V. Van Speybroeck

Silicalite-1 zeolite exhibits a characteristic pentasil framework vibration around 540−550 cm−1. In the initial stages of zeolite synthesis, however, this band is observed at much higher wavenumbers: literature shows this vibration to depend on particle size and to shift over 100 cm−1 with increasing condensation. In this work, the pentasil vibration frequency was derived from theoretical molecular dynamics simulations to obtain the correct IR band assignments for important nanoparticles. The IR spectroscopic fingerprint of oligomeric five-ring containing precursors proposed in the literature was computed and compared with experimental data. Our theoretical results show that, while isolated five-membered rings show characteristic vibrational bands around 650 cm−1, the combination of five-membered rings in the full MFI-type structure readily generates the bathochromic shift to the typical pentasil vibration around 550 cm−1. As opposed to what was previously believed, the IR band does not shift gradually as nanoparticle size increases, but it is highly dependent on the specific way structural units are added. The most important feature is the appearance of an additional band when double five-membered rings are included, which allows for a clear distinction between the key stages of early zeolite nucleation. Furthermore, the combination of the simulated spectra with the experimental observation of this spectral feature in nanoparticles extracted from silicalite-1 clear solutions supports their structured nature. The theoretical insights on the dependency of pentasil vibrations with the degree of condensation offer valuable support toward future investigations on the genesis of a zeolite crystal.

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.

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.