Molecular Simulation Study of Montmorillonite in Contact with Water

Diffusivity of CO2 in H2O: A Review of Experimental Studies and Molecular Simulations in the Bulk and in Confinement

An in-depth review of the available experimental and molecular simulation studies of CO2 diffusion in H2O, which is a central property in important industrial and environmental processes, such as carbon capture and storage, enhanced oil recovery, and in the food industry is presented. The cases of both bulk and confined systems are covered. The experimental and molecular simulation data gathered are analyzed, and simple and computationally efficient correlations are devised. These correlations are applicable to conditions from 273 K and 0.1 MPa up to 473 K and 45 MPa. The available experimental data for diffusion coefficients of CO2 in brines are also collected, and their dependency on temperature, pressure, and salinity is examined in detail. Other engineering models and correlations reported in literature are also presented. The review of the simulation studies focuses on the force field combinations, the data for diffusivities at low and high pressures, finite-size effects, and the correlations developed based on the Molecular Dynamics data. Regarding the confined systems, we review the main methods to measure and compute the diffusivity of confined CO2 and discuss the main natural and artificial confining media (i.e., smectites, calcites, silica, MOFs, and carbon materials). Detailed discussion is provided regarding the driving force for diffusion of CO2 and H2O under confinement, and on the role of effects such as H2O adsorption on hydrophilic confining media on the diffusivity of CO2. Finally, an outlook of future research paths for advancing the field of CO2 diffusivity in H2O at the bulk phase and in confinement is laid out.

Microscopic Characterization of Deformation Behavior during Kerogen Evolution: Effects of Maturity and Skeleton Moisture Content

The CH4 storage and seepage capacity of shale kerogen are the main controlling factors of the natural gas production rate, and the porosity and permeability of kerogen are greatly affected by kerogen deformation. Therefore, the study of the deformation rule and CH4 adsorption characteristics of kerogen at different maturities and skeleton moisture contents has an important impact on the proper understanding of the development potential of shale gas reservoirs. In this paper, kerogen maturity (II-A, II-B, II-C, and II-D) and skeleton moisture content (0.0, 0.6, 1.2, 1.8, and 2.4 wt %) were considered. The deformation of kerogen, the adsorption of CH4 after deformation, and the quadratic deformation induced by CH4 were studied by using Grand Canonical Monte Carlo (GCMC) and molecular dynamics (MD). The results show that the kerogen volume strain increases with increasing skeleton moisture content, following the order II-A < II-B < II-C < II-D for the same moisture content. The density of the kerogen matrix decreases, and porosity increases with rising moisture content. The void fraction of immature kerogen decreases with increasing water content, while the opposite is true for postmature kerogen. The presence of skeleton moisture decreases the CH4 adsorption capacity of immature kerogen and increases the CH4 adsorption capacity of postmature kerogen. The chemical structure of immature kerogen is relatively soft, making its volume more affected by CH4 adsorption compared with postmature kerogen. In the same water environment, postmature kerogen has greater CH4 storage, diffusion, and seepage capacity compared to those of immature kerogen, suggesting that reservoirs with high organic matter maturity should be prioritized for development.

The Influence Mechanism of Interfacial Characteristics between CSH and Montmorillonite on the Strength Properties of Cement-Stabilized Montmorillonite Soil

To elucidate the impact mechanism of the interfacial characteristics of Calcium Silicate Hydrate gel (CSH)-Montmorillonite (MMT) at the nanoscale on the strength of cement-stabilized montmorillonite soil, this paper begins by examining the interfacial energy. Through Molecular Dynamics (MD) simulation methods, the energy at the MMT and CSH binding interface is quantitatively calculated, and the correlation between the interfacial energy and macroscopic strength is determined in conjunction with grey relational analysis. Finally, based on the characterization results from X-ray diffraction (XRD), the accuracy and sources of deviation in the MD simulation results are discussed. The study shows the CSH-MMT interfacial energy is composed of van der Waals forces, hydrogen bond energy, and electrostatic interactions, which are influenced by the migration of cations; there is a good consistency between the CSH-MMT interfacial energy and the unconfined compressive strength (UCS) of cement-stabilized soil (cemented soil), with the interfacial energy decreasing as the number of water molecules increases and first decreasing then increasing as the number of MMT layers grows; by adjusting the mix proportions, the magnitude of the CSH-MMT interfacial energy can be altered, thereby optimizing the strength of the cemented soil.

Microscopic insights into the variations of antibiotics sorption to clay minerals

Understanding the adsorption behavior of antibiotic molecules on minerals is crucial for determining the environmental fate and transport of antibiotics in soils and waters. However, the microscopic mechanisms that govern the adsorption of common antibiotics, such as the molecular orientation during the adsorption process and the conformation of sorbate species, are not well understood. To address this gap, we conducted a series of molecular dynamics (MD) simulations and thermodynamics analyses to investigate the adsorption of two typical antibiotics, tetracycline (TET) and sulfathiazole (ST), on the surface of montmorillonite. The simulation results indicated that the adsorption free energy ranged from - 23 to - 32 kJ·mol^-1, and - 9 to - 18 kJ·mol^-1 for TET and ST, respectively, which was consistent with the measured difference of sorption coefficient (K_d) for TET-montmorillonite of 11.7 L·g^-1 and ST-montmorillonite of 0.014 L·g^-1. The simulations also found that TET was adsorbed through dimethylamino groups (85% in probability) with a molecular conformation vertical to the montmorillonite's surface, while ST was adsorbed through sulfonyl amide group (95% in probability) with vertical, tilted and parallel conformations on the surface. The results confirmed that molecular spatial orientations could affect the adsorption capacity between antibiotics and minerals. Overall, the microscopic adsorption mechanisms revealed in this study provide critical insights into the complexities of antibiotics adsorption to soil and facilitate the prediction of adsorption capacity of antibiotics on minerals and their environmental transport and fate. This study contributes to our understanding of the environmental impacts of antibiotic usage and highlights the importance of considering molecular-level processes when assessing the fate and transport of antibiotics in the environment.

Control of the Structural Charge Distribution and Hydration State upon Intercalation of CO2 into Expansive Clay Interlayers

Numerous experimental investigations indicated that expansive clays such as montmorillonite can intercalate CO₂ preferentially into their interlayers and therefore potentially act as a material for CO₂ separation, capture, and storage. However, an understanding of the energy-structure relationship during the intercalation of CO₂ into clay interlayers remains elusive. Here, we use metadynamics molecular dynamics simulations to elucidate the energy landscape associated with CO₂ intercalation. Our free energy calculations indicate that CO₂ favorably partitions into nanoconfined water in clay interlayers from a gas phase, leading to an increase in the CO₂/H₂O ratio in clay interlayers as compared to that in bulk water. CO₂ molecules prefer to be located at the centers of charge-neutral hydrophobic siloxane rings, whereas interlayer spaces close to structural charges tend to avoid CO₂ intercalation. The structural charge distribution significantly affects the amount of CO₂ intercalated in the interlayers. These results provide a mechanistic understanding of CO₂ intercalation in clays for CO₂ separation, capture, and storage.

Interaction between Hydrated Smectite Clay Particles as a Function of Salinity (0-1 M) and Counterion Type (Na, K, Ca)

Swelling clay minerals control the hydrologic and mechanical properties of many soils, sediments, and sedimentary rocks. This important and well-known phenomenon remains challenging to predict because it emerges from complex multiscale couplings between aqueous chemistry and colloidal interaction mechanics in nanoporous clay assemblages, for which predictive models remain elusive. In particular, the predominant theory of colloidal interactions across fluid films, the widely used Derjaguin-Landau-Verwey-Overbeek model, fails to predict the ubiquitous existence of stable swelling states at interparticle distances below 3 nm that are stabilized by specific inter-atomic interactions in overlapping electrical double layers between the charged clay surfaces. Atomistic simulations have the potential to generate detailed insights into the mechanisms of these interactions. Recently, we developed a metadynamics-based molecular dynamics simulation methodology that can predict the free energy of interaction between parallel smectite clay particles in a wide range of interparticle distances (from 0.3 to 3 nm) and salinities (from 0.0 to 1.0 M NaCl). Here, we extend this work by characterizing the sensitivity of interparticle interactions to counterion type (Na, K, Ca). We establish a detailed picture of the free energy of interaction of parallel clay particles across water films as the sum of five interaction mechanisms with different sensitivities to salinity, counterion type, and interparticle distance.

Batch and Column Adsorption of Phosphorus by Modified Montmorillonite

Phosphorus pollutants are a crucial component of water eutrophication. In this study, montmorillonite modified by Keggin Al13 and hexadecyltrimethyl ammonium (Al13-O-MMt) was used as an adsorbent to remove phosphorus from solutions and thus simulate the practice of a field trial, such as in wastewater. The ammonium molybdate spectrophotometric method was used to determine the concentrations of phosphorus in samples. In the batch experiment, phosphorus was adsorbed by original montmorillonite (MMt) and Al13-O-MMt at various pH values (6–9) to identify the effect of pH during the adsorption process. The batch adsorption results demonstrate that Al13-O-MMt can adsorb up to 93% of phosphorus at pH = 8. Six graduated amounts (0.01–0.25 g) of montmorillonite were tested at three different temperatures to determine the most suitable temperature and the minimum dosage of Al13-O-MMt needed for the adsorption of 200 mg/L phosphorus in a 30 mL solution, which was 0.1 g at 25 °C. Therefore, the adsorption capacity of Al13-O-MMt was found to be 60 mg/g. Subsequently, a column experiment was conducted. The results showed that the optimized dosage of Al13-O-MMt was 6.667 g for phosphorus adsorption with a concentration of 200 mg/L in 2000 mL solution, and the breakthrough time was 4794.67 min.

Clay Swelling: Role of Cations in Stabilizing/Destabilizing Mechanisms

The stepwise hydration of clay minerals has been observed repeatedly in studies, but the underlying mechanism remains unclear. Previous numerical studies confirmed the presence of one-water layer (1W) and two-water layer (2W) hydration states. However, the undisturbed transition between these hydration states has never been captured. Using molecular dynamics simulation, this study (i) simulated for the first time the free 1W-2W transition during clay hydration and (ii) identified the underlying mechanism to be the detachment of cations from the clay surface and the formation of a shell of water molecules around the cation. The swelling dynamics of clay was found to be affected by the clay charge, clay mineralogy, and counterions through complex cation-clay interactions, cation hydration capacity, and cation migration rate.

Understanding the sorption behaviors of heavy metal ions in the interlayer and nanopore of montmorillonite: A molecular dynamics study

The molecular-scale adsorption mechanism of heavy metal ions in the interlayer and nanopore regions of montmorillonite (MMT) were investigated by molecular dynamics simulations. Three typical heavy metals (zinc, cadmium, and lead) were selected as the model ions, and two types of MMT (Arizona and Wyoming) were considered. The results showed that Cd²⁺ and Pb²⁺ can form both inner- and outer-sphere complexes on Wyoming MMT, while Zn²⁺ only formed outer-sphere complex due to the stronger hydration interaction of Zn²⁺ than Cd²⁺ and Pb²⁺. For Arizona MMT, all of the three ions only formed outer-sphere complexes on its interlayer and external basal surface in which the cations remained a fully hydrated state. The calculated diffusion coefficients of three cations in interlayer and nanopore indicated that their diffusion abilities were significantly impaired, implying that MMT adsorbents have a strong ability to fix and retard heavy metal ions. The derived results and mechanisms are instrumental to a profound understanding of the transport and retention of heavy metal elements in subsurface environments, and provide guidance for the management of heavy metal pollution.

Complexation of Alkali and Alkaline-Earth Metal Cations at Spodumene-Saltwater Interfaces by Molecular Simulation: Impact on Oleate Adsorption

Spodumene, a lithium aluminum inosilicate, is recovered by froth flotation using surfactants, so-called collectors. Therefore, the behavior and properties of the water-mineral interface in saline solutions are central. Here, molecular dynamics simulations are used to study the adsorption of alkali and alkaline-earth metal cations from concentrated solutions on the weakest (110) surface plane of negatively-charged spodumene. Results include the envelope density function of inner-sphere complexes for each cation and the density of complexes according to their adsorption contacts. Visualization of complexes for each cation is also included. Once the structure of the cation layers adsorbed on the surface of spodumene is defined, its role as a catalyst or barrier for adsorption of the spodumene collector in flotation is evaluated. The collector studied is the typical sodium oleate. The results show that oleate adsorption is poor and that the few adsorption contacts are mainly via cation bridges. The findings here indicate that molecular simulation can facilitate the search for effective collectors for environmentally sustainable spodumene flotation processes in saltwater.

Bulk and Interfacial Properties of the Decane + Water System in the Presence of Methane, Carbon Dioxide, and Their Mixture

Molecular dynamics simulations are carried out to study the two-phase behavior of the n-decane + water system in the presence of methane, carbon dioxide, and their mixture under reservoir conditions. The simulation studies were complemented by theoretical modeling using the perturbed-chain statistical associating fluid theory (PC-SAFT) equation of state (EoS) and density gradient theory. Our results show that the presence of methane and carbon dioxide decreases the interfacial tension (IFT) of the decane + water system. In general, the IFT increases with increasing pressure and decreasing temperature for the methane + decane + water and carbon dioxide + decane + water systems, similar to what has been found for the corresponding decane + water system. The most important finding of this study is that the presence of carbon dioxide decreases the IFT of the methane + decane + water system. The atomic density profiles provide evidence of the local accumulation of methane and carbon dioxide at the interface, in most of the studied systems. The results of this study show the preferential dissolution in the water-rich phase and enrichment at the interface for carbon dioxide in the methane + carbon dioxide + decane + water system. This indicates the preferential interaction of water with carbon dioxide relative to methane and decane. Notably, there is an enrichment of the interface by decane at high mole fractions of methane in the methane/decane-rich or methane/carbon dioxide/decane-rich phase. Overall, the solubility of methane and carbon dioxide in the water-rich phase increases with increasing pressure and temperature. Additionally, we find that the overall performance of the PC-SAFT EoS and the cubic-plus-association EoS is similar with respect to the calculation of bulk and interfacial properties of these systems.

Sorption and Diffusion of Methane and Carbon Dioxide in Amorphous Poly(alkyl acrylates): A Molecular Simulation Study

Molecular simulations were carried out to understand the structural features and the sorption and diffusion behavior of methane and carbon dioxide in amorphous poly(alkyl acrylates) in the temperature range of 300-600 K. The hybrid Monte Carlo/molecular dynamics approach was employed to address the effects of polymer swelling and framework flexibility on the gas sorption. Simulations show that the glass-transition temperature decreases with the side-chain length of poly(alkyl acrylate), consistent with experiments. This is due to the fact that the shielding of the polar ester groups increases with the side-chain length. The simulated sorption isotherms for methane and carbon dioxide were in agreement with the experimental data. The polymer swelling becomes more pronounced, especially in the case of sorption of carbon dioxide. A significant swelling occurs, possibly because of the greater interaction between carbon dioxide and the polar ester groups in the polymers. The uptake of methane and carbon dioxide by poly(alkyl acrylates) generally increases with the side-chain length. Our simulations confirm the experimental findings that the diffusion coefficients of methane and carbon dioxide in poly(alkyl acrylates) increase with the side-chain length. Interestingly, the activation energies of gas diffusion decrease with the side-chain length. The diffusion coefficients of the penetrants have an exponential relationship with the void fraction, which is in agreement with the free volume theory.

Swelling suppression of black cotton soil by means of liquid immersion and surface modification

Although the surface organic modification of smectite has been investigated widely, the swelling behavior of clays has been scarcely studied with consideration of civil engineering applications. In this work a facile strategy of liquid-immersion (dilute H₂SO₄ aqeuous solution) was proposed, and the 3-aminopropyltrimethoxysilane (APS) was utilized as surface modifier to suppress expansibility of black cotton soil (BCS) which is a type of highly swelling soils in tropical areas. Factors such as the incorporation dosage of APS, surface characters of soil treated by solution of H₂SO₄ or Na₂CO₃, and reaction temperatures/time were investigated to get lower swelling ratios. The treatment of BCS by H₂SO₄ was found more effective in immobilizing APS molecules, and hydronium ions were suppressed after the APS modification. The free swelling index (FSI) of BCS was decreased from 120% to 15% after treatment with H₂SO₄ and appropriate amount of APS modification. The reaction can be completed within several hours at the room temperature to ~80 °C. The soil samples were characterized by different means including the X-ray diffraction, Fourier-transform infrared spectroscopy, scanning electron microscope, thermogravimetric analysis and Zeta potential measurements. The APS molecules were found to react with –OH groups of the clay, and the hydrophobic groups provide surface hydrophobicity, which prevents hydration of cations within clay minerals. The APS was indicated to re-constructed lamellar structures of smectites after H₂SO₄ treatment, which suppressed the intra-crystalline and the subsequent osmotic swelling. This research highlights the liquid immersion and surface modification is applicable in diminishing swelling ratios of highly expansive black cotton soil.

Bulk and interfacial properties of decane in the presence of carbon dioxide, methane, and their mixture

Molecular dynamics simulations were performed to study the bulk and interfacial properties of methane + n-decane, carbon dioxide + n-decane, and methane + carbon dioxide + n-decane systems under geological conditions. In addition, theoretical calculations using the predictive Peng-Robinson equation of state and density gradient theory are carried out to compare with the simulation data. A key finding is the preferential dissolution in the decane-rich phase and adsorption at the interface for carbon dioxide from the methane/carbon dioxide mixture. In general, both the gas solubility and the swelling factor increase with increasing pressure and decreasing temperature. Interestingly, the methane solubility and the swelling of the methane + n-decane system are not strongly influenced by temperature. Our results also show that the presence of methane increases the interfacial tension (IFT) of the carbon dioxide + n-decane system. Typically, the IFT of the studied systems decreases with increasing pressure and temperature. The relatively higher surface excess of the carbon dioxide + n-decane system results in a steeper decrease in its IFT as a function of pressure. Such systematic investigations may help to understand the behavior of the carbon dioxide-oil system in the presence of impurities such as methane for the design and operation of carbon capture and storage and enhanced oil recovery processes.

Distribution and Mobility of Crude Oil-Brine in Clay Mesopores: Insights from Molecular Dynamics Simulations

The value of crude oil accommodated in shale has been recognized and has attracted increasing attention from the academic and industrial society. The occurrence and mobility of crude oil in clay pores, therefore, become essential issues for evaluation and recovery of shale oil. The distribution, structure, and transport of the oil-brine mixture confined in a slit-shaped montmorillonite mesopore with different water amounts have been investigated using equilibrium molecular dynamics and nonequilibrium molecular dynamics (NEMD) simulations. A mimic model of crude oil, a mixture of 19 organic molecules, was employed, and thus the behavior of different organic molecules could be characterized in detail. A temperature of 410 K and a pressure of 300 atm corresponding to a buried depth of 3 km were employed. The simulations indicate that the water amount determines the distribution of crude oil. Water and metal ions prefer to cover on hydrophilic montmorillonite surfaces, while nonpolar hydrocarbons tend to be far away from clay surfaces. As the water amount is too low to completely cover the clay surfaces, some polar organic molecules will come into contact with the uncovered clay surface. Abundant organic acid molecules adsorb onto montmorillonite surfaces mainly through participating in the inner-sphere complexes of Na⁺ ions closely located at montmorillonite surfaces (i.e., Na⁺ cation bridge) and forming hydrogen bonds with water molecules in the vicinity. Carbazole molecules tend to aggregate together due to π-π stacking, while thioether molecules mix within alkane molecules and exhibit no characteristic distributions. The mobility of all oil components decreases with the decrease of the water amount, and the mobility of polar components (i.e., organic acid and carbazole) is relatively lower than that of nonpolar hydrocarbons. NEMD simulations clearly indicate that the transport velocity of crude oil markedly increases with the water amount under a specific pressure gradient. The brine covering on clay surfaces significantly weakens oil-clay interfacial interactions. Polar components, especially organic acid, exhibit relatively low transport velocity compared with nonpolar hydrocarbons. These findings highlight the understanding of physical-chemical behaviors of shale oil and provide atomistic information for technology development for enhancing oil recovery.