Simplified Protocol for the Calculation of Multiconformer Transition State Theory Rate Constants Applied to Tropospheric OH-Initiated Oxidation Reactions

Chemical kinetics plays a fundamental role in the understanding and modeling of tropospheric chemical processes, one of the most important being the atmospheric degradation of volatile organic compounds. These potentially harmful molecules are emitted into the troposphere by natural and anthropogenic sources and are chemically removed by undergoing oxidation processes, most frequently initiated by reaction with OH radicals, the atmosphere's "detergent". Obtaining the respective rate constants is therefore of critical importance, with calculations based on transition state theory (TST) often being the preferred choice. However, for molecules with rich conformational variety, a single-conformer method such as lowest-conformer TST is unsuitable while state-of-the-art TST-based methodologies easily become unmanageable. In this Feature Article, the author reviews his own cost-effective protocol for the calculation of bimolecular rate constants of OH-initiated reactions in the high-pressure limit based on multiconformer transition state theory. The protocol, which is easily extendable to other oxidation reactions involving saturated organic molecules, is based on a variety of freeware and open-source software and tested against a series of oxidation reactions of hydrofluoropolyethers, computationally very challenging molecules with potential environmental relevance. The main features, advantages and disadvantages of the protocol are presented, along with an assessment of its predictive utility based on a comparison with experimental rate constants.

Reactivity of α,ω-Dihydrofluoropolyethers toward OH Predicted by Multiconformer Transition State Theory and the Interacting Quantum Atoms Approach

We report rate constants for the tropospheric reaction between the OH radical and α,ω-dihydrofluoropolyethers, which represent a specific class of the hydrofluoropolyethers family with the formula HF₂C(OCF₂CF₂)_p(OCF₂)_qOCF₂H. Four cases were considered: p0q2, p0q3, p1q0, and p1q1 (pxqy denoting p = x and q = y) with the calculations performed by a cost-effective protocol developed for bimolecular hydrogen-abstraction reactions. This protocol is based on multiconformer transition state theory and relies on computationally accessible M08-HX/apcseg-2//M08-HX/pcseg-1 calculations. Within the protocol's approximations, the results show that (1) the calculated rate constants are within a factor of five of the experimental results (p1q0 and p1q1) and (2) the chain length and composition have a negligible effect on the rate constants, which is consistent with the experimental work. The interacting quantum atoms energy decomposition scheme is used to analyze the observed trends and extract chemical information related to the imaginary frequencies and barrier heights that are key parameters controlling the reactivity of the reaction. The intramolecular exchange-correlation contributions in the reactants and transition states were found to be the dominating factor.

DFT Study of Photochemical Properties and Radiative Forcing Efficiency Features of the Stereoisomers cis- and trans-CHCl═CH-CF3

The accurate assessment of radiative forcing efficiency (RFE) of a greenhouse gas is based on the precise knowledge of its structure and infrared absorption spectrum. The present work investigates the UV-vis absorption spectra and IR absorption spectra that are used for the determination of RFE of the short-lived compounds, cis- and trans-CHCl═CHCF₃ (CTFP). These investigations were carried out with six different density functional theory (DFT) methods B3LYP, CAM-B3LYP, M06, M06-2X, TPSS0, and ωB97X-D associated to the basis set 6-31G(3df). Therefrom, the relative populations of the two states cis and trans for temperatures over the range 220-370 K at 1 atm and along the atmospheric altitude were assessed. It turns out that trans-CTFP is the abundant component between the two states. This review reveals that B3LYP and M06 reproduce well the experimental results of UV-vis spectra of trans-CTFP. As for cis-CTFP, EOM-CCSD is not well fitted by DFT methods. The cis- to trans-CTFP isomerization leads to the red shift for DFT methods and to the blue shift in regard to EOM-CCSD and experimental results. The IR absorption spectra are well fitted by B3LYP over the range 500-1600 cm^-1 and TPSS0 over 1300-2000 cm^-1 for both stereoisomers. Moreover, the root-mean-square errors (RMSEs) of frequencies from experimental data are lower for B3LYP and TPSS0 for both systems. The computed IR absorption band strengths over 500-2000 cm^-1 for cis- and over 600-1800 cm^-1 for trans-CTFP are consistent with the experiment. The relevant descriptor RFEs of the climate effect were calculated using a narrow band model for a constant vertical profile and then corrected with a lifetime factor for different computational methods. The computed values correlate well with the experimental results for both stereoisomers except M06-2X and TPSS0. It is worth noting that, for both systems, the intense radiative forcing spectra are located at frequencies ranging in 1000-1200 cm^-1. The lower forcings of trans-CTFP lying in the atmospheric window region 800-1000 cm^-1 are greater than those of cis-CTFP. Therefore, RFE(trans-CTFP) = 1.127 RFE(cis-CTFP).