Autonomous phonon maser in levitated spin-mechanics

TL;DR

This work demonstrates that a microwave-driven and optically pumped NV center can act as a tunable gain medium for the ultra-low-frequency center-of-mass motion of a levitated nanodiamond, enabling autonomous phonon masing. Through adiabatic elimination in a separation-of-timescales regime, the authors derive a reduced mechanical master equation with closed-form, detuning-dependent rates, revealing a sharp lasing threshold governed by the sign of the phonon-number damping $\gamma_{ m eff}^{(n)}(\delta)$. The threshold and saturation follow Maxwell–Bloch predictions, and a semiclassical phase-space analysis shows a phase-diffusing limit cycle that appears as a ring with a circulating coherent component, observable via $g^{(2)}(0)$ and radial distributions. These results connect microscopic spin control to macroscopic nonequilibrium self-oscillations in a high-$Q$, thermally noisy platform and provide a concrete design framework for realizing tunable autonomous oscillators with levitated nanomechanics.

Abstract

Levitated nanodiamonds hosting a single nitrogen-vacancy (NV) center provide an ultra-low-frequency mechanical mode with widely tunable dissipation and spin backaction under microwave dressing and optical pumping. We demonstrate that the driven NV spin can be tuned to act as an inverted gain medium for the center-of-mass motion, thereby stabilizing an autonomous phonon maser. In the separation-of-timescales regime where spin dynamics is fast, adiabatic elimination yields a reduced mechanical master equation with closed-form, detuning-dependent transition rates and a sharp threshold given by the sign change of the phonon-number damping. For representative levitated-NV parameters, we find that a percent-level dressed-basis inversion is sufficient to reach the threshold, and the small-signal gain can exceed the intrinsic mechanical loss by orders of magnitude. Full master-equation simulations confirm above-threshold self-oscillation and a phase-diffusing, coherent steady state, whose saturation follows the Maxwell-Bloch prediction.

Figures (3)

• Figure 1: Schematic of the levitated NV spin--mechanics platform. A nanodiamond hosting a single NV center is levitated above a source of static magnetic-field gradient (red field lines), coupling the spin to the center-of-mass coordinate $z$. The COM motion along the trap axis is approximated as a harmonic oscillator of frequency $\omega_m$. A microwave field at frequency $\omega_d$ drives a selected effective two-level transition of bare splitting $\omega_0$ and dresses the spin; when the dressed splitting $\tilde{\omega} \simeq \omega_m$, the interaction becomes near resonant and enables coherent exchange between spin excitations and phonons.
• Figure 2: Gain, saturation, and phase-diffusing steady state of the levitated-NV phonon maser.(a) Spin-induced contribution to the phonon-number damping, $\Gamma_{\rm opt}(\delta)\equiv \Gamma_-(\delta)-\Gamma_+(\delta)$, normalized to the intrinsic loss $\gamma_m$, versus detuning $\delta/\gamma_2$ for several dressed inversions $S_z^{(0)}$. From Eq. \ref{['eq:Gammaopt']}, $\Gamma_{\rm opt}(\delta)$ is Lorentzian in $\delta$ and changes sign with $S_z^{(0)}$. The horizontal dashed line marks threshold, $\Gamma_{\rm opt}(0)=-\gamma_m$ (equivalently $\gamma_{\rm eff}^{(n)}=0$), giving $S_z^{(\mathrm{th})}=0.016$ for the parameters used. (b) Above threshold ($\delta=0$, $S_z^{(0)}=0.2$), Langevin simulations show exponential growth and saturation of the mean occupation $\langle n(t)\rangle$. The blue dashed line is the Maxwell--Bloch limit-cycle intensity $n_{\rm las}$ Eq. \ref{['app:nlas0']}; the red dashed line is the noise-corrected stationary mean from the associated Fokker--Planck theory, $\langle n\rangle_{\rm FP}=n_{\rm las}+\delta n$ with $\delta n=(\mathrm{Var}(z)+\mathrm{Var}(p))/2$. (c) Stationary radial distribution $P(r)$ of the phase-space amplitude $r=\sqrt{z^2+p^2}=\sqrt{2n}$, peaked near the predicted ring radius $r_0\simeq \sqrt{2n_{\rm las}}$ (vertical dashed line). (d) Equal-time intensity correlations $g^{(2)}(0)=(\langle n^2\rangle-\langle n\rangle)/\langle n\rangle^2$ versus $S_z^{(0)}$, showing the crossover from thermal bunching, $g^{(2)}(0)\approx 2$, below threshold to near-coherent statistics, $g^{(2)}(0)\to 1$, above threshold; vertical dashed line: $S_z^{(\mathrm{th})}$. Parameters are those of Table \ref{['tab:est_inputs']} with $g=\eta\omega_m$ and $\eta=0.1$; the diffusion term uses an effective pre-cooled occupation $\bar{n}_{\rm eff}=5$.
• Figure 3: Phase-diffusing maser steady state in phase space. Phase-space quasiprobability $W(z,p)$ reconstructed from the semiclassical Langevin dynamics associated with the reduced mechanical master equation \ref{['eq:MEred']}. Above threshold the dynamics approaches a limit cycle with a well-defined amplitude but a freely diffusing phase; the steady state is therefore phase-averaged and forms a ring in phase space. The dashed circle marks the Maxwell--Bloch prediction for the limit-cycle radius, $r_0\simeq \sqrt{2n_{\rm las}}$, with $n_{\rm las}$ from Eq. \ref{['app:nlas0']}. The ring thickness quantifies amplitude fluctuations induced by the mechanical bath and spin-induced diffusion, and accounts for the difference between $n_{\rm las}$ and the total steady-state occupation $\langle n\rangle$. Parameters are as in Fig. \ref{['fig2']}; here $\delta= \tilde{\omega}-\omega_m=0$ and $S_z^{(0)}=0.2$.