Journal of Magnetic Resonance 310 (2020) 106662
Contents lists available at ScienceDirect
Journal of Magnetic Resonance
journal homepage: www.elsevier.com/locate/jmr
Pulsed electron spin resonance spectroscopy in the Purcell regime
V. Ranjan a, S. Probst a, B. Albanese a, A. Doll b, O. Jacquot a, E. Flurin a, R. Heeres a, D. Vion a, ⇑ D. Esteve a, J.J.L. Morton c, P. Bertet a, a Quantronics Group, SPEC, CEA, CNRS, Université Paris-Saclay, CEA Saclay 91191 Gif-sur-Yvette Cedex, France b Laboratoire Nanomagnétisme et Oxydes, SPEC, CEA, CNRS, Université Paris-Saclay, CEA Saclay 91191 Gif-sur-Yvette Cedex, France c London Centre for Nanotechnology, University College London, London WC1H 0AH, United Kingdom article info abstract
Article history: In EPR, spin relaxation is typically governed by interactions with the lattice or other spins. However, it has Received 5 September 2019 recently been shown that given a sufficiently strong spin-resonator coupling and high resonator quality Revised 22 November 2019 factor, the spontaneous emission of microwave photons from the spins into the resonator can become the Accepted 25 November 2019 main relaxation mechanism, as predicted by Purcell. With increasing attention on the use of microres- Available online 4 December 2019 onators for EPR to achieve high spin-number sensitivity it is important to understand how this novel regime influences measured EPR signals, for example the amplitude and temporal shape of the spin- Keywords: echo. We study this regime theoretically and experimentally, using donor spins in silicon, under different Electron paramagnetic resonance conditions of spin-linewidth and coupling homogeneity. When the spin-resonator coupling is distributed Purcell effect Superconducting resonators inhomogeneously, we find that the effective spin-echo relaxation time measured in a saturation recovery Quantum limited amplifier sequence strongly depends on the parameters for the detection echo. When the spin linewidth is larger than the resonator bandwidth, the different Fourier components of the spin echo relax with different characteristic times – due to the role of the resonator in driving relaxation – which results in the temporal shape of the echo becoming dependent on the repetition time of the experiment. Ó 2019 Elsevier Inc. All rights reserved.
1. Introduction sequence, the spins are strongly out of equilibrium, so that before the next sequence can be started, a waiting time is needed for the
Pulsed magnetic resonance spectroscopy proceeds by applying longitudinal spin polarization Sz to relax back towards Sz0 by sequences of control pulses to an ensemble of electron or nuclear energy exchange between each spin and its environment in a char- spins via an electromagnetic resonator of frequency x0 (at micro- acteristic time T1. wave or radio frequency, respectively). Driven by these pulses, the In solids, the dominant energy exchange processes are usually spins undergo nutations on the Bloch sphere. Spins with identical spin-phonon or spin-spin interactions. For isotropic systems, T1 Larmor and Rabi frequencies (forming a spin packet) follow the depends only on global sample properties (temperature, concen- same trajectory. A prominent pulse sequence is the Hahn echo: a tration in magnetic species, ...); as a result, all spin packets con- first pulse imprints its phase coherence among all spin packets, tributing to the echo emission relax in the same way [2].In which quickly vanishes due to the inhomogeneity in Larmor pre- anisotropic systems, correlations may exist between T1 and the cession frequency. Coherence is restored by a second control pulse Larmor frequency of the spins [1,3]. Spin and spectral diffusion, applied after a delay s, which imposes a p phase shift to the spin as well as polarization transfer mechanisms, may also play a role packets leading to their collective rephasing after another delay and lead to non-exponential Sz relaxation [4,5]. s. This causes the buildup of a macroscopic oscillating magnetiza- Though much less common in EPR, spins can also, in principle, tion and the emission of a pulse known as the spin-echo, whose relax to equilibrium by exchanging energy with the radiation field. amplitude, shape, and time-dependence bear the desired informa- In free space, such radiative relaxation is negligibly slow for spins tion on spin characteristics and environment [1]. ( 1012 s at for electron spins at 9 GHz), however, it can be consid- The maximum spin-echo amplitude is governed by the equilib- erably accelerated by coupling the spin to a resonator. This is rium spin longitudinal magnetization Sz0. After each echo known as the Purcell effect, and arises due to the spatial and spec- tral confinement of the microwave field provided by the resonator [6–8]. The timescale for spin relaxation resulting from the Purcell ⇑ Corresponding author. effect is a function of the spin-resonator frequency detuning, d, E-mail address: [email protected] (P. Bertet). https://doi.org/10.1016/j.jmr.2019.106662 1090-7807/Ó 2019 Elsevier Inc. All rights reserved. 2 V. Ranjan et al. / Journal of Magnetic Resonance 310 (2020) 106662 the resonator energy damping rate j, and the spin-resonator cou- quantitative agreement with the simulation of the Bloch equations pling strength g: including explicitly the Purcell decay contribution. "# j 2d 2 T ¼ 1 þ : ð1Þ 2. Spin-echo in the Purcell regime: a simple model 1 4g2 j In this section we analyze a simple model which exemplifies the Here, g is equal to the half the Rabi frequency for a spin being driven consequences of Purcell relaxation by considering analytically- by a resonator field with an average energy of one photon, whose tractable limiting scenarios. amplitude we denote dB1, and j ¼ x0=Q is determined by the res- ¼ j= 2 onator frequency and total quality factor. At resonance, T1 4g , 2.1. System description and equations of motion which shows that Purcell relaxation rate is enhanced for resonators with high quality factor and small mode volume. In usual EPR spec- The system to be modelled is shown in Fig. 1. An ensemble of trometers, g=2p 10 3 Hz, and the Purcell time 108 s still much spins S ¼ 1=2 interact with the microwave field inside the res- longer than typical spin-lattice relaxation times. Recent experi- onator used for inductive detection, which is capacitively coupled ments have demonstrated that it is possible to considerably to a measurement line. The resonator is characterized by its fre- 2 3 increase g=2p up to 10 10 Hz by using the modes of supercon- quency x0, and total energy damping rate j. Energy is lost by leak- j ¼ x = ducting micro-resonators, in which the microwave field is confined age into the measurement line (with a rate c 0 Q c) and j ¼ x = j ¼ j þ j at the micron- scale around superconducting electrodes. The main internal losses (rate i 0 Q i), with c i. Control pulses motivation to do so is to enhance the spin detection sensitivity, at x0 are sent to the resonator input through the measurement 2 line, into which the subsequent spin echo signals are emitted then which scalesp withffiffiffiffiffiffi g and record spin-echo sensitivities of 102 103spin= Hz have been demonstrated in this way [9,10]. routed via a circulator towards the detection chain. ¼ ; ...; However, as seen from Eq. (1), radiative relaxation by the Purcell Each of the j 1 N spins is characterized by its Larmor fre- x d ¼ x x effect becomes unavoidably enhanced in these regimes up to the quency j (or equivalently the spin-cavity detuning j j 0) point where it may become the dominant relaxation channel for and its coupling to the resonator field gj. In the weak spin- j spins that have long enough spin-lattice relaxation times. This Pur- resonator coupling limit gj , the spin-resonator quantum corre- cell regime was reached recently for an ensemble of electron spins lations can be neglected. The dynamics is well described by equa- coupled to a superconducting micro-resonator at millikelvin tem- tions involving only the expectation value of the spin and peratures [11,12]. It was furthermore shown to enable electron spin resonator field operators [14], written in the frame rotating at x0 hyperpolarization by radiative cooling [13]. as Given the increasing use of micro-resonators in EPR which are 8 _ ðÞj ðÞj ðÞj ðÞj ðÞj > S ¼ d S þ 2g YS S =T capable of achieving substantial spin-resonator coupling, it is <> x j y j z x 2 important to understand any significant influences which the Pur- _ ðÞj ðÞj ðÞj ðÞj ðÞj Sy ¼ djS 2g XS S =T ð2Þ > x j z y 2 cell relaxation regime may impose on typical pulsed EPR measure- > : _ ðÞj ¼ ðÞj ðÞj ðÞj = ðÞj : ments. For example, one noticeable aspect of Eq. (1) is that Purcell Sz 2gjXSy 2gjYSx Sz Sz0 T1 relaxation is a resonant phenomenon, with a strong frequency j ðÞj dependence over a scale given by the resonator bandwidth Here, Sx;y;z is the expectation value of the corresponding dimension- ðÞ [11,12]. Another consequence of Eq. (1) is that T1 depends on the j ¼ r = ; r less spin j operator (with Sx;y;z x;y;z 2 x;y;z being the Pauli opera- position r of a given spin within the resonator mode, since g is pro- tors) and X; Y is the expectation value of the intra-resonator field portional to the field mode amplitude dB1ðÞr . Therefore, in the Pur- cell regime, spin packets with different detuning and Rabi frequency also have a different relaxation time T1, an unusual sit- uation in magnetic resonance. Because a spin echo is the sum of the contribution of all spin packets, spin-echo relaxation is expected to display a complex behavior, particularly when either the Larmor frequency or the Rabi frequency is inhomogeneously distributed, which happens with superconducting micro- resonators whose B1 field is often spatially-inhomogeneous. It is the purpose of this article to analyze the implications of these correlations between relaxation time, detuning and Rabi fre- quencies on the shape and time-dependence of the spin-echo in the Purcell regime. We first provide a simplified model that yields analytical results for the spin-echo amplitude and shape in the Pur- cell regime, and enables us to identify qualitatively novel effects. When the Rabi frequency is inhomogeneously distributed, the spin-echo amplitude is found to come back to equilibrium with an approximately exponential temporal dependence, but with a time constant that depends on the control pulse amplitude. When the Larmor frequency is inhomogeneously distributed over a fre- Fig. 1. A Sketch of electron-spin resonance spectroscopy in the Purcell regime. (a) quency range broader than the resonator bandwidth, the various An ensemble of spins is detected via a resonator of frequency x0. Each spin can d Fourier components of a spin echo relax with different time con- have a different coupling strength gj and frequency detuning j with the resonator, ðÞj stants, which also implies that the echo shape becomes dependent which leads to different individual Purcell relaxation time constants T1 . The j on the waiting time between consecutive sequences. The third sec- resonator has an internal energy decay rate i and is coupled to a drive and measurement line with an energy decay rate j . (b) A typical Hahn echo sequence tion describes the experimental setup and samples used to test c with rectangular pulses. Final echo contains contributions from different spin these effects, and the fourth section presents the measurements, packets. The relative amplitude of contributions depends on experimental param- their qualitative agreement with our simplified model, and their eters such as the pulse amplitude b and the repetition rate. V. Ranjan et al. / Journal of Magnetic Resonance 310 (2020) 106662 3 operators expressed in dimensionless units [14]. Here we use the 2.2. Hahn echo amplitude convention that X ¼ ðÞa þ ay =2 and Y ¼ iaðÞ ay =2; a (resp. ay) being the resonator annihilation (resp. creation) operator, which differs We now derive an approximate analytical expression for the from the one used in Refs. [14,15] by a factor 2. Phase coherence amplitude of a Hahn echo, based on Eqs. (1)–(4) and a number of ðÞj ðÞj simplifying assumptions. of spin j relaxes in a time T2 , and Sz relaxes towards the thermal ðÞ To simplify the discussion, we suppose that the control pulses equilibrium polarization S ¼ tanhðÞh x =2kT in a time T j . z0 0 s 1 generate a quasi-instantaneous intra-resonator field with a simple Because we are specifically interested in the impact of Purcell relax- rectangular time-dependence of duration dt. This can be achieved ation on spin-echo signals, we will in the following consider that if dt j 1, or by using shaped pulses that compensate for the the only decoherence mechanism is the Purcell relaxation, implying finite resonator response time [19,14,15]. We thus consider input ðÞ d 2 ðÞ ðÞ that T j ¼ j 1 þ 2 j and T j ¼ 2T j . b b 1 4g2 j 2 1 pulses on the in-phase X direction, with an amplitude X dur- j pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi b ¼ b ¼ = x ing dt while Y 0, related to the input power Pin as Pin h 0. Note that Eq. (2) are identical to the usual Blochqffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi equations, with pffiffiffiffiffi The corresponding intra-resonator field amplitude is a ¼ 2 jcb=j. the Rabi frequency given by 2g aðÞt , and aðÞ t XtðÞ2 þ YtðÞ2. Pur- j Under each controlqffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi pulse, spin j undergoes a Rabi nutation [1] cell relaxation appears as a T mechanism, and its particular fea- ÀÁ 1 a 2 þ d2 r a with a frequency 2gj j . We assume that d 2gj for tures arise from the dependence of T1 on the properties of a spin d packet. all spins so that the dependence of the Rabi frequency on j can q Using standard input-output theory [16], the resonator field in be neglected. Also, we assume that d is symmetric and is excited in its centre. the rotating frame at x0 obey Eqs. The Hahn-echo pulse sequence is shown in Fig. 1. A first pulse of 8 amplitude b=2 and duration dt is followed by a waiting time s, then > pffiffiffiffiffi XN b > _ j ðÞj by a second pulse of amplitude and same duration, and by a sec- > XtðÞ¼ j b ðÞ t XtðÞ g S <> c X 2 j y ond waiting time t. An echo is formed around t ¼ s because of the j¼1 ð3Þ refocusing of the spin packets (we assume dt s). We also use an > pffiffiffiffiffi XN > _ j ðÞj equivalent Hahn-echo sequence where the two pulses have the > YtðÞ¼ j b ðÞ t YtðÞþ g S ; : c Y 2 j x b = j¼1 same amplitude , but the first pulse duration is dt 2. The spin- echo originates solely from the y component of the magnetization, and the contribution of spin j can be shown [20,1] to be where b ; ðÞt are the in-phase and out of phase components of the ÀÁ X Y ðÞj ðÞ¼ s þ ðÞj 3 a d ðÞ s : ð Þ control field. We also need to compute the field leaking out of the Sy t Sz sin 2 gjdt cos j t 5 cavity, as it contains the spin free-induction-decay and echo signals. ðÞj where Sz is spin j longitudinal polarization at the time of the first The output field components XoutðÞt ; YoutðÞt are given by the input- output relations control pulse. The latter is not necessarily equal to Sz0 because the waiting time since the previous pulse sequence may not be suffi- ðÞ pffiffiffiffiffi ciently long. Importantly, because of Purcell relaxation, S j depends ðÞ¼ j ðÞ b ðÞ z Xout t cXt X t on g and d , which leads to novel effects as shown below. Note that pffiffiffiffiffi ð4Þ j j ðÞ¼ j ðÞ b ðÞ: Y out t cYt Y t we have, however, neglected the impact of spin relaxation during the Hahn echo sequence, because in most relevant scenarios The combination of Eqs. (1)–(4) is the theoretical framework s ðÞj that describes EPR spectroscopy in the Purcell regime. All the T1 . specific characteristics of the spin ensemble are provided by the To obtain simple expressions for the echo amplitude, we also q ðÞd consider that the resonator field dynamics adjusts adiabatically distributions of Larmor frequencyR d Rand coupling constant q ðÞ q ðÞ¼ dq ðÞ¼d to the spin operators, as would be the situation in the limit of g g , normalized such that dg g g d d 1. We also low resonator Q. Eqs. (3) and (4) then yield define rd as the standard deviation of qd, thus corresponding to pffiffiffiffi X j ðÞ the spin ensemble inhomogeneous linewidth. ðÞ¼ s þ c j ðÞs þ Xout t 2 j gjSy t Although this is not immediately apparent, we note that radia- j ð Þ pffiffiffiffi 6 tion damping is automatically included in the above equations, j RR ¼ 2 c gq ðÞg q ðÞd S ðÞg; d; t dddg; since they treat on an equal footing the intra-resonator field and j g d y the spin operators and thus incorporate all feedback effects of where we have taken the continuous limit in the last expression the radiation field on the spin dynamics [17,14]. This underlines ðÞ and defined S ðÞg; d; t as being equal to S j ðÞs þ t for a spin j having the distinction between radiative damping and Purcell relaxation, y y a coupling g ¼ g and detuning d ¼ d. Note also that we have which both are radiative effects but with different characteristics j j and impact on the spin dynamics. implicitly assumed that there is no correlation between the fre- quency of a given spin and its coupling constant g, an assumption Throughout this article, we make the extra simplifyingP hypoth- 2 that is not always verified [21]. esis that the spin-ensemble cooperativity C ¼ g =ðÞjrd verifies j This echo amplitude X ðÞt þ s depends in an intricate way on C 1. In this limit, the field generated by the spins is small com- out the pulse amplitude and on the shape of the inhomogeneous distri- pared to the intra-resonator field [18] so that radiation damping butions q and q . In the following we consider two limiting cases. can be entirely neglected. One can thus (1) compute the intra- g d resonator field ½ XtðÞ; YtðÞ using Eq. (3) with the last term neglected, (2) use it to solve the spin dynamics (Eq. (2)), and (3) 2.3. Narrow-line case compute the output field with Eqs. (3) and (4). This is the approach that is used to simulate numerically the spin-echo signals under Let us first assume that the spin ensemble has a linewidthÀÁ much 2 arbitrary control pulse sequences, and in the next sections to narrower than the resonator, i.e. rd j. Then, T1 j= 4g , derive approximate expressions for the echo amplitude. More implying that Sz and Sy do not depend any longer on d. Overall, details on the simulations can be found in the Appendix. the echo amplitude at t ¼ s becomes 4 V. Ranjan et al. / Journal of Magnetic Resonance 310 (2020) 106662 pjffiffiffiffiffi Z ðÞ¼s c q ðÞ 3 ðÞa : ð Þ sider situations encountered in recent experiments, we assume Xout 2 2 gSz g g sin 2 gdt dg 7 j that B1 is generated by a narrow wire deposited at the surface of a sample containing the spins. The spins may be distributed within To appreciate the impact of Purcell relaxation, we now focus on the a thin layer just below the sample surface (type B), or homoge- so-called saturation recovery sequence schematically depicted in neously in the bulk (type C). Assuming that dB ðÞ¼r l di=ðÞ2pr ; r Fig. 2c, which is commonly used to measure spin relaxation time. 1 0 being the spin-wire distance and di the quantum fluctuations of It consists of applying a saturating pulse at time t ¼ 0 (so that the ac current in the resonator, it is straightforward to see that S ðÞt ¼ 0 ¼ 0 for all spins), followed by a Hahn-echo sequence z q ðÞ¼g dðÞg g in type A, qðÞ¼g g =g2 in type B, and applied after a waiting time T at t ¼ T. Its maximum amplitude at g A B qðÞg ¼ g2 =g3 in type C. time t ¼ T þ 2s is denoted by XoutðÞT for simplicity in the following. C At the beginning of the Hahn-echo sequence, For completeness, it is interesting to note that the situation may ÀÁ be complicated by the presence of correlations between the spin T=T1 SzðÞ¼T Sz0 1 e : ð8Þ Larmor and Rabi frequency. This was reported in Refs. [9,11]. There, silicon donor spins were confined to a thin layer (100 nm) In the usual situation where spin relaxation is not Purcell- below the surface of a silicon sample, on top of which a thin-film limited or when it is governed by interactions with the lattice, resonator with a 5 lm-wide inductance was deposited. Because and where the system is moreover isotropic, all spins then relax of the mechanical strain exerted by the thin metallic film on the with the same time constant T . In Eq. (7), S can be factorized; 1 z silicon substrate [21] due to differential thermal contractions the echo amplitude is proportional to S ðÞT , and it follows the expo- z between the metal and the silicon, the spin hyperfine interaction nential dependence of S ðÞT , which enables to measure T . The z 1 becomes correlated with the lateral position relative to the wire, echo amplitude effectively measures the longitudinal polarization which also happens to be approximately correlated to B .Asa S ðÞT , which justifies its denomination as a detection echo; in partic- 1 z result, the Larmor and Rabi frequencies of the spins are correlated, ular, its parameters (pulse amplitude, duration, ...) have no impact and by properly choosing the biasing field B the system is much on the measured T . 0 1 better described by type A than type B. This approximate correla- In the Purcell regime, T is different for spin-packets with differ- 1 tion is valid only when the wire transverse dimensions are large ent couplings g; the detection echo amplitude X ðÞT has thus no out compared to the spin layer thickness, and was indeed no longer reason to be even exponential. Because of the strong g- found in Ref. [10] where the wire width was decreased to 500 nm. dependence of the integrand in Eq. (7) however, an approximate We will now compute the signal expected from a saturation exponential dependence is nevertheless recovered in a number of recovery sequence for three afore-mentioned types of spin distri- situations, but with an effective relaxation time that now depends bution, by inserting Eq. (8) into Eq. (7). on the parameters of the detection echo. In type A, we obtain straightforwardly We now make this reasoning explicit by considering three types pffiffiffiffiffi q j T of coupling constant distributions g shown in Fig. 2. The spins c 3 ðÞ ðÞ¼ ðÞa T1 gA : ð Þ Xout T 2 gA sin 2 gAdt Sz0 1 exp 9 may be located in an area where B1 is very homogeneous (type j A), or on the contrary very inhomogeneous (types B and C). To con- Because the coupling constant has a well-defined value gA, the spin relaxation time is identical for all measured spins and we are in the same situation as in usual magnetic resonance, where the detection echo amplitude relaxes exponentially with the Pur- ðÞ a cell relaxation time T1 gA , independently of . To deal with types B and C, it is useful to introduce the Rabi nutation angle w ¼ 2agdt. The T dependent echo amplitude can then be written as pffiffiffiffiffi Z hi j T c ðÞa;w ðÞ¼ T1 ðÞw w; ð Þ Xout T 2 j Sz0 1 exp f d 10
ðÞ¼a; w ja2dt2 ðÞ¼w wq 3 w=ðÞa ðÞw with T1 w2 and f w sin 2 dt . The function f indicates the relative contribution of spin packets to the total echo signal as a function of their Rabi angle w. As shown in Fig. 2 for types w p= B and C, it displays a maximum at a value 0 close to 2. In a crude ðÞ w dwðÞ w approximation, f 0 , so that pffiffiffiffiffi j T c ðÞa;w ðÞ¼ T1 0 : ð Þ Xout T 2 j Sz0 1 exp 11
Despite the broad coupling constant distribution, one thus recovers an approximate exponential dependence of the detection echo amplitude on the waiting time T. However, the effective measured relaxation time scales like a2, the square of the amplitude of the pulse used in the detection echo sequence. The physical interpreta- tion is straightforward: the detection echo signal mostly originates p= Fig. 2. Rabi angle selectivity in the narrow line case. (a) Types of spin implantation from the contribution of spins that undergo first a 2 pulse and in the substrate, type A: point distribution, type B: thin layer distribution, type C: then a p pulse during the refocusing step. When the pulse ampli- bulk distribution. Contours of rf magnetic field B1 created by the inductor are tude is varied, this amounts to selecting spins with different cou- sketched as dashed curves. (b) Calculated relative echo contribution f ðÞw from spin pling constants g, and therefore different relaxation times. packets with different Rabi angle w. (c) Saturation recovery of the echo amplitude To verify the validity of approximating f ðÞw by a d function, we calculated using the model. We have taken small Rabi angles 0 < w < p. Solid black ðÞa; w compute Eq. (10) numerically for types B and C. The results are curves show exponential recovery of polarization with a decay constant T1 0 , w ðÞw ðÞ where 0 is the value where f is maximum in panel (b). shown in Fig. 2(c) for the decay of the spin-echo amplitude AT V. Ranjan et al. / Journal of Magnetic Resonance 310 (2020) 106662 5
(open circles), compared to the single-exponential approximation icon at millikelvin temperature have an intrinsic relaxation time of Eq. (11) (solid lines, taking into account the different values of T1;int ¼ 1500 s, and are thus in the Purcell regime whenever w ¼ 0 for types B and C). The qualitative agreement demonstrates that T1 T1;int. This was shown in Ref. [11], where T1 1 s was reached Eq. (11) correctly captures the impact of Purcell relaxation on the when the bismuth donors were at resonance with a high-quality- effective relaxation time measured in a saturation-recovery factor superconducting micro-resonator. All the measurements sequence. reported in this article are also in this Purcell regime. We present data from three different devices. Each device is a 2.4. Narrow-coupling case silicon sample that was implanted with bismuth atoms close to its surface, and on which a discrete-element superconducting LC We now consider the case where the coupling constant is resonator was patterned. The resonators consist of an interdigi- single-valued (corresponding to type A in the previous section) tated capacitance shunted by a micron- or sub-micron-scale wire q ðÞ¼ dðÞ g g g gA . The pulse amplitude is chosen such that which plays the role of the inductance [Fig. 3 (a-c)]. The devices 2g adt ¼ p=2, so that the control pulses implement the ideal are mounted in a copper sample holder, and coupled capacitively A j Hahn-echo sequence. The Larmor frequency on the other hand is to a microwave antenna which determines the coupling rate c. broadly distributed, with rd j. In devices 1 and 2 [Fig. 3 (d-e)], the implantation depth is 100 In that limit, the Hahn echo amplitude Eq. (6) becomes nm, with a peak concentration of 8 1016 cm 3, the silicon sample pjffiffiffiffiffi Z is isotopically enriched in 28Si, and the resonator is made in alu- ðÞ¼s þ c q dðÞ s d: ð Þ l Xout t 2 j Sz d cos t d 12 minum. In device 3 [Fig. 3 (f)], the implantation depth is 1 m, with a smaller peak concentration of 1016 cm 3, the silicon is of Consider now that the echo sequences are repeated multiple times, natural isotopic abundance, and the resonator is made of niobium. with a waiting time T in-between two consecutive sequences, and Of particular importance for this work is the geometry of the res- let us assume that the spins are fully un-polarized at the immediate onator inductance, which strongly impacts the coupling constant end of a sequence, S ¼ 0. Then at the beginning of each echo, the q ðÞ l z distribution g g . It is 100 nm wide and 10 m long in device 1, longitudinal polarization Sz is given by Eq. (8). 500 nm wide and 100 lm long in device 2, and 2 lm wide and In a usual situation where T1 is not correlated with the spin Lar- mor frequency, Sz can be factored out of the integral in Eq. (12), and the echo temporal shape is simply given by the Fourier transform of the Larmor frequency distribution qdðÞd . In the Purcell regime however, we get that ÀÁ = ðÞd T T1 Sz ¼ Sz0 1 e ; ð13Þ so that pffiffiffiffiffi Z ÀÁ jc = ðÞd ðÞ¼s þ T T1 q ðÞd dðÞ s d; Xout t 2 j Sz0 1 e d cos t d ð14Þ whose Fourier transform is pffiffiffiffiffi ÀÁ jc = ðÞd ðÞ¼d; T T1 q ðÞd : ð Þ Xout T 2 j Sz0 1 e d 15