M-F. Reyniers

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.

An efficient procedure for normal-mode analysis of extended systems, such as zeolites, is developed and illustrated for the physisorption and chemisorption of n-octane and isobutene in H-ZSM-22 and H-FAU using periodic DFT calculations employing the Vienna Ab Initio Simulation Package. Physisorption and chemisorption entropies resulting from partial Hessian vibrational analysis (PHVA) differ at most 10 J mol−1 K−1 from those resulting from full Hessian vibrational analysis, even for PHVA schemes in which only a very limited number of atoms are considered free. To acquire a well-conditioned Hessian, much tighter optimization criteria than commonly used for electronic energy calculations in zeolites are required, i.e., at least an energy cutoff of 400 eV, maximum force of 0.02 eV/Å, and self-consistent field loop convergence criteria of 10−8 eV. For loosely bonded complexes the mobile adsorbate method is applied, in which frequency contributions originating from translational or rotational motions of the adsorbate are removed from the total partition function and replaced by free translational and/or rotational contributions. The frequencies corresponding with these translational and rotational modes can be selected unambiguously based on a mobile block Hessian−PHVA calculation, allowing the prediction of physisorption entropies within an accuracy of 10−15 J mol−1 K−1 as compared to experimental values. The approach presented in this study is useful for studies on other extended catalytic systems.