M. Waroquier

Primary free radical formations in fructose single crystals X-irradiated at 10 K were investigated at the same temperature using X-band Electron Paramagnetic Resonance (EPR), Electron Nuclear Double Resonance (ENDOR) and ENDOR induced EPR (EIE) techniques. ENDOR angular variations in the three principal crystallographic planes and a fourth skewed plane allowed the unambiguous determination of five proton hyperfine coupling tensors. From the EIE studies, these hyperfine interactions were assigned to three different radicals, labeled T1, T1* and T2. For the T1 and T1* radicals, the close similarity in hyperfine coupling tensors suggests that they are due to the same type of radical stabilized in two slightly different geometrical conformations. Periodic density functional theory calculations were used to aid the identification of the structure of the radiation-induced radicals. For the T1/T1* radicals a C3 centered hydroxyalkyl radical model formed by a net H abstraction is proposed. The T2 radical is proposed to be a C5 centered hydroxyalkyl radical, formed by a net hydrogen abstraction. For both radicals, a very good agreement between calculated and experimental hyperfine coupling tensors was obtained.

Pauwels, E., Declerck, R., Van Speybroeck, V. and Waroquier, M. Evidence for a Grotthuss-Like Mechanism in the Formation of the Rhamnose Alkoxy Radical Based on Periodic DFT Calculations. Radiat. Res. 169, 8-18 (2008). Molecular modeling adopting a periodic approach based on density functional theory (DFT) indicates that a Grotthuss-like mechanism is active in the formation of the radiation-induced alkoxy radical in alpha-l-rhamnose. Starting from an oxidized crystal structure, a hydroxyl proton is transferred along an infinite hydrogen bond chain pervading the entire crystal. The result of this proton shuttling mechanism is a stable radical species dubbed RHop. Only after several reorientations of crystal waters and hydroxyl groups, the more stable radical form RO4 is obtained, which differs in structure from the former by the absence of only one hydrogen bond. Calculations of the energetics associated with the mechanism as well as simulated spectroscopic properties reveal that different variants of the rhamnose alkoxy radical can be observed depending on the temperature of irradiation and consecutive EPR measurement. Cluster calculations on both radical variants provide hyperfine coupling and g tensors that are in good agreement with two independent experimental measurements at different temperatures.

A consistent set of group additive values Delta GAV 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 Delta GAV degrees s in the broad range of 300-1300K is limited to +/- 4.5 kJmol(-1) on activation energies and to +/- 0.4 on logA (A: pre-exponential factor) for 90% of the groups. Validation of the Delta GAV 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.

In an earlier paper, the authors have developed a new method, the mobile block Hessian (MBH), to accurately calculate vibrational modes for partially optimized molecular structures [J. Chem. Phys. 2007, 126 (22), 224102]. The proposed procedure remedies the artifact of imaginary frequencies, occurring in standard frequency calculations, when parts of the molecular system are optimized at different levels of theory. Frequencies are an essential ingredient in predicting reaction rate coefficients due to their input in the vibrational partition functions. The question arises whether the MBH method is able to describe the chemical reaction kinetics in an accurate way in large molecular systems where a full quantum chemical treatment at a reasonably high level of theory is unfeasible due to computational constraints. In this work, such a validation is tested in depth. The MBH method opens a lot of perspectives in predicting accurate kinetic parameters in chemical reactions where the standard full Hessian procedure fails.

Radical polymerization processes occur through a complex network of many different reactions. It is well-known that the polymerization rate is directly related to the monomer structure. The experimental polymerizability behavior is expressed as kp/kt1/2, where kp is the rate coefficient of propagation and kt is the rate coefficient of termination. In this study, the reactivity of a series of acrylates and methacrylates is modeled in order to understand the effect of the pendant group size, the polarity of a pendant group, and the nature of the pendant group (linear vs cyclic) on their polymerizability. The geometries and frequencies are calculated with the B3LYP/6-31+G(d) methodology whereas the energetics and kinetics of these monomers have been studied using the two-component BMK/6-311+G(3df,2p)//B3LYP/6-31+G(d) level of theory. For rotations about forming/breaking bonds in the transition state, an uncoupled scheme for internal rotations has been applied with potentials determined at the B3LYP/6-31+G(d) level. Generally the rate constants for propagation kp mimic the qualitative polymerization trend of the monomers modeled and can be used with confidence in predicting the polymerizability behavior of acrylates. However in the case of 2-dimethylaminoethyl acrylate, chain transfer is found to play a major role in inhibiting the polymerization. On the other hand, the disproportionation reaction turns out to be too slow to be taken into consideration as a termination reaction.