pKa prediction of per- and polyfluoroalkyl acids in water using in silico gas phase stretching vibrational frequencies and infrared intensities

Abstract
To successfully understand and model the environmental fate of per- and polyfluoroalkyl substances (PFAS), it is necessary to know key physicochemical properties (PChPs) such as pKa; however, measured PChPs of PFAS are scarce and of uncertain reliability. In this study, quantitative structure–activity relationships (QSARs) were developed by correlating calculated (M062-X/aug-cc-pVDZ) vibrational frequencies (VF) and corresponding infrared intensities (IRInt) to the pKa of carboxylic acids, sulfonic acids, phosphonic acids, sulfonamides, betaines, and alcohols. Antisymmetric stretching VF of the anionic species were used for all subclasses except for alcohols where the OH stretching VF performed better. The individual QSARs predicted the pKa for each subclass mostly within 0.5 pKa units from the experimental values. The inclusion of IRInt as a pKa predictor for carboxylic acids improved the results by decreasing the root-mean-square error from 0.35 to 0.25 (n > 100). Application of the developed QSARs to estimate the pKa of PFAS within each subclass revealed that the length of the perfluoroalkyl chain has minimal effect on the pKa, consistent with other models but in stark contrast with the limited experimental data available.
Description
This article was originally published in Physical Chemistry Chemical Physics. The version of record is available at: https://doi.org/10.1039/D3CP01390A. This article will be embargoed until 09/01/2024.
Keywords
Citation
J. Murillo-Gelvez, O. Dmitrenko, T. Lorena Torralba Sánchez, P. G. Tratnyek and D. Di Toro, Phys. Chem. Chem. Phys., 2023, DOI: 10.1039/D3CP01390A.