A Quantum Chemistry Based Force Field for Perfluoroalkanes and Poly(tetrafluoroethylene)

In silico prediction of the interaction of legacy and novel per- and poly-fluoroalkyl substances (PFAS) with selected human transporters and of their possible accumulation in the human body

Per- and poly-fluorinated compounds constitute a wide group of fluorocarbon chemicals with widespread industrial applications, ranging from non-stick coating in cookware to water surfactants, from fire-fighting foams to water-repellent coatings on textiles. Presently, over 12,000 PFAS are known worldwide. In recent years, extensive research has focused on investigating the biological effects of these molecules on various organisms, including humans. Here, we conducted in silico simulations to examine the potential binding of a representative selection of PFAS to various human proteins known to be involved in chemical transportation and accumulation processes. Specifically, we targeted human serum albumin (HSA), transthyretin (TTR), thyroxine binding protein (TBG), fatty acid binding proteins (FABPs), organic anion transporters (OATs), aiming to assess the potential for bioaccumulation. Molecular docking simulations were employed for this purpose, supplemented by molecular dynamics (MD) simulations to account for protein flexibility, when necessary. Our findings indicate that so-called "legacy PFAS" such as PFOA or PFOS exhibit a higher propensity for interaction with the analysed human protein targets compared to newly formulated PFAS, characterised by higher branching and hydrophilicity, and possibly a higher accumulation in the human body.

Insights into PFAS environmental fate through computational chemistry: A review

Per- and polyfluoroalkyl substances (PFAS) are widely used chemicals that exhibit exceptional chemical and thermal stability. However, their resistance to degradation has led to their widespread environmental contamination. PFAS also negatively affect the environment and other organisms, highlighting the need for effective remediation methods to mitigate their presence and prevent further contamination. Computational chemistry methods, such as Density Functional Theory (DFT) and Molecular Dynamics (MD) offer valuable tools for studying PFAS and simulating their interactions with other molecules. This review explores how computational chemistry methods contribute to understanding and tackling PFAS in the environment. PFAS have been extensively studied using DFT and MD, each method offering unique advantages and computational limitations. MD simulates large macromolecules systems however it lacks the ability model chemical reactions, while DFT provides molecular insights however at a high computational cost. The integration of DFT with MD shows promise in predicting PFAS behavior in different environments. This work summarizes reported studies on PFAS compounds, focusing on adsorption, destruction, and bioaccumulation, highlighting contributions of computational methods while discussing the need for continued research. The findings emphasize the importance of computational chemistry in addressing PFAS contamination, guiding risk assessments, and informing future research and innovations in this field.

Phase transitions of fluorotelomer alcohols at the water|alkane interface studied via molecular dynamics simulation

Fluorosurfactants are long-lasting environmental pollutants that accumulate at interfaces ranging from aerosol droplet surfaces to cell membranes. Modeling of adsorption-based removal technologies for fluorosurfactants requires accurate simulation methods which can predict their adsorption isotherm and monolayer structure. Fluorotelomer alcohols with one or two methylene groups adjacent to the alcohol (7 : 1 FTOH and 6 : 2 FTOH, respectively) are investigated using the OPLS-AA force field at the water|hexane interface, varying the interfacial area per surfactant. The acquired interfacial pressure isotherms and monolayer phase behavior are compared with previous experimental results. The results are consistent with the experimental data inasmuch as, at realistic adsorption densities, only 7 : 1 FTOH shows a phase transition between liquid-expanded (LE) and 2D crystalline phases. Structures of the LE and crystalline phases are in good agreement with the sticky disc and Langmuir defective crystal models, respectively, used previously to interpret experimental data. Interfacial pressure of the LE phase agrees well with experiment, and sticky disc interaction parameters indicate no 2D LE-gas transition is present for either molecule. Conformation analysis reveals 7 : 1 FTOH favors conformers where the OH dipole is perpendicular to the molecular backbone, such that the crystalline phase is stabilized when these dipoles align.

A ReaxFF-based molecular dynamics study of the destruction of PFAS due to ultrasound

This work uses a computational approach to provide a mechanistic explanation for the experimentally observed destruction of per- and polyfluoroalkyl substances (PFAS) in water due to ultrasound. The PFAS compounds have caused a strong public and regulatory response due to their ubiquitous presence in the environment and toxicity to humans. In this research, ReaxFF -based Molecular Dynamics simulation under several temperatures ranging from 373 K to 5,000 K and different environments such as water vapor, O2, N2, and air were performed to understand the mechanism of PFAS destruction. The simulation results showed greater than 98% PFAS degradation was observed within 8 ns under a temperature of 5,000 K in a water vapor phase, replicating the observed micro/nano bubbles implosion and PFAS destruction during the application of ultrasound. Additionally, the manuscript discusses the reaction pathways and how PFAS degradation evolves providing a mechanistic basis for the destruction of PFAS in water due to ultrasound. The simulation showed that small chain molecules C1 and C2 fluoro-radical products are the most dominant species over the simulated period and are the impediment to an efficient degradation of PFAS. Furthermore, this research confirms the empirical findings observations that the mineralization of PFAS molecules occurs without the generation of byproducts. These findings highlight the potential of virtual experiments in complementing laboratory experiments and theoretical projections to enhance the understanding of PFAS mineralization during the application of ultrasound.

Controlling the self-assembly of perfluorinated surfactants in aqueous environments

Surface active per- and polyfluoroalkyl substances (PFAS) released in the environment generate great concern in the US and worldwide. The sequestration of PFAS amphiphiles from aqueous media can be limited by their strong tendency to form micelles that plug the pores in the adsorbent material, rendering most of the active surface inaccessible. A joint experimental and simulation approach has been used to investigate the structure of perfluorooctanoate ammonium (PFOA) micelles in aqueous solutions, focusing on the understanding of ethanol addition on PFOA micelle formation and structure. Structurally compact and slightly ellipsoidal in shape, PFOA micelles in pure water become more diffuse with increasing ethanol content, and break into smaller PFOA clusters in aqueous solutions with high ethanol concentration. A transition from a co-surfactant to a co-solvent behavior with the increase of ethanol concentration has been observed by both experiments and simulations, while the latter also provide insight on how to achieve co-solvent conditions with other additives. An improved understanding of how to modulate PFAS surfactant self-assembly in water can inform the fate and transport of PFAS in the environment and the PFAS sequestration from aqueous media.

Role of Chemical Dynamics Simulations in Mass Spectrometry Studies of Collision-Induced Dissociation and Collisions of Biological Ions with Organic Surfaces

In this article, a perspective is given of chemical dynamics simulations of collisions of biological ions with surfaces and of collision-induced dissociation (CID) of ions. The simulations provide an atomic-level understanding of the collisions and, overall, are in quite good agreement with experiment. An integral component of ion/surface collisions is energy transfer to the internal degrees of freedom of both the ion and the surface. The simulations reveal how this energy transfer depends on the collision energy, incident angle, biological ion, and surface. With energy transfer to the ion's vibration fragmentation may occur, i.e. surface-induced dissociation (SID), and the simulations discovered a new fragmentation mechanism, called shattering, for which the ion fragments as it collides with the surface. The simulations also provide insight into the atomistic dynamics of soft-landing and reactive-landing of ions on surfaces. The CID simulations compared activation by multiple "soft" collisions, resulting in random excitation, versus high energy single collisions and nonrandom excitation. These two activation methods may result in different fragment ions. Simulations provide fragmentation products in agreement with experiments and, hence, can provide additional information regarding the reaction mechanisms taking place in experiment. Such studies paved the way on using simulations as an independent and predictive tool in increasing fundamental understanding of CID and related processes.

Molecular Dynamics Simulation of the Superspreading of Surfactant-Laden Droplets. A Review

Superspreading is the rapid and complete spreading of surfactant-laden droplets on hydrophobic substrates. This phenomenon has been studied for many decades by experiment, theory, and simulation, but it has been only recently that molecular-level simulation has provided significant insights into the underlying mechanisms of superspreading thanks to the development of accurate force-fields and the increase of computational capabilities. Here, we review the main advances in this area that have surfaced from Molecular Dynamics simulation of all-atom and coarse-grained models highlighting and contrasting the main results and discussing various elements of the proposed mechanisms for superspreading. We anticipate that this review will stimulate further research on the interpretation of experimental results and the design of surfactants for applications requiring efficient spreading, such as coating technology.

Mechanistic details of energy transfer and soft landing in ala2-H(+) collisions with a F-SAM surface

Previous chemical dynamics simulations (Phys. Chem. Chem. Phys., 2014, 16, 23769-23778) were analyzed to delineate atomistic details for collision of N-protonated dialanine (ala2-H(+)) with a C8 perfluorinated self-assembled monolayer (F-SAM) surface. Initial collision energies Ei of 5-70 eV and incident angles θi of 0° and 45°, with the surface normal, were considered. Four trajectory types were identified: (1) direct scattering; (2) temporary sticking/physisorption on top of the surface; (3) temporary penetration of the surface with additional physisorption on the surface; and (4) trapping on/in the surface, by physisorption or surface penetration, when the trajectory is terminated. Direct scattering increases from 12 to 100% as Ei is increased from 5 to 70 eV. For the direct scattering at 70 eV, at least one ala2-H(+) heavy atom penetrated the surface for all of the trajectories. For ∼33% of the trajectories all eleven of the ala2-H(+) heavy atoms penetrated the F-SAM at the time of deepest penetration. The importance of trapping decreased with increase in Ei, decreasing from 84 to 0% with Ei increase from 5 to 70 eV at θi = 0°. Somewhat surprisingly, the collisional energy transfers to the F-SAM surface and ala2-H(+) are overall insensitive to the trajectory type. The energy transfer to ala2-H(+) is primarily to vibration, with the transfer to rotation ∼10% or less. Adsorption and then trapping of ala2-H(+) is primarily a multi-step process, and the following five trapping mechanisms were identified: (i) physisorption-penetration-physisorption (phys-pen-phys); (ii) penetration-physisorption-penetration (pen-phys-pen); (iii) penetration-physisorption (pen-phys); (iv) physisorption-penetration (phys-pen); and (v) only physisorption (phys). For Ei = 5 eV, the pen-phys-pen, pen-phys, phys-pen, and phys trapping mechanisms have similar probabilities. For 13.5 eV, the phys-pen mechanism, important at 5 eV, is unimportant. The radius of gyration of ala2-H(+) was calculated once it is trapped on/in the F-SAM surface and trapping decreases the ion's compactness, in part by breaking hydrogen bonds. The ala2-H(+) + F-SAM simulations are compared with the penetration and trapping dynamics found in previous simulations of projectile + organic surface collisions.

PFC and Triglyme for Li-Air Batteries: A Molecular Dynamics Study

In this work, we present an all-atom molecular dynamics (MD) study of triglyme and perfluorinated carbons (PFCs) using classical atomistic force fields. Triglyme is a typical solvent used in nonaqueous Li-air battery cells. PFCs were recently reported to increase oxygen availability in such cells. We show that O2 diffusion in two specific PFC molecules (C6F14 and C8F18) is significantly faster than in triglyme. Furthermore, by starting with two very different initial configurations for our MD simulation, we demonstrate that C8F18 and triglyme do not mix. The mutual solubility of these molecules is evaluated both theoretically and experimentally, and a qualitative agreement is found. Finally, we show that the solubility of O2 in C8F18 is considerably higher than in triglyme. The significance of these results to Li-air batteries is discussed.

Self-Consistent Determination of Atomic Charges of Ionic Liquid through a Combination of Molecular Dynamics Simulation and Density Functional Theory

A self-consistent scheme is developed to determine the atomic partial charges of ionic liquid. Molecular dynamics (MD) simulation was conducted to sample a set of ion configurations, and these configurations were subject to density functional theory (DFT) calculations to determine the partial charges. The charges were then averaged and used as inputs for the subsequent MD simulation, and MD and DFT calculations were repeated until the MD results are not altered any more. We applied this scheme to 1,3-dimethylimidazolium bis(trifluoromethylsulfonyl) imide ([C1mim][NTf2]) and investigated its structure and dynamics as a function of temperature. At convergence, the average ionic charges were ±0.84 e at 350 K due to charge transfer among ions, where e is the elementary charge, while the reduced ionic charges do not affect strongly the density of [C1mim][NTf2] and radial distribution function. Instead, major effects are found on the energetics and dynamics, with improvements of the overestimated heat of vaporization and the too slow motions of ions observed in MD simulations using commonly used force fields.

Chemical dynamics simulations of energy transfer, surface-induced dissociation, soft-landing, and reactive-landing in collisions of protonated peptide ions with organic surfaces

Different simulation approaches like MM, QM + MM, and QM/MM, were used to study surface-induced dissociation, soft-landing, and reactive-landing for the peptide-H⁺ + surface collisions.

Smoothing of contact lines in spreading droplets by trisiloxane surfactants and its relevance for superspreading

Superspreading, the greatly enhanced spreading of aqueous solutions of trisiloxane surfactants on hydrophobic substrates, is of great interest in fundamental physics and technical applications. Despite numerous studies in the last 20 years, the superspreading mechanism is still not well understood, largely because the molecular scale cannot be resolved appropriately either experimentally or using continuum simulations. The absence of molecular-scale knowledge has led to a series of conflicting hypotheses based on different assumptions of surfactant behavior. We report a series of large-scale molecular dynamics simulations of aqueous solutions of superspreading and non-superspreading surfactants on different substrates. We find that the transition from the liquid-vapor to the solid-liquid interface is smooth for superspreading conditions, allowing direct adsorption through the contact line. This finding complements a study [Karapetsas et al., J. Fluid Mech., 2011, 670, 5-37], which predicts that superspreading can occur if this adsorption path is possible. Based on the observed mechanism, we provide plausible explanations for the influence of the substrate hydrophobicity, the surfactant chain length, and the surfactant concentration on the superspreading phenomenon. We also briefly address that the observed droplet shape is a mechanism to overcome the Huh-Scriven paradox of infinite viscous dissipation at the contact line.

Perfluoroalkane force field for lipid membrane environments

In this work, we present atomic parameters of perfluoroalkanes for use within the CHARMM force field. Perfluorinated alkanes represent a special class of molecules. On the one hand, they are considerably more hydrophobic than lipids, but on the other hand, they are not lipophilic either. Instead, they represent an independent class of philicity, enabling a whole portfolio of applications within both materials science and biochemistry. We performed a thorough parametrization of all bonded and nonbonded parameters with a particular focus on van der Waals parameters. Here, the general framework of the CHARMM and CGenFF force fields has been followed. The van der Waals parameters have been fitted to experimental densities over a wide range of temperatures and pressures. This newly parametrized class of molecules will open the gate for a variety of simulations of biologically relevant systems within the CHARMM force field. A particular perspective for the present work is the influence of polyphilic transmembrane molecules on membrane properties, aggregation phenomena, and transmembrane channels.

Dynamics of energy transfer and soft-landing in collisions of protonated dialanine with perfluorinated self-assembled monolayer surfaces

Chemical dynamics simulations are reported which provide atomistic details of collisions of protonated dialanine, ala2-H(+), with a perfluorinated octanethiolate self-assembled monolayer (F-SAM) surface. The simulations are performed at collision energies Ei of 5.0, 13.5, 22.5, 30.00, and 70 eV, and incident angles 0° (normal) and 45° (grazing). Excellent agreement with experiment (J. Am. Chem. Soc., 2000, 122, 9703-9714) is found for both the average fraction and distribution of the collision energy transferred to the ala2-H(+) internal degrees of freedom. The dominant pathway for this energy transfer is to ala2-H(+) vibration, but for Ei = 5.0 eV ∼20% of the energy transfer is to ala2-H(+) rotation. Energy transfer to ala2-H(+) rotation decreases with increase in Ei and becomes negligible at high Ei. Three types of collisions are observed in the simulations: i.e. those for which ala2-H(+) (1) directly scatters off the F-SAM surface; (2) sticks/physisorbs on/in the surface, but desorbs within the 10 ps numerical integration of the simulations; and (3) remains trapped (i.e. soft-landed) on/in the surface when the simulations are terminated. Penetration of the F-SAM by ala2-H(+) is important for the latter two types of events. The trapped trajectories are expected to have relatively long residence times on the surface, since a previous molecular dynamics simulation (J. Phys. Chem. B, 2014, 118, 5577-5588) shows that thermally accommodated ala2-H(+) ions have an binding energy with the F-SAM surface of at least ∼15 kcal mol(-1).

The elastic modulus of isolated polytetrafluoroethylene filaments

We report vibrational Raman spectra of small extended perfluoro-n-alkanes (C{n}F{2n+2} with n=6, 8-10, 12-14) isolated in supersonic jet expansions and use wavenumbers of longitudinal acoustic vibrations to extrapolate the elastic modulus of cold, isolated polytetrafluoroethylene filaments. The derived value E = 209(10) GPa defines an upper limit for the elastic modulus of the perfectly crystalline, non-interacting polymer at low temperatures and serves as a benchmark for quantum chemical predictions.

Atomistic potentials for trisiloxane, alkyl ethoxylate, and perfluoroalkane-based surfactants with TIP4P/2005 and application to simulations at the air-water interface

The mechanism of superspreading, the greatly enhanced spreading of water droplets facilitated by trisiloxane surfactants, is still under debate, largely because the role and behavior of the surfactants cannot be sufficiently resolved by experiments or continuum simulations. Previous molecular dynamics studies have been performed with simple model molecules or inaccurate models, strongly limiting their explanatory power. Here we present a force field dedicated to superspreading, extending existing quantum-chemistry-based models for the surfactant and the TIP4P/2005 water model ( Abascal et al. J. Chem. Phys. , 2005 , 123 , 234505 ). We apply the model to superspreading trisiloxane surfactants and nonsuperspreading alkyl ethoxylate and perfluoroalkane surfactants at various concentrations at the air-water interface. We show that the developed model accurately predicts surface tensions, which are typically assumed important for superspreading. Significant differences between superspreading and traditional surfactants are presented and their possible relation to superspreading discussed. Although the force field has been developed for superspreading problems, it should also perform well for other simulations involving polymers or copolymers with water.

The elastic modulus of isolated polytetrafluoroethylene filaments

We report vibrational Raman spectra of small extended perfluoro-n-alkanes (C{n}F{2n+2} with n=6, 8-10, 12-14) isolated in supersonic jet expansions and use wavenumbers of longitudinal acoustic vibrations to extrapolate the elastic modulus of cold, isolated polytetrafluoroethylene filaments. The derived value E = 209(10) GPa defines an upper limit for the elastic modulus of the perfectly crystalline, non-interacting polymer at low temperatures and serves as a benchmark for quantum chemical predictions.

Hydrogenated/fluorinated catanionic surfactants as potential templates for nanostructure design

The structure and physicochemical properties of the nanoparticles spontaneously formed within aqueous mixtures of the hydrogenated/fluorinated catanionic surfactant cetyltrimetylammonium perfluorooctanoate in the absence of counterions as a function of its concentration are investigated by a combined experimental/computational study at room temperature. Apparent molar volumes, isentropic apparent molar compressibilities, and dynamic light scattering measurements together with transmission and cryo-scanning electron as well as confocal laser microscopy images, and computational molecular dynamics simulations indicate that a variety of structures of different sizes coexist in solution with vesicles of ∼160 nm diameter. Interestingly, the obtained nanostructures were observed to self-assemble from a random distribution of monomers in a time scale easily accessible by atomistic classical molecular dynamics simulations, allowing to provide a comprehensive structural and dynamic characterization of the surfactant molecules at atomic level within the different aggregates. Overall, it is demonstrated that the use of mixed fluorinated hydrogenated surfactant systems represents an easy strategy for the design of specific nanoscale structures. The detailed structural analysis provided in the present work is expected to be useful as a reference to guide the design of new nanoparticles based on different hydrogenated/fluorinated catanionic surfactants.

Langmuir monolayers of a hydrogenated/fluorinated catanionic surfactant: from the macroscopic to the nanoscopic size scale

Langmuir monolayers of the hydrogenated/fluorinated catanionic surfactant cetyltrimethylammonium perfluorooctanoate at the air/water interface are studied at room temperature. Excess Gibbs energies of mixing, DeltaG(E), as well as transition areas and pressures, were obtained from the surface pressure-area isotherm. The DeltaG(E) curve indicates that tail-tail interactions are more important than head-head interactions at low pressures and vice versa. Atomic force microscopy and molecular dynamics simulations allowed a fine characterization of the monolayer structure as a function of the area per molecule at mesoscopic and nanoscopic size scales, respectively. A combined analysis of the techniques allow us to conclude that electrostatic interactions between the ionic head groups are dominant in the monolayer while hydrophobic parts are of secondary importance. Overall, results obtained from the different techniques complement to each other, giving a comprehensive characterization of the monolayer.

Experimental pKa determination for perfluorooctanoic acid (PFOA) and the potential impact of pKa concentration dependence on laboratory-measured partitioning phenomena and environmental modeling

An accurately measured equilibrium acid dissociation constant (pKa) is essential for understanding and predicting the fate of perfluorocarboxylic acids (PFCAs) in the environment. The aqueous pKa of perfluorooctanoic acid (PFOA) has been determined potentiometrically using a standard water-methanol mixed solvent approach and was found to be 3.8 +/- 0.1. The acidity of PFOA is thus considerably weaker than its shorter-chain PFCA homologues. This was attributed to differences in molecular and electronic structure, coupled with solvation effects. The pKa of PFOA was suppressed to approximately 2.3 at higher concentrations because of the aggregation of perfluorooctanoate (PFO). Often, PFCA partion coefficients are determined at concentrations above those found in the environment. Thus, it was suggested that a pKa correction factor, which accounts for this concentration-dependent shift in acid/base equilibrium, should be applied to PFCA partition efficients before they are implemented in environmental fate models. A pKa of 3.8 +/- 0.1 suggests that a considerable concentration of the PFCA exists as the neutral species in the aqueous environment for example, in typical Ontario rainwater, it is approximately 17%. Transport, fate, and partitioning models have often ignored the presence this species completely. The environmental dissemination of PFCAs could, in part, be explained by considering the role of the neutral species.

Structural and phase properties of tetracosane (C24H50) monolayers adsorbed on graphite: an explicit hydrogen molecular dynamics study

We discuss molecular dynamics (MD) computer simulations of a tetracosane (C24H50) monolayer physisorbed onto the basal plane of graphite. The adlayer molecules are simulated with explicit hydrogens, and the graphite substrate is represented as an all-atom structure having six graphene layers. The tetracosane dynamics modeled in the fully atomistic manner agree well with experiment. The low-temperature ordered solid organizes into a rectangularly centered structure that is not commensurate with underlying graphite. Above T=200 K, as the molecules start to lose their translational and orientational order via gauche defect formation a weak smectic mesophase (observed experimentally but never reproduced in united atom (UA) simulations) appears. The phase behavior of the adsorbed layer is critically sensitive to the way the electrostatic interactions are included in the model. If the electrostatic charges are set to zero (as for a UA force field), then the melting temperature increases by approximately 70 K with respect to the experimental value. When the nonbonded 1-4 interaction is not scaled, the melting temperature decreases by approximately 90 K. If the scaling factor is set to 0.5, then melting occurs at T=350 K, in very good agreement with experimental data.