Effect of layer charge density and cation concentration on sorption behaviors of heavy metal ions in the interlayer and nanopore of montmorillonite: A molecular dynamics simulation

Soil and ground water pollution by heavy metal is an identified global environment concern. Montmorillonite is widely adopted as an adsorbent used in retardation of heavy metals due to its excellent adsorption capacity, conveniency, and stability. However, the adsorption mechanism of heavy metal ions in the interlayer and nanopore of montmorillonite needs to be further studied. The classical grand canonical monte carlo (GCMC) and molecular dynamics (MD) simulations were performed to quantify the structural and dynamic properties of heavy metal ions in the interlayer and nanopore of montmorillonite (MMT). Two layer charge densities of MMT, two cation concentrations, and three typical heavy metal ions (Pb2+, Ba2+, and Cs+) were considered. With the increase of layer charge density, outer-sphere adsorption (OSA) was transformed into inner-sphere adsorption (ISA) for three heavy metal ions. The adsorption quantity increased with the increase of layer charge density and cation concentration. The calculated diffusion coefficients of three heavy metal caions in montmorillonite was bulk solution > nanopore > interlayer. The results and adsorption mechanism are helpful for understanding the natural processes of heavy metal ions in environment and developing more efficient remediation materials.