Articles | Volume 2, issue 1
Magn. Reson., 2, 33–48, 2021

Special issue: Geoffrey Bodenhausen Festschrift

Magn. Reson., 2, 33–48, 2021

Research article 17 Feb 2021

Research article | 17 Feb 2021

Room-temperature hyperpolarization of polycrystalline samples with optically polarized triplet electrons: pentacene or nitrogen-vacancy center in diamond?

Room-temperature hyperpolarization of polycrystalline samples with optically polarized triplet electrons: pentacene or nitrogen-vacancy center in diamond?
Koichiro Miyanishi1, Takuya F. Segawa2,3, Kazuyuki Takeda4, Izuru Ohki5, Shinobu Onoda6,7, Takeshi Ohshima6,7, Hiroshi Abe6,7, Hideaki Takashima8, Shigeki Takeuchi8, Alexander I. Shames9, Kohki Morita5, Yu Wang4, Frederick T.-K. So2,6, Daiki Terada2,6, Ryuji Igarashi6,7,10, Akinori Kagawa1,10,11, Masahiro Kitagawa1,11, Norikazu Mizuochi5, Masahiro Shirakawa2,6, and Makoto Negoro6,10,11 Koichiro Miyanishi et al.
  • 1Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
  • 2Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Nishikyo-Ku, Kyoto 615-8510, Japan
  • 3Laboratory for Solid State Physics, ETH Zurich, 8093 Zurich, Switzerland
  • 4Division of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan
  • 5Institute for Chemical Research, Kyoto University, Gokasho, Uji, Kyoto 611-0011, Japan
  • 6Institute for Quantum Life Science, National Institutes for Quantum and Radiological Science and Technology, 4-9-1, Anagawa, Inage-Ku, Chiba 263-8555, Japan
  • 7Takasaki Advanced Radiation Research Institute, National Institutes for Quantum and Radiological Science and Technology, 1233 Watanuki, Takasaki, Gunma 370-1292, Japan
  • 8Department of Electronic Science and Engineering, Kyoto University, Nishikyo-ku, Kyoto 615-8510, Japan
  • 9Department of Physics, Ben-Gurion University of the Negev, 8410501 Beer-Sheva, Israel
  • 10JST, PRESTO, Kawaguchi, Japan
  • 11Center for Quantum Information and Quantum Biology, Institute for Open and Transdisciplinary Research Initiatives, Osaka University, 1-2 Machikaneyama, Toyonaka, Osaka 560-0043, Japan

Correspondence: Koichiro Miyanishi (, Takuya F. Segawa (, and Makoto Negoro (


We demonstrate room-temperature 13C hyperpolarization by dynamic nuclear polarization (DNP) using optically polarized triplet electron spins in two polycrystalline systems: pentacene-doped [carboxyl-13C] benzoic acid and microdiamonds containing nitrogen-vacancy (NV) centers. For both samples, the integrated solid effect (ISE) is used to polarize the 13C spin system in magnetic fields of 350–400 mT. In the benzoic acid sample, the 13C spin polarization is enhanced by up to 0.12 % through direct electron-to-13C polarization transfer without performing dynamic 1H polarization followed by 1H−13C cross-polarization. In addition, the ISE has been successfully applied to polarize naturally abundant 13C spins in a microdiamond sample to 0.01 %. To characterize the buildup of the 13C polarization, we discuss the efficiencies of direct polarization transfer between the electron and 13C spins as well as that of 13C−13C spin diffusion, examining various parameters which are beneficial or detrimental for successful bulk dynamic 13C polarization.

1 Introduction

Nuclear magnetic resonance (NMR) spectroscopy and magnetic resonance imaging (MRI) are powerful tools covering fields from physics, chemistry, and biology to medicine. However, the poor sensitivity remains the Achilles heel to all magnetic resonance experiments. Dynamic nuclear polarization (DNP), a technique to transfer spin polarization from electrons to nuclei, has been extensively studied since its early discovery (Overhauser1953) and has seen an impressive revival in the current century (Ardenkjær-Larsen et al.2003; Maly et al.2008; Lesage et al.2010; Nelson et al.2013). In DNP using unpaired electrons as the sources of polarization, the polarization enhancement factor ϵ is limited to γe/γn, where γe(n) are the gyromagnetic ratios of the electron (nuclear) spins. To attain nuclear polarization of as high as ∼10 %, DNP needs to be performed at very low temperature (<20 K) and in high magnetic fields (>3 T).

Conversely, the spins of optically polarized electrons can have much higher polarization than the thermal equilibrium value. DNP using such attractive sources of polarization leads to nuclear hyperpolarization beyond the limit of DNP using thermal electron polarization. Moreover, the spin polarization of the optically polarized electrons does not depend on the temperature and the magnetic field, except that the magnetic field happens to be such that level (anti)crossings take place (He et al.1993; Corval et al.1994; Sosnovsky and Ivanov2019). It follows that nuclear hyperpolarization is expected even at ambient temperatures and in relatively low magnetic fields. In times of global liquid helium shortage, this could become an elegant method without requiring liquid helium for sample cooling. Microwave-induced DNP using optically polarized electron spins in the triplet state, triplet DNP, was demonstrated in a single crystal of p-dibromobenzene doped with p-dichlorobenzene (Deimling et al.1980), fluorene doped with phenanthrene (van Kesteren et al.1982, 1985), and naphthalene doped with pentacene (Henstra et al.1990; Stehlik and Vieth1992), and this DNP technique using pentacene achieved 1H polarization of 34 % at room temperature in 0.4 T (Tateishi et al.2014). DNP using pentacene has been extended to polycrystalline samples (Takeda et al.2001), and dissolution DNP using pentacene-doped powder samples (Negoro et al.2018; Kagawa et al.2019) has been implemented. Furthermore, the range of molecules, hyperpolarized by triplet DNP, was successfully expanded (Kagawa et al.2018; Tateishi et al.2019; Nishimura et al.2020). In these DNP experiments, the integrated solid effect (ISE) was used as a method for transferring the polarization. The ISE employs microwave irradiation and external magnetic-field sweep, so that the Hartmann–Hahn matching is implemented between the electron spins in the rotating frame and the nuclear spins in the laboratory frame (Henstra et al.1990).

In proton-rich organic solids or frozen solutions, a demonstrated strategy to polarize dilute or low-γ nuclear spins is first to transfer the electron polarization to the protons, then to let the protons undergo spin diffusion to spatially transport the enhanced polarization away from the source electron spins, and finally to perform local polarization transfer from the protons to the target spin species. The protons in rigid solids are known to be efficient carriers of spin polarization. However, when, as is often the case, it is not the protons but other dilute or low-γ spin species that are of NMR-spectroscopic interest, 1H hyperpolarization can be costly, because the relatively large heat capacity of the 1H reservoir consumes considerable resources of electron polarization. Then, a question arises: would direct dynamic polarization of dilute/low-γ spins using the electrons in the triplet state be feasible without the protons being involved in the polarization-buildup process, and if it is, how efficient would that be? In this context, dynamic polarization of dilute/low-γ spins in proton-free inorganic solids is also of interest.

Here, we study dynamic 13C polarization using optically polarized electron spins in the triplet state through direct electron-to-13C polarization transfer, i.e., without being mediated by the 1H spins. In this work, we deal with pentacene doped into a host matrix and negatively charged nitrogen-vacancy (NV) color centers in diamond. In both cases, the optically polarized triplet electron spins are the source for DNP. The electronic structures for these systems are shown in Fig. 1 (Takeda2009; Doherty et al.2013; Rogers et al.2008; Acosta et al.2010; Thiering and Gali2018). In the case of the NV center, laser irradiation excites the electronic state from the ground triplet state 3A2 to the excited triplet state 3E (Fig. 1), which then undergoes intersystem crossing (ISC) to the singlet state 1A1. Importantly, this transition is highly spin selective: while the mS=0 state has a low ISC probability, the mS=±1 state has a high ISC probability. The singlet state 1A1 further decays into a meta-stable singlet state 1E1, which preferentially decays into the mS=0 spin sublevel of the ground triplet state 3A2, hyperpolarizing the mS=0 state. In the case of pentacene, the ground singlet state S0 is photo-excited to the excited singlet state S1. The subsequent transition to the triplet state T3 is caused by ISC, and then to the lowest triplet state T1 by internal conversion, where the spin-orbit coupling selectively populates the triplet sublevels. Moreover, the triplet population depends on the host molecule.

Figure 1Schematic diagrams of the electronic structure for NV and pentacene. The optical transitions are denoted by solid straight arrows, the transitions caused by intersystem crossing are denoted by wavy arrows, and the transitions caused by internal conversion are denoted by dashed straight arrows.


Even though the relevant electron spins are in the triplet state (S=1) for both the pentacene molecules and the NV centers, they differ in many aspects. For the NV centers, the ground state is the triplet state, while in the case of pentacene, the ground state is the singlet state (S=0). The zero-field splitting (ZFS) parameters D and E are 2870 MHz and ∼0 MHz for the NV center and ∼1350 and ∼42 MHz for pentacene in the excited state, the latter of which may vary depending on the host molecule (Yu et al.1984). The differences between these two systems are summarized in Table 1.

Microwave-free optical hyperpolarization for a bulk ensemble of 13C nuclear spins using NV centers in a diamond single crystal was first demonstrated at cryogenic temperature in a magnetic field of 9.4 T (King et al.2010). Optically pumped electron spins were used for nuclear hyperpolarization by exploiting an excited-state level anticrossing at 50 mT, followed by sample shuttling to a magnetic field of 4.7 T for NMR detection (Fischer et al.2013). The method was generalized to a broader range of magnetic fields and correspondingly different orientations of the NV center by adding microwave irradiation (Álvarez et al.2015). For a single crystal of diamond, the 13C polarization of 6 % at room temperature was achieved via DNP as a combination of the thermal mixing and the solid effect (King et al.2015). Recently, DNP using the frequency-swept ISE on NV in powdered microdiamonds in a magnetic field of as low as ca. 30 mT has been reported by Ajoy et al. (2018a, b), who took advantage of the reduced width of the anisotropic electron paramagnetic resonance (EPR) powder pattern of the NV centers.

In the following, we report on optical DNP of 13C spins in microdiamonds containing NV centers and in pentacene-doped [carboxyl-13C] benzoic acid (PBA) microcrystals. While for the diamond samples the 13C isotope is naturally abundant, the carboxyl carbons are isotopically 13C-enriched in PBA microcrystals. We study the behavior of the buildup of the 13C polarization in terms of the efficiency of the polarization transfer from the electron to nuclear spins and that of 13C spin diffusion. As demonstrated below, we obtained 13C polarization of 0.01 % in the microdiamonds and 0.12 % in PBA at room temperature in a magnetic field of ∼0.4 T by the ISE scheme.

In a slightly different context, both the pentacene and NV systems can be used for single-spin optically detected magnetic resonance (ODMR). Pentacene in p-terphenyl was the first single molecule which was detected by magnetic resonance (Köhler et al.1993; Wrachtrup et al.1993). While these experiments were performed at cryogenic temperatures, the NV center in diamond opened the door for single-spin experiments at room temperature (Gruber et al.1997; Jelezko et al.2004). In this article, we will focus on the conventional induction detection of EPR/NMR on ensembles of spins.

Table 1Differences between the two samples, NV in microdiamonds and pentacene in [carboxyl-13C] benzoic acid (α13CBA).

Download Print Version | Download XLSX

2 Theory

Dynamic nuclear polarization using the integrated solid effect

Let us suppose that pulsed laser irradiation creates hyperpolarized electron spins in the triplet state. In the ISE, microwave irradiation and magnetic-field sweep are applied in such a way that the resonances of the individual electron spin packets between two triplet sublevels are adiabatically swept over and the Hartmann–Hahn condition is fulfilled at some point during the sweep, causing the contact between the electron spins in the rotating frame and the nuclear spins in the laboratory frame (Henstra et al.1990). The Hartmann–Hahn condition is expressed as

(1) ω eff , e = ω 0 , C .

Here, ωeff,e represents the electron nutation frequency around the effective field and is given by ω1,e2+(Δωe)2, where ω1,e is the microwave intensity and Δωe is the time-dependent resonance-offset frequency. ω0,C is the Larmor precession frequency of the nuclear spin species of interest, which is 13C in the present case. The exchange of the spin states and thereby of spin polarization is driven by the dipolar interaction between them through the same mechanism as that of cross-polarization (Mehring1983; Peng et al.2006). Thus, the source electron spins need to be locked along the effective field for polarization transfer to take place. This is ensured when the adiabatic sweep starts from far-off resonance, so that the effective field is initially aligned nearly along the static field and gradually tilts.

Since the polarization transfer is driven by the local dipolar interaction, direct hyperpolarization is limited to those nuclei which happen to be located in the vicinity of the electrons. Nevertheless, such a locally hyperpolarized 13C spin state can be transported away by spin diffusion. Thus, by repeating the ISE sequence as described in Fig. 2, the 13C polarization can be accumulated until the DNP buildup and the nuclear spin-lattice relaxation balance and a steady state is established.

Figure 2The pulse sequence for the ISE followed by NMR detection. A laser pulse is used for photoexcitation and optical polarization of the triplet electron spins. Then, a microwave pulse is applied together with a magnetic-field sweep. The ISE pulse sequence is repeated before the enhanced 13C magnetization is detected by applying a radiofrequency π/2 pulse. During the acquisition of the 13C signal, 1H decoupling is to be applied in the case of the benzoic acid sample. The definitions of the symbols used in the figure are as follows. Plaser: laser-beam intensity, ω1,e: microwave intensity, tMW: microwave-pulse width, Bsweep: field-sweep width, R: ISE repetition rate, and N: ISE repetition number.


In general, both the direct electron-to-13C polarization transfer by the ISE and spin diffusion among the 13C spins contribute to the overall buildup of the bulk 13C polarization. Conversely, when 13C spin diffusion is fast such that the 13C polarization is smoothed out within the interval of the repeated ISE sequences, the buildup behavior of the bulk 13C spin polarization PC becomes independent of the spin-diffusion rate. In this rapid-diffusion limit (Takeda2009), the time evolution of PC is governed by

(2) d d t P C R ξ η t ρ e ρ C ( η p P e - P C ) - 1 T 1 , C P C ,

where R is the repetition rate of the ISE sequence, ηt is the fraction of the triplet electrons in two of the triplet sublevels and is obtained as 23-Pe by considering that the populations in the mS=+1 and mS=-1 states are the same, ξ is the exchange probability, i.e., the probability of the spin states being exchanged between the electron and the nuclear spins during a single ISE sequence, ρe is the density of the electrons in the triplet state, ρC is the density of 13C spins, and T1,C is the longitudinal relaxation time of the 13C spins. ηp(≤1), called the active spin-packet fraction, is introduced to deal with such situations where not all but some of the electron spin packets participate in the process of polarization transfer. This is indeed the case for polycrystalline samples, where the anisotropy of the ZFS tensor causes significant broadening to such an extent that the resonance line can only partly be excited. Pe is the electron polarization between two of the triplet sublevels averaged over the ISE sweep time tMW, which is given as

(3) P e = 1 t MW 0 t MW P e exp ( - t / T 1 , e ) d t .

Here, Pe is the initial electron polarization between the relevant two triplet sublevels. T1,e is the time constant introduced to take either spin-lattice relaxation or lifetime decay of the triplet state into account.

Neglecting the thermal polarization of 13C spins, we obtain the solution of this differential equation as


Here, Pfin is the maximum attainable nuclear spin polarization and Tb is the DNP buildup time constant. We note that the initial buildup rate, given as

(6) d d t P C ( t = 0 ) = R ρ e ρ C η t η p ξ P e ,

is independent of the longitudinal relaxation time T1,C.

In the following, we will use Eqs. (4)–(6) for data analysis.

The rapid-diffusion limit can be made valid by performing experiments with a sufficiently slow repetition rate R. Buildup experiments with the slow repetition rate, albeit not advantageous for attaining the highest possible nuclear polarization, give helpful data that allow us to characterize the parameters.

3 Experimental

3.1 Sample preparation

We used the following three samples: (i) 70 mg of 500 µm microdiamonds, with 8.9×1017 cm−3 (5 ppm in atomic ratio) NV concentration and 4.6±0.1×1018 cm−3 (26 ppm) P1 center concentration. So-called “P1 centers” are another type of electron spin defects in diamond with S=1/2 and represent single substitutional nitrogen atoms. (ii) 5.6 mg of 100 nm nanodiamonds, with 5.3×1017 cm−3 (3 ppm) NV concentration, 1.2×1019 cm−3 (70 ppm) P1 center concentration, and a concentration of 3.1×1019 cm−3 (177 ppm) for the total S=1/2 electron spins (P1 centers, dangling bonds, etc.); 100 nm sized nanodiamonds “MICRON+ MDA M0.10” were purchased from Element Six, UK, and electron irradiated at a fluence of 1019 e/cm2 at room temperature to form vacancies. To create the NV centers, the nanodiamond samples were annealed at 800 C under vacuum (<10-4 Pa). The temperature was increased to 400 C over 1 h and kept for 4 h (“baking”). After that, the temperature was increased to 800 C over 11 h, where the sample was annealed for 2 h. The temperature was decreased to 350 C in 1 h and then down to room temperature. These “dark” nanodiamonds (due to amorphous sp2 carbon on the surface) were “made white” by oxidation in air at 550 C for 2 h, followed by boiling acid cleaning in H2SO4/HNO3 (3:1, 125 C) for 3 d (Terada et al.2019). Both diamond samples contain 1.1 % natural abundance 13C. (iii) 4 mg of PBA powder, with 0.04 mol% (2.6×1018 cm−3) concentration of pentacene. The [carboxyl-13C] benzoic acid was purchased from Cambridge Isotope Laboratories. The characteristics of the samples are summarized in Table 2. All the samples were packed in a glass tube, and the pictures of these samples are shown in the Supplement.

Table 2Summary of sample conditions and the optimized parameters of the DNP experiments.

Download Print Version | Download XLSX

3.2 Experimental setup

The procedure and experimental setup of triplet DNP used in this work are similar to those described by Tateishi et al. (2014).

A solid-state laser with a wavelength of 527 nm, pulse length of 200 ns, and pulse energy of 30 mJ was used as a light source for the excitation of NV centers, and a dye laser with a wavelength of 594 nm, pulse length of 200 ns, and pulse energy of 6 mJ was used as a light source for the excitation of pentacene. A static magnetic field from 0.3 to 0.5 T was generated by an electromagnet. All experiments were carried out at room temperature.

The EPR experiments were done with a home-built spectrometer, similar to that reported by Yap et al. (2015), where an RF pulse at 400 MHz generated by an arbitrary waveform generator was converted to 11.6 GHz in a superheterodyne architecture. The microwave pulse was amplified to 1 W and led into the cavity.

4 Results and discussion

4.1 EPR measurements

Using the pulse sequence depicted in Fig. 3a, magnetic-field dependencies of the amplitude of the spin echoes of the electrons in the triplet state were measured for the three samples described above.

Figure 3(a) EPR echo pulse sequence used in the measurement. τlaser is the delay time used for the lifetime decay measurement and τecho is the echo time. (b–f) EPR spectra of polycrystalline samples under study obtained by experiment and simulation: (b) optically polarized EPR powder spectrum of microdiamonds with MW frequency of 11.63 GHz. (c) EPR powder spectrum of microdiamonds with MW frequency of 11.63 GHz without the laser pulse irradiation. (d) Optically polarized EPR powder spectrum of nanodiamonds with MW frequency of 11.52 GHz. (e) EPR powder spectrum of powder BA doped with 0.04 mol% pentacene with MW frequency of 11.66 GHz. (f) Simulated EPR powder spectra of thermal state triplet electron in NV center and pentacene at 11.6 GHz. The shaded area is the magnetic field-sweep range Bsweep for each system.


The optically excited NV spectrum of microdiamonds in Fig. 3b does not show a smooth polycrystalline pattern (powder spectrum) but rather a partially averaged crystalline pattern due to the average diameter of the particles around 500 µm and the number of microdiamonds in the sample tube being only ≃200. The echo time τecho was set to 2.3 µs and the measurement took 3 h.

By comparing the signal intensity with and without laser irradiation (see the Supplement), we estimate an enhancement of ca. 170. This corresponds to an electron spin polarization of ca. 16 %, derived from a thermal polarization of ca. 0.095 %.

The three lines in the center of the field (around g=2) do not stem from the NV center, but from the P1 centers and other S=1/2 electron spin defects. The hyperfine structure observed arises from the coupling of the electron spin to the adjacent 14N nucleus in the P1 centers. This signal is not optically polarized and is also present “in the dark” as shown in Fig. 3c for the case of microdiamonds.

Figure 3d shows the optically polarized NV powder spectrum of the nanodiamonds. The echo time τecho was set to 400 ns, and the measurement took 7 h. The line exhibited a typical dipolar powder pattern (Pake pattern), where the two “horns” are separated by the zero-field splitting parameter D and the two “shoulders” by 2 D. The horns come from the crystallites with such orientation that the N-V axis is perpendicular to the magnetic field, whereas the shoulders correspond to those in which the N-V axis is along the field. The result that one of the transitions is inverted is ascribed to the optical polarization mechanism, which exclusively populates the electron spins of the ground triplet state in the NV center into the mS=0 state. To the best of our knowledge, this is the first report on the full powder EPR spectrum of optically polarized nanodiamonds observed with the conventional (i.e., not optical) detection. A similar spectrum from a large size distribution of nano-/micro-diamonds (up to ca. 1.5 µm) recently showed an enhancement of ca. 6-fold (Mindarava et al.2020). Importantly, an optically polarized NV signal was not seen for the “dark” nanodiamond powder before the chemical surface cleaning (see the Supplement). We assume that the visible laser light was absorbed on the particle surface.

Figure 3e shows the powder spectrum of the photo-excited triplet state of pentacene. In the measurement, the echo time τecho was set to 2.3 µs and the measurement took 5 h. The shape is also a dipolar powder pattern, where the two “horns” are separated by the zero-field splitting parameter D and the two “shoulders” by 2 D. The zero-field splitting parameter D for pentacene is only about half of that for NV, which explains half the powder “linewidth” of pentacene compared to NV. In the case of pentacene, the part of the “horn” is broadened due to the finite ZFS parameter E, which is zero for NV. Figure 3f shows EPR spectra of the thermally populated NV center and pentacene simulated using EasySpin, a Matlab package (Stoll and Schweiger2006).

4.2 Hyperpolarized 13C NMR

4.2.1 NV-containing microdiamonds

We performed hyperpolarization of the 13C spins in microdiamonds using the ISE pulse sequence in a magnetic field of 0.36 T, which corresponds to the position of the low-field “horn” in the spectrum of Fig. 3b. To find the experimental parameters that maximize the efficiency of DNP, we varied the range of the magnetic-field sweep BsweepNV, the width tMWNV, and the amplitude ω1,eNV of the microwave pulse and examined the enhanced 13C magnetization. As demonstrated in the upper part of Fig. 4a–c, the optimal conditions were found to be BsweepNV5 mT, tMWNV1250µs, and ω1,eNV3.74 MHz.

Figure 4(a–c) Optimization of the ISE pulse sequence for microdiamonds (upper row, filled circles) and pentacene-doped benzoic acid (PBA) (lower row, open circles) through adjustment of (a) the field-sweep width Bsweepα, (b) the microwave duration tMWα, and (c) the microwave intensity ω1,eα, where α=NV,P for microdiamonds and PBA, respectively. The values of the fixed parameters are indicated in each graph. The gray lines in (c) indicate the 13C Larmor frequency for the magnetic field that corresponds to the low-field horn of the EPR spectrum. (d) Longitudinal relaxation (upper graph) and lifetime decay (lower graph) of the electron spins in the triplet state of NV centers in microdiamonds and pentacene, respectively.


Then, adopting these parameters, we performed dynamic 13C polarization by repeating the ISE sequence at a rate of 100 Hz for 240 s, measured 13C NMR at a Larmor frequency ω0,C of 3.85 MHz, and successfully obtained a hyperpolarized 13C spectrum of the microdiamonds, as demonstrated in Fig. 5a. The result of the identical 13C measurement except that DNP was not performed, also plotted in Fig. 5a for comparison, did not show any appreciable sign of the signal above the noise level. Figure 5b shows buildup curves of 13C polarization with various ISE repetition rates R ranging from 10 to 100 Hz. For R of up to 60 Hz, the buildup rate and the finally attained 13C polarization increased with R, whereas they saturated for R>60 Hz. The maximum 13C polarization was 0.01 %, corresponding to 324-fold enhancement of 13C polarization compared to that in thermal equilibrium.

Figure 5Optical hyperpolarization in microdiamonds. (a) NMR spectra of 13C in microdiamonds averaged over four times. The solid line shows the 13C NMR signals with DNP and the broken line shows that without DNP. The spectrum was obtained with a laser repetition frequency of 100 Hz and 240 s of DNP buildup time. (b) Buildup curves of 13C polarization in microdiamonds obtained with ISE repetition rates ranging from 10 to 100 Hz. (c) The relaxation curve for the 13C NMR signal of microdiamonds. (d) The DNP and EPR field profile around 0.36 T.


Figure 5c shows the profile of 13C depolarization after hyperpolarization, from which the time constant T1,CDia of 13C longitudinal relaxation was determined to be 99±14 s. Here, 13C relaxation is mainly caused by the P1 centers (Ajoy et al.2019) creating fluctuating local fields at the 13C atomic sites. Figure 5d compares the magnetic-field dependences of the enhanced 13C magnetization and of the EPR signal obtained in the microdiamond sample at around 0.36 T. The result that they showed similar profiles indicates that the underlying polarization mechanism is indeed the ISE.

To determine the exchange probability ξNV appearing in Eq. (2) under the optimized ISE conditions, it is necessary to estimate the DNP-active spin-packet fraction ηpNV, the triplet fraction ηtNV of the NV in the microdiamonds, and the electron polarization PeNV averaged over the ISE sweep time. These parameters can be extracted from the experimental powder EPR spectra shown in Fig. 3 and the calculated orientational dependence of the resonance magnetic field shown in Fig. 6. The optimized range of the magnetic-field sweep, indicated by the shaded areas in Fig. 3b and f, corresponds to the region between the black lines in Fig. 6a. From the area of this region compared to that of the whole solid angle, the active spin packet fraction ηpNV was estimated to be 0.204.

Figure 6Orientational dependence of the resonant magnetic field for the electron spins in the triplet state of (a) the NV center and (b) pentacene doped in benzoic acid at a microwave frequency of 11.66 GHz. The x, y, and z axes represent the principal axes of the ZFS tensor. (c) Orientational dependence of the polarization of the electron spins in the photo-excited triplet state of pentacene doped in benzoic acid in 0.39 T calculated from the zero-field population (Yu et al.1984). The upper and lower rows correspond to the EPR transitions between the mS=+1 and mS=0 states and between the mS=0 and mS=-1 states, respectively. The regions between the black lines correspond to the orientations for which the crystallites experience polarization transfer by the ISE using the optimized field-sweep range found in the experiments demonstrated in Fig. 4.


Having estimated the electron polarization PeNV of the NV center between the mS=0 and mS=-1 states to be 16 %, we then obtained the triplet fraction ηtNV=2/(3-PeNV)=0.7. From the measured longitudinal relaxation time T1,eNV of 2.5 ms (Fig. 4d), we obtained the averaged electron polarization PeNV over the optimal pulse width tMWNV=1250µs as 0.125.

In evaluating the exchange probability ξNV, it is important to use such experimental data as were taken under the rapid-diffusion regime that validates Eq. (4). Accordingly, the buildup curve with the ISE repetition rate of 10 Hz, which was rather low and not suitable for attaining the highest polarization, was adopted for the purpose of data fitting to Eq. (4). From the final polarization Pfin=0.0026 %, the time constant Tb=47.4 s of the buildup curve at the repetition rate R of 10 Hz, the electron polarization averaged over the ISE sweep time PeNV0.125, the active spin-packet fraction ηpNV0.204, the triplet fraction ηtNV0.7, the density of triplet electrons in the microdiamonds ρeNV8.9×1017 cm−3, and the 13C spin density ρCDia1.9×1021 cm−3, we obtained the exchange probability ξNV as 0.0067.

4.2.2 Pentacene-doped benzoic acid

We also successfully hyperpolarized the 13C spins in PBA through direct electron-to-13C polarization transfer by the ISE. Our results contrast with the previous works that relied on prior 1H polarization followed by 1H−13C cross-polarization (Kagawa et al.2019). The lower parts of Fig. 4a–c show the dependence of the efficiency of dynamic 13C polarization by the ISE on the field-sweep width BsweepP, the ISE duration tMWP, and the microwave intensity ω1,eP, experimentally examined at 0.39 T, which corresponds to the EPR at the position where the low-field “horn” appears in the spectrum of Fig. 3e. The optimal conditions were found to be BsweepP20 mT, tMWP30µs, and ω1,eP3.44 MHz.

The optimal values for the field-sweep width BsweepP and the microwave intensity ω1,eP differed from those (BsweepNV, ω1,eNV) in the case of the microdiamond sample just by small factors, whereas the optimal microwave-pulse width tMWP was shorter than tMWNV by more than an order of magnitude. Such a striking difference can be ascribed to the sources of polarization for PBA being the electron spins in the metastable, photo-excited triplet state of pentacene, while those for the diamond sample are the electron spins in the persistent, ground triplet state. For the PBA sample, the contact time is limited by the lifetime of the photo-excited triplet state of pentacene, whereas for the microdiamond sample it is by the electron spin-lattice relaxation time T1ρ in the rotating frame. As we shall discuss below, it is not those 13C spins adjacent to the NV centers but those located at a moderate distance that need to receive the polarization of the electron spins of the NV center, in order for the enhanced polarization to be eventually transported away by 13C spin diffusion. Since the interaction between such 13C spins and the NV centers is expected to be relatively weak, direct polarization transfer by the ISE needs to be performed for a relatively longer time duration, as long as the electron magnetization is retained along the effective field in the rotating frame.

13C NMR measurements were performed at 4.19 MHz. An enhanced 13C NMR spectrum, in comparison with that obtained without performing the ISE, is demonstrated in Fig. 7a. Figure 7b shows 13C polarization buildup behaviors for various ISE repetition frequencies. With the highest experimentally feasible ISE repetition rate R of 90 Hz, the 13C polarization finally reached 0.12 %, i.e., ∼3600 times the thermal 13C polarization. The 13C longitudinal relaxation time T1,CBA was determined to be 474±30 s (Fig. 7c). We estimated the active spin-packet fraction ηpPBA and the triplet fraction ηtP of pentacene in a similar way to the case of the microdiamond sample. From comparison of the optimal magnetic-field sweep range Bsweep indicated by the shaded region in Fig. 3e and f with the area between the black lines in Fig. 6b and c, ηpPBA was estimated to be 0.565. The average electron polarization PeP between the two triplet sublevels over the field-sweep time was estimated using the populations over the zero-field eigenstates of the triplet state of pentacene doped in benzoic acid, which is known to be Px:Py:Pz=0.44:0.34:0.22 (Yu et al.1984). Assuming that the external magnetic field is nearly perpendicular to the z axis of the principal axis system of the ZFS tensor for the relevant electron-spin packets, we calculated the populations over the triplet sublevels in the magnetic field and obtained PeP=0.105 and ηtP=2/(3-PeP)=0.69. Then, taking account of the lifetime decay with the time constant τeP=9µs (lower row of Fig. 4d), we determined PeP to be 0.03.

Using Eq. (4), we performed curve fitting of the buildup data experimentally obtained with the ISE repetition rate R of 10 Hz, at which the rapid-diffusion limit is expected to be valid, and obtained PfinPBA(R=10Hz)=0.023 % and TbPBA(R=10Hz)=141 s, whence, with PeP0.03, ηpPBA0.535, ηtP0.69, ρeP2.6×1018 cm−3, and ρCPBA6.26×1021 cm−3, we determined ξPBA to be 0.035.

For the microdiamond and PBA samples, the exchange probabilities were found to be ξPBA=0.035 and ξNV=0.0067, respectively. To account for such a significant difference in the probability of the spin states being transferred between the electron and the 13C spins in a single shot of the ISE sequence, we again note the different types of the sources of polarization. In the case of the microdiamond sample, the relevant electron spins are in the ground triplet state, persistently causing significant local fields at the 13C sites in the vicinity. They would create the 13C spin-diffusion barriers (Wenckebach2016), in which the 13C spins would not be able to transport their polarization to other 13C spins via the mutual spin flip-flop process. Unfortunately, those 13C spins that are most likely to receive spin polarization from the NV centers are inside the barrier and thus are least likely to distribute the enhanced polarization away.

Nevertheless, the experimental result that the bulk enhancement of the 13C polarization in the microdiamonds was indeed realized indicates the presence of those 13C spins in the sample that are capable of both receiving the polarization from the NV center and passing it over other 13C spins through spin diffusion. Such 13C spins ought to be located at a moderate distance from the NV center just outside the barrier. That is, if the 13C spins are too close to the NV center, the interaction with the electron spin would overwhelm the dipolar interaction among the 13C spins, hindering the flip-flop transitions between the 13C spins. The relatively long distance between such mediating 13C spins and the NV centers would result in the low exchange probability ξNV.

Conversely, the 13C spins in PBA do not suffer from the local fields produced by the electron spins in the triplet state for most of the time during the buildup experiment, because of the transient nature of the excited triplet state. Indeed, the lifetime of the triplet scale, found to be ca. 9µs (Fig. 4d), is 3 orders of magnitude shorter than the time interval between the ISE sequence even with its highest feasible repetition rate. For the 13C spins located relatively close to the pentacene molecules, the probability of receiving the polarization from the electron spins in the triplet state is expected to be relatively high. Even if they ought to temporarily feel the strong local field and are trapped inside the spin-diffusion barrier during the photo-excitation cycles, they are allowed to undergo spin diffusion for most of the time when the pentacene molecules are in the ground, diamagnetic state. Hence, the exchange probability is expected to be higher.

Figure 7(a) 13C spectra of PBA with and without DNP. The enhanced spectrum (solid line) was obtained after performing the ISE sequence at a repetition frequency of 90 Hz for 600 s. (b) Buildup curves of 13C polarization in PBA for ISE repetition rates ranging from 10 to 90 Hz. Each data point represents the 13C polarization derived from the intensity of the NMR signal accumulated over 16 times. (c) Dependence of the 13C magnetization in PBA on the time interval between DNP and NMR measurement.


4.313C spin diffusion

4.3.1 In diamond

In addition to the efficiency ξ of direct polarization transfer from the electron spins in the triplet state to the 13C spins, spin diffusion among the 13C spins is another important factor that affects the overall enhancement of polarization. Spin diffusion is driven by the flip-flop component of homonuclear dipolar interactions, and in turn, the dipolar interaction is determined by the geometrical configuration of the relevant spins. It follows that spin diffusion is characterized by the geometrical configuration. According to Lowe and Gade (Lowe and Gade1967), the (α,β)-component (α,β=x,y,z) of the spin-diffusion tensor 𝒟αβ is given by

(7) D α β = π 2 i A 0 i α 0 i β 0 i Δ 0 i - 1 2 .

Here, 0 in the subscripts represents one of the 13C sites arbitrarily chosen to be the origin of the coordinate system, and the sum is taken over all 13C sites. αij(βij) is the α(β) component of the internuclear vector rij between sites i and j. Aij and Δij are expressed as


where μ0=4π×10-7 Hm−1 is the vacuum permeability, γ is the gyromagnetic ratio of 13C, and θij is the angle between the internuclear vector rij and the external static field.

To estimate the 13C spin-diffusion coefficient in the diamond sample outside the diffusion barrier where the effect of the hyperfine couplings is negligibly small, we made a table of vectors representing the carbon atomic sites in diamond. We then randomly picked up the 13C sites assuming the natural abundance of 1 %. Using the Lowe–Gade formula, we calculated the isotropic spin-diffusion coefficient D=(Dxx+Dyy+Dzz)/3 by taking the 13C sites inside the spherical region with various radii into account. As shown in Fig. 8, the calculated value of 𝒟 increased with the size of the region that we considered and asymptotically converged to 7.28×10-18 m2 s−1.

Figure 813C spin-diffusion coefficients calculated for diamond crystal using the Lowe–Gade formula (Eq. 7) by taking account of naturally abundant 13C sites within spherical regions with various radii; 750 13C sites are involved within the sphere with a radius of 4 nm.


4.3.2 In PBA

Unlike the case of diamond, where all atomic sites other than 13C are magnetically inert, estimation of the 13C spin-diffusion coefficient in PBA by the Lowe–Gade formula cannot be made in a straightforward way, because of the presence of the abundant 1H spins causing the considerable dipolar fields at the 13C sites. The effect of the 1H−13C dipolar interaction is to lift up the degeneracy in the energy levels of the 13C spin packets, so that the flip-flop process among the 13C spins tends to not conserve the energy. Thus, the 13C diffusion rate ought to be slower than in the case if it were not for the proton spins.

In this work, we estimated the 13C spin-diffusion coefficient in [carboxyl-13C] benzoic acid from the experimental repetition-rate dependence of the initial buildup rate of the 13C polarization (Kagawa et al.2009; Takeda2009). With relatively low ISE repetition rates for which the rapid-diffusion condition is valid, the buildup rate is proportional to the repetition rate. As the repetition rate is increased, spin diffusion would no longer be able to transport the polarization completely during the time interval of ISE repetition, and the repetition-rate dependence of the buildup rate begins to saturate. The saturation of the initial buildup rate of 13C polarization was observed in the case of PBA (Fig. 9a). Since the initial buildup rate is independent of relaxation, one can estimate the spin-diffusion rate from the data plotted in Fig. 9a without having to take the effect of the former into account.

In order to reproduce the profile of the experimentally obtained initial buildup rate that increased with the repetition rate R and exhibited saturation at relatively higher values of R, we considered a cubic region with a side length of 10.4 nm, which includes one single pentacene molecule on average, and supposed that the position-dependent 13C polarization evolves in time according to the diffusion equation with a periodic boundary condition. We also assumed that the timescale of spin diffusion is much longer than the time interval tMW of the ISE sequence, so that the point source inside the cubic region instantly creates the 13C polarization PeP=0.030 at each moment when the ISE sequence is implemented, with a probability ξPBA=0.035.

The profiles of numerically simulated time evolution of the net polarization for various spin-diffusion coefficients 𝒟, ranging from 8.0×10-20 to 1.5×10-19 m2 s−1, are plotted in Fig. 9a. We found that D=9.75×10-20 m2 s−1 is the most likely value for the 13C spin-diffusion coefficient that minimized the residual sum of squares in the case of our 13C-labeled benzoic acid sample (Fig. 9b).

Interestingly, the 13C spin-diffusion coefficient in PBA was estimated to be 2 orders of magnitude smaller than that in diamond, despite the 13C enrichment in the former sample. This is ascribed to the presence of the 1H spins that slow down the flip-flop process between the 13C spins in PBA and to the dense packing of the carbon atoms in the diamond crystal. From the DNP point of view, 13C spin diffusion should be as fast as possible. One way to make it is to continuously apply 1H decoupling throughout the buildup experiment (Negoro et al.2010). This strategy, however, is not practical because of the complexity of the hardware that realizes simultaneous application of radiofrequency and microwave irradiation and of the serious heating of the circuit. Deuteration of the sample can be an alternative way to make 13C spin diffusion faster in PBA.

Figure 9(a) Dependence of the initial buildup rate as a function of the ISE repetition rate R obtained for the PBA sample. The data points indicated by the circles were obtained from the slopes of the buildup curves shown in Fig. 7b at time zero. Solid lines represent the R dependence simulated for various spin-diffusion coefficients 𝒟 according to the model described in the text. (b) A plot of the residual sum of squares calculated for various spin-diffusion coefficients, which gave the minimum for D=9.75×10-20 m2 s−1.


4.4 The behavior of 13C polarization buildup

It is the contribution of both the direct polarization transfer from the electrons in the triplet state to the 13C spins and 13C spin diffusion that eventually leads to bulk 13C hyperpolarization. In reality, spin-lattice relaxation tends to drag the spin system back toward thermal equilibrium. The balance between the buildup and relaxation processes determines the profile of the 13C polarization-buildup curve, the finally attainable bulk 13C polarization Pfin, and the time required to attain Pfin.

The 13C longitudinal relaxation time T1,CDia for the microdiamond sample was 99±14 s, whereas T1,CBA=474±30 s for the PBA sample was longer than the former by a factor of ca. 5. The presence/absence of the paramagnetic electrons in the dark state can be a factor making the difference in the relaxation time. Table 3 summarizes the exchange probability ξ, 13C spin-diffusion coefficient 𝒟, and the spin-lattice relaxation time T1,C. As discussed above, in PBA, the average electron polarization Pe and the 13C spin-diffusion coefficient were found to be lower than those in microdiamonds. Nevertheless, the higher exchange probability ξPBA and longer longitudinal relaxation time T1,C more than compensate for the lower electron polarization and the slower 13C spin diffusion, resulting in the higher final 13C polarization.

In order to examine the possibility that the spin polarization is leaking into the proton spin system in PBA, we implemented the ISE sequence in PBA with the parameters that we found to be optimal for polarizing the 13C spins and then examined whether the 1H magnetization was enhanced. We found that the 1H polarization was indeed built up to 0.16 % (not shown), higher than the final 13C polarization. Thus, deuteration of the sample, also suggested above regarding the acceleration of 13C spin diffusion, would improve the efficiency of dynamic 13C polarization using the electron spin in the photo-excited triplet state without being mediated by the 1H spins.

Since only a part of the broad EPR lines is excited in triplet DNP in the polycrystalline samples, the buildup of nuclear polarization takes place only in those crystallites that happen to be oriented in such a way that the resonant field is within the region of the adiabatic sweep. In other words, other crystallites with their resonant field outside the region of the adiabatic sweep would not experience nuclear polarization. The size of the crystallites is such that spin diffusion across the boundaries of adjacent crystallites is negligibly slow. It follows that those inert crystallites continue to be inert, and polarization speckle would develop in the sample. In this respect, the relatively compact EPR powder spectrum of pentacene compared to the broader NV powder spectrum is an advantage for the former, as the larger number of the crystallites can be excited. We estimate that the fraction of the crystallites contributing to the DNP process is 0.21 for microdiamonds and 0.565 for PBA. Considering this, the local 13C polarization within the DNP-active crystallites approaches 0.05 % for microdiamonds and 0.21 % for PBA, whereas the 13C NMR transitions in the other inert crystallites remain unpolarized.

4.5 Summary and prospects

We have demonstrated room-temperature optical 13C DNP using the NV centers and the pentacene molecules. The 13C polarization was successfully built up in microdiamonds to 0.01 % and in PBA to 0.12 %. The difference between these polarizations is due to a multitude of favoring/disfavoring factors.

In the diamond sample, the contact time tMW in the ISE sequence can be set much longer than that in PBA, thanks to the persistent nature of the ground triplet state of the NV center. However, the permanent nature of the paramagnetic NV system turns into a disadvantage for the next, spin-diffusion step: paramagnetic electrons of the NV centers cause a significant local field around the 13C spins in the vicinity, creating spin-diffusion barriers around them. Unfortunately, those 13C spins that are most likely to be directly polarized are inside the barrier, so that they are least likely to distribute the polarization away outside the barrier. The 13C spins just outside the barrier can both receive the electron polarization and undergo spin diffusion. However, their relatively long distance to the NV center leads to the low exchange probability ξNV≃0.0067. Moreover, many other paramagnetic electron defects in diamond, such as P1 centers, play the same detrimental role in the efficient spread of the nuclear polarization.

The PBA system is free from the problem of diffusion barriers due to the transient nature of the paramagnetic electrons in the triplet state, which decays to the “dark”, diamagnetic, ground state within the lifetime on the order of several microseconds. Even though the contact time is limited by the rather short lifetime, the efficiency of the direct, electron-to-13C polarization transfer was found to be higher than that for the case of the microdiamond sample.

However, the effect of abundant 1H spins reduces both the efficiency of the direct polarization transfer and of 13C spin diffusion. Even though one carefully optimizes the microwave amplitude in favor of the 13C spins, a considerable amount of electron polarization was found to leak to the 1H spins.

In addition, 1H−13C dipolar broadening (see the Supplement) lifts up the degeneracy in the energy levels of the 13C spin packets, reducing the 13C spin-diffusion rate and thereby the buildup of overall 13C polarization. The 13C spin-diffusion coefficient in PBA evaluated from the R dependence of the initial buildup rate in comparison with numerical simulation was found to be 2 orders of magnitude lower than the spin-diffusion coefficient for the naturally abundant 13C spins in diamond.

The optical bulk DNP experiments in NV-containing diamonds reported so far all used continuous-wave laser light sources, whereas we have demonstrated a pulsed optical excitation for the first time. The optical polarization of NV centers in microdiamonds was 16 %. The polarization of the electrons in the ground triplet state created through the pulsed photo-excitation can be lower than that obtained through continuous photo-excitation, in particular when the pulse length is shorter than the decay time of the meta-stable singlet state 1E to the triplet ground state 3A2, which is ca. 300–400 ns (Doherty et al.2013) (see Fig. 1). On the other hand, the pulsed excitation can be beneficial in preventing NV from photo-ionization to NV0 (S=1/2) (Loretz et al.2017; Chen et al.2017).

How could we improve optical DNP experiments, with these results in mind? Ajoy et al. (2018a) showed a smart way to strongly increase the number of excited crystallites, by choosing a low magnetic field, where the total broadening is <2D (a range where the Zeeman term is smaller than the ZFS term) and where a 13C polarization of 0.25 % (detected in a high-field NMR spectrometer) was achieved (Ajoy et al.2018a). A direct comparison of the results is practically impossible, due to differences in diamond samples (varying concentrations in 13C, NV, and P1 centers) (Parker et al.2019). It would be interesting to repeat DNP experiments of both samples in this low-field setup. At higher fields, a frequency-swept (Scheuer et al.2016; Can et al.2017) (in contrast to our field-swept) version combined with ultra-wide bandwidth MW chirp pulses (Segawa et al.2015) could increase the number of excited crystallites. Also, DNP schemes using the double-quantum transition could enhance the polarization by exciting more crystallites at different orientations (Chen et al.2015). From the diamond engineering side, highly pure crystals with a low concentration of paramagnetic defects (P1 centers and others) would be desirable, extending both electron and nuclear spin relaxation times. Diamonds fabricated using chemical vapor deposition (CVD) may fulfill this criterion, albeit at a very high price. We could not perform DNP experiments on our 100 nm nanodiamonds, since the sample amount was too low. The best DNP enhancement for 100 nm nanodiamonds reported in the literature is only ϵ=3 (compared with the thermal signal intensity at 7 T) (Ajoy et al.2020). The difficulty in achieving DNP in real nano-particles is directly linked to the drastic reduction of nuclear longitudinal relaxation times T1,C, finally arriving at values shorter than 1 s for 5 nm nanodiamonds (Casabianca et al.2011). Methods to eliminate paramagnetic surface defects, known as dangling bonds, as well as ways to transfer the polarization outside the diamond crystals (Abrams et al.2014), are yet to be experimentally demonstrated.

In this respect, pentacene-doped organic host crystals can highly compete; 1H spins in 150 nm size crystals of pentacene in p-terphenyl were optically polarized at room temperature in water, achieving a polarization of 0.086 % (Nishimura et al.2019). Optically hyperpolarized benzoic acid crystals doped with pentacene were dissolved and a DNP-enhanced liquid-state 1H NMR signal of benzoic acid was measured (Negoro et al.2018). These examples show impressive advantages on the material side of organic crystals with a controlled pentacene doping. However, organic crystals suffer from melting or damage when they are irradiated with high-power laser beams, whereas this is not a big problem for diamonds. The quest for new defects with ODMR signals, which has just started, will also lead to a variety of new candidates for optical triplet DNP (Segawa and Shames2020). The field of semiconductor spintronics presents a rich source of inspiration (Awschalom and Flatté2007). Hopefully, our comparative approach serves as a guideline to select or even engineer improved systems for room-temperature hyperpolarization of dilute/low-γ spins.

Table 3Comparison of parameters of the DNP experiment.

Download Print Version | Download XLSX

Data availability

Experimental data are available upon request from the corresponding authors.


The supplement related to this article is available online at:

Author contributions

KoMi, TFS, and MN designed the experiment and KoMi carried out the experiments under the supervision of TFS, AK, MK, and MN. SO, TO, HT, ST, FTKS, DT, RI, and MS prepared nanodiamonds containing NV centers. IO, HA, and NM prepared microdiamonds containing NV centers. YW and KT performed simulations estimating diffusion coefficients. AIS estimated the electron spin concentrations of the nanodiamonds from CW EPR spectra. KoMo estimated the P1 center concentration of the microdiamonds under the supervision of IO and NM. KoMi, TFS, KT, and MN analysed and interpreted the data and KoMi, TFS, KT, and MN wrote the manuscript with input from all the authors.

Competing interests

The authors declare that they have no conflict of interest.

Special issue statement

This article is part of the special issue “Geoffrey Bodenhausen Festschrift”. It is not associated with a conference.


This work is supported by MEXT Quantum Leap Flagship Program (MEXT Q-LEAP), PRESTO, and Osaka University's International Joint Research Promotion Program. Takuya F. Segawa is grateful to the Branco Weiss Fellowship – Society in Science, administered by the ETH Zurich, and Masahiro Shirakawa for hosting him as a Guest Research Associate at Kyoto University. Koichiro Miyanishi is supported by JSPS KAKENHI and the Program for Leading Graduate Schools: Interactive Materials Science Cadet Program. Kazuyuki Takeda is supported by JST CREST.

Financial support

This research has been supported by the Ministry of Education, Culture, Sports, Science and Technology (grant nos. JPMXS0120330644, JPMXS0118067634, and JPMXS0118067395), the Precursory Research for Embryonic Science and Technology (grant nos. JPMJPR1666, JPMJPR18G5, and PMJPR18G1), the Branco Weiss Fellowship – Society in Science, the Japan Society for the Promotion of Science (JSPS) KAKENHI (grant no. 19J10976), and the Core Research for Evolutional Science and Technology (JST CREST) (grant no. JPMJCR1873).

Review statement

This paper was edited by Konstantin Ivanov and reviewed by Jeffrey Reimer and one anonymous referee.


Abrams, D., Trusheim, M. E., Englund, D. R., Shattuck, M. D., and Meriles, C. A.: Dynamic nuclear spin polarization of liquids and gases in contact with nanostructured diamond, Nano Lett., 14, 2471–2478,, 2014. a

Acosta, V., Jarmola, A., Bauch, E., and Budker, D.: Optical properties of the nitrogen-vacancy singlet levels in diamond, Phys. Rev. B, 82, 201202,, 2010. a

Ajoy, A., Liu, K., Nazaryan, R., Lv, X., Zangara, P. R., Safvati, B., Wang, G., Arnold, D., Li, G., Lin, A., Raghavan, P., Druga, E., Dhomkar, S., Pagliero, D., Reimer, J. A., Suter, D., Meriles, C. A., and Pines, A.: Orientation-independent room temperature optical 13C hyperpolarization in powdered diamond, Sci. Adv., 4, eaar5492,, 2018a. a, b, c

Ajoy, A., Nazaryan, R., Liu, K., Lv, X., Safvati, B., Wang, G., Druga, E., Reimer, J. A., Suter, D., Ramanathan, C., Meriles, C. A., and Pines, A.: Enhanced dynamic nuclear polarization via swept microwave frequency combs, P. Natl. Acad. Sci. USA, 115, 10576–10581,, 2018b. a

Ajoy, A., Safvati, B., Nazaryan, R., Oon, J. T., Han, B., Raghavan, P., Nirodi, R., Aguilar, A., Liu, K., Cai, X., Lv, X., Druga, E., Ramanathan, C., Reimer, J. A., Meriles, C. A., Suter, D., and Pines, A.: Hyperpolarized relaxometry based nuclear T1 noise spectroscopy in diamond, Nat. Commun., 10, 5160,, 2019. a

Ajoy, A., Nazaryan, R., Druga, E., Liu, K., Aguilar, A., Han, B., Gierth, M., Oon, J. T., Safvati, B., Tsang, R., Walton, J. H., Suter, D., Meriles, C. A., Reimer, J. A., and Pines, A.: Room temperature “optical nanodiamond hyperpolarizer”: Physics, design, and operation, Rev. Sci. Instrum., 91, 023106,, 2020. a

Álvarez, G. A., Bretschneider, C. O., Fischer, R., London, P., Kanda, H., Onoda, S., Isoya, J., Gershoni, D., and Frydman, L.: Local and bulk 13C hyperpolarization in nitrogen-vacancy-centred diamonds at variable fields and orientations, Nat. Commun., 6, 8456,, 2015. a

Ardenkjær-Larsen, J. H., Fridlund, B., Gram, A., Hansson, G., Hansson, L., Lerche, M. H., Servin, R., Thaning, M., and Golman, K.: Increase in signal-to-noise ratio of > 10 000 times in liquid-state NMR, P. Natl. Acad. Sci. USA, 100, 10158–10163,, 2003. a

Awschalom, D. D. and Flatté, M. E.: Challenges for semiconductor spintronics, Nat. Phys., 3, 153–159,, 2007. a

Can, T. V., Weber, R. T., Walish, J. J., Swager, T. M., and Griffin, R. G.: Frequency-Swept Integrated Solid Effect, Angew. Chem. Int. Edit., 56, 6744–6748,, 2017. a

Casabianca, L. B., Shames, A. I., Panich, A. M., Shenderova, O., and Frydman, L.: Factors affecting DNP NMR in polycrystalline diamond samples, J. Phys. Chem. C, 115, 19041–19048,, 2011. a

Chen, J., Lourette, S., Rezai, K., Hoelzer, T., Lake, M., Nesladek, M., Bouchard, L.-S., Hemmer, P., and Budker, D.: Optical quenching and recovery of photoconductivity in single-crystal diamond, Appl. Phys. Lett., 110, 011108,, 2017. a

Chen, Q., Schwarz, I., Jelezko, F., Retzker, A., and Plenio, M. B.: Optical hyperpolarization of 13C nuclear spins in nanodiamond ensembles, Phys. Rev. B, 92, 184420,, 2015. a

Corval, A., Kryschi, C., Astilean, S., and Trommsdorff, H. P.: Resonant intersystem crossing in pentacene, J. Phys. Chem., 98, 7376–7381,, 1994. a

Deimling, M., Brunner, H., Dinse, K., Hausser, K., and Colpa, J.: Microwave-induced optical nuclear polarization (MI-ONP), J. Magn. Reson., 39, 185–202,, 1980. a

Doherty, M. W., Manson, N. B., Delaney, P., Jelezko, F., Wrachtrup, J., and Hollenberg, L. C.: The nitrogen-vacancy colour centre in diamond, Phys. Rep., 528, 1–45,, 2013. a, b

Fischer, R., Bretschneider, C. O., London, P., Budker, D., Gershoni, D., and Frydman, L.: Bulk Nuclear Polarization Enhanced at Room Temperature by Optical Pumping, Phys. Rev. Lett., 111, 057601,, 2013. a

Gruber, A., Dräbenstedt, A., Tietz, C., Fleury, L., Wrachtrup, J., and von Borczyskowski, C.: Scanning confocal optical microscopy and magnetic resonance on single defect centers, Science, 276, 2012–2014,, 1997. a

He, X.-F., Manson, N. B., and Fisk, P. T.: Paramagnetic resonance of photoexcited N-V defects in diamond, I. Level anticrossing in the 3A ground state, Phys. Rev. B, 47, 8809–8815,, 1993. a

Henstra, A., Lin, T.-S., Schmidt, J., and Wenckebach, W.: High dynamic nuclear polarization at room temperature, Chem. Phys. Lett., 165, 6–10,, 1990. a, b, c

Jelezko, F., Gaebel, T., Popa, I., Gruber, A., and Wrachtrup, J.: Observation of Coherent Oscillations in a Single Electron Spin, Phys. Rev. Lett., 92, 076401,, 2004. a

Kagawa, A., Murokawa, Y., Takeda, K., and Kitagawa, M.: Optimization of 1H spin density for dynamic nuclear polarization using photo-excited triplet electron spins, J. Magn. Reson., 197, 9–13,, 2009. a

Kagawa, A., Negoro, M., Ohba, R., Ichijo, N., Takamine, K., Nakamura, Y., Murata, T., Morita, Y., and Kitagawa, M.: Dynamic Nuclear Polarization using Photoexcited Triplet Electron Spins in Eutectic Mixtures, J. Phys. Chem. A, 122, 9670–9675,, 2018. a

Kagawa, A., Miyanishi, K., Ichijo, N., Negoro, M., Nakamura, Y., Enozawa, H., Murata, T., Morita, Y., and Kitagawa, M.: High-field NMR with dissolution triplet-DNP, J. Magn. Reson., 309, 106623,, 2019. a, b

King, J. P., Coles, P. J., and Reimer, J. A.: Optical polarization of 13C nuclei in diamond through nitrogen vacancy centers, Phys. Rev. B, 81, 073201,, 2010. a

King, J. P., Jeong, K., Vassiliou, C. C., Shin, C. S., Page, R. H., Avalos, C. E., Wang, H.-J., and Pines, A.: Room-temperature in situ nuclear spin hyperpolarization from optically pumped nitrogen vacancy centres in diamond, Nat. Commun., 6, 8965,, 2015. a

Köhler, J., Disselhorst, J., Donckers, M., Groenen, E., Schmidt, J., and Moerner, W.: Magnetic resonance of a single molecular spin, Nature, 363, 242–244,, 1993. a

Lesage, A., Lelli, M., Gajan, D., Caporini, M. A., Vitzthum, V., Miéville, P., Alauzun, J., Roussey, A., Thieuleux, C., Mehdi, A., Bodenhausen, G., Copéret, C., and Emsley, L.: Surface enhanced NMR spectroscopy by dynamic nuclear polarization, J. Am. Chem. Soc., 132, 15459–15461,, 2010. a

Loretz, M., Takahashi, H., Segawa, T. F., Boss, J. M., and Degen, C. L.: Optical hyperpolarization of nitrogen donor spins in bulk diamond, Phys. Rev. B, 95, 064413,, 2017. a

Lowe, I. J. and Gade, S.: Density-Matrix Derivation of the Spin-Diffusion Equation, Phys. Rev., 156, 817–825,, 1967. a

Maly, T., Debelouchina, G. T., Bajaj, V. S., Hu, K.-N., Joo, C.-G., Mak-Jurkauskas, M. L., Sirigiri, J. R., van der Wel, P. C., Herzfeld, J., Temkin, R. J., and Griffin, R. G.: Dynamic nuclear polarization at high magnetic fields, J. Chem. Phys., 128, 052211,, 2008. a

Mehring, M.: Principles of High Resolution NMR in Solids, Springer, Berlin, Heidelberg, Germany, 344 pp.,, 1983. a

Mindarava, Y. L., Blinder, R., Liu, Y., Scheuer, J., Lang, J., Agafonov, V. N., Davydov, V. A., Laube, C., Knolle, W., Abel, B., Naydenov, B., and Jelezko, F.: Synthesis and coherent properties of 13C enriched sub-micron diamond particles with nitrogen vacancy color centers, Carbon, 165, 395–403,, 2020. a

Negoro, M., Nakayama, K., Tateishi, K., Kagawa, A., Takeda, K., and Kitagawa, M.: 2H-decoupling-accelerated 1H spin diffusion in dynamic nuclear polarization with photoexcited triplet electrons, J. Chem. Phys., 133, 154504,, 2010. a

Negoro, M., Kagawa, A., Tateishi, K., Tanaka, Y., Yuasa, T., Takahashi, K., and Kitagawa, M.: Dissolution dynamic nuclear polarization at room temperature using photoexcited triplet electrons, J. Phys. Chem. A, 122, 4294–4297,, 2018. a, b

Nelson, S. J., Kurhanewicz, J., Vigneron, D. B., Larson, P. E., Harzstark, A. L., Ferrone, M., Van Criekinge, M., Chang, J. W., Bok, R., Park, I., Reed, G., Carvajal, L., Small, E. J., Munster, P., Weinberg, V. K., Ardenkjaer-Larsen, J. H., Chen, A. P., Hurd, R. E., Odegardstuen, L.-I., Robb, F. J., Tropp, J., and Murray, J. A.: Metabolic imaging of patients with prostate cancer using hyperpolarized [1−13C] pyruvate, Sci. Transl. Med., 5, 198ra108,, 2013. a

Nishimura, K., Kouno, H., Tateishi, K., Uesaka, T., Ideta, K., Kimizuka, N., and Yanai, N.: Triplet dynamic nuclear polarization of nanocrystals dispersed in water at room temperature, Phys. Chem. Chem. Phys., 21, 16408–16412,, 2019. a

Nishimura, K., Kouno, H., Kawashima, Y., Orihashi, K., Fujiwara, S., Tateishi, K., Uesaka, T., Kimizuka, N., and Yanai, N.: Materials chemistry of triplet dynamic nuclear polarization, Chem. Commun., 56, 7217–7232., 2020. a

Overhauser, A. W.: Polarization of Nuclei in Metals, Phys. Rev., 92, 411–415,, 1953. a

Parker, A. J., Jeong, K., Avalos, C. E., Hausmann, B. J., Vassiliou, C. C., Pines, A., and King, J. P.: Optically pumped dynamic nuclear hyperpolarization in 13C-enriched diamond, Phys. Rev. B, 100, 041203,, 2019. a

Peng, W. K., Takeda, K., and Kitagawa, M.: A new technique for cross polarization in solid-state NMR compatible with high spinning frequencies and high magnetic fields, Chem. Phys. Lett., 417, 58–62,, 2006. a

Rogers, L., Armstrong, S., Sellars, M., and Manson, N.: Infrared emission of the NV centre in diamond: Zeeman and uniaxial stress studies, New J. Phys., 10, 103024,, 2008. a

Scheuer, J., Schwartz, I., Chen, Q., Schulze-Sünninghausen, D., Carl, P., Höfer, P., Retzker, A., Sumiya, H., Isoya, J., Luy, B., Plenio, M. B., Naydenov, B., and Jelezko, F.: Optically induced dynamic nuclear spin polarisation in diamond, New J. Phys., 18, 013040,, 2016. a

Segawa, T. F. and Shames, A. I.: How to Identify, Attribute, and Quantify Triplet Defects in Ensembles of Small Nanoparticles, J. Phys. Chem. Lett., 11, 7438–7442,, 2020. a

Segawa, T. F., Doll, A., Pribitzer, S., and Jeschke, G.: Copper ESEEM and HYSCORE through ultra-wideband chirp EPR spectroscopy, J. Chem. Phys., 143, 044201,, 2015. a

Sosnovsky, D. V. and Ivanov, K. L.: Magnetic field dependence of triplet-state ONP: theoretical analysis in terms of level anti-crossings, Mol. Phys., 117, 2740–2755,, 2019. a

Stehlik, D. and Vieth, H.-M.: Time evolution of electron-nuclear cross-polarization in radiofrequency induced optical nuclear spin polarization, in: Pulsed Magnetic Resonance: NMR, ESR, and Optics: A Recognition of Hahn, E. L., edited by: Bagguley, D. M. S., Oxford University Press, Oxford, UK, 446–477, 1992. a

Stoll, S. and Schweiger, A.: EasySpin, a comprehensive software package for spectral simulation and analysis in EPR, J. Magn. Reson., 178, 42–55,, 2006. a

Takeda, K.: Triplet State Dynamic Nuclear Polarization: Basics, Concepts, Methods, VDM Verlag Dr. Mueller, Saarbrücken, Germany, 172 pp., 2009. a, b, c

Takeda, K., Takegoshi, K., and Terao, T.: Dynamic nuclear polarization by photoexcited-triplet electron spins in polycrystalline samples, Chem. Phys. Lett., 345, 166–170,, 2001.  a

Tateishi, K., Negoro, M., Nishida, S., Kagawa, A., Morita, Y., and Kitagawa, M.: Room temperature hyperpolarization of nuclear spins in bulk, P. Natl. Acad. Sci. USA, 111, 7527–7530,, 2014. a, b

Tateishi, K., Negoro, M., Nonaka, H., Kagawa, A., Sando, S., Wada, S., Kitagawa, M., and Uesaka, T.: Dynamic nuclear polarization with photo-excited triplet electrons using 6,13diphenylpentacene, Phys. Chem. Chem. Phys., 21, 19737–19741,, 2019. a

Terada, D., Segawa, T. F., Shames, A. I., Onoda, S., Ohshima, T., Osawa, E., Igarashi, R., and Shirakawa, M.: Monodisperse Five-Nanometer-Sized Detonation Nanodiamonds Enriched in Nitrogen-Vacancy Centers, ACS Nano, 13, 6461–6468,, 2019. a

Thiering, G. and Gali, A.: Theory of the optical spin-polarization loop of the nitrogen-vacancy center in diamond, Phys. Rev. B, 98, 085207,, 2018. a

van Kesteren, H. W., Wenckebach, W. T., Schmidt, J., and Poulis, N.: Dynamic nuclear polarization of proton spins via photoexcited triplet states: the system phenanthrene in fluorene, Chem. Phys. Lett., 89, 67–70,, 1982. a

van Kesteren, H. W., Wenckebach, W. T., and Schmidt, J.: Production of High, Long-Lasting, Dynamic Proton Polarization by Way of Photoexcited Triplet States, Phys. Rev. Lett., 55, 1642–1644,, 1985. a

Wenckebach, W. T.: Essentials of dynamic nuclear polarization, Spindrift Publications, Burgh-Haamstede, The Netherlands, 296 pp., 2016. a

Wrachtrup, J., von Borczyskowski, C., Bernard, J., Orrit, M., and Brown, R.: Optical detection of magnetic resonance in a single molecule, Nature, 363, 244–245,, 1993. a

Yap, Y. S., Tabuchi, Y., Negoro, M., Kagawa, A., and Kitagawa, M.: A Ku band pulsed electron paramagnetic resonance spectrometer using an arbitrary waveform generator for quantum control experiments at millikelvin temperatures, Rev. Sci. Instrum., 86, 063110,, 2015. a

Yu, H.-L., Lin, T.-S., Weissman, S., and Sloop, D. J.: Time resolved studies of pentacene triplets by electron spin echo spectroscopy, J. Chem. Phys., 80, 102–107,, 1984. a, b, c

Short summary
We show 13C spin hyperpolarization using laser irradiation at room temperature in two different systems: 1. pentacene doped in organic crystals of benzoic acid and 2. nitrogen-vacancy (NV) color centers in microdiamonds. The 13C NMR signal was enhanced by a factor of > 3000 using pentacene and > 300 using NV centers as polarization sources. We analyze the pros and cons of these two systems and discuss their prospects for room-temperature enhancement of NMR signals.