Understanding the structure, distribution, and retention of nanoplastics in montmorillonite nanopore by multi-scale computational simulations

The interfacial adsorption, aggregation and deposition processes of nanoplastics (NPs) on clay mineral surfaces critically regulate their environmental mobility, transformation pathways, and ecotoxicological risks in aquatic ecosystems. A quantitative understanding of the nanoscale interfacial processes is essential. This study employs molecular dynamics (MD) simulations and density functional theory (DFT) calculations to elucidate the aggregation and deposition mechanisms of three types of NPs in their pristine and aged states in the nanopore solution of montmorillonite (Mt). In the wet environment, NPs tend to form aggregates in the nanopore and migrate in solution, increasing environmental risk, while in the dry environment, NPs are more likely to deposit on the basal surface to form larger aggregates, consequently reducing their mobility. Results show hydrophobic interactions play as the primary driving force for the aggregation of pristine NPs, and both hydrophilic and hydrophobic interactions contribute to the aggregation of aged NPs. Aged NPs exhibit stronger binding affinity to Mt through mechanism such as Ca²⁺ bridging and hydrogen bonding, compared to their pristine counterparts. DFT calculations further reveal the formation of hydrogen bonds between the hydroxyl groups of aged NPs and the tetrahedral oxygen atoms in Mt. Through atomic-level characterization of interfacial processes, this work establishes a predictive framework for NP environmental behavior by resolving migration dynamics and retention processes in nanopore water.