Abstract:During the last four decades, numerous studies have been directed to the swelling smectite-rich clays in the context of high-level radioactive waste applications and waste-liners for contaminated sites. The swelling properties of clay mineral particles arise due to hydration of the interlayer cations and the diffuse double layers formed near the negatively charged montmorillonite (MMT) surfaces. To accurately study the cation hydration in the interlayer nanopores of MMT, solvent-solute and solvent-clay surface interactions (i.e., the solvation effects and the shape effects) on the atomic level should be taken into account, in contrast to many recent electric double layer based methodologies using continuum models. Therefore, in this research we employed fully atomistic simulations using classical molecular dynamics (MD) simulations, the software package GROMACS along with the CLAYFF forcefield and the SPC/E water model. We present the ion distributions and the deformation of the hydrated coordination structures, i.e., the hydration shells of Na+ and Ca2+ in the interlayer, respectively, for MMT in the first-layer, the second-layer, the third-layer, the fourth-layer, and the fifth-layer (1W, 2W, 3W, 4W, and 5W) hydrate states. Our MD simulations show that Na+ in Na-MMT nanopores have an affinity to the ditrigonal cavities of the clay layers and form transient inner-sphere complexes at about 3.8 Å from clay midplane at water contents less than the 5W hydration state. However, these phenomena are not observed in Ca-MMT regardless of swelling states. For Na-MMT, each Na+ is coordinated to four water molecules and one oxygen atom of the clay basal-plane in the first hydration shell at the 1W hydration state, and with five to six water molecules in the first hydration shell within a radius of 3.1 Å at all higher water contents. In Ca-MMT, however each Ca2+ is coordinated to approximately seven water molecules in the first hydration shell at the 1W hydration state and about eight water molecules in the first hydration shell within a radius of 3.3 Å at all higher hydration states. Moreover, the MD results show that the complete hydration shells are nearly spherical with an orthogonal coordination sphere. They could only be formed when the basal spacing d001 ≥ 18.7 Å, i.e., approximately, the interlayer separation h ≥ 10 Å. Comparison between DFT and MD simulations shows that DFT failed to reproduce the outer-sphere complexes in the Stern-layer (within ∼5.0 Å from the clay basal-plane), observed in the MD simulations.

Modelling aqueous solubility of sodium chloride in clays at thermodynamic conditions of hydraulic fracturing by molecular simulations

Concentrated aqueous sodium chloride solution in clays at thermodynamic conditions of hydraulic fracturing: Insight from molecular dynamics simulations

Numerical simulation of formation water salinity redistribution in fractured shale reservoirs during hydraulic fracturing

Molecular Dynamic Simulations of Montmorillonite–Organic Interactions under Varying Salinity: An Insight into Enhanced Oil Recovery

Molecular Simulation of Carbon Dioxide, Brine, and Clay Mineral Interactions and Determination of Contact Angles

Pore‐Scale Modeling of Water and Ion Diffusion in Partially Saturated Clays

The solubility of H 2 in NaCl brine at high pressures and high temperatures: Molecular simulation study and thermodynamic modeling

The flowback behavior of salt in hydraulically fractured shale under multi-phase flow conditions: Modelling, simulation and application

Effects of Salt Ions on the Methane Hydrate Formation and Dissociation in the Clay Pore Water and Bulk Water

Mathematical Modeling and Numerical Simulation of Water‐Rock Interaction in Shale under Fracturing‐fluid Flowback Conditions

Molecular simulation of adsorption behaviors of methane and carbon dioxide on typical clay minerals

Molecule Simulation of CH4/CO2 Competitive Adsorption and CO2 Storage in Shale Montmorillonite

Molecular simulation on the effect of formation depth on methane adsorption by clay minerals

Investigation of Clay Type on Low Salinity Water Flooding Using a Glass Micromodel

Atomistic simulations of cation hydration in sodium and calcium montmorillonite nanopores

Modelling of the Phase-Partitioning Behaviors for CO2-Brine System at Geological Conditions

Coupling between clay swelling/collapse and cationic partition

Molecular Dynamics Simulations of Carbon Dioxide Intercalation in Hydrated Na-Montmorillonite

Simulations of vapor–liquid phase equilibrium and interfacial tension in the CO2–H2O–NaCl system

Multiphase Monte Carlo and Molecular Dynamics Simulations of Water and CO2 Intercalation in Montmorillonite and Beidellite

Distribution and mobility of crude oil-brine in clay mesopores: Insights from molecular dynamics simulations.