Hyperpolarized relaxometry based nuclear T1 noise spectroscopy in diamond

The origins of spin lifetimes in quantum systems is a matter of importance in several areas of quantum information. Spectrally mapping spin relaxation processes provides insight into their origin and motivates methods to mitigate them. In this paper, we map nuclear relaxation in a prototypical system of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{13}{\rm{C}}$$\end{document}13C nuclei in diamond coupled to Nitrogen Vacancy (NV) centers over a wide field range (1 mT-7 T). Nuclear hyperpolarization through optically pumped NV electrons allows signal measurement savings exceeding million-fold over conventional methods. Through a systematic study with varying substitutional electron (P1 center) and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{13}{\rm{C}}$$\end{document}13C concentrations, we identify the operational relaxation channels for the nuclei at different fields as well as the dominant role played by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{13}{\rm{C}}$$\end{document}13C coupling to the interacting P1 electronic spin bath. These results motivate quantum control techniques for dissipation engineering to boost spin lifetimes in diamond, with applications including engineered quantum memories and hyperpolarized \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{13}{\rm{C}}$$\end{document}13C imaging.

T he power of quantum technologies, especially those for information processing and metrology, relies critically on the ability to preserve the fragile quantum states that are harnessed in these applications 1 . Indeed noise serves as an encumbrance to practical implementations, causing both decoherence as well as dissipation of the quantum states 2,3 . Precise spectral characterization of the noise opens the door to strategies by which it can be effectively suppressed 4,5 -case in point being the emergence of dynamical decoupling techniques that preserve quantum coherence by periodic driving 6 . In these cases, quantum control sets up a filter that decouples components of noise except those resonant with the exact filter period 7 , allowing spectral decomposition of the dephasing noise afflicting the system. Experimentally implemented in ion traps 8 , superconducting qubits 9 , and solid-state NMR 10 , this has spurred development of Floquet engineering to enhance T 2 decoherence times by over an order of magnitude in these physical quantum device manifestations [11][12][13] .
Methods that analogously spectrally fingerprint T 1 relaxation processes, on the other hand, are more challenging to implement experimentally. If possible, however, they could reveal the origins of relaxation channels, and foster means to suppress them. Applications to real-world quantum platforms are pressing: relaxation in Josephson junctions and ion-trap qubits, for instance, occur due to often incompletely understood interactions with surface paramagnetic spins 14 . Relaxation studies are also important in the context of coupled quantum systems, such as those built out of electronic and nuclear spins. In the case of diamond Nitrogen Vacancy (NV) center electronic qubits coupled to 13 C nuclei 15 , for instance, a detailed understanding of nuclear relaxation can have important implications for quantum sensing 16 : engineered NV-13 C clusters form building blocks of quantum networks 17 , are the basis for spin gyroscopes 18 , and are harnessed as quantum memories in high-resolution nano-MRI probes 19 . Nuclear T 1 lifetimes are not dominated by phonon interactions, but instead are set by couplings with the intrinsic electronic spin baths themselves-a complex dynamics that is often difficult to probe experimentally. Indeed only a small proportion of 13 C spins can be addressed or readout via the NV centers, as also the direct inductive readout of these spins suffer from extremely weak signals. Moreover, as opposed to T 2 noise spectroscopy carried out in the rotating frame 13 , probing of T 1 processes have to be performed in the laboratory frame. This necessitates the ability to probe relaxation behavior while subjecting samples to widely varying magnetic field strengths.
In this paper, we develop a method of hyperpolarized relaxometry that overcomes these instrumentational and technical challenges. We measure T 1 relaxation rates of 13 C spins in diamond samples relevant for quantum sensing with a high density of NV centers. Our T 1 noise spectroscopy proceeds with high resolution and over four decades of noise spectral frequency, revealing the physical origins of the relaxation processes. While experiments herein are demonstrated for 13 C spins in electronrich diamond, these results are potentially more widely reflective of relaxation processes operational in other systems, including Si: P 20 , wide bandgap materials such as SiC 21,22 , and diamond-based quantum simulator platforms constructed out of 2D materials such as graphene and hBN [23][24][25] . These results are also pertinent for producing and maintaining polarization in hyperpolarized solids, for applications employing hyperpolarized nanoparticles of Si or diamond as MRI tracers 26,27 , and in the relayed optical DNP of liquids mediated through nanodiamonds 28 , since in these applications T 1 relaxation bounds the achievable polarization levels.
Results 13 C Hyperpolarized relaxometry. Key to our technique is the hyperpolarization of 13 C nuclei at room temperature, allowing the rapid and direct measurement of nuclear spin populations via bulk NMR 28 . Dynamic nuclear polarization (DNP) is carried out by optical pumping and polarizing the NV electrons to large values (>10% 31,32 ) and subsequently transferring this to 13 C nuclei (Fig. 1a). This routinely leads to nuclear polarization levels ≳ 0.5%. In a high-field (7 T) NMR detection spectrometer, for instance, the signals are enhanced by factors exceeding ε $ 300-800 times the Boltzmann value 28 , boosting measurement times by 10 5 -10 6 , and resulting in high single-shot detection SNRs. This permits T 1 spectroscopy experiments that would have otherwise been intractable. Hyperpolarization is equally efficiently generated in single crystals as well as randomly oriented diamond powders, and both at natural abundance as well as enriched 13 C concentrations. The hyperpolarized samples are interfaced to a home built field cycler instrument 33 (see Fig. 1d and video in ref. 34 ) that is capable of rapid and high-precision changes in magnetic field over a wide 1 mT-7 T range (extendable in principle from 1 nT-7 T), opening a unique way to peer into the origins of nuclear spin relaxation. Figure 1d, e schematically describe the experiment. Hyperpolarization in the 13 C nuclei is affected by optical pumping at low fields, typically B pol $ 40 mT, followed by rapid transfer to the intermediate field B relax where the spins are allowed to thermalize (see Fig. 1c), and subsequent bulk inductive measurement at 7 T. Experimentally varying B relax allows one to probe field-dependent lifetimes T 1 ðB relax Þ, and through them noise sources perpendicular to B relax and resonant with the nuclear Larmor frequency γ n B relax (Fig. 1b). Here, γ n ¼ 10:7 MHz T À1 is the 13 C gyromagnetic ratio. This allows the spectral decomposition of noise processes that spawn T 1 relaxation. For instance pairs of substitutional nitrogen impurities (P1 centers) undergoing flipflops (Fig. 1c) can apply on the 13 C nuclei a stochastic spinflipping field that constitutes a relaxation process.
Optical excitation for hyperpolarization involves 520-nm irradiation at low power (~80 mW mm −2 ) applied continuously for~40 s. Microwave (MW) sweeps, simultaneously applied across the NV center ESR spectrum, transfer this polarization to the 13 C spins (see Fig. 2a) 28,35 . To prevent excessive heating during optical irradiation, samples are immersed under water, with the powdered diamonds precipitated in the NMR tube. The hyperpolarization is carried out at room temperature, and at the relatively low optical power densities we employ, we do not observe sample heating (temperature < 100 C). DNP occurs in a manner that is completely independent of crystallite orientation. All parts of the underlying NV ESR spectrum produce hyperpolarization, with intensity proportional to the underlying electron density of states. The polarization sign depends solely on the direction of MW sweeps through the NV ESR spectrum (see Fig. 2a inset). Physically, hyperpolarization arises from partly adiabatic traversals of a pair of Landau-Zener (LZ) crossings in the rotating frame that are excited by the swept MWs. For a more detailed exposition of the DNP mechanism, we point the reader to Ref. 36 .
Low-field DNP. Low-field hyperpolarization is hence excited independent of the fields B relax under which relaxation dynamics is to be studied. There is significant acceleration in acquisition time since optical DNP obviates the need to thermalize spins at high fields where T 1 times can be long (for some samples >30 min). Gains averaging time are % ε 2 T 1 ð7TÞ T 1 ðB pol Þ , which in our experimental conditions exceeds five orders of magnitude. In Fig. 2a for instance on a 10% enriched single crystal, we obtain large DNP enhancements ε ¼ 380, and high single shot SNR % 400. It also reflects the inherently high DNP efficiency: every NV center has surrounding it $10 5 nuclear spins, which we polarize to a bulk value (averaged over all 13 C nuclei) of 0.37% employing just 3000 MW sweeps, indicating a transfer efficiency of %12.3% per sweep per fully polarized nuclear spin, assuming full NV polarization. Harnessing this large signal gain allows us to perform relaxometry at a density of field points that are about two orders of magnitude greater than previous efforts [37][38][39] . Such highresolution spectral mapping (for instance 55 field points in Fig. 2) can transparently reveal the underlying processes driving nuclear relaxation. Indeed, in future work, use of small flip-angle pulses might allow one to obtain the entire relaxation curve with a single round of DNP, and thus the ultrafast relaxometry of the nuclei. Our experiments are also aided by technological attributes of the DNP mechanism. DNP is carried out under low fields and laser and MW powers, and allows construction of a compact hyperpolarizer device that can accessorize a field cycling instrument 40 (see 41 for video of hyperpolarizer operation). The wide range (1 mT-7 T) field cycler is constructed over the 7 T detection magnet, and affects rapid magnetic field changes by physically transporting the sample in the axial fringe-field environment of the magnet 33 . This is accomplished by a fast (2 m s −1 ) conveyor belt actuator stage (Parker HMRB08) that shuttles the sample via a carbon fiber rod (see video in Ref. 34 ). The entire sample (field) trajectory can be programmed, allowing implementation of the polarization, relaxation, and detection periods as in Fig. 1e. Transfer times at the maximum travel range were measured to be 648 ± 4 ms (see Supplementary Fig. 2), short in comparison with the T 1n lifetimes we probe. High positional resolution (50 µm) allows access to field steps at high precision (Supplementary Fig. 3 shows full field-position map). The field is primarily in theẑ direction (parallel to the detection magnet), since sample transport occurs centrally, and the diameter of the shuttling rod (8 mm) is small in comparison with the magnet bore (54 mm). This ensures that the spatial gradients of the magnetic field, both longitudinally, as well as transverse, are negligible (field variation under 1% of the bias field) (see Supplementary Note 2).
Relaxation rates. In this paper, we perform T 1 noise spectroscopy on 13 C nuclei in a variety of diamond samples outlined in Table 1. Figure 2 shows representative results, considering here a 10% enriched single crystal (Sample 6). This intriguing data can be visualized in several complementary ways. First, considering relaxation rate R 1 ¼ 1=T 1 (Fig. 2b), the high-resolution data allow us to clearly discern three regimes: a steep narrow R 1  Table 1). a Signal gains due to hyperpolarization under optimal conditions at B pol % 36 mT. Red line shows a single-shot hyperpolarized signal (SNR % 400) after 60 s of optical pumping. Blue line is the 7 T thermal signal after 20 averages, allowing us to quantify signal enhancement from DNP %372 over 7 T, a time saving by %10 6 for equal SNR. Inset: Exemplary signals at B pol % 36 mT under low-to-high (high-to-low) frequency sweeps leading to positive (negative) 13 C hyperpolarization. b Relaxation rate R 1 ¼ 1=T 1 obtained from relaxometry over a wide-field range 20 mT-1.5 T. We observe a rapid growth in relaxation rate below a knee field of 0.5 T, and saturation at higher fields. Inset: Data can be fit to two Tsallian functions, which we ascribe to be originating from inter-13 C couplings and interactions to the P1 spin bath. c Spin lifetimes as a function of field, showing significant boost in nuclear T 1 beyond the knee field, approaching a lifetime % 2.1 min. Inset: Typical relaxation data at two representative fields showing monoexponential character. d Logarithmic scale data visualization, displaying a more equanimous sampling of experimental points, and the knee fields inflection points B ð1;2Þ K . Inset: Decomposition into the constituent Tsallians. Error bars in all panels are obtained from monoexponential fits. Table 1 Summary of samples and defining characteristics. We consider single-crystal (S.C.) samples both at natural abundance and with 13 C enrichment, as well as diamond particulate (P) samples. capture the decay rates at low and moderate fields, and a constant offset at high field (see Fig. 2b insets).
Field dependence. A second viewpoint of the data, presented in Fig. 2c, is of the T 1 relaxation times and highlights its highly nonlinear field dependence. There is a step-like behavior in T 1 ðB relax Þ, and an inflection point (knee field) %100 mT beyond which the T 1 's saturate. We can roughly define the knee-field value, B K , as the B relax at which the relaxation rate is twice the saturation R 1 , although verifying this claim is difficult due to the inherently high relaxometry error at high field. This somewhat counterintutive dependence has significant technological implications. (i) Long 13 C lifetimes can be fashioned even at relatively modest fields at room temperature. This adds value in the context of 13 C hyperpolarized nanodiamonds as potential MRI tracers 44 , since it provides enough time for the circulation and binding of surface functionalized particles to illuminate disease conditions. (ii) The step-behavior in Fig. 2c also would prove beneficial for 13 C hyperpolarization storage and transport. Exceedingly long lifetimes can be obtained by simply translating polarized diamond particles to modest~100 mT fields-low enough to be produced by simple permanent magnets 40 .
Finally, while the visualizations in Fig. 2b, c cast light on the low and high-field behaviors, respectively, the most natural representation of the wide-field data is on a logarithmic scale (Fig. 2d). The high-density data now unravel the rich relaxation behavior at play in the different field regimes. We discern an additional second inflection point B ð2Þ K at lower magnetic fields below which there is a sudden increase in the relaxation rates. The inset in Fig. 2d shows the decomposition into constituent Tsallian fits with a narrow and broad widths.
Microscopic origins of this relaxation behavior can be understood by first considering the diamond lattice to consist of three disjoint spin reservoirs-electron reservoirs of NV centers, P1 centers, and the 13 C nuclear spin reservoir. P1 centers arise predominantly during NV center production on account of finite conversion efficiency in the diamond lattice. Indeed the P1 centers are typically at 10-100 times higher concentration than NV centers; with typical lattice concentrations of NVs, P1s, and 13 C nuclei, respectively, P NV $ 1 ppm, P e $ 10-100 ppm, and P C $ 10 4 η ppm, where η is the 13 C lattice enrichment level. At any nonzero field of interest, B relax , the electron, and nuclear reservoirs are centered at widely disparate frequencies and do not overlap. There are also negligible effects of level-anticrossings of NVs to 13 C spins, since in all our experiments, the N-V axes are misaligned with the magnetic field. We can separate the relaxation processes in different field regimes to be driven, respectively, by (i) couplings of 13 C nuclei to pairs (or generally the reservoir) of P1 centers. This leads to the B ð1Þ K feature at moderate fields in Fig. 2d; (ii) 13 C spins interacting with individual P1 or NV centers undergoing lattice-driven relaxation (T 1e processes); (iii) inter-nuclear couplings within the 13 C reservoir that convert Zeeman order to dipolar order. Both of the latter processes contribute to the low field B ð2Þ K features in Fig. 2d; and finally, (iv) a high-field process >1 T that shows a slowly varying (approximately constant) field profile. We ascribe this to arise directly or indirectly (via electrons) from two-phonon Raman processes. Since these individual mechanisms are independent, the overall relaxation rate is obtained through a sum, 1 (shown in the inset of Fig. 2d).
Effect of electronic spin bath. Let us first experimentally consider the relaxation process stemming from 13 C spins coupling to the interacting P1 reservoir. While other lattice paramagnetic defects, e.g., V − , N3, and H3 centers 45,46 are also likely to be present, their concentration is much lower than those of the P1 centers, and their effects to 13 C relaxation are indistinguishable from the P1s. In Fig. 3, we consider single-crystal samples of natural 13 C abundance grown under similar conditions but with different nitrogen concentrations (Samples 1-2 in Table 1). Their P1 electron concentrations are P e ¼ 17 ppm and P e ¼ 48 ppm, respectively, and were measured from X-band ESR 47 shown in Fig. 3a (inset). To obtain the data with high density of field points, hyperpolarized relaxometry measurements are taken by an accelerated data collection strategy (outlined in Supplementary Note 3) over a ultra-wide-field range from 1 mT-7 T, with DNP being excited at B pol = 36 mT. For relaxometry at fields below B pol , we employ rapid current switching of Helmholtz coils within the hyperpolarizer device. Both the range of fields, as well as the density of field points being probed are significantly higher than previous studies 37,48 . This aids in quantitatively unraveling the underlying physics of the relaxation processes. We note that probing relaxation behavior below~1 mT in our experiments is currently limited by the finite sample shuttling time, which becomes of the order of the T 1 's being probed.
Experimental results in Fig. 3 reveal a remarkably sharp R 1 dependence, best displayed in Fig. 3a on a logarithmic scale, showing variation in relaxation rate over four orders of magnitude. Each curve fits to a sum of two Tsallian functions (solid line), and reveals the B ð1Þ K inflection point (closely resembling Fig. 2b) beyond which the lifetimes saturate. The second knee field B ð2Þ K at ultralow fields can also be discerned, although determining its exact position is difficult without relaxation data approaching truly zero field. Comparing the two samples (Fig. 3a), we observe a clear correlation in the B ð1Þ K kneefield values shifting to higher fields at higher electron concentration P e . The high-field relaxation rates, highlighted in Fig. 3b, increase with P e . Interestingly at low-fields (see Fig. 3c), the diamond with lower P e (Sample 1) has an enhanced relaxation rate, yielding an apparent "cross-over" in the relaxation data between the two samples at %50 mT.
While we have focused here on single crystals, we observe quantitatively identical relaxation behavior also for microdiamond powders (e.g., Samples 3-4 in Table 1), down to 5-µm sizes (see Fig. 4). Indeed, hyperpolarized particulate samples have a potentially wider application space, both as polarization agents and contrast agents in MRI, and the results indicate that the random orientations of the crystallites play no significant role in the dominant P1-driven nuclear relaxation process. We do expect, however, that for nanodiamond particles <100 nm, surface electronic spins will cause an additional relaxation channel.
Let us now develop a simple model to quantify this P1dominated relaxation process. Given the low relative density of the NV centers and consequently weak NV-NV couplings, to a good approximation they play no role except to inject polarization into the 13 C nuclei. Consider the Hamiltonian of the system, assumed for simplicity to be a single 13 C spin, and the environment-interacting bath of P1 centers surrounding it, H ¼ H S þ H E þ H SE þ H EE where, the first two terms capture the Zeeman parts, the third term is the coupling between reservoirs, and the last term captures the inter-electron dipolar couplings within the P1 bath. Specifically, where I (and S) refer to spin-1 2 Pauli operators on the nuclei (electrons), respectively, and A j zx the pseudo-secular hyperfine interaction that can drive nuclear spin-flips on the 13 C nuclei. For simplicity, we neglect here the effect of the P1 hyperfine couplings to host 14 N nuclei. In principle, they just split the electronic reservoirs seen by the 13 C nuclei into three manifolds separated by the large hyperfine coupling A P1 k % 114 MHz. In the rotating frame at H E , and going into an interaction picture with respect to H EE , the Hamiltonian becomes, Here, A zx h i is the total effective P1-13 C hyperfine interaction, and the norm k H EE k is set by the average dipolar interaction between electronic spins in the bath, henceforth d ee h i. We now make a semi-classical approximation, promotingŜ z to s z ðtÞ, a variable that represents a classical stochastic  process seen by the 13 C nuclear spins 10,49 , In summary therefore, a spin-flipping term I x is tethered to a stochastic variable sðtÞ and this serves as noise on the 13 C spins, flipping them at random instances, resulting in nuclear relaxation upon a time (or ensemble) average. Interestingly, this noise process arises due to electronic flip-flops in the remote P1 reservoir that are widely separated in frequency from 13 C spins. In a simplistic picture, shown in Fig. 1c, relaxation originates from pairs of P1 centers in the same 14 N nuclear manifold (energy-mismatched by δ) undergoing spin flip-flop processes, and flipping a 13 C nuclear spin (when ω L % δ) in order to make up the energy difference. In reality, the overall relaxation is constituted out of the several such processes over the entire P1 electronic spectrum.
Let us now assume the stochastic process s z ðtÞ is Gaussian with zero mean and an autocorrelation function gðτÞ ¼ expðÀτ=τ c Þ with correlation time τ c ¼ 1= d ee h i. The spectral density function SðωÞ ¼ 1 ffiffiffi ffi 2π p R 1 1 gðτÞe Àiωτ dτ that quantifies the power of the spin-flipping noise components at various frequencies is then a Lorentzian, SðωÞ ¼ 2τ c =ð1 þ ω 2 τ 2 c Þ. Going further now into an interaction picture respect to ω L I z , H ðIÞ I ¼ A zx h is z ðtÞ e Àiω L I z t 0 I x e iω L I z t 0 À Á : The survival probability of the spin is, pðtÞ ¼ 1 2 TrfI z e iH ðIÞ I t I z e ÀiH ðIÞ I t g $ e ÀχðtÞ where in an average Hamiltonian approximation, retaining effectively time-independent terms, the effective relaxation rate χðtÞ % R 1 t just can be obtained by sampling of the spectral density resonant with the nuclear Larmor frequency ω L at each field point. This is the basis behind noise spectroscopy of the underlying T 1 process 50 . We recover then the familiar Bloembergen-Purcell-Pound (BPP) result 51,52 , where the relaxation rate, The inter-spin couplings can be estimated from the typical inter-spin distance r e h i ¼ 3=4π  Fig. 3e). Moreover, we expect that these turning points (scaling / P e ) are independent of 13 C enrichment η, in agreement with the data in Fig. 2 (see also Fig. 5).
From lattice considerations (see Supplementary Note 4), we can also estimate the value of the effective hyperfine coupling A zx h i in Eq. (3), which we expect to grow slowly with P e . We make the assumption that there is barrier of r 0 % 2.15 nm around every P1 center in which the 13 C spins are unobservable because their hyperfine shifts exceed the measured 13 with N being the relative number of 13 C spins per P1 spin, and θ j the angle between the P1-13 C axis and the magnetic field, and index j runs over the region between neighboring P1 spins. This gives, For the two samples, we have r e h i ¼ 4.8 nm and 3.39 nm, respectively, giving rise to the effective P1- 13 Figure 3d calculates the resulting relaxation rates from this process R 1 ðω L Þ in a logarithmic plot. It shows good semi-quantitative agreement with the data in Fig. 3a and captures the experimental observation that the rates of the two samples cross over at a particular field. It is instructive to represent the data in terms of effective phase noise (see Fig. 3e), denoted logarithmically as, , where ω 0 ! 0 represents the relaxation rates approaching zero field. Figure 3e shows this for the two samples, employing ω 0 ¼1 mT, and with the estimated field-linewidths displayed by the dashed lines. This makes evident that the high field spin-flipping noise seen by the 13 C nuclei is about 15 dB lower in Sample 1.
While Eq. (3) is the dominant relaxation mechanism operational at moderate fields, let us now turn our attention to the the behavior at ultralow fields in Fig. 3. Equation (2) provides the framework to consider the effect of single P1 and NV electrons to the relaxation of 13 C nuclei. In this case, the stochastic process s z ðtÞ arises not on account of inter-electron couplings, but due to individual T 1e processes operational on the electrons, for instance coupling to lattice phonons. The width of the spectral density is then given by T 1e , While T 1e is also field-dependent, and dominated by two-phonon Raman processes at moderate-to-high field, typical values of T 1e $ 1 ms 53 , give rise to Lorentzian relaxometry widths of % 1 kHz, corresponding to field turning points of B Let us finally comment on the role NV centers themselves as a source of 13 C relaxation. Since NV concentrations are low for samples under study (see Table 1), with NV-NV coupling strengths ≲25 kHz, relaxation processes stemming from pairs of NV centers have a negligible contribution in comparison to single NV-driven processes. Moreover, close-shell 13 C nuclei that are most sensitive to this relaxation channel are unobservable in our experiments since they are hyperfine shifted by δ % A zx Δ sin θ γ e BþΔ cos θ beyond the NMR detection band width %20 kHz, where θ is the angle of the N-V axis to the B = 7 T detection field. Error bars obtained from confidence intervals of exponential decay curves for the 10% and 100% samples, whereas for the 3% sample the error bars emerge from the accelerated data strategy outlined in Supplementary Fig. 4. Characteristic knee-field B ð1Þ K (dashed vertical line) at moderate fields is independent of enrichment, evident in the inset. Knee field at ultralow fields B ð2Þ K is qualitatively indicated by the dashed line that serves as a guide to the eye. Inset: Second derivative of the fitted lines, showing the knee fields at the zero-crossings. c DNP polarization buildup curves also reflect differences in the nuclear spin lifetimes, displaying saturation at much shorter times upon increasing enrichment. DNP in all curves are performed at 36 mT sweeping the entire m s ¼ þ1 manifold. d Polarization buildup times extracted from the data showing that faster nuclear spin relaxation limits the final obtained hyperpolarization enhancements in highly enriched samples. e Spin diffusion constant and diffusion length for 13 C nuclei numerically estimated from the data as a function of lattice enrichment. Dashed line indicates the mean inter-electron distance r NV h i % 12 nm between NV centers at 1 ppm concentration, indicating that spin diffusion can homogeneously spread polarization in the lattice almost independent of 13 C enrichment.
Effect of 13 C enrichment. To systematically probe this low-field behavior as well as consider the effect of couplings within the 13 C reservoir, we consider in Fig. 5 diamond crystals with varying 13 C enrichment η and approximately identical NV and P1 concentrations (Samples 5-7). With increasing enrichment, a third relaxation mechanism becomes operational, wherein at low fields it becomes possible to dissipate Zeeman energy into the dipolar bath. The field dependence of this process is expected to be more Gaussian, centered at zero field and have a width $ d CC h i the mean inter-spin dipolar coupling between 13 C nuclei. We can estimate these couplings from the second moment, where in a lattice of size ', N ¼ N C ' 3 refers to the number of 13 C spins, and the spin density is the angle between the inter-nuclear vector and the direction of the magnetic field. In the numerical simulations outlined in Supple , for instance %39 µT for natural abundance samples, but scaling to %0.46 mT in case of the 100% enriched Sample 7. In real experiments, it is difficult to distinguish between this process and that arising directly from single electrons in Eq. (5), and hence we assign the same label to this field turning point.
Performing hyperpolarized relaxometry (see Fig. 5), we observe that increasing enrichment leads to a fall in nuclear T 1 s, evident both at low (Fig. 5a) and high (Fig. 5b) fields. R 1 rates for the diamonds with 13 C concentrations of 10% and 100% (Samples 6 and 7) are obtained by taking the full relaxation decay curves at every field point, while for the 3% enriched diamond (Sample 5) we use an accelerated data collection strategy (see Supplementary Note 3B) on account of the inherently long T 1 lifetimes. On a logarithmic scale (Fig. 5b), we observe the knee field B ð1Þ K is virtually identical across all the samples, indicating it is a feature independent of 13 C enrichment, originating from interactions with the electronic spin bath. This is in good agreement with the model in Eq. (3). A useful means to evaluate the inflection points arises from the zeros of the second derivative of the Tsallian fits, as indicated in the inset of Fig. 5a. Moreover, the lower inflection field B ð2Þ K scales to higher fields with increasing enrichment η, pointing to its origin from internuclear dipolar effects. At the low fields, we also notice that the samples with lower enrichment have higher relaxation rates, and with steeper field-profile slopes (Fig. 5b). This is once again consistent with the model that the spectral density height and width being probed scales with d CC h i. Changes in the nuclear lifetimes are also reflected directly in the DNP polarization buildup curves, shown in Fig. 5c. We perform here hyperpolarization of all samples under the same conditions, sweeping the entire m s = +1 manifold at B pol = 36 mT, sweeping over the full NV ESR spectrum. We note that polarization buildup is predominantly monoexponential (with time constants shown as dashed lines in Fig. 5c), except for at natural abundance 13 C, where biexponential growth (solid lines) is indicative of nuclear spin diffusion. Data demonstrate that highly enriched samples have progressively smaller polarization buildup times (see Fig. 5d) on account of limited nuclear lifetimes at B pol .
Moreover, the experimental data allow us to quantify the homogenization of polarization in the lattice. We assign a spin diffusion coefficient D ¼ r n h i 2 30T 2n (see Fig. 5e), where the T 2n are evaluated here by only taking the dipolar contribution to the linewidth, T 2n % 1= d CC h i 54 . Given a total time bounded by T 1 , we can calculate the rms overall diffusion length 55 as σ ¼ ffiffiffiffiffiffiffiffiffiffiffi 2DT 1 p that is displayed as the blue points in Fig. 5e. Also for reference is plotted the mean NV-NV distance %12 nm at 1 ppm concentration (dashed region in Fig. 5e), indicating that to a good approximation that the optically pumped polarization reaches to all parts of the diamond lattice between the NV centers.
We comment finally that determining the origins of 13 C relaxation in enriched samples can have several technological applications. Enrichment provides an immediate means to realize quantum registers and sensing modalities constructed out of hybrid NV-13 C spin clusters, and as such ascertaining nuclear relaxation profiles is of practical importance for such applications. Low η ( 3%) naturally engender NV-13 C pairs that can form quantum registers [56][57][58] . The nuclear spin can serve as an ancillary quantum memory that, when employed in magnetometry applications, can provide significant boosts in sensing resolution 19,59 . With increasing 13 C concentrations η ≳ 10% a single NV center can be coupled to several 13 C nuclei forming natural nodes for a quantum information processor, and where the nuclear spins can be actuated directly by hyperfine couplings to the NV electron 60,61 . Approaching full enrichment levels (η ¼ 100%), internuclear couplings become significant, permitting hybridized nuclear spin states and decoherence-protected subspaces 62 for information storage. In bulk quantum sensing too, for instance applied to diamond-based gyroscopes 18,63 , the high density of 13 C sensor spins ($10 22 cm −3 ), as much as >10 5 times the number of NV centers, can be harnessed to increase sensitivity.

Discussion
Experimental results in Fig. 3 and Fig. 5 substantiate the 13 C relaxation pathways operational at different field regimes, and highlight the important role played by the electronic reservoir toward setting the spin lifetimes. To independently affirm this connection, we perform in Fig. 6 relaxometry under blue (≲ 495 nm) wavelengths where the P1 electrons ionize strongly. In these experiments, we operate at the relatively low excitation powers 240 mW mm −2 due to technical limitations related to sample heating. We observe a comparative decrease in nuclear T 1 with respect to decay in the dark (see Fig. 6). In contrast, we do not observe significant change in the lifetimes under 520 nm excitation. We believe that we are accessing an intermediate regime where the P1 ionization rate is slower than the interelectron flipflop rate by d ee h i. As a result, ionization is not rapid enough to decouple each spin flip-flop event between neighboring P1 centers. On the contrary, upon ionization and subsequent recapture, the P1 electrons can now cause spin relaxation over a lattice length scale that is longer than purely given by dipole-dipole interactions. This stirring of the electronic spin bath through blue irradiation results in an increased spectral density component at the nuclear Larmor frequency and consequently an increased relaxation rate.
Our experimental observations open the door to a number of intriguing future directions. First, they suggest the prospect of increasing nuclear lifetimes by raising the NV center conversion efficiency 64 . More generally, they point to the efficacy of materials science approaches toward reducing paramagnetic impurities in the lattice. Finally, it opens the possibility of employing coherent quantum control for dissipation engineering, to manipulate the spectral density profile seen by the nuclei and consequently lengthen their T 1 . Applying a pulse sequence to increase T 1 has been a longstanding goal in magnetic resonance 65,66 , but is typically intractable because of inability to coherently control broadspectrum phonon interactions. Instead here since the nuclear T 1 stems from electronic T 2e processes, these can be echoed out; In particular, the application of electron decoupling (such as WAHUHA 67 or Lee-Goldburg 68 decoupling) on the P1 spin bath would suppress the inter-electron flip-flops, narrow the noise spectral density, and consequently shift the knee field B ð1Þ K to lower fields. Such T 1 gains just by spin driving at room temperature and without the need for cryogenic cooling, and consequent boosts in the hyperpolarization enhancements-scaling by the decoupling factor-will have far-reaching implications for the optical DNP of liquids under ambient conditions. Given the multi-frequency microwave control driving each of the 14 N manifolds would entail 69 , an attractive alternate all-optical means is via the optical ionization of P1 centers faster than their flip-flop rate. The exact interplay between optical ionization and recapture rates required for T 1 suppression will be the subject of future work.
Employing hyperpolarized relaxometry, we have mapped the 13 C nuclear spin lifetimes in a prototypical diamond quantum system over a wide-field range, in natural abundance and enriched 13 C samples, and for both single crystals as well as powders. We observe a dramatic and intriguing field dependence, where spin lifetimes fall rapidly below a knee field of~100 mT. The results indicate that the spin lifetimes predominantly arise from nuclear flip processes mediated by the P1 center electronic spin bath, and immediately opens the compelling possibility of boosting nuclear lifetimes by quantum control or optically induced electronic ionization. This has significant implications in quantum sensing, in building longer lived quantum memories, and in practically enhancing the 13 C hyperpolarization efficiency in diamond, with applications to hyperpolarized imaging of surface functionalized nanodiamonds and for the DNP of liquids brought in contact with high surface area diamond particles.

Methods
Materials. Table 1 summarizes the particulars of the samples we employ in this study. 13 C-enriched diamonds (Samples 5-7) used to conduct experiments in Fig. 2 and Fig. 5 were grown through chemical vapor deposition using a 13 C enrichment mixture of methane and nitrogen (660 ppm, Applied Diamond Inc) as precursor followed by 13 C enrichments of 10%, 25%, 50%, and 100% to produce the respective percent-enriched diamonds 30 . To produce an NV-concentration of 1-10 ppm, the enriched samples were irradiated with 1 MeV electrons at a fluence of 10 18 cm −2 (Prism Gem LLC) then annealed for 2 h at 800°C. The natural abundance samples used in Fig. 3 (Samples 1-2) were grown under synthetic highpressure, high-temperature conditions (Element 6, Sumitomo) 29 then annealed for 1 hour at 850°C. The NV and P1 concentration were measured to be 1.4 ± 0.02 ppm and 17 ± 2 ppm for the first sample and 6.9 ± 0.8 ppm and 48 ± 6 ppm for the second sample, respectively. The microdiamond powders in Fig. 4 (Samples 3-4), produced by HPHT techniques, were acquired from Element 6 and Columbus Nanoworks respectively.

Data availability
The data necessary to replicate the results of this paper are available from the corresponding author upon reasonable request.  c Exemplary decay curves obtained at 20 mT. d Relaxation rates on linear field scale. We observe that the blue irradiation leads to a decrease in 13 C nuclear lifetimes, which we hypothesize arises from fluctuations introduced in the electronic spin bath upon recapture after P1 center ionization. This illustrates that the electronic spin spectral density can be optically manipulated, and potentially ultimately also narrowed under sufficiently high-power ionization irradiation.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/ licenses/by/4.0/.