A Coarse-Grained Model for the Mechanical Behavior of Na-Montmorillonite Clay

Sodium montmorillonite (Na-MMT) is one of the most commonly found swelling clay minerals with diverse engineering and technological applications. The nanomechanical properties of this mineral have been extensively investigated computationally utilizing molecular dynamics (MD) simulations to portray the molecular-level changes at different environmental conditions. As the environmentally found Na-MMT clays are generally sized within hundreds of nanometers, all-atomistic (AA) MD simulations of clays within such size range are particularly challenging due to computational inefficiency. Informed from atomistic modeling, a coarse-grained (CG) modeling technique can be employed to overcome the spatiotemporal limitation. The current study presents a modeling strategy to develop a computationally efficient model of Na-MMT clay with a typical size over ≃100 nm by shrinking the atomistic platelet thickness and reducing the number of center-layer atoms. Using the "strain-energy conservation" approach, the force field parameters for the CG model are obtained and the developed CG model can well preserve in-plane tension, shear, and bending behaviors of atomistic counterparts. Remarkably, the CG tactoid model of Na-MMT, a hierarchical multilayer structure, can reproduce the interlayer shear and adhesion as well as d-spacing among the clay sheets as of atomistic one to a good approximation while gaining significantly improved computational speed. Our study demonstrates the efficacy of the CG modeling framework, paving the way for the bottom-up multiscale prediction of mechanical behaviors of clay and related minerals.