Modeling the thermochemistry of nitrogen-containing compounds via group additivity.

First-principles based kinetic modeling is essential to gain insight into the governing chemistry of nitrogen-containing compounds over a wide range of technologically important processes, e.g. pyrolysis, oxidation and combustion. It also enables the development of predictive, fundamental models key to improving understanding of the influence of nitrogen-containing compounds present as impurities or process additives, considering safety, operability and quality of the product streams. A prerequisite for the generation of detailed fundamental kinetic models is the availability of accurate thermodynamic properties. To address the scarcity of thermodynamic properties for nitrogen-containing compounds, a consistent set of 91 group additive values and three non-nearest-neighbor interactions has been determined from a dataset of CBS-QB3 calculations for 300 species, including 104 radicals. This dataset contains a wide range of nitrogen-containing functionalities, i.e. imine, nitrile, nitro, nitroso, nitrite, nitrate and azo functional groups. The group additivity model enables the approximation of the standard enthalpy of formation and standard entropy at 298 K as well as the standard heat capacities over a large temperature range, i.e. 300-1500 K. For a test set of 27 nitrogen-containing compounds, the group additivity model succeeds in approximating the ab initio calculated values for the standard enthalpy of formation with a MAD of 2.3 kJ mol-1. The MAD for the standard entropy and heat capacity is lower than 4 and 2 J mol-1 K-1, respectively. For a test set of 11 nitrogen-containing compounds, the MAD between experimental and group additivity approximated values for the standard enthalpy of formation amounts to 2.8 kJ mol-1.

Modeling the kinetics of hydrogen abstraction reactions in nitrogen-containing compounds via group additivity.

New group additivity values are presented to enable the modeling of a broad range of intermolecular hydrogen abstraction reactions involving nitrogen-containing compounds. From a dataset of 316 reaction rate coefficients calculated at the CBS-QB3 level of theory in the high-pressure limit, 76 group additivity values and 14 resonance corrections have been estimated. The influence of substituents on both the attacked hydrogen and attacking radical, being a carbon or nitrogen atom, has been investigated systematically. The new group additivity models can be applied to approximate the Arrhenius parameters of hydrogen abstraction reactions of nitrogen-containing compounds by hydrogen atoms, carbon-centered and nitrogen-centered radicals in the 300-1800 K temperature range. Complementary to the group additivity model, correlations for the tunneling coefficients, which depend on both the temperature and the activation energy of the reaction in the exothermic direction, have been generated. The good performance of the new group additivity schemes has been demonstrated using a test set of reactions. At 1000 K, the rate coefficients for all test set reactions are approximated on average within a factor of 1.45, 1.47 and 1.34, for the hydrogen abstractions with a reactive center of the type H-H-N, N-H-N and C-H-N respectively.

Improved Approach for ab Initio Calculations of Rate Coefficients for Secondary Reactions in Acrylate Free-Radical Polymerization.

Secondary reactions in radical polymerization pose a challenge when creating kinetic models for predicting polymer structures. Despite the high impact of these reactions in the polymer structure, their effects are difficult to isolate and measure to produce kinetic data. To this end, we used solvation-corrected M06-2X/6-311+G(d,p) ab initio calculations to predict a complete and consistent data set of intrinsic rate coefficients of the secondary reactions in acrylate radical polymerization, including backbiting, ß-scission, radical migration, macromonomer propagation, mid-chain radical propagation, chain transfer to monomer and chain transfer to polymer. Two new approaches towards computationally predicting rate coefficients for secondary reactions are proposed: (i) explicit accounting for all possible enantiomers for reactions involving optically active centers; (ii) imposing reduced flexibility if the reaction center is in the middle of the polymer chain. The accuracy and reliability of the ab initio predictions were benchmarked against experimental data via kinetic Monte Carlo simulations under three sufficiently different experimental conditions: a high-frequency modulated polymerization process in the transient regime, a low-frequency modulated process in the sliding regime at both low and high temperatures and a degradation process in the absence of free monomers. The complete and consistent ab initio data set compiled in this work predicts a good agreement when benchmarked via kMC simulations against experimental data, which is a technique never used before for computational chemistry. The simulation results show that these two newly proposed approaches are promising for bridging the gap between experimental and computational chemistry methods in polymer reaction engineering.

Connecting Gas-Phase Computational Chemistry to Condensed Phase Kinetic Modeling: The State-of-the-Art.

In recent decades, quantum chemical calculations (QCC) have increased in accuracy, not only providing the ranking of chemical reactivities and energy barriers (e.g., for optimal selectivities) but also delivering more reliable equilibrium and (intrinsic/chemical) rate coefficients. This increased reliability of kinetic parameters is relevant to support the predictive character of kinetic modeling studies that are addressing actual concentration changes during chemical processes, taking into account competitive reactions and mixing heterogeneities. In the present contribution, guidelines are formulated on how to bridge the fields of computational chemistry and chemical kinetics. It is explained how condensed phase systems can be described based on conventional gas phase computational chemistry calculations. Case studies are included on polymerization kinetics, considering free and controlled radical polymerization, ionic polymerization, and polymer degradation. It is also illustrated how QCC can be directly linked to material properties.

Ab initio derived group additivity model for intramolecular hydrogen abstraction reactions.

A set of group additivity values for intramolecular hydrogen abstraction reactions of alkanes, alkenes and alkynes is reported. Calculating 448 reaction rate coefficients at the CBS-QB3 level of theory for 1-2 up to 1-7 hydrogen shift reactions allowed the estimation of ΔGAV° values for 270 groups. The influence of substituents on (1) the attacking radical, (2) the attacked carbon atom, and (3) the carbon chain between the attacking and attacked reactive atom has been systematically studied. Substituents have been varied between hydrogen atoms and sp3, sp2 and sp hybridized carbon atoms. It has been assumed that substituents further away from the reactive atoms or their connecting carbon chain have negligible influences on the kinetics. This group additivity model is applicable to a wide variety of reactions in the 300-1800 K temperature range. Correlations for tunneling coefficients have been generated which are complementary to the ΔGAV°'s to obtain accurate rate coefficients without the need for imaginary frequencies or electronic energies of activation. These correlations depend on the temperature and activation energy of the exothermic step. The group additivity model has been successfully applied to a test set of reactions also calculated at the CBS-QB3 level of theory. A mean absolute deviation of 1.18 to 1.71 has been achieved showing a good overall accuracy of the model.

Computational Investigation of the Aminolysis of RAFT Macromolecules.

This work presents a detailed computational study and kinetic analysis of the aminolysis of dithioates, dithiobenzoates, trithiocarbonates, xanthates, and thiocarbamates, which are frequently used as chain-transfer agents for reversible addition-fragmentation chain-transfer (RAFT) polymerization. Rate coefficients were obtained from ab initio calculations, taking into account a diffusional contribution according to the encounter pair model. A kinetic model was constructed and reveals a reaction mechanism of four elementary steps: (i) formation of a zwitterionic intermediate, (ii) formation of a complex intermediate in which an assisting amine molecule takes over the proton from the zwitterionic intermediate, (iii) breakdown of the complex into a neutral tetrahedral intermediate with release of the assisting amine molecule, and (iv) amine-assisted breakdown of the neutral intermediate to the products. Furthermore, a comparative analysis indicates that the alkanedithioates and dithiobenzoates react the fastest, followed, respectively, by xanthates and trithiocarbonates, which react almost equally fast, and dithiocarbamates, which are not reactive at typical experimentally used conditions.

Computational Study and Kinetic Analysis of the Aminolysis of Thiolactones.

The aminolysis of three differently α-substituted γ-thiolactones (C4H5OSX, X = H, NH2, and NH(CO)CH3) is modeled based on CBS-QB3 calculated free energies corrected for solvation using COSMO-RS. For the first time, quantitative kinetic and thermodynamic data are provided for the concerted path and the stepwise path over a neutral tetrahedral intermediate. These paths can take place via an unassisted, an amine-assisted, or a thiol-assisted mechanism. Amine assistance lowers the free energy barriers along both paths, while thiol assistance only lowers the formation of the neutral tetrahedral intermediate. Based on the ab initio calculated rate coefficients, a kinetic model is constructed that is able to reliably describe experimental observations for the aminolysis of N-acetyl-dl-homocysteine thiolactone with n-butylamine in THF and CHCl3. Reaction path analysis shows that for all conditions relevant for applications in polymer synthesis and postpolymer modification, an assisted stepwise mechanism is operative in which the formation of the neutral tetrahedral intermediate is rate-determining and which is mainly amine-assisted at low conversions and thiol-assisted at high conversions.

Group Additive Kinetics for Hydrogen Transfer Between Oxygenates.

Hydrogen abstraction reactions involving oxygenates in gaseous phase play an important role in many biomass-related conversion processes. In this work, group additivity is used to provide Arrhenius parameters in a temperature range of 300-2500 K for hydrogen abstractions between oxygenate compounds such as alcohols, ethers, esters, acids, ketones, diketones, aldehydes, hydroxyperoxides, alkyl peroxides, and unsaturated ethers and ketones. The group additive values for Arrhenius parameters of hydrogen transfer reactions of the type O--H--C and O--H--O are derived from CBS-QB3 calculations in the high-pressure limit. From a total set of 118 reactions, 43 group additivity values are determined. Inclusion of an additional 37 corrections accounting for cross-resonance effects in the transition state further improves the accuracy of the model. For a set of 25 ab initio calculated and 60 experimental rate coefficients, group additive modeling reproduces rate coefficients within a mean factor of deviation of ∼3. Hence, the developed group additive models can be reliably used for an accurate and fast prediction of the kinetics of hydrogen abstractions involving oxygenates.

Thermodynamic study of benzene and hydrogen coadsorption on Pd(111).

Periodic density functional theory (DFT) has been used to study the coadsorption of hydrogen and benzene on Pd(111). The most stable coverages are predicted by constructing the thermodynamic phase diagram as a function of gas-phase temperature and pressure. The common approximation that neglects vibrational contributions to the surface Gibbs free energy, using the PW91 functional, is compared to the one that includes vibrational contributions. Higher coverages are predicted to be thermodynamically the most stable including vibrational frequencies, mainly due to the different entropy contributions. The first approach is also compared to the one using a (optPBE-vdW) vdW-DF functional without vibrational contributions, which predicts higher benzene coverages for benzene adsorption, and lower hydrogen coverages for hydrogen adsorption and coadsorption with a fixed benzene coverage. Inclusion of vibrational contributions using the vdW-DF method has not been implemented due to computational constraints. However, an estimate of the expected result is proposed by adding PW91 vibrational contributions to the optPBE-vdW electronic energies, and under typical hydrogenation conditions high coverages of about θH = 0.89 are expected. Inclusion of vibrational contributions to the surface Gibbs free energy and a proper description of van der Waals interaction are recommended to predict the thermodynamically most stable surface coverage.

Kinetic modeling of α-hydrogen abstractions from unsaturated and saturated oxygenate compounds by hydrogen atoms.

Hydrogen-abstraction reactions play a significant role in thermal biomass conversion processes, as well as regular gasification, pyrolysis, or combustion. In this work, a group additivity model is constructed that allows prediction of reaction rates and Arrhenius parameters of hydrogen abstractions by hydrogen atoms from alcohols, ethers, esters, peroxides, ketones, aldehydes, acids, and diketones in a broad temperature range (300-2000 K). A training set of 60 reactions was developed with rate coefficients and Arrhenius parameters calculated by the CBS-QB3 method in the high-pressure limit with tunneling corrections using Eckart tunneling coefficients. From this set of reactions, 15 group additive values were derived for the forward and the reverse reaction, 4 referring to primary and 11 to secondary contributions. The accuracy of the model is validated upon an ab initio and an experimental validation set of 19 and 21 reaction rates, respectively, showing that reaction rates can be predicted with a mean factor of deviation of 2 for the ab initio and 3 for the experimental values. Hence, this work illustrates that the developed group additive model can be reliably applied for the accurate prediction of kinetics of α-hydrogen abstractions by hydrogen atoms from a broad range of oxygenates.

Kinetic modeling of α-hydrogen abstractions from unsaturated and saturated oxygenate compounds by carbon-centered radicals.

Hydrogen abstractions are important elementary reactions in a variety of reacting media at high temperatures in which oxygenates and hydrocarbon radicals are present. Accurate kinetic data are obtained from CBS-QB3 ab initio (AI) calculations by using conventional transition-state theory within the high-pressure limit, including corrections for hindered rotation and tunneling. From the obtained results, a group-additive (GA) model is developed that allows the Arrhenius parameters and rate coefficients for abstraction of the α-hydrogen from a wide range of oxygenate compounds to be predicted at temperatures ranging from 300 to 1500 K. From a training set of 60 hydrogen abstractions from oxygenates by carbon-centered radicals, 15 GA values (ΔGAV°s) are obtained for both the forward and reverse reactions. Among them, four ΔGAV°s refer to primary contributions, and the remaining 11 ΔGAV°s refer to secondary ones. The accuracy of the model is further improved by introducing seven corrections for cross-resonance stabilization of the transition state from an additional set of 43 reactions. The determined ΔGAV°s are validated upon a test set of AI data for 17 reactions. The mean absolute deviation of the pre-exponential factors (log A) and activation energies (E(a)) for the forward reaction at 300 K are 0.238 log(m(3) mol(-1) s(-1)) and 1.5 kJ mol(-1), respectively, whereas the mean factor of deviation <ρ> between the GA-predicted and the AI-calculated rate coefficients is 1.6. In comparison with a compilation of 33 experimental rate coefficients, the <ρ> between the GA-predicted values and these experimental values is only 2.2. Hence, the constructed GA model can be reliably used in the prediction of the kinetics of α-hydrogen-abstraction reactions between a broad range of oxygenates and oxygenate radicals.

Group additive values for the gas-phase standard enthalpy of formation, entropy and heat capacity of oxygenates.

A complete and consistent set of 60 Benson group additive values (GAVs) for oxygenate molecules and 97 GAVs for oxygenate radicals is provided, which allow to describe their standard enthalpies of formation, entropies and heat capacities. Approximately half of the GAVs for oxygenate molecules and the majority of the GAVs for oxygenate radicals have not been reported before. The values are derived from an extensive and accurate database of thermochemical data obtained by ab initio calculations at the CBS-QB3 level of theory for 202 molecules and 248 radicals. These compounds include saturated and unsaturated, α- and ß-branched, mono- and bifunctional oxygenates. Internal rotations were accounted for by using one-dimensional hindered rotor corrections. The accuracy of the database was further improved by adding bond additive corrections to the CBS-QB3 standard enthalpies of formation. Furthermore, 14 corrections for non-nearest-neighbor interactions (NNI) were introduced for molecules and 12 for radicals. The validity of the constructed group additive model was established by comparing the predicted values with both ab initio calculated values and experimental data for oxygenates and oxygenate radicals. The group additive method predicts standard enthalpies of formation, entropies, and heat capacities with chemical accuracy, respectively, within 4 kJ mol(-1) and 4 J mol(-1) K(-1) for both ab initio calculated and experimental values. As an alternative, the hydrogen bond increment (HBI) method developed by Lay et al. (T. H. Lay, J. W. Bozzelli, A. M. Dean, E. R. Ritter, J. Phys. Chem.- 1995, 99, 14514) was used to introduce 77 new HBI structures and to calculate their thermodynamic parameters (Δ(f)H°, S°, C(p)°). The GAVs reported in this work can be reliably used for the prediction of thermochemical data for large oxygenate compounds, combining rapid prediction with wide-ranging application.

Benzene adsorption on binary Pt3M alloys and surface alloys: a DFT study.

Benzene adsorption on Pt3M/Pt(111) surfaces and Pt3M(111) bulk alloys (M = Fe, Co, Ni, Cu, Pd, Ag, Au) is analyzed using density functional theory calculations on 4-layered slabs in the framework of catalyst development for aromatics hydrogenation. Segregation in the top layers was allowed for, accounting for the actual stoichiometric composition of the top layers rather than using simplified 'skin' or 'sandwich' structures. On the surfaces that do not segregate (M = Pd, Ag, Au), the preferred benzene adsorption site is the hollow Pt3-hcp(0) site. On antisegregated "Pt-skin" surfaces (M = Fe, Co, Ni, Cu, Pd), which have a top layer composed entirely of Pt, benzene prefers bridge sites with a maximized number of solute atoms M in the subsurface layers. Benzene adsorption is weaker than on pure Pt(111), by 0.1-0.5 eV on the surface alloys and by 0.6-1.0 eV on bulk alloys, except for Pt3Pd alloys, which behave similarly to pure Pt. On the fully segregated Pt3Ag and Pt3Au alloys, which have a Ag resp. Au monolayer on top, only physisorption occurs. Benzene adsorption does not change the segregation state of the catalyst. From various DOS-based catalyst descriptors, the occupied d-band center of the clean catalyst slab shows the best correlation with benzene adsorption energies, allowing the prediction of benzene adsorption energies on a range of other Pt-based bimetallic alloys.

Kinetics of α hydrogen abstractions from thiols, sulfides and thiocarbonyl compounds.

Hydrogen abstraction reactions involving organosulfur compounds play an important role in many industrial, biological and atmospheric processes. Despite their chemical relevance, little is known about their kinetics. In this work a group additivity model is developed that allows predicting the Arrhenius parameters for abstraction reactions of α hydrogen atoms from thiols, alkyl sulfides, alkyl disulfides and thiocarbonyl compounds by carbon-centered radicals at temperatures ranging from 300 to 1500 K. Rate coefficients for 102 hydrogen abstractions were obtained using conventional transition state theory within the high-pressure limit. Electronic barriers were calculated using the CBS-QB3 method and the rate coefficients were corrected for tunneling and hindered rotation about the transitional bond. Group additivity values for 46 groups are determined. To account for resonance and hyperconjugative stabilization in the transition state, 8 resonance corrections were fitted to a set of 32 reactions. The developed group additivity scheme was validated using a test set containing an additional 30 reactions. The group additivity scheme succeeds in reproducing the rate coefficients on average within a factor of 2.4 at 300 K and 1.4 at 1000 K. Mean absolute deviations of the Arrhenius parameters amount to, respectively, 2.5 kJ mol(-1) for E(a) and 0.13 for log A, both at 300 and 1000 K. This work hence illustrates that the recently developed group additivity methods for Arrhenius parameters extrapolate successfully to hetero-element containing compounds.

Modeling the gas-phase thermochemistry of organosulfur compounds.

Key to understanding the involvement of organosulfur compounds in a variety of radical chemistries, such as atmospheric chemistry, polymerization, pyrolysis, and so forth, is knowledge of their thermochemical properties. For organosulfur compounds and radicals, thermochemical data are, however, much less well documented than for hydrocarbons. The traditional recourse to the Benson group additivity method offers no solace since only a very limited number of group additivity values (GAVs) is available. In this work, CBS-QB3 calculations augmented with 1D hindered rotor corrections for 122 organosulfur compounds and 45 organosulfur radicals were used to derive 93 Benson group additivity values, 18 ring-strain corrections, 2 non-nearest-neighbor interactions, and 3 resonance corrections for standard enthalpies of formation, standard molar entropies, and heat capacities for organosulfur compounds and organosulfur radicals. The reported GAVs are consistent with previously reported GAVs for hydrocarbons and hydrocarbon radicals and include 77 contributions, among which 26 radical contributions, which, to the best of our knowledge, have not been reported before. The GAVs allow one to estimate the standard enthalpies of formation at 298 K, the standard entropies at 298 K, and standard heat capacities in the temperature range 300-1500 K for a large set of organosulfur compounds, that is, thiols, thioketons, polysulfides, alkylsulfides, thials, dithioates, and cyclic sulfur compounds. For a validation set of 26 organosulfur compounds, the mean absolute deviation between experimental and group additively modeled enthalpies of formation amounts to 1.9  kJ mol(-1). For an additional set of 14 organosulfur compounds, it was shown that the mean absolute deviations between calculated and group additively modeled standard entropies and heat capacities are restricted to 4 and 2 J mol(-1) K(-1), respectively. As an alternative to Benson GAVs, 26 new hydrogen-bond increments are reported, which can also be useful for the prediction of radical thermochemistry.

Modeling the influence of resonance stabilization on the kinetics of hydrogen abstractions.

Resonance stabilization of the transition state is one of the key factors in modeling the kinetics of hydrogen abstraction reactions between hydrocarbons. A group additive model is developed which allows the prediction of rate coefficients for bimolecular hydrogen abstraction reactions over a broad range of hydrocarbons and hydrocarbon radicals between 300 and 1300 K. Group additive values for 50 groups are determined from rate coefficients determined using the high level CBS-QB3 ab initio method, corrected for tunneling and the hindered internal rotation around the transitional bond. Resonance and hyperconjugative stabilization of the transition state is accounted for by introducing 4 corrections based on the structure of the reactive moiety of the transition state. The corrections, fitted to a set of 28 reactions, are temperature-independent and reduce the mean absolute deviation on E(a) to 0.7 kJ mol(-1) and to 0.05 for log A. Tunneling contributions are accounted for by using a fourth order polynomial in the activation energy. Final validation for 19 reactions yields a mean factor of deviation between group additive prediction and ab initio calculation of 2.4 at 300 K and 1.8 at 1000 K. In comparison with 6 experimental rate coefficients (600-719 K), the mean factor of deviation is less than 3.

Hydrogen radical additions to unsaturated hydrocarbons and the reverse beta-scission reactions: modeling of activation energies and pre-exponential factors.

The group additivity method for Arrhenius parameters is applied to hydrogen addition to alkenes and alkynes and the reverse beta-scission reactions, an important family of reactions in thermal processes based on radical chemistry. A consistent set of group additive values for 33 groups is derived to calculate the activation energy and pre-exponential factor for a broad range of hydrogen addition reactions. The group additive values are determined from CBS-QB3 ab-initio-calculated rate coefficients. A mean factor of deviation of only two between CBS-QB3 and experimental rate coefficients for seven reactions in the range 300-1000 K is found. Tunneling coefficients for these reactions were found to be significant below 400 K and a correlation accounting for tunneling is presented. Application of the obtained group additive values to predict the kinetics for a set of 11 additions and beta-scissions yields rate coefficients within a factor of 3.5 of the CBS-QB3 results except for two beta-scissions with severe steric effects. The mean factor of deviation with respect to experimental rate coefficients of 2.0 shows that the group additive method with tunneling corrections can accurately predict the kinetics and is at least as accurate as the most commonly used density functional methods. The constructed group additive model can hence be applied to predict the kinetics of hydrogen radical additions for a broad range of unsaturated compounds.

First principles based group additive values for the gas phase standard entropy and heat capacity of hydrocarbons and hydrocarbon radicals.

In this work a complete and consistent set of 95 Benson group additive values (GAVs) for standard entropies S(o) and heat capacities C(p)(o) of hydrocarbons and hydrocarbon radicals is presented. These GAVs include 46 groups, among which 25 radical groups, which, to the best of our knowledge, have not been reported before. The GAVs have been determined from a set of B3LYP/6-311G(d,p) ideal gas statistical thermodynamics values for 265 species, consistently with previously reported GAVs for standard enthalpies of formation. One-dimensional hindered rotor corrections for all internal rotations are included. The computational methodology has been compared to experimental entropies (298 K) for 39 species, with a mean absolute deviation (MAD) between experiment and calculation of 1.2 J mol(-1) K(-1), and to 46 experimental heat capacities (298 K) with a resulting MAD = 1.8 J mol(-1) K(-1). The constructed database allowed evaluation of corrections on S(o) and C(p)(o) for non-nearest-neighbor effects, which have not been determined previously. The group additive model predicts the S(o) and C(p)(o) within approximately 5 J mol(-1) K(-1) of the ab initio values for 11 of the 14 molecules of the test set, corresponding to an acceptable maximal deviation of a factor of 1.6 on the equilibrium coefficient. The obtained GAVs can be applied for the prediction of S(o) and C(p)(o) for a wide range of hydrocarbons and hydrocarbon radicals. The constructed database also allowed determination of a large set of hydrogen bond increments, which can be useful for the prediction of radical thermochemistry.

Carbon-centered radical addition and beta-scission reactions: modeling of activation energies and pre-exponential factors.

A consistent set of group additive values DeltaGAV degrees for 46 groups is derived, allowing the calculation of rate coefficients for hydrocarbon radical additions and beta-scission reactions. A database of 51 rate coefficients based on CBS-QB3 calculations with corrections for hindered internal rotation was used as training set. The results of this computational method agree well with experimentally observed rate coefficients with a mean factor of deviation of 3, as benchmarked on a set of nine reactions. The temperature dependence on the resulting DeltaGAV degrees s in the broad range of 300-1300 K is limited to +/-4.5 kJ mol(-1) on activation energies and to +/-0.4 on logA (A: pre-exponential factor) for 90 % of the groups. Validation of the DeltaGAV degrees s was performed for a test set of 13 reactions. In the absence of severe steric hindrance and resonance effects in the transition state, the rate coefficients predicted by group additivity are within a factor of 3 of the CBS-QB3 ab initio rate coefficients for more than 90 % of the reactions in the test set. It can thus be expected that in most cases the GA method performs even better than standard DFT calculations for which a deviation factor of 10 is generally considered to be acceptable.

Theoretical study of the thermodynamics and kinetics of hydrogen abstractions from hydrocarbons.

Thermochemical and kinetic data were calculated at four cost-effective levels of theory for a set consisting of five hydrogen abstraction reactions between hydrocarbons for which experimental data are available. The selection of a reliable, yet cost-effective method to study this type of reactions for a broad range of applications was done on the basis of comparison with experimental data or with results obtained from computationally demanding high level of theory calculations. For this benchmark study two composite methods (CBS-QB3 and G3B3) and two density functional theory (DFT) methods, MPW1PW91/6-311G(2d,d,p) and BMK/6-311G(2d,d,p), were selected. All four methods succeeded well in describing the thermochemical properties of the five studied hydrogen abstraction reactions. High-level Weizmann-1 (W1) calculations indicated that CBS-QB3 succeeds in predicting the most accurate reaction barrier for the hydrogen abstraction of methane by methyl but tends to underestimate the reaction barriers for reactions where spin contamination is observed in the transition state. Experimental rate coefficients were most accurately predicted with CBS-QB3. Therefore, CBS-QB3 was selected to investigate the influence of both the 1D hindered internal rotor treatment about the forming bond (1D-HR) and tunneling on the rate coefficients for a set of 21 hydrogen abstraction reactions. Three zero curvature tunneling (ZCT) methods were evaluated (Wigner, Skodje & Truhlar, Eckart). As the computationally more demanding centrifugal dominant small curvature semiclassical (CD-SCS) tunneling method did not yield significantly better agreement with experiment compared to the ZCT methods, CD-SCS tunneling contributions were only assessed for the hydrogen abstractions by methyl from methane and ethane. The best agreement with experimental rate coefficients was found when Eckart tunneling and 1D-HR corrections were applied. A mean deviation of a factor 6 on the rate coefficients is found for the complete set of 21 reactions at temperatures ranging from 298 to 1000 K. Tunneling corrections play a critical role in obtaining accurate rate coefficients, especially at lower temperatures, whereas the hindered rotor treatment only improves the agreement with experiment in the high-temperature range.