Molecular simulation of aqueous solutions at clay surfaces

A Coarse-Grained Interaction Model for Sodium Dominant Montmorillonite

Montmorillonite is the main crystalline mineral present in bentonite. It is an absorbent, swelling material; the physical chemistry underlying its ability to absorb water and swell occurs at the nanoscale, governed by electrical double-layer interactions. In turn, absorption and swelling lead to important changes in the macroscopic transport properties of the clay. Mesoscale models can help us establish a link between these nanoscale processes and macroscale properties, notably by providing a detailed description of its pore network. Models on the scale of hundreds to thousands of nanometers are required, which cannot realistically be handled using traditional all-atom molecular dynamics simulations. This work presents a coarse-grained (CG) mesoscale model of sodium montmorillonite. In our model, montmorillonite platelets are represented by two types of particles: central nonhydrogen-bonded particles and edge hydrogen-bonding particles. The particle interactions are described by two-body potentials, which were optimized based on all-atom molecular dynamics simulations. Specifically, several potential mean force calculations involving dry and hydrated montmorillonite were performed, using the ClayFF potential to calculate interatomic forces. The CG model was validated by testing the scalability of the model, testing its ability to reproduce potentials of mean force reported elsewhere in the literature, and by comparing the calculated elastic properties of a system containing 1000 Na montmorillonite platelets to experimentally measured elastic properties of bentonite. The simulated elastic properties obtained using our mesoscale model agree with these experimental values.

Molecular Dynamics Simulation of Ion Adsorption and Ligand Exchange on an Orthoclase Surface

Orthoclase (K-feldspar) is one of the natural inorganic materials, which shows remarkable potential toward removing heavy metal ions from aqueous solutions. Understanding the interactions of the orthoclase and metal ions is important in the treatment of saline wastewater. In this paper, molecular dynamics simulations were used to prove the adsorption of different ions onto orthoclase. The adsorption isotherms show that orthoclase has remarkable efficiency in the removal of cations at low ion concentrations. Aluminol groups are the preferential adsorption sites of cations due to higher negative charges. The adsorption types and adsorption sites are influenced by the valence, radius, and hydration stability of ions. Monovalent cations can be adsorbed in the cavities, whereas divalent cations cannot. The hydrated cation may form an outer-sphere complex or an inner-sphere complex in association with the loss of hydration water. Na+, K+, and Ca2+ ions mainly undergo inner-sphere adsorption and Mg2+ ions prefer outer-sphere adsorption. On the basis of simulation results, the mechanism of ion removal in the presence of orthoclase is demonstrated at a molecular level.

Diffusion of confined fluids in microporous zeolites and clay materials

Fluids exhibit remarkable variation in their structural and dynamic properties when they are confined at the nanoscopic scale. Various factors, including geometric restriction, the size and shape of the guest molecules, the topology of the host, and guest-host interactions, are responsible for the alterations in these properties. Due to their porous structures, aluminosilicates provide a suitable host system for studying the diffusion of sorbates in confinement. Zeolites and clays are two classes of the aluminosilicate family, comprising very ordered porous or layered structures. Zeolitic materials are important due to their high catalytic activity and molecular sieving properties. Guest molecules adsorbed by zeolites display many interesting features including unidimensional diffusion, non-isotropic rotation, preferred orientation and levitation effects, depending on the guest and host characteristics. These are useful for the separation of hydrocarbons which commonly exist as mixtures in nature. Similarly, clay materials have found application in catalysis, desalination, enhanced oil recovery, and isolation barriers used in radioactive waste disposal. It has been shown that the bonding interactions, level of hydration, interlayer spacing, and number of charge-balancing cations are the important factors that determine the nature of diffusion of water molecules in clays. Here, we present a review of the current status of the diffusion mechanisms of various adsorbed species in different microporous zeolites and clays, as investigated using quasielastic neutron scattering and classical molecular dynamics simulation techniques. It is impossible to write an exhaustive review of the subject matter, as it has been explored over several decades and involves many research topics. However, an effort is made to cover the relevant issues specific to the dynamics of different molecules in microporous zeolites and clay materials and to highlight a variety of interesting features that are important for both practical applications and fundamental aspects.

Sorption Mechanisms of Chemicals in Soils

Sorption of chemicals onto soil particle surfaces is an important process controlling their availability for uptake by organisms and loss from soils to ground and surface waters. The mechanisms of chemical sorption are inner- and outer-sphere adsorption and precipitation onto mineral surfaces. Factors that determine the sorption behavior are properties of soil mineral and organic matter surfaces and properties of the sorbing chemicals (including valence, electron configuration, and hydrophobicity). Because soils are complex heterogeneous mixtures, measuring sorption mechanisms is challenging; however, advancements analytical methods have made direct determination of sorption mechanisms possible. In this review, historical and modern research that supports the mechanistic understanding of sorption mechanisms in soils is discussed. Sorption mechanisms covered include cation exchange, outer-sphere adsorption, inner-sphere adsorption, surface precipitation, and ternary adsorption complexes.

XPS and two-dimensional FTIR correlation analysis on the binding characteristics of humic acid onto kaolinite surface

The stabilization and preservation of soil organic matter have been attributed to the strong reactive sites of mineral surfaces that cause physical isolation and chemical stabilization due to the organic-mineral interface. However, much of the micro-scale knowledge about interactions between organic ligands and minerals largely remains at the qualitative level, and neglects the heterogeneity of functional groups of organic matter. Here, we report the use of molecular-scale technologies of two-dimensional FTIR Correlation Spectroscopy (2D-FTIR-CoS) and X-ray Photoelectron Spectroscopy (XPS) to directly measure the binding processes of humic acid (JGHA) groups onto kaolinite surface. The spectroscopy results showed that the carboxylate groups, aliphatic OH and aromatic structure participate in the binding of JGHA on kaolinite surface. The carboxylic and phenolic hydroxyl interact with kaolinite surface through the interfacial COAl/Si bonds. Kaolinite prefers to adsorb C-groups at pH 4.0 and O-groups at pH 8.0. The interaction of COO^- group at 1566 cm^-1 of JGHA leads to the formation of inner-sphere complex first and then outer-sphere complex with increasing contact time. The interaction of COOH group at 1261 cm^-1 with the AlOH²⁺ of kaolinite was could be ascribed to ligand exchange and/or electrostatic attraction, whose contribution was evaluated to be 13.90%, 7.65% and 0% at pH 4.0, 6.0 and 8.0, respectively. These results of molecular binding provide quantitative mechanistic insights into organic-mineral interactions and expound the effect of functional groups of HA on binding mechanisms, and thus bring important clues for better understanding the mobility and transformation of land‑carbon including mineral-bound carbon.

Atomistic Structure of Mineral Nano-aggregates from Simulated Compaction and Dewatering

The porosity of clay aggregates is an important property governing chemical reactions and fluid flow in low-permeability geologic formations and clay-based engineered barrier systems. Pore spaces in clays include interlayer and interparticle pores. Under compaction and dewatering, the size and geometry of such pore spaces may vary significantly (sub-nanometer to microns) depending on ambient physical and chemical conditions. Here we report a molecular dynamics simulation method to construct a complex and realistic clay-like nanoparticle aggregate with interparticle pores and grain boundaries. The model structure is then used to investigate the effect of dewatering and water content on micro-porosity of the aggregates. The results suggest that slow dewatering would create more compact aggregates compared to fast dewatering. Furthermore, the amount of water present in the aggregates strongly affects the particle-particle interactions and hence the aggregate structure. Detailed analyses of particle-particle and water-particle interactions provide a molecular-scale view of porosity and texture development of the aggregates. The simulation method developed here may also aid in modeling the synthesis of nanostructured materials through self-assembly of nanoparticles.

Molecular Dynamics Simulation of the Crystal Nucleation and Growth Behavior of Methane Hydrate in the Presence of the Surface and Nanopores of Porous Sediment

The behavior of hydrate formation in porous sediment has been widely studied because of its importance in the investigation of reservoirs and in the drilling of natural gas hydrate. However, it is difficult to understand the hydrate nucleation and growth mechanism on the surface and in the nanopores of porous media by experimental and numerical simulation methods. In this work, molecular dynamics simulations of the nucleation and growth of CH4 hydrate in the presence of the surface and nanopores of clay are carried out. The molecular configurations and microstructure properties are analyzed for systems containing one H2O hydrate layer (System A), three H2O hydrate layers (System B), and six H2O hydrate layers (System C) in both clay and the bulk solution. It is found that hydrate formation is more complex in porous media than in the pure bulk solution and that there is cooperativity between hydrate growth and molecular diffusion in clay nanopores. The hydroxylated edge sites of the clay surface could serve as a source of CH4 molecules to facilitate hydrate nucleation. The diffusion velocity of molecules is influenced by the growth of the hydrate that forms a block in the throats of the clay nanopore. Comparing hydrate growth in different clay pore sizes reveals that the pore size plays an important role in hydrate growth and molecular diffusion in clay. This simulation study provides the microscopic mechanism of hydrate nucleation and growth in porous media, which can be favorable for the investigation of the formation of natural gas hydrate in sediments.

Promoting the Adsorption of Metal Ions on Kaolinite by Defect Sites: A Molecular Dynamics Study

Defect sites exist abundantly in minerals and play a crucial role for a variety of important processes. Here molecular dynamics simulations are used to comprehensively investigate the adsorption behaviors, stabilities and mechanisms of metal ions on defective minerals, considering different ionic concentrations, defect sizes and contents. Outer-sphere adsorbed Pb²⁺ ions predominate for all models (regular and defective), while inner-sphere Na⁺ ions, which exist sporadically only at concentrated solutions for regular models, govern the adsorption for all defective models. Adsorption quantities and stabilities of metal ions on kaolinite are fundamentally promoted by defect sites, thus explaining the experimental observations. Defect sites improve the stabilities of both inner- and outer-sphere adsorption, and (quasi) inner-sphere Pb²⁺ ions emerge only at defect sites that reinforce the interactions. Adsorption configurations are greatly altered by defect sites but respond weakly by changing defect sizes or contents. Both adsorption quantities and stabilities are enhanced by increasing defect sizes or contents, while ionic concentrations mainly affect adsorption quantities. We also find that adsorption of metal ions and anions can be promoted by each other and proceeds in a collaborative mechanism. Results thus obtained are beneficial to comprehend related processes for all types of minerals.

Structure of tetraalkylammonium ionic liquids in the interlayer of modified montmorillonite

We perform molecular dynamics simulations of tetraalkylammonium ionic liquids confined in the interlayer of montmorillonite (MMT). We study the structure and energetics of the systems, which consist of cations with two different alkyl chain lengths and several ionic liquid concentrations. The results we obtained for the structure, namely the presence of a strong layering in all systems and the formation of nonpolar domains with interdigitated alkyl chains in some cases, are largely consistent with previous surface force balance experiments performed on similar systems. Finally, we show that swelling of the organo-modified MMT by a large amount of ionic liquid seems energetically favorable in all cases.

Driving force for the hydration of the swelling clays: case of montmorillonites saturated with alkaline-earth cations

Important structural modifications occur in swelling clays upon water adsorption. The multi-scale evolution of the swelling clay structure is usually evidenced by various experimental techniques. However, the driving force behind such phenomena is still not thoroughly understood. It appears strongly dependent on the nature of the interlayer cation. In the case of montmorillonites saturated with alkaline cations, it was inferred that the compensating cation or the layer surface could control the hydration process and thus the opening of the interlayer space, depending on the nature of the interlayer cation. In the present study, emphasis is put on the impact of divalent alkaline-earth cations compensating the layer charge in montmorillonites. Since no experimental technique offers the possibility of directly determining the hydration contributions related to interlayer cations and layer surfaces, an approach based on the combination of electrostatic calculations and immersion data is developed here, as already validated in the case of montmorillonites saturated by alkaline cations. This methodology allows to estimate the hydration energy for divalent interlayer cations and therefore to shed a new light on the driving force for hydration process occurring in montmorillonites saturated with alkaline-earth cations. Firstly, the surface energy values obtained from the electrostatic calculations based on the Electronegativity Equalization Method vary from 450 mJ m(-2) for Mg-montmorillonite to 1100 mJ m(-2) for Ba-montmorillonite. Secondly, considering both the hydration energy for cations and layer surfaces, the driving force for the hydration of alkaline-earth saturated montmorillonites can be attributed to the interlayer cation in the case of Mg-, Ca-, Sr-montmorillonites and to the interlayer surface in the case of Ba-montmorillonites. These results explain the differences in behaviour upon water adsorption as a function of the nature of the interlayer cation, thereby allowing the macroscopic swelling trends to be better understood. The knowledge of hydration processes occurring in homoionic montmorillonites saturated with both the alkaline and the alkaline-earth cations may be of great importance to explain the behaviour of natural clay samples where mixtures of the two types of interlayer cation are present and also provides valuable information on the cation exchange occurring in the swelling clays.