Properties and Mechanisms for PFAS Adsorption to Aqueous Clay and Humic Soil Components.

The proliferation of poly- and perfluorinated alkyl substances (PFASs) has resulted in global concerns over contamination and bioaccumulation. PFAS compounds tend to remain in the environment indefinitely, and research is needed to elucidate the ultimate fate of these molecules. We have investigated the model humic substance and model clay surfaces as a potential environmental sink for the adsorption and retention of three representative PFAS molecules with varying chain length and head groups. Utilizing molecular dynamics simulation, we quantify the ability of pyrophyllite and the humic substance to favorably adsorb these PFAS molecules from aqueous solution. We have observed that the hydrophobic nature of the pyrophyllite surface makes the material well suited for the sorption of medium- and long-tail PFAS moieties. Similarly, we find a preference for the formation of a monolayer on the surface for long-chain PFAS molecules at high concentration. Furthermore, we discussed trends in the adsorption mechanisms for the fate and transport of these compounds, as well as potential approaches for their environmental remediation.

A Classical Molecular Dynamics Simulation Study of Interfacial and Bulk Solution Aggregation Properties of Dirhamnolipids.

The rhamnolipids are a unique class of biosurfactants produced by the bacteria Pseudomonas aeruginosa. These molecules display a high level of surface activity as well as biodegradability. In this study nonionic dirhamnolipid was investigated by utilizing molecular dynamics simulation at the air-water interface as well as in bulk water. Detailed structural analysis is presented for both the interfacial simulations and the simulations in solution. A systematic comparison was made between our previous work on the monorhamnolipid at the air-water interface and in bulk water. The presence of a second rhamnose group in dirhamnolipid did not show any significant change in the aggregation at the air-water interface. An increase in the molecular weight resulted in the larger surface area per monomer for dirhamnolipid compared to monorhamnolipid at the air-water interface. However, aggregation of dirhamnolipid in bulk water was affected by the presence of a second rhamnose group. Dirhamnolipid aggregates into micellar structure around ∼N22 which was lower than the monorhamnolipid aggregation number ∼N40. The hydrophobic component of the dirhamnolipid was enhanced to balance the higher hydrophilic component. An increase in alkyl chain length has shown that the enhanced hydrophobic component supports the formation of micellar aggregates up to ∼N30 and above, which was not observed in dirhamnolipid with a shorter alkyl chain length.

Unraveling the Differential Aggregation of Anionic and Nonionic Monorhamnolipids at Air-Water and Oil-Water Interfaces: A Classical Molecular Dynamics Simulation Study.

The molecular structure of a surfactant molecule is known to have a great effect on the interfacial properties. We employ molecular dynamics simulations for a detailed atomistic study of monolayers of the nonionic and anionic form of the most common congener of monorhamnolipids, α-rhamnopyranosyl-ß-hydroxydecanoyl-ß-hydroxydecanoate (( R, R)-Rha-C10-C10), at the air-water and oil-water interfaces. An atomistic-level understanding of monolayer aggregation is necessary to explain a recent experimental observation indicating that nonionic and anionic Rha-C10-C10 show surprisingly different surface area per molecule at the critical micelle concentration. Surface-pressure analysis, interface formation energy calculations, and mass density profiles of the monolayers at the air-water interface show similar properties between nonionic and anionic Rha-C10-C10 aggregation. It is found that there is a significant difference in the headgroup conformations of Rha-C10-C10 in the nonionic and anionic monolayers. Hydrogen bonding interactions between the Rha-C10-C10 molecules in the monolayers is also significantly different between nonionic and anionic forms. Representative snapshots of the simulated system at different surface concentrations show the segregation of molecular aggregates from the interface into the bulk water in the anionic Rha-C10-C10 monolayer at higher concentrations, whereas in the nonionic Rha-C10-C10 monolayer, the molecules are still located at the interface. The present work provides insight into the different aggregation properties of nonionic and anionic Rha-C10-C10 at the air-water interface. Further analyses were carried out to understand the aggregation behavior of nonionic and anionic Rha-C10-C10 at the oil-water interface. It is observed that the presence of oil molecules does not significantly influence the aggregation properties of Rha-C10-C10 as compared to those of the air-water interface.

Molecular Dynamics Simulation of the Oil Sequestration Properties of a Nonionic Rhamnolipid.

A detailed molecular dynamics simulation study is presented on the behavior of aggregates composed of the nonionic monorhamnolipid α-rhamnopyranosyl-ß-hydroxydecanoyl-ß-hydroxydecanoate (Rha-C10-C10) and decane in bulk water. A graph theoretical approach was utilized to characterize the size and composition of the many aggregates generated in our simulations. Overall, we observe that the formation of oil in Rha-C10-C10 aggregates is a favorable process. Detailed analysis on the surfactant/oil aggregate shows that larger aggregates are stable. The shape and size of the aggregates are widely distributed, with the majority of the aggregates preferring ellipsoidal or cylindrical structures. Irrespective of the decane concentration in the system, we did not observe free decane in any of the simulations. Further insights into the binding energy of decane were carried out using free-energy perturbation calculations. The results showed that the trapped decane molecules provide stability to the Rha-C10-C10 aggregates of size N = 50 which are shown to be unstable in our previous study and allow for the growth of larger aggregates than pure Rha-C10-C10 in water. The density profile plots show that decane molecules encapsulated inside the aggregate preferred to remain closer to the center of mass. This study points to the feasibility of using this biosurfactant as an environmental remediation agent.

Structural Properties of Nonionic Monorhamnolipid Aggregates in Water Studied by Classical Molecular Dynamics Simulations.

Molecular dynamics simulations were carried out to investigate the structure and stabilizing factors of aggregates of the nonionic form of the most common congener of monorhamnolipids, α-rhamnopyranosyl-ß-hydroxydecanoyl-ß-hydroxydecanoate (Rha-C10-C10), in water. Aggregates of size ranging from 5 to 810 monomers were observed in the simulation forming spherical and ellipsoidal structures, a torus-like structure, and a unilamellar vesicle. The effects of the hydrophobic chain conformation and alignment in the aggregate, role of monomer···monomer and monomer···water H-bonds, and conformations of monomers in the aggregate were studied in detail. The unilamellar vesicle is highly stable due to the presence of isolated water molecules inside the core adding to the binding energy. Dissociation of a monomer from a larger micellar aggregate is relatively easy compared to that from smaller aggregates as seen from potential of mean force calculations. This analysis also shows that monomers are held more strongly in aggregates of Rha-C10-C10 than the widely used surfactant sodium dodecyl sulfate. Comparisons between the aggregation behavior of nonionic and anionic forms of Rha-C10-C10 are presented.