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.

Understanding structural anisotropy and mechanical properties of Na-montmorillonite with crystalline swelling and uniaxial deformation under different hydration degrees

A Molecular Dynamics Study on the Structural and Mechanical Properties of Pyrophyllite and M-Montmorillonites (M = Na, K, Ca, and Ba)

A Molecular Dynamics Study on the Structural and Mechanical Properties of Hydrated Kaolinite System under Tension

Deformation and Failure Processes of Kaolinite under Tension:Insights from Molecular Dynamics Simulations

Atomistic study on thermo-mechanical behavior and structural anisotropy of montmorillonite under triaxial tension and compression

Effects of Temperature on Structural Properties of Hydrated Montmorillonite: Experimental Study and Molecular Dynamics Simulation

Investigating Temperature Effects on the Shear Behavior of Clays: Molecular Dynamics Simulations

Molecular Dynamics Simulation of Deformation and Failure Mechanism of Kaolinite at Different Temperatures

A Molecular Dynamics Study of the Mechanical Properties of Kaolinite under Uniaxial and Isothermal Compression at Various Temperatures

Deformation and failure processes of Na-montmorillonite under uniaxial compressive strain condition via molecular dynamic method

Study on Anisotropic Permeability Model for Mixed Kaolin-Montmorillonite Clays

Internal Forces Within the Layered Structure of Na-Montmorillonite Hydrates: Molecular Dynamics Simulation

Atomistic simulations of cation hydration in sodium and calcium montmorillonite nanopores

Molecular Dynamics Study of Crystalline Swelling of Montmorillonite as Affected by Interlayer Cation Hydration

Rheological Behaviors of Na-montmorillonite Considering Particle Interactions: A Molecular Dynamics Study

Insight on molecular interactions in shrinkage of Na-montmorillonite clay by molecular dynamics simulation

Deciphering anomalous dielectric traits of nanoconfined water and its roles in modeling crystalline swelling of bentonite

Evolution of Chemical Bonding and Crystalline Swelling-Shrinkage of Montmorillonite upon Temperature Changes Probed by in Situ Fourier Transform Infrared Spectroscopy and X-ray Diffraction

Nanoscale chemomechanical variations of montmorillonite induced by the specificity of counterions—An in situ XRD and AFM study

Early stage of swelling process of dehydrated montmorillonite through molecular dynamics simulation

Extracting Hydromechanical Properties of Montmorillonite Quasi-Crystals Via Nanoindentation of Oriented Thin Films