Spin dynamics calculations of electron and nuclear spin relaxation times in paramagnetic solutions

Nuclear spin relaxation in aqueous paramagnetic ion solutions

A Brownian shell model describing the random rotational motion of a spherical shell of uniform particle density is presented and validated by molecular dynamics simulations. The model is applied to proton spin rotation in aqueous paramagnetic ion complexes to yield an expression for the Larmor-frequency-dependent nuclear magnetic resonance spin-lattice relaxation rate T_{1}^{-1}(ω) describing the dipolar coupling of the nuclear spin of the proton with the electronic spin of the ion. The Brownian shell model provides a significant enhancement to existing particle-particle dipolar models without added complexity, allowing fits to experimental T_{1}^{-1}(ω) dispersion curves without arbitrary scaling parameters. The model is successfully applied to measurements of T_{1}^{-1}(ω) from aqueous manganese(II), iron(III), and copper(II) systems where the scalar coupling contribution is known to be small. Appropriate combinations of Brownian shell and translational diffusion models, representing the inner and outer sphere relaxation contributions, respectively, are shown to provide excellent fits. Quantitative fits are obtained to the full dispersion curve of each aquoion with just five fit parameters, with the distance and time parameters each taking a physically justifiable numerical value.

A method to calculate the NMR spectra of paramagnetic species using thermalized electronic relaxation

For paramagnetic species, it has been long understood that the hyperfine interaction between the unpaired electrons and the nucleus results in a nuclear magnetic resonance (NMR) peak that is shifted by a paramagnetic shift, rather than split by the coupling, due to an averaging of the electronic magnetic moment caused by electronic relaxation that is fast in comparison to the hyperfine coupling constant. However, although this feature of paramagnetic NMR has formed the basis of all theories of the paramagnetic shift, the precise theory and mechanism of the electronic relaxation required to predict this result has never been discussed, nor has the assertion been tested. In this paper, we show that the standard semi-classical Redfield theory of relaxation fails to predict a paramagnetic shift, as does any attempt to correct for the semi-classical theory using modifications such as the inhomogeneous master equation or Levitt-di Bari thermalization. In fact, only the recently-introduced Lindbladian theory of relaxation in magnetic resonance [J.Magn.Reson., 310, 106645 (2019)] is able to correctly predict the paramagnetic shift tensor and relaxation-induced linewidth in pNMR. Furthermore, this new formalism is able to predict the NMR spectra of paramagnetic species outside the high-temperature and weak-order limits, and is therefore also applicable to dynamic nuclear polarization. The formalism is tested by simulations of five case studies, which include Fermi-contact and spin-dipolar hyperfine couplings, g-anisotropy, zero-field splitting, high and low temperatures, and fast and slow electronic relaxation.

Novel blue-green emitting NLOphoric triphenylamine-imidazole based donor-π-acceptor compound: Solvatochromism, DFT, TD-DFT and non-linear optical studies

Novel Donor (D)-π-Acceptor (A) NLOphoric triphenylamine-imidazole based dye 9 was designed, synthesized, and confirmed by Mass, ¹³C NMR, and ¹H NMR analysis. Photophysical properties of 9 were studied in solvents of different polarities and compared with analogues compounds 7 and 8. Phenonthroline acceptor based dye 9 shows highly bathochromic shifted absorption and emission compared to dyes 7 and 8. Positive solvatochromism was noticed in 7, 8, and 9 which was supported by the linear (i.e. Lippert-Mataga and Mac-Rae polarity functions) and multi-linear (i.e. Kamlet-Taft and Catalan parameters) analysis. Moreover, solvent polarizability (d_SP) and solvent dipolarity (C_SdP) are the major factors responsible for red shift in absorption as well as in emission spectra. Charge transfer descriptors as well as the polarity graphs are in good relation with Generalized Mulliken-Hush (GMH) parameters. NLO properties of 7, 8, and 9 were studied by using solvatochromic and computational methods. The static first hyperpolarizability (β₀) and relevant microscopic parameters (μ,α₀,α,β,γ) were determined using DFT with B3LYP, BHHLYP, and CAM-B3LYP functionals. Third-order NLO properties of nitrogen containing phenanthroline based compound 9 were observed to be several times higher than those of the compounds 7 and 8, justify the design approach.

Paramagnetic NMR in solution and the solid state

The field of paramagnetic NMR has expanded considerably in recent years. This review addresses both the theoretical description of paramagnetic NMR, and the way in which it is currently practised. We provide a review of the theory of the NMR parameters of systems in both solution and the solid state. Here we unify the different languages used by the NMR, EPR, quantum chemistry/DFT, and magnetism communities to provide a comprehensive and coherent theoretical description. We cover the theory of the paramagnetic shift and shift anisotropy in solution both in the traditional formalism in terms of the magnetic susceptibility tensor, and using a more modern formalism employing the relevant EPR parameters, such as are used in first-principles calculations. In addition we examine the theory first in the simple non-relativistic picture, and then in the presence of spin-orbit coupling. These ideas are then extended to a description of the paramagnetic shift in periodic solids, where it is necessary to include the bulk magnetic properties, such as magnetic ordering at low temperatures. The description of the paramagnetic shift is completed by describing the current understanding of such shifts due to lanthanide and actinide ions. We then examine the paramagnetic relaxation enhancement, using a simple model employing a phenomenological picture of the electronic relaxation, and again using a more complex state-of-the-art theory which incorporates electronic relaxation explicitly. An additional important consideration in the solid state is the impact of bulk magnetic susceptibility effects on the form of the spectrum, where we include some ideas from the field of classical electrodynamics. We then continue by describing in detail the solution and solid-state NMR methods that have been deployed in the study of paramagnetic systems in chemistry, biology, and the materials sciences. Finally we describe a number of case studies in paramagnetic NMR that have been specifically chosen to highlight how the theory in part one, and the methods in part two, can be used in practice. The systems chosen include small organometallic complexes in solution, solid battery electrode materials, metalloproteins in both solution and the solid state, systems containing lanthanide ions, and multi-component materials used in pharmaceutical controlled-release formulations that have been doped with paramagnetic species to measure the component domain sizes.

Investigation of NLO Properties of Fluorescent BORICO Dyes: a Comprehensive Experimental and Theoretical Approach

BORICO dyes with N, N-diethyl as a strong donor and BF₂ complexed iminocoumarin six member core as strong acceptor are investigated as an efficient non linear optical chromophores. Extended π-conjugation over iminocoumarin moiety is useful to make ICT character of BORICO dyes more significant and is established on the scale of Generalised Mulliken Hush analysis scale. Bond length alternation and bond order alternation values for three BORICO chromophores estimates the cyanine like framework for optimal non linear optical response. The frontier molecular orbital diagrams obtained from density functional theory calculations shows that there is charge transfer from donor to accepter as well as effective overlap between them making the basis for optimal NLO response of BORICO chromophores. The theoretical values of linear and non linear optical responses for three BORICO NLOphores obtained by using three different functionals B3LYP, CAMB3LYP and BHandHLYP with 6-311+g(d,p) basis set are quite consistent for the values of static dipole moment (μ), linear polarizability (α) and first hyperpolarizability (β). However in case of the γ values calculation, compare to the similar values obtained by CAMB3LYP and BHandHLYP functionals, B3LYP overestimates the same. The vibrational motions play decisive role in the overall non linear optical properties of BORICO chromophores.

Transient versus Static Electron Spin Relaxation in Mn(2+) Complexes Relevant as MRI Contrast Agents

The zero-field splitting (ZFS) parameters of the [Mn(EDTA)(H2O)](2-)·2H2O and [Mn(MeNO2A)(H2O)]·2H2O systems were estimated by using DFT and ab initio CASSCF/NEVPT2 calculations (EDTA = 2,2',2″,2‴-(ethane-1,2-diylbis(azanetriyl))tetraacetate; MeNO2A = 2,2'-(7-methyl-1,4,7-triazonane-1,4-diyl)diacetate). Subsequent molecular dynamics calculations performed within the atom-centered density matrix propagation (ADMP) approach provided access to the transient and static ZFS parameters, as well as to the correlation time of the transient ZFS. The calculated ZFS parameters present a reasonable agreement with the experimental values obtained from the analysis of (1)H relaxation data. The correlation times calculated for the two systems investigated turned out to be very short (τc ∼ 0.02-0.05 ps), which shows that the transient ZFS is modulated by molecular vibrations. On the contrary, the static ZFS is modulated by the rotation of the complexes in solution, which for the small complexes investigated here is characterized by rotational correlation times of τR ∼ 35-60 ps. As a result, electron spin relaxation in small Mn(2+) complexes is dominated by the static ZFS.

Nuclear magnetic resonance relaxivities: investigations of ultrahigh-spin lanthanide clusters from 10 MHz to 1.4 GHz

Paramagnetic relaxation enhancement is often explored in magnetic resonance imaging in terms of contrast agents and in biomolecular nuclear magnetic resonance (NMR) spectroscopy for structure determination. New ultrahigh-spin clusters are investigated with respect to their NMR relaxation properties. As their molecular size and therefore motional correlation times as well as their electronic properties differ significantly from those of conventional contrast agents, questions about a comprehensive characterization arise. The relaxivity was studied by field-dependent longitudinal and transverse NMR relaxometry of aqueous solutions containing Fe(III)(10)Dy(III)(10) ultrahigh-spin clusters (spin ground state 100/2). The high-field limit was extended to 32.9 T by using a 24 MW resistive magnet and an ultrahigh-frequency NMR setup. Interesting relaxation dispersions were observed; the relaxivities increase up to the highest available fields, which indicates a complex interplay of electronic and molecular correlation times.

Jahn-Teller dynamics in a series of high-symmetry Co(II) chelates determine paramagnetic relaxation enhancements

NMR paramagnetic relaxation enhancements (PREs) of a series of structurally characterized, trigonal bis-trispyrazolylborate (Tp) chelates of high-spin Co(II), spanning 100-850 MHz in field, are reported. Prior knowledge of the metal-nucleus distances allows numerical extraction of position-dependent electron spin relaxation rates (τ(c)(-1)) from direct measurement of the individual PREs of the four symmetry distinct protons in Co(Tp)(2), using available closed-form expressions. The data for this electronically complex system where spin-orbit coupling defines the ground state electronic structure are analyzed in terms of the Solomon-Bloembergen-Morgan (SBM) relations, as well as available zero-field splitting limit theories. A simple angular correction is shown to be sufficient to reconcile the individual τ(c)(T) data for the four classes of protons. The data identify a previously unrecognized dynamic Jahn-Teller effect in these historically important complexes, with a barrier of ~230 cm(-1), pointing to a level of dynamics in trispyrazolylborate chemistry that has not been described before, and further show that it is the Jahn-Teller that is responsible for the PREs in fluid solution. A field-dependent component is also identified for the two protons nearest g(//), which is suggested to arise due to Zeeman mixing of excited state character into the ground level.

Nuclear quadrupole resonance lineshape analysis for different motional models: stochastic Liouville approach

A general theory of lineshapes in nuclear quadrupole resonance (NQR), based on the stochastic Liouville equation, is presented. The description is valid for arbitrary motional conditions (particularly beyond the valid range of perturbation approaches) and interaction strengths. It can be applied to the computation of NQR spectra for any spin quantum number and for any applied magnetic field. The treatment presented here is an adaptation of the "Swedish slow motion theory," [T. Nilsson and J. Kowalewski, J. Magn. Reson. 146, 345 (2000)] originally formulated for paramagnetic systems, to NQR spectral analysis. The description is formulated for simple (Brownian) diffusion, free diffusion, and jump diffusion models. The two latter models account for molecular cooperativity effects in dense systems (such as liquids of high viscosity or molecular glasses). The sensitivity of NQR slow motion spectra to the mechanism of the motional processes modulating the nuclear quadrupole interaction is discussed.

Implicit numerical schemes for the stochastic Liouville equation in Langevin form

We present and numerically test implicit as well as explicit numerical schemes for solving the Stochastic Liouville Equation in Langevin form. It is found that implicit schemes provide significant gain in robustness, for example, when nonsecular Hamiltonian terms cannot be ignored in electron and nuclear spin resonance. Implicit schemes open up several spectroscopic relaxation problems for direct interpretation using the Stochastic Liouville Equation. To illustrate the proposed numerical schemes, studies are presented for an electron paramagnetic resonance problem involving a coordinated copper complex and a fluorescence problem.

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.

Joint analysis of ESR lineshapes and 1H NMRD profiles of DOTA-Gd derivatives by means of the slow motion theory

The "Swedish slow motion theory" [Nilsson and Kowalewski, J. Magn. Reson. 146, 345 (2000)] applied so far to Nuclear Magnetic Relaxation Dispersion (NMRD) profiles for solutions of transition metal ion complexes has been extended to ESR spectral analysis, including in addition g-tensor anisotropy effects. The extended theory has been applied to interpret in a consistent way (within one set of parameters) NMRD profiles and ESR spectra at 95 and 237 GHz for two Gd(III) complexes denoted as P760 and P792 (hydrophilic derivatives of DOTA-Gd, with molecular masses of 5.6 and 6.5 kDa, respectively). The goal is to verify the applicability of the commonly used pseudorotational model of the transient zero field splitting (ZFS). According to this model the transient ZFS is described by a tensor of a constant amplitude, defined in its own principal axes system, which changes its orientation with respect to the laboratory frame according to the isotropic diffusion equation with a characteristic time constant (correlation time) reflecting the time scale of the distortional motion. This unified interpretation of the ESR and NMRD leads to reasonable agreement with the experimental data, indicating that the pseudorotational model indeed captures the essential features of the electron spin dynamics.

General treatment of paramagnetic relaxation enhancement associated with translational diffusion

A theory of nuclear spin relaxation in isotropic liquids for nuclear spins interacting with electron spins, residing in other molecules (the outer-sphere relaxation), is presented. The approach, valid outside of the Redfield limit for electron spin relaxation, is an extension of the Swedish slow motion theory [Benetis et al., Mol. Phys. 48, 329 (1983); Nilsson and Kowalewski, J. Magn. Reson. 146, 345 (2000)] for inner-sphere relaxation. It is demonstrated that the outer-sphere relaxation rate can be expressed as an integral of a product of a translational diffusion correlation function and a function analogous to the inner-sphere spectral density. A numerical implementation of the theory is described and applied to a large number of realistic parameter sets for S = 7/2 and S = 1, which may correspond to Gd(III) and Ni(II) systems. It is shown that the outer-sphere contribution is relevant and should be included into the analysis of nuclear magnetic relaxation dispersion relaxation profiles, especially for slow relative translational diffusion and fast molecular tumbling.

Relaxation theory of the electronic spin of a complexed paramagnetic metal ion in solution beyond the Redfield limit

The relaxation of the electronic spin S of a paramagnetic metal ion with fully quenched orbital angular momentum in its ground state is investigated in an external magnetic field through a systematic study of the time correlation functions governing the evolution of the statistical operator (density matrix). Let omega0 be the Larmor angular frequency of S. When the relaxation is induced by a time-fluctuating perturbing Hamiltonian hH1(t) of time correlation tauc, it is demonstrated that after a transient period the standard Redfield approximation is relevant to calculate the evolution of the populations of the spin states if parallelH1 parallel2tauc2/(1+omega0(2)tauc2)<<1 and that this transient period becomes shorter than tauc at sufficiently high field for a zero-field splitting perturbing Hamiltonian. This property, proven analytically and confirmed by numerical simulation, explains the surprising success of several simple expressions of the longitudinal electronic relaxation rate 1/T1e derived from the Redfield approximation well beyond its expected validity range parallelH1 paralleltauc<<1. It has favorable practical consequences on the interpretation of the paramagnetic relaxation enhancement of nuclei used for structural and dynamic studies.

Design of Gd(III)-Based Magnetic Resonance Imaging Contrast Agents: Static and Transient Zero-Field Splitting Contributions to the Electronic Relaxation and Their Impact on Relaxivity

A multiple-frequency (9.4-325 GHz) and variable-temperature (276-320 K) electron paramagnetic resonance (EPR) study on low molecular weight gadolinium(III) complexes for potential use as magnetic resonance imaging (MRI) contrast agents has been performed. Peak-to-peak linewidths Delta Hpp and central magnetic fields have been analyzed within the Redfield approximation taking into account the static zero-field splitting (ZFS) up to the sixth order and the transient ZFS up to the second order. Longitudinal electronic relaxation is dominated by the static ZFS contribution at low magnetic fields (B < 0.3 T) and by the transient ZFS at high magnetic fields (B > 1.5 T). Whereas the static ZFS clearly depends on the nature of the chelating ligand, the transient ZFS does not. For the relatively fast rotating molecules studied water proton relaxivity is mainly limited by the fast rotation and electronic relaxation has only a marked influence at frequencies below 30 MHz. From our EPR results we can conclude that electronic relaxation will have no influence on the efficiency of Gd(III)-based MRI contrast agents designed for studies at very high magnetic fields (B > 3T).

Direct calculation of1H2O T1NMRD profiles and EPR lineshapes for the electron spin quantum numbers S = 1, 3/2, 2, 5/2, 3, 7/2, based on the stochastic Liouville equation combined with Brownian dynamics simulation

Direct calculation of electron spin relaxation and EPR lineshapes, based on Brownian dynamics simulation techniques and the stochastic Liouville equation approach (SLE-L) [Mol. Phys., 2004, 102, 1085-1093], is here generalized to high spin systems with spin quantum number S = 3/2, 2, 5/2, 3 and 7/2. A direct calculation method is demonstrated for electron spin-spin and spin-lattice relaxation, S-, X- and Q-band EPR-lineshapes and paramagnetic enhanced water proton T(1)- NMRD profiles. The main relaxation mechanism for the electron spin system is a stochastic second rank zero field splitting (ZFS). Brownian dynamics simulation techniques are used in describing a fluctuating ZFS interaction which comprises two parts namely the "permanent" part which is modulated by isotropic reorientation diffusion, and the transient part which is modulated by fast local distortion, which is also modelled by the isotropic rotation diffusion model. The SLE-L approach present is applicable both in the perturbation (Redfield) regime as well as outside the perturbation regime, in the so called slow motion regime.

Four complementary theoretical approaches for the analysis of NMR paramagnetic relaxation

Four theoretical and computational approaches used at the University of Michigan to analyze NMR paramagnetic relaxation enhancement (NMR-PRE) are described. The primary objective of the theory is to describe the relationship of the NMR-PRE phenomenon to the electron spin hamiltonian and the spin energy level structure when zero field splitting interactions are significant. Four formulations of theory are discussed: (1) spin dynamics simulation; (2) the laboratory frame "constant H(S)" formulation; (3) the Molecular Frame "constant H(S)" formulation; and (4) the zfs-limit "constant H(S)" formulation. No single theoretical approach describes all important aspects of the relaxation mechanism in a fully satisfactory way. We use the four formulations in a complementary manner to provide as complete a picture of the relaxation mechanism as possible. We also discuss the integration of NMR-PRE theory and recently developed theory of electron spin relaxation which accounts for effects of the permanent zfs hamiltonian.

Electronic relaxation of paramagnetic metal ions and NMR relaxivity in solution: Critical analysis of various approaches and application to a Gd(III)-based contrast agent

The time correlation functions (TCFs) G(alphaalpha(t)[triple bond](Salpha(t)Salpha(0)) (alpha = x,y,z) of the electronic spin components of a complexed paramagnetic metal ion give information about the time fluctuations of its zero-field splitting (ZFS) Hamiltonian due to the random dynamics of the coordination polyhedron. These TCFs reflect the electronic spin relaxation which plays an essential role in the inner- and outer-sphere paramagnetic relaxation enhancements of the various nuclear spins in solution. When a static ZFS Hamiltonian is allowed by symmetry, its modulation by the random rotational motion of the complex has a great influence on the TCFs. We discuss several attempts to describe this mechanism and show that subtle mathematical pitfalls should be avoided in order to obtain a theoretical framework, within which reliable adjustable parameters can be fitted through the interpretation of nuclear-magnetic relaxation dispersion experimental results. We underline the advantage of the numerical simulation of the TCFs, which avoids the above difficulties and allows one to include the effect of the transient ZFS for all the relative magnitudes of the various terms in the electron-spin Hamiltonian and arbitrary correlation times. This method is applied for various values of the magnetic field taken to be along the z direction. At low field, contrary to previous theoretical expectations, if the transient ZFS has negligible influence, the longitudinal TCF GII(t) [triple bond] G(zz)(t) has a monoexponential decay with an electronic relaxation time T1e different from 1/(2D(r)), D(r) being the rotational diffusion coefficient of the complex. At intermediate and high field, the simulation results show that GII (t) still has a monoexponential decay with a characteristic time T1e, which is surprisingly well approximated by a simple analytical expression derived from the Redfield perturbation approximation of the time-independent Zeeman Hamiltonian, even in the case of a strong ZFS where this approximation is expected to fail. These results are illustrated for spins S = 1, 3/2, and 5/2 in axial and rhombic symmetries. Finally, the simulation method is applied to the reinterpretation of the water-proton relaxivity profile due to P760-Gd(III), an efficient blood pool contrast agent for magnetic-resonance imaging.

NMR paramagnetic relaxation due to the S=5∕2 complex, Fe(III)-(tetra-p-sulfonatophenyl)porphyrin: Central role of the tetragonal fourth-order zero-field splitting interaction

The metalloporphyrins, Me-TSPP [Me=Cr(III), Mn(III), Mn(II), Fe(III), and TSPP=meso-(tetra-p-sulfonatophenyl)porphyrin], which possess electron spins S=3/2, 2, 5/2, and 5/2, respectively, comprise an important series of model systems for mechanistic studies of NMR paramagnetic relaxation enhancement (NMR-PRE). For these S>1/2 spin systems, the NMR-PRE depends critically on the detailed form of the zero-field splitting (zfs) tensor. We report the results of experimental and theoretical studies of the NMR relaxation mechanism associated with Fe(III)-TSPP, a spin 5/2 complex for which the overall zfs is relatively large (D approximately = 10 cm(-1)). A comparison of experimental data with spin dynamics simulations shows that the primary determinant of the shape of the magnetic relaxation dispersion profile of the water proton R1 is the tetragonal fourth-order component of the zfs tensor. The relaxation mechanism, which has not previously been described, is a consequence of zfs-induced mixing of the spin eigenfunctions of adjacent Kramers doublets. We have also investigated the magnetic-field dependence of electron-spin relaxation for S=5/2 in the presence of a large zfs, such as occurs in Fe(III)-TSPP. Calculations show that field dependence of this kind is suppressed in the vicinity of the zfs limit, in agreement with observation.

NMR-Paramagnetic Relaxation Due to the High-Spin d3 Electron Configuration: Cr(III)−TSPP

Sulfonated metalloporphyrins (Me-TSPP, where Me = Cr(III), Mn(III), Fe(III), and Mn(II)) comprise a well-characterized series of water-soluble paramagnetic complexes with electron spins of S = 3/2, 2, 5/2, and 5/2, respectively, which provide important model systems for mechanistic studies of paramagnetic NMR relaxation in solution. Previous studies of Mn(III), Fe(III), and Mn(II)-TSPP have uncovered relaxation mechanisms which differ qualitatively from each other and exhibit numerous unexpected features. In this study, Cr(III)-TSPP was examined as a model system for the d3 S = 3/2 electron configuration. Magnetic relaxation dispersion (MRD) profiles of the water proton R1 were measured as a function of pH between pH 1 and pH 9. In acid samples, R1 results from acid-catalyzed prototropic chemical exchange involving the Cr(III)-TSPP x 2 H2O. In neutral and basic solution, this species deprotonates, and base-catalyzed prototropic exchange becomes important. The pH 1 data were analyzed quantitatively using theory that accounts for the role of the permanent zero field splitting (zfs) tensor and for the effects of Brownian reorientation. Two levels of theory were employed: (1) spin dynamics simulation, which accurately describes the effects of Brownian reorientation on the spin wave functions, and (2) the "constant H(S)" approximation, which incorporates the effects of multiexponential electron spin relaxation and facilitates the physical interpretation of the relaxation mechanism. It was found that neither level of theory alone provides a fully satisfactory quantitative description of the data due to the fact that both reorientational modulation of the spin wave functions and multiexponential electron spin relaxation are important. The zero field splitting parameter, D = 0.27 cm(-1), is well defined by the data and was measured.

Electron spin relaxation due to reorientation of a permanent zero field splitting tensor

Electron spin relaxation of transition metal ions with spin S> or =1 results primarily from thermal modulation of the zero field splitting (zfs) tensor. This occurs both by distortion of the zfs tensor due to intermolecular collisions and, for complexes with less than cubic symmetry, by reorientational modulation of the permanent zfs tensor. The reorientational mechanism is much less well characterized in previous work than the distortional mechanism although it is an important determinant of nuclear magnetic resonance (NMR) paramagnetic relaxation enhancement phenomena (i.e., the enhancement of NMR relaxation rates produced by paramagnetic ions in solution or NMR-PRE). The classical density matrix theory of spin relaxation does not provide an appropriate description of the reorientational mechanism at low Zeeman field strengths because the zero-order spin wave functions are stochastic functions of time. Using spin dynamics simulation techniques, the time correlation functions of the spin operators have been computed and used to determine decay times for the reorientational relaxation mechanism for S=1. In the zfs limit of laboratory field strengths (H(Zeem)<<H(zfs) (composite function)), when the zfs tensor is cylindrical, the spin decay is exponential, the spin relaxation time, tau(S) (composite function) approximately 0.53tau(R)((1)), where tau(R)((1)) is the reorientational correlation time of a molecule-fixed vector. The value of tau(S) (composite function) is independent of the magnitude of the cylindrical zfs parameter (D), but it depends strongly on low symmetry zfs terms (the E/D ratio). Other spin dynamics (SD) simulations examined spin decay in the intermediate regime of field strengths where H(Zeem) approximately H(zfs) (composite function), and in the vicinity of the Zeeman limit. The results demonstrate that the reorientational electron spin relaxation mechanism is often significant when H(zfs) (composite function)> or =H(Zeem), and that its neglect can lead to serious errors in the interpretation of NMR-PRE data.

On the philosophy of optimizing contrast agents. An analysis of 1H NMRD profiles and ESR lineshapes of the Gd(III)complex MS-325+HSA

A generalization of the modified SBM theory is developed in closed analytical form. The theory is applied to describe the paramagnetically enhanced water proton spin-lattice relaxation rates of the aqueous-systems containing a gadolinium(S=7/2) complex(MS-325) in the presence or absence of human serum albumin (HSA). MS-325 binds to HSA: in the absence of the protein the reorientational time, tauR, is short, but when HSA is added tauR becomes much longer. In this way, the effect of reorientational motion, static (Delta s), and transient (Delta t) zero-field splitting (ZFS) interactions on both the water proton relaxivity and the Gd ESR lineshapes are investigated. Two dynamic models of electron spin relaxation are presented, characterized by transient and static ZFS-interactions. X-, Q-, and W-bands ESR spectra of MS-325+HSA are analyzed in order to describe the effect on the electron spin system upon binding to a macromolecule. A computer program based on this theory is developed which calculates solvent water proton T1 NMRD profiles and the corresponding X-, Q-, U-, and W-bands ESR lineshapes.

Nuclear spin relaxation in paramagnetic systems (S>/=1) under fast rotation conditions

A new theoretical model for nuclear spin relaxation in paramagnetic systems in solution has been developed. Fast rotational motion is included in the model, both as a source of modulation of the static zero-field splitting, which provides a mechanism for electron spin relaxation, and as an origin of the stochastic variation of the electron spin-nuclear spin dipole-dipole interaction leading to nuclear spin relaxation. At the limit of low magnetic field, the model is essentially identical to the earlier formulations from our laboratory, but new closed-form expressions are given for the inner- and outer-sphere relaxation at the high-field limit. Numerical comparisons with a general theory are reported for the inner-sphere case. In addition, some nuclear magnetic relaxation dispersion (NMRD) profiles from the literature are considered for systems where experiments have been done with both low-molecular weight paramagnetic complexes and their adducts with proteins. Previously developed theories are used to interpret data for the slowly rotating protein adducts, and good fits of the fast-rotating counterparts are obtained by further adjustment of one or two additional parameters.

Multiexponential electronic spin relaxation and Redfield's limit in Gd(III) complexes in solution: consequences for 17O/1H NMR and EPR simultaneous analysis

Multiple experiments (17O NMR, 1H NMR, and EPR) have been performed in the past to understand the microscopic parameters that control the magnetic relaxation rate enhancement induced by paramagnetic molecules on neighboring water protons, the so-called relaxivity. The generally accepted theories of the electron spin relaxation of S = 7/2 ions such as Gd3+ (Solomon-Bloembergen-Morgan or simplified Hudson-Lewis) are unsatisfactory for a simultaneous analysis. Recently, an improved theory, where the electron spin relaxation is due to the combination of a static (thus explicitly linked to the molecular structure) and a dynamic zero field splitting, has been developed and tested on experimental EPR data. The model has also been extended beyond the electronic Redfield limit using Monte Carlo simulations. Using the aqua ion [Gd(H2O)8]3+ as a test case, we present here the first simultaneous analysis of 17O NMR, 1H NMR, and EPR relaxation data using this rigorous approach of the electron spin relaxation. We discuss the physical meaning of the calculated parameters. The consequences on future experiments are also considered, especially regarding the analysis of nuclear magnetic relaxation dispersion (NMRD) profiles in the study of Gd3+ complexes.

Closed-form expressions for level-averaged electron spin relaxation times outside the Zeeman limit: application to paramagnetic NMR relaxation

Paramagnetic enhancement of NMR relaxation (NMR-PRE) depends on thermal relaxation of the electron spin system. Most previous analyses of experimental NMR-PRE data have relied on Bloembergen--Morgan (B--M) theory to describe the magnetic field dependence of electron spin relaxation in liquid samples. However, B--M theory assumes a Zeeman-limit situation and is not physically appropriate to the common case of S > or = 1 transition metal ions which possess a permanent zero-field splitting (zfs) that is comparable to or larger than the Zeeman splitting. Theory has been needed which (1) includes the effects of the zfs interaction, thus providing a realistic description of the magnetic field dependence of the NMR-PRE outside the Zeeman limit, and (2) describes electron spin relaxation phenomena at a comparable level of complexity to that of B--M theory, i.e., with two magnetic field-dependent electron spin relaxation times, tau(S1) and tau(S2), defined in the laboratory coordinate frame. Theory of this kind is developed. Expressions derived in a previous study (R. R. Sharp and L. L. Lohr, J. Chem. Phys. 115, 5005 (2001).) for level-specific relaxation rates have been averaged over spin eigenstates to give level-averaged quantities, tau(S1,2). This kind of averaging leads to a great simplification in the mathematical form of the results. Simple zfs-limit molecular-frame and laboratory-frame expressions are given for electron spin S=1, 3/2, 2, and 5/2. General expressions, valid for S > or = 1 and for arbitrary magnitudes of the Zeeman and zfs energies, are derived for level-averaged electron spin relaxation times defined in both the laboratory- and the molecule-fixed coordinate frames. The new theory coincides with B--M theory in the Zeeman limit.

NMR paramagnetic relaxation enhancement: test of the controlling influence of zfs rhombicity for S = 1

Prior theoretical work has predicted that the NMR paramagnetic relaxation enhancement (NMR-PRE) produced by electron spin S = 1 ions is highly sensitive to orthorhombic terms in the static zero field splitting (zfs) tensor. Zfs orthorhombicity (which implies chemical inequivalence of the three principal directions of the zfs-principal axis system and is described by the zfs E-parameter) is predicted to suppress the NMR-PRE profoundly relative to the reference cylindrical zfs-limit situation. This expectation was tested experimentally by a comparison of the zfs-limit NMR-PRE produced by [Ni(II)(en)(3)](2+) (en = ethylenediamine), a trigonal complex which lacks zfs-rhombicity, with the zfs-limit NMR-PRE produced by two orthorhombic complexes, [Ni(II)(en)(2)(H(2)O)(2)](2+) and [Ni(II)(en)(H(2)O)(4)](2+). As predicted, the zfs-limit NMR-PRE produced by the orthorhombic complexes in the proton resonance of a dioxane probe species in the solvent was strongly suppressed (by factors of approximately 5 and 7, respectively) relative to the comparable measurement on the trigonal complex. The suppression of the NMR-PRE due to the orthorhombic zfs terms is counteracted by an applied Zeeman field, leading to a predicted rise in the NMR-PRE with increasing Zeeman field strength; this rise occurs when the Zeeman energy is comparable to the orthorhombic zfs splitting, 2E. This second prediction of theory was likewise confirmed: the expected rhombicity-induced magnetic field dependence in the NMR-PRE was observed for the orthorhombic complexes but not for the trigonal complex.

Slow-motion theory of nuclear spin relaxation in paramagnetic low-symmetry complexes: A generalization to high electron spin

The slow-motion theory of nuclear spin relaxation in paramagnetic low-symmetry complexes is generalized to comprise arbitrary values of S. We describe the effects of rhombic symmetry in the static zero-field splitting (ZFS) and allow the principal axis system of the static ZFS tensor to deviate from the molecule-fixed frame of the nuclear-electron dipole-dipole tensor. We show nuclear magnetic relaxation dispersion (NMRD) profiles for different illustrative cases, ranging from within the Redfield limit into the slow-motion regime with respect to the electron spin dynamics. We focus on S = 3/2 and compare the effects of symmetry-breaking properties on the paramagnetic relaxation enhancement (PRE) in this case with that of S = 1, which we have treated in a previous paper. We also discuss cases of S = 2, 5/2, 3, and 7/2. One of the main objectives of this investigation, together with the previous papers, is to provide a set of standard calculations using the general slow-motion theory, against which simplified models may be tested. Copyright 2000 Academic Press.