Observing hyperfine interactions of NV centers in diamond in an advanced quantum teaching lab

The negatively charged nitrogen-vacancy (NV$^-$) center in diamond is a model quantum system for university teaching labs due to its room-temperature compatibility and cost-effective operation. Based on the low-cost experimental setup that we have developed and described for the coherent control of the electronic spin (Sewani et al.), we introduce and explain here a number of more advanced experiments that probe the electron-nuclear interaction between the \nv electronic and the \NN~and \CC~nuclear spins. Optically-detected magnetic resonance (ODMR), Rabi oscillations, Ramsey fringe experiments, and Hahn echo sequences are implemented to demonstrate how the nuclear spins interact with the electron spins. Most experiments only require 15 minutes of measurement time and can, therefore, be completed within one teaching lab.


I. INTRODUCTION
The recent progress in the field of quantum technologies -quantum computing, quantum communications and quantum sensors -has led to considerable commercial interest from entities like the semiconductor industry, software companies and consulting firms. Quantum technologies are predicted to constitute a ∼$10 billion market within the next decade, 2 and the high demand for quantum-specialists, as reflected by the high number of current job advertisements, 3 corroborates this prediction. Universities need to cater to this need and provide their graduates with a thorough understanding of quantum mechanics and solid-state physics, including experimental experience in how to perform quantum measurements and interpret results. However, these skills are difficult to convey in a university teaching lab, due to the sensitivity of quantum systems to environmental disturbances, and the high cost of specialized measurement equipment.
As discussed amply in literature, the nitrogen vacancy (NV − ) centre in diamond possesses an electronic spin that can be initialized and readout optically. It displays microsecond coherence times at room temperature in ambient conditions. 4 This makes it a model system for highschool or university teaching labs to perform magnetometry [5][6][7] and coherent spin control experiments. 1,7-10 Based on the same low-cost experimental setup that we have developed in Ref. 1, in combination with an off-the-shelf, commercially available, high-quality chemical vapour deposition (CVD) diamond sample from Element Six, we expand the portfolio of experiments to include the characterization of electron-nuclear spin coupling dynamics.
This paper is organized as follows: In the proceeding section (Section II) we will give an overview of the learning outcomes that can be conveyed with the experiments in this paper. Section III will then introduce the NV − Hamiltonian with a focus on the terms describing the nuclear spin of the 14 N nucleus and its interaction with the electronic spin, and Section IV will introduce the Rabi formula that describes the time evolution of a spin during spin resonance. In Section V we will detail the changes made to the equipment with regards to Ref. 1, and in Section VI we will present the different experiments with a special emphasis on how the coupling of the electron spin to the 14 N nucleus manifests in the optically-detected magnetic resonance (ODMR) spectrum (Subsection VI B), Rabi oscillations (Subsection VI C), and Ramsey fringes (Subsection VI D). Finally, we show how the coupling to the bath of 13 C nuclei manifests in the Hahn echo decay (Subsection VI E).

II. LEARNING OUTCOMES
The experiments described below are a set of advanced experiments that build upon the basic ones described in Ref. 1. The basic experiments taught the students fundamental concepts such as the NV − center structure, the spin initialization and readout procedure, the features of the ODMR spectrum, development of pulse sequences, the rotating frame, magnetic resonance, two-axes control, the Bloch sphere, and dynamical decoupling. We will assume knowledge of these concepts for this manuscript, and build upon them for the students to: • Develop a detailed understanding of the system's Hamiltonian and the ability to map the observed experimental features to its properties.
• Comprehend the effect of experimental nonidealities like power-broadening and spectral detuning.
• Observe the interaction between the NV − electronic spin and the 14 N nuclear spin, mediated by the hyperfine interaction.
• Experience how the electron-nuclear interaction manifests in both the frequency domain and the time domain, and how those are linked by the Fourier transform.
The concepts demonstrated with these experiments will provide the students with an advanced understanding of quantum experiments and a detailed knowledge of electron-nuclear spin interactions. They are applicable to state-of-the-art research on NV − centers in diamond, [11][12][13] other types of color centers in diamond 14,15 and silicon carbide, 16 donor spins in silicon, [17][18][19]  These four orientations are in principle equivalent, but lead to different alignments of the NV − center axes with respect to an externally applied static or oscillating magnetic field (see Fig. 3). The natural abundance of 14 N is 99.6% and its nuclear spin number is I = 1, which means that the electron spin in almost all NV − centers is coupled to a nuclear spin. We therefore need to extend the ground state Hamiltonian of the NV − center to account for the coupled 14 N nuclear spin. The static Hamiltonian H 0 contains three main terms: where H S is the electron spin Hamiltonian, H I is the nuclear spin Hamiltonian, and H SI describes the interaction of the electron spin with the nuclear spin. The electron spin Hamiltonian is the corresponding term for the 14 N nuclear spin. P ≈ −5.01 MHz is the nuclear quadrupole moment of 14 N, I is the vector of nuclear spin matrices [I x , I y , I z ] for I = 1, and γ n = 3.077 MHz/T is the nuclear gyromagnetic ratio of 14 N. [21][22][23] Since γ n is roughly a factor 10, 000 smaller than γ e , the interaction between the nuclear spin and the external B 0 field is not observed in our experiments as it would require several 100 mT of field to be resolved. The final term in H 0 describes the interaction of the NV − electron spin with the 14 N nuclear spin: where A || = 2.14 MHz and A ⊥ = 2.7 MHz are the parallel and perpendicular components of the hyperfine tensor with respect to the NV − center axis. 21 In Fig. 1(a) we show a schematic of the energy levels of the system (not to scale) according to the Hamiltonian H 0 in Eq. 1. We start with the bare ground state energy level |g of the NV − center on the left side. Turning on the zero-field splitting D then lifts the degeneracy between the different S = 1 levels, by raising the m s = ±1 levels by D = 2.87 GHz above the m s = 0 level. Next, we include the electron Zeeman interaction as given by γ e B 0 · S. Assuming that B 0 is applied along the zaxis, the Zeeman interaction raises the m s = +1 level by 2γ e B z above the m s = −1 level. Finally, we include the 14 N nuclear spin. The quadrupole component P lifts the degeneracy between the different I = 1 nuclear spin levels by raising the m I = 0 level by 5.01 MHz above the m I = ±1 levels (see Eq. 3), and the hyperfine interaction between electron spin and nuclear spin further splits the m I = +1 and m I = −1 levels by |2A || | for the m s = ±1 electron spin states (see Eq. 4). The nuclear spin Zeeman interaction as defined by γ n B 0 · I in Eq. 3 is not included in this schematic as it is too small to be observed in our experiments.

IV. RABI FORMULA
One of the learning outcomes is the understanding of experimental non-idealities. A very common non-ideality is a driving field that is detuned in frequency from the spin resonance transition and that is strong in power. This will lead to non-perfect rotations of the spin, and the power-broadening of the ODMR resonance lines which prevents the observation of fine features of the spectrum at high driving powers. The effect of both can be easily understood via the Rabi formula, that describes the probability of flipping a spin (or any two-level system) as a function of pulse time when a coherent driving pulse is applied: 24 where t is the length of the driving pulse, ∆ = ν MW −ν res is the frequency detuning between the frequency of the driving field ν MW and the resonance frequency of the spin transition ν res , and Ω R is the rotation frequency of the spin while it is driven -the so-called Rabi frequency. 25 In general Ω R = 1 √ 2 γ e B 1 for a spin-1 system, where B 1 is the amplitude of the oscillating magnetic field provided Zero field splitting Zeeman interaction 14   by an antenna. 26 Note, that this is different to a spin-1/2 system for which Ω R = 1 2 γ e B 1 . Eq. 5 can be considered in two parts. The first part is the envelope, Ω 2 R Ω 2 R +∆ 2 that describes the amplitude of the Rabi oscillations. If the amplitude is plotted as a function of the frequency detuning ∆, the result is a Lorentzian lineshape that has a half-width at halfmaximum (HWHM) of ∆ HWHM = Ω R . This means that for a detuning of ∆ = Ω R the spin can still be rotated with 50% amplitude, clearly showing how a high Rabi frequency caused by a large B 1 amplitude leads to a large linewidth in the ODMR spectrum.
The second part of Eq. 5 is the time-dependent part sin 2 1 2 Ω 2 R + ∆ 2 2π t that describes the actual time evolution of the spin during the Rabi oscillations. The frequency of the Rabi oscillations is given by Ω eff R = Ω 2 R + ∆ 2 27 , which reduces to Ω eff R = Ω R for the resonant case with ∆ = 0. A profound consequence on experiments being performed at a frequency detuning ∆ = 0 is that the Rabi oscillations are faster and have a smaller amplitude than in the resonant case. This can be visualized by imagining a state on the Bloch sphere that starts at the south pole and rotates about different axes (see Fig. 2). For a rotation axis along the equator of the Bloch sphere, the path that the state would describe on the Bloch sphere would take it all the way to the north pole (i.e., large amplitude) as seen in Fig. 2(a). However, when the rotation axis cants closer to one of the poles, the path that the state describes on the Bloch sphere would not take it very far away from the south pole (i.e., small amplitude) as seen in Fig. 2(b). teractions in the experiments. Firstly, while a photodiode is sufficient to perform the measurements presented here, and was in fact used for Fig. 1(d), the signal to noise ratio can be improved using an avalanche photodiode (APD), e.g., Excelitas Technologies SPCM-AQRH.
As with the photodiode, the APD can simply be connected to the lock-in amplifier input (in voltage mode), where the low-pass filter will average over the short voltage pulses. When these experiments were implemented in the teaching labs at UNSW Sydney, the students performed all experiments using standard photodiodes. Secondly, as shown in Fig. 1(c), we have added a pair of homemade electromagnetic coils, powered by a DC lab power supply, to generate a well-defined and homogeneous B 0 field. Each coil is constructed from a reel and an L-shaped mount (see Fig. 8). A hole at the center of the reel, allows the addition of a core, in our case a supermendur cylinder, to further increase the magnetic field. The advantage over permanent magnets is that the magnetic field strength can be easily tuned by adjusting the current going through the coils, and the magnetic field direction is known as well. In App. A, we present more details on the design and geometry of the homemade electromagnets. A low-cost, commercial solution are the Adafruit Industries 3875 electromagnets (see e.g. Digi-Key 1528-2691-ND. An alternative way to measure a more homogeneous ensemble of NV − centers would be to interchange the multi-mode fiber for detection with a single-mode fiber to collect signal from a smaller sample volume. The sample used for the experiments presented here is a commercially available, off-the-shelf, 100 -oriented CVD diamond sample from Element Six , 28 and is used without any processing. The dimensions of the sample are 2.6 × 2.6 × 0.25 mm, and the nitrogen concentration is < 1 ppm. In Appendix C of Ref. 1, there is a comparison between this sample and a high-pressure, high-temperature (HPHT) diamond sample that was used for the experiments in Ref. 1. The CVD sample has a ∼ 1000 times lower PL signal, and while the spin signal decreases proportionally, we obtain high-quality data with only ∼ 10 times longer measurement times. In contrast, the CVD sample has a much longer coherence time than the HPHT sample, on account of the lower N concentration, that allows the observation of effects originating from coupling to nuclear spins. More recently, we have also started using a commercially available sample from Element Six with an NV center density of ∼ 300 ppb, optimized for quantum experiments. 29 This sample has a > 50× higher ODMR signal than the CVD sample, 28 and was used for the experiments in Fig. 4(a).

VI. EXPERIMENTS
All experiments described here are based on electron spin resonance (ESR), that allows us to coherently control the electron spin of the NV − center and perform different types of measurements depending on the exact pulse sequence. More precisely, the type of ESR experiment that we perform here is termed opticallydetected magnetic resonance (ODMR), since the electron spin state is read out (and initialized) via optical means, i.e., laser excitation and spin-dependent photon emission via photoluminescence (PL). The exact process was already described in Ref. 1 and many other publications. 4, 5,7,8 As introduced in Ref. 1 and can be seen in Fig. 1(a), there are a number of different ESR transitions that can be driven with an oscillating or rotating magnetic field B 1 . For B 0 = 0 T and neglecting the nuclear spin, all transitions are degenerate and centered around D = 2.87 GHz. When the magnetic field is turned on, the Zeeman splitting (Eq. 3) lifts this degeneracy and leads to different resonance frequencies for m s = 0 ↔ m s = −1 and m s = 0 ↔ m s = +1. When the interaction with the 14 N nuclear spin is included, there are six possible transitions -two each for m I = −1 (red arrows), m I = 0 (orange arrows), and m I = −1 (yellow arrows). These transitions will keep the nuclear spin state unchanged, but modify the electron spin state.

A. Magnetic field dependence
We start by measuring the ODMR spectrum of the diamond sample as a function of magnetic field B 0 .  Fig. 3). The B 1 field created by the printed-circuit-board antenna is applied along the [001] direction, and therefore perpendicular to the sample surface. The B 0 field created by the coils of the electromagnet is applied parallel to the sample surface. By rotating the CVD sample on the sample holder, different B 0 directions with respect to the crystallographic axes can be created. More specifically, when the sample is mounted in a square orientation [i.e., like 0, as shown in Fig. 3(a)] the B 0 field acts along the [100] direction. This results in all four NV − orientations experiencing the same effective angle between their NV − axes and B 0 and having the same two ODMR frequencies. However, when the sample is mounted in a diamond orientation [i.e., like &, as shown in Fig. 3 direction, however, we introduce a slight misalignment and let the students work out the exact orientation of the magnetic field. They can achieve this by recording B 0 -dependent ODMR spectra and mapping the observed splittings onto the calculated transition energies from Hamiltonian H S (Eq. 2) using Matlab. We plot such a magnetic field map in Fig. 1(d). Here, the magnitude of B 0 can be increased without altering the magnetic field direction by increasing the current through the electromagnet (see top x-axis) from 0 A to 1.0 A in 21 steps. For each B 0 magnitude an ODMR spectra was recorded to build up the 2D ODMR map. The students then fit the simple electron spin Hamiltonian H S to the data to work out the coil constant that describes the relationship between current and B 0 field. The dashed lines in Fig. 1(d) show the best fit for this data. The fit indicates that B 0 is applied along the [0.56 0.83 0.03] direction and the coil constant is 11.9 mT/A.
Note that the resistance of the electromagnetic coils is about 15 Ω and that it increases when the coils heat up at larger driving currents. It is therefore important to use a current source to control the magnetic field. Furthermore, the coils will also heat up the diamond sample and therefore the NV − center environment. At a high current, the ODMR spectrum will shift according to the temperature dependency of the zero-field-splitting dD/dT = −74.2 kHz/K. 30 For our experimental setup a current of 1 A heats up the diamond sample by ∼ 40 K and shifts the zero-field splitting by ∼ 3 MHz.
In the teaching labs at UNSW Sydney, the students recorded eleven ODMR spectra for electromagnet currents between 0 A and 0.5 A in 0.05 A steps. Each spectrum was measured from 2.65 GHz to 3.09 GHz in 221 steps with a wait time per point of 600 ms using the photodiode. This resulted in a total measurement time of ∼ 25 minutes. While the measurements were running, the students were given a prepared Matlab script to examine the relationship between the coil constant and the magnetic field orientation. The students were then able to iteratively extract the relevant B 0 parameters by fit- ting their values to the measurement data, comfortably within the allocated time frame of the lab. This exercise addresses the first learning outcome in Section II.
B. Spectral signature of hyperfine coupled 14 N Having mapped out the Zeeman splitting of the S = 1 NV − center electron spin in Fig. 1(d), we can now take a closer look at the ODMR spectrum to examine the spectral signature of the hyperfine coupled 14 N nuclear spins. Owing to the optical properties of the NV − center, the electron-nuclear spin interactions for the native 14 N nucleus can be seen directly in the ODMR spectrum. The 14 N hyperfine interaction tensor (see Eq. 4) splits each electron energy level [see Fig. 1(a)], causing the electron spin transitions to become nuclear spin dependent, with a 2.16 MHz splitting in between the different I z projections. The different colored arrows in Fig. 1(a)    ange) and m I = +1 (red). These splittings can be resolved with a high-quality, low concentration NV − center ensemble sample (≤ 1ppm 31 ), a high ODMR frequency resolution and a low microwave power. 32 In Fig. 3(c) Fig. 3(c) displays a triplet structure that corresponds to the three different nuclear spin orientations. In the following sections, we will examine closely what effect the presence of the hyperfine lines has in coherent control experiments such as Rabi oscillations, Ramsey oscillations and Hahn echoes.
From the Rabi formula (see Section IV), we know that a lower microwave power results in a narrower HWHM (Half Width at Half Maximum) of the peaks in the spectrum. We map out this dependency in Fig. 4(a), where we plot five ODMR spectra recorded for different MW powers (open circles), together with a fit to the sum of three Lorentzian peaks with equal width and area (solid lines). At the lowest MW power, the three hyperfine peaks are clearly separated, while the power broadening at the higher powers leads to an overlap of the individual peaks.
We therefore need to work with low MW powers if we want to resolve the hyperfine splitting without power broadening. This requires a balance between resolution and signal strength, and we noticed that it is advantageous to work at low magnetic fields where the quenching of the photoluminescence is insignificant and the electron spin transitions from two different NV − orientations with slightly different angles overlap. For all following experiments, however, we choose a MW power that leads to some overlap between peaks [similar to the +3 dBm curve in Fig. 4(a)] in order to observe effects of simultaneously driving neighboring transitions.
In principle, it is also possible to work with NV − centres of a single orientation. To achieve this, a slightly misaligned and higher magnetic field is needed to split the electron spin transitions associated with the four different NV − orientations. 33 Fig. 4(d) shows an example spectrum of the electron spin transitions between m s = 0 and m s = +1 at B 0 = 14.3 mT. The three different nuclear spin states can be clearly resolved and are labeled accordingly. We fit the spectrum with the sum of three Lorentzian peaks with equal amplitude and FWHM shown as black line, with the individual Lorentzian peaks show as red lines.

C. Rabi oscillations in the presence of 14 N
Applying a microwave pulse in resonance with one of the three hyperfine lines, caused by the 14 N nuclear spin, induces coherent electron spin oscillations between the m s = 0 and m s = ±1 spin eigenstates, known as Rabi oscillations. However, it is important to realize that the microwaves are simultaneously driving the detuned hyperfine transitions, due to the finite linewidth of the individual lines as well as the power broadening as described by the Rabi formula [e.g. see overlap of peaks in Fig. 4(a,c,d)]. Figure 4(b) shows a Rabi oscillation experiment, driven with the microwave frequency in resonance with the m I = 0 transition [middle peak of the left triplet in Fig. 4(a)], and the data is fitted to the sum of two exponentially decaying sinusoids A α sin(2πf α ) exp(−τ /T α ) + A β sin(2πf β ) exp(−τ /T β ) (black line). The frequencies for the two sinusoids correspond to the on-resonance drive of the m I = 0 transition and the detuned drive of the m I = ±1 transitions. The two oscillation frequencies become apparent when taking the fast Fourier transform (FFT) of the Rabi oscillation experiment data [see Fig. 4(b) -inset], showing peaks at Ω res R ≈ 0.7 MHz and Ω det R ≈ 2.3 MHz. According to the Rabi formula (Eq. 5), the Rabi oscillation frequency is given by: where Ω res R is the on-resonance Rabi oscillation frequency and ∆ is the detuning. Ω res R is given by: where B 1 is the oscillating magnetic field component perpendicular to the NV − axis delivered by the antenna. 34 Therefore, the Rabi frequency of Ω R = 0.667 ± 0.005 MHz in Fig. 4(b) corresponds to a field amplitude of B 1 = 0.034 mT, and we keep it fixed to this value for all coherent control experiments. The faster frequency of Ω det R = 2.35 ± 0.06 MHz in Fig. 4(b) corresponds to Rabi oscillations on the detuned hyperfine transition and agrees well with the one calculated using Eq. 6, where ∆ A = 2.16 MHz is the detuning between adjacent hyperfine peaks: Ω det R = (Ω res R ) 2 + ∆ 2 A = (0.667 MHz) 2 + (2.16 MHz) 2 = 2.26 MHz. The Rabi coherence time extracted from the on-resonance component is T R 2 = 3.96 µs, which is most likely limited by the inhomogeneity of the B 1 field over the sample area. 26 The effect of power broadening and cross-talk between the different hyperfine lines can be well visualized in a so-called Rabi chevron experiment, where the MW frequency and the pulse length are both swept -or in other words, a Rabi oscillation measurement is performed for various MW frequencies. Note that due the large number of data points, a Rabi chevron is usually outside of the scope of a teaching lab. Fig. 4(c), shows such a measurement for the same triplet of peaks as in Fig. 4(a,b). Each of the three hyperfine lines shows its own chevron fringes, overlapping with the neighbouring lines, and leading to the deviation from a single, decaying sinusoid in Fig. 4(b).
In principle, a simultaneous drive of all three hyperfine lines with similar strength can be achieved by applying a MW drive with stronger power to increase the Rabi oscillation frequency beyond the hyperfine coupling (2.16 MHz). 26 However, in our experiments the MW power was purposefully kept low, so that the effect of hyperfine transitions can be observed. This allows the students to develop an intuitive understanding of the consequences and origin of power broadening, fulfilling the second learning outcome in section II.
The Rabi experiment discussed in this section is the gateway between performing spin resonance spectroscopy and applying discrete gate operations in pulse sequences. Extracting the Rabi frequency Ω R = 0.667 MHz allows us to calibrate the lengths of our MW pulses to perform controlled rotations around the Bloch sphere. Here a π-pulse corresponds to a MW pulse of correct length τ π = 0.75 µs to rotate the electron spin 180 • around the Bloch sphere, e.g. from m s = 0 to m s = −1. A π/2-pulse of length τ π/2 = 0.38 µs rotates the electron spin by 90 • , e.g. from the m s = 0 state to the |+y = 1 √ 2 (|0 + i|−1 ) superposition state at the equator of the Bloch sphere. Once the pulse lengths have been calibrated, experiments using pulse sequences, such as the coherence times measurements in the next sections can be implemented. See also Ref. 1 for a more detailed discussion and visualizations on the spin rotations.  Fig. 6. Channel 0 (CH0) provides the reference for the lockin amplifier, channel 1 (CH1) modulates the laser output for the optical initialization and readout pulses, channel 2 (CH2) modulates the "I" output, and channel 3 (CH3) modulates the "Q" output of the MW IQ modulator. The second π/2-pulse is either applied along the X-axis to implement the standard Ramsey sequence X π/2 − τ − X π/2 or about the Y -axis to implement the modified Ramsey sequence X π/2 − τ − Y π/2 .

D. Ramsey fringes in the presence of 14 N
A Ramsey experiment measures the free-induction decay or T * 2 coherence time of a spin or a spin ensemble. The T * 2 time is affected not only by noise during the pulse sequence, but also by detunings between the transition and the MW frequency. This is because the sequence does not contain any refocussing or echo pulses, in contrast to the Hahn echo and Carr-Purcell-Meiboom-Gill (CPMG) sequences discussed in Section VI E. Figure 5 shows the pulse sequence that we use in our experiments and the outputs of the individual channels (CH0 -CH3) of the Pulse Blaster. The spins are initialized into the |0 state by the first laser pulse and read out by the second laser pulse (see CH1). The MW pulse sequence (See CH2 and CH3) is applied only during the first half-cycle for the lock-in amplifier (see CH0), so that the lock-in detection scheme is sensitive to changes in the photoluminescence caused by the MW pulse sequence (see also Ref. 1). The first microwave π/2-pulse, applied along the X-axis, rotates the NV − electron spins from |0 to the Y -axis of the Bloch sphere, corresponding to the |+y = 1 √ 2 (|0 + i|−1 ) superposition state in the rotating frame. We then let the spins dephase freely for a duration τ , before applying a second microwave π/2pulse along the X-axis to transfer the NV − spin polarization from |+y to |−1 for readout. A free induction decay curve can be recorded by varying the free precession time τ and the resulting data shows an exponential decay in |−1 population as a function of τ with decay time T * 2 caused by decoherence and inhomogeneities. 34 For NV − , the main inhomogeneities come from the variations in local strain and electric fields, as well as temporal and spatial fluctuations of the surrounding 13 C nuclear spin bath. 4 In addition, a constant frequency detuning between the spins and the MW drive will lead to a sinusoidal modulation of the data. This is because the spins precess at the speed of the detuning in the reference frame of the MW source, causing them to rotate from |+y towards |+x or |−x and onward, depending on the sign of the detuning. As the second π/2-pulse in the sequence maps |+y → |−1 , |+x → |+x , |−y → |0 , and |−x → |−x , the constant detuning manifests as a sinusoidal oscillation in the |−1 population at the end of the sequence. As a consequence, a Ramsey sequence is often used for the accurate estimation of shifts in the spin qubit resonance frequency, as is important in magnetometry experiments. 35 In Fig. 6(a) we show in blue the Ramsey data obtained for driving the NV − in resonance with the m I = 0 transition [same as the Rabi oscillations in Fig. 4(b)]. The data shows the exponential decay superimposed with a sinusoidal modulation. We fit the data to the function where n corresponds to the number of expected transitions. 36 Looking at the corresponding FFT of the data in Fig. 6(c), we can clearly see two peaks. The peak at low frequencies indicates a slight detuning of our MW source to the m I = 0 transition frequency. The smaller peak at ∼ 2 MHz corresponds to the signal from the detuned m I = ±1 hyperfine lines, which are detuned from the m I = 0 transition by ∆ A = 2.16 MHz.
In Fig. 6(b), we show very similar data with the Ramsey experiment performed in resonance with the m I = −1 transition. The data now shows two discrete modulation frequencies, one for the m I = 0 electron spin subsystem that is detuned by ∆ A and a second one for the m I = +1 subsystem that is detuned by 2∆ A . The two frequency components can be clearly seen as peaks in the Fouriertransform of the decay data in Fig. 6(d).
For both Ramsey experiments in Figs. 6(a) and (b) we also show data that is obtained by applying the second microwave π/2-pulse along the Y -axis (purple data). This modifies the mapping of the spin before readout to |y → |y , |x → |0 , |−y → |−y , and |−x → |−1 . As the projection pulse along the Y -axis (purple data) is shifted by π/2 in phase with respect to the projection pulse along the X-axis (blue date), the sinusoidal modulations are also shifted in phase by π/2 as seen in Figs. 6(a) and (b). The possibility to modify the projection axis by adjusting the phase of the MW source forms the basis for quantum state tomography 37 and allows the extension of the teaching labs to include further experiments.

E. Hahn echo in the presence of 13 C
The Hahn echo experiment is designed to cancel out constant detunings and frequency offsets by adding an extra π-pulse to the sequence. The π-pulse inverts the spins halfway through the precession and as such refo- The blue circles correspond to the standard Ramsey sequence X π/2 − τ − X π/2 , while the purple circles correspond to the modified Ramsey sequence X π/2 − τ − Y π/2 . The solid lines are fits to the data with T * 2 = 900 ± 30 ns and T * 2 = 1210 ± 60 ns for the two pulse sequences, respectively. (b) Same as (a) measured for mI = −1. The extracted coherence times are T * 2 = 620 ± 30 ns and T * 2 = 620 ± 30 ns for the two pulse sequences, respectively. cusses variations in their Larmor precession frequency. 38 We have already introduced this experiment in Ref. 1 for the purpose of dynamical decoupling, where it helps to improve the transverse relaxation time of the relevant spins, allowing the coherence time to be extended beyond the free induction decay time.
The MW pulse sequence we implement for the Hahn echo measurement is X π/2 − τ − X π − τ − X π/2 . In this sequence the electron spins are initialised in the |0 state and rotated into the superposition state |y = 1 √ 2 (|0 + i |−1 ) using a MW π 2 -pulse. The spins are then left to freely precess for a length of time τ , after which they are inverted on the Bloch sphere with a π-pulse. The spins are again left to precess for the same amount of time, causing them to refocus before being rotated back into the |0 state with another π 2 -pulse for optical readout. The spin signal is then measured as a function of the free precession time τ , and the T Hahn 2 coherence time is given as the time constant of the exponential decay of the signal.
The Hahn echo can also be used to reveal interactions between spins -in this case the hyperfine coupling to the 13 C nuclei in proximity to the NV spins that have a natural abundance of 1.1%. The randomly distributed 13 C nuclei of the diamond lattice produce a local magnetic field that is modulated by their Larmor precession frequencies, effectively resulting in an oscillating background magnetic field. Since the Hahn echo sequence can only remove slow fluctuations in detuning, the Larmor precession of the nuclei will lead to gradual decoherence of the electron spin. However, if the refocussing pulse occurs after exactly one period of nuclear spin precession, the background field modulation in both free precession intervals is exactly the same and the coherence reappears.
In Fig. 7(a), we present the Hahn echo signal measured at a magnetic field of B 0 = 1.3 mT. From the exponential decay, we extract T Hahn 2 = 11.7 ± 0.2 µs. Upon further investigation, a measurement at a higher magnetic field of B 0 = 5.84 mT reveals the presence of the abovementioned coherence revival within the Hahn echo signal, as shown in Fig. 7(b). The decay and revival in Fig. 7(b) is attributed to the weak hyperfine coupling to the 13 C spins near the NV center lattice sites. 39 The signal revival at 2τ ∼ 30 µs corresponds to a frequency of ∼ 66.7 kHz, which is consistent with the f C = γ C B 0 = 62.5 kHz Larmor precession frequency of the spin-1/2 nuclei with of an NV − ensemble can be as long as 600 µs at room temperature. However, the coherence time reduces when the magnetic field is not parallel to the NV − axis, as is the case in our samples and experiments.
As the coherence revival rate is linked to the Larmor precession frequency of the 13 C nuclei and therefore proportional to the applied magnetic field, this revival is not present in Fig. 7(a) where the lower magnetic field of B 0 = 1.3 mT would have led to a revival at 2τ = 144 µs, longer than the T Hahn 2 coherence time.

VII. CONCLUSION
Based on the cost effective, room temperature teaching setup described in Ref. 1 with some small improvements (electromagnet and lower NV − density diamond sample), we have presented here a set of more advanced experiments that aim to provide students with an in-depth understanding of physical concepts that are important for various quantum technologies. The experiments are designed to highlight the interaction of the NV − centers with surrounding nuclear spins via the hyperfine interactions. The experiments allow the students to gain a better understanding of experimental imperfections as well as experience the manifestation of spin-spin interactions with the intrinsic 14 N nuclear spin and the bath of surrounding 13 C nuclear spins. These interactions become visible in the presented Rabi oscillation, free induction decay, and Hahn echo experiments.
Despite the weaker signal intensity of the lower-density NV − diamond sample used here, the experiments work robustly on a standard desk in room light, and have partially been implemented in a fourth year undergraduate course at UNSW Sydney. Furthermore, the possibility for students to design their own pulse sequences and collect measurement data in real-time, makes this setup also an interesting option for more complex experiments in an undergraduate thesis.  ] using a stainless steel screw. In addition we have added 12 mm long supermendur cylinders with 4 mm diameter into the holes to slightly increase the magnetic field (optional). The 1 mm diameter hole (drilled at a 45 • angle) allows the inner end of the copper wire to pass through. The copper wire we used is American Wire Gauge (AWG) No.22 with 0.644 mm diameter and 52.96 mΩ resistance per meter. Several 1 mm thick and 5 mm long cooling fences allow for better heat dissipation. We chose brass as material due to its non-magnetic properties and ease of manufacturability. (b) Sketch of the L-shaped mount. The 5.3 mm diameter hole is used to mount the reel, while the 1 mm diameter trench allows the copper wire from the reel to pass through. The connected reel and mount are screwed onto the optical breadboard via two 6.4 mm diameter holes on the L-shaped mount. (c-e) Magnetic field simulation results. The results were obtained by calculating the magnetic field produced by the wire loops using the Biot-Savart law in matlab. The simulation is done under the condition that the distance between the two reels is 32 mm and a 1 A current is passing through the coils. The supermendur cylinders were not included in the simulations. Two black areas in (c) and (d) indicate the physical locations of the coils. (c) shows a 2D map of the total magnetic field along the radial and axial directions of the coils, with (d) showing the axial magnetic field component along the axial direction at r = 0 m and (e) along the radial direction at x = 0 m. The magnetic field at the center is 14.5 mT [see also dip in (d) and peak in (e)] and the simulated power consumption is 10.6 W. The real power consumption is slightly higher due to the heat generated by the coils causing an increase in the copper wire resistance. The negative values in (e) indicate that the axial component of the magnetic field is pointing along the negative x direction. The gradient of the magnetic field is smallest at the center between the two reels, providing the most homogeneous magnetic field in both strength and direction.