Initiation rate for shock-heated hydrogen-oxygen-carbon monoxide-argon mixtures as determined by OH induction time measurements

Reaction energy benchmarks of hydrocarbon combustion by Gaussian basis and plane wave basis approaches

The reaction energies of 275 elementary reactions from the hydrocarbon combustion model GRI-Mech 3.0 were evaluated by electronic structure calculations using both localized Gaussian basis and plane wave basis sets. In the Gaussian basis calculations, the d-polarization function on C, N, and O elements reduces the mean absolute deviation (MAD) from the experimental value by 53%, a significant improvement in computational accuracy. In the plane wave basis calculation using different exchange-correlation (XC) functionals, the MAD values were 0.316-0.426 eV when non-hybrid type XC functionals such as RPBE, PBE, PW91, revPBE, and PBEsol were used. On the other hand, hybrid functionals like B3LYP and HSE06 reduced the MAD values significantly down to 0.182 and 0.233 eV, respectively. The B3LYP results have 49% less MAD compared to the PBE results. These demonstrated the strong advantage of the hybrid functional for calculating gas-phase reaction energies. The present comprehensive benchmarks will be crucial for future microkinetics as well as machine learning studies on the catalytic reactions. © 2019 Wiley Periodicals, Inc.

Theoretical study of the reaction of carbon monoxide with oxygen molecules in the ground triplet and singlet delta states

Quantum chemical calculations were carried out to study the reaction of carbon monoxide with molecular oxygen in the ground triplet and singlet delta states. Transition states and intermediates that connect the reactants with products of the reaction on the triplet and singlet potential energy surfaces were identified on the base of coupled-cluster method. The values of energy barriers were refined by using compound techniques such as CBS-Q, CBS-QB3, and G3. The calculations showed that there exists an intersection of triplet and singlet potential energy surfaces. This fact leads to the appearance of two channels for the triplet CO+O(2)(X(3)Σ(g)(-)) reaction, which produces atomic oxygen in the ground O((3)P) and excited O((1)D) states. The appropriate rate constants of all reaction paths were estimated on the base of nonvariational transition-state theory. It was found that the singlet reaction rate constant is much greater than the triplet one and that the reaction channel CO+O(2)(a(1)Δ(g)) should be taken into consideration to interpret the experimental data on the oxidation of CO by molecular oxygen.

Oxidation of CO by SO2: A Theoretical Study

The elementary reaction SO(2) + CO --> CO(2) + SO((3)Sigma) (1) and the subsequent reaction SO((3)Sigma) + CO --> CO(2) + S((3)P) (2) have been studied by the application of the Gaussian-3//B3LYP quantum chemical approach to characterize the potential energy surfaces and transition state kinetic analysis to derive rate coefficients. Reaction 1 is found to take place via two transition states (TS), a cis-OSOCO TS and a trans-OSOCO TS. Reaction via the cis-TS is concerted and takes place on a singlet surface. Intersystem crossing to the final products occurs after passage through the barrier on the singlet surface. The trans-TS leads to a very weakly bound singlet OSOCO intermediate that then passes through a second TS (on the triplet surface) to form the products. Reaction 2 takes place on triplet surfaces. There is a concerted reaction through a cis-SOCO TS and a weakly bound trans-SOCO has also been identified. Reaction 2 is analogous to the reaction CO + O(2)((3)Sigma) --> CO(2) + O((3)P) (3), and this reaction has been reinvestigated at a similar level of theory and the rate coefficient derived by quantum chemistry is compared with experiment. The sensitive effects of trace impurities such as H(2), H(2)O, and hydrocarbons on the accurate experimental determination of the rate coefficient of reaction 3 is discussed. Using rate coefficients for reactions 1 and 2 obtained via quantum chemical calculations, we have been unable to model the extent of decomposition of SO(2) measured in a shock tube study of reaction between SO(2) and CO [Bauer, S. H.; Jeffers, P.; Lifshitz, A.; Yadava, B. P. Proc. Combust. Inst. 1971, 13, 417]. In light of the known sensitivity of reaction 3 to trace impurities, we have incorporated trace amounts of H(2), CH(4), or H(2)O, together with our rate coefficients for (1) and (2), in a kinetic model of Alzueta et al. [Combust. Flame 2001, 127, 2234], which is then shown to be able to substantially model the SO(2) data of Bauer et al. In the course of this modeling study we also computed heats of formation for a number of sulfur-containing small molecules: HS, HSO, HSOH, HOSO, HS(2), HSO(2), HOSO(2), HOSOH, and HOSHO.