M. Waroquier

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.

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.

Electron paramagnetic resonance (EPR), electron nuclear double resonance (ENDOR), and ENDOR-induced EPR (EIE) measurements on sucrose single crystals at 10 K after in situ X irradiation at this temperature reveal the presence of at least nine different radical species. Nine proton hyperfine coupling tensors were determined from ENDOR angular variations and assigned to six of these species (R1−R6) using EIE. Spectral simulations indicate that four of those (R1−R3 and R6) dominate the EPR absorption. Assisted by periodic density functional theory (DFT) calculations, R1 and R2 are identified as H-abstracted C1- and C5-centered radicals, R3 is tentatively assigned to an H-abstracted C6-centered radical, and R6 is identified as an alkoxy radical where the abstracted hydroxy proton has migrated to a neighboring OH group via intermolecular proton transfer. The latter radical had been characterized and identified in a previous study, but the present DFT calculations provide additional insight into its conformation and particular properties. This study provides the first step in unraveling the formation mechanism of the stable sucrose radicals detected after room-temperature irradiation and contributes to the understanding of the initial stages of radiation damage to solid-state carbohydrates.

Thermodynamic properties of phosphorus-containing compounds were investigated using high-level ab initio computations. An extended set of contemporary density functional theory (DFT) procedures was assessed for their ability to accurately predict bond dissociation energies of a set of phosphoranyl radicals. The results of meta- and double-hybrids as well as more recent methods, in particular M05, M05-2X, M06, and M06-2X, were compared with benchmark G3(MP2)-RAD values. Standard heats of formation, entropies, and heat capacities of a set of ten organophosphorus compounds were determined and the low-cost BMK functional was found to provide results consistent with available experimental data. In addition, bond dissociation enthalpies (BDEs) were computed using the BMK, M05-2X, and SCS-ROMP2 procedure. The three methods give the same stability trend. The BDEs of the phosphorus(III) molecules were found to be lower than their phosphorus(V) counterparts. Overall, the following ordering is found: BDE(P−OPh)