Gevin
von Witte
ab,
Aaron
Himmler
b,
Sebastian
Kozerke
a and
Matthias
Ernst
*b
aInstitute for Biomedical Engineering, University and ETH Zurich, 8092 Zurich, Switzerland
bDepartment of Chemistry and Applied Biosciences, ETH Zurich, 8093 Zurich, Switzerland. E-mail: maer@ethz.ch
First published on 7th February 2024
Dynamic nuclear polarization enables the hyperpolarization of nuclear spins beyond the thermal-equilibrium Boltzmann distribution. However, it is often unclear why the experimentally measured hyperpolarization is below the theoretically achievable maximum polarization. We report a (near-) resonant relaxation enhancement by microwave (MW) irradiation, leading to a significant increase in the nuclear polarization decay compared to measurements without MW irradiation. For example, the increased nuclear relaxation limits the achievable polarization levels to around 35% instead of hypothetical 60%, measured in the DNP material TEMPO in 1H glassy matrices at 3.3 K and 7 T. Applying rate-equation models to published build-up and decay data indicates that such relaxation enhancement is a common issue in many samples when using different radicals at low sample temperatures and high Boltzmann polarizations of the electrons. Accordingly, quantification and a better understanding of the relaxation processes under MW irradiation might help to design samples and processes towards achieving higher nuclear hyperpolarization levels.
Dynamic nuclear polarization (DNP) is a hyperpolarization method that is based on the transfer of polarization (spin magnetic moment) from unbound valence electrons (radicals, defects, paramagnetic centers) to hyperfine-coupled nuclear spins under microwave (MW) irradiation. DNP experiments at cryogenic temperatures offer the advantages of long electron relaxation time and increased Boltzmann polarization of the electrons to near unity. Typically, the electron spins are very dilute in large nuclear spins systems with on average hundreds of nuclear spins per electron. Such a configuration leads to a slow build up of hyperpolarization as successively more and more nuclear spins are polarized either directly from the electrons or indirectly through nuclear spin diffusion. Eventually a steady-state polarization (P0) is reached where nuclear relaxation and polarization injection compensate each other. The build-up of the polarization can often be described by a characteristic time constant τbup. Upon switching off the MW irradiation, the hyperpolarized state relaxes back to the thermal equilibrium, commonly called polarization decay and described by a characteristic time constant τdecay. Rate equation or kinetic models offer the possibility to qualitatively and even quantitatively describe the balance between hyperpolarization injection and relaxation.17
In this paper we use such rate equation models to determine the relaxation-rate constant of the nuclei during the DNP polarization build-up process. We find that in many cases the nuclear relaxation-rate constants during the hyperpolarization build-up under MW irradiation are much larger than the ones obtained without MW irradiation. Microwave irradiation was previously shown to shorten rotating-frame relaxation times (T1ρ) in cross-polarization experiments18 and also nuclear coherence life times (T′2)19 at liquid helium temperatures and 6.7 T field strength. In this work, the increased nuclear relaxation-rate constants under microwave irradiation during the build up of hyperpolarization is found to limit the achievable maximum polarization compared to the prediction when using the electron equilibrium polarization, nuclear relaxation-rate constants (without MW irradiation) and polarization injection rates.
Electron paramagnetic resonance (EPR) spectra were acquired with our in-house developed longitudinal-detection (LOD) EPR set-up at 7 T at a temperature of 5 K.22 A home-made coil was inserted along the main magnetic field (the NMR coil is perpendicular to it) with the MW source irradiating along the NMR coil's direction (same MW wave guide as for NMR). MW irradiation leads to a saturation of the electron polarization which is inductively detected with the longitudinal coil. Upon switching off the MW, the electron polarization returns to its thermal equilibrium value, again creating a detectable voltage. The LOD EPR and NMR coils cannot be mounted at the same time but only require minor changes to the cryostat inset with the MW set-up and cooling unaffected.
For all build-up curves (except where otherwise stated to check the off-resonant irradiation), the MW frequency was set to 197.05 GHz without using any MW frequency modulation. The MW power was set to the maximum output unless explicitly stated. For decay curves under MW irradiation, the MW frequency was switched to 197.28 GHz (zero crossing of the DNP profile) after the sample (nearly) reached its steady-state polarization using maximum output power unless otherwise stated.
All data processing was performed with in-house developed MATLAB (MathWorks Inc., Natick, Massachusetts, U.S.A) scripts. All experimental uncertainties in the processed experimental data result from the 95% fit intervals. Experimental instabilities like slight changes in the MW output or minor temperature fluctuations as well as uncertainties in the thermal equilibrium measurement were not included in the analysis.
![]() | ||
Fig. 1 Nuclear relaxation depending on the temperature and MW irradiation. Build-ups (a) and the respective decays (b) of the (second) natural abundance sample at 3.3 K (red and orange) and 2.0 K (light and darker blue) without correcting for the effects of radio-frequency (RF) pulses. The MW is set to the same power and frequency (197.05 GHz) during all build-ups. During the decays depicted in darker blue and red, the MW is switched off. For the decays depcited in orange (3.3 K) and light blue (2.0 K), the MW is switched to the central EPR frequency (197.28 GHz) resulting in zero DNP, as can be seen in (d). The measured polarization in (a) can be simulated (solid black line) with a single compartment rate equation model.17 The only model input parameters are the thermal electron polarization, measured build-up time and steady-state polarization. (c) Relaxation rates during MW-on decays (dots) and MW-off decays (diamonds) for the natural abundance samples (green) and partially deuterated (DNP juice) sample (brown-gray/taupe, see Methods for details). The data shown are without RF correction (explained in the main text) (d) Overlay of the longitudinal-detected EPR line at 5 K (7 T, blue) from22 and a DNP profile. For simplicity, all measurements presented in this work were recorded with continuous wave (CW) irradiation without MW modulation. |
The corresponding decays shown in Fig. 1b were recorded with the MW source switched off (dark blue and red curves) as well as with the MW source turned on (light blue and yellow curves) but tuned to the frequency of the zero crossing of the DNP profile (197.28 GHz, see Fig. 1d), minimizing the DNP enhancement during the decay measurements. There are several effects visible when comparing the decay curves with and without MW irradiation:
(i) Without MW irradiation, the decay of the nuclear polarization is slower (dark blue and red line). Switching the MW irradiation to the frequency of the zero crossing of the DNP profile significantly speeds up the nuclear spin–lattice (T1) relaxation (light blue and yellow line in Fig. 1b).
(ii) The decay curves with on-resonance MW irradiation (light blue and yellow lines in Fig. 1b) start at a higher value, the ones without MW irradiation (dark blue and red lines in Fig. 1b) at a lower value than the end point of the build-up curves. This can be explained by the presence of strongly hyperfine-coupled and, therefore, invisible (hidden or quenched) nuclear spins. Microwave irradiation will decouple some of these spins and make them visible and contribute to the observable nuclear signal intensity.23,24 Irradiation closer to the center of the EPR line would result in a stronger decoupling effect as more electrons are affected by the decoupling.
(iii) In the low temperature (2.0 K), MW-off decays corrected for the perturbations by RF pulses shown in Fig. S1 (dark blue) (ESI†), the observed signal first increases slightly, in contrast to the expected decay as observed in the uncorrected data (Fig. 1b and Fig. S2, ESI†). This can again be explained by invisible strongly hyperfine-shifted spins that are coupled to the visible spins by slow spin diffusion.25
In the following, we will not discuss the latter two points but focus on the first point. Such a dependence of nuclear decay rates on the MW irradiation at low temperatures has not yet been reported to the best of our knowledge. Similar and related effects, i.e., a decrease of rotating-frame relaxation times (T1ρ) in cross-polarization experiments under MW irradiation18 and a decrease of transverse nuclear coherence life times (T′2) in echo experiments19 with MW irradiation at liquid-He temperatures have been reported recently.
Fig. 1c shows the fitted decay-rate constants for three different samples with and without MW irradiation at three different temperatures. In all cases, a significant increase of the nuclear decay-rate constant is observed under MW irradiation. Under MW irradiation, the relaxation is nearly independent of the temperature between 2.0 and 3.9 K, contrasting the nearly one order of magnitude difference in decay-rate constants without MW irradiation over this range of temperatures.
To investigate this further, we measured the decay-rate constants as a function of the MW settings (power and frequency) which is summarized in Fig. 2a and b. The relaxation enhancement shows a frequency-dependent resonant behavior with stronger relaxation enhancement close to the zero crossing of the DNP profile which corresponds to the center of the EPR spectrum (Fig. 1d). For a DNP frequency of 197.28 GHz (zero crossing of the DNP profile), the power dependence of the relaxation enhancement can be described with a square-root scaling (Fig. 2b) and scales in the same way as the build-up time constant (see Fig. 2c) as a function of MW power. The steady-state polarization P0, however, shows a different behavior as shown in Fig. 2d: It is nearly independent of the power between 20 and 100% of the MW output.
To understand the almost MW-power independent steady-state polarization while the decay-rate as well as build-up time constant τbup follow a square root scaling, we use a single compartment model.17 For a single homogeneous compartment, the hyperpolarization build-up can be described using a first-order differential equation with a hyperpolarization injection rate constant kW and a relaxation-rate constant of the build-up kbupR
![]() | (1) |
τbup−1 = kW + kbupR | (2a) |
![]() | (2b) |
For the decay, kW would be set to zero (MW off), leading to τdecay−1 = τR−1 = kdecayR. The model can reproduce the experimental build-up curve as shown in Fig. 1a (black line). A different approach to measure the rate-equation parameters is discussed in Section S2 of the ESI.†
From the definition of the build-up time in eqn (2a) and the observed square root scaling of τdecay and τbup with MW power, the injection rate constant, kW will obey a power scaling similar to a square root. Since all rates in the steady-state polarization (cf.eqn (2b)) follow roughly a square root scaling, a power-independent steady-state polarization is expected. Therefore, higher MW power will lead to a faster build-up but not to a higher steady-state polarization.
However, this does not hold true for low MW powers as in this case the relaxation enhancement under MW irradiation is small compared to the thermal relaxation and the total relaxation-rate constant will not scale with the square root of the MW power. Hence, the total relaxation in the system competing with the DNP injection can be described as kR = kdecayR + kRMWsB1,MW with the MW-field amplitude (MW power PMW) and a scaling factor s. For small MW powers, the total relaxation in the system is dominated by thermal relaxation, explaining the increase of the steady-state polarization with MW power as only the injection-rate constant grows with the square root of the MW power (linear in B1,MW). For larger MW powers, MW induced relaxation governs the total relaxation and no further improvements of the steady-state polarization and only a faster build-up is expected (see Fig. 2). A shortening of the build-up time with a preserved enhancement might help to improve the NMR signal per unit time.26
Extracting the rate-equation parameters kR and kW directly from the experimental build-up time constant and steady-state polarization value using eqn (2) leads to relaxation-rate constants three to four times larger than the decay-rate constant measured without MW irradiation (compare Fig. 1c and Fig. S3c, ESI†). Using the injection-rate constant from the highest applied MW power and assuming the decay-rate constant measured without MW irradiation, the steady-state polarization would be around 60–70% instead of the experimentally measured 30–40% as shown in Fig. S4 (ESI†).
The injection rate as used in eqn (1) refers to the creation of polarization (normalized magnetization) and not to absolute magnetization which would describe the absolute number of effective electron-nuclear polarization transfers. In our experiments, all samples contained the same concentration of radicals (50 mM). The approximately 90% deuteration reduced the number of protons per radical. Hence, the similar (polarization) injection rates for deuterated and natural abundance samples (cf. Fig. S3c, ESI†) indicate that DNP (at least in the partially deuterated sample) is limited by transport from the strongly hyperfine coupled spins into the bulk rather than by the initial electron-to-nucleus polarization transfer: Similar steady-state polarizations are found for natural abundance and partially deuterated samples (cf. Fig. S4, ESI†) which requires around ten times more nuclear spins to be polarized per TEMPO radical in the the natural abundance sample. If the electron-to-nuclear spin angular momentum transfer (converting electron polarization/electron Zeeman energy into nuclear spin magnetic moment/nuclear Zeeman energy) would be the DNP limiting step, a ten times higher (polarization) injection rate would be expected for the deuterated sample.
Similar conclusions have been made for partially deuterated samples,27 termed “spin-diffusion limited” samples, with a somewhat more complicated model. For high 1H concentrations, ref. 27 ascribes the polarization to be “spin-exchange limited” as too many nuclei per radical need to be polarized, possibly supported by the slight decrease in kW for the natural abundance samples compared to the partially deuterated sample (cf. Fig. S3c, ESI†).
However, in the current work it remains unclear why the partially deuterated sample has higher relaxation rate constants during build-up and decay (cf.Fig. 1c, 2a, b and Fig. S4c, ESI†) and similar injection rate constant (Fig. S4c, ESI†) compared to the natural abundance samples. This higher relaxation results in a lower steady-state polarization.
To investigate whether an enhancement of the apparent decay rates under MW irradiation is a particular property of our experimental sample composition, we have reevaluated data available in the literature using the one-compartment model. Using experimental build-up and decay times as well as steady-state polarization levels allows us to calculate and compare the relaxation-rate constants during build up and decay.
c(Gd) | B 0 | T | c(e−) | P 0 | τ bup | T 1 | k W | k bupR | k decayR | k bupR/kdecayR | Ref. | ||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
mM | T | K | mM | % | s | s | (103 s)−1 | (103 s)−1 | (103 s)−1 | ||||
a The ratio kbupR/kdecayR = 1 indicates the absence of relaxation enhancement by MW irradiation. b The sample formulation was broken down to the radical together with its target spins and possible gadolinium (Gd) doping, although the exact sample formulations vary between the different works. The first three samples from this work refer to the natural abundance sample and the second three to the partially deuterated (DNP juice) sample (cf.Fig. 1 and 2). RF correction was applied by the authors to the decay of ref. 29 but not to any of the other reported measurements according to our knowledge. The TEMPO data set from ref. 21 (same sample formulation as for our natural abundance samples apart from the [13C]urea) was recorded in our lab previously although with an older MW set-up with an approximately eight times lower MW power at the sample (see Methods), possibly explaining the discrepancy to the measurements reported herein. | |||||||||||||
Trityl | 13C | 1.8 | 3.35 | 3.5 | 15 | 3.5 | 353 | 419 | 0.17 | 2.66 | 2.39 | 1.1 | 28 |
Trityl | 13C | 3.35 | 1.0 | 15 | 27 | 1200 | 10![]() |
0.23 | 0.60 | 0.10 | 6.0 | 29 | |
Trityl | 13C | 3.35 | 1.0 | 14.9 | 50 | 1500 | 9000 | 0.34 | 0.33 | 0.11 | 2.9 | 29 | |
Trityl | 13C | 4.64 | 1.0 | 14.3 | 64 | 3000 | 20![]() |
0.21 | 0.12 | 0.05 | 2.4 | 29 | |
Trityl | 13C | 4.64 | 1.0 | 18.5 | 58 | 1600 | 22![]() |
0.36 | 0.26 | 0.05 | 5.7 | 29 | |
Trityl | 13C | 1.5 | 4.64 | 1.0 | 45.4 | 23 | 475 | 9000 | 0.49 | 1.62 | 0.11 | 15 | 29 |
Trityl | 13C | 4.64 | 1.0 | 14.3 | 70 | 3000 | 14![]() |
0.23 | 0.10 | 0.07 | 1.4 | 29 | |
Trityl | 13C | 3.35 | 1.3 | 15 | 27 | 1000 | 6400 | 0.29 | 0.71 | 0.16 | 4.6 | 30,31 | |
Trityl | 13C | 6.7 | 1.3 | 15 | 65 | 3700 | 33![]() |
0.18 | 0.09 | 0.03 | 3.1 | 30,31 | |
Trityl | 13C | 6.7 | 1.3 | 30 | 70 | 1200 | 33![]() |
0.58 | 0.25 | 0.03 | 8.2 | 30,31 | |
Trityl | 13C | 10.1 | 1.3 | 15 | 70 | 26![]() |
90![]() |
0.03 | 0.01 | 0.01 | 1.0 | 30,31 | |
Trityl | 13C | 10.1 | 1.3 | 45 | 70 | 2200 | 90![]() |
0.32 | 0.14 | 0.01 | 12 | 30,31 | |
Trityl | 1H | 7 | 20 | 10 | 0.0 | 80 | 85 | 0.01 | 12.5 | 11.8 | 1.1 | 32 | |
Trityl | 1H | 7 | 20 | 40 | 0.3 | 40 | 50 | 0.29 | 24.7 | 20.0 | 1.2 | 32 | |
Trityl | 1H | 7 | 20 | 100 | 1.0 | 1.75 | 2.5 | 24.5 | 547 | 400 | 1.4 | 32 | |
TEMPO | 13C | 3.35 | 4.2 | 25 | 1.0 | 346 | 420 | 0.06 | 2.83 | 2.38 | 1.2 | 21 | |
TEMPO | 13C | 3.35 | 4.2 | 50 | 1.4 | 117 | 169 | 0.24 | 8.30 | 5.92 | 1.4 | 21 | |
TEMPO | 13C | 3.35 | 4.2 | 75 | 1.3 | 57 | 91 | 0.45 | 17.1 | 11.0 | 1.6 | 21 | |
TEMPO | 1H | 3.35 | 4.2 | 25 | 3.7 | 102 | 192 | 0.75 | 9.05 | 5.21 | 1.7 | 21 | |
TEMPO | 1H | 3.35 | 4.2 | 50 | 4.7 | 34 | 84 | 2.84 | 26.6 | 11.9 | 2.2 | 21 | |
TEMPO | 1H | 3.35 | 4.2 | 75 | 4.5 | 18 | 42 | 5.08 | 50.5 | 23.8 | 2.1 | 21 | |
TEMPO | 13C | 7 | 4.2 | 25 | 3.6 | 1271 | 1850 | 0.04 | 0.75 | 0.54 | 1.4 | 21 | |
TEMPO | 13C | 7 | 4.2 | 50 | 4.4 | 414 | 666 | 0.13 | 2.28 | 1.50 | 1.5 | 21 | |
TEMPO | 13C | 7 | 4.2 | 75 | 4.2 | 205 | 349 | 0.26 | 4.62 | 2.87 | 1.6 | 21 | |
TEMPO | 1H | 7 | 4.2 | 25 | 7.7 | 242 | 439 | 0.39 | 3.74 | 2.28 | 1.6 | 21 | |
TEMPO | 1H | 7 | 4.2 | 50 | 18.4 | 81 | 216 | 2.81 | 9.54 | 4.63 | 2.1 | 21 | |
TEMPO | 1H | 7 | 4.2 | 75 | 16.7 | 39 | 106 | 5.30 | 20.3 | 9.43 | 2.2 | 21 | |
TEMPO | 1H | 7 | 3.9 | 50 | 27 | 40 | 197 | 7.99 | 17.0 | 5.08 | 3.4 | This work | |
TEMPO | 1H | 7 | 3.3 | 50 | 37 | 54 | 340 | 7.63 | 10.9 | 2.94 | 3.7 | This work | |
TEMPO | 1H | 7 | 2.0 | 50 | 40 | 58 | 1200 | 7.06 | 10.2 | 0.83 | 12 | This work | |
TEMPO | 1H | 7 | 3.9 | 50 | 24 | 27 | 93 | 10.8 | 26.2 | 10.8 | 2.4 | This work | |
TEMPO | 1H | 7 | 3.3 | 50 | 25 | 29 | 149 | 9.75 | 24.7 | 6.71 | 3.7 | This work | |
TEMPO | 1H | 7 | 2.0 | 50 | 30 | 32 | 870 | 9.61 | 21.6 | 1.15 | 19 | This work |
The relaxation-rate constant extracted from the build-up kbupR (compare eqn (2)) and decay kdecayR are compared to estimate the relaxation enhancement by MW irradiation with a ratio of kbupR/kdecayR = 1 indicating no relaxation enhancement in the sample. The estimated ratios for the relaxation enhancements are sensitive to inaccuracies in the thermal equilibrium measurement as this affects the steady-state polarization quoted. The perturbing effects of monitoring RF pulses can be minimized by using few and small flip-angle pulses or correcting for the pulses.17
In the following, relaxation enhancement by MW irradiation depending on the experimental conditions (temperature and field) as well as sample formulation (radical, concentration and target spins) are discussed. Trends are first discussed for trityl-based and then for TEMPO-based samples.
For trityl radicals and direct DNP to 13C nuclei, the relaxation enhancement is more pronounced at lower temperatures and lower fields. At 3.35 T, the relaxation enhancement is nearly absent at 3.5 K, grows to more than fourfold at 1.3 K and sixfold at 1.0 K. Increasing the field strength reduces the relaxation enhancement in the range of 3.35 and 10.1 T and temperatures of 1.0–1.3 K. Higher trityl (electron) concentrations lead to pronounced increases of the relaxation enhancement. For example, at 10.1 T and 1.3 K, relaxation enhancement is absent for a 15 mM sample while a twelvefold relaxation enhancement is found for 45 mM of trityl – resulting in twelvefold shortening of τbup while preserving steady-state polarization P0. The relaxation enhancement at higher electron concentrations can be reduced by sample doping with gadolinium (Gd).
Raising the temperature to 20 K (and DNP to the 1H nuclei), the relaxation enhancement is nearly completely suppressed or masked by increased thermal relaxation even at very high trityl concentrations of 100 mM.
For TEMPO with DNP performed to 1H or 13C and 4.2 K, no field dependence is evident comparing measurements at 3.35 and 7 T. The elevated temperature complicates discrimination between temperature-related and radical-related (DNP mechanism) changes. The relaxation enhancement appears slightly larger for 1H compared to 13C although it is not increasing beyond a factor of 2.2.
Two general trends may be deduced from Table 1:
(i) The relaxation enhancement by MW irradiation requires elevated electron concentrations and electron dipolar couplings causing electron spectral diffusion (eSD). Thus, suppression of eSD should reduce the relaxation enhancement. This is in agreement with the results of Gd doping, shortening the electronic relaxation times,33,34 or higher magnetic fields, broadening the EPR line and shortening the electronic relaxation.34
(ii) Lower temperatures increase the relaxation enhancement, prolong the electronic relaxation times and increase the electron polarization. Together, these determine the nuclear relaxation times near liquid helium temperature as all nuclear relaxation can be assumed to originate from paramagnetic relaxation. With lower thermal relaxation, the relaxation enhancement by MW irradiation contributes increasingly more to the total relaxation.
In the following we will discuss possible explanations for the relaxation enhancement by MW irradiation. MW irradiation saturates parts of the electron line, creating a polarization difference between different parts of the electron line. At liquid-He temperatures, nuclear relaxation is most likely due to electron flip-flops.35 Cross effect and thermal mixing DNP convert this polarization difference into nuclear hyperpolarization. Assuming that the non-irradiated part of the electron line remains at the thermal electron polarization (A = Pethermal) while the irradiated part of the electron line is saturated to PeMW, nuclear relaxation for electrons with different polarization is proportional to 1 − PethermalPeMW35,36 (a similar expression is derived for direct nuclear (paramagnetic) relaxation in ref. 35). Hence, nuclear relaxation and DNP injection scale with the electron saturation, leading to a MW power independent steady-state polarization (for small enough thermal relaxation, cf.Fig. 2d). For MW irradiation at the center of the electron line, the higher electron spectral density and relative ease to fulfill the matching condition for triple spin flips results in a frequency dependence with the higher nuclear relaxation at the center of the electron line. Increased electron spectral diffusion increases the probability for electron flip-flops which cause nuclear relaxation.
For (water) ice the frequency dependence of electromagnetic absorption has been reported in the literature37 and we can assume it to be frequency independent over the relatively narrow experimental frequency range (192–198 GHz). Since the set-up does not involve any type of cavity, no additional frequency dependence except for the weak power dependence of the MW source is expected. Thus, a frequency-independent heating would occur if the local cooling is weaker than the absorbed MW power. A rough estimate for this is given in Section S5 of the ESI.† If a sample heating would occur, a higher sample temperature reduces the thermal electron polarization which increases the relaxation scaling as 1 − (Pethermal)2. This might explain the slightly decreasing steady-state polarization for the highest MW powers (see Fig. 2d) and enhanced relaxation for far off-resonant MW irradiation (see Fig. 2a) compared to the MW off situation. For elevated temperature, the thermal relaxation is following a 1 − (Pethermal)2 scaling and is large such that a relaxation enhancement by MW irradiation is difficult to observe.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3cp06025j |
This journal is © the Owner Societies 2024 |