M. Sabbe

First Principles Based Group Additive Values for the Gas Phase Standard Entropy and Heat Capacity of Hydrocarbons and Hydrocarbon Radicals

M. Sabbe, F. De Vleeschouwer, M-F. Reyniers, M. Waroquier, G.B. Marin
Journal of Physical Chemistry A
112 (47), 12235-12251
2008
A1

Abstract 

In this work a complete and consistent set of 95 Benson group additive values (GAVs) for standard entropies S° and heat capacities Cp° 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° and Cp° for non-nearest-neighbor effects, which have not been determined previously. The group additive model predicts the S° and Cp° within 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° and Cp° 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.

Theoretical Study of the Thermodynamics and Kinetics of Hydrogen Abstractions from Hydrocarbons

A.G. Vandeputte, M. Sabbe, M-F. Reyniers, V. Van Speybroeck, M. Waroquier, G.B. Marin
Journal of Physical Chemistry A
111 (46), 11771–11786
2007
A1

Abstract 

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.

arbon-Centered Radical Addition and β-Scission Reactions: Modeling of Activation Energies and Pre-exponential Factors

M. Sabbe, A.G. Vandeputte, M-F. Reyniers, V. Van Speybroeck, M. Waroquier, G.B. Marin
Journal of Physical Chemistry A
111 (34), 8416-8428
2007
A1

Abstract 

A consistent set of group additive values ΔGAV° for 46 groups is derived, allowing the calculation of rate coefficients for hydrocarbon radical additions and β-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 ΔGAV°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 ΔGAV°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.

Group additive values for the gas phase standard enthalpy of formation of hydrocarbons and hydrocarbon radicals

M. Sabbe, M. Saeys, M-F. Reyniers, G.B. Marin, V. Van Speybroeck, M. Waroquier
Journal of Physical Chemistry A
109 (33), 7466-7480
2005
A1

Abstract 

A complete and consistent set of 95 Benson group additive values (GAV) for the standard enthalpy of formation of hydrocarbons and hydrocarbon radicals at 298 K and 1 bar is derived from an extensive and accurate database of 233 ab initio standard enthalpies of formation, calculated at the CBS-QB3 level of theory. The accuracy of the database was further improved by adding newly determined bond additive corrections (BAC) to the CBS-QB3 enthalpies. The mean absolute deviation (MAD) for a training set of 51 hydrocarbons is better than 2 kJ mol-1. GAVs for 16 hydrocarbon groups, i.e., C(Cd)3(C), C−(Cd)4, C−(Ct)(Cd)(C)2, C−(Ct)(Cd)2(C), C−(Ct)(Cd)3, C−(Ct)2(C)2, C−(Ct)2(Cd)(C), C−(Ct)2(Cd)2, C−(Ct)3(C), C−(Ct)3(Cd), C−(Ct)4, C−(Cb)(Cd)(C)(H), C−(Cb)(Ct)(H)2, C−(Cb)(Ct)(C)(H), C−(Cb)(Ct)(C)2, Cd−(Cb)(Ct), for 25 hydrocarbon radical groups, and several ring strain corrections (RSC) are determined for the first time. The new parameters significantly extend the applicability of Benson's group additivity method. The extensive database allowed an evaluation of previously proposed methods to account for non-next-nearest neighbor interactions (NNI). Here, a novel consistent scheme is proposed to account for NNIs in radicals. In addition, hydrogen bond increments (HBI) are determined for the calculation of radical standard enthalpies of formation. In particular for resonance stabilized radicals, the HBI method provides an improvement over Benson's group additivity method.

Hydrogen Radical Additions to Unsaturated Hydrocarbons and the Reverse β-Scission Reactions: Modeling of Activation Energies and Pre-Exponential Factors

M. Sabbe, M-F. Reyniers, M. Waroquier, G.B. Marin
ChemPhysChem
11 (1), 195-210
2010
A1

Abstract 

The group additivity method for Arrhenius parameters is applied to hydrogen addition to alkenes and alkynes and the reverse β-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 β-scissions yields rate coefficients within a factor of 3.5 of the CBS-QB3 results except for two β-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.

Open Access version available at UGent repository

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

M. Sabbe, A.G. Vandeputte, M-F. Reyniers, M. Waroquier, G.B. Marin
Physical Chemistry Chemical Physics (PCCP)
12, 1278-1298
2010
A1

Abstract 

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 Ea 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.

Standard enthalpy of formation, entropy and heat capacity of hydrocarbons and hydrocarbon radicals: first principles group additive values

ISBN/ISSN:
Talk

Conference / event / venue 

Thermodynamics 2007
Rueil-Malmaison, France
Wednesday, 26 September, 2007 to Friday, 28 September, 2007
Subscribe to RSS - M. Sabbe