the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
The solid effect of dynamic nuclear polarization in liquids – accounting for g-tensor anisotropy at high magnetic fields
Danhua Dai
Thomas F. Prisner
In spite of its name, the solid effect of dynamic nuclear polarization (DNP) is also operative in viscous liquids, where the dipolar interaction between the polarized nuclear spins and the polarizing electrons is not completely averaged out by molecular diffusion on the timescale of the electronic spin–spin relaxation time. Under such slow-motional conditions, it is likely that the tumbling of the polarizing agent is similarly too slow to efficiently average the anisotropies of its magnetic tensors on the timescale of the electronic T_{2}. Here we extend our previous analysis of the solid effect in liquids to account for the effect of g-tensor anisotropy at high magnetic fields. Building directly on the mathematical treatment of slow tumbling in electron spin resonance (Freed et al., 1971), we calculate solid-effect DNP enhancements in the presence of both translational diffusion of the liquid molecules and rotational diffusion of the polarizing agent. To illustrate the formalism, we analyze high-field (9.4 T) DNP enhancement profiles from nitroxide-labeled lipids in fluid lipid bilayers. By properly accounting for power broadening and motional broadening, we successfully decompose the measured DNP enhancements into their separate contributions from the solid and Overhauser effects.
- Article
(5782 KB) - Full-text XML
- BibTeX
- EndNote
The sensitivity of NMR experiments is greatly increased by dynamic nuclear polarization (DNP)^{1}, where the much larger static polarization that is available to electronic spins is transferred to nuclear spins (Atsarkin, 2011; Wenckebach, 2016). For the transfer to take place, the electronic and nuclear spins should be able to flip simultaneously (Abragam and Goldman, 1978). Such concerted flips correspond to the zero-quantum (ZQ) and double-quantum (DQ) transitions of the electron–nucleus spin system, which are enabled by inter-spin interactions. Among the four DNP mechanisms, i.e., the Overhauser effect (OE), the solid effect (SE), the cross effect and thermal mixing, only the first two have been conclusively shown to be operative in the liquid state, where the spin–spin interactions change randomly in time due to the thermal motions of the molecules.
In OE-DNP, the ZQ and DQ transitions are in fact possible because the dipole–dipole and contact interactions are modulated by molecular motions. In SE-DNP, on the other hand, the ZQ and DQ transitions are driven directly by mw excitation, and the modulation of the dipolar interaction is detrimental because it constantly modifies the matching condition that the mw frequency should satisfy in order to resonantly drive these transitions.
The initial theoretical treatments of OE (Solomon, 1955) and SE (Abragam and Proctor, 1958) modeled the ZQ and DQ transitions by expressing the transition probabilities per unit time using Fermi's golden rule. As the mathematical description of (semi-classical) relaxation theory matured around the same time (Redfield, 1957; Abragam, 1961), Fermi's golden rule was promptly replaced in the theory of OE-DNP in liquids (Hausser and Stehlik, 1968) with the correlation function of the dipolar interaction (or its Laplace transform, which is known as spectral density). Because the time-domain description of relaxation leads to a correlation function in a very general way (Abragam, 1961), the same formalism works naturally with different spectral densities (e.g., for rotational or translational diffusion). As an example, the improved analytical treatment of isotropic translational diffusion achieved in 1975 was immediately applied to paramagnetic relaxation in liquids (Ayant et al., 1975; Hwang and Freed, 1975).^{2}
During the same time period, it also became possible to account for spin dephasing and relaxation beyond second order (Anderson, 1954; Kubo, 1954), which is important for understanding spectral line shapes outside the regime of fast averaging (Kubo, 1969). These initial ideas were transformed into a powerful tool for the calculation and analysis of slow-motional EPR spectra by Freed (Freed et al., 1971; Freed, 1976).
When first presented, Abragam's quantitative description of SE-DNP in terms of mixing of the Zeeman energy levels by the dipolar interaction (Abragam and Proctor, 1958) conclusively explained that the NMR signal is maximally enhanced when the mw frequency is shifted from the electronic resonance by ±ω_{I}, where ω_{I} is the Larmor frequency of the polarized nuclear spin. Abragam's perturbative analysis also correctly predicted that the effect should drop quadratically with the magnitude of the static magnetic field, which has lasting implications for SE-DNP at high magnetic fields. In spite of these successes, however, the perturbative approach to SE is practically impossible to integrate with other relevant spin phenomena whose mathematical treatment matured subsequently.
Recently, Sezer (2023a) presented a time-domain description of SE which, like semi-classical relaxation theory, allows for different dynamical processes to modulate the relevant spin interactions. By interfacing this description with the spectral density of isotropic translational diffusion (Ayant et al., 1975; Hwang and Freed, 1975), it was possible to treat SE-DNP in the presence of molecular translation as relevant to homogeneous liquids (Sezer, 2023b). The requirement that the dipolar interaction should not be completely averaged out by the molecular dynamics during the electronic T_{2} restricts liquid-state SE-DNP to viscous media, where the tumbling of the polarizing agent may similarly be too slow to average the anisotropies of its magnetic tensors. The current paper accounts for the effect of g-tensor anisotropy on SE in this slow-tumbling regime. To this end, the time-domain description of SE-DNP in liquids is interfaced here with the established mathematical treatment of slow-motional EPR spectra (Freed et al., 1971). For the illustrative purposes of the current paper, we only consider free (i.e., unrestricted) rotational diffusion with an isotropic diffusion coefficient. Nevertheless, the treatment can be analogously extended to anisotropic diffusion in an orienting potential by building on the general mathematical formalism of the MOMD (microscopic order macroscopic disorder) and SRLS (slowly relaxing local structure) models (Meirovitch et al., 1984; Polimeno and Freed, 1995).
To motivate the presented theoretical analysis, in Sect. 2 we formulate one specific practical problem that it addresses. There we also introduce the experimental EPR and DNP data that are analyzed subsequently in Sect. 5 using the developed theory. The needed background from Sezer (2023a, b) is presented in Sect. 3. Building on this, in Sect. 4 we adapt the slow-motional formalism of Freed et al. (1971) to the treatment of SE in the liquid state. Our conclusions are in Sect. 6, and several supporting figures are left to the Appendix.
DNP aims to increase the longitudinal nuclear magnetization, i_{z}, beyond its equilibrium Boltzmann value, ${i}_{z}^{\mathrm{eq}}$. This is done by doping the sample with unpaired electrons, whose spins are then subjected to near-resonance mw irradiation. In cw-DNP, which is the only variety that we consider here, a steady-state magnetization ${i}_{z}^{\mathrm{ss}}$ is reached after the microwaves have been applied for a sufficiently long time. The enhancement of i_{z} under such steady-state conditions is
where ϵ=0 corresponds to the absence of DNP.
In both OE and SE, ϵ is directly proportional to the ratio of the gyromagnetic factors of the electronic and nuclear spins, γ_{S} and γ_{I}. For OE (Hausser and Stehlik, 1968; Müller-Warmuth and Meise-Gresch, 1983),
where s, c and f are, respectively, the electronic saturation factor, the coupling factor and the leakage factor. The first is defined as
and reflects the deviation of the longitudinal electronic magnetization at steady state, ${s}_{z}^{\mathrm{ss}}$, from its equilibrium value, ${s}_{z}^{\mathrm{eq}}$. The other two factors, c and f, quantify the interaction between the electronic and nuclear spins. Specifically, the leakage factor
compares the nuclear T_{1}s in the presence (T_{1I}) and absence (${T}_{\mathrm{1}I}^{\mathrm{0}}$) of the polarizing agent. In DNP, T_{1I} is typically (much) shorter than ${T}_{\mathrm{1}I}^{\mathrm{0}}$ due to the elevated concentration of the electronic spins, and hence f≈1.
Similarly, the SE enhancement can be expressed as (Sezer, 2023a)
where $p=\mathrm{1}-s$ quantifies how “non-saturated” the electronic transition is, and the rate constants v_{+} and v_{−} are related to the ability of the microwaves to excite simultaneous flips of the electronic and nuclear spins. These concerted flips correspond to the “forbidden” ZQ and DQ transitions, which are enabled by the dipolar interaction. In fact,
where v_{0} and v_{2} denote, respectively, the ZQ and DQ transition rate constants. In liquids, where the dipolar interaction is partially averaged, the contribution of the mw excitation to the nuclear relaxation rate ${R}_{\mathrm{1}I}=\mathrm{1}/{T}_{\mathrm{1}I}$, which is quantified by v_{+}, is generally negligibly small. As a result, ${v}_{+}/{R}_{\mathrm{1}I}\ll \mathrm{1}$, and the expression in parentheses in Eq. (5) is essentially 1. Then the SE enhancement acquires the following multiplicative form:
which is analogous to ϵ_{OE} with the factors s, c and f being replaced with the factors p, v_{−} and T_{1I}, respectively. In the numerical work presented in Sect. 5 we use the approximation in Eq. (7). The condition ${v}_{+}{T}_{\mathrm{1}I}\ll \mathrm{1}$ is validated at the end of the analysis by comparing the estimated v_{+} to the measured T_{1I}.
In the current paper we study the dependence of the DNP enhancement on the displacement from the electronic resonance. Following Gizatullin et al. (2022), we call the profile of ϵ against the offset from resonance a “DNP spectrum”. Because DNP experiments in the liquid state are carried out with a mw resonator (Erb et al., 1958a, b; Leblond et al., 1971b; Neudert et al., 2016; Gizatullin et al., 2021a; Kuzhelev et al., 2022, 2023), off-resonance conditions are achieved by varying the stationary magnetic field at a constant mw frequency (i.e., field sweep). In theoretical analysis, however, it is more convenient to work with a fixed B_{0} and a variable mw frequency. Thus, when comparing calculations and experiments, we will convert the horizontal axis of the experiments from units of magnetic field to units of offset frequency.
In the case of ϵ_{OE} (Eq. 2), the entire offset dependence is due to the saturation factor s, as the factors c and f are practically constant over such a narrow frequency range. In the case of ϵ_{SE} (Eq. 5), both pv_{−} and v_{+} are functions of the offset. For a single, homogeneously broadened EPR line the saturation factor can be obtained in closed analytical form from the Bloch equations (as we review below in Sect. 3.1). Recently, Sezer (2023a) showed that the SE spin dynamics is described by two coupled Bloch equations, whose steady state can similarly be solved analytically to obtain closed-form expressions for the rate constants v_{±} (reviewed in Sect. 3.2). In liquids, where the random molecular motion modulates the dipolar interaction between the electronic and nuclear spins, these rate constants are no longer available analytically but can be calculated numerically for motional models with known dipolar spectral densities (Sezer, 2023b), as reviewed below in Sect. 3.3.
Liquid-state SE-DNP is restricted to viscous media, where the dipolar interaction is not averaged out completely on the decoherence timescale of the electronic spins. Under these conditions, the tumbling of the polarizing agent is also expected to be too slow to average the anisotropies of its magnetic tensors on the timescale of the electronic T_{2}. One thus expects substantial deviations from the Lorentzian EPR line shape of the Bloch equations. Such deviations are unavoidable in the case of nitroxide-based polarizing agents whose g and A tensors are rather anisotropic. A recent SE-DNP study at 9.4 T demonstrated that even the narrow-line radical trityl exhibited g-tensor broadening in liquid glycerol (Kuzhelev et al., 2023).
This paper extends the theoretical description of SE-DNP to the regime of slow radical tumbling, where the cw-EPR line shape is not Lorentzian. Given our long-standing efforts in liquid-state DNP at 9.4 T, here we focus on high magnetic fields, where the width of the EPR spectrum is dominated by the anisotropy of the g tensor. We will thus completely neglect the hyperfine tensor. This possibility greatly simplifies the needed adjustments to the Lorentzian case (Sect. 4).
To illustrate the practical problem that motivated this theoretical work, we now turn to the experimental data in Fig. 1. The characterized samples comprised liposomes of hydrated lipid bilayers composed of DOPC (1,2-dioleoyl-sn-glycero-3-phosphocholine) lipids. As the phase transition temperature of DOPC is about −17 ^{∘}C, the lipids were in their fluid, liquid-crystalline phase in the experiments at ≈320 K. The DOPC lipids were mixed at a ratio of 20:1 with PSPC lipids spin-labeled either at position 10 (1-palmitoyl-2-stearoyl-(10-doxyl)-sn-glycero-3-phosphocholine) or at position 16 along one of their aliphatic chains. Both the EPR spectra (Fig. 1a, b) and the DNP enhancements (Fig. 1c, d) were recorded in our home-built Fabry–Pérot resonator at 9.4 T equipped with a temperature control (Denysenkov et al., 2022). While the target temperature of the experiments was 320 K, an extra temperature rise of less than 10 ^{∘}C can be expected at the maximum mw power of 5.5 W that was used for DNP (Denysenkov et al., 2022). Details about the experiments and the sample preparation will be published elsewhere.
The cw-EPR spectrum of 10-Doxyl-PC in Fig. 1a (orange line) is seen to deviate substantially from (the derivative of) a Lorentzian line shape. At this high magnetic field, the EPR line width is expected to be dominated by the large anisotropy of the nitroxide g tensor, with a comparatively much smaller contribution from the nitroxide hyperfine tensor. (These expectations are tested and verified below in Sect. 5.1.) For comparison, in Fig. 1b we show the cw-EPR spectrum of the sample doped with 16-Doxyl-PC. Visually, this narrower spectrum more closely resembles a homogeneous Lorentzian line, although it still deviates from it (as discussed in Sect. 5.2).
In Fig. 1c we show the DNP spectrum (filled red circles) of the sample containing 10-Doxyl-PC as a polarizing agent. The enhanced NMR signal belongs to the acyl chain protons of the lipids. Thanks to the high magnetic field of the experiment, it was possible to resolve the NMR signal of these non-polar protons from the polar protons of water and of the lipid head groups. The DNP spectrum is seen to have a complex line shape, with positive enhancement values at offsets of about +400 MHz demonstrating a contribution from SE. At the same time, the comparatively larger negative enhancements in the vicinity of the electronic resonance (i.e., around 0 MHz) point to a contribution from OE. Such coexistence of SE and OE is well documented for nitroxide free radicals at the classical EPR fields of about 0.35 T (Leblond et al., 1971b; Neudert et al., 2017; Gizatullin et al., 2021a, b). Evidently, it also persists at 9.4 T. The DNP spectrum of 16-Doxyl-PC in Fig. 1d also exhibits a mixture of SE and OE.
More than half a century ago, Korringa and coworkers developed a rigorous theoretical framework to predict such mixed DNP spectra in viscous liquids (Papon et al., 1968; Leblond et al., 1971a). Likely because of its complexity and its neglect of translational diffusion, their formal analysis has not been applied to recent DNP data. As a simple and practical alternative, Neudert et al. (2017) disentangled the OE and SE components of such mixed DNP spectra using only the integral of the measured cw-EPR signal. Their approach is based on the following insightful observations: (i) up to an overall scaling factor, the EPR line shape is equal to the saturation factor and thus to the OE enhancement (Eq. 2); (ii) up to an overall scaling factor, the SE enhancement lines at ±ω_{I} are shifted versions (and flipped for the ZQ transition) of the same EPR line shape. One can thus identify the contributions of OE and SE to the DNP spectrum by placing the integrated cw-EPR spectrum at, respectively, zero and ±ω_{I} offsets and independently adjusting the magnitudes of the two components.
This approach is illustrated in Fig. 1c and 1d, where the dashed–dotted blue lines are the integrals of the cw-EPR spectra from Fig. 1a and 1b, respectively (flipped here to reflect the dipolar nature of OE), and the dashed green lines are the same EPR spectra but centered at −400 and +400 MHz. The sum of the OE and SE contributions determined in this way is shown with a dotted black line. This sum is seen to agree closely with the DNP spectrum of 10-Doxyl-PC (Fig. 1c) and to capture well the overall shape of the DNP spectrum of 16-Doxyl-PC (Fig. 1d).
In spite of the good general agreement between the experimental DNP spectra and the dotted black lines in Fig. 1c and 1d, some persistent differences remain. In particular, (i) the OE feature in the experiment appears to be consistently broader than the EPR line and (ii) the enhancement between the central OE feature and the negative SE feature is consistently larger than what is predicted by the overlap of the two copies of the EPR line shape. Both of these aspects are especially clear in the case of 16-Doxyl-PC (Fig. 1d). The theory presented below (Sect. 4) aims to address these deficiencies of the simple approach.
In fact, the first deficiency is easy to rationalize. Cw-EPR spectra are recorded at low mw power, and their widths reflect mechanisms contributing to the electronic T_{2} relaxation. The DNP spectrum, on the other hand, is recorded at high mw power, where the EPR line width experiences power broadening that also depends on the electronic T_{1} relaxation. That the OE-DNP spectrum “represents an indirect observation of the electron resonance when greatly saturated” was understood early on (Carver and Slichter, 1956, Fig. 6). To properly model the contribution of OE to mixed DNP spectra, therefore, it is necessary to calculate the cw-EPR spectrum under saturating conditions. How to rigorously do that in the regime of slow radical tumbling is known (Freed et al., 1971).
While power broadening affects OE, it is not immediately clear whether one should also take it into account when modeling SE. (We address this point in Sect. 4.4.) Even leaving power broadening aside, however, we know that in liquids the SE lines of the DNP spectrum should also be broader than the EPR line width because of the fluctuations of the dipolar interaction (Sezer, 2023b). Although Sezer (2023b) showed how to quantify this additional motional broadening in the case of translational molecular diffusion, the theoretical treatment there assumed a Lorentzian EPR line and is thus not directly applicable to the experiments in Fig. 1. In the current paper, we extend the formalism to slow radical tumbling and g-tensor anisotropy (Sect. 4). In Sect. 5 we apply the developed theory to the analysis of the experimental spectra in Fig. 1, disentangling the contributions of SE and OE to the observed DNP. The needed theoretical background from Sezer (2023a, b) is reviewed next.
The classical Bloch equations describe the dynamics of the electronic magnetization, including under saturating conditions. In Sect. 3.1 we recall the relationship between the steady-state solution of the Bloch equations and cw-EPR. Then, in Sect. 3.2, closed-form expressions are obtained for the rate constants of the forbidden transitions that are driven by the microwaves in SE-DNP. These expressions, derived in this form for the first time (Eqs. 33 and 34), are similar to the steady-state solutions of the Bloch equations but additionally contain (i) the strength of the electron–nucleus dipolar interaction and (ii) the Larmor frequency of the polarized nuclear spin (Sezer, 2023a). Finally, in Sect. 3.3 we remind the reader how these expressions should be modified in the presence of random modulation of the dipolar interaction, as relevant for liquids (Sezer, 2023b).
The reviewed results, which apply to a single Lorentzian line, will be extended in Sect. 4 to the regime of slow radical tumbling and an anisotropic g tensor. In the process, some of the scalar variables that appear below, like the offset frequency and the electronic relaxation rates, will be replaced with square matrices, as we explain in Sect. 4.1 and 4.2. The generalization of Sect. 3.1, 3.2 and 3.3 along these lines is carried out in, respectively, Sect. 4.3, 4.4 and 4.5.
3.1 Bloch equations
The evolution of the expectation values of the electronic spin operators S_{i} ($i=x,y,z$), which we denote by s_{i}, is described by the classical Bloch equations (in the rotating frame)
Here, the dot above the variable indicates differentiation with respect to time, R_{2} and R_{1} are the reciprocals of the electronic relaxation times T_{2} and T_{1}, respectively, and^{3}
is the offset between the Larmor frequency of the electronic spins, ω_{0}, and the (angular) frequency of the oscillating magnetic field, ω. In the case of an isotropic g factor, g_{0},
where μ_{B} is the Bohr magneton and ℏ is the reduced Planck constant.
At steady state,
Solving these algebraic equations for the variables ${s}_{i}^{\mathrm{ss}}$, one can calculate the cw-EPR spectrum and the electronic saturation profile. Making use of the zeros in the first and last rows of the Bloch matrix in Eq. (11), we first express ${s}_{x}^{\mathrm{ss}}$ and ${s}_{z}^{\mathrm{ss}}$ in terms of ${s}_{y}^{\mathrm{ss}}$:
The middle row of the matrix then yields
where we defined
The in-phase (absorptive) and out-of-phase (dispersive) components of the cw-EPR signal are then found to be
From the longitudinal component at steady state, we similarly find
which shows that the saturation factor is directly proportional to the absorptive EPR line shape. This proportionality holds for all mw powers, including the large powers used in DNP. In Sect. 4.3 we show that this also remains valid in the case of g-tensor anisotropy and isotropic rotational diffusion.
When generalizing the Bloch equations to an anisotropic g tensor, we will need to work with high-dimensional abstract vectors. To distinguish these vectors from the vectors in 3D space, we will denote the latter by placing an arrow above their symbols and will use bold symbols for the former. (A 3D unit vector will be indicated with a hat rather than an arrow.) Additionally, we will use capital hollow letters to denote 3×3 matrices that act on the 3D vectors. With this understanding, we will write the Bloch equations in Eq. (8) as
where
and ${\mathbb{B}}_{\mathrm{0}}=\mathbb{B}({\mathit{\omega}}_{I}=\mathrm{0})$. The iω_{I} that has been added to the main diagonal of the matrix 𝔹 will be needed for the dynamical description of the solid effect (see Sect. 3.2). The subscript of 𝔹_{0} is intended as a reminder that 𝔹 is evaluated at ω_{I}=0, where ω_{I} is the Larmor frequency of the polarized nuclear spin.
3.2 Solid effect in solids
SE relies on the dipolar interaction between the electronic and nuclear spins whose coupling is
Here ${D}_{\mathrm{dip}}=({\mathit{\mu}}_{\mathrm{0}}/\mathrm{4}\mathit{\pi})\mathrm{\hslash}{\mathit{\gamma}}_{S}{\mathit{\gamma}}_{I}$ is the dipolar constant, which equals approximately 2π (79 kHz nm^{3}) for protons, and $(r,\mathit{\theta},\mathit{\varphi})$ are the spherical polar coordinates of the inter-spin vector.
In liquids, A_{1} changes in time because of molecular diffusion. The treatment of SE-DNP for a time-dependent A_{1} in Sezer (2023b) was developed under the assumption that the nuclear T_{1} is orders of magnitude larger than the correlation time of the electron–nucleus dipolar interaction, which is practically always the case in liquids. For the same analysis to apply to solids, nuclear spin diffusion, which analogously to molecular diffusion in liquids spreads out the nuclear polarization across the sample, should be much faster than the nuclear T_{1}. Although this condition is not necessarily satisfied in the solid state, for the mathematical description in terms of a dipolar correlation function to apply, we will assume that spin diffusion is fast when referring to solids. Similarly, when accounting for g-tensor anisotropy below, we will assume that the tumbling of the radical is much faster than the nuclear T_{1}. This assumption is clearly violated in solids, where “tumbling” is infinitely slow. Nevertheless, for the purposes of comparison, we will refer in the following to “solids” with the understanding that the correlation time of the dipolar interaction is much shorter than the nuclear T_{1} (in order to treat nuclear spin diffusion on the level of a translational correlation function) but much longer than all other relaxation timescales (in order to treat the electron–nucleus dipolar interaction as constant). Because we will keep all other parameters, including the timescale of radical tumbling, the same when comparing “solids” and liquids, it should be kept in mind that our treatment is not a good model for the solid state (hence the quotation marks).
For SE-DNP, in addition to the Bloch equations, it is necessary to consider the following dynamical equations of the electron–nucleus coherences ${g}_{i}=\langle {S}_{i}{I}_{+}\rangle $ ($i=x,y,z$) (Sezer, 2023a):
Again, we are only interested in the steady state of the dynamics where
The rate constants pv_{−} and v_{+} needed to calculate the SE enhancement (Eq. 5) are determined from ${g}_{z}^{\mathrm{ss}}$ using the following equality, which combines Sezer (2023a, Eq. 31) and Sezer (2023b, Eq. 42):
(Re{} takes the real part of its argument.) The term proportional to ${R}_{\mathrm{1}I}^{A}$ on the right-hand side of Eq. (23) accounts for the contribution of the coherences g_{i} to the nuclear T_{1} relaxation in the absence of mw excitation. This contribution should be removed when calculating the mw-related rates v_{+} and pv_{−}.
To a good approximation, the electronic spin dynamics is independent of the dipolar interaction with the nuclear spins, as other mechanisms are more efficient at causing electronic relaxation, especially in liquids. As a result, the steady-state expressions from Sect. 3.1 can be used when solving Eq. (22) for ${g}_{z}^{\mathrm{ss}}$.
Inverting the matrix 𝔹 in Eq. (22) and using ${s}_{x,y}^{\mathrm{ss}}$ from before, we find
where [𝔹^{−1}]_{ij} is the ijth matrix element of 𝔹^{−1}. Substituting this ${g}_{z}^{\mathrm{ss}}$ into Eq. (23), we identify the desired SE rate constants
where
reflects the strength of the dipolar interaction. (Im{} takes the imaginary part of its argument.)
In liquids, where A_{1} is time-dependent, we will need to modify the matrix 𝔹^{−1} in Eq. (25) without changing the structure of these expressions (Sect. 3.3). In the case of solids (i.e., when A_{1} does not change with time), it is possible to carry out the inversion of 𝔹 by expressing ${g}_{x}^{\mathrm{ss}}$ and ${g}_{z}^{\mathrm{ss}}$ in terms of ${g}_{y}^{\mathrm{ss}}$, analogously to our treatment of the Bloch equations in the previous subsection.
From the upper and lower rows of 𝔹 in Eq. (22), we find
Substituting this ${g}_{z}^{\mathrm{ss}}$ into Eq. (23), we obtain
The first term on the right-hand side of Eq. (28) vanishes when ω_{1}=0. In contrast, the term in the second line is independent of ω_{1} and thus also contributes in the absence of mw excitation. We thus identify this second term with the thermal relaxation rate
Since we are not interested in this rate, the second summand in Eq. (28) can be dropped at this stage. The rate constants v_{+} and pv_{−} will thus be identified using only the first line in Eq. (28):
Substituting ${g}_{x}^{\mathrm{ss}}$ and ${g}_{z}^{\mathrm{ss}}$ from Eq. (27) into the middle equality of Eq. (22) and using the electronic steady state, we find
where
generalizes Eq. (14) such that ${P}_{\mathrm{0}}=P({\mathit{\omega}}_{I}=\mathrm{0})$. Finally, using this ${g}_{y}^{\mathrm{ss}}$ in Eq. (30), we obtain
and
In these expressions we have written down the combinations P and P_{0} explicitly in order to show in closed form how v_{+} and pv_{−} depend on all the parameters. For example, we immediately see that pv_{−} is odd in the offset Δ, while v_{+} is even. Because the SE-DNP enhancement is proportional to the ratio of these two rates (Eq. 5), it has the characteristic odd (i.e., antisymmetric) dependence on the offset from the electronic resonance.
When generalizing the SE spin dynamics to g-tensor anisotropy, we will write the dynamical equations (Eq. 21) as
with
3.3 Solid effect in liquids
The modulation of the dipolar interaction by translational diffusion was described in Sezer (2023b) on the level of the spectral density of the motional model, which was denoted by J_{11}(s) since this is the Laplace transform of the autocorrelation function of the dipolar interaction A_{1} (hence the double subscript of J). As an example, the spectral density of the FFHS model of translational diffusion is (Ayant et al., 1975; Hwang and Freed, 1975)
Here, the parameter
is the diffusive timescale of the model, which depends on the contact distance of the electronic and nuclear spins, b, and on the coefficient of their relative translational diffusion, D_{trans}, and
is the average of the dipolar interaction strength δ^{2} over the sample volume times the concentration of the electronic spins, N.
It is convenient to write J_{11}, which has units of angular frequency, as
where j_{11}(s) has units of time. This factorization confines the effect of the parameters N and b and the constant D_{dip} to the scaling factor 〈δ^{2}〉. The factor j_{11}(s) then fully accounts for the line shape of the SE-DNP spectrum, which results from the interplay between the offset frequency and the timescale of the translational motion.
According to Sezer (2023b), the modification from solids to liquids amounts to replacing the matrix 𝔹^{−1} in Eq. (25) with the matrix
and also replacing δ^{2} with 〈δ^{2}〉. The desired SE rate constants in liquids are thus
We now clarify the meaning of Eq. (41). Following the definition of the function of a matrix, one should first solve the eigenvalue problem of 𝔹, i.e., 𝔹𝕌=𝕌Λ, where the diagonal matrix $\mathrm{\Lambda}=\text{diag}({\mathit{\lambda}}_{\mathrm{1}},{\mathit{\lambda}}_{\mathrm{2}},{\mathit{\lambda}}_{\mathrm{3}})$ contains the three eigenvalues and the columns of 𝕌 contain the corresponding (right) eigenvectors. Then one should evaluate the spectral density at the three eigenvalues: ℓ_{n}=j_{11}(λ_{n}). Finally, one should form the diagonal matrix $\mathbb{L}=\text{diag}({\mathrm{\ell}}_{\mathrm{1}},{\mathrm{\ell}}_{\mathrm{2}},{\mathrm{\ell}}_{\mathrm{3}})$ and calculate $\mathbb{Q}=\mathbb{U}\mathbb{L}{\mathbb{U}}^{-\mathrm{1}}$. Comparing this expression of ℚ with ${\mathbb{B}}^{-\mathrm{1}}=\mathbb{U}\mathrm{\Lambda}{\mathbb{U}}^{-\mathrm{1}}$, we see that, in the transition from solids to liquids, where 𝔹^{−1} is replaced with ℚ, we essentially “process” the eigenvalues of 𝔹 with the spectral density function j_{11}. This step prevents us from eliminating the variables ${g}_{x,z}^{\mathrm{ss}}$ in the way we did previously for solids (Sect. 3.2). Because of that, the rate constants in liquids (Eq. 42) need to be calculated numerically.
Nonetheless, it is still possible to simplify the expression for ${R}_{\mathrm{1}I}^{A}$, since when ω_{1}=0, the zz component of 𝔹 is decoupled from the rest of the matrix. One then finds
Clearly, the time dependence of the dipolar interaction modifies all rate constants, including ${R}_{\mathrm{1}I}^{A}$ (cf. Eq. 29).
In this section we show how to account for g-tensor anisotropies when the tumbling of the radical is slow. Because our description of SE is built around the Bloch equations (Sezer, 2023a, b), we first adapt the treatment of isotropic rotational diffusion of Freed et al. (1971) to our needs (Sect. 4.1, 4.2 and 4.3) and then generalize it to SE-DNP (Sect. 4.4 and 4.5). If needed, further generalization to anisotropic diffusion and an orienting potential can be carried out analogously, following the mathematical treatment of the MOMD and SRLS models for slow-motional EPR (Meirovitch et al., 1984; Polimeno and Freed, 1995).
4.1 Stochastic Liouville equation for isotropic tumbling
Following Freed et al. (1971), we account for the effect of tumbling on the EPR spectrum using the SLE formalism (Anderson, 1954; Kubo, 1954). We describe the rotational state of the radical statistically with the probability density P(Ω,t), which quantifies the likelihood that at time t the molecular system of coordinates in which the g tensor is diagonal will have orientation Ω with respect to the laboratory system of axes defined by the magnetic fields B_{0} and B_{1}. In the case of isotropic rotation, this probability evolves with the Fokker–Planck equation
where D_{rot} is the rotational diffusion constant of the radical and the Laplace differential operator ${\mathrm{\nabla}}_{\mathrm{\Omega}}^{\mathrm{2}}$ acts on the orientation variable Ω. The operator
satisfies the following eigenvalue problem:
where the eigenfunctions ${\mathcal{D}}_{mn}^{\mathrm{\ell}}\left(\mathrm{\Omega}\right)$ are the Wigner rotation matrix elements, which are orthogonal to each other:
From Eq. (46) it is clear that the time derivative on the left-hand side of Eq. (44) vanishes for the equilibrium probability
In the presence of g-tensor anisotropy, the electronic Larmor frequency depends on the orientation Ω of the radical as follows:^{4}
with the angular frequencies (Freed et al., 1971)
These are formed from the components g_{xx}, g_{yy} and g_{zz} of the g tensor in the molecular frame. In Eq. (49), the first index in the subscripts of the Wigner rotation matrix elements refers to the molecular system of axes, while the second index refers to the laboratory system. The second indices are zero here because we only consider the secular terms, which are proportional to S_{z}.
Since the electronic Larmor frequency depends on Ω, the offset frequency Δ also becomes a function of the molecular orientation. As an example, for a fixed Ω, the Bloch equations (Eq. 17) should be modified as
where the orientation dependence is confined to the 3×3 matrix
(The matrix 𝔾 was introduced in Eq. 36.) It should be stressed, however, that Eq. (51) is not a legitimate equation of motion, as it does not account for the dynamics of the orientation Ω.
The SLE formalism remedies this deficiency by introducing the orientation-conditioned averages $\overrightarrow{s}(\mathrm{\Omega},t)$, whose spatial part evolves according to the Bloch equations (Eq. 51) and whose Ω dependence evolves according to the diffusion equation (Eq. 44):
Here 𝔼 is the 3×3 identity matrix in 3D space and E_{Ω} is the identity operator in the same abstract space as K_{Ω}. The outer product ⊗ is needed to create a combined operator that acts simultaneously in both of these spaces.
Since the functions ${\mathcal{D}}_{mn}^{\mathrm{\ell}}\left(\mathrm{\Omega}\right)$ form a complete set, we expand $\overrightarrow{s}(\mathrm{\Omega},t)$ as follows:
The coefficients ${\overrightarrow{s}}_{mn}^{\mathrm{\ell}}$, which contain the time dependence, can be obtained from $\overrightarrow{s}(\mathrm{\Omega},t)$ using the orthogonality of ${\mathcal{D}}_{mn}^{\mathrm{\ell}}\left(\mathrm{\Omega}\right)$ (Eq. 47):
Ultimately, the only property that we care about is the integral of the SLE variable $\overrightarrow{s}(\mathrm{\Omega},t)$ over all the orientations:
In that sense, the (vector) coefficient ${\overrightarrow{s}}_{\mathrm{00}}^{\mathrm{0}}\left(t\right)$ is the main object of interest, while all the other coefficients ${\overrightarrow{s}}_{mn}^{\mathrm{\ell}}\left(t\right)$ play an auxiliary, bookkeeping role.
Substituting $\overrightarrow{s}(\mathrm{\Omega},t)$ from Eq. (54) into Eq. (53), multiplying both sides by ${\mathcal{D}}_{MN}^{L*}\left(\mathrm{\Omega}\right)$ and integrating over Ω, we get
Clearly, the terms proportional to K_{Ω} and 𝔹_{0} in Eq. (53) do not mix coefficients ${\overrightarrow{s}}_{MN}^{L}$ with different values of L, M and N. In other words, these two operators are diagonal in the selected representation. The term proportional to 𝔽(Ω), on the other hand, mixes coefficients with different L and M (but not N, as we discuss below).
The integral in the last line of Eq. (57) contains the product of three Wigner rotation matrix elements. These can be expressed in terms of the Clebsch–Gordan coefficients ${C}_{{\mathrm{\ell}}_{\mathrm{1}}{m}_{\mathrm{1}}{\mathrm{\ell}}_{\mathrm{2}}{m}_{\mathrm{2}}}^{LM}$. Specifically, for the ${\mathcal{D}}_{K\mathrm{0}}^{\mathrm{2}}\left(\mathrm{\Omega}\right)$ in Eq. (52), we have
which leads to
In Eq. (59), the sum over ℓ mixes only expansion coefficients with $\mathrm{\ell}=L,L\pm \mathrm{2}$ (Freed et al., 1971) because all three Wigner rotation matrix elements in 𝔽 have L=2 (Eq. 52). Since we need ${\overrightarrow{s}}_{\mathrm{00}}^{\mathrm{0}}$ at the end, it is sufficient to consider only coefficients with even values of ℓ. Furthermore, as the Wigner rotation matrix elements in 𝔽 have $M=\mathrm{0}\pm \mathrm{2}$ and N=0, the sum over m mixes only coefficients whose values m are either equal to M or differ from it by two units, while the sum over n does not mix any coefficients with n different from N. These considerations imply that the triple sum in Eq. (59) will only go over ${\overrightarrow{s}}_{m\mathrm{0}}^{\mathrm{\ell}}$ with even ℓ and m. Finally, because the Wigner rotation matrix elements with M=2 and $M=-\mathrm{2}$ appear in a symmetrical way in 𝔽, it becomes possible to work with the symmetrized coefficients (Freed et al., 1971)
thus restricting M to non-negative values ($\mathrm{0}\le M\le L$). The lowest-order coefficients that are coupled by the SLE dynamics are thus ${\overrightarrow{s}}^{\mathrm{00}}$, ${\overrightarrow{s}}^{\mathrm{20}}$, ${\overrightarrow{s}}^{\mathrm{22}}$, ${\overrightarrow{s}}^{\mathrm{40}}$, ${\overrightarrow{s}}^{\mathrm{42}}$, ${\overrightarrow{s}}^{\mathrm{44}}$, ${\overrightarrow{s}}^{\mathrm{60}}$, etc.
4.2 Matrix representation of the SLE dynamics
While the above considerations greatly reduce the needed coefficients, there are still an infinite number left. In any practical work, this infinite set is truncated by selecting a maximum value of L to account for and setting to zero the coefficients with L>L_{max}. Since the total number of even L such that L≤L_{max} is ${n}_{L}={L}_{\mathrm{max}}/\mathrm{2}+\mathrm{1}$, the total number of remaining coefficients ${\overrightarrow{s}}^{LM}$ is ${n}_{\mathrm{tot}}={n}_{L}({n}_{L}+\mathrm{1})/\mathrm{2}={L}_{\mathrm{max}}^{\mathrm{2}}/\mathrm{8}+\mathrm{3}{L}_{\mathrm{max}}/\mathrm{4}+\mathrm{1}$. For the smallest non-trivial choice of L_{max}=2, n_{tot}=3 (with ${\overrightarrow{s}}^{\mathrm{00}}$, ${\overrightarrow{s}}^{\mathrm{20}}$ and ${\overrightarrow{s}}^{\mathrm{22}}$). The number of coefficients increases quadratically with L_{max} (e.g., ${n}_{\mathrm{tot}}=\mathrm{15},\mathrm{28},\mathrm{45}$ for ${L}_{\mathrm{max}}=\mathrm{8},\mathrm{12},\mathrm{16}$, respectively).
To compactly write down how these coefficients are mixed by the SLE dynamics, we introduce the following abstract vectors with n_{tot} elements:
where the former is needed for the first term on the right-hand side of Eq. (59). The SLE dynamics then becomes
where
is a 3n_{tot}×3n_{tot} matrix, and E, R_{1}, R_{2} and Δ are n_{tot}×n_{tot} matrices.
The first three of these sub-matrices are purely diagonal: E is the identity matrix and
with the diagonal elements of C_{D} being equal to L(L+1). For the simplest case of L_{max}=2 with only three coefficients (${\overrightarrow{s}}^{\mathrm{00}}$, ${\overrightarrow{s}}^{\mathrm{20}}$ and ${\overrightarrow{s}}^{\mathrm{22}}$),
In Eq. (63), the diagonal matrices R_{1,2} and E, which originate from the second line of Eq. (59), do not mix coefficients with different L and M. Only the sub-matrix Δ, which is of the form
mixes coefficients of different orders. In fact, the mixing is due to the matrices C_{0,2}, which modify the frequency offset Δ in proportion to the g-tensor anisotropies ${\mathit{\gamma}}_{\mathrm{0}}^{\mathrm{2}}$ and ${\mathit{\gamma}}_{\mathrm{2}}^{\mathrm{2}}$. For L_{max}=2,
(The factors of 2 in the last column of C_{2} arise from the fact that coefficients with M=0 pose an exception to the symmetrization in Eq. 60.) The matrix elements of these two matrices in the most general case are
where the summand proportional to δ_{M0} in the second line accounts for the factor of 2 that is needed by the coefficients ${\overrightarrow{s}}^{L\mathrm{0}}$.
Selecting L_{max}=0 in the above formalism amounts to retaining only the (3D vector) coefficient ${\overrightarrow{s}}^{\mathrm{00}}$. Then the matrix ℬ_{0} in Eq. (63) reduces to 𝔹_{0}, and Eq. (62) reduces to the classical Bloch equations for a homogeneous line. For L_{max}>0, the diagonal matrices R_{1} and R_{2} cause the coefficients ${s}_{z}^{LM}$ and ${s}_{x,y}^{LM}$, respectively, to decay exponentially, with those with larger L being suppressed more strongly by the tumbling. Analogously to the Bloch equations, the mw excitation mixes the y and z components of ${\overrightarrow{s}}^{LM}$ without mixing their LM dependence. The latter is mixed only by the offset matrix Δ, as elaborated above.
By building the SLE dynamics on top of the classical Bloch equations, we have arrived at a rather intuitive picture of how the g-tensor anisotropy is incorporated into the spin dynamics. Specifically, every element of the Bloch matrix 𝔹_{0} (Eq. 19 with ω_{I}=0) is replaced with a matrix in the space of LM indices (Eq. 63). In this replacement, all the elements except the frequency offset become diagonal matrices in the LM space, with the mixing in this space being entirely due to the offset. Since we describe the solid effect by two coupled Bloch equations, this intuition about the effect of g-tensor anisotropy on the spin dynamics will be helpful when adapting the approach to SE-DNP in Sect. 4.4 and 4.5.
4.3 EPR spectrum and saturation
The cw-EPR spectrum and the electronic saturation factor under g-tensor anisotropy are obtained from the steady state of Eq. (62),
which can be solved by inverting the 3n_{tot}×3n_{tot} matrix ℬ_{0} numerically. However, it is also possible to solve Eq. (69) by inverting a single matrix with dimensions that are 3 times smaller (i.e., n_{tot}×n_{tot}), as we show next.
First, taking advantage of the zeros in ℬ_{0} (Eq. 63), we express ${\mathit{s}}_{x}^{\mathrm{ss}}$ and ${\mathit{s}}_{z}^{\mathrm{ss}}$ in terms of ${\mathit{s}}_{y}^{\mathrm{ss}}$:
Because only the first element of 1^{00} is non-zero and the diagonal matrix R_{1} does not mix coefficients with different values of LM, the second equality in Eq. (70) becomes
For the 00th (i.e., first) element of ${\mathit{s}}_{z}^{\mathrm{ss}}$, we thus have ${s}_{z}^{\mathrm{00}}={\mathit{\omega}}_{\mathrm{1}}{T}_{\mathrm{1}}{s}_{y}^{\mathrm{00}}+{s}_{z}^{\mathrm{eq}}$, which is identical to the second equality in Eq. (12). We thus conclude that the proportionality between the electronic saturation factor and the in-phase EPR line shape (Eq. 16) is not limited to a homogenous line but also applies under g-tensor anisotropy, at least in the case of isotropic rotational diffusion.
Second, from the middle row of the matrix ℬ_{0} (Eq. 63) and after substituting ${\mathit{s}}_{x}^{\mathrm{ss}}$ and ${\mathit{s}}_{z}^{\mathrm{ss}}$ from Eq. (70), we find
where we have introduced the n_{tot}×n_{tot} matrix
and ${\mathbf{P}}_{\mathrm{0}}=\mathbf{P}({\mathit{\omega}}_{I}=\mathrm{0})$. The matrix P_{0} generalizes P_{0} (Eq. 14), and Eq. (72) generalizes Eq. (13) to the case of g-tensor anisotropy.
From the 00th components of ${\mathit{s}}_{y}^{\mathrm{ss}}$ and ${\mathit{s}}_{x}^{\mathrm{ss}}$, we find
where we used the fact that R_{2} is a diagonal matrix. These expressions generalize Eq. (15) to the case of g-tensor anisotropy. The corresponding saturation factor as a function of the offset is then (from Eqs. 16 and 74)
As claimed, to solve for the steady state numerically, we need to invert the matrix P_{0}, whose dimensions are 3 times smaller than those of ℬ_{0}. (The two matrix inversions needed to calculate P_{0} itself involve the diagonal matrices R_{1,2}.)
The cw-EPR spectrum in derivative mode can be calculated from the derivative of P_{0} with respect to the (scalar) frequency offset Δ:
The in-phase and out-of-phase derivative spectra are then obtained from the first (i.e., 00th) components of the vectors
These expressions are used in Sect. 5 to fit the experimental EPR spectra from Fig. 1.
In Fig. 2 we show examples of (integral) EPR spectra calculated using the presented approach for different tumbling times τ_{rot}. The different columns in the figure correspond to different choices of L_{max}. The g-tensor values used in the simulations are characteristic of nitroxide spin labels. We also selected a small mw magnetic field (B_{1}=0.02 G) to mimic the low-power conditions typical of cw-EPR. The main message of this figure is that slower tumbling requires larger L_{max}. At the same time, we see that L_{max}=8 is already good enough for τ_{rot}≤10 ns, which is the range of rotational timescales of relevance to our experimental data (Sect. 5). By selecting L_{max}=10, to be on the safe side, we only need to invert a 21×21 matrix at every frequency offset, which makes the calculation of g-broadened EPR spectra very fast. This allows us to perform an automated search over the various parameters and to fit the experimental cw-EPR spectra in less than a minute.
4.4 Solid effect in “solids”
Extending the above treatment to SE-DNP, we combine the spin dynamics in Eq. (35) with the rotational dynamics in Eq. (44) to form the following SLE:
As before, we introduce the expansion
where we have set n=0 from the start and find
Again, we switch to the symmetrized coefficients
and form the following three n_{tot}-dimensional vectors from the spatial components of the 3D vectors ${\overrightarrow{g}}^{LM}$:
The steady state of the resulting spin dynamics is then
where
generalizes the matrix ℬ_{0} from Eq. (63).
Our goal is to solve for ${\mathit{g}}_{z}^{\mathrm{ss}}$, since its 00th component should be used in Eq. (23) to calculate the rate constants pv_{−} and v_{+}. After inverting ℬ in Eq. (83), we find
Note that now [ℬ^{−1}]_{ij} denotes the n_{tot}×n_{tot} sub-matrix of ℬ^{−1} at position ij and not a scalar matrix element. Using ${\mathit{s}}_{x,y}^{\mathrm{ss}}$ from the previous subsection, we find that the first component of ${\mathit{g}}_{z}^{\mathrm{ss}}$ is
Substituting this result into Eq. (23), we obtain
These expressions, which require the inversion of the 3n_{tot}×3n_{tot} matrix ℬ, are directly generalizable to liquids (Sect. 4.5). In “solids”, it is possible to obtain alternative expressions that require the inversion of a smaller n_{tot}×n_{tot} matrix. To this end, we express ${\mathit{g}}_{x}^{\mathrm{ss}}$ and ${\mathit{g}}_{z}^{\mathrm{ss}}$ in terms of ${\mathit{g}}_{y}^{\mathrm{ss}}$ using the first and last rows of ℬ:
Substituting the first (i.e., 00th) component of ${\mathit{g}}_{z}^{\mathrm{ss}}$ into Eq. (23), we find
Because R_{1} is a diagonal matrix, this result is identical to Eq. (29), showing that ${R}_{\mathrm{1}I}^{A}$ is not affected by the anisotropy of the g tensor in the case of isotropic rotation.
Similarly, from the middle part of ℬ, we obtain
We first observe that $({\mathbf{R}}_{\mathrm{1}}+\text{i}{\mathit{\omega}}_{I}{)}^{-\mathrm{1}}{\mathbf{1}}^{\mathrm{00}}={\mathbf{1}}^{\mathrm{00}}({R}_{\mathrm{1}}+\text{i}{\mathit{\omega}}_{I}{)}^{-\mathrm{1}}$ because R_{1} is diagonal. Then we substitute ${\mathit{s}}_{x,y}^{\mathrm{ss}}$ from before to get
Finally, substituting the 00th element of ${\mathit{g}}_{y}^{\mathrm{ss}}$ into Eq. (30), we find
Observe how these expressions generalize Eqs. (33) and (34) to the case of g-tensor anisotropy.
In the last two rows of Fig. 3 we show ${v}_{+}/{\mathit{\delta}}^{\mathrm{2}}$ and $p{v}_{-}/{\mathit{\delta}}^{\mathrm{2}}$, which have units of time. Although the electronic non-saturation factor p and the rate constant v_{−} always appear together as pv_{−}, it is helpful to separate these two factors when rationalizing SE. We show p and ${v}_{-}/{\mathit{\delta}}^{\mathrm{2}}$ in the first two rows of Fig. 3. Note that ${v}_{+}/{\mathit{\delta}}^{\mathrm{2}}$ and $p{v}_{-}/{\mathit{\delta}}^{\mathrm{2}}$ were calculated directly from Eq. (92), whereas ${v}_{-}/{\mathit{\delta}}^{\mathrm{2}}$ was determined by dividing $p{v}_{-}/{\mathit{\delta}}^{\mathrm{2}}$ by $p=\mathrm{1}-{\mathit{\omega}}_{\mathrm{1}}^{\mathrm{2}}{T}_{\mathrm{1}}[{\mathbf{P}}_{\mathrm{0}}^{-\mathrm{1}}{]}_{\mathrm{11}}$ (Eq. 75).
The columns in Fig. 3 reveal the effect of the g-tensor anisotropy on the different factors relevant to SE. ${v}_{+}/{\mathit{\delta}}^{\mathrm{2}}$ in the third row of the figure is composed of two SE lines centered at −ω_{I} and +ω_{I}. At the fastest tumbling (leftmost column), each of these two lines is symmetric and approximately Lorentzian. When the tumbling slows down, each line broadens and becomes asymmetric. At the slowest tumbling rate (rightmost column), each line resembles a powder EPR spectrum with an anisotropic g tensor. We see that in the regime of slow tumbling the profile of ${v}_{+}/{\mathit{\delta}}^{\mathrm{2}}$ is no longer symmetric (i.e., even) with respect to the electronic resonance at zero offset frequency.
In the second row of Fig. 3 we show ${v}_{-}/{\mathit{\delta}}^{\mathrm{2}}$ (orange line), which is also composed of two SE lines centered at −ω_{I} and +ω_{I}, with the former flipped with respect to the horizontal axis. For comparison, in the second row we also plotted ${v}_{+}/{\mathit{\delta}}^{\mathrm{2}}$ and $-{v}_{+}/{\mathit{\delta}}^{\mathrm{2}}$ (dashed red lines). We see that, for all the tumbling rates, the two SE lines comprising ${v}_{-}/{\mathit{\delta}}^{\mathrm{2}}$ exactly match their counterparts in ${v}_{+}/{\mathit{\delta}}^{\mathrm{2}}$.
The first row of Fig. 3 shows the electronic saturation under g-tensor anisotropy (we actually plot the “non-saturation” $p=\mathrm{1}-s$). Because of the large B_{1} used in the calculations (B_{1}=5.5 G), appreciable electronic saturation is achieved for all the shown tumbling rates. From the perspective of the solid effect, it is noteworthy that the saturation is more localized to on-resonance conditions when the g-tensor anisotropy is averaged out by the tumbling and spreads to larger off-resonance frequencies when the tumbling slows down. This spread broadens the saturation profile and reduces its maximum. However, in spite of the substantial increase in the spectral width of the saturation when going from τ_{rot}=2 ns to τ_{rot}=20 ns, the maximum decreases only moderately, remaining close to 50 % at the slower tumbling rate.
Of course, the amplitude of the saturation profile depends not only on B_{1}, but also on the electronic T_{1} relaxation time. To illustrate this dependence, we recalculated all the curves in Fig. 3 after increasing T_{1} 5-fold to 500 ns. The result, which is shown in Fig. A1, demonstrates larger saturation for all tumbling rates. At the same time, ${v}_{-}/{\mathit{\delta}}^{\mathrm{2}}$ and ${v}_{+}/{\mathit{\delta}}^{\mathrm{2}}$ (second and third rows) remain entirely unaffected. This demonstrates that, in our case of a high constant magnetic field, the SE lines do not experience the power broadening that affects the EPR spectrum.
Finally, the last row of Fig. 3 shows $p{v}_{-}/{\mathit{\delta}}^{\mathrm{2}}$ (solid green line), which equals the product of the first and second rows. From Eq. (7), we know that $p{v}_{-}/{\mathit{\delta}}^{\mathrm{2}}$ basically gives the SE-DNP spectrum up to an overall scaling factor. Since pv_{−} is suppressed by the electronic saturation compared to v_{−}, we see that $p{v}_{-}/{\mathit{\delta}}^{\mathrm{2}}$ is somewhat reduced at offsets between the canonical SE positions ±ω_{I}. Because both the electronic saturation profile and the profile of v_{−} are asymmetric in the slow motional regime where the EPR line exhibits clear g broadening, the line shape of the SE-DNP spectrum (proportional to pv_{−}) is no longer antisymmetric (i.e., odd) with respect to the electronic resonance. This is most visible for the green line in the lower rightmost corner of Fig. 3.
4.5 Solid effect in liquids
In the light of Sect. 3.3, the generalization to liquids consists of calculating the matrix
and using it instead of ℬ^{−1} in Eq. (87):
Because the zz sub-block of ℬ is diagonal and does not couple to the rest when ω_{1}=0, we deduce that
which is identical to Eq. (43). Thus, as we already observed for “solids”, the expression for ${R}_{\mathrm{1}I}^{A}$ is not affected by the anisotropy of the g tensor and the slow tumbling of the radical, in the case of isotropic rotational diffusion.
In Fig. 4 we show the same properties as in Fig. 3 but now in the presence of translational diffusion treated by the FFHS model with the motional timescale τ_{ffhs}=6 ns. Several changes compared to “solids” (Fig. 3) are worth pointing out.
In line with our previous understanding (Sezer, 2023b), the SE lines comprising ${v}_{+}/\langle {\mathit{\delta}}^{\mathrm{2}}\rangle $ are broadened by the translational motion that modulates the dipolar interaction (red lines in the third row of Fig. 4). This motional broadening reduces their maximum intensities compared to “solids” (Fig. 3, third row). Previously, in the case of Lorentzian lines, the reduction in intensity in the transition from solids to liquids was dramatic, by more than a factor of 10 (Sezer, 2023b, Figs. 3, 4 and 5). In contrast, the reduction in the presence of g-tensor broadening is about a factor of 2 (compare the third rows of Figs. 3 and 4). This observation may help rationalize why the maximum SE-DNP enhancement in liquids, e.g., about 50 for trityl in glycerol at 320 K (Kuzhelev et al., 2023), is not negligibly smaller compared to the enhancements that are obtained in the solid state. We also point out that, while reducing the maximum SE intensities in the vicinity of ±ω_{I}, the motional broadening substantially increases the intensities at the smaller offsets around the electronic resonance.
Besides the motional broadening, the progression from left to right in the third row of Fig. 4 demonstrates additional g-tensor broadening, which was also present in “solids”. However, now the two SE lines are affected differently by the g-tensor anisotropy, making the profile of ${v}_{+}/\langle {\mathit{\delta}}^{\mathrm{2}}\rangle $ at slow tumbling rates rather irregular.
Moving on to the second row in Fig. 4, we see that the SE lines that make up ${v}_{-}/\langle {\mathit{\delta}}^{\mathrm{2}}\rangle $ (orange) are now completely different from their counterparts in ${v}_{+}/\langle {\mathit{\delta}}^{\mathrm{2}}\rangle $ (dashed red). The increased intensity in the vicinity of the electronic resonance due to motional broadening is also manifested by ${v}_{-}/\langle {\mathit{\delta}}^{\mathrm{2}}\rangle $. For the fastest tumbling in the figure (leftmost column), the fluctuations of the dipolar interaction not only broaden the SE lines, but also enable a new phenomenon, which is manifested as near-resonance peaks that are comparable in magnitude to the peaks at ±ω_{I} but are clearly distinct from them (orange line). These peaks reflect the multiplicative contribution of the dispersive EPR signal to v_{−} (Sezer, 2023a, b). For faster translational diffusion, the near-resonance peaks may become larger than the peaks at ±ω_{I}, as can be seen in the leftmost column of Fig. A3a (orange line). Because they are more strongly suppressed by the electronic saturation, however, these peaks do not exceed the SE peaks in the final enhancement profile (Fig. A3a, leftmost column, green line).
Up to an overall scaling factor, the green lines in the last row of Fig. 4 correspond to the SE-DNP enhancement profile. Because its middle part is suppressed by the electronic saturation, this profile in the presence of g-tensor broadening becomes very non-symmetric and responds sensitively to the tumbling of the polarizing agent. To further illustrate the influence of the electronic saturation on the SE-DNP spectrum, in Fig. A2 we show the same curves but calculated with a 5-fold longer electronic spin-lattice relaxation time (T_{1}=500 ns), which leads to larger saturation. Similarly, to illustrate the effect of translational diffusion, we recalculated the curves in Fig. 4 for τ_{ffhs}=3 ns (2 times faster) and τ_{ffhs}=12 ns (2 times slower). The results are presented in Fig. A3. These additional simulations show that the SE-DNP line shape is very sensitive to the timescales of molecular motion.
In the next section, we systematically vary the degrees of power broadening and motional broadening to match the experimental DNP profiles from Fig. 1.
Using the developed methodology, we now analyze the experiments from Fig. 1. In the light of Eqs. (2) and (7) for the OE and SE enhancements, we will identify the profile of the electronic saturation (Fig. 4, first row) with ϵ_{OE} and the profile of $p{v}_{-}/\langle {\mathit{\delta}}^{\mathrm{2}}\rangle $ (Fig. 4, last row, green line) with ϵ_{SE}. The tumbling times to be used in the DNP calculations will be obtained by fitting the experimental cw-EPR spectra. We start with 10-Doxyl-PC (Fig. 1a, c) as its experimental spectra were more amenable to unrestricted fits of all parameters.
5.1 Analysis of 10-Doxyl-PC
5.1.1 Fit to the cw-EPR spectrum
Derivative EPR spectra were calculated from the first (i.e., 00th) components of the expressions in Eq. (77) for different values of the fitting parameters. In the fit, we varied the timescale of tumbling, τ_{rot}, as well as the g-tensor anisotropies ${\mathit{\gamma}}_{\mathrm{0}}^{\mathrm{2}}$ and ${\mathit{\gamma}}_{\mathrm{2}}^{\mathrm{2}}$ (Eq. 50). As we have no precise knowledge of the field B_{0} at the sample, we freely shifted the calculated spectra along the horizontal axis to achieve the best match with the experiment. Since this leaves one of the g-tensor components undetermined, we took g_{zz}=2.0023, which is typical for nitroxides.
The numerical integrals of the derivative EPR spectra in Fig. 1a and 1b (dotted–dashed blue lines) do not come down exactly to zero at the end of the integration range at high-frequency offsets. This points to the possibility that the in-phase component, s_{y}, is mixed slightly with the out-of-phase component, s_{x}. To account for this possibility, we fitted the derivative EPR spectra by calculating
where the angle ϕ controlled the degree of mixing.
All in all, not counting the shift along the horizontal axis, we had four fitting parameters: ${\mathit{\gamma}}_{\mathrm{0}}^{\mathrm{2}}$, ${\mathit{\gamma}}_{\mathrm{2}}^{\mathrm{2}}$, τ_{rot} and ϕ. The best fit to the cw-EPR spectrum of 10-Doxyl-PC is shown in Fig. 5a. The corresponding fitting parameters are given in the upper half of Table 1.
Encouragingly, our fitted spectrum shows rather good agreement with the experiment, in spite of the simplifying assumptions of the theoretical model, i.e., isotropic rotational diffusion and the absence of hyperfine interaction. To check the effect of the latter on the cw-EPR spectrum, we used Easyspin (Stoll and Schweiger, 2006) to simulate spectra with our fitted parameters but now also including a nitroxide hyperfine tensor, $A=\text{diag}(\mathrm{14},\mathrm{14},\mathrm{90})$ MHz. The result is given in Fig. A4a. The modification due to the hyperfine interaction, although small as expected at high magnetic fields, is clearly visible. Nevertheless, the comparison of the integrals of the cw-EPR spectra in Fig. A4b suggests that the error made by neglecting the hyperfine interaction when calculating the DNP spectrum should be small.
Regarding the values of the fitted parameters, it was encouraging to see that the fit resulted in a negligibly small mixing angle of $\mathit{\varphi}=-\mathrm{1.3}{}^{\circ}$, indicating that the measured spectrum correctly reflects the in-phase EPR component. With B_{0}=9.4029 T and g_{zz}=2.0023, the fitted g-tensor anisotropies that are given in Table 1 implied that
These values are rather reasonable for a nitroxide spin label. Finally, the fitted timescale of rotational diffusion was τ_{rot}=5.2 ns. For comparison, the same timescale for the nitroxide free radical TEMPOL in water is about 20 ps (Sezer et al., 2009). However, unlike TEMPOL, our spin label is covalently attached to the lipid chain.
5.1.2 Fit to the DNP spectrum
Fixing the g-tensor components and the tumbling time to the values obtained from the fit to the cw-EPR spectrum, we proceeded to fit the DNP spectrum of 10-Doxyl-PC (Fig. 1c). In the calculations, we fixed the mw field to B_{1}=5.5 G, which is our best estimate for the home-built Fabry–Pérot resonator operating at maximum power (Denysenkov et al., 2022). During the fits, we again allowed for global shift of the calculation along the horizontal axis. In addition, we fitted the electronic T_{1} time, which has a direct effect on the electronic saturation profile, as well as the timescale of translational diffusion, τ_{ffhs}, which is responsible for the motional broadening of the SE lines.
In the fit, we calculated the electronic saturation factor (Eq. 75) and the timescale $p{v}_{-}\left(\mathrm{\Delta}\right)/\langle {\mathit{\delta}}^{\mathrm{2}}\rangle $ (last equality in Eq. 94) as functions of the offset frequency Δ. Up to unknown multiplicative factors, these correspond to, respectively, the OE and SE enhancement profiles (Eqs. 2 and 7). We then fit the experimental DNP spectrum by calculating
where the scaling parameters σ_{OE} and σ_{SE} were also allowed to vary freely. As a result, not counting the shift along the horizontal axis, our fit contained four fitting parameters: τ_{ffhs}, T_{1}, σ_{OE} and σ_{SE}. The best fit to the DNP spectrum of 10-Doxyl-PC is shown in Fig. 5b. It is noteworthy how the total DNP enhancement (dashed black line) emerges from the sum of the SE (green line) and OE (dotted–dashed blue line) contributions. The corresponding fitting parameters are given in the bottom half of Table 1.
In the case of 10-Doxyl-PC, the intuitive analysis of Neudert et al. (2017) for identifying the OE and SE components of a mixed DNP spectrum using the integrated cw-EPR line shape already performed very well (Fig. 1c). It is, therefore, not surprising that our analysis, which has more fitting parameters, agrees better with the experimental DNP spectrum (Fig. 5b). Both deficiencies of the intuitive approach, i.e., too narrow OE and SE contributions due to the lack of, respectively, power broadening and motional broadening, appear to be satisfactorily addressed.
On a more fundamental level, our simulation shows that, due to the simultaneous power and motional broadening, the OE and SE contributions to the DNP enhancement are not only rather asymmetric, but also overlap extensively. It should, therefore, be practically impossible to extract any molecular information from the mixed DNP spectrum without a complex, quantitative analysis. In our specific case, the fit resulted in a translational timescale τ_{ffhs}=6.4 ns and suggested that the electronic relaxation time should be about T_{1}=120 ns. At the high magnetic field of the experiment (B_{0}=9.4 T), this spin-lattice relaxation time is practically impossible to measure in the liquid state.
In addition to τ_{ffhs} and T_{1}, the fit to the DNP spectrum of 10-Doxyl-PC also produced the following numerical values for the two scaling parameters in Eq. (98): σ_{OE}=2.4 and σ_{SE}=1.5 ps^{−1}. These will be analyzed in Sect. 5.3 together with the corresponding values for 16-Doxyl-PC.
5.2 Analysis of 16-Doxyl-PC
Because the g-tensor anisotropies are largely averaged in the cw-EPR spectrum of 16-Doxyl-PC (Fig. 1b), we did not attempt to fit them. Instead, we fixed all three components to the values obtained from 10-Doxyl-PC. This left only the rotational time, τ_{rot}, and the mixing angle, ϕ, as fitting parameters, not counting the shift along the horizontal axis. As the automated fitting did not behave well, we varied these two parameters manually. One satisfactory fit, obtained with the parameters that are given in Table 1, is shown in Fig. 6a. We mention that the relative heights of the two lines in the calculation were slightly improved by using a small mixing angle of $\mathit{\varphi}=-\mathrm{2}{}^{\circ}$.
Although, overall, the fit is not bad, the middle part of the calculated spectrum changes too sharply, and its high-frequency line is too narrow compared to the experiment. We again used Easyspin to check whether these deficiencies are due to the lack of hyperfine interaction. The spectra for τ_{rot}=1.9 ns with and without hyperfine interaction are shown in Fig. A5a. As the whole spectrum is narrower than that of 10-Doxyl-PC, the effect of the hyperfine tensor is comparatively larger. Nonetheless, the integrated EPR lines in Fig. A5b show that the extra width due to the hyperfine tensor should not compromise our subsequent analysis of the DNP spectrum, which will experience additional power broadening and motional broadening.
Moving on to the DNP spectrum, we observed that the free fit of all the parameters resulted in a τ_{ffhs} that was more than 2 times larger than that of 10-Doxyl-PC, as we explain below. Considering this to be unrealistic, we fixed τ_{ffhs} to the value that was obtained from 10-Doxyl-PC. Thus, not counting the horizontal translation of the calculated DNP spectrum, our automated fit had three fitting parameters: T_{1}, σ_{OE} and σ_{SE}. The outcome is shown in Fig. 6b. The corresponding parameters are given in the second-last column of the lower half of Table 1.
At 9.4 T the electronic Larmor precession timescale is about 0.5 ps, which is 3 orders of magnitude less than the rotational timescales inferred from the cw-EPR spectra. On such sub-picosecond timescales, the local dynamics of the spin labels at positions 10 and 16 should not be very different from each other. Since the spin-lattice relaxation is determined by dynamics on the electronic Larmor timescale, we were satisfied that the fitted T_{1}=150 ns was close to that from 10-Doxyl-PC.
The performance of the simple analysis of Neudert et al. (2017) was poorer for 16-Doxyl-PC (Fig. 1d). Compared to it, our fit to the DNP enhancement profile is excellent (Fig. 6b). The only parts of the DNP spectrum that our calculation systematically underestimates are the five leftmost experimental points. Although there are other individual experimental points that lie further from the calculated spectrum, these five points are persistently lower by about 0.2 enhancement units.
Observe that the downward shift in the fifth experimental point (together with the first four points) produces an enhancement peak at around −400 MHz. The only way our automated fit can create a pronounced peak at this offset is by making the SE contribution (green line) more “solid-like”, i.e., by increasing τ_{ffhs} and reducing the motional broadening. (The lower-left corner of Fig. A3b provides an example of such a more solid-like SE line shape.) We thus identify the systematic displacement of the leftmost five points as responsible for the increase in τ_{ffhs} when it is allowed to vary freely during the fit.
The best fit that we obtained when τ_{ffhs} was included among the fitting parameters is shown in Fig. 7. (The resulting fit parameters are given in the last column of the lower half of Table 1.) Indeed, with τ_{ffhs}=15.3 ns, the SE lines (green) have become sharper, and a small enhancement peak at −400 MHz has emerged (dashed black line). Although the enhancement around +400 MHz has been compromised in the process, the overall fit to all the experimental points is improved compared to Fig. 6b.
The two alternative fits in Figs. 6b and 7 correspond to very different timescales of translational diffusion. Nevertheless, within the variability of the measurements, they both agree with the DNP data. Considering the experimental challenges of liquid-state DNP at such high magnetic fields and large mw powers, further decreasing the experimental variability will be very hard. It is, therefore, important to analyze together several different experimental constructs, like our 10- and 16-Doxyl-PC. The final decision of which fit to the DNP spectrum of 16-Doxyl-PC is “better” can only be based on the overall consistency of the fitted parameters across all the analyzed data. We return to this point in Sect. 5.3.
The other two parameters that emerged from the fit to the DNP spectrum of 16-Doxyl-PC were σ_{SE} and σ_{OE}. These determine the amplitudes of the SE contribution (solid green lines in Figs. 6b and 7) and OE contribution (dotted–dashed blue lines) to the DNP enhancement (dashed black lines). We now turn to the analysis of these scaling parameters.
5.3 Additional molecular parameters
Ultimately, the motivation to disentangle a mixed DNP spectrum into its OE and SE components lies in the desire to extract information about the molecular and spin properties that the respective DNP mechanism depends on. The main advantage of our procedure over the intuitive approach of Neudert et al. (2017) is that our decomposition produces physically interpretable parameters like τ_{ffhs} and T_{1}. In addition, our scaling parameters σ_{OE} and σ_{SE} multiply, respectively, the saturation factor and ${v}_{-}/\langle {\mathit{\delta}}^{\mathrm{2}}\rangle $, whose absolute magnitudes are part of the calculation (Fig. 4, vertical axes). Thus, we can extract further information from the fitted values of σ_{OE} and σ_{SE}. In contrast, because the simple approach rescales the integrated cw-EPR spectrum whose amplitude is arbitrary, the values of its scaling factors are not informative.
Using Eq. (2) for the OE enhancement, the coupling factor c is readily expressed in terms of σ_{OE}:
where the leakage factor f can be obtained by measuring the nuclear spin-lattice relaxation times (Eq. 4).
We measured the T_{1} values for the chain protons of DOPC (without spin-labeled lipids) at 310 and 330 K using the Fabry–Pérot probe. These are given in the ${T}_{\mathrm{1}I}^{\mathrm{0}}$ column of Table 2. Additionally, we measured the nuclear spin-lattice relaxation times in the presence of either 10- or 16-Doxyl-PC (column T_{1I} of Table 2). The target temperature of the DNP experiments (320 K) lies between the two temperatures at which the nuclear T_{1} times were measured. However, considering the possibility of mild temperature rise by several degrees, we expect the values at 330 K to closely reflect the DNP conditions. Nonetheless, we carry out the following analysis using the T_{1} values measured at both 310 and 330 K.
The leakage factors obtained from Eq. (4) are shown in column f of Table 2. Using the values of σ_{OE} from Table 1 in Eq. (99), we arrived at the coupling factors in column c of Table 2. In the case of 16-Doxyl-PC, the analysis was performed for the fit where τ_{ffhs} was fixed at 6.4 ns (denoted as 16 in Table 2) as well as for the fit where τ_{ffhs} was free to change (denoted as 16*). (These two alternatives correspond to the last two columns of Table 1.) For both choices, somewhat larger coupling factors were deduced for 16-Doxyl-PC compared to 10-Doxyl-PC. The estimated coupling factors are less than 2 times smaller than what we have obtained previously for TEMPOL in DMSO and about 4 times smaller than the coupling factors between TEMPOL and the protons of toluene (Prisner et al., 2016; Sezer, 2013; Küçük et al., 2015).
Turning now to SE, using the enhancement in Eq. (7), we express the unknown strength of the dipolar interaction in terms of the scaling parameter σ_{SE} as follows:
The values of 〈δ^{2}〉 T_{1I}, which were calculated from the right-hand side of Eq. (100), are about 2 ns^{−1} for 10, 16 and 16* (Table 2). Since ${v}_{+}/\langle {\mathit{\delta}}^{\mathrm{2}}\rangle $ is about 1 ps (Fig. 4, third row), we conclude that ${v}_{+}{T}_{\mathrm{1}I}\ll \mathrm{1}$, which justifies our use of the approximation in Eq. (7) throughout the analysis, including during the fit to the DNP spectra.
From the expression of 〈δ^{2}〉 (Eq. 39), we can write the contact distance of the translational FFHS model as
where N is the number density of the electronic spins. Since, in principle, all parameters on the right-hand side of Eq. (101) are measurable, we can determine b. To estimate N, we note that the molecular volume of DOPC is 1.3 nm^{3} (Greenwood et al., 2006). Since there are 20 unlabeled lipids for 1 labeled one, we estimate N = (20×1.3 nm^{3})^{−1}, which corresponds to a molar concentration of 64 mM. Using this number in Eq. (101), we obtained the values of b that are given in the second-last column of Table 2.
When the values of b are interpreted literally as the “contact distance” between the nitroxide spin label and the protons of the lipid chains, their substantial variation between 10- and 16-Doxyl-PC is disturbing. From that perspective, it is clear that the parameter b of the FFHS model, which we used to account for the fluctuations of the dipolar interaction due to molecular translations, cannot reflect the actual molecular distances of the closest approach.
Because b was obtained from the scaling parameter σ_{SE}, only information about the amplitude of the SE enhancement has been directly used in its estimate. In contrast, the motional timescale τ_{ffhs} (Table 1) encodes information about the line shape of the SE enhancement. From these complementary features of the SE contribution to the DNP spectrum, we have managed to determine both b and τ_{ffhs}. Having access to these two parameters, we can calculate the diffusion constant of the FFHS model from Eq. (38). The results are given in the last column of Table 2. To our surprise, we obtained very similar values for 10 and 16*, while the diffusion constant for 16 is 2-fold larger. (Given the variability in the experimental data and the fact that the fits to the DNP spectra are not unique, the differences between D_{ffhs} of 10 and 16* should not be seen as meaningful.)
In an effort to identify a potential candidate for the physical motion that the FFHS model emulates, we observed that the coefficients of lateral translational diffusion for DOPC in oriented bilayers are 20 nm^{2} µs^{−1} at 323 K and 26 nm^{2} µs^{−1} at 333 K (Filippov et al., 2003, Fig. 6a). These, we expect, will bracket the value under our DNP conditions. The diffusion in the FFHS model corresponds to the relative translation of the nuclear and electronic spins, i.e., ${D}_{\mathrm{ffhs}}={D}_{I}+{D}_{S}$. Assuming that the lateral diffusion of spin-labeled PSPC in a DOPC bilayer is similar to that of DOPC, from the measured values given above we would expect D_{ffhs} to be between 40 and 52 nm^{2} µs^{−1}. This range is surprisingly close to the estimates of 10 and 16* in the last column of Table 2, which suggests that the FFHS model in our analysis likely accounts for the lateral diffusion of the lipids in the plane of the bilayer.
Since it leads to a diffusion constant that is similar to (i) the known lateral diffusion of DOPC and (ii) the estimate obtained for 10-Doxyl-PC, we conclude that the fit to the DNP spectrum of 16-Doxyl-PC that is shown in Fig. 7 (i.e., the one that led to “unreasonably” large τ_{ffhs}) is more realistic than the one with fixed τ_{ffhs} (Fig. 6b). From the perspective of the diffusion constant, the longer translational timescale of 16* compared to 10, which resulted in a more solid-like SE line shape with less motional broadening, reflects the fact that the “contact distances” in the two cases are different. In retrospect, it is amazing how the independent estimates of b and τ_{ffhs} combine to yield practically identical diffusion constants for the two spin-labeling positions.
At the moment, it is not clear to us how to properly interpret the different values of b at positions 10 and 16. Atomistic molecular dynamics simulations (Oruç et al., 2016) could, in principle, be used to investigate whether these effective contact distances reflect differences in proton density along the normal of the lipid bilayer or arise for some other reason.
5.4 Limitations of the modeling
The calculated DNP spectra of 10-Doxyl-PC (Fig. 5b) and 16-Doxyl-PC (Fig. 7) agree well with the experiments in spite of our simplistic treatment of the quantum and classical dynamics. Specifically, when modeling the spin dynamics, (i) we completely neglected the hyperfine interaction with the nuclear spin of ^{14}N, which is present in nitroxide spin labels. In the case of the classical dynamics, (ii) we modeled the reorientation of the spin labels at positions 10 and 16 of the lipid chain as free, isotropic diffusion, and (iii) we modeled the dynamics of the acyl protons relative to the unpaired electron as isotropic translational diffusion that extends to infinity in all three spatial directions. We now comment on these deficiencies of the modeling.
Starting with the third point, it is clear that the translational diffusion of the polarized aliphatic protons (as well as that of the chain-attached spin labels) must be confined to the interior of the lipid bilayer and should not extend arbitrarily far along the direction perpendicular to the bilayer plane. In contrast, the FFHS model whose analytical correlation function we used in the calculations assumes isotropic diffusion in all spatial directions. To properly address this deficiency of the modeling, one would need to solve the diffusion equation with boundary conditions that reflect the confining planar geometry of the lipid bilayer and then calculate the dipolar correlation function for such confined diffusion (preferably in closed, analytical form). In the meantime, one could argue that, because the dipolar interaction drops rapidly with distance, an overwhelming contribution to the dipolar correlation function should come from configurations in which the electron and nucleus are close to each other. In that case, the unphysical configurations that place the acyl protons outside the plane of the lipid bilayer (but are allowed in the FFHS model) may contribute relatively little. To support this argument, we observe that the numerical values of the FFHS parameter b in Table 2 indicate that the shortest relevant distances for SE are about 0.6 nm (for 10-Doxyl-PC) and 0.9 nm (for 16-Doxyl-PC). These are 3 to 5 times smaller than the hydrophobic thickness of the DOPC lipid bilayer, which is about 3 nm (Kučerka et al., 2008).
Regarding the second deficiency, the problem here is that the nitroxide spin label is covalently fused to the lipid chain; thus, its possible orientations should reflect the preferred alignment of the chain in the hydrophobic core of the bilayer. Furthermore, the fused nitroxide is not expected to have identical diffusion rates for rotations about different spatial directions. Clearly, both of these aspects (i.e., the orientational preference and the anisotropy) are missing from the free, isotropic rotational diffusion that we implemented. It is, however, known how to account for them in a rigorous and efficient way. Indeed, the MOMD model from the Freed lab (Meirovitch et al., 1984) treats anisotropic rotational diffusion in a restoring potential. In fact, this model has been extensively used to simulate high-field cw-EPR spectra of lipid bilayers by Freed (Lou et al., 2001; Costa-Filho et al., 2003) and Marsh (Livshits et al., 2004, 2006). The studies of Marsh and colleagues have focused on DMPC lipid bilayers containing appreciable amounts of cholesterol, which puts them in a liquid-ordered phase. For DMPC with 40 mol % cholesterol at 30 ^{∘}C, 10-Doxyl-PC was deduced to be aligned with the director (i.e., the direction normal to the bilayer plane) with an order parameter S=0.67 (Livshits et al., 2004). If the orientational motion is imagined as being confined to a cone (Lipari and Szabo, 1982), this order parameter would correspond to a maximum possible deviation from the director of ${\mathit{\theta}}_{\mathrm{0}}=\mathrm{40}{}^{\circ}$ in all directions. The lipid bilayers in our experiments are composed of pure DOPC lipids and are in their liquid-crystalline phase, where the ordering is substantially reduced. The liquid-crystalline phase of pure DPPC lipid bilayers has been characterized in the studies of Freed and colleagues. The order parameter reported for 16-Doxyl-PC in pure DPPC at 50 ^{∘}C is S=0.16 (Costa-Filho et al., 2003). It corresponds to a maximum possible deviation from the director of ${\mathit{\theta}}_{\mathrm{0}}=\mathrm{75}{}^{\circ}$, assuming the diffusion is confined to a cone. Although 10-Doxyl-PC is expected to be more ordered than 16-Doxyl-PC, it is not clear how much smaller than ${\mathit{\theta}}_{\mathrm{0}}=\mathrm{75}{}^{\circ}$ its corresponding cone angle would be. (Because S is the expectation value of a rank-2 spherical harmonic, the free rotational diffusion that we use corresponds to ${\mathit{\theta}}_{\mathrm{0}}=\mathrm{90}{}^{\circ}$.) From these studies we conclude that the MOMD model (with an axial diffusion tensor) will likely improve our fits to the experimental cw-EPR spectra. Nevertheless, free rotation may still be a good first approximation to the orientational dynamics of 10- and 16-Doxyl-PC in the liquid-crystalline phase of our lipid bilayers.
We should emphasize that our aim in the current paper is to show how to account for the rotational dynamics of the polarizing agent in the calculation of SE-DNP. In this context, we observe that, while the cw-EPR spectra in derivative mode (Figs. A4a and A5a) are extremely sensitive to the details of the rotational motion of the radical, their integrals (Figs. A4b and A5b) are much more forgiving. When contributing to the DNP spectrum, these integrated EPR line shapes are additionally broadened by mw power (OE) and translational diffusion (SE) (Figs. 5b and 7). All these factors are expected to reduce the sensitivity of the DNP spectrum to the details of the radical tumbling (at least in comparison to the sensitivity of the cw-EPR line shape). We therefore think that, for the purposes of fitting the DNP spectrum, further improving the fit to the cw-EPR spectra at the cost of introducing more fitting parameters is not really justified. That being said, we stress that the formalism of Sect. 4 can be straightforwardly extended to anisotropic diffusion in an orienting potential (i.e., the MOMD model). This would lead to larger matrices R_{1}, R_{2} and Δ, whose matrix elements would be different than the expressions we gave in Sect. 4.2 for free, isotropic rotational diffusion.^{5} Once correctly formed, these three matrices can be directly used in Eqs. (92) and (94) to calculate the SE-DNP spectra in, respectively, “solids” and liquids.
Moving on to the first deficiency mentioned above, we remind the reader that we describe the SE spin dynamics in terms of two sets of Bloch equations that are connected in series (Sezer, 2023a). These are the classical Bloch equations with the Bloch matrix 𝔹_{0} (Sect. 3.1) and the “new Bloch equations” (Eq. 21) with the matrix $\mathbb{B}={\mathbb{B}}_{\mathrm{0}}+\text{i}{\mathit{\omega}}_{I}$ (Sect. 3.2). Because our description of SE-DNP is based on Bloch equations, in Sect. 4.1 and 4.2 we reformulated Freed's treatment of slow tumbling as a generalization of the classical Bloch equations, such that the scalar elements of 𝔹_{0} became matrices in the space of the angular-momentum indices LM. The result was the “expanded Bloch matrix” ℬ_{0} in Eq. (63). For this reformulation to work, however, we had to neglect the hyperfine interaction, which is in fact treated by Freed et al. (1971). As a result, our analysis is formally deficient for nitroxide radicals. Nevertheless, we reasoned that it should be possible to illustrate the theoretical formalism in its current form by focusing on nitroxides at high magnetic fields, where the hyperfine interaction is expected to be negligible compared to the anisotropy of the g tensor. From this perspective, it should be clear that the DNP experiments that we analyzed here had been carefully selected.
Figures A4 and A5 show our attempt to assess the contribution of the neglected hyperfine interaction to the (integrated) EPR spectra of 10- and 16-Doxyl-PC. A somewhat more detailed analysis is contained in our response to the reviewers, which is freely accessible online. There we observe that the hyperfine interaction slightly broadens the EPR line of 16-Doxyl-PC (which is also visible in Fig. A5a). Since the only mechanism of broadening in our case is the rotational tumbling, our choice of τ_{rot}=1.9 ns (Table 1) likely compensates for some of the “missing” hyperfine broadening. Such compensation does not appear to be happening in the case of 10-Doxyl-PC, where the hyperfine interaction changes the shape but not the width of the EPR line (Fig. A4a). Ultimately, for the theory to be applicable to SE-DNP with nitroxide polarizing agents at lower magnetic fields, like X band (Gizatullin et al., 2021a, b), the description of the spin dynamics will need to include the nuclear spin of ^{14}N. Since the dimension of the resulting Liouville space would need to increase by a factor of 9, it should be possible to preserve the two sets of connected Bloch equations after replacing each of their scalar matrix elements with a 9×9 matrix. Alternatively, the two sets of Bloch equations should be replaced with the corresponding equations of motion for the density matrix in Liouville space. However, considering the inherent experimental uncertainty of the DNP enhancements that we compare with (Fig. 1c and 1d) and the achieved agreement between simulation and experiment (Figs. 5b and 7), we believe that such more complex modeling is presently not justified.
Once the spin dynamics of the solid effect has been formulated in the time domain (Sezer, 2023a), it becomes possible to interface this quantum dynamics with various types of classical dynamics. The classical dynamics in Sezer (2023b) was the translational diffusion of the spins in a liquid; here we additionally included the rotational diffusion of the polarizing agent. To illustrate the practical utility of the resulting formalism, we analyzed either previously published (Sezer, 2023b) or previously unpublished (current paper) experimental DNP data on lipid bilayers. In our analysis, the treatment of molecular translation and rotation was limited to the simplest possible models of free, isotropic diffusion. Surprisingly, in spite of the spatial anisotropy that one expects for hydrated lipid bilayers, previously we found that isotropic translation, as described by the FFHS model, worked well for the free radical BDPA in DMPC bilayers (Sezer, 2023b; Kuzhelev et al., 2022). Similarly, in the current paper we found that the simplest treatment of free, isotropic rotation (together with FFHS translation) reproduced well the DNP field profiles of nitroxide-labeled lipids in DOPC bilayers.
DNP experiments with nitroxide free radicals in viscous liquids invariably manifest a mixture of SE and OE (Leblond et al., 1971b; Neudert et al., 2017). As these two DNP mechanisms are sensitive to molecular motions on vastly different timescales, it should be possible to obtain rich dynamical information by analyzing their contributions to the overall DNP enhancement. Disentangling the SE and OE contributions, however, has proven to be challenging (Leblond et al., 1971a). Here we fitted liquid-state DNP spectra by calculating enhancements that were affected by both the translational diffusion of the spins and the rotational diffusion of the free radical. Since different motions modify the amplitude and the shape of the DNP spectrum in a highly concerted manner, by fitting the entire line shape of the enhancement, we also gained access to the absolute magnitudes of the SE and OE contributions.
Our current treatment of SE-DNP in liquids uses only the correlation function of the dipolar interaction to describe the translational motion of the spins (Sezer, 2023b). This is formally correct only when the diffusion is much faster than the nuclear T_{1} relaxation. It should be possible to relax this condition and model slower spin diffusion, as relevant for SE in the solid state.
The analyzed experimental data and the code used to generate the figures in the paper are available at https://github.com/dzsezer/solidDNPliquids_g-tensor (last access: 23 October 2023) (https://doi.org/10.5281/zenodo.8360325, Sezer, 2023c).
TFP envisioned the presented high-field DNP experiments in lipid bilayers and acquired funding for their execution. DD carried out these EPR and NMR experiments, deconvoluted the measured NMR signals to calculate the presented DNP spectra for the lipid chain protons, and performed the reported Easyspin simulations. DS conceived of the presented analysis of the experimental data, developed the reported theoretical and numerical framework, analyzed the data, and wrote the manuscript with feedback from all the coauthors.
At least one of the (co-)authors is a member of the editorial board of Magnetic Resonance. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
We are grateful to Vasyl Denysenkov, without whom the high-field DNP experiments analyzed here would not have been possible. Andrei Kuzhelev was the first to realize that the DNP spectra presented in this paper manifest not only the Overhauser effect, but also the solid effect.
This research has been supported by the Deutsche Forschungsgemeinschaft (grant no. 405972957).
This open-access publication was funded by the Goethe University Frankfurt.
This paper was edited by Geoffrey Bodenhausen and reviewed by Gunnar Jeschke and two anonymous referees.
Abragam, A.: The Principles of Nuclear Magnetism, Oxford University Press, New York, ISBN 978 0 19 852014 6, 1961. a, b
Abragam, A. and Goldman, M.: Principles of dynamic nuclear polarisation, Rep. Prog. Phys., 41, 395, https://doi.org/10.1088/0034-4885/41/3/002, 1978. a
Abragam, A. and Proctor, W. G.: Une nouvelle méthode de polarisation dynamique des noyaux atomiques dans les solides, Compt. Rend., 246, 2253–2256, 1958. a, b
Anderson, P. W.: A Mathematical Model for the Narrowing of Spectral Lines by Exchange or Motion, J. Phys. Soc. Jpn., 9, 316–339, https://doi.org/10.1143/JPSJ.9.316, 1954. a, b
Atsarkin, V. A.: Dynamic nuclear polarization: Yesterday, today, and tomorrow, J. Phys. Conf. Ser., 324, 012003, https://doi.org/10.1088/1742-6596/324/1/012003, 2011. a
Ayant, Y., Belorizky, E., Alizon, J., and Gallice, J.: Calcul des densités spectrales résultant d'un mouvement aléatoire de translation en relaxation par interaction dipolaire magnétique dans les liquides, J. Phys. (Paris), 36, 991–1004, 1975. a, b, c
Carver, T. R. and Slichter, C. P.: Experimental Verification of the Overhauser Nuclear Polarization Effect, Phys. Rev., 102, 975–980, https://doi.org/10.1103/PhysRev.102.975, 1956. a
Costa-Filho, A. J., Shimoyama, Y., and Freed, J. H.: A 2D-ELDOR Study of the Liquid Ordered Phase in Multilamellar Vesicle Membranes, Biophys. J., 84, 2619–2633, https://doi.org/10.1016/S0006-3495(03)75067-X, 2003. a, b
Denysenkov, V., Dai, D., and Prisner, T. F.: A triple resonance (e, ^{1}H, ^{13}C) probehead for liquid-state DNP experiments at 9.4 Tesla, J. Magn. Reson., 337, 107185, https://doi.org/10.1016/j.jmr.2022.107185, 2022. a, b, c
Erb, E., Motchane, J.-L., and Uebersfeld, J.: Effet de polarisation nucléaire dans les liquides et les gaz adsorbés sur les charbons, Compt. Rend., 246, 2121–2123, 1958a. a
Erb, E., Motchane, J.-L., and Uebersfeld, J.: Sur une nouvelle méthode de polarisation nucléaire dans les fluides adsorbés sur les charbons. extension aux solides et en particulier aux substances organiques irradiées., Compt. Rend., 246, 3050–3052, 1958b. a
Filippov, A., Orädd, G., and Lindblom, G.: The Effect of Cholesterol on the Lateral Diffusion of Phospholipids in Oriented Bilayers, Biophys. J., 84, 3079–3086, https://doi.org/10.1016/S0006-3495(03)70033-2, 2003. a
Freed, J. H.: Theory of Slow Tumbling ESR Spectra for Nitroxides, in: Spin Labeling: Theory and Applications, edited by: Berliner, L. J., Academic Press, Inc., 53–132, https://doi.org/10.1016/B978-0-12-092350-2.50008-4, 1976. a
Freed, J. H., Bruno, G. V., and Polnaszek, C. F.: Electron spin resonance line shapes and saturation in the slow motional region, J. Phys. Chem., 75, 3385–3399, https://doi.org/10.1021/j100691a001, 1971. a, b, c, d, e, f, g, h, i, j, k, l
Gizatullin, B., Mattea, C., and Stapf, S.: Field-cycling NMR and DNP – A friendship with benefits, J. Magn. Reson., 322, 106851, https://doi.org/10.1016/j.jmr.2020.106851, 2021a. a, b, c
Gizatullin, B., Mattea, C., and Stapf, S.: Molecular Dynamics in Ionic Liquid/Radical Systems, J. Phys. Chem. B, 125, 4850–4862, https://doi.org/10.1021/acs.jpcb.1c02118, 2021b. a, b
Gizatullin, B., Mattea, C., and Stapf, S.: Three mechanisms of room temperature dynamic nuclear polarization occur simultaneously in an ionic liquid, Phys. Chem. Chem. Phys., 24, 27004–27008, https://doi.org/10.1039/D2CP03437A, 2022. a
Greenwood, A. I., Tristram-Nagle, S., and Nagle, J. F.: Partial molecular volumes of lipids and cholesterol, Chem. Phys. Lipids, 143, 1–10, https://doi.org/10.1016/j.chemphyslip.2006.04.002, 2006. a
Hausser, K. H. and Stehlik, D.: Dynamic Nuclear Polarization in Liquids, in: Advances in Magnetic Resonance, edited by: Waugh, J. S., vol. 3, 79–139 pp., Academic Press, https://doi.org/10.1016/B978-1-4832-3116-7.50010-2, 1968. a, b
Hwang, L.-P. and Freed, J. H.: Dynamic effects of pair correlation functions on spin relaxation by translational diffusion in liquids, J. Chem. Phys., 63, 4017–4025, 1975. a, b, c
Kubo, R.: Note on the Stochastic Theory of Resonance Absorption, J. Phys. Soc. Jpn., 9, 935–944, https://doi.org/10.1143/JPSJ.9.935, 1954. a, b
Kubo, R.: A Stochastic Theory of Line Shape, in: Advances in Chemical Physics, edited by: Shuler, K. E., vol. XV, chap. 6, John Wiley & Sons Inc., 101–127, https://doi.org/10.1002/9780470143605.ch6, 1969. a
Kučerka, N., Nagle, J. F., Sachs, J. N., Feller, S. E., Pencer, J., Jackson, A., and Katsaras, J.: Lipid Bilayer Structure Determined by the Simultaneous Analysis of Neutron and X-Ray Scattering Data, Biophys. J., 95, 2356–2367, https://doi.org/10.1529/biophysj.108.132662, 2008. a
Küçük, S. E., Neugebauer, P., Prisner, T. F., and Sezer, D.: Molecular simulations for dynamic nuclear polarization in liquids: a case study of TEMPOL in acetone and DMSO, Phys. Chem. Chem. Phys., 17, 6618–6628, https://doi.org/10.1039/C4CP05832A, 2015. a
Kuzhelev, A. A., Dai, D., Denysenkov, V., and Prisner, T. F.: Solid-like Dynamic Nuclear Polarization Observed in the Fluid Phase of Lipid Bilayers at 9.4 T, J. Am. Chem. Soc., 144, 1164–1168, https://doi.org/10.1021/jacs.1c12837, 2022. a, b
Kuzhelev, A. A., Denysenkov, V., Ahmad, I. M., Rogozhnikova, O. Y., Trukhin, D. V., Bagryanskaya, E. G., Tormyshev, V. M., Sigurdsson, S. T., and Prisner, T. F.: Solid-Effect Dynamic Nuclear Polarization in Viscous Liquids at 9.4 T Using Narrow-Line Polarizing Agents, J. Am. Chem. Soc., 145, 10268–10274, https://doi.org/10.1021/jacs.3c01358, 2023. a, b, c
Leblond, J., Papon, P., and Korringa, J.: Stochastic Theory of Dynamic Spin Polarization in Viscous Liquids with Anisotropic Electron-Spin Relaxation, Phys. Rev. A, 4, 1539–1549, https://doi.org/10.1103/PhysRevA.4.1539, 1971a. a, b
Leblond, J., Uebersfeld, J., and Korringa, J.: Study of the Liquid-State Dynamics by Means of Magnetic Resonance and Dynamic Polarization, Phys. Rev. A, 4, 1532–1539, https://doi.org/10.1103/PhysRevA.4.1532, 1971b. a, b, c
Lipari, G. and Szabo, A.: Model-free approach to the interpretation of nuclear magnetic resonance relaxation in macromolecules. 2. Analysis of experimental results, J. Am. Chem. Soc., 104, 4559–4570, https://doi.org/10.1021/ja00381a010, 1982. a
Livshits, V. A., Kurad, D., and Marsh, D.: Simulation Studies on High-Field EPR Spectra of Lipid Spin Labels in Cholesterol-Containing Membranes, J. Phys. Chem. B, 108, 9403–9411, https://doi.org/10.1021/jp035915p, 2004. a, b
Livshits, V. A., Kurad, D., and Marsh, D.: Multifrequency simulations of the EPR spectra of lipid spin labels in membranes, J. Magn. Reson., 180, 63–71, https://doi.org/10.1016/j.jmr.2006.01.005, 2006. a
Lou, Y., Ge, M., and Freed, J. H.: A Multifrequency ESR Study of the Complex Dynamics of Membranes, J. Phys. Chem. B, 105, 11053–11056, https://doi.org/10.1021/jp013226c, 2001. a
Meirovitch, E., Nayeem, A., and Freed, J. H.: Analysis of protein-lipid interactions based on model simulations of electron spin resonance spectra, J. Phys. Chem., 88, 3454–3465, https://doi.org/10.1021/j150660a018, 1984. a, b, c
Müller-Warmuth, W. and Meise-Gresch, K.: Molecular Motions and Interactions as Studied by Dynamic Nuclear Polarization (DNP) in Free Radical Solutions, in: Advances in Magnetic Resonance, edited by: Waugh, J. S., vol. 11, 1–45, Academic Press, https://doi.org/10.1016/B978-0-12-025511-5.50007-4, 1983. a, b
Neudert, O., Mattea, C., and Stapf, S.: A compact X-Band resonator for DNP-enhanced Fast-Field-Cycling NMR, J. Magn. Reson., 271, 7–14, https://doi.org/10.1016/j.jmr.2016.08.002, 2016. a
Neudert, O., Mattea, C., and Stapf, S.: Molecular dynamics-based selectivity for Fast-Field-Cycling relaxometry by Overhauser and solid effect dynamic nuclear polarization, J. Magn. Reson., 276, 113–121, https://doi.org/10.1016/j.jmr.2017.01.013, 2017. a, b, c, d, e, f
Oruç, T., Küçük, S. E., and Sezer, D.: Lipid bilayer permeation of aliphatic amine and carboxylic acid drugs: rates of insertion, translocation and dissociation from MD simulations, Phys. Chem. Chem. Phys., 18, 24511–24525, https://doi.org/10.1039/C6CP05278A, 2016. a
Papon, P., Motchane, J. L., and Korringa, J.: Unitary Theory of Dynamic Polarization of Nuclear Spins in Liquids and Solids, Phys. Rev., 175, 641–649, https://doi.org/10.1103/PhysRev.175.641, 1968. a
Polimeno, A. and Freed, J. H.: Slow Motional ESR in Complex Fluids: The Slowly Relaxing Local Structure Model of Solvent Cage Effects, J. Phys. Chem., 99, 10995–11006, https://doi.org/10.1021/j100027a047, 1995. a, b, c
Prisner, T., Denysenkov, V., and Sezer, D.: Liquid state DNP at high magnetic fields: Instrumentation, experimental results and atomistic modelling by molecular dynamics simulations, J. Magn. Reson., 264, 68–77, https://doi.org/10.1016/j.jmr.2015.11.004, 2016. a
Redfield, A. G.: On the Theory of Relaxation Processes, IBM J. Res. Dev., 1, 19–31, https://doi.org/10.1147/rd.11.0019, 1957. a
Schneider, D. J. and Freed, J. H.: Calculating Slow Motional Magnetic Resonance Spectra, in: Spin Labeling: Theory and Applications, edited by: Berliner, L. J. and Reuben, J., Springer US, Boston, MA, 1–76, https://doi.org/10.1007/978-1-4613-0743-3_1, 1989. a
Sezer, D.: Computation of DNP coupling factors of a nitroxide radical in toluene: seamless combination of MD simulations and analytical calculations, Phys. Chem. Chem. Phys., 15, 526–540, https://doi.org/10.1039/C2CP42430D, 2013. a
Sezer, D.: Non-perturbative treatment of the solid effect of dynamic nuclear polarization, Magn. Reson., 4, 129–152, https://doi.org/10.5194/mr-4-129-2023, 2023a. a, b, c, d, e, f, g, h, i, j, k, l, m
Sezer, D.: The solid effect of dynamic nuclear polarization in liquids, Magn. Reson., 4, 153–174, https://doi.org/10.5194/mr-4-153-2023, 2023b. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s, t
Sezer, D.: dzsezer/solidDNPliquids_g-tensor: v1.0.0 (v1.0.0), Zenodo [code and data set], https://doi.org/10.5281/zenodo.8360325, 2023c. a
Sezer, D., Gafurov, M., Prandolini, M. J., Denysenkov, V. P., and Prisner, T. F.: Dynamic nuclear polarization of water by a nitroxide radical: rigorous treatment of the electron spin saturation and comparison with experiments at 9.2 Tesla, Phys. Chem. Chem. Phys., 11, 6638–6653, https://doi.org/10.1039/B906719C, 2009. a
Solomon, I.: Relaxation Processes in a System of Two Spins, Phys. Rev., 99, 559–565, https://doi.org/10.1103/PhysRev.99.559, 1955. a
Stoll, S. and Schweiger, A.: EasySpin, a comprehensive software package for spectral simulation and analysis in EPR, J. Magn. Reson., 178, 42–55, https://doi.org/10.1016/j.jmr.2005.08.013, 2006. a, b
Wenckebach, T.: Essentials of Dynamic Nuclear Polarization, Spindrift Publications, the Netherlands, ISBN 9789075541182, 2016. a
Abbreviations used in the text: 1,3-bisdiphenylene-2-phenylallyl (BDPA), continuous wave (cw), 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), double quantum (DQ), dynamic nuclear polarization (DNP), electron paramagnetic resonance (EPR), force-free hard sphere (FFHS), microscopic order macroscopic disorder (MOMD), microwave (mw), nuclear magnetic resonance (NMR), Overhauser effect (OE), 1-palmitoyl-2-stearoyl-sn-glycero-3-phosphocholine (PSPC), solid effect (SE), stochastic Liouville equation (SLE), slowly relaxing local structure (SRLS), 4-Hydroxy-2,2,6,6-tetramethylpiperidin-1-oxyl (TEMPOL), zero quantum (ZQ).
Surprisingly, this improved treatment is not mentioned by Müller-Warmuth and Meise-Gresch (1983), who continue to use the older, deficient expression of spectral density for translational diffusion.
In Sezer (2023a, b), the frequency offset was denoted by Ω. Here we denote the offset by Δ and reserve the symbol Ω for the orientation of the polarizing agent (Sect. 4).
We follow Freed et al. (1971) and consider the effect of the g-tensor anisotropy only on the secular terms in the electronic spin Hamiltonian, i.e., those proportional to the spin operator S_{z}. The response of the non-secular terms to the g anisotropy is neglected.
The orienting potential will mix coefficients with different values of L, which will result in non-diagonal R_{1} and R_{2}. The anisotropic rotation will mix coefficients whose M indices differ by ±1, so odd values of M will also need to be included. Finally, since the potential is defined with respect to the director axis, which may differ from the axes of both the laboratory frame and the molecular frame, it will be necessary to consider coefficients ${\mathit{s}}_{MN}^{L}$ with non-zero index N. Clearly, for a given L_{max}, the resulting matrices R_{1}, R_{2} and Δ will be substantially larger. A detailed presentation can be found in Schneider and Freed (1989) and Polimeno and Freed (1995).
- Abstract
- Introduction
- Motivation
- Theoretical background
- Slow-motional EPR and DNP spectra for an anisotropic g tensor
- Disentangling the solid and Overhauser DNP effects
- Conclusion
- Appendix A: Additional figures
- Code and data availability
- Author contributions
- Competing interests
- Disclaimer
- Acknowledgements
- Financial support
- Review statement
- References
- Abstract
- Introduction
- Motivation
- Theoretical background
- Slow-motional EPR and DNP spectra for an anisotropic g tensor
- Disentangling the solid and Overhauser DNP effects
- Conclusion
- Appendix A: Additional figures
- Code and data availability
- Author contributions
- Competing interests
- Disclaimer
- Acknowledgements
- Financial support
- Review statement
- References