Molecular Dynamics Simulations of the Adsorption of Phthalate Esters on Smectite Clay Surfaces

The partitioning of organic contaminants between water and solid surfaces is a key process controlling their fate and transport in natural environments. A novel methodology was developed to predict the adsorption of organic contaminants by smectite clay minerals (high specific surface area adsorbents abundant in natural soils) using molecular dynamics (MD) simulations. The methodology models a stack of flexible Ca-montmorillonite lamellae in direct contact with a bulk aqueous reservoir and uses the metadynamics technique to facilitate the exploration of the free energy landscape. The methodology was tested and validated in the case of six phthalate esters, widely used chemical plasticizers with endocrine disrupting properties. Simulation predictions reveal strong phthalate adsorption, especially for the larger and more hydrophobic phthalates. Predicted partition coefficients (Kd values) are consistent with collected batch experimental data. Adsorption was observed on both the exterior basal surfaces and within the interlayer nanopore, with phthalate molecules predominately adopting a flat orientation on the clay surface. Intercalation was also detected in complementary X-ray diffraction (XRD) experiments. A strong inverse relationship between extent of adsorption and clay surface charge density was observed, as phthalate molecules preferentially occupied the more hydrophobic uncharged patches on each surface. Detailed analysis of the free energy of adsorption revealed that phthalate affinity for the clay surface results from a small favorable van der Waals contribution and a large favorable entropic contribution. Overall, this research demonstrates the substantial affinity of smectite clays for phthalate esters and establishes a computational methodology capable of predicting the water–clay partition coefficients of organic contaminants, a key parameter in environmental fate and transport models.