The solid effect of dynamic nuclear polarization in liquids – accounting for g-tensor anisotropy at high magnetic fields

Abstract 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 T2 . 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 , 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.

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). 2uring 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 slowmotional 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.

Motivation
DNP aims to increase the longitudinal nuclear magnetization, i z , beyond its equilibrium Boltzmann value, i eq z .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 ss z 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 ss z , from its equilibrium value, s eq z .
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 0 1I ) of the polarizing agent.In DNP, T 1I is typically (much) shorter than T 0 1I 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 = 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 1I = 1/T 1I , which is quantified by v + , is generally negligibly small.As a result, v + /R 1I 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 1I 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., 2022Kuzhelev et al., , 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 D. Sezer et al.: Solid effect in liquids 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.

Theoretical background
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.

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, and3 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, D. Sezer et al.: Solid effect in liquids Solving these algebraic equations for the variables s ss i , 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 ss x and s ss z in terms of s ss y : 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 abs = s ss y /s 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 where and B 0 = B(ω I = 0).The iω I that has been added to the main diagonal of the matrix B will be needed for the dynamical description of the solid effect (see Sect. 3.2).The subscript of B 0 is intended as a reminder that B is evaluated at ω I = 0, where ω I is the Larmor frequency of the polarized nuclear spin.

Solid effect in solids
SE relies on the dipolar interaction between the electronic and nuclear spins whose coupling is Here D dip = (µ 0 /4π ) γ S γ I is the dipolar constant, which equals approximately 2π (79 kHz nm 3 ) for protons, and (r, θ, φ) 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 electronnucleus 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 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 ss z 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 A 1I 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 ss z .Inverting the matrix B in Eq. ( 22) and using s ss x,y from before, we find where 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 B −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 B by expressing g ss x and g ss z in terms of g ss y , analogously to our treatment of the Bloch equations in the previous subsection.
From the upper and lower rows of B in Eq. ( 22), we find Substituting this g ss z 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 ss x and g ss z 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 0 = P (ω I = 0).Finally, using this g ss y 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

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 B −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 B, i.e., BU = U , where the diagonal matrix = diag(λ 1 , λ 2 , λ 3 ) contains the three eigenvalues and the columns of U 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 L = diag( 1 , 2 , 3 ) and calculate Q = ULU −1 .Comparing this expression of Q with B −1 = U U −1 , we see that, in the transition from solids to liquids, where B −1 is replaced with Q, we essentially "process" the eigenvalues of B with the spectral density function j 11 .This step prevents us from eliminating the variables g ss x,z 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 A 1I , since when ω 1 = 0, the zz component of B 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 A 1I (cf.Eq. 29).

Slow-motional EPR and DNP spectra for an anisotropic g tensor
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).

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 ∇ 2 acts on the orientation variable .The operator satisfies the following eigenvalue problem: where the eigenfunctions D mn ( ) are the Wigner rotation matrix elements, which are orthogonal to each other: From Eq. ( 46) it is clear that the time derivative on the lefthand 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: 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 G 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 − → s ( , t), whose spatial part evolves according to the Bloch equations (Eq.51) and whose dependence evolves according to the diffusion equation (Eq.44): Here E 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 D mn ( ) form a complete set, we expand − → s ( , t) as follows: The coefficients − → s mn , which contain the time dependence, can be obtained from − → s ( , t) using the orthogonality of D mn ( ) (Eq. 47): Ultimately, the only property that we care about is the integral of the SLE variable − → s ( , t) over all the orientations: In that sense, the (vector) coefficient − → s 0 00 (t) is the main object of interest, while all the other coefficients − → s mn (t) play an auxiliary, bookkeeping role.Substituting − → s ( , t) from Eq. ( 54) into Eq.( 53), multiplying both sides by D L * MN ( ) and integrating over , we get Clearly, the terms proportional to K and B 0 in Eq. ( 53) do not mix coefficients − → s L MN with different values of L, M and N .In other words, these two operators are diagonal in the selected representation.The term proportional to F( ), 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 LM 1 m 1 2 m 2 .Specifically, for the D 2 K0 ( ) in Eq. ( 52), we have In Eq. ( 59), the sum over mixes only expansion coefficients with = L, L±2 (Freed et al., 1971) because all three Wigner rotation matrix elements in F have L = 2 (Eq.52).Since we need − → s 0 00 at the end, it is sufficient to consider only coefficients with even values of .Furthermore, as the Wigner rotation matrix elements in F have M = 0 ± 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 − → s m0 with even and m.Finally, because the Wigner rotation matrix elements with M = 2 and M = −2 appear in a symmetrical way in F, it becomes possible to work with the symmetrized coefficients (Freed et al., 1971) thus restricting M to non-negative values (0 ≤ M ≤ L).The lowest-order coefficients that are coupled by the SLE dynamics are thus − → s 00 , − → s 20 , − → s 22 , − → s 40 , − → s 42 , − → s 44 , − → s 60 , etc.

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 For the smallest non-trivial choice of L max = 2, n tot = 3 (with − → s 00 , − → s 20 and − → s 22 ).The number of coefficients increases quadratically with L max (e.g., n tot = 15, 28, 45 for L max = 8, 12, 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 ( − → s 00 , − → s 20 and − → s 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 γ 2 0 and γ 2 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 − → s L0 .
Selecting L max = 0 in the above formalism amounts to retaining only the (3D vector) coefficient − → s 00 .Then the matrix B 0 in Eq. ( 63) reduces to B 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 LM z and s LM x,y , 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 − → 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 B 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 gtensor anisotropy on the spin dynamics will be helpful when adapting the approach to SE-DNP in Sect.4.4 and 4.5.

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 B 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 B 0 (Eq.63), we express s ss x and s ss z in terms of s ss y : 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 s ss z , we thus have s 00 z = ω 1 T 1 s 00 y + s eq z , 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 B 0 (Eq.63) and after substituting s ss x and s ss z from Eq. ( 70), we find where we have introduced the n tot × n tot matrix and P 0 = P(ω I = 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 s ss y and s ss x , 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 B 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.

Solid effect in "solids"
Extending the above treatment to SE-DNP, we combine the spin dynamics in Eq. ( 35  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 − → g LM : The steady state of the resulting spin dynamics is then where generalizes the matrix B 0 from Eq. ( 63).
Our goal is to solve for g ss z , since its 00th component should be used in Eq. ( 23) to calculate the rate constants pv − and v + .After inverting B in Eq. ( 83), we find Note that now [B −1 ] ij denotes the n tot × n tot sub-matrix of B −1 at position ij and not a scalar matrix element.Using s ss x,y from the previous subsection, we find that the first component of g ss z is Substituting this result into Eq.( 23), we obtain These expressions, which require the inversion of the 3n tot × 3n tot matrix B, 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 g ss x and g ss z in terms of g ss y using the first and last rows of B: Substituting the first (i.e., 00th) component of g ss z into Eq.( 23), we find Because R 1 is a diagonal matrix, this result is identical to Eq. ( 29), showing that R A 1I is not affected by the anisotropy of the g tensor in the case of isotropic rotation.
Similarly, from the middle part of B, we obtain We first observe that (R 1 + iω I ) −1 1 00 = 1 00 (R 1 + iω I ) −1 because R 1 is diagonal.Then we substitute s ss x,y from before to get Finally, substituting the 00th element of g ss y 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 + /δ 2 and pv − /δ 2 , which have units of time.Although the electronic nonsaturation 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 − /δ 2 in the first two rows of Fig. 3.Note that v + /δ 2 and pv − /δ 2 were calculated directly from Eq. ( 92), whereas v − /δ 2 was determined by dividing pv − /δ 2 by p = 1 − ω 2 1 T 1 [P −1 0 ] 11 (Eq.75).The columns in Fig. 3 reveal the effect of the g-tensor anisotropy on the different factors relevant to SE. v + /δ 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 + /δ 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 − /δ 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 + /δ 2 and −v + /δ 2 (dashed red lines).We see that, for all the tumbling rates, the two SE lines comprising v − /δ 2 exactly match their counterparts in v + /δ 2 .
The first row of Fig. 3 shows the electronic saturation under g-tensor anisotropy (we actually plot the "nonsaturation" p = 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 https://doi.org/10.5194/mr-4-243-2023Magn.Reson., 4, 243-269, 2023 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 − /δ 2 and v + /δ 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 pv − /δ 2 (solid green line), which equals the product of the first and second rows.From Eq. ( 7), we know that pv − /δ 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 pv − /δ 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.

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 B −1 in Eq. ( 87): Because the zz sub-block of B 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 A 1I 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 + / δ 2 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 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 gtensor 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 + / δ 2 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 − / δ 2 (orange) are now completely different from their counterparts in v + / δ 2 (dashed red).The increased intensity in the vicinity of the electronic resonance due to motional broadening is also manifested by v − / δ 2 .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 mag-nitude 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  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 pv − / δ 2 (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.

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 γ 2 0 and γ 2 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 highfrequency offsets.This points to the possibility that the inphase 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: γ 2 0 , γ 2 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 = diag(14, 14, 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 φ = −1.3• , 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.

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 pv − ( )/ δ 2 (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  1.
Table 1.Parameters obtained from the fits to the experimental data.B 1 = 0.02 G for EPR and 5.5 G for DNP.T 2 = 20 ns was used for both EPR and DNP.  (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.

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 φ = −2 • .
Although, overall, the fit is not bad, the middle part of the calculated spectrum changes too sharply, and its highfrequency 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 tenhttps://doi.org/10.5194/mr-4-243-2023Magn.Reson., 4, 243-269, 2023 sor 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.Neverthe-less, 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.

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 − / δ 2 , 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 0 1I 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   1. Observe that the OE contribution to the DNP spectrum (dotted-dashed blue line) is narrower (i.e., more "liquid-like") than that of 10-Doxyl-PC (Fig. 5b), while the SE contribution is more "solid-like" because the timescale τ ffhs is 2.4 times longer.
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 righthand side of Eq. ( 100), are about 2 ns −1 for 10, 16 and 16* (Table 2).Since v + / δ 2 is about 1 ps (Fig. 4, third row), we conclude that v + T 1I 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 (  ) and without (T 0 1I ) spin labels were measured at two different temperatures.These determine the leakage factor f .The coupling factor c is obtained from f and σ OE using Eq. ( 99).The magnitude of the dipolar interaction responsible for SE ( δ 2 ), obtained from σ SE using Eq. ( 100), provides information about the effective contact distance (b).Combining b with τ ffhs from (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 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.

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 correla-tion 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(Livshits et al., , 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 θ 0 = 40 • 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 θ 0 = 75 • , 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 θ 0 = 75 • 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 θ 0 = 90 • .)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 dynam-ics 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. 5Once 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 B 0 (Sect.3.1) and the "new Bloch equations" (Eq.21) with the matrix B = B 0 +iω 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 B 0 became matrices in the space of the angular-momentum indices LM.The result was the "expanded Bloch matrix" B 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.

Conclusion
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.

Figure 1 .
Figure 1.Experimental cw-EPR spectra (a, b) and DNP spectra (c, d) of spin-labeled lipids in DOPC lipid bilayers at 9.4 T and ≈ 320 K.The nitroxide spin label (Doxyl) is either at position 10 (a, c) or at position 16 (b, d) of the aliphatic lipid chain.The integrated cw-EPR spectra (dashed-dotted blue lines in panels a and b) are used to decompose the DNP spectra (c and d) into contributions from OE (dashed-dotted blue lines) and SE (dashed green lines).The sum of these two contributions is indicated with a dotted black line.

Figure 2 .
Figure 2. Cw-EPR spectra for g = diag(2.00755,2.00555, 2.0023) and different tumbling times τ rot .A larger L max is necessary for slower tumbling.The needed L max also depends on the anisotropies, which are (γ 0 , γ 2 ) = (−373, 107) MHz for B 0 = 9.403 T. Other simulation parameters were B 1 = 0.02 G, T 1 = 100 ns and T 2 = 20 ns.Note that these are the values of T 1 and T 2 that are introduced by hand.The rotational tumbling itself further shortens the actual spin-spin relaxation time T 2 .

Figure 3 .
Figure 3. Solid-effect rates v + (red line) and pv − (green line) calculated at high mw power (B 1 = 5.5 G) for static dipolar interaction (i.e., "solid") and several different rates of rotational tumbling.The factorization of pv − into p (blue line) and v − (orange line) is also shown.All other parameters are as in Fig. 2, and L max = 10.In particular, T 1 = 100 ns.

Figure 4 .
Figure 4. Same as Fig. 3 but for the FFHS model of translational diffusion with τ ffhs = 6 ns.

Figure 5 .
Figure 5. Fits to the experimental cw-EPR spectrum (a) and DNP spectrum (b) of 10-Doxyl-PC.In both cases, our best fits are shown with dashed black lines.The DNP spectrum in panel (b) is calculated by adding the contributions of SE (solid green line) and OE (dotted-dashed blue line), both of which are affected by the g-tensor anisotropy.The fit parameters are given inTable 1.

Figure 6 .
Figure 6.Same as Fig. 5 but for 16-Doxyl-PC.The fitted parameters are given in the second-last column of Table1.Because all fitted lines in Figs. 5, 6 and 7 are calculated only at the experimental offsets, the green SE lines are not perfectly smooth.

Figure 7 .
Figure 7. Same as Fig. 6b but also fitting τ ffhs .The fitted parameters are given in the last column of Table1.Observe that the OE contribution to the DNP spectrum (dotted-dashed blue line) is narrower (i.e., more "liquid-like") than that of 10-Doxyl-PC (Fig.5b), while the SE contribution is more "solid-like" because the timescale τ ffhs is 2.4 times longer.

Figure A3 .
Figure A3.Same as Fig. 4 but with (a) τ ffhs = 3 ns, i.e., 2-fold faster translational motion which broadens the SE lines to a larger extent, and (b) τ ffhs = 12 ns, i.e., more solid-like behavior.Observe how the predicted SE-DNP line shape (green line in the last row) responds sensitively to the timescale of the translational motion that is responsible for averaging the dipolar interaction.

Figure A4 .
Figure A4.Effect of the hyperfine tensor on the calculated EPR spectrum of 10-Doxyl-PC.Various derivative cw-EPR spectra (a) and their numerical integrals (b) are compared with each other.Our calculation (solid black line) agrees perfectly with the Easyspin(Stoll and Schweiger, 2006) simulation without a hyperfine tensor (dashed blue line).Including a hyperfine tensor with components (14, 14, 90) MHz in the Easyspin calculation (dotted-dashed orange line) leads to visible changes in the derivative cw-EPR spectrum.However, the difference of the integrated EPR lines with and without a hyperfine tensor in panel (b) should be negligible as far as the simulation of the DNP spectrum is concerned.

Figure A5 .
Figure A5.Same as Fig. A4 but for 16-Doxyl-PC.Because the cw-EPR spectrum is narrower to begin with, the relative contribution of the hyperfine tensor with components (14, 14, 90) MHz is larger than in the case of 10-Doxyl-PC.Considering that the EPR line will experience additional power broadening and motional broadening in DNP, it should still be possible to safely neglect the extra width that the hyperfine tensor brings to the integrated EPR line in panel (b).

Table 2 .
Analysis of the scaling parameters σ OE and σ SE .Nuclear spin-lattice relaxation times with (T 1I

Table 1 ,
we estimate the diffusion constant of the FFHS model (D ffhs ).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. the