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.

Hydration Characteristics of Ca2+ and Mg2+: A Density Functional Theory, Polarized Continuum Model and Molecular Dynamics Investigation

Second Hydration Shell of Mg 2+ : Competition between Ion-Water Interaction and Hydrogen Bonding Interaction

A density functional theory study of the hydration of calcium ions confined in the interlayer space of montmorillonites

Hydration numbers of biologically relevant divalent metal cations from ab initio molecular dynamics and continuum solvation methods

Molecular Dynamics Study of the Structure and Dynamics of the Hydration Shell of Alkaline and Alkaline-Earth Metal Cations

Molecular Dynamics Simulation of the Hydrated Magnesium Ion

Molecular Dynamics Simulation, Ab Initio Calculation, and Size-Selected Anion Photoelectron Spectroscopy Study of Initial Hydration Processes of Calcium Chloride

Hydration Energies of Divalent Beryllium and Magnesium Ions: An ab Initio Molecular Orbital Study

Hydration properties of alkali and alkaline earth metal ions in aqueous solution: A molecular dynamics study

Molecular dynamics study on ionic hydration

Structures, Stability, Vibration Entropy and Ir Spectra of Hydrated Calcium Ion Clusters [ca(h2o)(n)](2+) (N=1-20, 27): A Systematic Investigation by Density Functional Theory

Molecular Simulation of the Ions Ultraconfined in the Nanometer-Channel of Calcium Silicate Hydrate: Hydration Mechanism, Dynamic Properties, and Influence on the Cohesive Strength

Hydration shell structure of the calcium ion from simulations with ab initio effective pair potentials

Ab Initio Molecular Dynamics Study of Aqueous Solutions of Magnesium and Calcium Nitrates: Hydration Shell Structure, Dynamics and Vibrational Echo Spectroscopy.

Molecular dynamics simulations study for the effect of cations hydration on the surface tension of the electrolyte solutions

Ion Hydration and Associated Defects in Hydrogen Bond Network of Water: Observation of Reorientationally Slow Water Molecules Beyond First Hydration Shell in Aqueous Solutions of MgCl$_2$

Atomistic simulations of cation hydration in sodium and calcium montmorillonite nanopores

Hydration and Mobility of Interlayer Ions of (Na x , Ca y )-Montmorillonite: A Molecular Dynamics Study

Ground-State Structures of Hydrated Calcium Ion Clusters From Comprehensive Genetic Algorithm Search

Solvation Structure of Calcium and Magnesium Ions in Water with the Presence of Hydroxide: a Study by Deep Potential Molecular Dynamics