LAB - Quantum Sensing with Solid State Spins (Week 1)#
In this lab, you will explore using RF electronics to probe the microwave cavities and coupling the cavities to NV centers in diamond to enable quantum magnetometry.
Specific aims are:
1. Understand the optical setup
2. Characterize the microwave cavities and the Helmholtz coil
3. Simulate strong coupling between NV centers and the microwave resonator
4. Achieve strong coupling between NV centers and the microwave resonator
5. Improving the avoided-crossing contrast
Please note the questions throughout the lab instructions. You will be expected to address these in your lab writeup. You are also expected to keep lab notebooks (can be digital, paper, or both) for notes, data, and analysis.
Aim 0: Understand the optical setup#
Prelab#
In this lab, you will be using a home-built optical setup containing the following main components:
a 3W 532nm free space laser (this is a class-4 laser requiring laser safety!)
a magnetic shield
microwave cylindrical resonators
a loop coupler
a Helmholtz coil
diamond samples (~3mm x 3mm x 0.5mm) containing NV centers
cavity mounts (PTFE disks)
Not shown in the picture of the setup are the electronics and the RF components needed to probe the microwave cavity. You will be tasked to put together the RF circuitry for running the experiment.
As mentioned during the demo session, the experiment is based mainly on E. Eisenach et al, Cavity-enhanced microwave readout of a solid-state spin sensor, which you may find to be a valuable reference.
Other helpful references may be:
Aim 1: Characterize the microwave cavities and the Helmholtz coil#
Lab Exercise#
For the following steps, please take extra care when handling the microwave cavities and the loop coupler!#
i) Estimate the cavity’s resonance frequency#
There are two microwave resonators provided to you. Both of them should have identical geometries, except with slightly different surface qualities. Take one of them and measure its outer radius and length. Use the measured geometry parameters and estimate the cavity’s resonance frequency.
Can you measure the resonance frequency via the VNA, whose bandwidth covers 100 kHz to 6 GHz?
What would be the expected resonance frequency if you stack the two resonators on top of another, with their central holes as aligned as possible?
Can you measure the combined stack’s resonance frequency via the VNA?
ii) Build the RF circuitry for probing the cavity’s resonance frequency#
The VNA has a send port (Port 1) and a receive port (Port 2). Port 1 sends out a microwave tone by an internal signal generator, and port 2 measures the transmitted/reflected signal to compute the S-parameters. Use the provided circulator and note the numbers written on its, connect its ports to the VNA and the loop coupler appropriately such that the loop coupler can be used to probe the cavity’s signal. Check with the TA after you are done.
iii) Measure the microwave resonator’s frequency#
Now with the loop coupler properly connected to the VNA via the circulator, first put on gloves and grab the pair of tweezers. Carefully place the resonator stack inside the magnetic shield on the cavity mount, as shown in the figure below.
Mount the loop coupler such that it is directly above the resonator. You should use the right-angle clamp to adjust the loop coupler position, as shown in the figure below.
Open the libreVNA GUI and connect to the device. Make sure you are monitoring the “S21” window. Change the center frequency to what you estimated above for the stacked resonator. Change the span to 20 MHz. By adjusting the position of the loop coupler relative to the resonator, you should be able to observe the cavity resonance dip. Once you have done so, play around with the position of the loop coupler and observe how the cavity’s observed loaded \(Q\) and resonance frequency \(\omega_c\) change. Write down your observations and explanations to these changes as a function of position. In Aim 2, you will need to methodically measure the cavity resonance frequency and quality factor as a function of \(z\).
Try repeating the above procedure except now with a thick microscope glass slide in between the two resonators in the resonator stack. Is there a difference in the cavity response from before? Now switch to using a thin microscope glass slide. Any difference now?
iv) Characterize the Helmholtz coil#
The Helmholtz coil we will be using actually has two coils, each with two electrical leads. Connect the red and black alligator clips to the leads in both coils. Note that red typically indicates the positive end, and black indicates the negative end. There is a correct way of making the connection. Otherwise, you may not be able to generate much magnetic field. Check with the TA if you are unsure.
After you are done, grab the alligator clips’ BNC cable ends and connect them to a BNC-Tee connector. Then, attach the BNC-banana connector adapter to the BNC-Tee connector, and insert the adapter into Channel 1 of the R&S power supply (HMP4040).
Next, turn on the Gauss meter.
What is the intrinsic reading displayed on the Gauss meter? Why is it non-zero?
Set Channel 1’s voltage to 1 V, and click on “Output” on the top right of the front panel of the power supply. Put the Gauss meter’s probe in between the loops and note down the measured magnetic field.
Estimate the number of loops per coil. You may find this Helmholtz coil calculator useful.
Increase the voltage in increments of 1 V all the way to 20 V, and write down the measured magnetic field for each voltage setting. Please do not set the voltage beyond 20 V without consulting the TA
Generate below a plot of magnetic field in Gauss (y-axis) vs Channel 1’s voltage (x-axis).
Aim 2: Simulate strong coupling between NV centers and the microwave resonator#
Prelab#
Here are some questions to better prepare you for the lab exercise.
Consult this product link for one of the diamond types. Obtain the nitrogen concentration under “Materials Properties” and assume that all nitrogen defects in the crystal lattice are converted to NV centers. Estimate the density of NV centers in diamond in units of NV/m\(^3\).
Hint: you may need to look up the atomic density of diamond and convert to its units to number of atoms/m\(^3\).
Given our diamond dimensions, 3mm x 3mm x 0.5mm, what is the effective coupling strength \(g\)? Assume a single spin coupling strength is \(g_0=0.01\) Hz. You may find the paper by Eisenach et al. helpful.
Now you should be equiped to do the lab exercise!
Lab Exercise#
In this lab session, you will be strictly working on creating a Python script that simulates the strong coupling dynamics between NV centers and the microwave resonator. In scientific research, it is often extremely valuable to be able to predict how the experimental results may vary depending numerous parameters. By producing a digital twin, you will be able to better understand the feasibility of using such a solid state spin-based system for performing quantum sensing.
i) Estimate the microwave cavity’s intrinsic quality factor#
Based on the data you have acquired last week, plot out the loaded \(Q\) as a function of distance (in \(z\)) between the loop coupler and the top of the cavity. Rough estimation of the distance is fine. If you did not record the data from last week, you may repeat the same measurements. Be sure to have the loop coupler as directly above the microwave resonator to remove any \(x,y\) dependence.
Remember that the loaded \(Q\propto 1/\kappa_c\) and the coupling rate to the loop coupler \(\kappa_1\) are both functions of \(z\). Fit your plot with an exponential function, e.g. \(\propto 1-exp(-z)\), and estimate the intrinsic quality factor, \(Q_0\propto 1/\kappa_0\).
Then, generate a plot of \(\kappa_1\) as a function of \(z\).
ii) Plot the frequency difference between \(m_s=0\) and \(m_s=+1\) states as a function of magnetic field (in Gauss)#
iii) Generate a 1D plot of the cavity reflection coefficient#
Recall that the cavity reflection coefficient is defined as:
For simplicity, let’s first take \(g=0\). Generate a 1D plot of the real part of \(\Gamma(\omega)\) as a function of \(\omega\). Choose and define \(\omega_c,\kappa_c,\kappa_{c1}\) based on what you measured in Aim 1 and calculated above.
iv) Generate 2D plots of the cavity reflection coefficient#
Now, use the effective coupling strength \(g\) you computed in the Prelab, generate 2D plots of the following, with the \(x\)-axis being magnetic field strength and \(y\)-axis being frequency:
In-phase component of \(\Gamma(\omega)\)
Quadrature component of \(\Gamma(\omega)\)
Amplitude of \(\Gamma(\omega)\)
Phase of \(\Gamma(\omega)\)
Select \(\omega_c,\kappa_c,\kappa_{c1}\) appropriately to ensure you can observe avoided crossing similar to what was observed by Eisenach et al. (Figure 2). You may assume \(\alpha=0\). \(\kappa_s=2/T_2\) where \(T_2\) is the coherence time of the NV centers. You may assume \(T_2=20~\mu\)s.
You may find using the Python module ipywidgets
helpful. By creating a widget with scroll bars, you may be able to freely change the parameters without manually typing in different values. Specifically, the interact_manual
method (see documentation here).
If you are done with the simulator and playing with various parameters, please ask the TA to give you a demo on how to handle diamond pieces. For the following sessions, you will be picking up and placing small diamond samples (3mm x 3mm x 0.5mm). It may benefit you to get some practice on handling small diamonds without dropping them.
Practice picking up a diamond without losing grip. Take care to not exert too much strength on tweezing, as too much force would sometimes to lead to “propelling” the diamond.
Practice placing down the diamond and nudging it gently from the sides to adjust its lateral position.
Practice having the diamond stand on its 0.5mm-thick side.