G.B. Marin

Ab Initio Study of Radical Reactions: Cyclization Pathways for the Butylbenzene Radical (II)

V. Van Speybroeck, Y. Borremans, D. Van Neck, M. Waroquier, S. Wauters, M. Saeys, G.B. Marin
Journal of Physical Chemistry A
105 (32), 7713–7723
2001
A1

Abstract 

Ab initio density functional theory calculations are presented on some model reactions involved in coke formation during the thermal cracking of hydrocarbons. The reactions under consideration are different cyclization pathways for the butylbenzene radical, which can lead to a further growth of the coke layer. This study enables us to gain more microscopic insight into the mechanistic and kinetic aspects of the reactions. Special attention is paid to the exact treatment of internal rotations and their impact on the kinetic parameters. Pre-exponential factors are very sensitive to the accuracy of constructing the microscopic partition functions. In particular, the relative importance of cyclization toward five and six-membered rings is studied on the basis of the calculated rate constants and concentration profiles of the reactants. The influence of the size of the ring and of the relative stability of the primary and secondary butylbenzene radical on the cyclization reaction is discussed. The activation energy for the formation of six-membered rings is approximately 30 kJ/mol lower than that for five-ring formation. The predicted values for the kinetic parameters enable us to validate some basic assumptions on coke formation. The calculations as presented here are especially important for complex reaction schemes, for which experimental data are not always available.

Ab initio study of radical addition reactions: Addition of a primary ethylbenzene radical to ethene (I)

V. Van Speybroeck, D. Van Neck, M. Waroquier, S. Wauters, M. Saeys, G.B. Marin
Journal of Physical Chemistry A
104 (46), 10939–10950
2000
A1

Abstract 

Ab initio density functional theory calculations have been carried out on a model reaction involved in coke formation during the thermal cracking of hydrocarbons, namely, the addition of the ethylbenzene radical to ethene. This study enables one to get more microscopic insight into the mechanistic and kinetic aspects of the reaction. A profound ab initio conformational analysis of the formed products, reactants, and transition states is made. The impact of internal rotations on the two kinetic parameters deduced from transition state theory (TST), the activation energy and the preexponential factor, has been studied in detail. Furthermore, we report on the various components that govern the kinetic parameters. Preexponential factors are very sensitive to the accuracy of constructing the microscopic partition functions. Internal rotations play a dominant role in the reaction mechanism, and their impact on the preexponential factor is large. Hence, a very accurate handling of internal rotations is of crucial importance. We present a new algorithm to extract exactly on a quantum mechanical basis the partition functions of the internal rotations. The calculations as presented here are especially important for complex reaction schemes, for which experimental data are not always available.

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.

The remarkable catalytic activity of the saturated metal organic framework V-MIL-47 in the cyclohexeneoxidation

K. Leus, I. Muylaert, M. Vandichel, G.B. Marin, M. Waroquier, V. Van Speybroeck, P. Van der Voort
Chemical Communications
46, 5085-5087
2010
A1

Abstract 

The remarkable catalytic activity of the saturated metal organic framework MIL-47 in the epoxidation of cyclohexene is elucidated by means of both experimental results and theoretical calculations.

Open Access version available at UGent repository

First principle kinetic studies of zeolite-catalyzed methylation reactions

V. Van Speybroeck, J. Van der Mynsbrugge, M. Vandichel, K. Hemelsoet, D. Lesthaeghe, A. Ghysels, G.B. Marin, M. Waroquier
JACS (Journal of the American Chemical Society)
133 (4), 888–899
2011
A1

Abstract 

Methylations of ethene, propene, and butene by methanol over the acidic microporous H-ZSM-5 catalyst are studied by means of state of the art computational techniques, to derive Arrhenius plots and rate constants from first principles that can directly be compared with the experimental data. For these key elementary reactions in the methanol to hydrocarbons (MTH) process, direct kinetic data became available only recently [J. Catal.2005, 224, 115−123; J. Catal.2005, 234, 385−400]. At 350 °C, apparent activation energies of 103, 69, and 45 kJ/mol and rate constants of 2.6 × 10−4, 4.5 × 10−3, and 1.3 × 10−2 mol/(g h mbar) for ethene, propene, and butene were derived, giving following relative ratios for methylation kethene/kpropene/kbutene = 1:17:50. In this work, rate constants including pre-exponential factors are calculated which give very good agreement with the experimental data: apparent activation energies of 94, 62, and 37 kJ/mol for ethene, propene, and butene are found, and relative ratios of methylation kethene/kpropene/kbutene = 1:23:763. The entropies of gas phase alkenes are underestimated in the harmonic oscillator approximation due to the occurrence of internal rotations. These low vibrational modes were substituted by manually constructed partition functions. Overall, the absolute reaction rates can be calculated with near chemical accuracy, and qualitative trends are very well reproduced. In addition, the proposed scheme is computationally very efficient and constitutes significant progress in kinetic modeling of reactions in heterogeneous catalysis.

Normal Mode Analysis in Zeolites: Toward an Efficient Calculation of Adsorption Entropies

B. De Moor, A. Ghysels, M-F. Reyniers, V. Van Speybroeck, M. Waroquier, G.B. Marin
Journal of Chemical Theory and Computation (JCTC)
7(4), 1090-1101
2011
A1

Abstract 

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.

Catalytic and molecular separation properties of Zeogrids and Zeotiles

J.A. Martens, J.W. Thybaut, J.F.M. Denayer, S. Pulinthanathu Sree, A. Aerts, M-F. Reyniers, V. Van Speybroeck, M. Waroquier, A. Buekenhoudt, I. Vankelecom, W. Buijs, J. Persoons, G.V. Baron, S. Bals, G. Van Tendeloo, G.B. Marin, P.A. Jacobs, C. Kirschhock
Catalysis Today
168, 17-27
2011
A1

Abstract 

Zeogrids and Zeotiles are hierarchical materials built from assembled MFI zeolite precursor units. Permanent secondary porosity in these materials is obtained through self assembly of nanoparticles encountered in MFI zeolite synthesis in the presence of supramolecular templates. Hereon, the aggregated species are termed nanoslabs. Zeogrids are layered materials with lateral spacings between nanoslabs creating galleries qualifying as supermicropores. Zeotiles present a diversity of tridimensional nanoslab assemblies with mesopores. Zeotile-1, -4 and -6 are hexagonal mesostructures. Zeotile-1 has triangular and hexagonal channels; Zeotile-4 has hexagonal channels interconnected via slits. Zeotile-2 has a cubic structure with gyroid type mesoporosity. The behavior of Zeogrids and Zeotiles in adsorption, membrane and chromatographic separation and catalysis has been characterized and compared with zeolites and mesoporous materials derived from unstructured silica sources. Shape selectivity was detected via adsorption of n- and iso-alkanes. The mesoporosity of Zeotiles can be exploited in chromatographic separation of biomolecules. Zeotiles present attractive separation properties relevant to CO2 sequestration. Because of its facile synthesis procedure without hydrothermal steps Zeogrid is convenient for membrane synthesis. The performance of Zeogrid membrane in gas separation, nanofiltration and pervaporation is reported. In the Beckmann rearrangement of cyclohexanone oxime Zeogrids and Zeotiles display a catalytic activity characteristic of silicalite-1 zeolites. Introduction of acidity and redox catalytic activity can be achieved via incorporation of Al and Ti atoms in the nanoslabs during synthesis. Zeogrids are active in hydrocracking, catalytic cracking, alkylation and epoxidation reactions. Zeogrids and Zeotiles often behave differently from ordered mesoporous materials as well as from zeolites and present a valuable extension of the family of hierarchical silicate based materials.

The coordinatively saturated vanadium MIL-47 as a low leaching heterogeneous catalyst in the oxidation of cyclohexene

K. Leus, M. Vandichel, Y-Y Liu, I. Muylaert, J. Musschoot, H. Vrielinck, F. Callens, G.B. Marin, C. Detavernier, Y.Z. Khimyak, M. Waroquier, V. Van Speybroeck, P. Van der Voort
Journal of Catalysis
285 (1) 196-207
2012
A1

Abstract 

A Metal Organic Framework, containing coordinatively saturated V+IV sites linked together by terephthalic linkers (V-MIL-47), is evaluated as a catalyst in the epoxidation of cyclohexene. Different solvents and conditions are tested and compared. If the oxidant TBHP is dissolved in water, a significant leaching of V-species into the solution is observed, and also radical pathways are prominently operative leading to the formation of an adduct between the peroxide and cyclohexene. If, however, the oxidant is dissolved in decane, leaching is negligible and the structural integrity of the V-MIL-47 is maintained during successive runs. The selectivity toward the epoxide is very high in these circumstances. Extensive computational modeling is performed to show that several reaction cycles are possible. EPR and NMR measurements confirm that at least two parallel catalytic cycles are co-existing: one with V+IV sites and one with pre-oxidized V+V sites, and this is in complete agreement with the theoretical predictions.

Open Access version available at UGent repository

Pages

Subscribe to RSS - G.B. Marin