Thermochemistry of Radicals and Molecules Relevant to Atmospheric Chemistry: Determination of Group Additivity Values using G3//B3LYP Theory

Comprehensive Theoretical Study on Four Typical Intramolecular Hydrogen Shift Reactions of Peroxy Radicals: Multireference Character, Recommended Model Chemistry, and Kinetics

Intramolecular hydrogen shift reactions in peroxy radicals (RO₂^• → ^•QOOH) play key roles in the low-temperature combustion and in the atmospheric chemistry. In the present study, we found that a mild-to-moderate multireference character of a potential energy surface (PES) is widely present in four typical hydrogen shift reactions of peroxy radicals (RO₂^•, R = ethyl, vinyl, formyl methyl, and acetyl) by a systematic assessment based on the T₁ diagnostic, %TAE diagnostic, M diagnostic, and contribution of the dominant configuration of the reference CASSCF wavefunction (C₀²). To assess the effects of these inherent multireference characters on electronic structure calculations, we compared the PESs of the four reactions calculated by the multireference method CASPT2 in the complete basis set (CBS) limit, single-reference method CCSD(T)-F12, and single-reference-based composite method WMS. The results showed that ignoring the multireference character will introduce a mean unsigned deviation (MUD) of 0.46-1.72 kcal/mol from CASPT2/CBS results by using the CCSD(T)-F12 method or a MUD of 0.49-1.37 kcal/mol by WMS for three RO₂^• reactions (R = vinyl, formyl methyl, and acetyl) with a stronger multireference character. Further tests by single-reference Kohn-Sham (KS) density functional theory methods showed even larger deviations. Therefore, we specifically developed a new hybrid meta-generalized gradient approximation (GGA) functional M06-HS for the four typical H-shift reactions of peroxy radicals based on the WMS results for the ethyl peroxy radical reaction and on the CASPT2/CBS results for the others. The M06-HS method has an averaged MUD of 0.34 kcal/mol over five tested basis sets against the benchmark PESs, performing best in the tested 38 KS functionals. Last, in a temperature range of 200-3000 K, with the new functional, we calculated the high-pressure-limit rate coefficients of these H-shift reactions by the multi-structural variational transition-state theory with the small-curvature tunneling approximation (MS-CVT/SCT) and the thermochemical properties of all of the involved key radicals by the multi-structural torsional (MS-T) anharmonicity approximation method.

Thermochemistry of Oxygen-Containing Organosilane Radicals and Uncertainty Estimations of Organosilane Group-Additivity Values

Si-C-H-O-containing radicals are important intermediates during the combustion and pyrolysis of precursors applied for the gas-phase synthesis of silica nanoparticles. Despite the industrial importance of silica nanoparticles, a comprehensive thermodynamics database of organosilane species is still missing. This work presents thermochemical data of 91 Si-C-H-O radical species. Quantum-chemical calculations and isodesmic reaction schemes are used to determine the standard enthalpy of formation, entropy, and heat capacities covering the 298-2000 K temperature range. In addition, 90 group-additivity values (GAVs) are calculated, which cover all relevant group increments. A combinatorial approach is used to ensure that all possible group increments are considered. The theoretically calculated species are used as a training set to derive 90 GAVs of Si-C-H-O radical species for the first time. In addition, uncertainty contributions of GAVs were estimated. These uncertainty estimates also comprise GAVs that were previously derived to compute thermochemical data of stable Si-C-H species and radicals as well as stable Si-C-H-O compounds. Therefore, uncertainty contributions of GAVs for a whole set of 243 group increments used to predict thermochemical data of Si-organic species are reported.

Low Temperature Oxidation of Cyclohexane: Uncertainty of Important Thermo-Chemical Properties

The study of the standard formation enthalpy, entropy, and heat capacity for key species relevant to the low-temperature combustion of cyclohexane has been performed by applying the group additivity method of Benson. The properties of 18 Benson groups (8 of them for the first time), and 10 ring correction factors for cyclic species were estimated through different empirical and semi-empirical methods. The method validation proceeded through comparison of predicted values for certain number of newly estimated groups and available literature data derived from quantum chemistry estimations. Further validations of the estimated properties of groups have been provided by comparing estimated properties of test species with data in literature and kinetic databases. Also the standard deviation between prediction and reported values has been evaluated for each validation case. A similar approach has been applied for validation of the estimated ring correction groups. For selected well-studied cyclic molecules the predicted values and the literature data have been compared with each other, and the standard deviations have been also reported. The evaluated properties of the cyclohexane relevant species were also compared with similar ones available in other kinetic models and in databases. At the end the estimated properties have been presented in a tabulated form of NASA polynomial coefficients with extrapolation up to 3500 K.