Abstract:The edges of montmorillonite (MMT) react strongly with metals and organic matter, but the atomic structure of the edge and its surface complexes are not unambiguous since the experimental isolation of the edge is challenging. In this study, we introduce an atomistic model of a Na MMT edge that is suitable for classical molecular dynamics (MD) simulations, in particular for the B edge, a representative edge surface of 2:1 phyllosilicates. Our model possesses the surface groups identified through density functional theory (DFT) geometry optimizations performed with variation in the structural charge deficit and Mg substitution sites. The edge structure of the classical MD simulations agreed well with previous DFT-based MD simulation results. Our MD simulations revealed an extensive H-bond network stabilizing the Na-MMT edge surface, which required an extensive simulation trajectory. Some Na counter ions formed inner-sphere complexes at two edge sites. The stronger edge site coincided with the exposed vacancy in the dioctahedral sheet; a weaker site was associated with the cleaved hexagonal cavity of the tetrahedral sheet. The six-coordinate Na complexes were not directly associated with the Mg edge site. Our simulations have demonstrated the heterogeneous surface structures, the distribution of edge surface groups, and the reactivity of the MMT edge.