Global Master Equation Analysis of Rate Data for the Reaction C2H4 + H ⇄ C2H5: ΔfH0C2H5

J Phys Chem A. 2021 Nov 4;125(43):9548-9565. doi: 10.1021/acs.jpca.1c05911. Epub 2021 Oct 27.

Abstract

While forward and reverse rate constants are frequently used to determine enthalpies of reaction and formation, this process is more difficult for pressure-dependent association/dissociation reactions, especially since the forward and reverse reactions are usually studied at very different temperatures. The problems can be overcome by using a data-fitting procedure based on a master equation model. This approach has been applied to existing experimental pressure-dependent forward and reverse rate coefficients for the reaction C2H4 + H ⇄ C2H5 (k1, k-1) using the MESMER code to determine ΔfH0C2H5 from the enthalpy of the reaction. New measurements of k1, k-1 were included in analysis. They are based on laser flash photolysis with direct observation of H atom time profiles by vacuum ultraviolet laser-induced fluorescence under conditions where the approach to equilibrium could be observed. Measurements were made over the temperature range 798-828 K and with [He] from 2.33 to 7.21 × 1018 molecule cm-3. These data were then combined with a wide range of existing experimental data with helium as the bath gas (112 measurements of k1 and k-1, covering the temperature range 285-1094 K, and [He] = 7.1 × 1015-1.9 × 1019 molecule cm-3) and fitted using the master equation solver MESMER. The required vibrational frequencies and rotational constants of the system were obtained from ab initio calculations, and the activation threshold for association (ΔEthresh), enthalpy of reaction (ΔrH0), imaginary frequency (υimag), and helium energy-transfer parameters (⟨ΔEd,298(T/298)n) were optimized. The resulting parameters (errors are 2σ) are ΔEthresh = 11.43 ± 0.34 kJ mol-1, ΔrH0 = -145.34 ± 0.60 kJ mol-1, υimag = 730 ± 130 cm-1, ⟨ΔEd,298 = 54.2 ± 7.6 cm-1, and n = 1.17 ± 0.12. A value of ΔfH298.15(C2H5) = 120.49 ± 0.57 kJ mol-1 is obtained by combining ΔrH0 with standard enthalpies of formation for H and C2H4 and making the appropriate temperature corrections. The dependence of these parameters on how the internal rotor and CH2 inversion modes are treated has been explored. The experimental data for other bath gases have been analyzed, and data sets compatible with the potential energy surface parameters determined above have been identified. The parameters are virtually identical but with slightly smaller error limits. Parameterization of k1, k-1 using the Troe formalization has been used to investigate competition between ethyl decomposition and reaction with oxygen under combustion conditions.