Single-Spin Magnetomechanics with Levitated Micromagnets

J. Gieseler,1 A. Kabcenell,1 E. Rosenfeld,1 J. D. Schaefer,1 A. Safira,1 M. J. A. Schuetz,1 C. Gonzalez-Ballestero,2, 3 C. C. Rusconi,4 O. Romero-Isart,2, 3 and M. D. Lukin1 1Physics Department, Harvard University, Cambridge, MA 02318, USA. 2Institute for Quantum Optics and Quantum Information of the Austrian Academy of sciences, A-6020 Innsbruck, Austria. 3Institute for Theoretical Physics, University of Innsbruck, A-6020 Innsbruck, Austria. 4Max-Planck-Institut für Quantenoptik, Hans-Kopfermann-Strasse 1, 85748 Garching, Germany. (Dated: December 24, 2019) We demonstrate a new mechanical transduction platform for individual spin qubits. In our ap- proach, single micro-magnets are trapped using a type-II superconductor in proximity of spin qubits, enabling direct magnetic coupling between the two systems. Controlling the distance between the magnet and the superconductor during cooldown, we demonstrate three dimensional trapping with quality factors around one million and kHz trapping frequencies. We further exploit the large mag- netic moment to mass ratio of this mechanical oscillator to couple its motion to the spin degree of freedom of an individual nitrogen vacancy center in diamond. Our approach provides a new path towards interfacing individual spin qubits with mechanical motion for testing quantum mechanics with mesoscopic objects, realization of quantum networks, and ultra-sensitive metrology.

Realizing coherent coupling between individual spin in the solid state [3], quantum sensing [20–22] and quan- degrees of freedom and massive mechanical modes is an tum communication [23]. outstanding challenge in quantum science and engineer- Before proceeding, we note that low dissipation me- ing. Such a coupling could be used to create mechani- chanical resonators based on magnetic levitation have cal systems with a strong quantum non-linearity, allow- been explored previously [24, 25]. However, experiments ing preparation of macroscopic quantum states of motion with superconducting levitation have so far been limited [1]. In addition, mechanical systems can be used to medi- to millimeter-scale magnets [24, 26–28]. A recent ex- ate effective spin-spin interactions between distant spin- periment [29] demonstrated superconducting levitation qubits [2], enabling applications ranging from quantum of micro-magnets, but without spin-mechanical coupling information processing [3] and sensing [4–6] to tests of and with much lower frequencies and Q-factors than fundamental physics [7,8]. One particularly promising shown here. Levitated spin-mechanical systems in which approach is to engineer a strong spin-mechanical coupling the spin is hosted inside the resonator have been imple- via magnetic field gradients [9–13]. Achieving strong mented with nano-diamonds containing NV-center de- spin-resonator coupling requires a combination of high fects trapped in optical tweezers [30, 31], Paul traps [32– quality mechanical resonators, strong magnetic field gra- 34] and magneto-gravitational traps [35, 36]. Nonethe- dients, and spin qubits with very long spin coherence less, the challenge remains to integrate these systems times. with strong magnetic field gradients, long coherence NV- In this Letter, we propose and demonstrate a new centers and operation under ultra-high vacuum condi- platform for strong spin-mechanical coupling based on tions. Our approach fulfills all these criteria simultane- levitated microscopic magnets coupled to the electronic ously (Fig.1). ground state manifold of a single nitrogen vacancy (NV) Levitation of micromagnets We levitate single hard center in diamond (Fig.1). The key idea is to utilize a micro-magnets with a thin film of the type-II supercon- levitated magnet that is localized in free space by elec- ductor (sc), yttrium barium copper oxide (YBCO) (c.f. tromagnetic fields. In such a system, dissipation is min- Fig.1). Since we do not apply additional magnetic fields, imized since there is no direct contact with the environ- the micromagnet is the only magnetic field source. Thus,

arXiv:1912.10397v1 [quant-ph] 22 Dec 2019 ment. Specifically, we make use of a levitating mechan- the magnetic flux through the YBCO film is determined ical resonator based on magnetostatic fields. This ap- only by the distance hcool between the magnet and the proach not only avoids clamping losses, but also circum- YBCO film, the orientation of the magnet θcool, the mag- vents photon recoil and heating associated with optical netization of the magnet, and its radius a. After cooling levitation [14–16] and is therefore predicted to yield large the YBCO film below its critical temperature Tc ≈ 90K, mechanical quality factors [17, 18]. it becomes superconducting and magnetic flux that pen- In addition, the levitated magnet naturally generates etrates the film is frozen in. As a consequence, below the strong magnetic field gradient that is required for Tc, motion of the magnet induces currents in the super- spin-mechanical coupling. We specifically demonstrate conductor that counteract changes in the magnetic field. the coupling to an individual NV-center, one of the most This allows for stable 3D trapping using a procedure il- studied color centers in diamond [19]. Besides optical lustrated in (Fig.2a). The levitation height (hlev) and initialization and readout, the NV-center features long trapping frequencies ωj (j = x, y, z) depend on the con- coherence times even at room temperature, which makes ditions during cooldown, which we can control by adjust- it an attractive candidate for scalable quantum networks ing the relative distance between the superconductor and 2 a) b) a) b) camera nitrogen-vacany center

hlev hcool diamond silicon pocket cryostat magnetic field set height cool down approach & levitate membrane c) 100 100 superconductor ) piezo stack 1 2 101 1

10 m 10

( 0

2 10

10 2 x y y 10 2 1

FIG. 1. Concept of the experiment. (a) Individual micro- ) 10 z 0 200 400 H 3 3 scopic hard ferromagnets are isolated in microfabricated pock- / 10 10

2 time (s)

ets and levitated with a type II superconductor, exploiting its m ( flux trapping properties. The magnet’s stray field allows for 4 4

D 10 10 S

efficient coupling between the magnet’s motional degrees of P freedom and bulk nitrogen vacancy centers in the nearby dia- 10 5 z 10 5 mond, which generally feature long spin coherence times. (b) The coupling λ depends on the relative orientation of the NV- 10 6 10 6 center, ns, the magnetic moment nµ, the direction of motion of the magnet nm, the distance between the magnet and the 10 7 10 7 0 0.805 0.815 0.83 0.84 0.85 2.25 2.30 NV-center r , the magnet radius a and the frequencies of the frequency (kHz) motion. d) e) 1.2 z y x 1.0 103 ) the particle during cooldown (hcool). 6 0 1

( 0.8 We observe the levitated magnet through a long work- r o t ing distance microscope objective that is positioned out- c a

f 0.6 side the vacuum chamber, and record its motion with - Q a fast camera (Fig.2b). From the video frames, we ex- frequency (Hz) 102 0.4 tract time-traces of the particle’s center-of-mass position 23.2 m 15.5 m in the camera coordinate system (xc(t), yc(t)) and cal- 0.2 culate their power spectral densities (PSD). The distinct 2 3 4 5 6 7 1 2 peaks in the PSD shown correspond to the three center- hlev/a frequency (kHz) of-mass modes ωj/2π (Fig.2c). Fig.2d displays the center-of-mass frequencies as a FIG. 2. Mechanical properties of levitated micro-magnets. ¯ (a) To levitate the magnets, we adjust the magnet-sc dis- function of the normalized levitation height hlev = hlev/a tance above the critical temperature Tc. Then we cool the for two particles with radius a1 = 23.2 ± 0.7µm and sc below Tc to freeze in the magnetic flux from the mag- a2 = 15.5 ± 0.3µm, respectively. The lines are a fit net. After cooldown, we reduce the distance until the mag- ¯ ¯−β to a power law f(hlev) = fmaxhlev , where fmax is the net begins to levitate. (b) We observe the magnet motion frequency in the limit when the gap between the parti- through a microscope objective and record its motion with cle and the superconductor goes to zero. In our experi- a video camera. (c) Power spectral density extracted from ment, we find fmax = (2.3 ± 0.4, 2.4 ± 0.4, 5.6 ± 1.0)kHz video analysis of the magnet motion in the vertical (hori- and β = (1.9 ± 0.1, 2.1 ± 0.1, 2.0 ± 0.1) for particle 1 zontal) direction, corresponding to the data shown as circles in (e). The inset shows a typical ringdown measurement of and fmax = (8.8 ± 1.1, 9.5 ± 1.1, 25.2 ± 3.3)kHz and β = (2.1±0.1, 2.1±0.1, 2.3±0.1) for particle 2, which is in the y-mode (ωy = 0.839kHz) from which we extract the Q- good agreement with the expected value of β = 2.5 from factor. (d) Frequencies of center-of-mass motion as a function of the levitation height normalized to magnet radius. Dashed a simple dipole model [SI]. The measured center-of-mass (solid) lines show the three center-of-mass modes for mag- frequencies are comparable to those achieved with optical net 1 (2). (e) Q-factor as a function of trap frequency for levitation [37] and significantly exceed motional frequen- magnet 2. Each symbol corresponds to a different levitation cies in Paul traps [38] and magneto-gravitational traps experiment and the different colors corresponds to the three [35, 36]. However, the observed dependence of the maxi- different center-of-mass modes. mum frequency on the particle radius is stronger than the dp dipole model’s prediction of fmax ∝ 1/a. We attribute this to the multi-domain nature of our particles and note Fig.2e shows the Q-factors of magnet 1 for three dif- that spherical particles as large as a ∼ 1µm, which we ferent levitation heights for all three translational modes. expect to be achievable with this technique, can be single We measure the dissipation with ring down measure- domain [39]. ments, exciting one mode with an AC magnetic field 3 and observing its energy decay (Fig.2c inset). From an measure the magnet’s motion with the NV, we excite exponential fit of the energy decay we extract the de- one of its modes with a broadband fluctuating magnetic cay time 1/γj and the Q-factor Qj = ωj/γj for each field, which drives it into a quasi-thermal state and al- mode. The measured Q-factors are around one million lows us to observe it as a peak in the power spectral and depend only weakly on the trapping frequencies and density of the NV PL counts (Fig.3f). We confirm that thus on the levitation height. Air damping can be ruled the peak is due to the moving magnet with the camera out at our experimental conditions with pressures below (Fig.3d). The camera and NV measurements are taken 10−5 mBar and the most likely source of dissipation is sequentially. Notably, we observe a small systematic fre- the magnet-superconductor interaction. Note that, even quency shift of ≈ 1 Hz between the NV and camera mea- though this Q-factor is somewhat lower than what has surements, which is likely due to the laser field turned on been demonstrated with non-magnetic optically levitated during the NV measurement. 2 [16] and nano-fabricated mechanical resonators [40–43], it The mean spectral power in the NV peak is hcNVi = 2 2 represents the state of the art for magnetized resonators s hδωNVi, where s is the slope of the ODMR signal at [44, 45] and the ultimate limit, in particular for magnets the microwave frequency, which we measure by apply- with a < µm, is still an open question. ing a calibration tone to the microwave signal. The 2 Coupling to NV-center Next, we demonstrate cou- mean spectral power in the camera measurement, hx i, pling the motion of a levitated micromagnet to the allows us then to extract the spin-mechanical coupling as p 2 2 electronic spin associated with an individual negatively λg = xzp hδωNVi /hx i. To measure the coupling and charged NV-center. In our sample, NV-centers are hosted confirm the thermal character of the driven mode, we inside a diamond slab and implanted dimpl ∼ 40 nm be- consider the area under the PSD integrated over a time low the diamond surface. The diamond replaces the glass interval much shorter than 1/γ, and construct its distri- slide of the previous experiment and is placed across bution over repeated measurements. For both the cam- the pocket that contains the magnet. The pocket is era and the NV measurements, the distribution agrees ∼ 80µm deep and the magnet radius is a3 = 15.1±0.1µm. with an exponential distribution P (E) = β exp(−βE), We levitate the magnet with the superconductor zmd = where the decay constant β is the inverse of the vari- 2 44 ± 5µm below the diamond using the same method as ances hx i and hδωNVi (Fig.3e,g). The resulting cou- before. The NV-center is located at (xd, yd) = (83, 29) ± pling strength is 48 ± 2mHz, in satisfactory agreement 5µm with respect to the magnet center such that the with the theoretical value for the gradient coupling to a 3 µ0ρµa distance between the center of the magnet and the NV- dipole λg = γe r04 xzpfg(θ) = 2π × (18 ± 3)mHzfg(θ) 0 p 2 2 2 center is |r | = (zmd + a3) + xd + yd = 99 ± 5µm. (Fig.4b). Here fg(θ) is a function on the order of 1 that The NV-center’s electronic ground state has spin S = 1 depends on the relative position and orientation of the with the lower-energy |ms = 0i level separated from the NV-center and the magnet. |ms = ±1i levels by a zero-field splitting Dzf/(2π) ≈ Discussion We now discuss the prospects of using this 2.87 GHz and its symmetry axes ns aligned along one of system to achieve strong coupling. The minimal sepa- min 0 four crystallographic orientations set by the tetrahedral ration dq = |r | − a between magnet and NV-center symmetry of the diamond lattice. A microwave (MW) is given by the NV implantation depth and the onset of signal at the transition frequency ωMW, drives the tran- strong attractive surface forces that will make the magnet sition |ms = 0i → |ms = ±1i which results in a decrease stick to the diamond surface. For a given separation and −n of the PL signal. The spin-dependent PL of the NV de- assuming that the frequency scales as ωj /2π = αa , min fect is due to a non-radiative intersystem crossing decay the radius a = (n + 3)/(5 − n)dq yields the maximum pathway, which also allows for efficient spin-polarization gradient coupling for a dipolar particle. A conservative min in the |ms = 0i spin sublevel through optical pumping gap dq = 0.25µm, α = 15kHz µm, corresponding to our [46]. The magnetic field dependent PL can therefore be observations in Fig.2, and n = 1 for the dipole model, re- used to optically detect magnetic fields [47], which we sults in a = 0.25µm and λg /2π ∼ 2.6kHz. Since the will use to sense the motion of the magnet (Fig.3b). motional frequency can be reduced by adjusting the levi- Fig.3c shows the optically detected magnetic reso- tation height, one can even reach the elusive ultra-strong nance (ODMR) spectrum of the NV-center with a fit to coupling regime λg > ωj [48]. 2 ˜ a Lorentzian corresponding to the |+1i transition. The The cooperativity C = λ QT2~ /(2πkBT ) > 1 marks spectrum is measured with a home-built fluorescence mi- the onset of coherent quantum effects in a coupled croscope that we integrated with the cryostat, replacing spin-phonon system. With a mechanical Q-factor of the long working distance objective outside the chamber 108, which has been demonstrated in levitated systems with a high NA objective inside the chamber to maxi- [49], the coupling exceeds the thermal decoherence rate mize the collection efficiency of the fluorescence photons Γth/2π = kBT/(2π~Q) = 0.8kHz at T = 4K. NV-centers (Fig.3a). In the presence of a microwave tone, the spin in bulk diamond, such as the sample used in our exper- mechanical coupling λg causes a variation in PL, since a iment, can have up to second long extended coherence displacement x of the magnet shifts the electron spin res- times T˜ at these temperatures using spin manipulation p 2 onance by δωNV = (λg/xxp)x, where xzp = ~ /2mωx such as multi-pulse dynamical decoupling sequences lim- is the zero point motion, in this case of 24 ± 1fm. To ited only by pulse errors [50, 51]. The minimum spin 4

a) b) d) e)

APD 532nm dichroic mirror

flip-mirror camera microwave

MW generator cryostat

c) f) g)

FIG. 3. Measurement of spin-motional coupling. (a) Schematic of the experiment. (b) Nitrogen-vacancy ground state spin levels. Resonant microwaves drive transitions between the ms = 0 and ms = ±1 states. The transition frequency depends on the magnetic field, which depends on the magnet-NV distance. The ms = 0 is brighter than the ms = ±1 states. Hence, for near-resonant MW driving, the fluorescence intensity depends on the magnet position. (c) ODMR spectrum. The magnet motion is measured by recording a timetrace of the fluorescence photons while applying a MW tone at the steepest slope of the transition (2.918 GHz) (solid red line). We measure the slope by frequency modulating the MW tone (red dashed lines). (d) Power spectral density of the video timetrace and (f) fluorescence timetrace. The gray area corresponds to the numerical 2 2 2 integration of the variance hx i, hcNVi and hccali. The second sharp peak is the calibration peak due to the MW frequency modulation. The black dashed line represents the photon shot noise. Histograms of the variances reveal the thermal character of the motional state during both the (e) camera and (g) NV measurement which we fit (red line) to extract the variances hx2i 2 and hδωNVi, respectively.

manipulation (and therefore sensing) frequency in such (spin) density. Since the Larmor frequency ωL = γ0B0, sequences is set by the power spectrum of the noise. It is ωl is higher than the translational mode frequencies even typically a few kHz for bulk diamond NV-centers, which without additional fields and can be tuned with moder- is within reach for the mechanical frequencies in our cur- ate magnetic fields ∼ 10mT to MHz frequencies. The rent geometry. Hence, this system can reach the high co- high frequencies of the librational modes make them in- operativity (C > 1) and even the strong coupling regime accessible to our current detection based on video anal- (λ > 2π/T˜2, Γth). Such a strong coupling enables ground- ysis and DC magnetometry. Future work will explore state cooling, quantum-by-quantum generation of arbi- these modes using optical interferometry and SQUID trary states of motion [9], and spin-spin entanglement [53] or NV-AC magnetometry [4, 38]. The rotational [2]. modes couple to the NV-center with a dipole-dipole 3 coupling λ = γ µ m a f (θ, ϕ), where f (θ, ϕ) is Besides the translational degrees of freedom, levitated dp e 0 zp r03 dp dp a function ∼ 1 that depends on the relative position particles are free to rotate. For the hard magnets used and orientation of the NV-center and the magnet and in our experiment, the anisotropy energy strongly cou- m = p γ ρ /2V is the zero point magnetization ples the particle orientation to the magnetization axis. zp ~ e µ mag of the Kittel magnon [54]. The weaker distance depen- The coupling leads to hybrid magneto-rotational modes, dence yields λ = 0.4kHz for a a = 5µm magnet at which correspond to a librational mode at frequency ω = dp √ l 5µm distance from the NV-center (Fig.4b). This param- ω ω that precesses around the local magnetic field B L I 0 eter regime is readily accessible with the experimental at the Einstein-deHaas frequency ω = ρ V /(I γ ) I µ mag 0 0 approach presented here, and it is sufficient to proba- due to the intrinsic spin angular momentum of the po- bilistically cool the librational mode near its ground state larized electrons in the magnet. Here, B is the sum 0 [55] (based on the requirements for a harmonic oscillator of the field due to the sc and additional external fields 6 2 mode ω0/2π ∼ MHz, λ/2π ∼ 100Hz, Q > 10 ). and I0 = 2ρmVmaga /5 is the moment of inertia, Vmag being the volume of the magent and ρm (ρµ) its mass 5

even for magnets with a ∼ 5µm. Magnets of this a) b) 105 106 size can be levitated with the experimental technique f0 introduced here. The combination of high mechanical f1 3 f2 10 Q-factor, strong spin-mechanical coupling, and long 5 10 f3 spin-coherence is key for a range of applications such f4 101 as magnetometers, accelerometers and gyroscopes [53], where the magnet is the sensor which is read out 104

coupling (Hz) through the NV-center [4]. Furthermore, it may enable frequency (Hz) 10 1 the exploration of new phenomena, including dynamics

dp/2 3 between a levitated nanomagnet and a single flux 10 gr/2 10 3 vortex[57–59], precession of a non-rotating magnet due 10 1 100 101 10 1 100 101 to its intrinsic spin angular momentum [60], preparation radius a ( m) radius a ( m) of non-Gaussian quantum states [9], mechanically mediated quantum networks[3, 61], detection of dark FIG. 4. Prospect for frequencies and couplings. (a) Mode matter[56], and measuring the magnets internal degrees frequency following dipole model [52] for hlev/a = 3. The of freedom [54]. lowest (solid) lying modes correspond to the center-of-mass and the highest (dashed) lying modes to rotational motion. Data points are the experimental frequencies for particle 1 This work was supported by the NSF, the Center Fig.2d. MHz frequencies are predicted for sub-µm particles. for Ultracold Atoms (CUA), the ONR MURI Quantum (b) Gradient (solid) and dipole-dipole (dashed) couplings, re- Opto-Mechanics with Atoms and Nanostructured Dia- spectively. Straight lines show couplings for constant gap- particle size ratio d¯ = d/a = 5.5 (color) and d¯ = 1 (black). mond (QOMAND), the Vannever Bush Faculty Fellow- Curved lines (gray) show coupling for constant gap d = 250nm ship, and the Moore Foundation. This material is based with a maximum coupling at a = d. The black data point is upon work supported by the National Science Foundation the experimental gradient coupling for particle 3. Graduate Research Fellowship Program under Grant No. DGE1144152 and DGE1745303. Any opinions, findings, and conclusions or recommendations expressed in this Outlook These considerations indicate that our material are those of the author(s) and do not necessar- approach is a promising platform for quantum nanome- ily reflect the views of the National Science Foundation. chanics. Our experimental technique also allows us JG was supported by the European Union (SEQOO, to achieve levitation with the superconductor in the H2020-MSCA-IF-2014, no. 655369). AK and AS were Meissner state, and thus presents a path forward to supported by the Department of Defense (DoD) through observe precession due to the intrinsic spin angular the National Defense Science and Engineering Graduate momentum of the magnet with applications in highly Fellowship (NDSEG) Program. C. G. B. was supported sensitive magnetometry [56]. Since our mechanical by the European Union (PWAQUTEC, H2020-MSCA- resonator is all magnetic, we maximize the spin to mass IF-2017 no. 796725). This work was performed in part √ ratio ρµ/ ρm for a given magnetic material, which at the Center for Nanoscale Systems (CNS), a member of maximizes the spin-mechanical coupling. This leads the National Nanotechnology Coordinated Infrastructure to strong spin-mechanical coupling even for moderate Network (NNCI), which is supported by the National Sci- experimental parameters. In addition, this system ence Foundation under NSF award no. 1541959. CNS is features libration modes, which are expected to reach part of Harvard University. We gratefully acknowledge an unprecedented spin-mechanical parameter regime Frank Zhao for assistance with sample magnetization.

[1] P Rabl. Cooling of mechanical motion with a two-level with a Single-Spin Qubit. Science, 335(6076):1603–1606, system: The high-temperature regime. Physical Review March 2012. B, 82(16):165320, October 2010. [5] S D Bennett, N Y Yao, J Otterbach, P Zoller, P Rabl, [2] M J A Schuetz, G Giedke, L M K Vandersypen, and J I and Mikhail D Lukin. Phonon-Induced Spin-Spin In- Cirac. High-fidelity hot gates for generic spin-resonator teractions in Diamond Nanostructures: Application to systems. Physical Review A, 95(5):052335–26, May 2017. Spin Squeezing. Physical Review Letters, 110(15):156402, [3] P Rabl, S J Kolkowitz, F H L Koppens, J G E Harris, April 2013. P Zoller, and Mikhail D Lukin. A quantum spin trans- [6] Michael S J Barson, Phani Peddibhotla, Preeti ducer based on nanoelectromechanical resonator arrays. Ovartchaiyapong, Kumaravelu Ganesan, Richard L Tay- Nature Physics, 6(8):602–608, May 2010. lor, Matthew Gebert, Zoe Mielens, Berndt Koslowski, [4] S Kolkowitz, A C Bleszynski Jayich, Q P Unterreith- David A Simpson, Liam P McGuinness, Jeffrey McCal- meier, S D Bennett, P Rabl, J G E Harris, and Mikhail D lum, Steven Prawer, Shinobu Onoda, Takeshi Ohshima, Lukin. Coherent Sensing of a Mechanical Resonator Ania C Bleszynski Jayich, Fedor Jelezko, Neil B Manson, 6

and Marcus W Doherty. Nanomechanical Sensing Using P Cappellaro, L Jiang, M V Gurudev Dutt, E Togan, Spins in Diamond. Nano Letters, 17(3):1496–1503, March A S Zibrov, A Yacoby, R L Walsworth, and M D Lukin. 2017. Nanoscale magnetic sensing with an individual electronic [7] J van Wezel and T H Oosterkamp. A nanoscale exper- spin in diamond. Nature, 455(7213):644–647, October iment measuring gravity’s role in breaking the unitarity 2008. of quantum dynamics. Proceedings of the Royal Soci- [22] S Kolkowitz, A Safira, A A High, R C Devlin, S Choi, ety A: Mathematical, Physical and Engineering Sciences, Q P Unterreithmeier, D Patterson, A S Zibrov, V E 468(2137):35–56, November 2011. Manucharyan, H Park, and M D Lukin. Probing John- [8] William Marshall, Christoph Simon, Roger Penrose, and son noise and ballistic transport in normal metals with a Dik Bouwmeester. Towards quantum superpositions of single-spin qubit. Science, 347:1129–1132, March 2015. a mirror. Physical Review Letters, 91:130401, September [23] B Hensen, H Bernien, A E Dréau, A Reiserer, N Kalb, 2003. M S Blok, J Ruitenberg, R F L Vermeulen, R N Schouten, [9] P Rabl, P Cappellaro, M Dutt, L Jiang, J Maze, and C Abellán, W Amaya, V Pruneri, M W Mitchell, Mikhail D Lukin. Strong magnetic coupling between an M Markham, D J Twitchen, D Elkouss, S Wehner, T H electronic spin qubit and a mechanical resonator. Physi- Taminiau, and R Hanson. Loophole-free Bell inequality cal Review B, 79(4):041302, January 2009. violation using electron spins separated by 1.3 kilometres. [10] O Arcizet, V Jacques, A Siria, P Poncharal, P Vin- Nature, 526(7575):682–686, October 2015. cent, and S Seidelin. A single nitrogen-vacancy defect [24] H Barowski, K M Sattler, and W Schoepe. Static and dy- coupled to a nanomechanical oscillator. Nature Physics, namic forces on a permanent magnet levitating between 7(11):879–883, September 2011. superconducting surfaces. Journal of Low Temperature [11] B Pigeau, S Rohr, L Mercier de L eacute pinay, Physics, 93(1):85–100, 1993. A Gloppe, V Jacques, and O Arcizet. Observation [25] K Gloos, J H KOIVUNIEMI, W SCHOEPEN, J T of a phononic Mollow triplet in a multimode hybrid Simola, and J T TUORINIEMI. VISCOMETER UTI- spin-nanomechanical system. Nature Communications, LIZING A FLOATING CHARGED MAGNETIC PAR- 6:8603, October 2015. TICLE. Physica B: Condensed . . . , pages 119–120, July [12] Florian Mintert and Christof Wunderlich. Ion-Trap 1990. Quantum Logic Using Long-Wavelength Radiation. [26] J Druge, C Jean, O Laurent, M-A Méasson, and I Favero. Physical Review Letters, 87(25):257904–5, November Damping and non-linearity of a levitating magnet in ro- 2001. tation above a superconductor. New Journal of Physics, [13] Sungkun Hong, Michael S Grinolds, Patrick Maletinsky, page 075011, July 2014. Ronald L Walsworth, Mikhail D Lukin, and Amir Ya- [27] V V Nemoshkalenko, E H Brandt, A A Kordyuk, and coby. Coherent, Mechanical Control of a Single Electronic B G Nikitin. Dynamics of a permanent magnet levitating Spin. Nano Letters, 12(8):3920–3924, August 2012. above a high-Tc superconductor. Physica C: Supercon- [14] Oriol Romero-Isart, Mathieu L Juan, Romain Quidant, ductivity, 170(5-6):481–485, October 1990. and J Ignacio Cirac. Toward quantum superposition of [28] John R Hull and Ahmet Cansiz. Vertical and lat- living organisms. New Journal of Physics, 12(3):033015, eral forces between a permanent magnet and a high- March 2010. temperature superconductor. Journal of Applied Physics, [15] Darrick E Chang, C A Regal, S B Papp, D J Wil- 86(11):6396–6404, December 1999. son, J Ye, O Painter, H Jeff Kimble, and P Zoller. [29] Tao Wang, Sean Lourette, Sean R O’Kelley, Metin Kayci, Cavity opto-mechanics using an optically levitated Y B Band, Derek F Jackson Kimball, Alexander O nanosphere. Proceedings of the National Academy of Sci- Sushkov, and Dmitry Budker. Dynamics of a Ferromag- ences, 107(3):1005–1010, January 2010. netic Particle Levitated over a Superconductor. Physical [16] Vijay Jain, Jan Gieseler, Clemens Moritz, Christoph Del- Review Applied, 11(4):044041, April 2019. lago, Romain Quidant, and Lukas Novotny. Direct Mea- [30] Levi P Neukirch, Jan Gieseler, Romain Quidant, Lukas surement of Photon Recoil from a Levitated Nanoparti- Novotny, and A Nick Vamivakas. Observation of nitrogen cle. Physical Review Letters, 116(24):243601, June 2016. vacancy photoluminescence from an optically levitated [17] O Romero-Isart, L Clemente, C Navau, A Sanchez, and nanodiamond. Optics letters, 38(16):2976–2979, January J I Cirac. Quantum Magnetomechanics with Levitating 2013. Superconducting Microspheres. Physical Review Letters, [31] Thai M Hoang, Jonghoon Ahn, Jaehoon Bang, and Tong- 109(14):147205–5, October 2012. cang Li. Electron spin control of optically levitated nan- [18] M Cirio, G K Brennen, and J Twamley. Quantum Mag- odiamonds in vacuum. Nature Communications, 7:1–8, netomechanics: Ultrahigh- Q-Levitated Mechanical Os- July 2016. cillators. Physical Review Letters, 109(14):147206–5, Oc- [32] Alexander Kuhlicke, Andreas W Schell, Joachim Zoll, tober 2012. and Oliver Benson. Nitrogen vacancy center fluores- [19] Marcus W Doherty, Neil B Manson, Paul Delaney, Fe- cence from a submicron diamond cluster levitated in dor Jelezko, Jörg Wrachtrup, and Lloyd C L Hollen- a linear quadrupole ion trap. Applied Physics Letters, berg. Physics Reports. Physics Reports, 528(1):1–45, 105(7):073101, August 2014. July 2013. [33] T Delord, L Nicolas, L Schwab, and G Hétet. Electron [20] J M Taylor, P Cappellaro, L Childress, L Jiang, D Bud- spin resonance from NV centers in diamonds levitating ker, P R Hemmer, A Yacoby, R Walsworth, and M D in an ion trap. New Journal of Physics, 19(3):033031–11, Lukin. High-sensitivity diamond magnetometer with March 2017. nanoscale resolution. Nature Physics, 4(10):810–816, [34] I Alda, J Berthelot, R A Rica, and Romain Quidant. September 2008. Trapping and manipulation of individual nanoparti- [21] J R Maze, P L Stanwix, J S Hodges, S Hong, J M Taylor, cles in a planar Paul trap. Applied Physics Letters, 7

109(16):163105–5, October 2016. [49] Jan Gieseler, Lukas Novotny, and Romain Quidant. [35] Jen-Feng Hsu, Peng Ji, Charles W Lewandowski, and Thermal nonlinearities in a nanomechanical oscillator. Brian D’Urso. Cooling the Motion of Diamond Nanocrys- Nature Physics, 9(12):806–810, November 2013. tals in a Magneto- Gravitational Trap in High Vacuum. [50] M H Abobeih, J Cramer, M A Bakker, N Kalb, Scientific Reports, 6:30125, July 2016. M Markham, D J Twitchen, and T H Taminiau. One- [36] M C O’Brien, S Dunn, J E Downes, and J Twamley. second coherence for a single electron spin coupled to a Magneto-mechanical trapping of micro-diamonds at low multi-qubit nuclear-spin environment. Nature Commu- pressures. Applied Physics Letters, 114(5):053103–6, nications, 9:2552, June 2018. February 2019. [51] N Bar-Gill, L M Pham, A Jarmola, D Budker, and R L [37] Jan Gieseler, B Deutsch, Romain Quidant, and Lukas Walsworth. Solid-state electronic spin coherence time ap- Novotny. Subkelvin Parametric Feedback Cooling of a proaching one second. Nature Communications, 4:1743– Laser-Trapped Nanoparticle. Physical Review Letters, 6, April 2013. 109(10):103603, 2012. [52] Alexander A Kordyuk. Magnetic levitation for hard su- [38] P Huillery, T Delord, L Nicolas, M Van Den Bossche, perconductors. Journal of Applied Physics, 83(1):610– M Perdriat, and G. Hetet. Spin-mechanics with levitating 612, June 1998. ferromagnetic particles. arXiv.org. [53] J Prat-Camps, C Teo, C C Rusconi, W Wieczorek, and [39] R Skomski. Nanomagnetics. Journal of Physics: Con- O Romero-Isart. Ultrasensitive Inertial and Force Sen- densed Matter, 15(20):R841–R896, May 2003. sors with Diamagnetically Levitated Magnets. Physical [40] R A Norte, J P Moura, and S Gröblacher. Mechan- Review Applied, 8(3):034002, September 2017. ical Resonators for Quantum Optomechanics Experi- [54] Carlos Gonzalez-Ballestero, Jan Gieseler, and Oriol ments at Room Temperature. Physical Review Letters, Romero-Isart. Quantum Acoustomechanics with a Mi- 116(14):021001–6, April 2016. cromagnet. arXiv.org, July 2019. [41] Christoph Reinhardt, Tina Müller, Alexandre Bourassa, [55] D D Bhaktavatsala Rao, S Ali Momenzadeh, and Jörg and Jack C Sankey. Ultralow-Noise SiN Trampoline Res- Wrachtrup. Heralded Control of Mechanical Motion by onators for Sensing and Optomechanics. Physical Review Single Spins. Physical Review Letters, 117(7):077203–8, X, 6(2):021001–8, April 2016. August 2016. [42] Y Tsaturyan, A Barg, E S Polzik, and A Schliesser. Ul- [56] Derek F Jackson Kimball, Alexander O Sushkov, and tracoherent nanomechanical resonators via soft clamping Dmitry Budker. Precessing Ferromagnetic Needle Mag- and dissipation dilution. 12(8):776–783, June 2017. netometer. Physical Review Letters, 116(19):190801, May [43] A H Ghadimi, S A Fedorov, N J Engelsen, M J Bereyhi, 2016. R Schilling, D J Wilson, and Tobias J Kippenberg. Elas- [57] L Thiel, D Rohner, M Ganzhorn, P Appel, E Neu, tic strain engineering for ultralow mechanical dissipation. B Müller, R Kleiner, D Koelle, and P Maletinsky. Quan- Science, 360(6390):764–768, May 2018. titative nanoscale vortex imaging using a cryogenic quan- [44] Jonilyn G Longenecker, H J Mamin, Alexander W Senko, tum magnetometer. 11(8):677 –681, May 2016. Lei Chen, Charles T Rettner, Daniel Rugar, and John A [58] Matthew Pelliccione, Alec Jenkins, Preeti Marohn. High-Gradient Nanomagnets on Cantilevers for Ovartchaiyapong, Christopher Reetz, Eve Emmanouili- Sensitive Detection of Nuclear Magnetic Resonance. ACS dou, Ni Ni, and Ania C Bleszynski Jayich. Scanned Nano, 6(11):9637–9645, November 2012. probe imaging of nanoscale magnetism at cryogenic [45] R Fischer, D P McNally, C Reetz, G G T Assumpcão, temperatures with a single-spin quantum sensor. Nature T Knief, Y Lin, and C A Regal. Spin detection with Nanotechnology, 11(8):700 –705, May 2016. a micromechanical trampoline: towards magnetic reso- [59] Schlussel Yechezkel, Lenz Till, Rohner Dominik, Bar- nance microscopy harnessing cavity optomechanics. New Haim Yaniv, Bougas Lykourgos, Kieschnick Michael, Journal of Physics, 21(4):043049–14, April 2019. Rozenberg Evgeny, Thiel Lucas, Waxman Amir, [46] J-P Tetienne, L Rondin, P Spinicelli, M Chipaux, T De- Maletinsky Patrick, Budker Dmitry, and Folman Ron. buisschert, J-F Roch, and V Jacques. Magnetic-field- Wide-Field Imaging of Superconductor Vortices with dependent photodynamics of single NV defects in dia- Electron Spins in Diamond. Physical Review Applied, mond: an application to qualitative all-optical magnetic 10(1):1, September 2018. imaging. New Journal of Physics, 14(10):103033–16, Oc- [60] C C Rusconi, V Pöchhacker, K Kustura, J I Cirac, and tober 2012. O Romero-Isart. Quantum Spin Stabilized Magnetic Lev- [47] L Rondin, J-P Tetienne, T Hingant, J-F Roch, itation. Physical Review Letters, 119(16):167202, Octo- P Maletinsky, and V Jacques. Magnetometry with ber 2017. nitrogen-vacancy defects in diamond. Reports on [61] C C Rusconi, M J A Schuetz, Jan Gieseler, M D Lukin, Progress in Physics, 77(5):056503, May 2014. and O Romero-Isart. Hybrid architecture for engineer- [48] P Forn-Díaz, L Lamata, E Rico, J Kono, and E Solano. ing magnonic quantum networks. Physical Review A, Ultrastrong coupling regimes of light-matter interaction. 100(2):022343, August 2019. Reviews of Modern Physics, 91(2):025005, June 2019.