Benchmarking Compound Methods (CBS-QB3, CBS-APNO, G3, G4, W1BD) against the Active Thermochemical Tables: Formation Enthalpies of Radicals

Basis Set Effects in the Description of the Cl-O Bond in ClO and XClO/ClOX Isomers (X = H, O, and Cl) Using DFT and CCSD(T) Methods

The performance of a group of density functional methods of progressive complexity for the description of the ClO bond in a series of chlorine oxides was investigated. The simplest ClO radical species and the two isomeric structures XClO/ClOX for each X = H, Cl, and O were studied using the PW91, TPSS, B3LYP, PBE0, M06, M06-2X, BMK, and B2PLYP functionals. Geometry optimizations and reaction enthalpies and enthalpies of formation for each species were calculated using Pople basis sets and the (aug)-cc-pVnZ Dunning sets, with n = D, T, Q, 5, and 6. For the calculation of enthalpies of formation, atomization and isodesmic reactions were employed. Both the precision of the methods with respect to the increase of the basis sets, as well as their accuracy, were gauged by comparing the results with the more accurate CCSD(T) calculations, performed using the same basis sets as for the DFT methods. The results obtained employing composite chemical methods (G4, CBS-QB3, and W1BD) were also used for the comparisons, as well as the experimental results when they are available. The results obtained show that error compensation is the key for successful description of molecular properties (geometries and energies) by carefully selecting the method and basis sets. In general, expansion of the one-electron basis set to the limit of completeness does not improve results at the DFT level, but just the opposite. The enthalpies of formation calculated at the CCSD(T)/aug-cc-pV6Z for the species considered are generally in agreement with experimental determinations and the most accurate theoretical values. Different sources of error in the calculations are discussed in detail.

Mechanism and kinetics of the oxidation of dimethyl carbonate by hydroxyl radical in the atmosphere

The mechanism and kinetics for the reaction of dimethyl carbonate (DMC) with OH radical have been studied by using quantum chemical methods. Four reaction pathways were identified for the initial reaction. In the first two pathways, hydrogen atom abstraction is taking place and alkyl radical intermediate is formed with the energy barrier of 6.4 and 7.9 kcal/mol. In the third pathway, OH addition reaction to the carbonyl carbon (C2) atom of DMC and intermediate, I2, is formed with an energy barrier of 11.9 kcal/mol. In the fourth pathway, along with CH₃O^●, methyl hydrogen carbonate is formed. For this C-O bond breaking and O-H addition reaction, the energy barrier is 27 kcal/mol. The calculated enthalpy and Gibbs energy values show that the studied initial reactions are exothermic and exoergic except the OH addition reaction. For the initial reactions, the rate constants were calculated by using canonical variational transition state theory (CVT) with small curvature tunneling (SCT) correction over the temperature range of 278-1200 K. At 298 K, the calculated rate coefficient for the in-plane and out-of-plane hydrogen atom abstraction reaction pathway is 2.30 × 10^-13 and 0.02 × 10^-13 cm³ molecule^-1 s^-1. Further, the reaction between alkyl radical intermediate formed from the first pathway and O₂ is studied. The reaction of alkyl peroxy radical intermediate with atmospheric oxidants, HO₂, NO, and NO₂ is also studied. It was found that the formic (methyl carbonic) anhydride is the end product formed from the atmospheric oxidation and secondary reactions of DMC.

Theoretical Investigation on the Mechanism and Kinetics of Atmospheric Reaction of Methyldichloroacetate with Hydroxyl Radical

The atmospheric reaction of methyldichloroacetate (MDCA) with OH radical is studied using electronic structure calculations. Five different pathways were considered for the initial reactions, which results in the formation of alkyl radical of MDCA along with H₂O, HOCl, and CH₃O^•. Among the five pathways studied, the α-carbon atom (-CHCl₂ site) H atom abstraction reaction, which leads to the formation of the alkyl radical intermediate ^•CCl₂C(O)OCH₃ (I1) is found to be more favorable with an energy barrier of 7.3 kcal/mol, and Cl-atom abstraction reaction is having high energy barrier of 21.3 kcal/mol at M06-2X/6-311++G(2df,2p) level. The calculated thermochemical parameters show that except Cl-atom abstraction channel the other initial reaction channels are highly exothermic. The rate constant is calculated for the initial H atom abstraction reactions using canonical variational transition state theory over the temperature range of 278 to 350 K. The Arrhenius plot shows positive temperature dependence for both the reactions. The results from the calculated thermochemical parameters and rate constants show that the formation of the alkyl radical intermediate (I1) is more favorable with the rate constant of 2.07 × 10^-13 cm³ molecule^-1 s^-1 at 298 K. The calculated atmospheric lifetime of MDCA is 28 days at normal atmospheric OH concentration. The results obtained from secondary reactions show that the major product formed from the oxidation chemistry of MDCA is methyl-2-chloro-2-oxoacetate (or) methyl oxalyl chloride.

Extensive Theoretical Study of the Thermochemical Properties of Unsaturated Hydrocarbons and Allylic and Super-Allylic Radicals: The Development and Optimization of Group Additivity Values

In this study, the thermochemistry of C₂-C₇ unsaturated hydrocarbons (22 alkene and 6 diene molecules) and 16 allylic and 5 super-allylic radicals is determined using high-accuracy quantum chemistry calculations. In addition, the group additivity values (GAVs) of a total of 19 relevant groups are systematically optimized on the basis of the calculated thermochemistry of species clusters. The M06-2X method using the 6-311++G(d,p) basis set is used for the geometry optimizations, vibrational frequency calculations, and internal rotation scans for lower-frequency modes. The composite compound methods, CBS-APNO, G3, and G4, are utilized to derive the average atomization formation enthalpies. The entropy and temperature-dependent heat capacity values of all species are calculated using statistical thermodynamics in MultiWell. These results are in good agreement with literature data. A GAVs optimization is performed on the basis of a statistical analysis: a Bland-Altman plot, which is employed to visualize the agreement between the results from the quantum chemical calculations and the GA method. It is found that the 298 K entropies of the CD/C2, C/CD2/H2, C/C/CD2/H, and C/CD3/H groups disagree by more than 5 cal K^-1 mol^-1 compared to existing values, while the values for the ALLYLS and ALLYLT radical groups also differ by ∼2.4 and 4.1 cal K^-1 mol^-1, respectively. The 298 K formation enthalpies of the C/CD2/H2, C/C/CD2/H, C/CD3/H, and ALLYLT groups are modified by more than 1 kcal mol^-1, compared to existing values. The updated GAVs can be used with increased confidence to estimate the thermochemical properties of combustion-relevant unsaturated hydrocarbon molecules and their radicals which are critical for the development of accurate chemical kinetic models describing the pyrolysis and oxidation of hydrocarbon and oxygenated hydrocarbon fuels.

Ground-State Gas-Phase Structures of Inorganic Molecules Predicted by Density Functional Theory Methods

We tested a battery of density functional theory (DFT) methods ranging from generalized gradient approximation (GGA) via meta-GGA to hybrid meta-GGA schemes as well as Møller-Plesset perturbation theory of the second order and a single and double excitation coupled-cluster (CCSD) theory for their ability to reproduce accurate gas-phase structures of di- and triatomic molecules derived from microwave spectroscopy. We obtained the most accurate molecular structures using the hybrid and hybrid meta-GGA approximations with B3PW91, APF, TPSSh, mPW1PW91, PBE0, mPW1PBE, B972, and B98 functionals, resulting in lowest errors. We recommend using these methods to predict accurate three-dimensional structures of inorganic molecules when intramolecular dispersion interactions play an insignificant role. The structures that the CCSD method predicts are of similar quality although at considerably larger computational cost. The structures that GGA and meta-GGA schemes predict are less accurate with the largest absolute errors detected with BLYP and M11-L, suggesting that these methods should not be used if accurate three-dimensional molecular structures are required. Because of numerical problems related to the integration of the exchange-correlation part of the functional and large scattering of errors, most of the Minnesota models tested, particularly MN12-L, M11, M06-L, SOGGA11, and VSXC, are also not recommended for geometry optimization. When maintaining a low computational budget is essential, the nonseparable gradient functional N12 might work within an acceptable range of error. As expected, the DFT-D3 dispersion correction had a negligible effect on the internuclear distances when combined with the functionals tested on nonweakly bonded di- and triatomic inorganic molecules. By contrast, the dispersion correction for the APF-D functional has been found to shorten the bonds significantly, up to 0.064 Å (AgI), in Ag halides, BaO, BaS, BaF, BaCl, Cu halides, and Li and Na halides and hydrides. These results do not agree well with very accurate structures derived from microwave spectroscopy; we therefore believe that the dispersion correction in the APF-D method should be reconsidered. Finally, we found that inaccurate structures can easily lead to errors of few kcal/mol in single-point energies.

Theoretical Kinetics Analysis for Ḣ Atom Addition to 1,3-Butadiene and Related Reactions on the Ċ4H7 Potential Energy Surface

The oxidation chemistry of the simplest conjugated hydrocarbon, 1,3-butadiene, can provide a first step in understanding the role of polyunsaturated hydrocarbons in combustion and, in particular, an understanding of their contribution toward soot formation. On the basis of our previous work on propene and the butene isomers (1-, 2-, and isobutene), it was found that the reaction kinetics of Ḣ-atom addition to the C═C double bond plays a significant role in fuel consumption kinetics and influences the predictions of high-temperature ignition delay times, product species concentrations, and flame speed measurements. In this study, the rate constants and thermodynamic properties for Ḣ-atom addition to 1,3-butadiene and related reactions on the Ċ₄H₇ potential energy surface have been calculated using two different series of quantum chemical methods and two different kinetic codes. Excellent agreement is obtained between the two different kinetics codes. The calculated results including zero-point energies, single-point energies, rate constants, barrier heights, and thermochemistry are systematically compared among the two quantum chemical methods. 1-Methylallyl (Ċ₄H₇1-3) and 3-buten-1-yl (Ċ₄H₇1-4) radicals and C₂H₄ + Ċ₂H₃ are found to be the most important channels and reactivity-promoting products, respectively. We calculated that terminal addition is dominant (>80%) compared to internal Ḣ-atom addition at all temperatures in the range 298-2000 K. However, this dominance decreases with increasing temperature. The calculated rate constants for the bimolecular reaction C₄H₆ + Ḣ → products and C₂H₄ + Ċ₂H₃ → products are in excellent agreement with both experimental and theoretical results from the literature. For selected C₄ species, the calculated thermochemical values are also in good agreement with literature data. In addition, the rate constants for H atom abstraction by Ḣ atoms have also been calculated, and it is found that abstraction from the central carbon atoms is the dominant channel (>70%) at temperatures in the range of 298-2000 K. Finally, by incorporating our calculated rate constants for both Ḣ atom addition and abstraction into our recently developed 1,3-butadiene model, we show that laminar flame speed predictions are significantly improved, emphasizing the value of this study.

Modeling Nitrogen Species as Pollutants: Thermochemical Influences

To simulate emissions of nitrogen-containing compounds in practical combustion environments, it is necessary to have accurate values for their thermochemical parameters, as well as accurate kinetic values to describe the rates of their formation and decomposition. Significant disparity is observed in the literature for the former, and we therefore present herein high-accuracy ab initio gas-phase thermochemistry for 60 nitrogenous compounds, many of which are important in the formation and consumption chemistry of NOx species. Several quantum-chemical composite methods (CBS-APNO, G3, and G4) were utilized to derive enthalpies of formation via the atomization method. Entropies and heat capacities were calculated from traditional statistical thermodynamics, with oscillators treated as anharmonic based on ro-vibrational property analyses carried out at the B3LYP/cc-pVTZ level of theory. The use of quantum chemical methods, along with the treatments of anharmonicities and hindered rotors, ensures accurate enthalpy of formation, entropy, and heat capacity values across the temperature range 298.15-3000 K. The implications of these results for atmospheric and combustion modeling are discussed.

A Database of Formation Enthalpies of Nitrogen Species by Compound Methods (CBS-QB3, CBS-APNO, G3, G4)

Accurate thermochemical data for compounds containing C/H/N/O are required to underpin kinetics simulation and modeling of the reactions of these species in different environments. There is a dearth of experimental data so computational quantum chemistry has stepped in to fill this breach and to verify whether particular experiments are in need of revision. A number of composite model chemistries (CBS-QB3, CBS-APNO, G3, and G4) are used to compute theoretical atomization energies and hence enthalpies of formation at 0 and 298.15 K, and these are benchmarked against the best available compendium of values, the Active Thermochemical Tables or ATcT. In general the agreement is very good for some 28 species with the only discrepancy being for hydrazine. It is shown that, although individually the methods do not perform that well, collectively the mean unsigned error is <1.7 kJ mol(-1); hence, this approach provides a useful tool to screen published values and validate new experimental results. Using multiple model chemistries does have some drawbacks but can produce good results even for challenging molecules like HOON and CN2O2. The results for these smaller validated molecules are then used as anchors for determining the formation enthalpies of larger species such as methylated hydrazines and diazenes, five- and six-membered heterocyclics via carefully chosen isodesmic working reactions with the aim of resolving some discrepancies in the literature and establishing a properly validated database. This expanded database could be useful in testing the performance of computationally less-demanding density function methods with newer functionals that have the capacity to treat much larger systems than those tested here.