Swelling Phenomena of the Nonswelling Clay Induced by CO2 and Water Cooperative Adsorption in Janus-Surface Micropores.

With the development of microscopy and sensor techniques, it becomes evident that nonswelling clays show swelling behavior under CO2-water mixture environments at high pressures and temperatures. The examples include Illite, muscovite, and kaolinite-rich rock samples. Here, we investigated the underlying mechanisms of kaolinite swelling induced by CO2 and water using molecular simulations and low-pressure gas adsorption experiments. The results suggest the cooperative adsorption behavior of CO2 and water on contact with kaolinite micropores, which have distinct wettabilities on the two adjoining interlayer surfaces. Even if clay-bound water exists, CO2 can enter the micropores to induce swelling. The measured micropore volume, simulated equilibrium stable interlayer distance with pure water, and that with CO2-water mixture were used in the swelling estimation, which shows good agreement with our experiments. The CO2 and water molecule distributions inside the interlayer micropores verify the importance of the wettabilities of the kaolinite surfaces in this cooperative adsorption behavior. The result extends the traditional understanding of the swelling mechanism, i.e., cation hydration and subsequent osmotic processes. In addition to earlier observations of kaolinite swelling behavior with potassium acetate, our study indicates the significance of the subtle balance of the noncovalent interactions between CO2, water, and the kaolinite Janus surfaces.

Modeling CO2-Water-Mineral Wettability and Mineralization for Carbon Geosequestration.

Carbon dioxide (CO2) capture and storage (CCS) is an important climate change mitigation option along with improved energy efficiency, renewable energy, and nuclear energy. CO2 geosequestration, that is, to store CO2 under the subsurface of Earth, is feasible because the world's sedimentary basins have high capacity and are often located in the same region of the world as emission sources. How CO2 interacts with the connate water and minerals is the focus of this Account. There are four trapping mechanisms that keep CO2 in the pores of subsurface rocks: (1) structural trapping, (2) residual trapping, (3) dissolution trapping, and (4) mineral trapping. The first two are dominated by capillary action, where wettability controls CO2 and water two-phase flow in porous media. We review state-of-the-art studies on CO2/water/mineral wettability, which was found to depend on pressure and temperature conditions, salt concentration in aqueous solutions, mineral surface chemistry, and geometry. We then review some recent advances in mineral trapping. First, we show that it is possible to reproduce the CO2/water/mineral wettability at a wide range of pressures using molecular dynamics (MD) simulations. As the pressure increases, CO2 gas transforms into a supercritical fluid or liquid at ∼7.4 MPa depending on the environmental temperature. This transition leads to a substantial decrease of the interfacial tension between CO2 and reservoir brine (or pure water). However, the wettability of CO2/water/rock systems depends on the type of rock surface. Recently, we investigated the contact angle of CO2/water/silica systems with two different silica surfaces using MD simulations. We found that contact angle increased with pressure for the hydrophobic (siloxane) surface while it was almost constant for the hydrophilic (silanol) surface, in excellent agreement with experimental observations. Furthermore, we found that the CO2 thin films at the CO2-hydrophilic silica and CO2-H2O interfaces displayed a linear correlation, which can in turn explain the constant contact angle on the hydrophilic silica surface. In view of the literature and our study results, a few recommendations seem necessary to construct a molecular system suitable to study wettability with MD simulations. Future work should be conducted to determine the influence of brine salinity on the wettability of minerals with high cation exchange capacity. Mineral trapping is believed to be an extremely slow process, likely taking thousands of years. However, a recent pilot study demonstrated that CO2 mineralization occurs within 2 years in highly reactive basalt reservoirs. A first-principles MD study has also shown that carbonation reactions occur rapidly at the surface oxygen sites of a reactive mineral. We observed carbonate ions on both a newly cleaved quartz surface (without hydrolysis), and a basalt andesine surface after hydrolysis in a CO2-rich environment. Future work should consider the influence of water, gas impurities, and mineral cation type on carbonation.

Elasticity and Stability of Clathrate Hydrate: Role of Guest Molecule Motions.

Molecular dynamic simulations were performed to determine the elastic constants of carbon dioxide (CO2) and methane (CH4) hydrates at one hundred pressure-temperature data points, respectively. The conditions represent marine sediments and permafrost zones where gas hydrates occur. The shear modulus and Young's modulus of the CO2 hydrate increase anomalously with increasing temperature, whereas those of the CH4 hydrate decrease regularly with increase in temperature. We ascribe this anomaly to the kinetic behavior of the linear CO2 molecule, especially those in the small cages. The cavity space of the cage limits free rotational motion of the CO2 molecule at low temperature. With increase in temperature, the CO2 molecule can rotate easily, and enhance the stability and rigidity of the CO2 hydrate. Our work provides a key database for the elastic properties of gas hydrates, and molecular insights into stability changes of CO2 hydrate from high temperature of ~5 °C to low decomposition temperature of ~-150 °C.

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.

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.

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.

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.

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.

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.

Comparison of thermodynamic stabilities and mechanical properties of CO2, SiO2, and GeO2 polymorphs by first-principles calculations.

The recent discovery that molecular CO(2) transforms under compression into carbon four-coordinated, 3-dimensional network solid phases has generated considerable interests on possible new phases in the fourth-main-group elemental oxides. Based on density-functional theory calculations, we have investigated the thermodynamic stability, mechanical properties and electronic structure of proposed guest-free clathrates, quartz and cristobalite phases for CO(2), SiO(2), and GeO(2), and the dry ice phase for CO(2). It was predicted that a GeO(2) clathrate, likely a semiconductor, could be synthesized presumably with some suitable guest molecules. The hypothetical CO(2) guest-free clathrate phase was found hardly to be formed due to the large energy difference with respect to the other polymorphs. This phase is unstable at all pressures, which is also implied by its different electronic structure in comparison with SiO(2) and GeO(2). Finally, the SiO(2) clathrate presents a uniquely high bulk modulus, which is higher than that of quartz and three times of the experimental data, might not be a weak point of ab-initio calculations such as pseudopotentials, correlation functional etc., instead it can be readily understood by the constraint as imposed by the high symmetry. Either temperature or an "exhausted" relaxation (without any symmetry constraint) can remedy this problem.

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.