Yunfeng Liang

Self-accumulation of aromatics at the oil-water interface through weak hydrogen bonding.

It is well-known that the amphiphilic solutes are surface-active and can accumulate at the oil-water interface. Here, we have investigated the water and a light-oil model interface by using molecular dynamic simulations. It was found that aromatics concentrated in the interfacial region, whereas the other hydrocarbons were uniformly distributed throughout the oil phase. Similar to previous studies, such concentrations were not observed at pure aromatics-water interfaces. We show that the self-accumulation of aromatics at the oil-water interface is driven by differences in the interfacial tension, which is lower for aromatics-water than between the others. The weak hydrogen bonding between the aromatic rings and the water protons provides the mechanism for lowering the interfacial tension.

Molecular dynamics study of salt–solution interface: Solubility and surface charge of salt in water

The NaCl salt-solution interface often serves as an example of an uncharged surface. However, recent laser-Doppler electrophoresis has shown some evidence that the NaCl crystal is positively charged in its saturated solution. Using molecular dynamics (MD) simulations, we have investigated the NaCl salt-solution interface system, and calculated the solubility of the salt using the direct method and free energy calculations, which are kinetic and thermodynamic approaches, respectively. The direct method calculation uses a salt-solution combined system. When the system is equilibrated, the concentration in the solution area is the solubility. In the free energy calculation, we separately calculate the chemical potential of NaCl in two systems, the solid and the solution, using thermodynamic integration with MD simulations. When the chemical potential of NaCl in the solution phase is equal to the chemical potential of the solid phase, the concentration of the solution system is the solubility. The advantage of using two different methods is that the computational methods can be mutually verified. We found that a relatively good estimate of the solubility of the system can be obtained through comparison of the two methods. Furthermore, we found using microsecond time-scale MD simulations that the positively charged NaCl surface was induced by a combination of a sodium-rich surface and the orientation of the interfacial water molecules.

Lattice Boltzmann simulations of supercritical CO2-water drainage displacement in porous media: CO2 saturation and displacement mechanism.

CO2 geosequestration in deep aquifers requires the displacement of water (wetting phase) from the porous media by supercritical CO2 (nonwetting phase). However, the interfacial instabilities, such as viscous and capillary fingerings, develop during the drainage displacement. Moreover, the burstlike Haines jump often occurs under conditions of low capillary number. To study these interfacial instabilities, we performed lattice Boltzmann simulations of CO2-water drainage displacement in a 3D synthetic granular rock model at a fixed viscosity ratio and at various capillary numbers. The capillary numbers are varied by changing injection pressure, which induces changes in flow velocity. It was observed that the viscous fingering was dominant at high injection pressures, whereas the crossover of viscous and capillary fingerings was observed, accompanied by Haines jumps, at low injection pressures. The Haines jumps flowing forward caused a significant drop of CO2 saturation, whereas Haines jumps flowing backward caused an increase of CO2 saturation (per injection depth). We demonstrated that the pore-scale Haines jumps remarkably influenced the flow path and therefore equilibrium CO2 saturation in crossover domain, which is in turn related to the storage efficiency in the field-scale geosequestration. The results can improve our understandings of the storage efficiency by the effects of pore-scale displacement phenomena.

Asymmetric orientation of toluene molecules at oil-silica interfaces.

The interfacial structure of heptane and toluene at oil-silica interfaces has previously been studied by sum frequency generation [Z. Yang et al., J. Phys. Chem. C. 113, 20355 (2009)]. It was found that the toluene molecule is almost perpendicular to the silica surface with a tilt angle of about 25°. Here, we have investigated the structural properties of toluene and heptane at oil-silica interfaces using molecular dynamics simulations for two different surfaces: the oxygen-bridging (hydrophobic) and hydroxyl-terminated (hydrophilic) surfaces of quartz (silica). Based on the density profile, it was found that both heptane and toluene oscillate on silica surfaces, with heptane showing more oscillation peaks. Furthermore, the toluene molecules of the first layer were found to have an asymmetric distribution of orientations, with more CH(3) groups pointed away from the silica surface than towards the silica surface. These findings are generally consistent with previous experiments, and reveal enhanced molecular structures of liquids at oil-silica interfaces.

Ion Distribution and Hydration Structure in the Stern Layer on Muscovite Surface

Based on molecular dynamics simulations of eight ions (Na+, K+, Rb+, Cs+, Mg2+, Ca2+, Sr2+, and Ba2+) on muscovite mica surfaces in water, we demonstrate that experimental data on the muscovite mica surface can be rationalized through a unified picture of adsorption structures including the hydration structure, cation heights from the muscovite surface, and state stability. These simulations enable us to categorize the inner-sphere surface complex into two different species: an inner-sphere surface complex in a ditrigonal cavity (IS1) and that on top of Al (IS2). By considering the presence of the two inner-sphere surface complexes, the experimental finding that the heights of adsorbed cations from the muscovite surface are proportional to the ionic radius for K+ and Cs+ but inversely proportional to the ionic radius for Ca2+ and Ba2+ was explained. We find that Na+, Ca2+, Sr2+, and Ba2+ can form both IS1 and IS2; K+, Rb+, and Cs+ can form only IS1; and Mg2+ can form only IS2. It is suggested that the formation of IS1 and IS2 is governed by the charge density of the ions. Among the eight ions, we also find that the hydration structure for the outer-sphere surface complexes of divalent cations differs from that of the monovalent cations by one adsorbed water molecule (i.e., a water molecule located in a ditrigonal cavity).

Molecular Dynamics Simulation of Atomic Force Microscopy at the Water–Muscovite Interface: Hydration Layer Structure and Force Analysis

With the development of atomic force microscopy (AFM), it is now possible to detect the buried liquid-solid interfacial structure in three dimensions at the atomic scale. One of the model surfaces used for AFM is the muscovite surface because it is atomically flat after cleavage along the basal plane. Although it is considered that force profiles obtained by AFM reflect the interfacial structures (e.g., muscovite surface and water structure), the force profiles are not straightforward because of the lack of a quantitative relationship between the force and the interfacial structure. In the present study, molecular dynamics simulations were performed to investigate the relationship between the muscovite-water interfacial structure and the measured AFM force using a capped carbon nanotube (CNT) AFM tip. We provide divided force profiles, where the force contributions from each water layer at the interface are shown. They reveal that the first hydration layer is dominant in the total force from water even after destruction of the layer. Moreover, the lateral structure of the first hydration layer transcribes the muscovite surface structure. It resembles the experimentally resolved surface structure of muscovite in previous AFM studies. The local density profile of water between the tip and the surface provides further insight into the relationship between the water structure and the detected force structure. The detected force structure reflects the basic features of the atomic structure for the local hydration layers. However, details including the peak-peak distance in the force profile (force-distance curve) differ from those in the density profile (density-distance curve) because of disturbance by the tip.

Mechanisms for Enhanced Hydrophobicity by Atomic-Scale Roughness.

It is well known that the close-packed CF3-terminated solid surface is among the most hydrophobic surfaces in nature. Molecular dynamic simulations show that this hydrophobicity can be further enhanced by the atomic-scale roughness. Consequently, the hydrophobic gap width is enlarged to about 0.6 nm for roughened CF3-terminated solid surfaces. In contrast, the hydrophobic gap width does not increase too much for a rough CH3-terminated solid surface. We show that the CF3-terminated surface exists in a microscopic Cassie–Baxter state, whereas the CH3-terminated surface exists as a microscopic Wenzel state. This finding elucidates the underlying mechanism for the different widths of the observed hydrophobic gap. The cage structure of the water molecules (with integrated hydrogen bonds) around CH3 terminal assemblies on the solid surface provides an explanation for the mechanism by which the CH3-terminated surface is less hydrophobic than the CF3-terminated surface.

Microscopic Origin of Strain Hardening in Methane Hydrate

It has been reported for a long time that methane hydrate presents strain hardening, whereas the strength of normal ice weakens with increasing strain after an ultimate strength. However, the microscopic origin of these differences is not known. Here, we investigated the mechanical characteristics of methane hydrate and normal ice by compressive deformation test using molecular dynamics simulations. It is shown that methane hydrate exhibits strain hardening only if the hydrate is confined to a certain finite cross-sectional area that is normal to the compression direction. For normal ice, it does not present strain hardening under the same conditions. We show that hydrate guest methane molecules exhibit no long-distance diffusion when confined to a finite-size area. They appear to serve as non-deformable units that prevent hydrate structure failure, and thus are responsible for the strain-hardening phenomenon.

Molecular Dynamics Study of Aqueous NaCl Solutions: Flash Crystallization Caused by Solution Phase Change

Solutions under sub- or supercritical conditions receive much attention because of their significance in geology and industrial applications. One of the most important properties is their very low salt solubility, which leads to rapid crystallization. The morphology of sodium chloride crystals formed in supercritical fluids has been shown to be affected by the presence of the high-concentration liquid phase and low-concentration vapor phase [Armellini and Tester, J. Supercritical Fluid 4, 254–264 (1991)]. However, because of the short time scales involved, experimental observation of the underlying mechanism is difficult. In the present study, microsecond-scale molecular dynamics, which provide insight with picosecond resolution, were conducted for the NaCl salt–solution interfacial system at sub- and super-critical conditions. We propose the utilization of the correlation between the number density of ions and water to parameterize the solution phase. This correlation and the two-dimensional number densities provide insight into flash crystallization at the atomic scale. Vapor–liquid phase coexistence was found at high pressure and the diffusive high-concentration liquid phase could transport sodium and chloride ions to form a compact solid phase. In contrast, an isolated crystal forms at low pressure owing to the rapid volume expansion of the vapor phase.

Microscopic mechanisms of the pressure-induced amorphization of SiO2

When the low-pressure tetrahedral phases of SiO2 are compressed at room temperature they retain their structures at pressures well above their stability domain, but they all collapse into denser structures when pressure reaches about 20 GPa. Depending on the experimental conditions, pressure-induced densification can be accompanied by amorphization; by the formation of crystalline, metastable polymorphs; and can be preceded by the appearance of intermediate phase. Based on molecular dynamic simulations, we show that this rich phenomenology can be rationalized through a unified theoretical framework of the atomistic pathways leading to densification. The model emphasizes the role played by the oxygen sublattice, which transforms from a bcc-like order into close-packed arrangements in the denser structures, through a ferroelastic instability of martensitic nature.