I. INTRODUCTION
Perfluorononanoic acid (PFNA or C9HF17O2, IUPAC name: heptadecafluorononanoic acid) is a compound belonging to the family of per- and polyfluoroalkyl substances (PFAS), characterized by alkyl chains decorated by fluorine atoms. Due to widespread industrial usage in diverse applications (surfactants, surface treatment agents for textiles, fire-fighting foams, and production of other fluorinated compounds, among others) and the strong nature of carbon-fluorine bonds making the compounds highly resistant to degradation, PFAS have become ubiquitous environmental contaminants (Dadashi Firouzjaei et al., Reference Dadashi Firouzjaei, Zolghadr, Ahmadalipour, Taghvaei, Afkhami, Nejati and Elliott2022; Evich et al., Reference Evich, Davis, McCord, Acrey, Awkerman, Knappe, Lindstrom, Speth, Tebes-Stevens, Strynar, Wang, Weber, Henderson and Washington2022). PFAS have been implicated in numerous health concerns such as liver and kidney disease, reproductive and developmental issues, and cancer (Fenton et al., Reference Fenton, Ducatman, Boobis, DeWitt, Lau, Ng, Smith and Roberts2021), but due to the variety of different compounds (numbering in the thousands) and complexity of human population studies, much still remains unknown about their health effects (Kirk et al., Reference Kirk, Smurthwaite, Braunig, Trevenar, D'Este, Lucas, Lal, Korda, Clements, Mueller and Armstrong2018). PFNA is prominently included in a toxicology literature database with 29 common PFAS compounds (Pelch et al., Reference Pelch, Reade, Kwiatkowski, Merced-Nieves, Cavalier, Schultz, Wolffe and Varshavsky2022).
In this work, synchrotron powder X-ray diffraction (PXRD) data were obtained for PFNA and used to solve the crystal structure. Subsequent Rietveld refinement strongly suggested the presence of disorder which we were unable to model. A structureless Le Bail refinement was used to extract a complete Bragg reflection list unbiased by the structural model for inclusion in the Powder Diffraction File (Kabekkodu et al., Reference Kabekkodu, Dosen and Blanton2024).
II. EXPERIMENTAL
A specimen of PFNA (Sigma-Aldrich, 97% purity) was lightly ground with a mortar and pestle, then inserted into a 1.2 mm inner diameter glass capillary for data collection.
PXRD data were collected using the Brockhouse X-ray Diffraction and Scattering sector Wiggler Low Energy (WLE) beamline (Leontowich et al., Reference Leontowich, Gomez, Moreno, Muir, Spasyuk, King, Reid, Kim and Kycia2021) at the Canadian Light Source (CLS). WLE is an in-vacuum wiggler beamline which employs a Si (111) single-side bounce crystal monochromator. Data were collected at a photon energy of ~15.14 keV (calibrated wavelength of 0.81906Å), using a series of eight Dectris Mythen2 X series 1K strip detectors on the end of the goniometer arm (sample-detector distance of 767 mm). The precise wavelength and instrument resolution function were determined using Rietveld refinement of a pattern obtained from a lanthanum hexaboride (LaB6) standard reference material (NIST SRM 660a LaB6) using GSASII (Toby and Von Dreele, Reference Toby and Von Dreele2013). The sample was spun at 2 Hz during data acquisition.
Indexing of the PFNA pattern with DICVOL06 (Boultif and Louër, Reference Boultif and Louër2004) suggested a monoclinic unit cell with lattice parameters a = 26.1585 Å, b = 5.6307 Å, c = 10.9423 Å, β = 98.759°, and a cell volume of 1592.9 Å3 (M 20 = 77.9, F 20 = 418.7). The cell volume was suitable for four formula units and space group determination with ChekCell (Laugier and Bochu, Reference Laugier and Bochu2000) suggested P2 1/c as the most likely space group. A PFNA molecule (minus the hydrogen atom) was obtained from the Cambridge Structural Database (Groom et al., Reference Groom, Bruno, Lightfoot and Ward2016) entry KATVUY (Motreff et al., Reference Motreff, da Costa, Allouchi, Duttine, Mathonière, Duboc and Vincent2012), and then converted to a Fenske–Hall Z-matrix with Open Babel (O'Boyle et al., Reference O'Boyle, Banck, James, Morley, Vandermeersch and Hutchison2011). The Fenske–Hall Z-matrix with the PFNA molecule was uploaded into FOX (Favre-Nicolin and Černý, Reference Favre-Nicolin and Černý2002) for structure solution with parallel tempering. Initial parallel tempering sets employing least-squares refinement every 250 000 trials were observed to significantly distort the PFNA molecule into unrealistic configurations. Subsequent parallel tempering sets employed rigid body constraints on the molecule without least-squares refinements, in order to find possible molecular packing arrangements without distorting the molecule. A number of different structures, including both dimers and columnar structures supporting zig-zag chains of hydrogen bonds, were investigated with Rietveld refinement. Ultimately, a dimer structure was chosen based on the best fit to the data.
Rietveld refinement of the crystal structures was performed with the GSASII program (Toby and Von Dreele, Reference Toby and Von Dreele2013) using a Thompson–Cox–Hastings modified pseudo-Voigt peak shape function (Thompson et al., Reference Thompson, Cox and Hastings1987). The background was modeled using a 16-term log interpolation function. Positional parameters were refined with restraints on bond distances and angles for all non-hydrogen atoms in the PFNA molecule obtained using Mogul (Bruno et al., Reference Bruno, Cole, Kessler, Luo, Motherwell, Purkis, Smith, Taylor, Cooper, Harris and Orpen2004), with the mean and standard deviation quantities used to prepare the restraints. Individual isotropic displacement parameters were constrained to be the same for each element, with the value for the hydrogen atom constrained to be 1.3 times the value of the oxygen atoms. A fourth-order spherical harmonic texture model (Von Dreele, Reference Von Dreele1997) was used to correct for preferred orientation (PO), utilizing eight refined parameters.
The crystal data, data collection, and refinement details are summarized in Table I.
The crystal data, data collection, and refinement parameters obtained for monoclinic PFNA
Density functional theory (DFT) geometry optimization was performed starting from the Rietveld-refined PFNA structure using CRYSTAL17 (Dovesi et al., Reference Dovesi, Erba, Orlando, Zicovich-Wilson, Civalleri, Maschio, Rérat, Casassa, Baima, Salustro and Kirtman2018). Basis sets were obtained from the literature for the C, H, O (Gatti et al., Reference Gatti, Saunders and Roetti1994), and F atoms (Vilela Oliveira et al., Reference Vilela Oliveira, Laun, Peintinger and Bredow2019). The calculations used the B3LYP functional (Becke, Reference Becke1993) and D3 dispersion correction (Grimme et al., Reference Grimme, Antony, Ehrlich and Krieg2010) with eight k-points.
III. RESULTS AND DISCUSSION
The final Rietveld refinement obtained for the PFNA pattern is illustrated in Figure 1 (top), with the final Le Bail refinement (Figure 1, bottom) plotted below for comparison. While the weighted R-factor is relatively low (wR = 3.14%, which includes contributions of both the raw data points and the restraints) and the main features of the pattern are well described, discrepancies between observed and calculated peak intensities can be observed in the minor Bragg reflections. The weighted R-factor for the Le Bail refinement (wR = 1.24%) is significantly lower, and visually the Le Bail fit is superior. The spherical harmonic PO correction yielded a texture index of 1.201(4) and improved the fit, but did not eliminate the discrepancies between the observed and calculated peak intensities in the Rietveld refinement.
Plots illustrating the final Rietveld refinement (wR = 3.14%, top) and Le Bail refinement (wR = 1.24%, bottom) obtained for PFNA with GSASII. The region above 12° is magnified by 6× for each refinement to aid visualization.
The final Rietveld-refined structure (red) is compared to the DFT-optimized structure (blue) for PFNA in Figure 2. The structures are dimers where the carboxylic acid groups form pairs of hydrogen bonds between molecules. The fundamental dimer structure, monoclinic unit cell, and space group are reasonably consistent with the hydrogenated analog of PFNA, nonanoic acid (C9H18O2), which was previously determined using single crystal data (Bond, Reference Bond2004). While the Rietveld-refined and DFT-optimized structures are similar, there is clearly a shift in the position of the DFT-optimized molecule, which is more pronounced toward the trifluoromethyl end. The root-mean-squared (RMS) Cartesian displacement between the Rietveld-refined and the DFT-optimized structures for the non-hydrogen atoms is 1.12 Å, which falls outside the range generally expected for correct crystal structures using powder diffraction data (Van de Streek and Neumann, Reference Van de Streek and Neumann2014). The RMS displacement obtained using the molecular overlay function of Mercury (Macrae et al., Reference Macrae, Bruno, Chisholm, Edgington, McCabe, Pidcock, Rodriguez-Monge, Taylor, Streek and Wood2008) is 0.515 Å. Mercury performs a least-squares alignment of the two molecules, resulting in the improved agreement compared with the absolute Cartesian displacement. The molecular overlay of the Rietveld-refined and DFT-optimized molecules is illustrated in Figure 3 and differences are observed, particularly with respect to the trifluoromethyl and carboxylic acid orientations at the ends of the molecules. Despite restraints, the intermolecular hydrogen bonds formed between the carboxylic acid groups are shorter in the Rietveld-refined structure (H29⋯O2 = 1.539(6) Å) than in the DFT-optimized structure (H29⋯O2 = 1.60 Å). The eight-atom ring formed by the two hydrogen-bonded carboxylic acid groups is also perfectly planar in the DFT-optimized structure, but slightly irregular in the Rietveld-refined structure, as illustrated in Figure 4. Arguably, the most unusual features of the experimental structure are the torsion angles about the C7–C10 axis (four torsion angles based on two fluorine atoms per carbon atom). These torsion angles fall outside the distributions given by Mogul for similar geometries. Calculation of the powder pattern using the DFT-optimized structure results in significantly different pattern intensities than the experimental pattern, suggesting it is somewhat different than the true experimental structure. The Rietveld-refined atomic coordinates are given in Table II, and the atomic labels used in the refinement are illustrated in Figure 5. The largest peak and hole in the difference Fourier map for the Rietveld-refined structure were 0.67 and −0.83 eÅ−3, respectively.
A comparison of the Rietveld-refined (red) and the DFT-optimized (blue) structures of PFNA, viewed along the b-axis (top) and the c-axis (bottom). The figure was prepared with Mercury (Macrae et al., Reference Macrae, Bruno, Chisholm, Edgington, McCabe, Pidcock, Rodriguez-Monge, Taylor, Streek and Wood2008).
A comparison of the molecular overlay of the Rietveld-refined (red) and the DFT-optimized (blue) structures of PFNA obtained with Mercury (Macrae et al., Reference Macrae, Bruno, Chisholm, Edgington, McCabe, Pidcock, Rodriguez-Monge, Taylor, Streek and Wood2008).
The Rietveld-refined (top) and DFT-optimized (bottom) crystal structure of PFNA, viewed along the b-axis. The atom types can be identified by color including carbon (black), fluorine (green), oxygen (red), and hydrogen (pink). Hydrogen bonds are represented by the dotted lines. The unit cell is outlined in black. The figure was prepared with VESTA (Momma and Izumi, Reference Momma and Izumi2011).
The Rietveld-refined crystal structure of PFNA with lattice parameters a = 26.172(1) Å, b = 5.6345(2) Å, c = 10.9501(4) Å, and β = 98.752(2)°.
All atoms belong on Wyckoff 4e sites and were refined with fixed site occupancies of 1.
The atom labels used for the PFNA structure.
Given the plausibility of the refined structure and the relatively large isotropic displacement parameters observed (see Table II), disorder was considered as a possible reason for the discrepancies between the Rietveld and Le Bail refinements. Disorder can be observed in long-chain molecular structures; for example, the structure of Teflon, (C2F4)n, can be disordered (Sprik et al., Reference Sprik, Roethlisberger and Klein1999). In an (unsuccessful) attempt to model the disorder, a structure solution was performed with FOX using two PFNA molecules with half occupancy to ascertain if a plausible overlapping structure with partial occupancies could be found. Despite running over 50 sets of parallel tempering, no plausible structures were observed with this approach, but this may be asking too much of the data.
In addition to the dimer structure, Rietveld refinement was performed on a columnar structure composed of double columns of PFNA molecules running parallel to the b-axis, with the molecules oriented perpendicular to the b-axis. The double columns were held together by a zig-zag ladder of hydrogen bonds. While the structure appeared plausible, the refinement was inferior to that obtained for the final dimer structure, and a DFT-optimization of the Rietveld-refined columnar structure failed to converge.
The most similar-length fluorinated compound to PFNA observed in the Cambridge Structural Database (CSD) is perfluorooctanoic acid (PFOA or C8HF15O2) contained in CSD entry KEQREG (Omorodion et al., Reference Omorodion, Palenzuela, Ruether, Twamley, Platts and Baker2018). PFOA crystallizes in a significantly different structure with triclinic symmetry, composed of single columns of PFOA molecules running parallel to the a-axis. Each column is held together by a chain of hydrogen bonds connecting carboxylic acid groups of adjacent molecules in the column.
The Hirshfeld surface (Hirshfeld, Reference Hirshfeld1977; Spackman and Byrom, Reference Spackman and Byrom1997) of the PFNA molecule obtained using Crystal Explorer (Spackman et al., Reference Spackman, Turner, McKinnon, Wolff, Grimwood, Jayatilaka and Spackman2021) is illustrated in Figure 6. The volume enclosed within the Hirshfeld surface (395.08 Å3) represents 99.0% of one-fourth of the unit cell volume, suggesting that the molecules are not tightly packed.
The Hirshfeld surface of PFNA. Regions in red represent intermolecular contacts shorter than the sum of the van der Waals radii, while regions in blue represent longer contacts and regions in white represent contacts equal to the sum of the van der Waal radii. The figure was prepared with Crystal Explorer (Spackman et al., Reference Spackman, Turner, McKinnon, Wolff, Grimwood, Jayatilaka and Spackman2021).
A Bragg reflection list for PFNA based on Le Bail refinement was prepared by summing reflections closer than 0.02°, 2θ as multiple reflections and assigning a weighted average reflection position, then including all reflections with relative integrated intensities of 0.2% or greater, up to 25° in 2θ.
IV. DEPOSITED DATA
Individual Crystallographic Information Framework (CIF) files containing the results of the final Rietveld refinement (crystal structure), Le Bail refinement (raw data and Bragg reflection list), and DFT-optimized structure were deposited with ICDD. The data can be requested at pdj@icdd.com.
ACKNOWLEDGEMENTS
JR would like to thank Jim Kaduk and Denis Spasyuk for helpful discussions regarding the crystal structure and suggestions which significantly improved the paper. MP and TM would like to acknowledge support from: (i) Gregory Holmbeck at Idaho National Laboratory through a Laboratory Directed Research & Development (LDRD) Program under U.S. Department of Energy (DOE) Idaho Operations Office Contract DE-AC07-05ID14517; and (ii) the U. S. DOE Office of Science Basic Energy Sciences program under Award Number DE-SC0023248.
Research described in this paper was performed using the WLE beamline at the Canadian Light Source, which is supported by the Canadian Foundation for Innovation, the Natural Sciences and Engineering Research Council of Canada, the National Research Council Canada, the Canadian Institutes of Health Research, the Government of Saskatchewan, Western Economic Diversification Canada, and the University of Saskatchewan.
CONFLICTS OF INTEREST
The authors have no conflicts of interest to declare.
Open access
