Combined electron magnetic resonance and density functional theory study of 10 K X-irradiated beta-D-fructose single crystals

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.     

Q-band EPR and ENDOR of low temperature X-irradiated beta-D-fructose single crystals

Beta-D-fructose single crystals were in situ X-irradiated at 80 K and measured using electron paramagnetic resonance (EPR), electron nuclear double resonance (ENDOR) and ENDOR-induced EPR (EIE) techniques at Q-band (34 GHz) microwave frequencies. The measurements revealed the presence of at least four carbon-centered radicals stable at 80 K. By means of ENDOR angular variations in the three principal crystallographic planes, six proton hyperfine coupling tensors could be determined and were assigned to four different radicals by the aid of EIE. Two of the radicals exhibit only beta-proton hyperfine couplings and reveal almost identical EIE spectra. For the other two radicals, the major hyperfine splitting originates from a single alpha-proton hyperfine coupling and their EIE spectra were also quite similar. The similarity of the EIE spectra and hyperfine tensors led to the assumption that there are only two essentially different radical structures. The radical exhibiting only beta-proton hyperfine couplings was assigned to a C3 centered radical arising from H3 abstraction and the other radical suggested to be an open-ring species with a disrupted C2-C3 bond and a double C2-O2 bond. A possible formation mechanism for the latter open-ring radical is presented. By means of cluster density functional theory (DFT) calculations, the structures of the two radicals were determined and a fairly good agreement between the calculated and experimental hyperfine tensors was found.     

Radiation-induced radicals in alpha-D-glucose: Comparing DFT cluster calculations with magnetic resonance experiments

Using density functional theory (DFT) calculations, an enhanced theoretical examination was made of the radiation-induced radicals in alpha-d-glucose. For the carbon-centred radicals in this sugar, the effect of the model space on the radical geometry as well as on the calculated radical hyperfine coupling tensors was examined. The findings were compared with previously published tensors, as determined by electron paramagnetic resonance (EPR) experiments and single molecule DFT calculations. A cluster approach was adopted, in which intermolecular interactions (predominantly hydrogen bonds) between the radical species and its environment were explicitly incorporated. This substantially improved the correspondence with experimental findings in comparison with single molecule calculations of an earlier examination. In a direct comparison between both computational methods for the glucose radicals, it was shown that the extent of the model space plays an important part in the determination of the radical geometry. Furthermore, the model space also has an impact on the calculated hyperfine coupling tensors. Full cluster EPR calculations, in which the paramagnetic properties are calculated for the entire model space of the cluster, give an excellent agreement with the experimental EPR measurements.     

Study of rhamnose radicals in the solid state adopting a density functional theory cluster approach

A theoretical study is performed on the radiation-induced radicals in crystalline alpha-l-rhamnose, using density functional theory (DFT) calculations. Irrespective of earlier structural assignments, a host of possible radical models is examined in search for a structure that accurately reproduces experimental electron paramagnetic resonance (EPR) properties. A cluster approach is followed, incorporating all hydrogen bond interactions between radical and crystalline environment. Hyperfine coupling tensors as well as g tensors are determined and a comparison is made with available experimental data. Three carbon-centered hydroxyalkyl radicals are validated, in accordance with experimental suggestions for their structure. The occurrence of a carbon-centered oxygen anion radical for one of the radical species is rejected on theoretical grounds, and instead an altered hydroxyalkyl structure is suggested. Our cluster calculations are able to determine g and hyperfine tensors for the oxygen-centered alkoxy radical in rhamnose, in accordance with one of the two measurements for this species. For all radical models, quantitative agreement with experimental hyperfine tensors is obtained by performing full cluster DFT calculations. The inclusion of the molecular environment for the determination of this EPR property proved to be essential.     

Comparative density-functional study of the electron paramagnetic resonance parameters of amavadin

Electronic g tensors and hyperfine coupling tensors have been calculated for amavadin, an unusual eight-coordinate vanadium(IV) complex isolated from Amanita muscaria mushrooms. Different density-functional methods have been compared, ranging from local via gradient-corrected to hybrid functionals with a variable Hartree-Fock exchange admixture. For both electron paramagnetic resonance (EPR) properties, hybrid functionals with an appreciable exact-exchange admixture provide the closest agreement with experimental data. Second-order spin-orbit corrections provide non-negligible contributions to the 51V hyperfine tensor. The orientation of g and A tensors relative to each other also depends on spin-orbit corrections to the A tensor. A rationalization for the close resemblance of the EPR parameters of amavadin to those of the structurally rather different vanadyl complexes is provided, based on the nature of the relevant frontier orbitals.     

ESR study of the carbon 13 hyperfine coupling tensors of the (COOH)H(OH) and COOH(OH2) radicals

(COOH)(H)OH and (COOH)(OH)₂ radicals trapped in single crystals of glycolic and mesoxalic acids are studied by EPR. ¹³C and hydrogen coupling hyperfine tensors are determined. It is concluded from results relative to ¹³C couplings that both radicals are planar.     

Automatic fitting to 'powder' EPR spectra of coupled paramagnetic species employing Feynman's theorem

A previous automatic fitting procedure of EPR spectra has been extended with the purpose to characterise coupled paramagnetic complexes in powders and frozen solutions. The theoretical EPR spectra were obtained by matrix diagonalization of a general spin Hamiltonian. A least-squares fitting procedure using analytical derivatives of the calculated spectrum with respect to the spectroscopic, fine structure, nuclear quadrupole, electron-electron, and hyperfine coupling tensors was used to refine those parameters. The powder spectra of matrix isolated *CF3 and RCF2CF2* radicals, previously measured at low temperature, were reanalysed with this method. A theoretically modeled complex consisting of a Cu2+ ion, featuring an axially symmetric g-tensor and 63Cu hyperfine structure anisotropy, and a free radical located at different orientations, with respect to the symmetry axis of the Cu2+ ion, was examined in order to investigate the possibility to recover the magnetic parameters of the separate units and the magnetic couplings between them.     

HYSCORE and Davies ENDOR study of irradiated ultra high molecular weight polyethylene

Ultra high molecular weight polyethylene (UHMWPE) has been studied with different magnetic resonance techniques to elicit information on the nature and the location of radicals generated during high energy irradiation. Field swept electron paramagnetic resonance, pulsed Davies electron nuclear double resonance and hyperfine sublevel correlation spectroscopic measurements allowed extracting for the first time the full ¹H hyperfine coupling tensors of the most abundant radical, i.e. a secondary alkyl radical and to ascertain the formation of allyl radicals in the first stages of the irradiation process. The ¹H hyperfine coupling tensors are analogous to those reported for single crystal irradiated polyethylene, suggesting that radicals generated in UHMWPE are located in the crystalline region of the polymer. Copyright © 2012 John Wiley & Sons, Ltd.     

Determining the geometry and magnetic parameters of fluorinated radicals by simulation of powder ESR spectra and DFT calculations: the case of the radical RCF2CF2* in nafion perfluorinated ionomers

The ESR spectrum of the chain-end radical RCF2CF2* detected in Nafion perfluorinated membranes exposed to the photo-Fenton reagent was accurately simulated by an automatic fitting procedure, using as input the hyperfine coupling tensors of the two F alpha and two F beta nuclei as well as the corresponding directions of the principal values from density functional theory (DFT) calculations. An accurate fit was obtained only for different orientations of the hyperfine coupling tensors for the two F alpha nuclei, indicating a nonplanar structure about the C alpha radical center. The fitted isotropic hyperfine splittings for the two F beta nuclei in the Nafion radical, 24.9 and 27.5 G, are significantly larger than those for the chain-end radical in Teflon (15 G), implying different radical conformations in the two systems. The excellent fit indicated that the geometry and electronic structure of free radicals can be obtained not only from single-crystal ESR spectroscopy, but also, in certain cases, from powder spectra, by combination with data from DFT calculations. The optimized structures obtained by DFT calculations for the CF3CF2CF2CF2* or CF3OCF2CF2* radicals as models provided additional support for the pyramidal structure determined from the spectral fit. Comparison and analysis of calculated and fitted values for the hyperfine splittings of the two F beta nuclei suggested that the radical detected by ESR in Nafion is ROCF2CF2*, which originates from attack of oxygen radicals on the Nafion side chain. The combination of spectrum fitting and DFT is considered important in terms of understanding the hyperfine splittings from 19F nuclei and the different conformations of fluorinated chain-end-type radicals RCF2CF2* in different systems, and also for elucidating the mechanism of Nafion fragmentation when exposed to oxygen radicals in fuel cell conditions.     

Radiation products at 77 K in trehalose single crystals: EMR and DFT analysis

The radicals obtained in trehalose dihydrate single crystals after 77 K X-irradiation have been investigated at the same temperature using X-band electron paramagnetic resonance (EPR), electron nuclear double resonance (ENDOR), and ENDOR-induced EPR (EIE) techniques. Five proton hyperfine coupling tensors were unambiguously determined from the ENDOR measurements and assigned to three carbon-centered radical species (T1, T1*, and T2) based on the EIE spectra. EPR angular variations revealed the presence of four additional alkoxy radical species (T3 to T6) and allowed determination of their g tensors. Using periodic density functional theory (DFT) calculations, T1/T1*, T2, and T3 were identified as H-loss species centered at C4, C1', and O2', respectively. The T4 radical is proposed to have the unpaired electron at O4, but considerable discrepancies between experimental and calculated HFC values indicate it is not simply the (net) H-loss species. No suitable models were found for T5 and T6. These exhibit a markedly larger g anisotropy than T3 and T4, which were not reproduced by any of our DFT calculations.     

Electronic structure and hyperfine interactions for hydrazyl radicals

Some aspects of electronic structure and hyperfine interactions for α,α-diphenyl-β-picrylhydrazyl and its derivatives have been studied. High resolution EPR spectra have been obtained and analyzed by the least-quares method in Fourier space. The results obtained are confirmed by NMR experiments. Anisotropic EPR spectra of hydrazyl radicals in frozen solutions and glasses have been recorded. Parameters of hyperfine tensors and g-tensor have been determined by the least-squares method. Spin distributions in radical molecules are calculated by the unrestricted Hartree—Fock method and calculations of isotropic and anisotropic hyperfine interactions have been carried out.     

Computational studies of EPR parameters for paramagnetic molybdenum complexes. II. Larger MoV systems relevant to molybdenum enzymes

The careful validation of modern density functional methods for the computation of electron paramagnetic resonance (EPR) parameters in molybdenum complexes has been extended to a number of low-symmetry MoV systems that model molybdoenzyme active sites. Both g and hyperfine tensors tend to be reproduced best by hybrid density functionals with about 30-40% exact-exchange admixture, with no particular spin contamination problems encountered. Spin-orbit corrections to hyperfine tensors are mandatory for quantitative and, in some cases, even for qualitative agreement. The g11 (g||) component of the g tensor tends to come out too positive when spin-orbit coupling is included only to leading order in perturbation theory. Compared to single-crystal experiments, the calculations reproduce both g- and hyperfine-tensor orientations well, both relative to each other and to the molecular framework. This is significant, as simulations of the EPR spectra of natural-abundance frozen-solution samples frequently do not allow a reliable determination of the hyperfine tensors. These may now be extracted based on the quantum-chemically calculated parameters. In a number of cases, revised simulations of the experimental spectra have brought theory and experiment into substantially improved agreement. Systems with two terminal oxo ligands, and to some extent with an oxo and a sulfido ligand, have been confirmed to exhibit particularly large negative Deltag33 shifts and thus large g anisotropies. This is discussed in the context of the experimental data for xanthine oxidase.     

Identification of primary free radicals in trehalose dihydrate single crystals X-irradiated at 10 K

Primary free radical formation in trehalose dihydrate single crystals X-irradiated at 10 K was investigated at the same temperature using X-band Electron Paramagnetic Resonance (EPR), Electron Nuclear Double Resonance (ENDOR) and ENDOR-induced EPR (EIE) techniques. The ENDOR results allowed the unambiguous determination of six proton hyperfine coupling (HFC) tensors. Using the EIE technique, these HF interactions were assigned to three different radicals, labeled R1, R2 and R3. The anisotropy of the EPR and EIE spectra indicated that R1 and R2 are alkyl radicals (i.e. carbon-centered) and R3 is an alkoxy radical (i.e. oxygen-centered). The EPR data also revealed the presence of an additional alkoxy radical species, labeled R4. Molecular modeling using periodic Density Functional Theory (DFT) calculations for simulating experimental data suggests that R1 and R2 are the hydrogen-abstracted alkyl species centered at C5' and C5, respectively, while the alkoxy radicals R3 and R4 have the unpaired electron localized mainly at O2 and O4'. Interestingly, the DFT study on R4 demonstrates that the trapping of a transferred proton can significantly influence the conformation of a deprotonated cation. Comparison of these results with those obtained from sucrose single crystals X-irradiated at 10 K indicates that the carbon situated next to the ring oxygen and connected to the CH(2)OH hydroxymethyl group is a better radical trapping site than other positions.     

Computational studies of electron paramagnetic resonance parameters for paramagnetic molybdenum complexes. 1. Method validation on small and medium-sized systems

A variety of density functional methods have been evaluated in the computation of electronic g-tensors and molybdenum hyperfine couplings for systems ranging from the Mo atom through MoIIIN, [MoVOCl4]-, and [MoVOF5]2- to two larger MoV complexes MoXLCl2 (X=O, S; L=tris(3,5-dimethylpyrazolyl)hydroborate anion). In particular, the influence of the molybdenum basis set and of various exchange-correlation functionals with variable admixtures of Hartree-Fock exchange on the computed EPR parameters have been evaluated in detail. Careful basis-set studies have provided a moderate-sized 12s6p5d all-electron basis on molybdenum that gives hyperfine tensors in excellent agreement with much larger basis sets and that will be useful for calculations on larger systems. The best agreement with experimental data for both hyperfine and g-tensors is obtained with hybrid functionals containing approximately 30-40% Hartree-Fock exchange. Only for MoSLCl2 does increasing spin contamination with increasing exact-exchange admixture restrict the achievable computational accuracy. In all cases, spin-orbit corrections to the hyperfine tensors are sizable and have to be included in accurate calculations. Scalar relativistic effects enhance the isotropic Mo hyperfine coupling by approximately 15-20%. Two-component g-tensor calculations with variational inclusion of spin-orbit coupling show that the Deltag parallel components in [MoVOCl4]- and [MoVOF5]2- depend on higher-order spin-orbit contributions and are thus described insufficiently by the usual second-order perturbation approaches. Computed orientations of g- and hyperfine tensors relative to each other and to the molecular framework for the MoXLCl2 complexes provide good agreement between theory and single-crystal electron paramagnetic resonance experiments. In these cases, the hyperfine tensor orientations are influenced only slightly by spin-orbit effects.     

Role of the geometry, restricted rotations and solvents on the computed 2,2'-diphenyl-1-picrylhydrazyl hyperfine tensors

The nuclear hyperfine tensor (A) components of the 2,2'-diphenyl-1-picrylhydrazyl neutral radical are computed using the UB1LYP hybrid density functional method. Solvent interactions via hydrogen bonding are found to play a crucial role in the position of the two phenyl rings relative to the picryl moiety. Under these conditions, the calculated isotropic hyperfine tensor components of the N 1 and N 2 hydrazyl backbone are within approximately 1.3 Gauss (G) of the experimental values determined by EPR and ENDOR spectroscopy. Just as important are the effects of restricted rotations of the phenyl rings on these tensors. Rotational averaging using a Maxwell-Boltzmann type distribution improves the agreement between theory and experiment to less than 1.0 G. In addition, rotational averaging of the twelve isotropic proton coupling constants has also been performed. They come within 0.3 G of the experimental values. Thus, for the first time, all the nuclear hyperfine tensor components of this large class of molecules are accurately calculated without resorting to post Hartree-Fock techniques.     

Covalent bonding in copper(II) complexes with hetarylformazans according to EPR spectra

Covalent bonding in a number of copper(II) complexes with hetarylformazans that have pseudotetrahedral or square-planar symmetry of the nearest metal environment was analyzed from EPR spectra. The dependence of the unpaired electron delocalization on the pseudotetrahedral distortion of the coordination polyhedron was determined. A change in the Zeeman coupling parameters was interpreted. Various contributions to the components of the hyperfine coupling (HFC) and ligand hyperfine coupling (LHFC) tensors were calculated. pd-Mixing of the AO of the copper ion was found to have a slight effect on the HFC parameters. In the components of the LHFC tensor, the contribution from isotronic LHFC is decisive.     

Magnetic properties of single crystal alpha-benzoin oxime: An EPR study

The electron paramagnetic resonance (EPR) spectra of gamma irradiated single crystals of alpha-benzoinoxime (ABO) have been examined between 120 and 440 K. Considering the dependence on temperature and the orientation of the spectra of single crystals in the magnetic field, we identified two different radicals formed in irradiated ABO single crystals. To theoretically determine the types of radicals, the most stable structure of ABO was obtained by molecular mechanic and B3LYP/6-31G(d,p) calculations. Four possible radicals were modeled and EPR parameters were calculated for the modeled radicals using the B3LYP method and the TZVP basis set. Calculated values of two modeled radicals were in strong agreement with experimental EPR parameters determined from the spectra. Additional simulated spectra of the modeled radicals, where calculated hyperfine coupling constants were used as starting points for simulations, were well matched with experimental spectra.     

Structure and electron paramagnetic resonance parameters of the manganese site of concanavalin A studied by density functional methods

The EPR parameters of the manganese site in the saccharide-binding protein concanavalin A have been studied by density functional methods, with an emphasis on metal (⁵⁵Mn) and ligand (¹H and ¹⁷O) hyperfine couplings, in comparison with high-field EPR and ENDOR data. Results for gradient-corrected and hybrid functionals with different exact-exchange admixture have been compared with experiment for the ⁵⁵Mn and the ¹H ligand hyperfine coupling and have been predicted for ¹⁷O hyperfine coupling based on comparison with experiment for the related [Mn(H₂O)₆]²⁺. Appreciable exact-exchange admixture in the hybrid functional is needed to obtain an adequate spin-density distribution and thus near-quantitative agreement with experimental EPR parameters. The common use of experimental proton hyperfine coupling tensors together with the point-dipole approximation for determination of bond lengths is evaluated by explicit calculations.     

Simplification of complex EPR spectra by cepstral analysis

As the Fourier transform of time-series data is known as the spectrum, the Fourier transform of the logarithm of the time-series data is called the cepstrum of the data. When cepstral analysis is applied to free induction decay signals of free radicals showing first-order EPR spectra, the identification of nuclear hyperfine coupling constants becomes simple. In a systematic manner, we have examined how the technique of cepstral analysis is affected by the presence of aliasing, noise, uncertainty in the time origin of the free induction decay, the presence of second-order hyperfine couplings, and the applications of various apodization methods. This technique was then applied to analyze the EPR spectrum of anthraquinone anion radical, and anion radicals of porphycene and tetrapropyl-porphycene, and the hyperfine coupling constants thus obtained were compared with published data. A good agreement was always found. We make a case for the usefulness of cepstral analysis in determining the hyperfine coupling constants of complex EPR spectra of organic free radicals.     

EPR and DFT study on the stabilization of radiation-generated methyl radicals in dehydrated Na-A zeolite

Electron paramagnetic resonance (EPR) spectroscopy was applied to study paramagnetic species stabilized in Na-A zeolite exposed to gaseous methane and gamma-irradiated at 77 K. Two types of EPR spectra were recorded during thermal annealing of zeolite up to room temperature. Owing to the results for the zeolite exposed to (13)CH(4) the multiplet observed at 110 K was assigned to a (.-)CH(3)...Na(+) complex. After decay of the multiplet, the isotropic quartet of methyl radical was recorded in the temperature range of 170-280 K. On the basis of the EPR parameters it is postulated that (.-)CH(3) radicals in this temperature region are able to freely rotate inside the zeolite cage. The structures of the (.-)CH(3)...Na(+) adsorption complex and respective hyperfine coupling constants were calculated by applying DFT quantum chemical methods. Two different models were applied to represent the zeolite framework: the 6T structure of one six-membered ring and the 3T cluster. The hyperfine coupling constants calculated for the (.-)CH(3)...Na(+) adsorption complex for both applied models show very good agreement with those obtained experimentally.