Ab initio modeling of CW-ESR spectra of the double spin labeled peptide Fmoc-(Aib-Aib-TOAC)2-Aib-OMe in acetonitrile

In this work we address the interpretation, via an ab initio integrated computational approach, of the CW-ESR spectra of the double spin labeled, 310-helical, peptide Fmoc-(Aib-Aib-TOAC)2-Aib-OMe dissolved in acetonitrile. Our approach is based on the determination of geometric and local magnetic parameters of the heptapeptide by quantum mechanical density functional calculations taking into account solvent and, when needed, vibrational averaging contributions. The system is then described by a stochastic Liouville equation for the two electron spins interacting with each other and with two 14N nuclear spins, in the presence of diffusive rotational dynamics. Parametrization of the diffusion rotational tensor is provided by a hydrodynamic model. CW-ESR spectra are simulated with minimal resorting to fitting procedures, proving that the combination of sensitive ESR spectroscopy and sophisticated modeling can be highly helpful in providing 3D structural and dynamic information on molecular systems.     

On the interpretation of continuous wave electron spin resonance spectra of tempo-palmitate in 5-cyanobiphenyl

Electron spin resonance (ESR) measurements are highly informative on the dynamic behavior of molecules, which is of fundamental importance to understand their stability, biological functions and activities, and catalytic action. The wealth of dynamic information which can be extracted from a continuous wave electron spin resonance (cw-ESR) spectrum can be inferred by a basic theoretical approach defined within the stochastic Liouville equation formalism, i.e., the direct inclusion of motional dynamics in the form of stochastic (Fokker-Planck/diffusive) operators in the super Hamiltonian H governing the time evolution of the system. Modeling requires the characterization of magnetic parameters (e.g., hyperfine and Zeeman tensors) and the calculation of ESR observables in terms of spectral densities. The magnetic observables can be pursued by the employment of density functional theory which is apt, provided that hybrid functionals are employed, for the accurate computation of structural properties of molecular systems. Recently, an ab initio integrated computational approach to the in silico interpretation of cw-ESR spectra of multilabeled systems in isotropic fluids has been discussed. In this work we present the extension to the case of nematic liquid crystalline environments by performing simulations of the ESR spectra of the prototypical nitroxide probe 4-(hexadecanoyloxy)-2,2,6,6-tetramethylpiperidine-1-oxy in isotropic and nematic phases of 5-cyanobiphenyl. We first discuss the basic ingredients of the integrated approach, i.e., (1) determination of geometric and local magnetic parameters by quantum-mechanical calculations, taking into account the solvent and, when needed, the vibrational averaging contributions; (2) numerical solution of a stochastic Liouville equation in the presence of diffusive rotational dynamics, based on (3) parameterization of diffusion rotational tensor provided by a hydrodynamic model. Next we present simulated spectra with minimal resorting to fitting procedures, proving that the combination of sensitive ESR spectroscopy and sophisticated modeling can be highly helpful in providing three-dimensional structural and dynamic information on molecular systems in anisotropic environments.     

Ab initio study of the ground state properties of molecular oxygen

The electric and magnetic properties of the ground state of oxygen molecule are calculated by multiconfiguration self-consisted field (MCSCF) method and compared with experimental data: the quadrupole moment, polarizability, the 17O nuclear quadrupole coupling constant, magnetizability tensor, nuclear spin-rotation coupling constant and rotational g factor. The last two constants are calculated for all possible isotope modifications. The rotational, ESR and NMR spectra are discussed. Fermi-contact hyperfine coupling parameter is additionally estimated by different methods. The NMR chemical shielding tensor for 17O16O species at high temperature limit (without electron spin contribution) is predicted. Potential energy curves for 10 excited bound states and the internuclear distance dependence of the studied properties are also presented.     

Calculating slow-motional electron paramagnetic resonance spectra from molecular dynamics using a diffusion operator approach

A number of groups have utilized molecular dynamics (MD) to calculate slow-motional electron paramagnetic resonance (EPR) spectra of spin labels attached to biomolecules. Nearly all such calculations have been based on some variant of the trajectory method introduced by Robinson, Slutsky and Auteri (J. Chem. Phys. 1992,96, 2609-2616). Here we present an alternative approach that is specifically adapted to the diffusion operator-based stochastic Liouville equation (SLE) formalism that is also widely used to calculate slow-motional EPR line shapes. Specifically, the method utilizes MD trajectories to derive diffusion parameters such as the rotational diffusion tensor, diffusion tilt angles, and expansion coefficients of the orienting potential, which are then used as direct inputs to the SLE line shape program. This approach leads to a considerable improvement in computational efficiency over trajectory-based methods, particularly for high frequency, high field EPR. It also provides a basis for deconvoluting the effects of local spin label motion and overall motion of the labeled molecule or domain: once the local motion has been characterized by this approach, the label diffusion parameters may be used in conjunction with line shape analysis at lower EPR frequencies to characterize global motions. The method is validated by comparison of the MD predicted line shapes to experimental high frequency (250 GHz) EPR spectra.     

On the growth and shape of sodium taurodeoxycholate micellar aggregates: a spin-label and quasielastic light scattering investigation

Electron spin resonance (ESR) and quasielastic laser scattering (QELS) measurements have been carried out on sodium taurodeoxycholate (NaTDC) micellar aqueous solutions. Computer simulation of the ESR line shape has been used to quantitatively analyze the rotational dynamics of the cholestan-spin label (CSL) dissolved by the NaTDC micellar aggregates as a function of temperature and NaCl concentration. The local reorientation of CSL has been accounted for motionally-averaged g- and A-tensors assuming fast oscillation around the spin-probe long molecular axis. The overall Brownian tumbling of CSL-micelle complexes has been modeled by an axially symmetric rotational tensor. Good agreement with experimental spectra is obtained. Best-fit rotational parameters and QELS data suggest that, in the circumstance of large aggregation, NaTDC micelles have cylindrical shape and micellar growth occurs along the cylinder axis.     

Microphase separation in poly(acrylonitrile–butadiene–styrene) (ABS) studied with paramagnetic spin probes. II. Simulation of electron spin resonance spectra

We recently presented electron spin resonance spectra of poly(acrylonitrile–butadiene–styrene) (ABS) doped with 10‐doxylnonadecane (10DND) and 5‐doxyldecane (5DD) as spin probes. The spectra were measured in three types of ABS that differed in their butadiene contents and methods of preparation. Results for the ABS polymers were evaluated by comparison with similar studies on the homopolymers polybutadiene (PB) and polystyrene (PS) and the copolymers poly(styrene‐co‐acrylonitrile) (SAN) and poly(styrene‐co‐butadiene) (SB). Only one spectral component was detected for 10DND in PB, PS, SAN, and SB. In contrast, two spectral components differing in their dynamic properties were detected in the ABS samples and were assigned to spin probes located in butadiene‐rich domains (the fast component) and SAN‐rich domains (the slow component). The presence of two spectral components was taken as an indication of microphase separation. In this study, we present details on the dynamics and microphase separation by simulating spectra of 10DND in ABS, PB, PS, and SAN. The simulations are based on a dynamic model defined by the components of the rotational diffusion tensor and the diffusion tilt angle between the symmetry axis of the rotational diffusion tensor and the direction of the nitrogen 2p_z atomic orbital. The jump diffusion model led to good agreement with experimental spectra. In this model, the spin probe has a fixed orientation for a given time and then jumps instantaneously to a new orientation. The temperature variation of the rotational correlation time in PB and PS consisted of two dynamic regimes, with different activation energies. The transition temperature at which the change in dynamics occurs (T_tr) is 380 K for PS and 205 K for PB, essentially the same as the corresponding glass‐transition temperatures measured by differential scanning calorimetry. We suggest that T_tr is a better indicator of the glass transition than the temperature at which the total spectral width is 50 G, especially for large probes. The simulation program allowed the determination of the relative intensities of the fast and slow spectral components as a function of temperature; this information was used to clarify the redistribution of the probe above the glass transition of the SAN‐rich component in ABS systems. © 2002 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 40: 424–433, 2002; DOI 10.1002/polb.10110     

Stochastic modeling of CW-ESR spectroscopy of [60]fulleropyrrolidine bisadducts with nitroxide probes

In this work, we address the interpretation of continuous wave electron spin resonance (CW-ESR) spectra of fulleropyrrolidine bisadducts with nitroxide addends. Our approach is based on a definition of the spin Hamiltonian which includes exchange and dipolar interactions and on a complete numerical solution of the resulting stochastic Liouville equation, with inclusion of diffusive rotational dynamics. CW-ESR spectra are simulated for a series of C60 bisadducts made up of four trans isomers and the equatorial isomer. A nonlinear least-squares fitting procedure allows extraction directly from the available experimental spectra of a wide range of parameters, namely interprobe relative distances, diffusion tensors, and values of the exchange parameter J. Results are in good agreement with previous, more phenomenological estimates, proving that the combination of sensitive ESR spectroscopy based on multiple spin labeling with nitroxide radicals and sophisticated modeling can be highly helpful in providing structural and dynamic information on molecular systems.     

Molecular motion of a nickel-bis(dithiolato) complex in solution

The molecular motion of the planar bis(maleonitriledithiolato)nickel anion, Ni(mnt)(2)(-), has been studied as a function of temperature using electron spin resonance (ESR) in several polar solvents; they are ethyl alcohol, eugenol, dimethyl phthalate, tri-n-butyl phosphate, tris(2-ethyl-hexyl)phosphate, diglyme, and a dimethylformamide-chloroform mixed solvent. Calculated spectra in agreement with the experimental X-band spectra are obtained using axially symmetric reorientation when the long in-plane axis is the unique (parallel) axis of the rotational diffusion tensor with D parallel/D perpendicular = 3.0-4.0; D parallel and D perpendicular are the diffusion constants for reorientation about the parallel and perpendicular axes, respectively. The reorientational model required for the simulations is either in or close to the Brownian rotational diffusion limit. In the slow motional (low temperature) region, the spectra can be simulated using the glassy g values. As the temperature increases, however, agreement is obtained only if the intermediate g factor, g(y), for the non-axially symmetric Zeeman interaction increases while g(x), g(z), and the motional model remain unchanged; this scheme and others for which gx and g(z) are possibly temperature-dependent are discussed. The values of D perpendicular from the simulations are in general agreement with those from earlier analyses of the width of the central spectral feature. The simulations and width analyses indicate (as do electrochemical, conductivity, and vapor-phase osmometry data) that the paramagnetic species reorienting in solution has a shape similar to that of the Ni(mnt)(2)(-) ion.     

Optical detection of spin multipole order in the ground state of alkali atoms

Using low light intensities, only spin multipole order of rank k = 1 and 2 (orientation and alignment) can be created directly by optical pumping with σ ⁺ and π polarized light. In order to create and observe higher rank multipole order (k > 2) under weak pumping conditions, however, optical pumping and radio-frequency multiple pulse techniques can be applied. We observed all allowed ranks of multipole order for a given spin system in the single-quantum Zeeman spectra of ground state alkali atoms by resorting to the technique of multiple pulse radio frequency-optical double resonance. The theoretical approach to model the experimental results is based on irreducible tensor operators describing both the atomic states and the light fields involved, and results are in good agreement with the experimental observations for Rb and Cs atoms.     

Simulating electron spin resonance spectra of nitroxide spin labels from molecular dynamics and stochastic trajectories

Simulating electron spin resonance spectra of nitroxide spin labels from motional models is necessary for the quantitative analysis of experimental spectra. We present a framework for modeling the spin label dynamics by using trajectories such as those from molecular dynamics (MD) simulations combined with stochastic treatment of the global protein tumbling. This is achieved in the time domain after two efficient numerical integrators are developed: One for the quantal dynamics of the spins and the other for the classical rotational diffusion. For the quantal dynamics, we propagate the relevant part of the spin density matrix in Hilbert space. For the diffusional tumbling, we work with quaternions, which enables the treatment of anisotropic diffusion in a potential expanded as a sum of spherical harmonics. Time-averaging arguments are invoked to bridge the gap between the smaller time step of the MD trajectories and the larger time steps appropriate for the rotational diffusion and/or quantal spin dynamics.     

A general approach for prediction of motional EPR spectra from Molecular Dynamics (MD) simulations: application to spin labelled protein

A general approach for the prediction of EPR spectra directly and completely from single dynamical trajectories generated from Molecular Dynamics (MD) simulations is described. The approach is applicable to an arbitrary system of electron and nuclear spins described by a general form of the spin-Hamiltonian for the entire motional range. It is shown that for a reliable simulation of motional EPR spectra only a single truncated dynamical trajectory generated until the point when correlation functions of rotational dynamics are completely relaxed is required. The simulation algorithm is based on a combination of the propagation of the spin density matrix in the Liouville space for this initial time interval and the use of well defined parameters calculated entirely from the dynamical trajectory for prediction of the evolution of the spin density matrix at longer times. A new approach is illustrated with the application to a nitroxide spin label MTSL attached to the protein sperm whale myoglobin. It is shown that simulation of the EPR spectrum, which is in excellent agreement with experiment, can be achieved from a single MD trajectory. Calculations reveal the complex nature of the dynamics of a spin label which is a superposition of the fast librational motions within dihedral states, of slow rotameric dynamics among different conformational states of the nitroxide tether and of the slow rotational diffusion of the protein itself. The significance of the slow rotameric dynamics of the nitroxide tether on the overall shape of the EPR spectrum is analysed and discussed.     

Rotational spectra of rare isotopic species of bromofluoromethane: determination of the equilibrium structure from ab initio calculations and microwave spectroscopy

Guided by theoretical predictions, the rotational spectra of the mono- and bideuterated species of bromofluoromethane, CDH(79)BrF, CDH(81)BrF, CD(2) (79)BrF, and CD(2) (81)BrF, have been recorded for the first time. Assignment of a few hundred rotational transitions led to the accurate determination of the ground-state rotational constants, all of the quartic and most of the sextic centrifugal distortion constants, as well as the full bromine quadrupole-coupling tensor for both (79)Br and (81)Br, in good agreement with corresponding theoretical predictions based on high-level coupled-cluster calculations. The rotational spectra of the (13)C containing species (13)CH(2) (79)BrF and (13)CH(2) (81)BrF have been observed in natural abundance and have been assigned, thus allowing the determination of the rotational and centrifugal distortion constants as well as the bromine quadrupole-coupling tensor. Furthermore, empirical equilibrium structures have been obtained within a least-squares fit procedure using the available experimental ground-state rotational constants for various isotopic species. Vibrational effects have been accounted for in the analysis using vibration-rotation interaction constants derived from anharmonic force fields computed at the second-order Moller-Plesset perturbation theory as well as coupled-cluster (CC) levels. The empirical equilibrium geometries obtained in this way agree well with the corresponding theoretical predictions obtained from CC calculations [at the CCSD(T) level] after extrapolation to the complete basis set limit and inclusion of core-valence correlation corrections and relativistic effects.     

Carbon-13 NMR relaxation study of 1,8-bis(dimethylamino)naphthalene in isotropic solution

Carbon-13 nuclear spin relaxation in 1,8-bis(dimethylamino)naphthalene (DMAN) was investigated in a dimethylformamide- d 7 solution. In addition, the chemical shielding tensors were measured in the crystalline powder. Detailed analysis of (13)C longitudinal relaxation in this molecule yielded its rotational diffusion tensor. Comparison to the protonated form of DMAN, DMANH(+), leads to conclusions concerning interaction of the latter with its counterion.     

Carbon-13 spin–lattice relaxation of tropane alkaloids. Anisotropic rotational motion of tropine and pseudotropine

The ¹³C spin–lattice relaxation times of tropine and pseudotropine have been measured in CDCl₃ as a function of concentration. The same relative increase in concentration serves to increase the relaxation rates much less in the region 0.9–5.0 wt.% than in the region 5.0–14.3 wt.%. The rotational diffusion coefficients have been calculated from the relaxation data using Woessner's anisotropic rotational diffusion model. Reorientation of both molecules is shown to be moderately anisotropic. The principal axes of the rotational diffusion tensor in the symmetry plane of both molecules are rotationally shifted from the principal axes of the moment of inertia tensor of the free molecules, and the main rotational axis is parallel with a line passing through the centre of mass of the molecule and the nitrogen atom.     

Cu[CC6H11O)2PS2]2配合物单晶的电子顺磁共振

本文报道了Cu[(C6H11O)2PS2]2配合物单晶在X波段室温下的电子顺磁共振研究. 电子顺磁共振谱显示出由^6^3Cu和^6^5Cu的磁性核引起的超精细结构以及由配体^3^1P的磁性核引起的配体超精细结构. 用非同轴的g张量和A张量系统的最小二乘拟合技术, 严格地计算了自旋Hamiltonian参数. g张量的主值表征, Cu^2^+处在由四个配体S形成的平行四方形的中心, 具有四角对称性, 但是由于配体中两个P的影响, 在CuS4平面上A张量的主值出现较大的各向异性. g张量和A张量有一个主轴是共轴的, 它们与CuS4平面垂直. 实验上观察到电子自旋与配体中^3^1P的相互作用是各向同性的, 并获得相应的配体超精细耦合常数A^p值.     

Ab initio and DFT studies of the spin-orbit and spin-spin contributions to the zero-field splitting tensors of triplet nitrenes with aryl scaffolds

Spin-orbit and spin-spin contributions to the zero-field splitting (ZFS) tensors (D tensors) of spin-triplet phenyl-, naphthyl-, and anthryl-nitrenes in their ground state are investigated by quantum chemical calculations, focusing on the effects of the ring size and substituted position of nitrene on the D tensor. A hybrid CASSCF/MRMP2 approach to the spin-orbit term of the D tensor (D(SO) tensor), which was recently proposed by us, has shown that the spin-orbit contribution to the entire D value, termed the ZFS parameter or fine-structure constant, is about 10% in all the arylnitrenes under study and less depends on the size and connectivity of the aryl groups. Order of the absolute values for D(SO) can be explained by the perturbation on the energy level and spatial distributions of π-SOMO through the orbital interaction between SOMO of the nitrene moiety and frontier orbitals of the aryl scaffolds. Spin-spin contribution to the D tensor (D(SS) tensor) has been calculated in terms of the McWeeny-Mizuno equation with the DFT/EPR-II spin densities. The D(SS) value calculated with the RO-B3LYP spin density agrees well with the D(Exptl) -D(SO) reference value in phenylnitrene, but agreement with the reference value gradually becomes worse as the D value decreases. Exchange-correlation functional dependence on the D(SS) tensor has been explored with standard 23 exchange-correlation functionals in both RO- and U-DFT methodologies, and the RO-HCTH/407 method gives the best agreement with the D(Exptl) -D(SO) reference value. Significant exchange-correlation functional dependence is observed in spin-delocalized systems such as 9-anthrylnitrene (6). By employing the hybrid CASSCF/MRMP2 approach and the McWeeny-Mizuno equation combined with the RO-HCTH/407/EPR-II//U-HCTH/407/6-31G* spin densities for D(SO) and D(SS), respectively, a quantitative agreement with the experiment is achieved with errors less than 10% in all the arylnitrenes under study. Guidelines to the putative approaches to D(SS) tensor calculations are given.     

An investigation of the segmental flexibility of spin labelled synthetic macromolecules in solutions

Within the framework of a model, which takes into account fast anisotropic rotation of the spin label and slow diffusion motion of the polymer chain segment, the ESR spectra of 1% solutions of spin labelled poly(4-vinylpyridine) have been analysed. It has been shown that the model describes the observed spectra satisfactorily. The rotational correlation times for the segment together with its effective dimensions (～ 12 monomer units) have been calculated. The dependence of segment correlation time on temperature is described by the Arrhenius equation with E = 4 kcal/mol. τ_o = 2·10^?13 sec.     

A 13C-Labeled Triarylmethyl Radical as an EPR Spin Probe Highly Sensitive to Molecular Tumbling

A stable triarylmethyl spin probe whose electron paramagnetic resonance (EPR) spectrum is highly sensitive to molecular tumbling is reported. The strong anisotropy of the hyperfine coupling tensor with the central carbon of a ¹³C₁-labeled triarylmethyl radical enables the measurement of the probe rotational correlation time with applications to measure microviscosity and molecular dynamics.     

A 13C‐Labeled Triarylmethyl Radical as an EPR Spin Probe Highly Sensitive to Molecular Tumbling

A stable triarylmethyl spin probe whose electron paramagnetic resonance (EPR) spectrum is highly sensitive to molecular tumbling is reported. The strong anisotropy of the hyperfine coupling tensor with the central carbon of a ¹³C₁‐labeled triarylmethyl radical enables the measurement of the probe rotational correlation time with applications to measure microviscosity and molecular dynamics.