Effect of interlayer cations on montmorillonite swelling: Comparison between molecular dynamic simulations and experiments

The understanding of the swelling phenomenon of montmorillonite is essential to predict the physical and chemical behavior of clay-based barriers in geoenvironmental engineering. This study investigated the crystalline swelling behavior of montmorillonite with different interlayer cations including monovalent Na, K, and Cs and divalent Ca and Sr by molecular dynamics simulations and experimental measurements including X-ray diffraction and water vapor adsorption. The stepwise swelling behavior for Na- and Ca/Sr-montmorillonite demonstrated good agreement between molecular dynamics (MD) simulated and experimental results. The swelling curves of K- and Cs-montmorillonite remaining in the 1-hydrated state even at high relative humidity showed systematic discrepancies between MD simulations and experimental results. The comparative discussion between the experimental and MD derived swelling curves offer insights that may help to understand the swelling mechanism of montmorillonite and the dependence on interlayer cations. The water content in the interlayer and external surface in five homoionic montmorillonite was strongly correlated to the hydration number of each ion. On the basis of these comparisons, the hydration free energy of interlayer cations seems to contribute greatly to the total driving force for the crystalline swelling of montmorillonite. An additional key factor is the preference of an outer- or inner-sphere complex of interlayer cations and the distribution of cations in the interlayer space. This effect has a considerable impact in the case of monovalent cations and results in different swelling behavior between the outer-sphere Na and inner-sphere K and Cs ions.