Accurate Calculation of the Energy Barriers and Rate Constants of the Large-size Molecular Reaction System for Abstraction from Alkyl Hydroperoxides

Chen Fangfang,Sun Xiaohui,Yao Qian,Li Zerong,Wang Jingbo,Li Xiangyuan
DOI: https://doi.org/10.6023/a18010015
2018-01-01
Acta Chimica Sinica
Abstract:The reaction class of a free radical with a molecule are non-elementary reactions with negative activation energies and they are usually proceeded through two reaction steps with the first step being a reactant complex formation. This class of reactions are widespread in the atmospheric chemistry and the mechanism of hydrocarbon fuel combustion, so they are extensively studied in the theoretical calculation and experimental studies. The reaction class of alpha-H abstraction from alkyl hydroper-oxides (ROOH) by hydroxyl radicals, which are important in the mechanism of hydrocarbon fuel combustion, are chosen as the object of this study. The regularity of this reaction class are revealed by quantum chemical calculations and their kinetic parameters are accurately calculated. When the standard molar Gibbs free energy change of the formation of the reactant complex in the first step is equal to zero, the corresponding temperature is defined as the conversion temperature T-c in this study, and it is shown that a steady state approximation method are applicable for this kind of reaction system to calculate the overall reaction rate constants when the temperature is much higher than the T-c. Geometric optimization and frequency analysis for all species were conducted at the BHandHLYP/6-311G(d,p) level. Five reactions are chosen as the representative for the reaction class and their single point energies are calculated using the method of CCSD(T)/CBS and it is shown that the highest conversion temperature for the five reactions is 195.17 K, far below usual modeling lowest temperature of the hydrocarbon fuel combustion, and therefore, the steady state approximation method is reasonable. It is also shown that the reaction-center geometries of the transition states are conserved, and thus the isodesmic reaction method is applicable to this reaction class to correct the energy barriers and rate constants at low-level BHandHLYP method. The obtained energy barriers are compared with the results using high-level ab initio CCSD(T)/CBS method and it is shown that the maximum absolute deviation of reaction energy barriers can be reduced from 19.99 kJ.mol(-1) before correction to 1.47 kJ.mol(-1) after correction, indicating that the isodesmic reaction method are applicable for the accurate calculation of the kinetic parameters for large-size molecular systems with the negative activation energy reaction. Finally, energy barriers for 20 reactions in the class are calculated with the isodesmic reaction method, and then based on steady state approximation, the rate constants for the overall reactions are calculated using the transition state theory in combination with the isodesmic correction scheme. It is shown that the negative activation energy relationship for the reaction class only exists in the low temperature region.
What problem does this paper attempt to address?