In this work, we demonstrate the use of frequency-tunable superconducting NbTiN coplanar waveguide microresonators for multi-frequency pulsed electron spin resonance (ESR) experiments. By applying a bias current to the center pin, the resonance frequency (∼7.6 GHz) can be continuously tuned by as much as 95 MHz in 270 ns without a change in the quality factor of 3000 at 2 K. We demonstrate the ESR performance of our resonators by measuring donor spin ensembles in silicon and show that adiabatic pulses can be used to overcome magnetic field inhomogeneities and microwave power limitations due to the applied bias current. We take advantage of the rapid tunability of these resonators to manipulate both phosphorus and arsenic spins in a single pulse sequence, demonstrating pulsed double electron-electron resonance. Our NbTiN resonator design is useful for multi-frequency pulsed ESR and should also have applications in experiments where spin ensembles are used as quantum memories.

Superconducting microresonators have dramatically enhanced the detection sensitivity of conventional electron spin resonance (ESR). Recently, single-shot detection of 107 spins has been demonstrated at 2 K using coplanar waveguide (CPW) resonators,1,2 and further improvements in the readout at mK temperatures have enabled detection of 103 spins3,4 using lumped-element resonators. As a result, these resonators are routinely used for the manipulation of spin ensembles and have also been proposed as buses for coherently coupling superconducting qubits with spins.5,6 The resonance frequencies of CPW and lumped-element resonators are typically fixed by fabrication, and conventional methods used to tune the resonance frequency post-fabrication using superconducting quantum interference devices are incompatible with the high magnetic fields that are typically necessary for X-band ESR. A method that enables frequency tunability of superconducting microresonators at high magnetic fields is thus desirable for multi-frequency ESR. The ability to tune the resonance frequency on-demand also allows greater control of spin dynamics7,8 as well as the coupled spin-cavity dynamics in the strong-coupling regime where spin ensembles are used as quantum memories.6,9 In this work, we fabricate frequency-tunable CPW resonators that are compatible with high magnetic fields following recently developed techniques.10–12 We show multi-frequency ESR with 31P donors in 28Si, finding no effect on spin coherence when the resonators are operated at different frequencies. We show that the resonance frequency can be tuned in no more than ∼270 ns. As a practical application, we take advantage of the rapid tunability of our resonators to manipulate 31P and 75As donor spin ensembles that are 33 MHz apart in a single pulse sequence.

The resonators are fabricated from a 20 nm NbTiN thin film (TC = 13.8 K) on a c-axis sapphire substrate using optical lithography followed by reactive-ion etching with an SF6/Ar plasma. Devices are wirebonded to a copper printed circuit board equipped with microwave connectors and attached to a rotatable sample holder which enables in situ alignment of the device with its surface parallel to the externally applied magnetic field.

An optical micrograph of a device is shown in Fig. 1(a). The structure is a photonic bandgap (PBG) resonator with two Bragg mirrors on either side of a cavity which defines the resonance frequency. A key feature of the PBG resonator is that the center pin is continuous throughout the device, enabling the application of bias currents. Further details of the PBG resonator design have been discussed previously.2,13 The Bragg mirrors are constructed using four periods of stepped-impedance waveguides, where a period is defined as a length of low-impedance (35 Ω) followed by a length of high-impedance (137 Ω) transmission line. The top inset of Fig. 1(a) shows a magnified view near one of the impedance steps. Three devices are fabricated where the width of the center pin of the cavity is 1.5, 2.5, and 4 μm. Care is taken to prevent current crowding14,15 at both ends of the cavity by tapering from the impedance steps, as shown in the bottom inset of Fig. 1(a).

FIG. 1.

(a) Optical micrograph of a tunable superconducting coplanar photonic bandgap microresonator. The device is patterned from a 20 nm NbTiN thin film on a c-axis sapphire substrate. The photonic bandgap is defined by the two Bragg mirrors on either side of the cavity. The Bragg mirrors are implemented using stepped impedances (top inset) and the resonance frequency is determined by the length of the cavity. The stepped impedances are tapered at both ends of the cavity to avoid current crowding (bottom inset). The center pin is continuous throughout the device from the input port to the output port. Three such devices are investigated with cavity center pin widths of 1.5 μm, 2.5 μm, and 4 μm. (b) The resonance frequency of the 4 μm device shifts by δf = 52 MHz for a bias current of 5 mA, with no change in the quality factor of 3000. (c) Dependence of δf on bias current for the three different devices. The maximum value of δf before the bias current exceeds the critical current is seen to be 95 MHz, 78 MHz, and 52 MHz for the 1.5 μm, 2.5 μm, and 4 μm devices, respectively. The fits (red lines) are generated using Eq. (1).

FIG. 1.

(a) Optical micrograph of a tunable superconducting coplanar photonic bandgap microresonator. The device is patterned from a 20 nm NbTiN thin film on a c-axis sapphire substrate. The photonic bandgap is defined by the two Bragg mirrors on either side of the cavity. The Bragg mirrors are implemented using stepped impedances (top inset) and the resonance frequency is determined by the length of the cavity. The stepped impedances are tapered at both ends of the cavity to avoid current crowding (bottom inset). The center pin is continuous throughout the device from the input port to the output port. Three such devices are investigated with cavity center pin widths of 1.5 μm, 2.5 μm, and 4 μm. (b) The resonance frequency of the 4 μm device shifts by δf = 52 MHz for a bias current of 5 mA, with no change in the quality factor of 3000. (c) Dependence of δf on bias current for the three different devices. The maximum value of δf before the bias current exceeds the critical current is seen to be 95 MHz, 78 MHz, and 52 MHz for the 1.5 μm, 2.5 μm, and 4 μm devices, respectively. The fits (red lines) are generated using Eq. (1).

Close modal

Microwave transmission through each device is monitored at a temperature of 2 K while a bias current is applied to the center pin. With no bias current, the resonance frequency of the 4 μm device is 7636.6 MHz [Fig. 1(b)] with a temperature-limited loaded quality factor of ∼3000 and a coupling coefficient of 0.6. With 5 mA of current applied to this device, the resonance frequency shifts to 7584.3 MHz, indicating a frequency shift, δf = 52 MHz. The device maintains its quality factor as the bias current is increased from 0 mA to 5 mA. The critical current of the device is measured to be 5.014 mA, in agreement with previously reported values for similar NbTiN wires where current crowding was minimized.14 

Figure 1(c) shows the dependence of δf as a function of the applied bias current for all three devices. The 1.5 μm device shows the largest change in resonance frequency, δf = 95 MHz, while the 2.5 μm device shows a change of δf = 78 MHz before the applied current exceeds the critical current of the superconductor.

The resonance frequency shift follows a well-known dependence on the applied bias current.10,16 The total inductance of our transmission lines is composed of both geometric and kinetic inductance. The application of a bias current, i, modulates the kinetic inductance, resulting in a change from the zero-current resonance frequency, f0, according to the expression

δ f ( i ) f 0 = [ ( i I 2 * ) 2 + ( i I 4 * ) 4 ] ,
(1)

where we note that the response is always negative, independent of the sign of the current.16 Here, I 2 * and I 4 * are parameters that set the scale of the current non-linearity of the kinetic inductance. We fit the measured values of δf for our three devices and extract I 2 * = { 28.9 , 50.1 , 62.5 } mA and I 4 * = { 13.9 , 26.3 , 35.7 } mA for the { 1.5 , 2.5 , 4 } μ m devices, respectively. The 1.5 μm device shows the largest response to applied current since it has the largest kinetic inductance fraction of the three devices. While the quartic component has typically been small in previous reports,12,17 we find that it is necessary when biasing near the critical current. Such dependence has been reported in devices with large kinetic inductance fractions.11,18

Next, we evaluate the ESR performance of our devices using a sample consisting of a 2 μm epitaxial layer of isotopically enriched 28Si grown on high-resistivity p-type Si. The epi-layer is bulk-doped with 31P donors to a concentration of 5 × 1015/cm3 and implanted with 75As donors at a nominal depth of 200 nm to a concentration of 5 × 1016 donors/cm3 as discussed in Ref. 2. We clip the sample on top of the 4 μm resonator with a phosphor bronze spring. We choose this resonator over the 1.5 μm and 2.5 μm devices because of better expected homogeneity of the microwave magnetic field seen by donor spins in the 2 μm epi-layer.

The device and sample are cooled to 2 K in a magnetic field of B0 = 275 mT oriented parallel to the plane of the superconductor. This orientation maintains the loaded quality factor of 3000 by avoiding dissipation due to the motion of trapped flux vortices in the superconductor under the influence of microwave currents.19,20 We measure echo-detected field sweeps of 31P donors with and without bias currents as shown in Fig. 2. In these measurements, the intensity of a Hahn echo [ π / 2 ( + x ) τ π ( + y ) τ echo ] is monitored with τ = 60 μs as the external magnetic field, B0, is swept. Figure 2(a) shows the results of field sweeps measured without bias current, where the mI = –1/2 31P donor hyperfine line appears at 274.78 mT. We repeated this experiment using both rectangular (π/2-pulse = 200 ns, π-pulse = 400 ns) and adiabatic pulses [BIR4-WURST20 (Ref. 1), 10 μs, chirp frequency ±2 MHz] with microwave powers of −29 dBm and −32 dBm, respectively, at the input port of the resonator. Adiabatic pulses have previously been shown to compensate for rotation errors arising from microwave magnetic field inhomogeneities in the bulk-doped sample.1 As shown in Fig. 2(a), the signals obtained using rectangular and adiabatic pulses are comparable, indicating good homogeneity of the microwave magnetic field in the 2 μm epi-layer.

FIG. 2.

Echo-detected field sweeps of the mI = –1/2 ESR transition of 31P donor spins in 28Si with 0 mA (a) or 4 mA (b) of bias current applied to the center pin of the resonator. Rectangular and adiabatic pulses are tested in both cases, and are found to produce similar signal intensities when no current flows through the center pin. On the other hand, rectangular pulses show reduced signal intensity when current flows through the center pin. The nonlinearity of the superconductor limits the maximum microwave power that can be applied at the input port of the resonator in the presence of a bias current. The reduced echo intensity for rectangular pulses is due to a combination of the reduced microwave power and magnetic field inhomogeneities introduced by the applied bias current. Adiabatic pulses are not power-limited and nearly recover the full signal intensity. (c) Coherence times are found to be the same with and without bias current, indicating that the bias current noise does not introduce additional decoherence to the spins at timescales shorter than ∼1 ms.

FIG. 2.

Echo-detected field sweeps of the mI = –1/2 ESR transition of 31P donor spins in 28Si with 0 mA (a) or 4 mA (b) of bias current applied to the center pin of the resonator. Rectangular and adiabatic pulses are tested in both cases, and are found to produce similar signal intensities when no current flows through the center pin. On the other hand, rectangular pulses show reduced signal intensity when current flows through the center pin. The nonlinearity of the superconductor limits the maximum microwave power that can be applied at the input port of the resonator in the presence of a bias current. The reduced echo intensity for rectangular pulses is due to a combination of the reduced microwave power and magnetic field inhomogeneities introduced by the applied bias current. Adiabatic pulses are not power-limited and nearly recover the full signal intensity. (c) Coherence times are found to be the same with and without bias current, indicating that the bias current noise does not introduce additional decoherence to the spins at timescales shorter than ∼1 ms.

Close modal

We performed the field sweep experiment with a bias current of i = 4 mA applied to the center pin of the resonator throughout the entire pulse sequence. The 31P signal is shifted to 273.72 mT in agreement with the resonance frequency change of 33 MHz. The center pin of the device is aligned with the external magnetic field, B0, such that the inhomogeneous magnetic field generated by the applied bias current, Bi, mostly adds to B0 in quadrature. However, a small misalignment (B0 in the plane of NbTiN but misaligned with the long axis of the resonator) leads to a component of Bi along B0 that broadens the observed ESR transition for spins directly above the center pin. Such a misalignment does not affect spins above the gap between the center pin and the ground plane where Bi remains orthogonal to B0 with the misalignment. From simulations that take into account the full distribution of Bi and the microwave magnetic field, B1, we estimate that a broadening comparable to the ESR linewidth can be caused by a misalignment of ∼4.7°. Additionally, we observed that microwave powers larger than −32 dBm affected the resonance lineshape with a bias current of 4 mA due to the non-linearity of the superconductor.21,22 The reduced signal intensity for rectangular pulses is due to a combination of reduced microwave power and broadening due to misalignment. By using adiabatic pulses which are not power-limited at −32 dBm, we nearly recover the full signal intensity as shown in Fig. 2(b).

In order to evaluate the effect of the applied bias currents on the coherence of the spin ensemble, we measure T2 with and without bias currents as shown in Fig. 2(c). In both cases, we use adiabatic pulses and measure the intensity of the Hahn echo as the interpulse delay, τ, is increased. We find T2 = 448 ± 5.2 μs without bias current and T2 = 450± 9.1 μs with 4.9 mA bias current. The agreement of these two values indicates that there is no effect on spin coherence due to noise in the applied bias current at timescales shorter than ∼1 ms. It is conceivable that the stability of the current source will become important at longer timescales.

In Fig. 3, we measure how quickly we can tune the resonator to a new frequency by applying a current pulse to the center pin and monitoring microwave transmission through the device. As the resonance frequency shifts in response to the current pulse at t = 0, it coincides with the frequency of the probe tone at a later time, t = T, enabling transmission through the resonator and resulting in a peak in the measurement at time T. We show the results of this measurement for a 3.9 mA current pulse in which the resonator reaches the target δf = 31.2 MHz in 270 ns. The inset compiles all measured tuning times for the 4 μm and 1.5 μm devices showing no frequency dependence with an average delay of 270 ns. In our experiments, we used resistive power combiners and DC blocks to combine the microwaves with the bias currents. The measured tuning times are likely limited by the time constant of our current-biasing circuit rather than the response of the superconductor.

FIG. 3.

Dynamics of the resonance frequency change in response to a 3.9 mA current pulse at time t = 0 measured by monitoring microwave transmission through the device. The inset shows tuning times for various values of the target δf for the 4 μm and 1.5 μm devices. The measured tuning times are independent of the frequency change.

FIG. 3.

Dynamics of the resonance frequency change in response to a 3.9 mA current pulse at time t = 0 measured by monitoring microwave transmission through the device. The inset shows tuning times for various values of the target δf for the 4 μm and 1.5 μm devices. The measured tuning times are independent of the frequency change.

Close modal

We demonstrate the use of this rapid tunability in a double electron-electron resonance (DEER) experiment.23,24 Here, we monitor the Hahn echo intensity from 31P donors while also applying a π pulse to 75As donors (hereafter referred to as πAs). The pulse sequence is shown in the inset of Fig. 4. The effect of πAs is to flip 75As donors such that the local dipolar field which is seen by each 31P donor spin in the first free evolution period (τ = 34 μs) is different from that of the second. This results in additional dipole-induced dephasing of the 31P donor spin ensemble due to the random distribution of 75As donors; the effect is maximized when πAs is applied close to the 31P donor π-pulse.25 

FIG. 4.

Double electron-electron resonance with 31P and 75As spin ensembles. The Hahn echo from 31P spins is monitored while a π-pulse (labeled πAs) is applied to the 75As spins. The position t of this pulse is swept and its effect on the 31P spin echo is measured. The two spin transitions are 33 MHz apart. Details about the pulse sequence are described in the main text. When πAs is resonant with the 75As spins, dipole-induced dephasing of the 31P spins is enhanced and the echo intensity is reduced by ∼20%.

FIG. 4.

Double electron-electron resonance with 31P and 75As spin ensembles. The Hahn echo from 31P spins is monitored while a π-pulse (labeled πAs) is applied to the 75As spins. The position t of this pulse is swept and its effect on the 31P spin echo is measured. The two spin transitions are 33 MHz apart. Details about the pulse sequence are described in the main text. When πAs is resonant with the 75As spins, dipole-induced dephasing of the 31P spins is enhanced and the echo intensity is reduced by ∼20%.

Close modal

The πAs pulse is applied at time t after the π/2 pulse which tips the 31P donors. The minimum value of t is 6 μs, arising from timing delays that protect the detection circuitry of our spectrometer and a 1 μs wait time for the resonator to settle to the frequency of 75As donors before the application of πAs. The results are shown in Fig. 4, where πAs is either on or off resonance with 75As donors, keeping the rest of the experimental parameters constant. When πAs is resonant, a reduction in the Hahn echo intensity of ∼20% at t = 20 μs is observed due to additional dephasing in the 31P donors arising from changes in the local magnetic field environment when the 75As donor spins are flipped.

We note here that a typical fixed-frequency resonator used for DEER requires a low Q-factor of ∼230 in order to accommodate the two donor frequencies that are 33 MHz apart. This is disadvantageous because the use of a low-Q resonator requires high-power pulses to invert both donors, and also leads to poor detection sensitivity. In contrast, our tunable resonator maintains a Q-factor of 3000 at both donor frequencies. The use of a fixed-frequency resonator would lead to a reduction in the detected echo signal by a factor of ∼10 meaning an increase in averaging time by a factor of ∼100 to maintain a similar signal-to-noise ratio. Therefore, our tunable resonator enables high-sensitivity DEER with significantly reduced signal-averaging time.

We remark here on practical considerations when using these devices. (1) The large kinetic inductance non-linearity of NbTiN results in an undesirable limitation on the maximum microwave power that can be applied for spin manipulation. For our 4 μm and 2.5 μm devices, microwave powers at the input port of the resonator beyond –15 and –23 dBm (B1 = 0.2 and 0.09 μT), respectively, were seen to affect the resonance lineshape without bias currents. This behavior has been reported previously21,22 and is a signature of the strong non-linearity of the superconductor. (2) In our DEER experiments, we applied two bias current pulses—one before, and another after the π pulse on 31P donors—to cancel phase accumulation from the inhomogeneous magnetic fields generated by the bias currents. We have found that bipolar current pulses can also be used to compensate for the unwanted phases since the negative currents reverse phase accumulation due to the positive currents. Bipolar pulses are advantageous because they only need to be applied on one side of the π pulse. (3) Pulsed DEER experiments for biochemical samples frequently must flip broad lines.26 Recent experiments have used chirped adiabatic pulses,27 but very low-Q resonators are then required, resulting in poor sensitivity. The tunable high-Q resonator frequency could be varied during the pulse to follow the chirped pulse frequency, thus allowing broad lines to be inverted without compromising the sensitivity of the experiment.

Our frequency-tunable devices may also be useful in experiments where a spin ensemble that is strongly coupled to the resonator is used as a quantum memory.6 This enables storage of quantum information for timescales of order T 2 * of the spin ensemble28 and can, in principle, be extended to T2, which can exceed seconds29 in silicon, by the use of refocusing π pulses. However, the application of a π pulse also inverts the spin ensemble into an unstable state that can emit the stored photons into the cavity by superradiance.9,28,30 In order to avoid uncontrolled emission of photons into the cavity, the resonator can be detuned away from the ESR frequency after the π pulse.30 In this work, we have demonstrated the desired frequency tunability. Further work is needed to elucidate the effect of magnetic field inhomogeneities introduced by the bias currents on the coupled spin-cavity dynamics.

In summary, we have demonstrated frequency-tunable ESR microresonators fabricated from NbTiN by taking advantage of the non-linear kinetic inductance of the superconductor. We have shown the ability to rapidly tune the resonance frequency of these devices by as much as δf = 95 MHz within 270 ns while maintaining a quality factor of 3000 at 7.6 GHz in a magnetic field of 275 mT at 2 K. We evaluated the performance of these resonators for ESR, finding that adiabatic (BIR4-WURST20) pulses can be used to overcome magnetic field inhomogeneities and microwave power limitations due to the applied bias current. We find that the applied bias currents have no effect on spin coherence at timescales up to 1 ms. We take advantage of the rapid frequency tunability to demonstrate high-sensitivity double electron-electron resonance in which we address both 31P and 75As donor spin transitions that are separated by 33 MHz in a single pulse sequence. Our devices enable multi-frequency electron spin resonance with CPW microresonators and may also be used in a hybrid quantum computation platform where spin ensembles that are strongly coupled to our resonator serve as quantum memories.

We thank Dr. Robin Cantor for preparation of NbTiN thin films. Devices were fabricated in the Princeton Institute for the Science and Technology of Materials Micro/Nano Fabrication Laboratory and the Princeton University Quantum Device Nanofabrication Laboratory. Work at Princeton was supported by the NSF through the Princeton MRSEC (Grant No. DMR-01420541) and by the ARO (Grant No. W911NF-13-1-0179). Work at LBNL was performed under the auspices of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.

1.
A. J.
Sigillito
,
H.
Malissa
,
A. M.
Tyryshkin
,
H.
Riemann
,
N. V.
Abrosimov
,
P.
Becker
,
H.-J.
Pohl
,
M. L. W.
Thewalt
,
K. M.
Itoh
,
J. J. L.
Morton
,
A. A.
Houck
,
D. I.
Schuster
, and
S. A.
Lyon
,
Appl. Phys. Lett.
104
,
222407
(
2014
).
2.
A. J.
Sigillito
,
A. M.
Tyryshkin
,
T.
Schenkel
,
A. A.
Houck
, and
S. A.
Lyon
, e-print arXiv:1701.06650 [quant-ph].
3.
C.
Eichler
,
A. J.
Sigillito
,
S. A.
Lyon
, and
J. R.
Petta
,
Phys. Rev. Lett.
118
,
037701
(
2017
).
4.
A.
Bienfait
,
J. J.
Pla
,
Y.
Kubo
,
M.
Stern
,
X.
Zhou
,
C. C.
Lo
,
C. D.
Weis
,
T.
Schenkel
,
M. L. W.
Thewalt
,
D.
Vion
,
D.
Esteve
,
B.
Julsgaard
,
K.
Mølmer
,
J. J. L.
Morton
, and
P.
Bertet
,
Nat. Nanotechnol.
11
,
253
(
2016
).
5.
K. D.
Petersson
,
L. W.
McFaul
,
M. D.
Schroer
,
M.
Jung
,
J. M.
Taylor
,
A. A.
Houck
, and
J. R.
Petta
,
Nature
490
,
380
(
2012
).
6.
Y.
Kubo
,
C.
Grezes
,
A.
Dewes
,
T.
Umeda
,
J.
Isoya
,
H.
Sumiya
,
N.
Morishita
,
H.
Abe
,
S.
Onoda
,
T.
Ohshima
,
V.
Jacques
,
A.
Dréau
,
J.-F.
Roch
,
I.
Diniz
,
A.
Auffeves
,
D.
Vion
,
D.
Esteve
, and
P.
Bertet
,
Phys. Rev. Lett.
107
,
220501
(
2011
).
7.
H.
Kaupp
,
T.
Hümmer
,
M.
Mader
,
B.
Schlederer
,
J.
Benedikter
,
P.
Haeusser
,
H.-C.
Chang
,
H.
Fedder
,
T. W.
Hänsch
, and
D.
Hunger
,
Phys. Rev. Appl.
6
,
054010
(
2016
).
8.
A.
Bienfait
,
J. J.
Pla
,
Y.
Kubo
,
X.
Zhou
,
M.
Stern
,
C. C.
Lo
,
C. D.
Weis
,
T.
Schenkel
,
D.
Vion
,
D.
Esteve
,
J. J. L.
Morton
, and
P.
Bertet
,
Nature
531
,
74
(
2016b
).
9.
C.
Grezes
,
B.
Julsgaard
,
Y.
Kubo
,
M.
Stern
,
T.
Umeda
,
J.
Isoya
,
H.
Sumiya
,
H.
Abe
,
S.
Onoda
,
T.
Ohshima
,
V.
Jacques
,
J.
Esteve
,
D.
Vion
,
D.
Esteve
,
K.
Mølmer
, and
P.
Bertet
,
Phys. Rev. X
4
,
021049
(
2014
).
10.
A. J.
Annunziata
,
D. F.
Santavicca
,
L.
Frunzio
,
G.
Catelani
,
M. J.
Rooks
,
A.
Frydman
, and
D. E.
Prober
,
Nanotechnology
21
,
445202
(
2010
).
11.
M. R.
Vissers
,
J.
Hubmayr
,
M.
Sandberg
,
S.
Chaudhuri
,
C.
Bockstiegel
, and
J.
Gao
,
Appl. Phys. Lett.
107
,
062601
(
2015
).
12.
A. A.
Adamyan
,
S. E.
Kubatkin
, and
A. V.
Danilov
,
Appl. Phys. Lett.
108
,
172601
(
2016
).
13.
Y.
Liu
and
A. A.
Houck
,
Nat. Phys.
13
,
48
(
2017
).
14.
H. L.
Hortensius
,
E. F. C.
Driessen
,
T. M.
Klapwijk
,
K. K.
Berggren
, and
J. R.
Clem
,
Appl. Phys. Lett.
100
,
182602
(
2012
).
15.
J. R.
Clem
and
K. K.
Berggren
,
Phys. Rev. B
84
,
174510
(
2011
).
16.
J.
Zmuidzinas
,
Annu. Rev. Condens. Matter Phys.
3
,
169
(
2012
).
17.
J.
Luomahaara
,
V.
Vesterinen
,
L.
Grönberg
, and
J.
Hassel
,
Nat. Commun.
5
,
4872
(
2014
).
18.
A.
Kher
,
P. K.
Day
,
B. H.
Eom
,
J.
Zmuidzinas
, and
H. G.
Leduc
,
J. Low Temp. Phys.
184
,
480
(
2016
).
19.
V.
Ginzburg
and
L.
Landau
,
Zh. Eksp. Teor. Fiz.
20
,
1064
(
1950
).
20.
H.
Malissa
,
D. I.
Schuster
,
A. M.
Tyryshkin
,
A. A.
Houck
, and
S. A.
Lyon
,
Rev. Sci. Instrum.
84
,
025116
(
2013
).
21.
B.
Abdo
,
E.
Segev
,
O.
Shtempluck
, and
E.
Buks
,
J. Phys.: Conf. Ser.
43
,
1346
(
2006
).
22.
L. J.
Swenson
,
P. K.
Day
,
B. H.
Eom
,
H. G.
Leduc
,
N.
Llombart
,
C. M.
McKenney
,
O.
Noroozian
, and
J.
Zmuidzinas
,
J. Appl. Phys.
113
,
104501
(
2013
).
23.
A.
Milov
,
K.
Salikhov
, and
M.
Shirov
,
Fiz. Tverd. Tela
23
,
975
(
1981
).
24.
A.
Milov
,
A.
Ponomarev
, and
Y. D.
Tsvetkov
,
Chem. Phys. Lett.
110
,
67
(
1984
).
25.
A. D.
Milov
,
A. G.
Maryasov
, and
Y. D.
Tsvetkov
,
Appl. Magn. Reson.
15
,
107
(
1998
).
26.
O.
Schiemann
and
T. F.
Prisner
,
Q. Rev. Biophys.
40
,
1
(
2007
).
27.
I.
Niemeyer
,
J. H.
Shim
,
J.
Zhang
,
D.
Suter
,
T.
Taniguchi
,
T.
Teraji
,
H.
Abe
,
S.
Onoda
,
T.
Yamamoto
,
T.
Ohshima
,
J.
Isoya
, and
F.
Jelezko
,
New J. Phys.
15
,
033027
(
2013
).
28.
B. C.
Rose
,
A. M.
Tyryshkin
,
H.
Riemann
,
N. V.
Abrosimov
,
P.
Becker
,
H.-J.
Pohl
,
M. L. W.
Thewalt
,
K. M.
Itoh
, and
S. A.
Lyon
, e-print arXiv:1702.00504 [quant-ph].
29.
A. M.
Tyryshkin
,
S.
Tojo
,
J. J. L.
Morton
,
H.
Riemann
,
N. V.
Abrosimov
,
P.
Becker
,
H.-J.
Pohl
,
T.
Schenkel
,
M. L. W.
Thewalt
,
K. M.
Itoh
, and
S. A.
Lyon
,
Nat. Mater.
11
,
143
(
2012
).
30.
B.
Julsgaard
and
K.
Mølmer
,
Phys. Rev. A
88
,
062324
(
2013
).