The dissociation/recombination reaction CH4 (+M) ⇔ CH3 + H (+M): A case study for unimolecular rate theory

Reconciling Accuracy and Feasibility for Barrierless Reaction Steps by the PCS/DDCI/MC-PDFT Protocol: Methane and Ethylene Dissociations as Case Studies

Several enhancements have been introduced into state-of-the-art computational protocols for the treatment of barrierless reaction steps in the framework of variable reaction coordinate variational transition state theory. The first step is the synergistic integration of the Iterative Difference Dedicated Configuration Interaction (I-DDCI) and Pisa Composite Scheme, which defines a reduced cost, yet very accurate, computational workflow. This approach provides a near black box tool for obtaining 1D reference potentials. Then, a general strategy has been devised for tuning the level of theory used in Monte Carlo (MC) sampling, employing Multiconfiguration Pair Density Functional Theory (MC-PDFT) with dynamically adjusted Hartree-Fock exchange. Concurrently, partial geometry optimizations during the MC simulations account for the coupling between the reaction coordinates and conserved modes. The protocol closely approaches full size consistency and yields highly accurate results, with several test computations suggesting rapid convergence of the I-DDCI correction with the basis set dimensions. The capabilities of the new platform are illustrated by two case studies (the hydrogen dissociation from CH4 and C2H4), which highlight its flexibility in handling different carbon hybridizations (sp3 and sp2). The remarkable accuracy of the computed rate constants confirms the robustness of the proposed method. Together with their intrinsic interest, these results pave the way for systematic investigations of complex gas-phase reactions through a reliable, user-friendly tool accessible to specialists and nonspecialists alike.

Quantum versus classical unimolecular fragmentation rate constants and activation energies at finite temperature from direct dynamics simulations

In the present work, we investigate how nuclear quantum effects modify the temperature dependent rate constants and, consequently, the activation energies in unimolecular reactions. In the reactions under study, nuclear quantum effects mainly stem from the presence of a large zero point energy. Thus, we investigate the behavior of methods compatible with direct dynamics simulations, the quantum thermal bath (QTB) and ring polymer molecular dynamics (RPMD). To this end, we first compare them with quantum reaction theory for a model Morse potential before extending this comparison to molecular models. Our results show that, in particular in the temperature range comparable with or lower than the zero point energy of the system, the RPMD method is able to correctly capture nuclear quantum effects on rate constants and activation energies. On the other hand, although the QTB provides a good description of equilibrium properties including zero-point energy effects, it largely overestimates the rate constants. The origin of the different behaviours is in the different distance distributions provided by the two methods and in particular how they differently describe the tails of such distributions. The comparison with transition state theory shows that RPMD can be used to study fragmentation of complex systems for which it may be difficult to determine the multiple reaction pathways and associated transition states.

A theoretical study of the addition of CH2OO to hydroxymethyl hydroperoxide and its implications on SO3 formation in the atmosphere

The reaction of hydroxymethyl hydroperoxide (HMHP, HOCH₂OOH) with the simplest Criegee intermediate, CH₂OO, has been examined using quantum chemical methods with transition state theory. Geometry optimization and IRC calculations were performed using the M06-2X, MN15-L, and B2PLYP-D3 functionals in conjunction with the aug-cc-pVTZ basis set. Single point energy calculations using QCISD(T) and BD(T) with the same basis set have been performed to determine the energy of reactants, reactive complexes, transition states, and products. Rate coefficients have been obtained using variational transition state theory. The addition of CH₂OO on the three different oxygen atoms in HMHP has been considered and the ether oxide forming channel, CH₂OO + HOCH₂OOH → HOCH₂O(O)CH₂OOH (channel 2), is the most favorable. The best computed standard enthalpy of reaction (ΔH) and zero-point corrected barrier height are -20.02 and -6.33 kcal mol^-1, respectively. The reaction barrier is negative and our results suggest that both the inner and outer transition states contribute to the corresponding overall reactive flux in the tropospheric temperature range (220 K to 320 K). A two-transition state model has been used to obtain reliable rate coefficients at the high-pressure limit. The pressure-dependent rate coefficient calculations using the SS-QRRK theory have shown that this channel is pressure-dependent. Moreover, our investigation has shown that the ether oxide formed may rapidly react with SO₂ at 298 K to form SO₃, which can, in turn, react with water to form atmospheric H₂SO₄. A similar calculation has been conducted for the reaction of HMHP with OH, suggesting that the titled reaction may be a significant sink of HMHP. Therefore, the reaction between CH₂OO and HOCH₂OOH could be an indirect source for generating atmospheric H₂SO₄, which is crucial to the formation of clouds, and it might relieve global warming.

High-Temperature Unimolecular Decomposition of Diethyl Ether: Shock-Tube and Theory Studies

The unimolecular decomposition of diethyl ether (DEE; C₂H₅OC₂H₅) is considered to be initiated via a molecular elimination and a C-O and a C-C bond fission step: C₂H₅OC₂H₅ → C₂H₄ + C₂H₅OH (1), C₂H₅OC₂H₅ → C₂H₅ + C₂H₅O (2), and C₂H₅OC₂H₅ → CH₃ + C₂H₅OCH₂ (3). In this work, two shock-tube facilities were used to investigate these reactions via (a) time-resolved H-atom concentration measurements by H-ARAS (atomic resonance absorption spectrometry), (b) time-resolved DEE-concentration measurements by high repetition-rate time-of-flight mass spectrometry (HRR-TOF-MS), and (c) product-composition measurements via gas chromatography/MS (GC/MS) after quenching the test gas. The experiments were conducted at temperatures ranging from 1054 to 1505 K and at pressures between 1.2 and 2.5 bar. Initial DEE mole fractions between 0.4 and 9300 ppm were used to perform the kinetics experiments by H-ARAS (0.4 ppm), GC/MS (200-500 ppm), and HRR-TOF-MS (7850-9300 ppm). The rate constants, k_total (k_total = k₁ + k₂ + k₃) derived from the GC/MS and HRR-TOF-MS experiments agree well with each other and can be described by the Arrhenius expression, k_total(1054-1467 K; 1.3-2.5 bar) = 10^12.81±0.22 exp(-240.27 ± 5.11 kJ mol^-1/RT) s^-1. From the H-ARAS experiments, overall rate constants for the bond fission channels, k₂₊₃ = k₂ + k₃ have been extracted. The k₂₊₃ data can be well described by the Arrhenius equation, k₂₊₃(1299-1505 K; 1.3-2.5 bar) = 10^14.43±0.33 exp(-283.27 ± 8.78 kJ mol^-1/RT) s^-1. A master-equation analysis was performed using CCSD(T)/aug-cc-pvtz//B3LYP/aug-cc-pvtz and CASPT2/aug-cc-pvtz//B3LYP/aug-cc-pvtz molecular properties and energies for the three primary thermal decomposition processes in DEE. The derived experimental data is very well reproduced by the simulations with the mechanism of this work. With regard to the branching ratios between bond fissions and elimination channels, uncertainties remain.

Shock Wave and Theoretical Modeling Study of the Dissociation of CH2F2. I. Primary Processes

The unimolecular dissociation of CH₂F₂ leading to CF₂ + H₂, CHF + HF, or CHF₂ + H is investigated by quantum-chemical calculations and unimolecular rate theory. Modeling of the rate constants is accompanied by shock wave experiments over the range of 1400-1800 K, monitoring the formation of CF₂. It is shown that the energetically most favorable dissociation channel leading to CF₂ + H₂ has a higher threshold energy than the energetically less favorable one leading to CHF + HF. Falloff curves of the dissociations are modeled. Under the conditions of the described experiments, the primary dissociation CH₂F₂ → CHF + HF is followed by the reaction CHF + HF → CF₂ + H₂. The experimental value of the rate constant for the latter reaction indicates that it does not proceed by an addition-elimination process involving CH₂F₂* intermediates, as assumed before.

Shock wave and modeling study of the reaction CF4 (+M) ⇔ CF3 + F (+M)

The thermal decomposition of CF4 (+Ar) → CF3 + F (+Ar) was studied in shock waves over the temperature range 2000-3000 K varying the bath gas concentration [Ar] between 4 × 10(-6) and 9 × 10(-5) mol cm(-3). It is shown that the reaction corresponds to the intermediate range of the falloff curve. By combination with room temperature data for the reverse reaction CF3 + F (+He) → CF4 (+He) and applying unimolecular rate theory, falloff curves over the temperature range 300-6000 K are modeled. A comparison with the reaction system CH4 (+M) ⇔ CH3 + H (+M) is made.

“Plug and play” full-dimensional ab initio potential energy and dipole moment surfaces and anharmonic vibrational analysis for CH4–H2O

The first full-dimensional potential energy surface of CH₄–H₂O dimer is presented, and vibrational analysis of this dimer is performed.

A theoretical study of three gas-phase reactions involving the production or loss of methane cations

Hydrocarbon ions are important species in flames, spectroscopy and the interstellar medium. Their importance is reflected in the extensive body of literature on the structure and reactivity of carbocations. However, the geometry, electronic structure and reactivity of carbocations are difficult to assess. This study aims to contribute to the current knowledge of this subject by presenting a quantum mechanics description of methane cation dissociation using multiconfigurational methods. The geometric and electronic parameters of the minimum structure were determined for three main reaction paths: the dissociation CH4(+)→ CH2(+) + H2 and the dissociation-recombination processes CH4(+)↔ CH3(+) + H. The electronic and energetic effects of these reactions were analyzed, and it was found that each reaction path has a strong dependence on the methodology used as well as a strong multiconfigurational character during dissociation. The first doublet excited states are inner-shell excited states and may correspond to the ions that are expected to be formed after electron detachment. The rate coefficient for each reaction path was determined using variational transition state theory and RRKM/master equation calculations. The major dissociation paths, with their rate coefficients at the high-pressure limit, are CH4(+)(X(~)(2)B1) → CH3(+)(A(2)A1') + H((2)S) (k∞(T) = 1.42 × 10(+14) s(-1) exp(-37.12/RT)) and CH4(+)(X(~)(2)B1) → CH2(+)(A(2)A1) + H2((2)Σg(+)) (k∞(T) = 9.18 × 10(+14) s(-1) exp(-55.77/RT)). Our findings help to explain the abundance of ions formed from CH4 in the interstellar medium and to build models of chemical evolution.

Direct measurements of rate constants for the reactions of CH3 radicals with C2H6, C2H4, and C2H2 at high temperatures

The shock tube technique has been used to study the reactions CH3 + C2H6 → C2H4 + CH4 + H (1), CH3 + C2H4 → Products + H (2), and CH3 + C2H2 → Products + H (3). Biacetyl, (CH3CO)2, was used as a clean high temperature thermal source for CH3-radicals for all the three reactions studied in this work. For reaction 1, the experiments span a T-range of 1153 K ≤ T ≤ 1297 K, at P ~ 0.4 bar. The experiments on reaction 2 cover a T-range of 1176 K ≤ T ≤ 1366 K, at P ~ 1.0 bar, and those on reaction 3 a T-range of 1127 K ≤ T ≤ 1346 K, at P ~ 1.0 bar. Reflected shock tube experiments performed on reactions 1-3, monitored the formation of H-atoms with H-atom Atomic Resonance Absorption Spectrometric (ARAS). Fits to the H-atom temporal profiles using an assembled kinetics model were used to make determinations for k1, k2, and k3. In the case of C2H6, the measurements of [H]-atoms were used to derive direct high-temperature rate constants, k1, that can be represented by the Arrhenius equation k1(T) = 5.41 × 10(-12) exp(-6043 K/T) cm(3) molecules(-1) s(-1) (1153 K ≤ T ≤ 1297 K) for the only bimolecular process that occurs, H-atom abstraction. TST calculations based on ab initio properties calculated at the CCSD(T)/CBS//M06-2X/cc-pVTZ level of theory show excellent agreement, within ±20%, of the measured rate constants. For the reaction of CH3 with C2H4, the present rate constant results, k2', refer to the sum of rate constants, k(2b) + k(2c), from two competing processes, addition-elimination, and the direct abstraction CH3 + C2H4 → C3H6 + H (2b) and CH3 + C2H4 → C2H2 + H + CH4 (2c). Experimental rate constants for k2' can be represented by the Arrhenius equation k2'(T) = 2.18 × 10(-10) exp(-11830 K/T) cm(3) molecules(-1) s(-1) (1176 K ≤ T ≤ 1366 K). The present results are in excellent agreement with recent theoretical predictions. The present study provides the only direct measurement for the high-temperature rate constants for these channels. Lastly, measurements of H-atoms from the reaction of CH3 with C2H2 provided direct unambiguous determinations of the rate constant for the dominant process under the present experimental conditions, the addition-elimination, CH3 + C2H2 → p-C3H4 + H (3b). Experimental rate constants for k(3b) can be represented by the Arrhenius equation k(3b)(T) = 5.16 × 10(-13) exp(-3852 K/T) cm(3) molecules(-1) s(-1) (1127 K ≤ T ≤ 1346 K). The present determinations for k(3b) represent the only direct measurements for this reaction and are also in good agreement with recent theoretical predictions. The present experimental k(3b) values were also used to derive rate constants, k(-3b), for the more extensively studied back-process, the reaction of H-atoms with propyne. The best fit Arrhenius equation, combining the presently derived k(-3b) values with a recent experimental determination for k(-3b), can be represented by k(-3b)(T) = 3.87 × 10(-11) exp(-1313 K/T) cm(3) molecules(-1) s(-1) (870 K ≤ T ≤ 1346 K). The present studies represent a novel implementation of the sensitive H-ARAS technique to measure rate constants for poorly characterized and difficult to isolate "slow" CH3-radical reactions with stable C2 hydrocarbons.