The influence of molecular shape and charge distribution on the molecular approach in solution. NMR dipolar relaxation studies of the many-body effects

The paramagnetic dipolar spin-lattice relaxation of nuclear spins is expressed in terms of the geometrical and dynamical parameters of the molecules and used to analyze the intra- and intermolecular microdynamics. The possible conformational changes of a TEMPOL nitroxide-labelled sugar in deuterated chloroform are studied and then the intermolecular processes are considered. The influence of the discrete polar natures of acetonitrile, chloroform, and methanol on the TEMPOL solvation dynamics is recognized. The importance of the water structure on the collisional approach of two attractive ions is emphasized. Finally, the attractive effects of the dipole-dipole electrostatic forces are shown for the water/TEMPOL pair. Diffusion and jump analytical theories, integral equations of statistical mechanics, and Monte-Carlo simulations of the diffusion are combined to treat the anisotropic interactions and the many-body effects.     

Surface nuclear spin relaxation of 199Hg

We report the results of a detailed study of surface relaxation of 199Hg nuclear spins in paraffin coated cells. From measurements of the magnetic field and temperature dependence of the spin relaxation rates we determine the correlation time for magnetic fluctuations and the surface adsorption energy. The data indicate that surface relaxation is caused by dipolar coupling to paramagnetic sites on the surface. We also observe changes in the spin relaxation rate caused by ultraviolet radiation resonant with the 254 nm transition in Hg.     

Two-dimensional measurements of the solvent structural relaxation dynamics in dipolar solvation

Resonant-pump polarizability response spectroscopy (RP-PORS) is based on an optical heterodyne detected transient grating (OHD-TG) method with an additional resonant pump pulse. In RP-PORS, the resonant pump pulse excites the solute-solvent system and the subsequent relaxation of the solute-solvent system is monitored by the OHD-TG spectroscopy. RP-PORS is shown to be an excellent experimental tool to directly measure the solvent responses in solvation. In the present work, we extended our previous RP-PORS (Park et al., Phys. Chem. Chem. Phys., 2011, 13, 214-223) to measure time-dependent transient solvation polarizability (TSP) spectra with Coumarin153 (C153) in acetonitrile. The time-dependent TSP spectra showed how the different solvent intermolecular modes were involved in different stages of the solvation process. Most importantly, the inertial and diffusive components of the solvent intermolecular modes in solvation were found to be spectrally and temporally well-separated. In a dipolar solvation of C153, high-frequency inertial solvent modes were found to be driven instantaneously and decay on a subpicosecond timescale while low-frequency diffusive solvent modes were induced slowly and decayed on a picosecond timescale. Our present result is the first experimental manifestation of frequency-dependent solvent intermolecular response in a dipolar solvation.     

Molecular properties determined from the relaxation of long-lived spin states

The populations of long-lived spin states, in particular, populations of singlet states that are comprised of antisymmetric combinations of product states, |alpha(I)beta(S)> - |beta(I)alpha(S)>, are characterized by very long lifetimes because the dipole-dipole interaction between the two "active" spins I and S that are involved in such states is inoperative as a relaxation mechanism. The relaxation rate constants of long-lived (singlet) states are therefore determined by the chemical shift anisotropy (CSA) of the active spins and by dipole-dipole interactions with passive spins. For a pair of coupled spins, the singlet-state relaxation rate constants strongly depend on the magnitudes and orientations of the CSA tensors. The relaxation properties of long-lived states therefore reveal new information about molecular symmetry and structure and about spectral density functions that characterize the dynamic behavior.     

Slow relaxation of longitudinal multispin orders in weakly and strongly coupled two‐spin systems

Longitudinal multispin order (LOMO) corresponds to a nonequilibrium population distribution in spin systems that exhibit scalar (J), dipolar, or quadrupolar coupling. We investigated the relaxation of longitudinal two‐spin order (2‐LOMO) in systems that had either weakly or strongly J‐coupled spins. Our results indicated longer relaxation times for the 2‐LOMO state compared with the corresponding longitudinal single‐spin state (1‐LOMO). Accessing nuclear spin states that have relaxation times longer than T₁, without the use of external contrast agents, is potentially useful for in vivo imaging and also for studying systems using dynamically hyperpolarized nuclear spins where longer life times are sought to increase the time available to study (bio)chemical events. Copyright © 2012 John Wiley & Sons, Ltd.     

Dynamic nuclear polarization assisted spin diffusion for the solid effect case

The dynamic nuclear polarization (DNP) process in solids depends on the magnitudes of hyperfine interactions between unpaired electrons and their neighboring (core) nuclei, and on the dipole-dipole interactions between all nuclei in the sample. The polarization enhancement of the bulk nuclei has been typically described in terms of a hyperfine-assisted polarization of a core nucleus by microwave irradiation followed by a dipolar-assisted spin diffusion process in the core-bulk nuclear system. This work presents a theoretical approach for the study of this combined process using a density matrix formalism. In particular, solid effect DNP on a single electron coupled to a nuclear spin system is considered, taking into account the interactions between the spins as well as the main relaxation mechanisms introduced via the electron, nuclear, and cross-relaxation rates. The basic principles of the DNP-assisted spin diffusion mechanism, polarizing the bulk nuclei, are presented, and it is shown that the polarization of the core nuclei and the spin diffusion process should not be treated separately. To emphasize this observation the coherent mechanism driving the pure spin diffusion process is also discussed. In order to demonstrate the effects of the interactions and relaxation mechanisms on the enhancement of the nuclear polarization, model systems of up to ten spins are considered and polarization buildup curves are simulated. A linear chain of spins consisting of a single electron coupled to a core nucleus, which in turn is dipolar coupled to a chain of bulk nuclei, is considered. The interaction and relaxation parameters of this model system were chosen in a way to enable a critical analysis of the polarization enhancement of all nuclei, and are not far from the values of (13)C nuclei in frozen (glassy) organic solutions containing radicals, typically used in DNP at high fields. Results from the simulations are shown, demonstrating the complex dependences of the DNP-assisted spin diffusion process on variations of the relevant parameters. In particular, the effect of the spin lattice relaxation times on the polarization buildup times and the resulting end polarization are discussed, and the quenching of the polarizations by the hyperfine interaction is demonstrated.     

Oxygen accessibility to ribonuclease a: quantitative interpretation of nuclear spin relaxation induced by a freely diffusing paramagnet

The nuclear spin relaxation induced by a freely diffusing paramagnetic center provides a direct measure of intermolecular accessibility. A number of factors are involved in a quantitative interpretation of relaxation data including excluded volume effects, solvation differences, and the details of the electron spin relaxation in the paramagnetic center. In the case where the electron relaxation time is short compared with correlation times describing the electron-nuclear coupling, the nuclear spin relaxation rates may be related to the effective local concentration of the paramagnetic center at different locations about the solute of interest. The local concentrations may in turn be related to differences in the local free energies of interaction between the diffusing paramagnet and the cosolute. We demonstrate this approach for the case of ribonuclease A and deduce surface free energy differences for a large number of protein proton sites. We find that the oxygen accessibility is poorly represented by hard-sphere models such as computed solvent or steric accessibility. There is a distribution of local intermolecular interactions with a width of the order of RT that dominates the report of the intermolecular exploration of the protein by this simple solute.     

About long-lived nuclear spin states involved in para-hydrogenated molecules

This study deals with a spin system constituted of three nonequivalent protons, two of them originating from para-hydrogen (p-H(2)) after a hydrogenation reaction carried out in the earth magnetic field. It is shown that three singlet states are created provided indirect (J) couplings exist between the three spins, implying hyperpolarization transfer toward the third spin. Upon insertion of the sample in the NMR (Nuclear Magnetic Resonance) high field magnet, the following events occur: (i) the longitudinal two-spin orders which are parts of the singlet states survive; (ii) the other two terms (of these singlet states) tend to be destroyed by magnetic field gradients but at the same time are partly converted into differences of longitudinal polarizations. Nuclear spin relaxation is studied by appropriate NMR measurements when evolution takes place in the high field magnet or in the earth field. In the former case, relaxation is classical although complicated by numerous relaxation rates associated with both longitudinal two-spin orders and longitudinal polarizations. In the latter case, an equilibration between the singlet states first occur, their disappearance being thereafter driven by relaxation rates which remain very small because of the absence of any dipolar contribution. Thus, even in the case of a three-spin system, long-lived states exist; this unexpected property could be very useful for many applications.     

Behavior of Two Almost Identical Spins during the CPMG Pulse Sequence

Multiple‐spin‐echo experiments have found wide use in nuclear magnetic resonance spectroscopy. In particular, the Carr–Purcell–Meiboom–Gill (CPMG) pulse sequence is used to determine transverse relaxation times T₂. Herein it is demonstrated, both theoretically and experimentally, that for a pair of almost identical spins‐1/2 the experimental setup can have a profound effect on the observed spin dynamics. It is shown that, in the case of dipolar relaxation, the measured T₂ values can roughly vary between the limits of identical and unlike spins, just depending on the repetition rate of π pulses with respect to chemical shift separation. Such an effect can, in the extreme narrowing regime, amount to a 50 % difference.     

Nonideality in diffusion of ionic and hydrophobic solutes and pair dynamics in water-acetone mixtures of varying composition

We have performed a series of molecular dynamics simulations of water-acetone mixtures containing either an ionic solute or a neutral hydrophobic solute to study the extent of nonideality in the dynamics of these solutes with variation of composition of the mixtures. The diffusion coefficients of the charged solutes, both cationic and anionic, are found to change nonmonotonically with the composition of the mixtures showing strong nonideality of their dynamics. Also, the extent of nonideality in the diffusion of these charged solutes is found to be similar to the nonideality that is observed for the diffusion and orientational relaxation of water and acetone molecules in these mixtures which show a somewhat similar changes in the solvation characteristics of charged and dipolar solutes with changes of composition of water-acetone mixtures. The diffusion of the hydrophobic solute, however, shows a monotonic increase with increase of acetone concentration showing its different solvation characteristics as compared to the charged and dipolar solutes. The links between the nonideality in diffusion and solvation structures are further confirmed through calculations of the relevant solute-solvent and solvent-solvent radial distribution functions for both ionic and hydrophobic solutes. We have also calculated various pair dynamical properties such as the relaxation of water-water and acetone-water hydrogen bonds and residence dynamics of water molecules in water and acetone hydration shells. The lifetimes of both water-water and acetone-water hydrogen bonds and also the residence times of water molecules are found to increase steadily with increase in acetone concentration. No maximum or minimum was found in the composition dependence of these pair dynamical quantities. The lifetimes of water-water hydrogen bonds are always found to be longer than that of acetone-water hydrogen bonds in these mixtures. The residence times of water molecules are also found to follow a similar trend.     

Magnetic bistability of isolated giant-spin centers in a diamagnetic crystalline matrix

Polynuclear single-molecule magnets (SMMs) were diluted in a diamagnetic crystal lattice to afford arrays of independent and iso-oriented magnetic units. Crystalline solid solutions of an Fe(4) SMM and its Ga(4) analogue were prepared with no metal scrambling for Fe(4) molar fractions x down to 0.01. According to high-frequency EPR and magnetic measurements, the guest SMM species have the same total spin (S=5), anisotropy, and high-temperature spin dynamics found in the pure Fe(4) phase. However, suppression of intermolecular magnetic interactions affects magnetic relaxation at low temperature (40 mK), where quantum tunneling (QT) of the magnetization dominates. When a magnetic field is applied along the easy magnetic axis, both pure and diluted (x=0.01) phases display pronounced steps at evenly spaced field values in their hysteresis loops due to resonant QT. The pure Fe(4) phase exhibits additional steps which are firmly ascribed to two-molecule QT transitions. Studies on the field-dependent relaxation rate showed that the zero-field resonance sharpens by a factor of five and shifts from about 8 mT to exactly zero field on dilution, in agreement with the calculated variation of dipolar interactions. The tunneling efficiency also changes significantly as a function of Fe(4) concentration: the zero-field resonance is significantly enhanced on dilution, while tunneling at ±0.45 T becomes less efficient. These changes were rationalized on the basis of a dipolar shuffling mechanism and transverse dipolar fields, whose effect was analyzed by using a multispin model. Our findings directly prove the impact of intermolecular magnetic couplings on SMM behavior and disclose the magnetic response of truly isolated giant spins in a diamagnetic crystalline environment.     

Nuclear magnetic resonance proton dipolar order relaxation in thermotropic liquid crystals: a quantum theoretical approach

By means of the Jeener-Broekaert nuclear magnetic resonance pulse sequence, the proton spin system of a liquid crystal can be prepared in quasiequilibrium states of high dipolar order, which relax to thermal equilibrium with the molecular environment with a characteristic time (T1D). Previous studies of the Larmor frequency and temperature dependence of T1D in thermotropic liquid crystals, that included field cycling and conventional high-field experiments, showed that the slow hydrodynamic modes dominate the behavior of T1D, even at high Larmor frequencies. This noticeable predominance of the cooperative fluctuations (known as order fluctuations of the director, OFD) could not be explained by standard models based on the spin-lattice relaxation theory in the limit of high temperature (weak order). This fact points out the necessity of investigating the role of the quantum terms neglected in the usual high temperature theory of dipolar order relaxation. In this work, we present a generalization of the proton dipolar order relaxation theory for highly correlated systems, which considers all the spins belonging to correlated domains as an open quantum system interacting with quantum bath. As starting point, we deduce a formulation of the Markovian master equation of relaxation for the statistical spin operator, valid for all temperatures, which is suitable for introducing a dipolar spin temperature in the quantum regime, without further assumptions about the form of the spin-lattice Hamiltonian. In order to reflect the slow dynamics occurring in correlated systems, we lift the usual short-correlation-time assumption by including the average over the motion of the dipolar Hamiltonian together with the Zeeman Hamiltonian into the time evolution operator. In this way, we calculate the time dependence of the spin operators in the interaction picture in a closed form, valid for high magnetic fields, bringing into play the spin-spin interactions within the microscopic time scale. Then, by adopting the spin-temperature density operator to represent the collective state of the spin system, and removing the traditional hypothesis of high temperature, we deduce an expression for the first order quantum contribution to T1D (-1), in terms of spectral densities, with coefficients in form of spin traces. The properties that distinguish our result from the high-temperature T1D (-1) are as follows. (a) It is exclusively associated to cooperative fluctuations. (b) Because of its quantum character, it relies on both considering the lattice degrees of freedom quantum mechanically and including the spin-spin interactions in the microscopic time scale. With regard to the average dipolar Hamiltonian, only the nonsecular part plays a relevant role. (c) Associated with the structure of the spin operator involved in the quantum contribution, a term arises which is proportional to the number of spins in the correlated molecular domains, showing that the quantum contribution may be of macroscopic size in highly correlated systems. When applied to nematic liquid crystals, the new term exhibits the typical nu(-1/2) Larmor frequency dependence through the spectral density of the OFD, in consistence with the experimental results.     

Electron spin relaxation of N@C60 in CS2 in CS2

We examine the temperature dependence of the electron spin relaxation times of the molecules N@C60 and N@C70 (which comprise atomic nitrogen trapped within a carbon cage) in liquid CS2 solution. The results are inconsistent with the fluctuating zero-field splitting (ZFS) mechanism, which is commonly invoked to explain electron spin relaxation for S> or =1 spins in liquid solution, and is the mechanism postulated in the literature for these systems. Instead, we find an Arrhenius temperature dependence for N@C60 , indicating the spin relaxation is driven primarily by an Orbach process. For the asymmetric N@C70 molecule, which has a permanent ZFS, we resolve an additional relaxation mechanism caused by the rapid reorientation of its ZFS. We also report the longest coherence time (T2) ever observed for a molecular electron spin, being 0.25 ms at 170 K.     

Longitudinal nuclear spin relaxation of ortho- and para-hydrogen dissolved in organic solvents

The longitudinal relaxation time of ortho-hydrogen (the spin isomer directly observable by NMR) has been measured in various organic solvents as a function of temperature. Experimental data are perfectly interpreted by postulating two mechanisms, namely intramolecular dipolar interaction and spin-rotation, with activation energies specific to these two mechanisms and to the solvent in which hydrogen is dissolved. This permits a clear separation of the two contributions at any temperature. Contrary to the self-diffusion coefficients at a given temperature, the rotational correlation times extracted from the dipolar relaxation contribution do not exhibit any definite trend with respect to solvent viscosity. Likewise, the spin-rotation correlation time obeys Hubbard's relation only in the case of hydrogen dissolved in acetone-d6, yielding in that case a spin-rotation constant in agreement with literature data. Concerning para-hydrogen, which is NMR-silent, the only feasible approach is to dissolve para-enriched hydrogen in these solvents and to follow the back-conversion of the para-isomer into the ortho-isomer. Experimentally, this conversion has been observed to be exponential, with a time constant assumed to be the relaxation time of the singlet state (the spin state of the para-isomer). A theory, based on intermolecular dipolar interactions, has been worked out for explaining the very large values of these relaxation times which appear to be solvent-dependent.     

Solvent response and dielectric relaxation in supercooled butyronitrile

We have measured the dynamics of solvation of a triplet state probe, quinoxaline, in the glass-forming dipolar liquid butyronitrile near its glass transition temperature T(g)=95 K. The Stokes shift correlation function displays a relaxation time dispersion of considerable magnitude and the optical linewidth changes along the solvation coordinate in a nonmonotonic fashion. These features are characteristic of solvation in viscous solvents and clearly indicate heterogeneous dynamics, i.e., spatially distinct solvent response times. Using the dielectric relaxation data of viscous butyronitrile as input, a microscopic model of dipolar solvation captures the relaxation time, the relaxation dispersion, and the amplitude of the dynamical Stokes shift remarkably well.     

Long‐Lived States of Magnetically Equivalent Spins Populated by Dissolution‐DNP and Revealed by Enzymatic Reactions

Hyperpolarization by dissolution dynamic nuclear polarization (D ‐DNP) offers a way of enhancing NMR signals by up to five orders of magnitude in metabolites and other small molecules. Nevertheless, the lifetime of hyperpolarization is inexorably limited, as it decays toward thermal equilibrium with the nuclear spin‐lattice relaxation time. This lifetime can be extended by storing the hyperpolarization in the form of long‐lived states (LLS) that are immune to most dominant relaxation mechanisms. Levitt and co‐workers have shown how LLS can be prepared for a pair of inequivalent spins by D ‐DNP. Here, we demonstrate that this approach can also be applied to magnetically equivalent pairs of spins such as the two protons of fumarate, which can have very long LLS lifetimes. As in the case of para‐hydrogen, these hyperpolarized equivalent LLS (HELLS) are not magnetically active. However, a chemical reaction such as the enzymatic conversion of fumarate into malate can break the magnetic equivalence and reveal intense NMR signals.     

Nonadiabatic molecular dynamics simulations of correlated electrons in solution. 1. Full configuration interaction (CI) excited-state relaxation dynamics of hydrated dielectrons

The hydrated dielectron is composed of two excess electrons dissolved in liquid water that occupy a single cavity; in both its singlet and triplet spin states there is a significant exchange interaction so the two electrons cannot be considered to be independent. In this paper and the following paper,we present the results of mixed quantum/classical molecular dynamics simulations of the nonadiabatic relaxation dynamics of photoexcited hydrated dielectrons, where we use full configuration interaction (CI) to solve for the two-electron wave function at every simulation time step. To the best of our knowledge, this represents the first systematic treatment of excited-state solvation dynamics where the multiple-electron problem is solved exactly. The simulations show that the effects of exchange and correlation contribute significantly to the relaxation dynamics. For example, spin-singlet dielectrons relax to the ground state on a time scale similar to that of single electrons excited at the same energy, but spin-triplet dielectrons relax much faster. The difference in relaxation dynamics is caused by exchange and correlation: The Pauli exclusion principle imposes very different electronic structure when the electrons' spins are singlet paired than when they are triplet paired, altering the available nonadiabatic relaxation pathways. In addition, we monitor how electronic correlation changes dynamically during nonadiabatic relaxation and show that solvent dynamics cause electron correlation to evolve quite differently for singlet and triplet dielectrons. Despite such differences, our calculations show that both spin states are stable to excited-state dissociation, but that the excited-state stability has different origins for the two spin states. For singlet dielectrons, the stability depends on whether the solvent structure can rearrange to create a second cavity before the ground state is reached. For triplet dielectrons, in contrast, electronic correlation ensures that the two electrons do not dissociate, even if the dielectron is artificially kept from reaching the ground state. In addition, both singlet and triplet dielectrons change shape dramatically during relaxation, so that linear response fails to describe the solvation dynamics for either spin state. In the following paper (Larsen, R. E.; Schwartz, B. J. J. Phys. Chem. B 2006, 110, 9692), we use these simulations to calculate the pump-probe spectroscopic signal expected for photoexcited hydrated dielectrons and to predict an experiment to observe hydrated dielectrons directly.     

Probing side-chain dynamics in the proteasome by relaxation violated coherence transfer NMR spectroscopy

A pair of experiments is presented for measuring intra-methyl 1H-1H dipolar cross-correlated spin relaxation rates in highly deuterated, methyl protonated proteins with significantly improved sensitivity relative to previously developed experiments that measure dynamics via 1H spin relaxation. In applications to proteins with correlation times in the macromolecular limit, these cross-correlation rates are related directly to order parameters, characterizing the amplitude of motion of methyl-containing side-chains. The experimental approach is validated by comparing extracted order parameters with those obtained via 2H and 13C spin relaxation methods for both protein L (7.5 kDa) and malate synthase G (82 kDa), with excellent correlations obtained. The methodology is applied to study Ile, Leu, and Val side-chain dynamics in a 360 kDa "half-proteasome" complex. In particular, order parameters obtained from the WT complex and from a second complex where the proteasome gating residues are deleted establish that the relative levels of dynamics in each of the two molecules are very similar. It thus becomes clear that there is no communication between gating residues and other regions of the molecule involving pico- to nanosecond time-scale dynamics of these methyl-containing side-chains.     

Early Time Solvation Dynamics Probed by Spectrally Resolved Degenerate Pump-Probe Spectroscopy

The dynamic role of solvent in influencing the rates of physico-chemical processes (for example, polar solvation and electron transfer) has been extensively studied using time-resolved fluorescence spectroscopy. Here we study ultrafast excited state relaxation dynamics of three different fluorescent probes (DNTTCI, IR-140 and IR-144) in two polar solvents, ethanol and ethylene glycol, using spectrally resolved degenerate pump-probe spectroscopy. We discuss how time-resolved emission spectra can be directly used for constructing relaxation correlation function, obviating spectral reconstruction and estimation of time-zero spectrum in non-polar solvents. We show that depending on the specific probe used, the relaxation dynamics is governed either by intramolecular vibrational relaxation (for IR140) or by intermolecular solvation (for DNTTCI) or by both (for IR144). We further show (using DNTTCI as a probe) that major differences in solvation by ethanol and ethylene glycol is contributed by early time (<1 ps) dynamics.     

A new method to reconstruct three-dimensional spatial distribution function from radial distribution function in solvation structure

Three-dimensional spatial distribution function (SDF) of solvent is a fundamental quantity for analysis of solvation. However, its calculation has been very limited because long computational time is required. We here developed a novel and robust method to construct approximated SDFs of solvent sites from radial distribution functions. In this method, the expansion of SDFs in real solid harmonics around atoms of solute leads to a linear equation, from which SDFs are evaluated with reasonable computational time. This method is applied to the analysis of the solvation structure of liquid water, as an example. The successful results clearly show that this method is very powerful to investigate solvation structure.