Atmospheric Degradation of Two Hydrofluoroketones: Theoretical Rate Constants for the Gas-Phase OH-Oxidation of HFK-447mcc and HFK-465mc

Abstract

Accurate calculation of rate constants for gas-phase OH-oxidation reactions of fluorinated compounds is crucial for the understanding of atmospheric processes that are subject of the Kigali Agreement. Here, we have determined two such rate constants for two hydrofluoroketones, HFK-447mcc and HFK-465mc. The calculations were performed with a cost-effective multiconformer transition state theory protocol coupled with the constrained transition state randomization sampling method. The calculated rate constants of $k({\mathrm{HO}}^{•}+\mathrm{HFK}-447\mathrm{mcc})=3.1\times {10}^{-15}{\mathrm{cm}}^3{\mathrm{molecule}}^{-1}{\mathrm{s}}^{-1}$ and $k({\mathrm{HO}}^{•}+\mathrm{HFK}-465\mathrm{mc})=3.2\times {10}^{-14}{\mathrm{cm}}^3{\mathrm{molecule}}^{-1}{\mathrm{s}}^{-1}$ at 298.15 K imply an atmospheric lifetime of 10 years and 1 year, respectively. To our knowledge, these rate constants have never been determined experimentally or theoretically, and the similarity between the ratios of these two rate constants and of the well-studied acetone and diethyl ketone suggest the validity of our approach toward obtaining accurate rate constants and branching ratios.

1. Introduction

After it became clear that CFC compounds were directly connected to ozone depletion and global warming [1,2,3,4,5], theoretical and experimental efforts in the search for environmental-friendly CFC replacements have been vastly increased. This field of research gained extra relevance after October 2016, when almost 200 nations adopted the Kigali Amendment to the Montreal Protocol [6]. The ratification of this protocol urges countries to phase down production and the use of hydrofluorocarbons (HFCs), thus reinforcing the urgency of finding environmentally friendly alternatives, as these second-generation replacements have [7] elevated values of global warming potential (GWP). One of the most important parameters defining the GWP of a compound is the value of its atmospheric lifetime [8]. One potential solution for third-generation replacements in industrial and commercial applications [9] is provided by fluorine-containing oxygenated volatile organic compounds (OVOCs), since their atmospheric lifetimes are frequently decreased due to the presence of oxygen-containing functional groups. These functional groups promote long-range activating effects with respect to hydrogen-atom abstraction at sites remote from the substituent group, increasing the atmospheric reactivity when compared with the corresponding non-oxygenated fluorinated compounds [10,11,12]. Compounds such as fluorinated alcohols, aldehydes, ethers and ketones have a variety of industrial uses [9] and can gradually substitute HFCs, such as required by the Kigali Amendment to the Montreal Protocol.

Specifically, the atmospheric chemistry of ketones is important because they can be emitted to the atmosphere from anthropogenic (e.g., chemical industries) and biogenic (e.g., plants) sources, while they can also be generated in situ by the oxidation of other VOCs in the atmosphere [11,12]. The main atmospheric degradation processes of ketones are photolysis (ketones absorb radiation in the ultraviolet region around 280 nm due to the $n\to {\pi }^{*}$ transition of the carbonyl group) and reaction with OH radicals, with the former being more relevant at higher altitudes [10,12]. Tropospheric oxidation reactions with OH radicals typically proceed via hydrogen atom abstraction [10,12,13,14], with the reactivity of alkyl $\mathrm{C}-\mathrm{H}$ bonds at the $\beta$-position being enhanced with respect to the $\mathrm{C}-\mathrm{H}$ bonds at the $\alpha$-position due to the deactivation/activation of the hydrogen atom by the carbonyl group at the $\alpha$/$\beta$ positions, respectively [15,16]. In other words, the rate coefficients for the ${\mathrm{C}}_{\mathrm{n}}$ alkanones generally exceed those of the equivalent ${\mathrm{C}}_{\mathrm{n}-1}$ alkanes [11]. Large alkyl groups will lead to hydrogen abstraction from the various sites with the respective alkoxy radicals reacting with ${\mathrm{O}}_2$ and NO and generating a large variety of products [10,12].

Because of the recognized environment-friendly effect of the carbonyl group on the atmospheric lifetimes of ketones, hydrofluoroketones (HFKs) are considered to be potential third-generation replacements for CFCs and have been used as industrial synthesis reagents and fire extinguishing agents [9]. Naturally, for a particular HFK to be seriously considered as a potential replacement, it is necessary to determine if it has the required physicochemical properties and technical specifications. An example of such analysis has been performed on four HFKs [17]: HFK-447mcc, HFK-465mc, HFK-354pc and HFK-549mccc. The mentioned study concluded that HFK-354pc is thermally unstable, therefore rendering this particular HFK non-viable as a third-generation replacement option. Moreover, and to our knowledge, there are no experimental or theoretical investigations concerning the atmospheric lifetime of any of the remaining three HFKs. Such studies would provide an indication of the atmospheric impact of these three HFKs, and, assuming that they could all be efficiently used for the same purpose, it could also provide crucial information to potential stakeholders on which is the greener option (greener in the sense of Green Chemistry principle number 10 [18,19], design for degradation). Consequently, it is therefore important to gather relevant information about the atmospheric lifetime of these HFKs.

Considering their two main degradation pathways, the global lifetime of ketones can be approximated as [8,20]

(1)$\frac{1}{{\tau }_{\mathrm{HFK}}^{\mathrm{global}}}=\frac{1}{{\tau }_{\mathrm{HFK}}^{\mathrm{photo}.}}+\frac{1}{{\tau }_{\mathrm{HFK}}^{\mathrm{OH}}},$

where ${\tau }_{\mathrm{HFK}}^{\mathrm{photo}.}$ represents the HFK atmospheric lifetime resulting from the photolytic process with rate constant J and ${\tau }_{\mathrm{HFK}}^{\mathrm{OH}}$ is the atmospheric lifetime associated with a pseudo-first-order loss process with a corresponding $k_{\mathrm{OH}}\left[\mathrm{OH}\right]$ rate constant. However, in this study, we will discard the calculation of ${\tau }_{\mathrm{HFK}}^{\mathrm{photo}.}$ as we will be focusing only on the reaction between HFKs and the OH radical, for which the respective atmospheric lifetimes can be obtained by

(2)${\tau }_{\mathrm{HFK}}^{\mathrm{OH}}=\frac{1}{k_{\mathrm{OH}}\left[\mathrm{OH}\right]}.$

Although photolysis is generally considered to be more relevant at higher altitudes [10,12] (larger overlap between the $n\to {\pi }^{*}$ absorption band and the available actinic flux [21]), the complete approach would require the calculation of the photolysis rate constant J in order to determine ${\tau }_{\mathrm{HFK}}^{\mathrm{photo}.}$ and solve Equation (1), which is outside the scope of this investigation. Alternatively, we will provide a brief overview of what we can expect from the photolysis process. For example, the ${\tau }^{\mathrm{OH}}$ atmospheric lifetimes for 2-butanone, 2-pentanone and 2-hexanone range between 1.5 and 11 days [22], while the respective ${\tau }^{\mathrm{photo}.}$ lifetimes associated with $J_{\max }$ values (overhead Sun, corresponding to a solar zenith angle of 0 degrees, see Figure VIII-D-8 of Ref. [22]) range from 2.5 to 5.5 days. Under the mentioned conditions, the photolysis and OH-oxidation processes clearly compete. Interestingly, the fluorination of ketones leads to a bathochromic shift in the UV absorption bands of the carbonyl groups [23,24], enabling their overlap with actinic radiation at lower altitudes. Recent photolysis studies of several fluorinated ketones [23,24] under outdoor atmospheric conditions point to ${\tau }^{\mathrm{photo}.}$ lifetimes in the range of a few days to two weeks.

In this study, our calculations will focus solely on HFK-447mcc ($\mathrm{C}{\mathrm{F}}_3\mathrm{C}{F}_2\mathrm{C}{\mathrm{F}}_2\mathrm{C}\left(\mathrm{O}\right)\mathrm{C}{\mathrm{H}}_3$) and HFK-465mc ($\mathrm{C}{\mathrm{F}}_3\mathrm{C}{\mathrm{F}}_2\mathrm{C}\left(\mathrm{O}\right)\mathrm{C}{\mathrm{H}}_2\mathrm{C}{\mathrm{H}}_3$), excluding HFK-549mccc ($\mathrm{C}{\mathrm{F}}_3\mathrm{C}{\mathrm{F}}_2\mathrm{C}{\mathrm{F}}_2\mathrm{C}{\mathrm{F}}_2\mathrm{C}\left(\mathrm{O}\right)\mathrm{C}{\mathrm{H}}_3$). The reason for this is that the size of the C_xF_2x+1 group usually has little to no effect on the reactivity of halogenated OVOCs toward OH radicals [11], and therefore, HFK-447mcc and HFK-549mccc should have very similar rate constants. The main goal of this study will then be to calculate the k_OH rate constants for the oxidation reactions between two HFKs (HFK-447mcc and HFK-465mc) and OH radicals

(3)$\mathrm{H}\mathrm{F}\mathrm{K}+{\mathrm{H}\mathrm{O}}^{•}\stackrel{k_{\mathrm{O}\mathrm{H}}}{\to }{\mathrm{H}\mathrm{F}\mathrm{K}}^{•}+{\mathrm{H}}_2\mathrm{O}.$

This will be completed by using multiconformer transition state theory (MC-TST) [25,26,27,28,29,30,31,32] coupled with a newly developed transition state sampling method denominated constrained transition state randomization (CTSR) [33]. Our theoretical approach will be summarized in the next section, which is followed by a discussion of the results. The conclusions are gathered in the final section.

2. Methodology

2.1. Theoretical Background

The reaction between HFKs and ${\mathrm{HO}}^{•}$ is mechanistically similar to other radical-molecule reactions of atmospheric interest [34], including reactions between OVOCs and ${\mathrm{HO}}^{•}$ [10,13,14]. Such reaction paths feature loose complexes in both the reactant and product channels connected by a hydrogen-abstraction transition state (TS) with the underlying mechanism consisting of three steps: (1) formation of the prereaction complex (PRC) from the isolated OVOC + OH radical reactants, (2) formation of the product complex (PC) from the PRC, and (3) formation of the final products (OVOC radical and water) from the PC. For the specific case of an HFK, one obtains

(4)$\mathrm{step}1:\mathrm{HFK}+{\mathrm{HO}}^{•}\underset{k_{-1}}{\overset{k_1}{\rightleftharpoons }}\left[\mathrm{HFK}\cdots \mathrm{HO}\right]{}^{•}$

(5)$\mathrm{step}2:\left[\mathrm{HFK}\cdots \mathrm{HO}\right]{}^{•}\stackrel{k_2}{\to }\left[\mathrm{HFK}\cdots \mathrm{HOH}\right]{}^{•}$

(6)$\mathrm{step}3:\left[\mathrm{HFK}\cdots \mathrm{HOH}\right]{}^{•}\stackrel{k_3}{\to }{\mathrm{HFK}}^{•}+{\mathrm{H}}_2\mathrm{O}.$

If one assumes that the PRC concentration remains constant over time (steady-state assumption [35], with $k_{-1}+k_2\gg k_1$), the effective rate constant for reaction (3) can be calculated as [36,37,38,39]

(7)$k_{\mathrm{OH}}=\frac{k_1{k}_2}{k_{-1}+k_2}$

In the context of tropospheric reactions, it is usually found that the $k_{-1}$ reverse rate constants for step 1 are much larger than the $k_2$ rate constant values for the irreversible unimolecular step 2 (the conditions for which $k_{-1}\gg k_2$ are discussed in detail elsewhere [30,36,37,38,39,40,41]), which simplifies the overall rate constant to $k_{\mathrm{OH}}=k_1{k}_2/k_{-1}$. Under such conditions, the system is assumed to be fully equilibrated and thermalized at the high-pressure limit (equilibrium-state assumption [40]), and therefore, the reactants, the PRC and the TS of step 2 can be characterized by Boltzmann equilibrium distributions. These conditions should be checked before performing the calculations, particularly if the hydrogen-abstraction TSs are submerged below the reactants. This can be achieved by calculating the Gibbs free energies along the reaction path and verifying their relative values. In many cases, the Gibbs free energy profile along the reaction coordinate will resemble the profile represented in Figure 5b of Ref. [40], providing a strong indication that $k_{-1}$ will be considerably larger than $k_1$ and $k_2$ and that both the steady-state and equilibrium assumptions will be satisfied for the reaction under study [40]. In that case, $k_1/k_{-1}$ is identified as the $K_{\mathrm{eq}}$ thermal equilibrium constant for the fast pre-equilibrium between the isolated reactants and the PRC, which can be determined by employing the usual statistical thermodynamic tools [36,37,38,39]. Similarly, $k_2$ can also be calculated by resorting to thermal TST, simplifying our TST-based rate coefficient to $k_{\mathrm{OH}}=K_{\mathrm{eq}}{k}_2$. This approach, which we will utilize in this investigation, is frequently used to describe reactions between OVOCs and ${\mathrm{HO}}^{•}$, and it has the advantage of avoiding the need to deal with the PRCs, since their partition functions and energies cancel out in the $K_{\mathrm{eq}}{k}_2$ product.

Alternatively, our recent version of the MC-TST/CTSR protocol [33,41] considers the approach followed by Petit and Harvey [27], where both the numerator and denominator of Equation (7) are multiplied by $K_{\mathrm{eq}}$, yielding

(8)$k_{\mathrm{OH}}=\frac{k_1{k}_2{K}_{\mathrm{eq}}}{k_1+k_2{K}_{\mathrm{eq}}}=\frac{k_1{k}_{\mathrm{OH}}^{\mathrm{MC}-\mathrm{TST}}}{k_1+k_{\mathrm{OH}}^{\mathrm{MC}-\mathrm{TST}}}\equiv k_{\mathrm{OH}}\left(\mathrm{calc}\right).$

Here, $k_1$ is the forward rate coefficient for the formation of the PRC and $k_{\mathrm{OH}}^{\mathrm{MC}-\mathrm{TST}}$ can be identified as our working equation for the calculation of the MC-TST rate constants. Naturally, the $k_1$ rate constant is unique for each specific reaction under analysis. However, since these $k_1$ rate coefficients are almost independent of the collider, we used [27] a value of $k_1$ for collision with acetone of ($k_1=(2.67\pm 0.15)\times {10}^{-11}{\mathrm{cm}}^3{\mathrm{molecule}}^{-1}{\mathrm{s}}^{-1}$) obtained via experimental rate constants [42,43,44,45] for the quenching of vibrationally excited ${\mathrm{HO}}^{•}$.

2.2. Computational Methods

We will now proceed to the calculation of $k_{\mathrm{OH}}^{\mathrm{MC}-\mathrm{TST}}$. Our cost-effective protocol [30,33,41] considers N conformers for each HFK reactant and M possible hydrogen abstraction TSs, with the rate coefficient being written as

(9)$k_{\mathrm{OH}}^{\mathrm{MC}-\mathrm{TST}}=\kappa \left(T\right)\frac{k_{\mathrm{B}}T}{c^{\circ }h{Q}_{\mathrm{OH}}}\frac{{\sum }_i^M\frac{{\alpha }_i{\omega }_{{\mathrm{TS}}_i}{Q}_{{\mathrm{TS}}_i}}{{\sigma }_{\mathrm{rot},{\mathrm{TS}}_i}}}{{\sum }_j^N\frac{{\alpha }_j{\omega }_{{\mathrm{HFK}}_j}{Q}_{{\mathrm{HFK}}_j}}{{\sigma }_{\mathrm{rot},{\mathrm{HFK}}_j}}}{\mathrm{e}}^{-V^{‡}/k_{B}T}.$

Here, $c^{\circ }=p^{\circ }/k_{\mathrm{B}}T$ is the standard concentration at a temperature of 298.15 K and a pressure of 1 atmosphere, while $k_{\mathrm{B}}$ and h are the Boltzmann and Planck’s constants, respectively. The Q values are the total partition functions for ${\mathrm{HO}}^{•}$, the ith HFK reactant and jth TS, which are evaluated as the product of the translational, electronic, vibrational and rotational partition functions ($Q=q_{\mathrm{t}}{q}_{\mathrm{e}}{q}_{\mathrm{v}}{q}_{\mathrm{r}}$) under the rigid rotor and harmonic oscillator approximation. For each HFK or a TS conformer, every total partition function is weighted by its respective ${\omega }_{{\mathrm{HFK}}_j}$ or ${\omega }_{{\mathrm{TS}}_i}$ thermal weight factor, which depends on the relative $E_{{\mathrm{HFK}}_j}$ and $E_{{\mathrm{TS}}_i}$ zero-point corrected energies of conformations ${\mathrm{HFK}}_j$ and ${\mathrm{TS}}_i$, respectively. In addition, ${\sigma }_{\mathrm{rot},{\mathrm{HFK}}_j}$, ${\sigma }_{\mathrm{rot},{\mathrm{TS}}_i}$, ${\alpha }_j$ and ${\alpha }_i$ are the rotational symmetry numbers and reaction path degeneracy parameters of conformations ${\mathrm{HFK}}_j$ and ${\mathrm{TS}}_i$, respectively. The calculation of the rotational symmetry numbers is performed by first determining the respective point group symmetry with a brute force symmetry analyzer [46] and then associating the obtained group symmetry with the correct factor by examining Table 2 of Ref. [26]. The $V^{‡}$ parameter is calculated [26,28,29,30] with the formula $V^{‡}=E_{{\mathrm{TS}}_0}-(E_{{\mathrm{HFK}}_0}+E_{\mathrm{OH}})$ where the subscript 0 indicates a lowest energy conformer. The $\kappa \left(T\right)$ tunneling correction is calculated [28,29,30] via the Eckart formula [47] using both the imaginary frequency characterizing ${\mathrm{TS}}_0$ and the zero-point corrected energies of the PRC and PC minima that are connected to ${\mathrm{TS}}_0$ through the respective intrinsic reaction coordinate path.

The above parameters are obtained through our cost-effective MC-TST/CTSR protocol [30,33,41], which we will briefly recapitulate. The protocol is initiated with the calculation of the reactant conformers with the use of confab [48], which is a knowledge-based conformer generation tool that utilizes force fields (for our protocol, we are employing the MMFF94 [49] force field) and a database of allowed torsion angles. These conformations are then minimized with the same force field as in the confab calculation by employing the Obminimize program. Both programs are available in the Open Babel [50] chemical toolbox. The minimized structures are then used as starting geometries for optimization at the M08-HX/pcseg-2 [51,52] level of theory. A similarity evaluation (SE) between all generated conformers is then performed with a superimposing algorithm [53]. The purpose of the SE is to identify identical structures and enantiomers in our pool of optimized reactants, which we can then safely eliminate, with the remaining structures constituting the final set of N distinguishable reactants.

The next step concerns the TS sampling, which was performed with the new CTSR method [33]. The CTSR code first generates random OH-molecule structures which are then subjected to constraints targeting an optimization of TSs at the PM3 level of theory. The three constraints for the HFK + ${\mathrm{HO}}^{•}$ TSs were the following: (1) the distance between the ${\mathrm{HO}}^{•}$ oxygen atom and the hydrogen atom under attack must be in the interval [1.2, 1.4] Å; (2) the distance between the ${\mathrm{HO}}^{•}$ oxygen atom and all reactants’ non-hydrogen atoms must be greater than 1.70 Å and (3) the distance between the ${\mathrm{HO}}^{•}$ hydrogen atom and all reactants’ hydrogen atoms must be greater than 0.70 Å. Guess structures were then generated by the CTSR code and subsequently used as input geometry in PM3 saddle-point optimizations. The converged PM3 calculations were then analyzed, with the coordinates corresponding to the TSs yielding an imaginary frequency within a specific interval being subjected to an SE step. The surviving geometries were used as guess structures for saddle-point optimizations at the M08-HX/pcseg-2 level, with the geometries of the successful optimizations being again analyzed in terms of their imaginary frequencies. The geometries filtered by the previous step were then put through another SE, so that the final set of unique TSs was obtained. It should be stated that the M08-HX functional was chosen due to its quality in describing hydrogen abstraction barrier heights, yielding, respectively, mean unsigned errors of 1.00 and 0.73 $\mathrm{kcal}{\mathrm{mol}}^{-1}$ for the DBH76 and HTBH38 data sets involving this particular type of barrier heights [51]. In this study, M08-HX was the best performing functional, surpassing M06-2X, which is routinely used in atmospheric chemistry studies of the same type. Moreover, we recently [54] made a rough estimate of 0.6 $\mathrm{kcal}{\mathrm{mol}}^{-1}$ for the uncertainty associated with $V^{‡}$, which led to rate constants and error bars with the approximate values of $k\pm 0.98k$. It should also be noted that although our protocol relies on many approximations and on error cancellation [30,41], the quality of its calculated rate constants is quite encouraging [33,41] when taking into account the complexity of the studied systems. The electronic structure calculations necessary to the determination of the parameters entering Equation (9) were performed with the GAMESS package [55]. Finally, the value of $V^{‡}$ was refined with single-point energy calculations at the M08-HX/apcseg-2 level.

3. Results and Discussion

We start by showing Figure 1, where we present the simulated IR spectra of both reactants, HFK-447mcc and HFK465mc, based on the M08-HX/pcseg-2 harmonic calculations and including the contribution of all conformers through the respective Boltzmann weighting at 298.15 K, which was calculated from the Gibbs energies. Each spectrum was obtained by convolution of the calculated vibrational frequencies ($\nu$) with a Lorentzian profile given by

(10)$I\left(\nu \right)=\sum _j^{3N-6}\left[\frac{I\left({\nu }_j\right)}{1+4{[(\nu -{\nu }_j)/W]}^2}\right],$

where N is the number of atoms in each molecule, ${\nu }_j$ is the calculated vibrational frequency with a corresponding $I\left({\nu }_j\right)$ intensity and W is the full width at half-maximum, which we empirically selected as $50{\mathrm{cm}}^{-1}$.

For saturated aliphatic ketones, the $\mathrm{C}=\mathrm{O}$ and $\mathrm{C}-\mathrm{H}$ stretch bands are usually located at 1715 and [2800,3000] ${\mathrm{cm}}^{-1}$, respectively. For the studied ketones, our calculations predict these bands to be found at slightly higher wavenumbers, as represented in Figure 1. However, it should be noted that although the M08-HX/pcseg-2 model chemistry is a good choice for MC-TST rate constants [30,33,41,56], it might not be the best theoretical method to employ if one is pursuing an accurate prediction of the vibrational frequencies [57].

3.1. OH Oxidation of HFK-447mcc

Starting the analysis of our results with the $\mathrm{C}{\mathrm{F}}_3\mathrm{C}{\mathrm{F}}_2\mathrm{C}{\mathrm{F}}_2\mathrm{C}\left(\mathrm{O}\right)\mathrm{C}{\mathrm{H}}_3\mathrm{H}\mathrm{F}\mathrm{K}$, our confab-based [48] conformer search [30,41] relying on the MMFF94 [49] force field and followed by M08-HX/pcseg-2 optimizations found 20 unique conformers of $C_1$ symmetry (each one belonging to an enantiomeric pair by definition), with the 15 lowest energy ones having a relative energy in the interval [0.00,0.45] $\mathrm{kcal}{\mathrm{mol}}^{-1}$ and representing approximately 90% of the total fractional population of reactants. Application of the CTSR sampling procedure yielded 90 distinguishable hydrogen abstraction TSs of $C_1$ symmetry, with the $V^{‡}$ barrier having a value of 4.2 $\mathrm{kcal}{\mathrm{mol}}^{-1}$. At this point, it becomes advisable to test the validity of our theoretical kinetic approach by investigating if the steady-state and equilibrium assumptions [40] are applicable to the reaction under scrutiny. For this, we have gathered the Gibbs free energies for the reactants, ${\mathrm{TS}}_0$, and the associated PRC and PC minima at the M08-HX/pcseg-2 level of theory. The $\Delta G$ results for the reactants, PRC, ${\mathrm{TS}}_0$ and PC, are, respectively, 0.0, 4.7, 12.6 and $-19.1\mathrm{kcal}{\mathrm{mol}}^{-1}$. These values indicate that the Gibbs free energy profile along the reaction coordinate is similar to the one shown in Figure 5b of Ref. [40], indicating that $k_{-1}$ will most likely be considerably larger than $k_1$ and $k_2$, guaranteeing that both the steady-state and equilibrium assumptions are satisfied for the reaction between HFK-447mcc and ${\mathrm{HO}}^{•}$ at 298.15 K.

The Eckart tunneling correction yielded a value of 6.3, with the final $k_{\mathrm{OH}}\left(\mathrm{calc}\right)$ rate constant for the reaction between HFK-447mcc and ${\mathrm{HO}}^{•}$ being $k_{\mathrm{OH}}^{\mathrm{HFK}-447\mathrm{mcc}}\left(\mathrm{calc}\right)=3.1\times {10}^{-15}{\mathrm{cm}}^3{\mathrm{molecule}}^{-1}{\mathrm{s}}^{-1}$. Such a rate constant implies a ${\tau }_{\mathrm{HFK}}^{\mathrm{OH}}$ atmospheric lifetime of approximately 10 years, considering a value of $1\times {10}^6\mathrm{molecules}{\mathrm{cm}}^{-3}$ for the global average concentration of OH radicals [58]. Additionally, we have extended our calculations to a temperature interval relevant for tropospheric chemistry [21,59,60], namely $230\le T/\mathrm{K}\le 320$, and fitted the corresponding data to an Arrhenius–Kooij equation [59,61,62] to give $k\left(T\right)=2.61\times {10}^{-18}exp(2126/T){(T/300)}^{10.42}{\mathrm{cm}}^3{\mathrm{molecule}}^{-1}{\mathrm{s}}^{-1}$. Cartesian coordinates for reactants and TSs can be found in the Supplementary Information, while Figure 2 shows ball and stick drawings of the lowest energy HFK-447mcc and corresponding TS (${\mathrm{HFK}}_0$ and ${\mathrm{TS}}_0$, respectively).

3.2. OH Oxidation of HFK-465mc

Proceeding to the $\mathrm{C}{\mathrm{F}}_3\mathrm{C}{\mathrm{F}}_2\mathrm{C}\left(\mathrm{O}\right)\mathrm{C}{\mathrm{H}}_2\mathrm{C}{\mathrm{H}}_3\mathrm{H}\mathrm{F}\mathrm{K}$, our calculations found 16 unique conformers of $C_1$ symmetry, with the seven lowest energy ones having a relative energy in the interval [0.00,0.75] $\mathrm{kcal}{\mathrm{mol}}^{-1}$ and representing slightly more than 90% of the total fractional population of reactants. Application of the CTSR sampling procedure yielded 52 distinguishable hydrogen abstraction TSs of $C_1$ symmetry, with the $V^{‡}$ barrier being calculated as 2.1 $\mathrm{kcal}{\mathrm{mol}}^{-1}$. For this reaction, the relative Gibbs free energies for the reactants, PRC, ${\mathrm{TS}}_0$ and PC minima at the M08-HX/pcseg-2 level of theory are, respectively: 0.0, 4.2, 10.8 and $-26.1\mathrm{kcal}{\mathrm{mol}}^{-1}$. As in the case above, this strongly suggests that both the steady-state and equilibrium assumptions are satisfied for the reaction between HFK-465mc and ${\mathrm{HO}}^{•}$ at 298.15 K.

The Eckart tunneling correction yielded a value of 2.7, with the final $k_{\mathrm{OH}}\left(\mathrm{calc}\right)$ rate constant for the reaction between HFK-465mc and ${\mathrm{HO}}^{•}$ being $k_{\mathrm{OH}}^{\mathrm{HFK}-465\mathrm{mc}}\left(\mathrm{calc}\right)=3.2\times {10}^{-14}{\mathrm{cm}}^3{\mathrm{molecule}}^{-1}{\mathrm{s}}^{-1}$. Such a rate constant implies a ${\tau }_{\mathrm{HFK}}^{\mathrm{OH}}$ atmospheric lifetime of approximately 1 year. Moreover, the fit of the theoretical data points in the $230\le T/\mathrm{K}\le 320$ interval to the Arrhenius–Kooij equation leads to $k\left(T\right)=3.14\times {10}^{-16}exp(1389/T){(T/300)}^{6.83}{\mathrm{cm}}^3{\mathrm{molecule}}^{-1}{\mathrm{s}}^{-1}$. Cartesian coordinates for reactants and TSs can be found in the Supplementary Information, while Figure 3 shows ball and stick drawings of the lowest energy HFK-465mc and corresponding TS (${\mathrm{HFK}}_0$ and ${\mathrm{TS}}_0$, respectively).

In Figure 4, we present the $k_{\mathrm{OH}}\left(\mathrm{calc}\right)$ points calculated in the $230\le T/\mathrm{K}\le 320$ temperature interval together with the Arrhenius–Kooij fits to those points for both reactions. Although both curves exhibit some curvature, it is clearly noticeable that the curvature of the HFK-447mcc + OH reaction is more pronounced mainly because of the differences in the tunneling correction factors, which increase with decreasing temperature: ${\kappa }^{\mathrm{HFK}-447\mathrm{mcc}}\left(320\mathrm{K}\right)=2.05{\kappa }^{\mathrm{HFK}-465\mathrm{mc}}\left(320\mathrm{K}\right)$ and ${\kappa }^{\mathrm{HFK}-447\mathrm{mcc}}\left(230\mathrm{K}\right)=5.24{\kappa }^{\mathrm{HFK}-465\mathrm{mc}}\left(230\mathrm{K}\right)$.

The reaction between HFK-465mc and ${\mathrm{HO}}^{•}$ can proceed through two different abstraction sites, as the OH radical can abstract a hydrogen atom in two distinct locations of the $\mathrm{C}{\mathrm{F}}_3\mathrm{C}{\mathrm{F}}_2\mathrm{C}\left(\mathrm{O}\right)\mathrm{C}{\mathrm{H}}_2\mathrm{C}{\mathrm{H}}_3\mathrm{H}\mathrm{F}\mathrm{K}$: at the $-\mathrm{C}{\mathrm{H}}_2$ group ($\alpha$-position) and at the $-\mathrm{C}{\mathrm{H}}_3$ group ($\beta$-position). The calculated branching ratios for the two abstraction sites are: ${\mathsf{\Gamma }}_{\mathrm{C}{\mathrm{H}}_2}=64\%$ and ${\mathsf{\Gamma }}_{\mathrm{C}{\mathrm{H}}_3}=36\%$. Interestingly, these values coincide with the branching ratios calculated for the corresponding unsubstituted diethyl ketone (DEK, $\mathrm{C}{\mathrm{H}}_3\mathrm{C}{\mathrm{H}}_2\mathrm{C}\left(\mathrm{O}\right)\mathrm{C}{\mathrm{H}}_2\mathrm{C}{\mathrm{H}}_3$) with the use of structure–activity relationships (SARs). Using the $k_{\left(\mathrm{CHx}\right)}$ additivity method [10,11,16,63] with the recommended values of $k_{\left(\mathrm{CHx}\right)}$ for the reaction of OH with ketones (see Table 4 of Ref. [10] or Table V-B-20 of Ref. [11]), we can determine $k_{\mathrm{OH}}^{\mathrm{DEK}}$ as

(11)$k_{\mathrm{OH}}^{\mathrm{DEK}}=2k_{\beta \left(-{\mathrm{CH}}_3-\right)}+2k_{\alpha \left(-{\mathrm{CH}}_2-\right)},$

which yields a value of $k_{\mathrm{OH}}^{\mathrm{DEK}}=2.14\times {10}^{-12}{\mathrm{cm}}^3{\mathrm{molecule}}^{-1}{\mathrm{s}}^{-1}$ and corresponding branching ratios of ${\mathsf{\Gamma }}_{\mathrm{C}{\mathrm{H}}_2}=65\%$ and ${\mathsf{\Gamma }}_{\mathrm{C}{\mathrm{H}}_3}=35\%$. Note that the recommended value [11] for this rate coefficient is $2\times {10}^{-12}{\mathrm{cm}}^3{\mathrm{molecule}}^{-1}{\mathrm{s}}^{-1}$ with an uncertainty of $\pm 30\%$.

Moreover, it is also observed that the ratio between the calculated OH-oxidation rate constants of HFK-465mc and HFK-447mcc is approximately 10, which reflects the different reactivity of the $-\mathrm{C}{\mathrm{H}}_3$ and $-{\mathrm{C}}_2{\mathrm{H}}_5$ alkyl groups. An almost equivalent ratio is obtained if the same type of comparison is performed with unsubstituted ketones: the ratio between the recommended values [11] for the OH-oxidation rate constants of DEK and acetone ($\mathrm{C}{\mathrm{H}}_3\mathrm{C}\left(\mathrm{O}\right)\mathrm{C}{\mathrm{H}}_3$, with a rate constant of $1.8\times {10}^{-13}{\mathrm{cm}}^3{\mathrm{molecule}}^{-1}{\mathrm{s}}^{-1}$) is approximately 11. In the absence of experimental rate coefficients concerning the HFK-447mcc and HFK-465mc molecules, such an agreement may be interpreted as an indication that despite all approximations, the MC-TST/CTSR protocol [30,33,41] correctly captures essential reactivity trends.

4. Conclusions

In this study, we have calculated the rate constants for the reactions between two HFKs (HFK-447mcc and HFK-465mc) and the OH radical, which is the most important oxidant in the atmosphere. The results, which were obtained via theoretical calculations employing our MC-TST/CTSR protocol, show that the rate constant for the reaction involving HFK-465mc is ten times larger than the rate constant involving HFK-447mcc, which results in an atmospheric lifetime of approximately 1 and 10 years for HFK-465mc and HFK-447mcc, respectively. These results indicate that photolysis will probably be the main atmospheric sink for HFK-447mcc and HFK-465mc. Although there are no experimental studies to compare our results to, the observed reactivity pattern and branching ratios resemble the ones observed in well-tested systems with the same reactive alkyl groups, such as acetone and diethyl ketone. Despite the fact that these are good indications and that our MC-TST protocol consistently yields accurate rate constants for reactions between fluorinated compounds and the OH radical, a full assessment of these two rate constants can only be made in the future if and when the corresponding experimental rate coefficients are determined.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available alongside the article and the Supplementary Materials.

Acknowledgments

L.P.V. acknowledges Dorothee Birke for proofreading this manuscript.

Conflicts of Interest

The author declares no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures

Figure 1. Boltzmann weighted (298.15 K) harmonic M08-HX/pcseg-2 calculated IR spectra of HFK-447mcc and HFK-465mc, including all obtained conformers for both species.

Figure 2. Ball and stick drawings for the lowest energy HFK-447mcc reactant and respective [Forumla omitted. See PDF.] for reaction (3).

Figure 3. Ball and stick drawings for the lowest energy HFK-465mc reactant and respective [Forumla omitted. See PDF.] for reaction (3).

Figure 4. Arrhenius–Kooij fits (solid lines) and [Forumla omitted. See PDF.] data points (solid circles) for both reactions in the [Forumla omitted. See PDF.] temperature interval.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/atmos13081256/s1, Cartesian coordinates for all unique conformers used in the calculation of both MC-TST rate constants present in the manuscript.

References

1. Molina, M.J.; Rowland, F.S. Stratospheric Sink for Chlorofluoromethanes: Chlorine Atom-Catalysed Destruction of Ozone. Nature; 1974; 249, pp. 810-812. [DOI: https://dx.doi.org/10.1038/249810a0]

2. Farman, J.D.; Gardiner, B.G.; Shanklin, J.D. Large Losses of Total Ozone in Antarctica Reveal Seasonal ClO_x/NO_x Interaction. Nature; 1985; 315, pp. 207-210. [DOI: https://dx.doi.org/10.1038/315207a0]

3. United Nations: Montreal Protocol on Substances that Deplete the Ozone Layer-Adjustments and Amendment. Int. Leg. Mater.; 1993; 32, pp. 874-887. [DOI: https://dx.doi.org/10.1017/S0020782900016338]

4. Lazarou, Y.G.; Papagiannakopoulos, P. Theoretical Investigation of the Thermochemistry of Hydrofluoroethers. Chem. Phys. Lett.; 1999; 301, pp. 19-28. [DOI: https://dx.doi.org/10.1016/S0009-2614(98)01455-9]

5. McGivern, B.P. Conference of the Parties to the Framework Convention on Climate Change: Kyoto Protocol. Int. Leg. Mater.; 1998; 37, pp. 22-43. [DOI: https://dx.doi.org/10.1017/S0020782900019409]

6. Heath, E.A. Amendment to the Montreal Protocol on Substances that Deplete the Ozone Layer (Kigali Amendment). Int. Leg. Mater.; 2017; 56, pp. 193-205. [DOI: https://dx.doi.org/10.1017/ilm.2016.2]

7. UNEP. HFCs: A Critical Link in Protecting Climate and the Ozone Layer; United Nations Environment Programme (UNEP): 2011; 36p. Available online: https://www.unep.org/resources/report/hfcs-critical-link-protecting-climate-and-ozone-layer (accessed on 6 August 2022).

8. Hodnebrog, Ø.; Etminan, M.; Fuglestvedt, J.S.; Marston, G.; Myhre, G.; Nielsen, C.J.; Shine, K.P.; Wallington, T.J. Global Warming Potentials and Radiative Efficiencies of Halocarbons and Related Compounds: A Comprehensive Review. Rev. Geophys.; 2013; 51, pp. 300-378. [DOI: https://dx.doi.org/10.1002/rog.20013]

9. Wallington, T.J.; Andersen, M.P.S.; Nielsen, O.J. Atmospheric Chemistry of Halogenated Organic Compounds. Advances in Atmospheric Chemistry; World Scientific: Singapore, 2017; Volume 1, pp. 305-402.

10. Mellouki, A.; Le Bras, G.; Sidebottom, H. Kinetics and Mechanisms of the Oxidation of Oxygenated Organic Compounds in the Gas Phase. Chem. Rev.; 2003; 103, pp. 5077-5096. [DOI: https://dx.doi.org/10.1021/cr020526x]

11. Calvert, J.G.; Mellouki, A.; Orlando, J.J.; Pilling, M.J.; Wallington, T.J. The Mechanisms of Atmospheric Oxidation of the Oxygenates; 1st ed. Oxford University Press: New York, NY, USA, 2011.

12. Mellouki, A.; Wallington, T.J.; Chen, J. Atmospheric Chemistry of Oxygenated Volatile Organic Compounds: Impacts on Air Quality and Climate. Chem. Rev.; 2015; 115, pp. 3984-4014. [DOI: https://dx.doi.org/10.1021/cr500549n]

13. Vereecken, L.; Francisco, J.S. Theoretical Studies of Atmospheric Reaction Mechanisms in the Troposphere. Chem. Soc. Rev.; 2012; 41, pp. 6259-6293. [DOI: https://dx.doi.org/10.1039/c2cs35070j]

14. Vereecken, L.; Glowacki, D.R.; Pilling, M.J. Theoretical Chemical Kinetics in Tropospheric Chemistry: Methodologies and Applications. Chem. Rev.; 2015; 115, pp. 4063-4114. [DOI: https://dx.doi.org/10.1021/cr500488p] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25843545]

15. Atkinson, R.; Aschmann, S.M.; Carter, W.P.L.; Pitts, J.N., Jr. Rate Constants for the Gas-Phase Reaction of OH Radicals with a Series of Ketones at 299 ± 2 K. Int. J. Chem. Kinet.; 1982; 14, pp. 839-847. [DOI: https://dx.doi.org/10.1002/kin.550140804]

16. Le Calvé, S.; Hitier, D.; Le Bras, G.; Mellouki, A. Kinetic Studies of OH Reactions with a Series of Ketones. J. Phys. Chem. A; 1998; 102, pp. 4579-4584. [DOI: https://dx.doi.org/10.1021/jp980848y]

17. Otake, K.; Yasumoto, M.; Yamada, Y.; Murata, J.; Urata, S. Critical Parameters and Vapor Pressure Measurements of Potential Replacements for Chlorofluorocarbons - Four Hydrofluoroketones and a Hydrofluoroamine. J. Chem. Eng. Data; 2003; 48, pp. 1380-1383. [DOI: https://dx.doi.org/10.1021/je020210m]

18. Anastas, P.T.; Warner, J.C. Green Chemistry: Theory and Practice; Oxford University Press: Oxford, UK, 1998.

19. Erythropel, H.C.; Zimmerman, J.B.; de Winter, T.M.; Petitjean, L.; Melnikov, F.; Lam, C.H.; Lounsbury, A.W.; Mellor, K.E.; Janković, N.Z.; Tu, Q. et al. The Green ChemisTREE: 20 Years After Taking Root With the 12 Principles. Green Chem.; 2018; 20, pp. 1929-1961. [DOI: https://dx.doi.org/10.1039/C8GC00482J]

20. Kurylo, M.J.; Orkin, V.L. Determination of Atmospheric Lifetimes via the Measurement of OH Radical Kinetics. Chem. Rev.; 2003; 103, pp. 5049-5076. [DOI: https://dx.doi.org/10.1021/cr020524c]

21. Finlayson-Pitts, B.J.; Pitts, J.N., Jr. Chemistry of the Upper and Lower Atmosphere: Theory, Experiments, and Applications; 1st ed. Academic Press: San Diego, CA, USA, 2000.

22. Calvert, J.G.; Orlando, J.J.; Stockwell, W.R.; Wallington, T.J. The Mechanisms of Reactions Influencing Atmospheric Ozone; 1st ed. Oxford University Press: New York, NY, USA, 2015.

23. Díaz-de-Mera, Y.; Aranda, A.; Notario, A.; Rodríguez, A.; Rodríguez, D.; Bravo, I. Photolysis Study of Fluorinated Ketones Under Natural Sunlight Conditions. Phys. Chem. Chem. Phys.; 2015; 17, pp. 22991-22998. [DOI: https://dx.doi.org/10.1039/C5CP03527A]

24. Ren, Y.; Bernard, F.; Daële, V.; Mellouki, A. Atmospheric Fate and Impact of Perfluorinated Butanone and Pentanone. Environ. Sci. Technol.; 2019; 53, pp. 8862-8871. [DOI: https://dx.doi.org/10.1021/acs.est.9b02974]

25. Vereecken, L.; Peeters, J. The 1,5-H-Shift in 1-Butoxy: A Case Study in the Rigorous Implementation of Transition State Theory for a Multirotamer System. J. Chem. Phys.; 2003; 119, pp. 5159-5170. [DOI: https://dx.doi.org/10.1063/1.1597479]

26. Fernández-Ramos, A.; Ellingson, B.A.; Meana-Pañeda, R.; Marques, J.M.C.; Truhlar, D.G. Symmetry Numbers and Chemical Reaction Rates. Theor. Chem. Acc.; 2007; 118, pp. 813-826. [DOI: https://dx.doi.org/10.1007/s00214-007-0328-0]

27. Petit, A.S.; Harvey, J.N. Atmospheric Hydrocarbon Activation by the Hydroxyl Radical: A Simple yet Accurate Computational Protocol for Calculating Rate Coefficients. Phys. Chem. Chem. Phys.; 2012; 14, pp. 184-191. [DOI: https://dx.doi.org/10.1039/C1CP21367A] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22073400]

28. Rissanen, M.P.; Kurtén, T.; Sipilä, M.; Thornton, J.A.; Kangasluoma, J.; Sarnela, N.; Junninen, H.; Jørgensen, S.; Schallhart, S.; Kajos, M.K. et al. The Formation of Highly Oxidized Multifunctional Products in the Ozonolysis of Cyclohexene. J. Am. Chem. Soc.; 2014; 136, pp. 15596-15606. [DOI: https://dx.doi.org/10.1021/ja507146s] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25283472]

29. Møller, K.H.; Otkjaer, R.V.; Hyttinen, N.; Kurtén, T.; Kjaergaard, H.G. Cost-Effective Implementation of Multiconformer Transition State Theory for Peroxy Radical Hydrogen Shift Reactions. J. Phys. Chem. A; 2016; 120, pp. 10072-10087. [DOI: https://dx.doi.org/10.1021/acs.jpca.6b09370] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27992191]

30. Viegas, L.P. Exploring the Reactivity of Hydrofluoropolyethers toward OH Through a Cost-Effective Protocol for Calculating Multiconformer Transition State Theory Rate Constants. J. Phys. Chem. A; 2018; 122, pp. 9721-9732. [DOI: https://dx.doi.org/10.1021/acs.jpca.8b08970]

31. Novelli, A.; Vereecken, L.; Bohn, B.; Dorn, H.; Gkatzelis, G.L.; Hofzumahaus, A.; Holland, F.; Reimer, D.; Rohrer, F.; Rosanka, S. et al. Importance of Isomerization Reactions for OH Radical Regeneration from the Photo-Oxidation of Isoprene Investigated in the Atmospheric Simulation Chamber SAPHIR. Atmos. Chem. Phys.; 2020; 20, pp. 3333-3355. [DOI: https://dx.doi.org/10.5194/acp-20-3333-2020]

32. Vereecken, L.; Nozière, B. H Migration in Peroxy Radicals Under Atmospheric Conditions. Atmos. Chem. Phys.; 2020; 20, pp. 7429-7458. [DOI: https://dx.doi.org/10.5194/acp-20-7429-2020]

33. Viegas, L.P. Gas-Phase OH-Oxidation of 2-Butanethiol: Multiconformer Transition State Theory Rate Constant with Constrained Transition State Randomization. Chem. Phys. Lett.; 2022; 803, 139829. [DOI: https://dx.doi.org/10.1016/j.cplett.2022.139829]

34. Hansen, J.C.; Francisco, J.S. Radical-Molecule Complexes: Changing Our Perspective on the Molecular Mechanisms of Radical-Molecule Reactions and Their Impact on Atmospheric Chemistry. ChemPhysChem; 2002; 3, pp. 833-840. [DOI: https://dx.doi.org/10.1002/1439-7641(20021018)3:10<833::AID-CPHC833>3.0.CO;2-0]

35. Raines, R.T.; Hansen, D.E. An Intuitive Approach to Steady-State Kinetics. J. Chem. Educ.; 1988; 65, pp. 757-759. [DOI: https://dx.doi.org/10.1021/ed065p757]

36. Singleton, D.L.; Cvetanović, R.J. Temperature Dependence of the Reaction of Oxygen Atoms with Olefins. J. Am. Chem. Soc.; 1976; 98, pp. 6812-6819. [DOI: https://dx.doi.org/10.1021/ja00438a006]

37. Alvarez-Idaboy, J.R.; Mora-Diez, N.; Vivier-Bunge, A. A Quantum Chemical and Classical Transition State Theory Explanation of Negative Activation Energies in OH Addition To Substituted Ethenes. J. Am. Chem. Soc.; 2000; 122, pp. 3715-3720. [DOI: https://dx.doi.org/10.1021/ja993693w]

38. Galano, A.; Alvarez-Idaboy, J.R.; Bravo-Pérez, G.; Ruiz-Santoyo, M.E. Gas Phase Reactions of C₁–C₄ Alcohols with the OH Radical: A Quantum Mechanical Approach. Phys. Chem. Chem. Phys.; 2002; 4, pp. 4648-4662. [DOI: https://dx.doi.org/10.1039/B205630E]

39. Cruz-Torres, A.; Galano, A. On the Mechanism of Gas-Phase Reaction of C₁–C₃ Aliphatic Thiols + OH Radicals. J. Phys. Chem. A; 2007; 111, pp. 1523-1529. [DOI: https://dx.doi.org/10.1021/jp0673415] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17279738]

40. Xiao, F.; Sun, X.; Li, Z.; Li, X. Theoretical Study of Radical-Molecule Reactions with Negative Activation Energies in Combustion: Hydroxyl Radical Addition to Alkenes. ACS Omega; 2020; 5, pp. 12777-12788. [DOI: https://dx.doi.org/10.1021/acsomega.0c00400]

41. Viegas, L.P. Simplified Protocol for the Calculation of Multiconformer Transition State Theory Rate Constants Applied to Tropospheric OH-Initiated Oxidation Reactions. J. Phys. Chem. A; 2021; 125, pp. 4499-4512. [DOI: https://dx.doi.org/10.1021/acs.jpca.1c00683]

42. Smith, I.W.M. Collisional Energy Transfer, Intramolecular Vibrational Relaxation and Unimolecular Reactions. J. Chem. Soc. Faraday Trans.; 1997; 93, pp. 3741-3750. [DOI: https://dx.doi.org/10.1039/a704144f]

43. Smith, I.W.M.; Ravishankara, A.R. Role of Hydrogen-Bonded Intermediates in the Bimolecular Reactions of the Hydroxyl Radical. J. Phys. Chem. A; 2002; 106, pp. 4798-4807. [DOI: https://dx.doi.org/10.1021/jp014234w]

44. McCabe, D.C.; Brown, S.S.; Gilles, M.K.; Talukdar, R.K.; Smith, I.W.M.; Ravishankara, A.R. Kinetics of the Removal of OH(v = 1) and OD(v = 1) by HNO₃ and DNO₃ from 253 to 383 K. J. Phys. Chem. A; 2003; 107, pp. 7762-7769. [DOI: https://dx.doi.org/10.1021/jp0346413]

45. Talukdar, R.K.; Gierczak, T.; McCabe, D.C.; Ravishankara, A.R. Reaction of Hydroxyl Radical with Acetone. 2. Products and Reaction Mechanism. J. Phys. Chem. A; 2003; 107, pp. 5021-5032. [DOI: https://dx.doi.org/10.1021/jp0273023]

46. Patchkovskii, S. Brute Force Symmetry Analyzer. Available online: https://github.com/nquesada/symmetry (accessed on 6 August 2022).

47. Eckart, C. The Penetration of a Potential Barrier by Electrons. Phys. Rev.; 1930; 35, pp. 1303-1309. [DOI: https://dx.doi.org/10.1103/PhysRev.35.1303]

48. O’Boyle, N.M.; Vandermeersch, T.; Flynn, C.J.; Maguire, A.R.; Hutchison, G.R. Confab-Systematic Generation of Diverse Low-Energy Conformers. J. Cheminf.; 2011; 3, 8. [DOI: https://dx.doi.org/10.1186/1758-2946-3-8] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21410983]

49. Halgren, T.A. Merck Molecular Force Field. I. Basis, Form, Scope, Parameterization, and Performance of MMFF94. J. Comput. Chem.; 1996; 17, pp. 490-519. [DOI: https://dx.doi.org/10.1002/(SICI)1096-987X(199604)17:5/6<490::AID-JCC1>3.0.CO;2-P]

50. O’Boyle, N.M.; Banck, M.; James, C.A.; Vandermeersch, T.; Hutchison, G.R. Open Babel: An Open Chemical Toolbox. J. Cheminf.; 2011; 3, 33. [DOI: https://dx.doi.org/10.1186/1758-2946-3-33] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21982300]

51. Zhao, Y.; Truhlar, D.G. Exploring the Limit of Accuracy of the Global Hybrid Meta Density Functional for Main-Group Thermochemistry, Kinetics, and Noncovalent Interactions. J. Chem. Theory Comput.; 2008; 4, pp. 1849-1868. [DOI: https://dx.doi.org/10.1021/ct800246v] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26620329]

52. Jensen, F. Unifying General and Segmented Contracted Basis Sets. Segmented Polarization Consistent Basis Sets. J. Chem. Theory Comput.; 2014; 10, pp. 1074-1085. [DOI: https://dx.doi.org/10.1021/ct401026a]

53. Marques, J.M.C.; Llanio-Trujillo, J.L.; Abreu, P.E.; Pereira, F.B. How Different Are Two Chemical Structures?. J. Chem. Inf. Model.; 2010; 50, pp. 2129-2140. [DOI: https://dx.doi.org/10.1021/ci100219f]

54. Viegas, L.P.; Jensen, F. Reactivity of α,ω-Dihydrofluoropolyethers toward OH Predicted by Multiconformer Transition State Theory and the Interacting Quantum Atoms Approach. J. Phys. Chem. A; 2020; 124, pp. 3460-3470. [DOI: https://dx.doi.org/10.1021/acs.jpca.0c02911]

55. Schmidt, M.W.; Baldridge, K.K.; Boatz, J.A.; Elbert, S.T.; Gordon, M.S.; Jensen, J.H.; Koseki, S.; Matsunaga, N.; Nguyen, K.A.; Su, S. et al. General Atomic and Molecular Electronic Structure System. J. Comput. Chem.; 1993; 14, pp. 1347-1363. [DOI: https://dx.doi.org/10.1002/jcc.540141112]

56. Viegas, L.P. Assessment of Model Chemistries for Hydrofluoropolyethers: A DFT/M08-HX Benchmark Study. Int. J. Quantum Chem.; 2017; 117, e25381. [DOI: https://dx.doi.org/10.1002/qua.25381]

57. Palafox, M.A. DFT Computations on Vibrational Spectra: Scaling Procedures to Improve the Wavenumbers. Phys. Sci. Rev.; 2018; 3, 0184. [DOI: https://dx.doi.org/10.1515/psr-2017-0184]

58. Prinn, R.G.; Huang, J.; Weiss, R.F.; Cunnold, D.M.; Fraser, P.J.; Simmonds, P.G.; McCulloch, A.; Harth, C.; Salameh, P.; O’Doherty, S. et al. Evidence for Substantial Variations of Atmospheric Hydroxyl Radicals in the Past Two Decades. Science; 2001; 292, pp. 1882-1888. [DOI: https://dx.doi.org/10.1126/science.1058673] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11337586]

59. McGillen, M.R.; Carter, W.P.L.; Mellouki, A.; Orlando, J.J.; Picquet-Varrault, B.; Wallington, T.J. Database for the Kinetics of the Gas-Phase Atmospheric Reactions of Organic Compounds. Earth Syst. Sci. Data; 2020; 12, pp. 1203-1216. [DOI: https://dx.doi.org/10.5194/essd-12-1203-2020]

60. Altarawneh, M. A Closer Look Into the Contribution of Atmospheric Gas-Phase Pathways in the Formation of Perfluorocarboxylic Acids. Atmos. Pol. Res.; 2021; 12, 101255. [DOI: https://dx.doi.org/10.1016/j.apr.2021.101255]

61. Kooij, D.M. Über die Zersetzung des Gasförmigen Phosphorwasserstoffs. Z. Phys. Chem.; 1893; 12, pp. 155-161. [DOI: https://dx.doi.org/10.1515/zpch-1893-1214]

62. Laidler, K.J. The Development of the Arrhenius Equation. J. Chem. Ed.; 1984; 61, pp. 494-498. [DOI: https://dx.doi.org/10.1021/ed061p494]

63. Ren, Y.; Baramoussi, E.M.E.; Daële, V.; Mellouki, A. Atmospheric Chemistry of Ketones: Reaction of OH Radicals with 2-Methyl-3-Pentanone, 3-Methyl-2-Pentanone and 4-Methyl-2-Pentanone. Sci. Total Environ.; 2021; 780, 146249. [DOI: https://dx.doi.org/10.1016/j.scitotenv.2021.146249]