Molecular Dynamics Simulation Prediction of the Partitioning Constants (KH, Kiw, Kia) of 82 Legacy and Emerging Organic Contaminants at the Water-Air Interface.

The tendency of organic contaminants (OCs) to partition between different phases is a key set of properties that underlie their human and ecological health impacts and the success of remediation efforts. A significant challenge associated with these efforts is the need for accurate partitioning data for an ever-expanding list of OCs and breakdown products. All-atom molecular dynamics (MD) simulations have the potential to help generate these data, but existing studies have applied these techniques only to a limited variety of OCs. Here, we use established MD simulation approaches to examine the partitioning of 82 OCs, including many compounds of critical concern, at the water-air interface. Our predictions of the Henry's law constant (KH) and interfacial adsorption coefficients (Kiw, Kia) correlate strongly with experimental results, indicating that MD simulations can be used to predict KH, Kiw, and Kia values with mean absolute deviations of 1.1, 0.3, and 0.3 logarithmic units after correcting for systematic bias, respectively. A library of MD simulation input files for the examined OCs is provided to facilitate future investigations of the partitioning of these compounds in the presence of other phases.

Molecular dynamics simulation of the adsorption of per- and polyfluoroalkyl substances (PFASs) on smectite clay.

Molecular dynamics (MD) simulations are used to predict the partitioning of per- and polyfluoroalkyl substances (PFASs) to smectite clay, a high surface area adsorbent ubiquitous in temperate soils. Simulated systems model a stack of flexible smectite lamellae in contact with a bulk aqueous reservoir containing PFAS molecules. Perfluorobutanesulfonic acid (PFBS), perfluorohexanesulfonic acid (PFHxS), and perfluorooctanesulfonic acid (PFOS) are simulated at various aqueous chemistry conditions to examine the effect of PFAS size, salinity, and coordinating cation type (K+, Na+, and Ca2+) on adsorption. The metadynamics technique is employed to facilitate the exploration of the simulation cell and to reconstruct the underlying free energy landscape. Adsorption is favorable on the hydrophobic domains of the external basal surfaces with the fluorinated chain adopting a flat orientation on the surface. Analysis of the adsorption energetics reveals large favorable entropic contributions to adsorption. The enthalpy of adsorption is unfavorable, though much less so in the presence of Ca2+ due to stabilizing 'lateral cation bridging' interactions between divalent cations and PFAS sulfonate head groups. Overall, this research advances the mechanistic understanding of PFAS-smectite interactions and provides new insights that could help inform fate and transport models and the development of adsorbents and remediation techniques.