E. Pauwels

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 method for the calculation of hyperfine parameters in extended systems under periodic boundary conditions is presented, using the Gaussian and augmented-plane-wave density functional method, and implemented in QUICKSTEP. In order to increase the efficiency in larger systems, a hybrid scheme is proposed, in which an all-electron treatment for the nuclei of interest and a pseudopotential approximation for the remaining atoms in the simulation cell are combined. The method is validated first by comparing the hyperfine parameters for a selection of atoms and small molecules (using a supercell technique) with other theoretical methods and experimental data from literature. As a typical example of a periodic system where our hybrid method can be applied, the hyperfine parameters of the well-characterized R2 L-α-alanine derived radical are evaluated, yielding results in excellent agreement with the available experimental data.

Monoatomic X- (X = O, S) chalcogen centers in MZ (M = Na, K, Rb and Z = Cl, Br, I) alkali halide lattices are investigated within the framework of density functional theory with the principal aim to establish defect models. In electron paramagnetic resonance (EPR) experiments, X- defects with tetragonal, orthorhombic, and monoclinic g-tensor symmetry have been observed. In this paper, models in which X- replaces a single halide ion, with a next nearest neighbor and a nearest neighbor halide vacancy, are validated for the X- centers with tetragonal and orthorhombic symmetry, respectively. As such defect models are extended, the ability to reproduce experimental data is a stringent test for various computational approaches. Cluster in vacuo and embedded cluster schemes are used to calculate energy and EPR parameters for the two vacancy configurations. The final assignment of a defect structure is based on the qualitative and quantitative reproduction of experimental g and (super)hyperfine tensors.

Density functional theory techniques are used to investigate the defect structure of X- (X = O, S, Se) ions in MZ (M = Na, K, Rb and Z = Cl, Br) alkali halides which exhibit monoclinic-I g-tensor symmetry, using cluster in vacuo, embedded cluster, and periodic embedding schemes. Although a perturbed interstitial defect model was suggested from electron paramagnetic resonance experiments (EPR), the nature of the perturbation is still unknown. An appropriate defect model is developed theoretically by comparing structural and energetical properties of various defect configurations. Further validation is achieved by cross referencing experimental and computed EPR data. On the basis of the computational results, the following defect model is proposed:  the X- ion is located interstitially with a charge compensating halide vacancy in its first coordination shell.

With the aid of density functional theory calculations, all conformers of several single-chain alcohols, thiols, ethers, and sulfides are investigated. Starting from earlier computational works on n-alkanes, we construct an extended set of general rules for predicting the number and occurrence of conformers in these oxygen- or sulfur-containing compounds. In alcohols and thiols, it is found that only the conformers generated by internal rotations in the HXCH2CH2CH2 (X = O or S) top are distinctive from those in n-alkanes. In ethers and sulfides, the primary influence of the heteroelement also extends up to three internal rotations, but many more conformers are possible. However, a number of double gauche sequences are forbidden, and therefore, several conformers can be eliminated. These exclusions in particular make up a set of rules for eventually deducing all possible conformers. Furthermore, on the basis of only an exact calculation of these gg conformations in addition to single gauche conformers, it is possible to make an accurate estimate of the relative energy. This two-dimensional approximation scheme constitutes an effective tool for adequately describing the relative energies of all possible conformers at a minimal computational cost.