skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Thermochemistry of per‐ and polyfluoroalkyl substances
Abstract The determination of gas phase thermochemical properties of per‐ and polyfluoroalkyl substances (PFAS) is central to understanding the long‐range transport behavior of PFAS in the atmosphere. Prior gas‐phase studies have reported the properties of perfluorinated sulfonic acid (PFOS) and perfluorinated octanoic acid (PFOA). Here, this study reports the gas phase enthalpies of formation of short‐ and long‐chain PFAS and their precursor molecules determined using density functional theory (DFT) and ab initio approaches. Two density functionals, two ab initio methods and an empirical method were used to compute enthalpies of formation with the total atomization approach and an isogyric reaction. The performance of the computational methods employed in this work were validated against the experimental enthalpies of linear alkanoic acids and perfluoroalkanes. The gas‐phase determinations will be useful for future studies of PFAS in the atmosphere, and the methodological choices will be helpful in the study of other PFAS.  more » « less
Award ID(s):
2107155
PAR ID:
10458833
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Journal of Computational Chemistry
Volume:
44
Issue:
4
ISSN:
0192-8651
Page Range / eLocation ID:
570 to 580
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Per- and polyfluoroalkyl substances (PFAS) are a group of synthetic organofluorine surfactants that are resistant to typical methods of degradation. Thermal techniques along with other novel, less energy-intensive techniques are currently being investigated for the treatment of PFAS-contaminated matrices. Non-equilibrium plasma is one technique that has shown promise for the treatment of PFAS-contaminated water. To better tailor non-equilibrium plasma systems for this application, knowledge of the energy required for mineralization, and in turn the roles that plasma reactive species and heat can play in this process, would be useful. In this study, fundamental thermodynamic equations were used to estimate the enthalpies of reaction (480 kJ/mol) and formation (−4640 kJ/mol) of perfluorooctanoic acid (PFOA, a long-chain legacy PFAS) in water. This enthalpy of reaction estimate indicates that plasma reactive species alone cannot catalyze the reaction; because the reaction is endothermic, energy input (e.g., heat) is required. The estimated enthalpies were used with HSC Chemistry software to produce a model of PFOA defluorination in a 100 mg/L aqueous solution as a function of enthalpy. The model indicated that as enthalpy of the reaction system increased, higher PFOA defluorination, and thus a higher extent of mineralization, was achieved. The model results were validated using experimental results from the gliding arc plasmatron (GAP) treatment of PFOA or PFOS-contaminated water using argon and air, separately, as the plasma gas. It was demonstrated that PFOA and PFOS mineralization in both types of plasma required more energy than predicted by thermodynamics, which was anticipated as the model did not take kinetics into account. However, the observed trends were similar to that of the model, especially when argon was used as the plasma gas. Overall, it was demonstrated that while energy input (e.g., heat) was required for the non-equilibrium plasma degradation of PFOA in water, a lower energy barrier was present with plasma treatment compared to conventional thermal treatments, and therefore mineralization was improved. Plasma reactive species, such as hydroxyl radicals (⋅OH) and/or hydrated electrons (e−(aq)), though unable to accelerate an endothermic reaction alone, likely served as catalysts for PFOA mineralization, helping to lower the energy barrier. In this study, the activation energies (Ea) for these species to react with the alpha C–F bond in PFOA were estimated to be roughly 1 eV for hydroxyl radicals and 2 eV for hydrated electrons. 
    more » « less
  2. ABSTRACT In this work, we benchmark the accuracy of the density‐functional tight‐binding (DFTB) method, namely the long‐range corrected second‐order (LC‐DFTB2) and third‐order (DFTB3) models, for predicting energetics of imidazolium‐based ionic liquid (IL) ion pairs. We compare the DFTB models against popular density functionals such as LC‐ωPBE and B3LYP, using ab initio domain‐based local pair‐natural orbital coupled cluster (DLPNO‐CC) energies as reference. Calculations were carried out in the gas phase, as well as in aqueous solution using implicit solvent methods. We find that the LC‐DFTB2 model shows excellent performance in the gas phase and agrees well with reference energies in implicit solvent, often outperforming DFTB3 predictions for complexation energetics. Our study identifies a range of opportunities for use of the LC‐DFTB method and quantifies its sensitivity to protonation states and the types of chemical interactions between ion pairs. 
    more » « less
  3. Poly- and Perfluorinated alkyl substances (PFAS) pose environmental and public health concerns. While incineration remains the most common PFAS remediation method, the complete combustion and pyrolysis mechanism of PFAS is unknown. This study aims to expand our understanding of the kinetics of gas-phase PFAS incineration by measuring the effect of difluoromethane (CHF) on propane ignition delay times (IDTs). The ignition delay times were measured by OH* emission and end-wall pressure time histories behind the reflected shock wave. Different concentrations of CH2F2 were mixed with fuel-lean propane-oxygen mixtures diluted in argon. Experiments were conducted at a nominal reflected shock pressure of P5 = 1 atm and reflected shock temperatures of 1200 < T5 < 1800 K. A new detailed chemical kinetic mechanism is presented. 135 new rate constants were computed using RRKM/ME theory, based upon stationary points computed using ANL0. The new mechanism is in excellent agreement with the measured ignition delay time. A novel sensitivity analysis helps to explain the elementary steps by which CH2F2 increases the ignition delay time. 
    more » « less
  4. Although most manufacturers stopped using long-chain per- and polyfluoroalkyl substances (PFASs), including perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS), short-chain PFASs are still widely employed. Short-chain PFASs are less known in terms of toxicity and have different adsorption behavior from long-chain PFASs. Previous studies have shown electrostatic interaction with the adsorbent to be the dominant mechanism for the removal of short-chain PFASs. In this study, we designed a high charge density cationic quaternized nanocellulose (QNC) to enhance the removal of both short- and long-chain PFASs from contaminated water. Systematic batch adsorption tests were conducted using the QNC adsorbent to compare its efficiency against PFASs with varying chain lengths and functional groups. From the kinetic study, PFBA (perfluorobutanoic acid), PFBS (perfluorobutanesulfonic acid) and PFOS showed rapid adsorption rates, which reached near equilibrium values (>95% of removal) between 1 min to 15 min, while PFOA required a relatively longer equilibration time of 2 h (it obtained 90% of removal within 15 min). According to the isotherm results, the maximum adsorption capacity ( Q m ) of the QNC adsorbent exhibited the following trend: PFOS ( Q m = 559 mg g −1 or 1.12 mmol g −1 ) > PFOA ( Q m = 405 mg g −1 or 0.98 mmol g −1 ) > PFBS ( Q m = 319 mg g −1 or 1.06 mmol g −1 ) > PFBA ( Q m = 121 mg g −1 or 0.57 mmol g −1 ). This adsorption order generally matches the hydrophobicity trend among four PFASs associated with both PFAS chain length and functional group. In competitive studies, pre-adsorbed short-chain PFASs were quickly desorbed by long-chain PFASs, suggesting that the hydrophobicity of the molecule played an important role in the adsorption process on to QNC. Finally, the developed QNC adsorbent was tested to treat PFAS-contaminated groundwater, which showed excellent removal efficiency (>95%) for long-chain PFASs (C7–C9) even at a low adsorbent dose of 32 mg L −1 . However, short-chain PFASs ( i.e. , PFBA and perfluoropentanoic acid (PFPeA)) were poorly removed by the QNC adsorbent (0% and 10% removal, respectively) due to competing constituents in the groundwater matrix. This was further confirmed by controlled experiments that revealed a drop in the performance of QNC to remove short-chain PFASs at elevated ionic strength (NaCl), but not for long-chain PFASs, likely due to charge neutralization of the anionic functional group of PFASs by inorganic cations. Overall, the QNC adsorbent featured improved PFAS adsorption capacity at almost two-fold of PFAS removal by granular activated carbons, especially for short-chain PFASs. We believe, QNC can complement the use of common treatment methods such as activated carbon or ionic exchange resin to remove a wide range of PFAS pollutants, heading towards the complete remediation of PFAS contamination. 
    more » « less
  5. The hydrated electron is of interest to both theorists and experimentalists as a paradigm solution-phase quantum system. Although the bulk of the theoretical work studying the hydrated electron is based on mixed quantum/classical (MQC) methods, recent advances in computer power have allowed several attempts to study this object using ab initio methods. The difficulty with employing ab initio methods for this system is that even with relatively inexpensive quantum chemistry methods such as density functional theory (DFT), such calculations are still limited to at most a few tens of water molecules and only a few picoseconds duration, leaving open the question as to whether the calculations are converged with respect to either system size or dynamical fluctuations. Moreover, the ab initio simulations of the hydrated electron that have been published to date have provided only limited analysis. Most works calculate the electron’s vertical detachment energy, which can be compared to experiment, and occasionally the electronic absorption spectrum is also computed. Structural features, such as pair distribution functions, are rare in the literature, with the majority of the structural analysis being simple statements that the electron resides in a cavity, which are often based only on a small number of simulation snapshots. Importantly, there has been no ab initio work examining the temperature-dependent behavior of the hydrated electron, which has not been satisfactorily explained by MQC simulations. In this work, we attempt to remedy this situation by running DFT-based ab initio simulations of the hydrated electron as a function of both box size and temperature. We show that the calculated properties of the hydrated electron are not converged even with simulation sizes up to 128 water molecules and durations of several tens of picoseconds. The simulations show significant changes in the water coordination and solvation structure with box size. Our temperature-dependent simulations predict a red-shift of the absorption spectrum (computed using TD-DFT with an optimally tuned range-separated hybrid functional) with increasing temperature, but the magnitude of the predicted red-shift is larger than that observed experimentally, and the absolute position of the calculated spectra are off by over half an eV. The spectral red-shift at high temperatures is accompanied by both a partial loss of structure of the electron’s central cavity and an increased radius of gyration that pushes electron density onto and beyond the first solvation shell. Overall, although ab initio simulations can provide some insights into the temperature-dependent behavior of the hydrated electron, the simulation sizes and level of quantum chemistry theory that are currently accessible are inadequate for correctly describing the experimental properties of this fascinating object. 
    more » « less