Source-linked AI summary
Unraveling the PFAS helix: A statistical approach
Pranoy Ray, Haden Cavalli, Gashaw Bizana, Andrew R. Castillo, Shubham Vyas, Ronald L. Siefert, Surya R. Kalidindi, Manoj Kolel-Veetil
TL;DR
PFAS helicity has lacked a continuous quantitative metric despite its relevance to environmental persistence and degradation. The paper derives three statistical descriptors from carbon-backbone coordinates and validates them against geometric and spectroscopic benchmarks. Persistent-homology scores provide the strongest continuous descriptor, with agreement extending across several PFAS subfamilies.
Problem
PFAS helicity is typically described with isolated dihedral angles, which do not capture global fold topology, and a robust continuous quantitative metric is absent.
Source: §1 · p. 2Method
Three statistical descriptors derived from sorted carbon-backbone coordinates are cross-validated against dihedral benchmarks and DFT-calculated VCD spectra.
Source: §1 · p. 3, §2.2.1 · p. 5, §2.2.3 · p. 6Results
PH PC1 is the strongest continuous helicity descriptor, and geometric and topological measures agree across PFCA, PFSA, PFHS, FTOH, and PFOA/PFOS subfamilies.
Source: §1 · p. 3, §4 · p. 17Takeaways & Limitations
The validated continuous descriptor can be incorporated into QSPR models of PFAS environmental fate and degradation reactivity.
Source: §4 · p. 17The descriptors are evaluated on single minimum-energy conformers; Boltzmann-weighted conformer ensembles are a stated extension for flexible or partially fluorinated chains.
Source: §4 · p. 18
Abstract
from arXiv · showhide
The extreme persistence of per- and polyfluoroalkyl substances (PFAS) in the environment is rooted in their three-dimensional molecular conformation. The helical twist adopted by perfluoroalkyl chains due to hyperconjugation involving their recalcitrant C-F bonds governs their resistance to degradation, yet a quantitative, continuous metric for helicity has remained absent. Here we introduce a data-driven framework that quantifies backbone helicity using three statistical descriptors: void-state (binary occupancy) spatial autocorrelations, local backbone principal component analysis, and persistent homology. We further validate each statistical descriptor against geometric dihedral benchmarks and DFT-calculated Vibrational Circular Dichroism (VCD) spectra. The framework is initially established on a homologous series of perfluorocarboxylic acids (FC2 through FC16) and their hydrogenated analogues, then extended across perfluorosulfonic acids, fluorotelomer alcohols, polyfluoroalkyl hexanoic acid analogues, and longer-chain PFOA and PFOS analogues. Agreement between statistical, geometric, and spectroscopic definitions of helicity is established for PFCAs and extended to structurally distinct subfamilies, including PFOA analogues of varying fluorine content and perfluorosulfonic acids, demonstrating that the descriptors are robust to changes in headgroup chemistry and fluorine substitution pattern. The framework provides a chemistry-agnostic toolset for incorporating backbone conformation into predictive models of PFAS environmental fate and degradation reactivity.
1. INTRODUCTION
PFAS helicity is linked to backbone stability and proposed degradation pathways, but has lacked a robust continuous measure. The paper develops three backbone-based statistical descriptors and validates them against geometric and spectroscopic benchmarks.
1. INTRODUCTION
- PFAS persistence is attributed primarily to stable C-F bonds and fluorine shielding, while helical twisting is stabilized by dipole interactions and hyperconjugation.The paper notes that nascent helical character appears in chains as short as three carbons, with pronounced compact twists beyond five carbons.Source: §1 · pp. 1–2
- The missing robust, continuous helicity metric limits quantitative characterization beyond isolated dihedral angles; the study tests coordinate-based descriptors against hydrogenated controls.The comparison uses PFCAs FC2–FC16 and corresponding hydrogenated n-alkyl carboxylic acids to isolate fluorination-related structural differences.Source: §1 · p. 2
- The three statistical approaches use sorted carbon-backbone coordinates; their relevance to degradation is framed alongside a proposed, but not directly confirmed, linearization-facilitated pathway.The proposed pathway links electron capture to helical-to-linear transition, while the radical-anion intermediate remains an open experimental challenge.Source: §1 · p. 2
- Persistent-homology principal-component scores are the strongest continuous helicity descriptor, resolving chain-length progression and validated against dihedral benchmarks and DFT-calculated VCD spectra.The descriptors are compared with C-C-C-C backbone dihedral deviation from planarity and median F-C-C-F dihedral distributions.Source: §1 · p. 3
2. METHODOLOGY
The methodology applies statistical descriptors to ordered carbon-backbone coordinates from DFT-optimized PFAS structures, then cross-validates them with geometric dihedral benchmarks and VCD spectra. Its six stages establish controlled molecular geometries, derive global, local, and topological features, and define independent geometric checks.
2.1 DFT CALCULATIONS
- Single minimum-energy conformers of FC2–FC16 perfluorinated acids and matched hydrogenated acids were optimized and verified as true minima by frequency calculations.The paired series holds headgroup chemistry constant while varying chain length and fluorination.Source: §2.1 · p. 3
2.2 FEATURE ENGINEERING
- All descriptors use the ordered carbon backbone, isolated from each molecule; statistical descriptors are cross-validated against dihedral benchmarks and DFT-calculated VCD spectra.The workflow derives three complementary features from the sorted backbone coordinates.Source: §2.2 · pp. 3, 5, Figure 1 · p. 4
2.2.1. Voxelization & Spatial Autocorrelations
- Binary void-state voxelization and spatial autocorrelation produce a global shape descriptor; globally fitted PCA uses Spatial PC2 to separate fluorinated from hydrogenated chains.Spatial PC1 is dominated by chain length; Spatial PC2 explains 5.9% of variance and captures cross-sectional morphology.Source: §2.2.1 · p. 5
2.2.2. Local Backbone Principal Component Analysis
- Per-molecule PCA treats PC1 as the chain axis, PC2 as in-plane bending, and Local PC3 EVR as a scalar measure of out-of-plane backbone twist.A strictly linear chain has Local PC3 EVR near zero.Source: §2.2.2 · pp. 5–6
2.2.3. Persistent Homology Fingerprints
- Persistent homology converts pose-normalized backbone geometry into persistence-image fingerprints, with PH PC1 retained as a per-molecule descriptor.Arc-length normalization and canonical alignment reduce sensitivity to rigid-body pose and local inter-carbon length scale while preserving backbone-shape information.Source: §2.2.3 · p. 6
2.2.4. Dihedral Geometric Benchmarks
- The independent geometric benchmarks are C-C-C-C dihedral deviation from planarity and median F-C-C-F, or H-C-C-H, dihedral distributions; neither is used as an input feature.Zero deviation indicates an all-anti planar conformation; helical perfluoroalkyl chains show F-C-C-F values around 60°–80°, versus near-180° anti values in hydrogenated analogues.Source: §2.2.4 · p. 6
3. RESULTS & DISCUSSION
Across PFCA homologues, statistical descriptors distinguish fluorinated from hydrogenated backbones and resolve a continuous increase in helicity with chain length. Geometric and VCD benchmarks support these descriptors, which also track fluorine-substitution differences in fixed-length PFOA and PFOS analogues.
3.1. Separation of Fluorinated and Hydrogenated Classes by Backbone Shape
- Spatial PC2 separates fluorinated from hydrogenated chains, Local PC3 EVR distinguishes their out-of-plane structure, and PH PC1 orders fluorinated acids by increasing chain-length helicity.Local PC3 EVR is not monotonic within the fluorinated series because chain elongation captures a growing fraction of coordinate variance; short chains remain close to hydrogenated analogues in PH PC1.Source: §3.1 · pp. 7–8, Figure 2 · p. 7
3.2. Geometric Validation of Statistical Descriptors
- PH PC1 is the strongest descriptor against PFCA geometric benchmarks, correlating with C-C-C-C deviation from planarity (r = +0.84) and median F-C-C-F dihedral (r = +0.69).The planarity benchmark saturates around 17°-19° and has limited within-series range; pooled correlations also include class-separation effects.Source: §3.2 · p. 10, Figure 4 · p. 11
- Independent experimental observations, including FC9 crystal structure and solid-state 19F MAS NMR, are consistent with ordered helical perfluoroalkyl conformations.The FC9 crystal study also benchmarked its refined geometry against a DFT-optimized structure.Source: §3.2 · p. 11
3.3. Spectroscopic validation with Vibrational Circular Dichroism (VCD)
- DFT VCD activity is negligible for hydrogenated analogues and increases from clearly nonzero levels at FC8 through FC16; PH PC1 correlates more tightly with VCD PC1 than Local PC3 EVR (r = -0.95 versus r = -0.82).The computed single-sense VCD magnitude is a chiroptical probe, not a prediction of net signal from a racemic ensemble, whose signal cancels by symmetry.Source: §3.3 · pp. 12, 14, Figure 6 · p. 14
- Across the PFCA series, PH PC1 is the strongest descriptor against both VCD summaries, followed by Local PC3 EVR and Spatial PC2.The benchmarks include integrated absolute VCD intensity in the 1000-1500 cm-1 C-F stretching band and VCD PC1.Source: Figure 4 · p. 11, §3.3 · p. 14
3.4. Generalizability across fixed-chain-length subfamilies with varying fluorine content
- In fixed-length PFOA and PFOS analogues, higher fluorine content is associated with larger helicity-descriptor values, and PH PC1 best tracks dihedral and VCD benchmarks across the tested subfamilies.Spatial PC1, rather than PC2, captures substitution-driven structural variation in the PFOA set; the PFOS set shows the same overall pattern.Source: §3.4 · pp. 15–16
4. CONCLUSIONS
The framework quantifies PFAS helicity as a continuum, with agreement across geometric and topological measures in multiple PFAS subfamilies. PH PC1 is the strongest continuous descriptor across geometric and spectroscopic benchmarks, while the current evaluation uses single minimum-energy conformers.
4. CONCLUSIONS
- PH PC1 is the most physically faithful continuous helicity descriptor, with the strongest correlations to median F-C-C-F dihedrals and both spectroscopic benchmarks across PFCA and expanded multi-class datasets.Its correlation with backbone dihedral deviation from planarity is comparable to Local PC3 EVR.Source: §4 · p. 17
- Helical character forms a continuum: statistical descriptors resolve increasing twist with chain length, while clearly nonzero VCD activity begins at FC8 and fluorine substitution alone distinguishes fixed-chain-length analogues.Agreement spans PFCA, PFSA, PFHS, FTOH, and PFOA/PFOS subfamilies.Source: §4 · p. 17
- Continuous helicity descriptors can be incorporated into QSPR models linking backbone conformation to environmental fate, bioavailability, and degradation reactivity.Using PH PC1 to screen substrates and reaction conditions for a helical-to-linear transition preceding defluorination is motivated by a computational hypothesis, not an established mechanism.Source: §4 · p. 17
- The feature-engineering approach characterizes structure–property relationships without relying on traditional human-defined chemical descriptors and may extend to other systems where backbone conformation governs fate or reactivity.The authors identify fluorinated chains and synthetic helical polymers as examples.Source: §4 · p. 17
- The descriptors are evaluated on single minimum-energy conformers; Boltzmann-weighted conformer ensembles are a proposed extension for flexible or partially fluorinated chains with multiple thermally accessible conformers.Source: §4 · p. 18