Molecular-Level Understanding of Surface Roughness Boosting Segregation Behavior at the ZIF-8/Ionic Liquid Interfaces

Here, we perform a series of classical molecular dynamics simulations for two different [HEMIM][DCA] and [BMIM][BF₄] ionic liquids (ILs) on the ZIF-8 surface to explore the interfacial properties of metal-organic framework (MOFs)/IL composite materials at the molecular level. Our simulation results reveal that the interfacial structures of anions and cations on the ZIF-8 surface are dominated by the surface roughness due to the steric hindrance, which is extremely different from the driving mechanism based on solid-ion interactions of ILs on flat solid surfaces. At the ZIF-8/IL interfaces, the open sodalite (SOD) cages of the ZIF-8 surface can block most of the large-size cations outside and significantly boost the segregation behavior of anions and cations. In comparison with the [BMIM][BF₄] IL, the [HEMIM][DCA] IL has much more anions entering into the open SOD cages owing to the combination of stronger ZIF-8-[DCA]^- interactions and more ordered arrangement of [DCA]^- anions on the ZIF-8 surface. Furthermore, more and stronger ZIF-8-[BF₄]^- hydrogen bonds (HBs) are found to exist on the cage edges than the ZIF-8-[DCA]^- HBs, further preventing [BF₄]^- anions from entering into SOD cages. By more detailed analyses, we find that the hydrophobic interaction has an important influence on the interfacial structures of the side chains of [HEMIM]⁺ and [BMIM]⁺ cations, while the π-π stacking interaction plays a key role in determining the interfacial structures of the imidazolium rings of both cations. Our simulation results in this work provide a molecular-level understanding of the underlying driving mechanism on segregation behavior at the ZIF-8/IL interfaces.

Local Structural Fingerprints of Nanoparticle-Bound Polymer Layers

The molecular structure of bound layers at attractive polymer-nanoparticle interfaces strongly influences the properties of nanocomposites. Thus, a unifying theoretical framework that can provide insights into the correlations between the molecular structure of the bound layers, their thermodynamics, and macroscopic properties is highly desirable. In this work, molecular dynamics simulations were used in combination with local fingerprint analysis of configurational entropy and interaction energy at the segmental scale, with the goal to establish such physical grounds. The thickness of bound polymer layers is found to be independent of the polymer chain length, as deduced from density oscillations at the surface of a nanotube. The local configurational entropy of layers is estimated from pair correlations in equilibrium structures. By plotting mean layer entropy vs internal energy on a phase diagram, a one-to-one equivalence is established between the local structures of layers and their thermodynamic properties. Moreover, a gradient in local dynamics of segments in bound layers under equilibrium is observed normal to the nanoparticle surface. The relaxation times of individual layers show correspondence to their phase diagram fingerprints, thus suggesting that a unified perspective can be envisioned for such materials built on the grounds of locally heterogeneous interfaces.

Single-Particle Dynamics at the Intrinsic Surface of Aqueous Alkali Halide Solutions

The distribution of ions in the proximity of the liquid-vapor interface of their aqueous solution has been the subject of an intense debate during the last decade. The effects of ionic polarizability have been one of its salient aspects. Much less has been said about the corresponding dynamical properties, which are substantially unexplored. Here, we investigate the single-particle dynamics at the liquid-vapor interface of several alkali halide solutions, using molecular dynamics simulations with polarizable and nonpolarizable force fields and intrinsic surface analysis. We analyze the diffusion coefficient, residence time, and velocity autocorrelation function of water and ions and investigate how these properties depend on the molecular layer where they reside. While anions are found in the first molecular layer for relatively long times, cations are only making quick excursions into it, thanks to thermal fluctuations. The in-layer residence time of ions and their molar fraction in the layer turned out to be linearly dependent on each other. We interpret this unexpected result using a simple two-state model. In addition, we found that, unlike water and other neat molecular liquids that show a different diffusion mechanism at the surface than in the bulk of their liquid phase, ions do not enjoy enhanced mobility in the surface layer of their aqueous solution. This result indicates that ions in the surface layer are shielded by their nearest water neighbors from being exposed to the vapor phase as much as possible. Such positions are available for the ions at the negatively curved troughs of the molecularly rugged liquid surface.

Finite-size corrections for the static structure factor of a liquid slab with open boundaries

The presence of a confining boundary can modify the local structure of a liquid markedly. In addition, small samples of finite size are known to exhibit systematic deviations of thermodynamic quantities relative to their bulk values. Here, we consider the static structure factor of a liquid sample in slab geometry with open boundaries at the surfaces, which can be thought of as virtually cutting out the sample from a macroscopically large, homogeneous fluid. This situation is a relevant limit for the interpretation of grazing-incidence diffraction experiments at liquid interfaces and films. We derive an exact, closed expression for the slab structure factor, with the bulk structure factor as the only input. This shows that such free boundary conditions cause significant differences between the two structure factors, in particular, at small wavenumbers. An asymptotic analysis of this result yields the scaling exponent and an accurate, useful approximation of these finite-size corrections. Furthermore, the open boundaries permit the interpretation of the slab as an open system, supporting particle exchange with a reservoir. We relate the slab structure factor to the particle number fluctuations and discuss conditions under which the subvolume of the slab represents a grand canonical ensemble with chemical potential μ and temperature T. Thus, the open slab serves as a test-bed for the small-system thermodynamics in a μT reservoir. We provide a microscopically justified and exact result for the size dependence of the isothermal compressibility. Our findings are corroborated by simulation data for Lennard-Jones liquids at two representative temperatures.

Amphiphilic Dimers at Liquid-Liquid Interfaces: A Density Functional Approach

We apply density functional theory to study the structure of dimers at the interface between two partially miscible symmetric liquids. The dimers are built of two tangentially jointed spheres and do not solve the coexisting liquids. The interactions in the system are modeled using Lennard-Jones potentials with different interactions between segments of the dimers and the liquid components. We study how asymmetry of the interactions between dimers and molecules of the liquid, i.e., the degree of dimer amphiphilicity, influences the interfacial structure. Two unexpected phenomena have been found. First, for some systems, the liquid-liquid interface is able to accommodate only a finite amount of dimers. If the amount of added dimers is larger than a threshold value, a part or all of the dimers move to the interior one of the coexisting phase, forming an insoluble sheet inside it, or the initial interface splits into separate parts. The second is a peculiar behavior of the dependence of the interfacial width with an increase of the amount of added dimers. In this case, we observe a discontinuous jump that is connected with reorientation of dimers with respect to the interface.

Characterizing the bifurcating configuration of hydrogen bonding network in interfacial liquid water and its adhesion on solid surfaces

The interfacial structures of liquid water molecules adjacent to a solid surface contribute significantly to the interfacial properties of aqueous solutions, and are of prime importance in a wide spectrum of applications. In this work, we use molecular dynamics (MD) simulations to explore the interfacial structures, mainly in term of hydrogen bonding network, of a liquid water film interacting intimately with solid surfaces, which are composed of [100] face centered cubic (FCC) lattices. We disclose the formation of a bifurcating configuration of hydrogen bonds in interfacial liquid water and ascribe its occurrence to the collective effects of water density depletion, hydrogen bonds and local polarization. Such bifurcating configuration of interfacial water molecules consists of repetitive layer by layer water sheets with intra-layer hydrogen bonding network being formed in each layer, and inter-layer defects, i.e., hydrogen bonds formed between two neighboring layers of interfacial water. A lower bound of 2.475 for the average number of hydrogen bonds per interfacial water molecule is expected. Our MD study on the interfacial configuration of water on solid surfaces reveals a quadratic dependence of adhesion on the solid-liquid affinity, bridging the gap between the macroscopic interfacial property W adh and the microscopic parameter ε SL of the depth of the Lennard-Jones solid-liquid potential.

Effect of Polymer/Solid and Polymer/Vapor Instantaneous Interfaces on the Interfacial Structure and Dynamics of Polymer Melt Systems

Polymers are used in a wide range of applications that involve chemical and physical processes taking place at surfaces or interfaces which influence the interaction between the polymer material and the substance that comes into contact with it. Polymer surfaces are usually modified either chemically or physically for specific applications such as facilitating wetting, reducing friction, and enhancing adhesion. The variety and complexity of surface and interfacial processes requires a molecular-level understanding of the structural and dynamical properties of the surface/interface layer to help in the design of materials with desired functional properties. Using molecular dynamics (MD) simulations, we investigate the structure and dynamics at the surface of polymer films. We find that the density profiles of the films as a function of distance relative to an instantaneous surface have a structure indicative of a layering at the polymer/vapor interface similar to the typical layered structure observed at the polymer/substrate interface. However, the interfacial molecules at the polymer/vapor interface have a higher mobility compared to that in the bulk while the mobility of the molecules is lower at the polymer/substrate interface. Time correlation of the instantaneous polymer/vapor interface shows that surface fluctuations are strongly temperature dependent and are directly related to the mobility of polymer chains near the interface.

Properties of the liquid–vapor interface of acetone–methanol mixtures, as seen from computer simulation and ITIM surface analysis

The intrinsic surface of acetone–methanol mixtures is studied by computer simulation and ITIM analysis.

Two-dimensional percolation at the free water surface and its relation with the surface tension anomaly of water

The percolation temperature of the lateral hydrogen bonding network of the molecules at the free water surface is determined by means of molecular dynamics computer simulation and identification of the truly interfacial molecules analysis for six different water models, including three, four, and five site ones. The results reveal that the lateral percolation temperature coincides with the point where the temperature derivative of the surface tension has a minimum. Hence, the anomalous temperature dependence of the water surface tension is explained by this percolation transition. It is also found that the hydrogen bonding structure of the water surface is largely model-independent at the percolation threshold; the molecules have, on average, 1.90 ± 0.07 hydrogen bonded surface neighbors. The distribution of the molecules according to the number of their hydrogen bonded neighbors at the percolation threshold also agrees very well for all the water models considered. Hydrogen bonding at the water surface can be well described in terms of the random bond percolation model, namely, by the assumptions that (i) every surface water molecule can form up to 3 hydrogen bonds with its lateral neighbors and (ii) the formation of these hydrogen bonds occurs independently from each other.

First principles nonequilibrium plasma mixing

We have performed nonequilibrium classical and quantum-mechanical molecular dynamics simulations that follow the interpenetration of deuterium-tritium (DT) and carbon (C) components through an interface initially in hydrostatic and thermal equilibrium. We concentrate on the warm, dense matter regime with initial densities of 2.5-5.5 g/cm3 and temperatures from 10 to 100 eV. The classical treatment employs a Yukawa pair-potential with the parameters adjusted to the plasma conditions, and the quantum treatment rests on an orbital-free density functional theory at the Thomas-Fermi-Dirac level. For times greater than about a picosecond, the component concentrations evolve in accordance with Fick's law for a classically diffusing fluid with the motion, though, described by the mutual diffusion coefficient of the mixed system rather than the self-diffusion of the individual components. For shorter times, microscopic processes control the clearly non-Fickian dynamics and require a detailed representation of the electron probability density in space and time.

Calculation of the intrinsic solvation free energy profile of an ionic penetrant across a liquid-liquid interface with computer simulations

We introduce the novel concept of an intrinsic free energy profile, allowing one to remove the artificial smearing caused by thermal capillary waves, which renders difficulties for the calculation of free energy profiles across fluid interfaces in computer simulations. We apply this concept to the problem of a chloride ion crossing the interface between water and 1,2-dichloroethane and show that the present approach is able to reveal several important features of the free energy profile which are not detected with the usual, nonintrinsic calculations. Thus, in contrast to the nonintrinsic profile, a free energy barrier is found at the aqueous side of the (intrinsic) interface, which is attributed to the formation of a water "finger" the ion pulls with itself upon approaching the organic phase. Further, by the presence of a nonsampled region, the intrinsic free energy profile clearly indicates the coextraction of the first hydration shell water molecules of the ion when entering the organic phase.

Molecular level properties of the free water surface and different organic liquid/water interfaces, as seen from ITIM analysis of computer simulation results

Molecular dynamics simulations of the interface of water with four different apolar phases, namely water vapour, liquid carbon tetrachloride, liquid dichloromethane (DCM) and liquid dichloroethane (DCE) are performed on the canonical ensemble at 298 K. The resulting configurations are analysed using the novel method of identification of the truly interfacial molecules (ITIM). Properties of the first three molecular layers of the liquid phases (e.g. width, spacing, roughness, extent of the in-layer hydrogen bonding network) as well as of the molecules constituting these layers (e.g., dynamics, orientation) are investigated in detail. In the analyses, particular attention is paid to the effect of the polarity of the non-aqueous phase and to the length scale of the effect of the vicinity of the interface on the various properties of the molecules. The obtained results show that increasing polarity of the non-aqueous phase leads to the narrowing of the interface, in spite of the fact that, at the same time, the truly interfacial layer of water gets somewhat broader. The influence of the nearby interface is found to extend only to the first molecular layer in many respects. This result is attributed to the larger space available for the truly interfacial than for the non-interfacial molecules (as the shapes of the two liquid surfaces are largely independent of each other, resulting in the presence of voids between the two phases), and to the fact that the hydrogen bonding interaction of the truly interfacial water molecules with other waters is hindered in the direction of the interface.

Structure and charging of hydrophobic material/water interfaces studied by phase-sensitive sum-frequency vibrational spectroscopy

We have studied the hydrophobic water/octadecyltrichlorosilane (OTS) interface by using the phase-sensitive sum-frequency vibrational spectroscopy (PS-SFVS), and we obtained detailed structural information of the interface at the molecular level. Excess ions emerging at the interface were detected by changes of the surface vibrational spectrum induced by the surface field created by the excess ions. Both hydronium (H(3)O(+)) and hydroxide (OH(-)) ions were found to adsorb at the interface, and so did other negative ions such as Cl(-). By varying the ion concentrations in the bulk water, their adsorption isotherms were measured. It was seen that among the three, OH(-) has the highest adsorption energy, and H(3)O(+) has the lowest; OH(-) also has the highest saturation coverage, and Cl(-) has the lowest. The result shows that even the neat water/OTS interface is not neutral, but charged with OH(-) ions. The result also explains the surprising observation that the isoelectric point appeared at approximately 3.0 when HCl was used to decrease the pH starting from neat water.