Shelvingstyle QND phononnumber detection in quantum optomechanics
Abstract
We propose a new method for optomechanical quantum nondemolition (QND) detection of phonon number, based on a “shelving” style measurement. The scheme uses a twomode optomechanical system where the frequency splitting of the two photonic modes is nearresonant with the mechanical frequency. The combination of a strong optical drive and the underlying nonlinear optomechanical interaction gives rise to spinlike dynamics which facilitate the measurement. This approach allows phonon number measurement to be accomplished parametrically faster than in other schemes which are restricted to weak driving. The ultimate power of the scheme is controlled by the size of the single photon optomechanical cooperativity.
I Introduction
Quantum optomechanical systems, where photons interact with mechanical motion, have been the recent subject of intense experimental and theoretical activity (1). They have enabled the exploration of a wide range of effects, including the cooling of macroscopical mechanical modes to near their quantum ground states (2); (3); (4) and the generation of squeezed states of light (5); (6); (7) and mechanical motion (8); (9); (10).
A key unrealized goal in optomechanics is the quantum nondemolition (QND) measurement of mechanical phonon number, something that would directly reveal the energy quantization of the mechanics (11); (12). The most promising proposals involve optomechanical systems where a single mechanical resonator couples to two photonic modes, as depicted schematically in Fig. 1a. Such systems can in principle be used for phononnumber detection both in the regime where the mechanical frequency is much smaller than the splitting between the optical modes (13); (14); (15); (16); (17), and in the regime where it is comparable to this splitting (18); (19); (20). Twocavity optomechanical systems have been realized experimentally using a variety of platforms (13); (21); (22).
In this work, we revisit the twomode system of Fig. a, and focus on the regime where the mechanics are almost resonant with the frequency difference of the photonic modes, and where one photonic mode is strongly driven. This strong drive generates a coherent, manyphoton optomechanical interaction which cannot be treated perturbatively. We show that this regime allows a new approach to phonon number measurement, one which offers advantages over previous proposals (namely a much faster measurement). The regime we study differs from that of Ludwig et al. (18), who considered more nonresonant interactions and weak driving, such that a pertubative treatment was possible. It also contrasts from the work of Kómár et al. (19), who considered a perfectly resonant system and weak optical driving, but did not consider phonon number measurement physics.
Our measurement scheme is sketched in \creffig:shelvingscheme, and is conceptually similar to the highly successful electron shelving technique used in trappedion systems (23); (24); (25), as explained in \crefsec:shelvinganalogy. The measurement starts by turning on the strong optical drive on the “primary” mode. This induces coherent oscillations between phonons and photons in the undriven “auxiliary” mode, oscillations which necessarily conserve the total number of auxmode and mechanical excitations. The remaining nonlinear interaction then allows one to use the output field of the driven optical mode to measure without changing its value  a true QND measurement. This is then effectively a measurement of mechanical phonon number, as without spurious dissipative effects, is simply equal to the initial phonon number.
As with other schemes, our approach still requires a relatively strong singlephoton optomechanical coupling (i.e. it must be large compared to the geometric mean of the two photonic damping rates) to resolve mechanical Fock states. However, it allows measurement rates that are parametrically larger than the perturbative regime considered in Ref. (18), or in the adiabatic regime of a positionsquared measurement considered in Ref. (13); (14); (15). This faster rate could be of significant utility in new experimental designs of twocavity optomechanical systems where one photonic mode is explicitly engineered to have extremely low dissipation (21); (22). In such systems, the ultimate limit to Fock state detection comes from the requirement that the measurement rate be faster than the heating rate of the mechanics due to their intrinsic dissipation. As we show, the relevant parameter controlling the utility of our scheme in this limit is the singlephoton cooperativity (26),
The remainder of the paper is arranged as follows. In \Crefsec:model we introduce the theoretical model of our twomode optomechanical system, and discuss the parameter regime of interest. In \Crefsec:pseudospin we introduce a mapping to an effective spin system that provides a convenient way to describe the system dynamics; we focus here on the ideal case where the only dissipation is due to the coupling to the input/output port used for the measurement. In \Crefsec:meas we discuss the measurement protocol. Finally, in \Crefsec:gamma we study the effects of the unwanted dissipation both in the aux mode and in the mechanics, and the limits that they impose on QND measurement.
Ii Model
The optomechanical system of interest (c.f. Fig. 1a) consists of two photonic modes and a single mechanical mode. Its coherent dynamics is described by the Hamiltonian
(1) 
The first term describes the two photonc eigenmodes in the absence of optomechanical coupling, and is ( throughout)
(2) 
where are annihilation operators for the two photonic modes (frequencies and respectively). We call these the primary () and auxiliary () modes.
The second term in Eq. (1) describes a mechanical mode with frequency and annihilation operator :
(3) 
We have introduced the parameter which represents the mismatch between the photonic mode splitting and the mechanical frequency. Similar to reference (18), we will be interested in the regime where the mechanical frequency is similar in magnitude to , and hence .
Finally, the last term in \crefeq:Htot describes the optomechanical coupling, which takes the form of a mechanicallymediated tunneling interaction between the two modes:
(4) 
where is the singlephoton optomechanical coupling amplitude. For the specific case of a membraneinthemiddle style optomechanical system (13), a detailed derivation of this model is given in Refs. (14); (27); its limitations are discussed in Refs. (27); (28).
We next use a standard treatment to describe the effects of dissipation on the system (29); (30). We assume that each optical mode is coupled to an independent zerotemperature Markovian reservoir, giving rise to an amplitude damping rate of for , respectively. The mechanics are similarly coupled to a thermal Markovian reservoir characterized by a thermal occupancy , leading to an amplitude damping rate of for . Note that for a membraneinthemiddle type system, treating the optical modes as seeing independent reservoirs can be problematic (17), as it neglects a possible dissipationinduced coupling of the optical normal modes. Such a coupling will not play a role in our scheme, as it is suppressed by the optical mode splitting , while the main optomechanical coupling is resonant.
Our measurement scheme will involve applying a strong coherent drive to the primary, , mode near its resonance frequency . For simplicity, we assume the case of a perfectly resonant drive, as qualitatively similar results are obtained for a detuned drive. In this case, it is useful to move to a rotating frame, and also make a standard displacement transformation on the optical modes:
(5) 
Here, is average primarymode amplitude induced by the drive in the absence of optomechanical coupling, and is the displaced annihilation operator for this mode in the rotating frame.
In this frame we find where
(6)  
(7) 
Finally, we assume that the optical mode gap is sufficiently large that we can safely make a rotating wave approximation and drop the nonresonant interactions described by . This requires in practice . We are left with the interactions depicted in Fig. 1a: a beam splitter interaction between the mechanics and the auxiliary mode, and a nonlinear interaction involving all three modes.
Note that in the absence of driving, i.e. , the photonic modes are only driven by vacuum noise, and the nonlinear optomechanical interaction in \crefeq:Heff effectively vanishes. Therefore, without driving the mechanics are completely decoupled from the photonic modes. We imagine that before the measurement is turned on by driving the primary photonic mode, the mechanics are first prepared close to the ground state, e.g. via standard cavity cooling techniques. The scheme we describe in what follows then allows one to measure the mechanical phonon number.
Iii Dynamics and effective spin description
iii.1 Conserved excitation number
The nonlinear interaction in \crefeq:Heff has the form of a standard threewave mixing Hamiltonian, describing e.g. a nondegenerate parametric amplifier with a dynamical pump mode. It has three conserved excitation numbers (typically known as ManleyRowe constants of motion (32)) reflecting the SU(1) and SU(1,1) symmetries of the interaction. Including the beamsplitter interaction reduces the symmetry of the Hamiltonian, and has only a single conserved excitation number :
(8) 
Note that as is independent of the primary mode, it continues to be conserved even if this mode is damped.
In what follows, we first analyze the ideal case where the only dissipation in the system is due to its coupling to the inputoutput port used for the measurement. This corresponds to , while . Starting with this limit will let us present a clear picture of how our proposed phonon number measurement operates: the conservation of in this limit is at the heart of the shelving scheme depicted in \creffig:diagramb. In this ideal case, we find it is always possible to determine the initial mechanical phonon number by simply measuring for long enough; this is akin to standard shelving measurements.
The regime is also relevant experimentally. While mechanical damping rates are typically orders of magnitude smaller than photonic damping rates, it might first seem surprising that one could achieve . However, recent experiments with microwavecircuit optomechanics have achieved precisely this kind of regime through clever experimental design.
We analyze the effects of nonzero in \crefsec:meas. As expected, they limit the maximum possible measurement time, putting constraints on whether phonon number detection is achievable.
iii.2 Mapping to an effective spin system
To understand the dynamics of our system, it is helpful to make the conservation of explicit by representing the mechanical and auxmodes by an effective spin (a Schwinger boson representation in reverse). We thus introduce the effective (dimensionless) angular momentum operators via
(9) 
These operators satisfy standard angular momentum commutation relations. A direct calculation shows:
(10) 
The conservation of now just corresponds to our effective spin having a fixed length.
With these operators, the effective Hamiltonian is
(11) 
where the effective magnetic field is
(12) 
with the unit vectors in the respective directions. As the mode is damped, we see that takes the form of a generalized spinboson model. The quadratic terms in the original version of now correspond to a Zeeman field on our spin. The order nonlinear term in now corresponds to a systembath coupling where a bath boson can be created while adding angular momentum in the direction to the spin.
The HeisenbergLangevin equations of motion take the form
(13)  
(14) 
where . The operator describes the vacuum fluctuations entering the system from the inputoutput line coupled to the principal mode (i.e. from the inputoutput waveguide that will be used for measurement); it has zero mean and correlation functions
(15) 
When the optomechanical drive is initially turned on, , where is the initial mechanical phonon number. From Eq. (10), we see that the phonon number manifests itself in the dynamics by directly setting the size of the effective spin . The nonlinear interaction creates photons in the principal mode in a way that depends on , and hence on the magnitude of . Thus, by monitoring the output field from this mode, one can determine the initial phonon number. We stress that the “backaction” from the this measurement only alters the direction our effective spin points in, but not its size. The measurement is thus QND.
For easy reference, a synopsis of the mapping from bosons to spins is provided in \creftab:spincomp.
Spin term  Optomechanical Equivalent  

The projection of the effective spin represents the difference between the number of phonons and number of auxmode quanta.  
,  ,  
The manyphoton optomechanical interaction exchanges photons and phonons, and hence acts like a magnetic field in direction.  
The mismatch between the photonic mode splitting and the mechanical frequency acts as an field. 
iii.3 Dynamics in the limit
We consider the system dynamics in the typical to experiment, . We begin by solving \crefeq:eomspinc to find
(16) 
where .
We see that depends on the behaviour of at earlier times, and that the range of this memory effect is . We can solve for the spin dynamics during this short interval to a good approximation by ignoring the effects of , as
(17) 
As long as is at most , we can drop the term as it is proportional to the small parameter
(18) 
With this approximation, the solution for takes the simpler form:
(19) 
Here, the nonHermitian operator is linear in the components of , and given by
(20) 
The initial transient behaviour of is described by .
eq:ctchi implies that the primary mode (and its output optical field) will be sensitive to components of in a manner that depends on both and . In the extreme Markovian limit we find , implying that the is able to measure both transverse components of our effective spin. In contrast, in the limit , the mode is only able to respond to the nonrotating component of the spin, and hence .
It is also instructive to consider the dynamics in terms of the reduced density matrix describing the spin. The interaction picture density matrix is defined in terms of the Schrödingerpicture density matrix via
(21) 
We again consider the limit, and derive a weakcoupling master equation, retaining terms to leading order in . For simplicity, we also consider the typical limit where , allowing us to drop rapidly oscillating terms (i.e. make a secular approximation).
We use to denote the projection of the spin along the axis of the effective magnetic field , and use to denote raising and lowering operators for . Defining as the standard Lindblad superoperator, we find
(22) 
This has the form of a standard master equation for a spin coupled to a bath. The first term describes a modification of the coherent dynamics, described by
(23) 
The remaining terms describe dephasing (rate ), and spin raising / lowering transitions along the axis (rates ). The dissipative rates are
(24) 
Note that unless , the vacuum noise driving the mode is able to cause spin flips upwards in energy: this is possible because we are strongly driving the mode, and can be viewed as being analogous to the phenomenon of “quantum heating” (33); (34); (35). Note also that we always assume that the drive frequency is much larger than .
From \crefeq:Jmaster we find that the steady state of the spin is simply a thermal state,
(25) 
with an effective temperature
(26) 
Again, unless , the steady state will not correspond to the spin being in its ground state (due to the previously mentioned quantum heating physics). Note that in the case where the effective field is completely in the direction (i.e. the detuning vanishes), our master equation describes an effective infinite temperature, and the spin will be completely depolarized in the steady state.
fig:Jbehave illustrates the dissipative spin dynamics described above (see \crefsec:numerics for details on numerics).
iii.4 QND nature of the measurement
As is a conserved quantity, it is clear that our protocol allows for a QND measurement of : this quantity commutes with the system Hamiltonian, even when we include the coupling to the inputoutput port used for the measurement.
Perhaps less obvious is that this also corresponds to a QND measurement of the phonon number in the mechanics. As already discussed, when the control drive is first turned on, the value of corresponds to the initial mechanical phonon number . For a true QND measurement of phonon number, one would like the final phonon number to return to this value after the measurement is complete. This is easily achieved by turning off the large control drive applied to the primary mode. Once this drive is off, , the nonlinear interaction (last term in \crefeq:Heff) keeps acting until , i.e. until all the original phonons are returned to the mechanical resonator.
Iv Description of the measurement output
iv.1 Homodyne measurement and measurement time in the absence of spurious dissipation
Having discussed its dynamics, we now turn to how information on the spin system (and hence the initial photon number) emerges in the output field leaving the primary phononic mode . \crefeq:ctchi tells us that the output field will be sensitive to particular components of the spin vector . Recall that in our original lab frame, the primary mode is driven at its resonance frequency . Information on our effective spin will be optimally encoded in one quadrature of the primary mode output field. We thus consider a homodyne measurement which directly probes this optimal output field quadrature, again starting with the ideal case where .
From standard inputoutput theory (and working in the interaction picture defined by Eq. (5)), the cavity output field leaving the mode is given by
(27) 
and the phonon number can be obtained from measuring an output quadrature,
(28) 
where the angle parameterizes the choice of quadrature.
Using the result in \crefeq:ctchi, we find that the average value of this output quadrature is given by:
(29) 
In the steady state, our spin relaxes to a thermal state, and in general the value of will be nonzero and reflect the overall size of the spin. It will thus be related to the initial mechanical phonon number, as illustrated in \creffig:chiJ.
One finds from \crefeq:chi,eq:rhoJss that in the steadystate is purely real, implying that the optimal choice of quadrature corresponds to an angle . This corresponds to measuring the phase quadrature of the output light leaving the primary mode.
The simplest way to extract the phonon number is to first turn on the measurement by driving the primary photonic mode, and then integrating the output homodyne current (which is proportional to ) for a time . For sufficiently long integration times, the signal of phonon number in the average homodyne current will be resolvable above the noise. For weak couplings, the noise in the homodyne current is simply due to vacuum noise and is Gaussian. As a result, standard calculations (30); (14) indicate that the mechanical Fock state can be distinguished from Fock state when the measurement time is sufficiently long, namely
(30) 
where
(31) 
and the prefactor is determined by
(32) 
As one can see from Fig. 3, for a broad range of , the prefactor is of order unity, and the timescale of the measurement is set by as given in \crefeq:tmeas.
We stress that the time scale required to distinguish between different mechanical Fock states in our setup is parametrically shorter than that obtained in the proposal of Ludwig et al. (18). Their scheme requires a measurement time which scales as to resolve two adjacent Fock states. However, the setup requires an interaction that is sufficiently nonresonant such that hybridization of phonons and photons is weak, and such that a perturbative treatment is applicable. In practice, this requires both the conditions and , yielding a measurement time which scales inversely with a small parameter. Our approach is not constrained in the same way, and can operate much closer to perfect resonance, where . As a result, for the same value of , our approach yields greatly enhanced measurement speed. This in turn implies that an equally rapid measurement can be achieved with a smaller optomechanical coupling: where Ludwig et al. use to resolve phonon number states after (see Ref. (18) Fig. 3), \crefeq:tmeas implies our shelving protocol could do the same in a system where .
iv.2 Optimal parameters
Our system has a fairly large number of tunable parameters. In particular, both nonzero components of the effective field can be selected for optimal performance. We remind the reader that is the manyphoton coupling strength and the mismatch between the mechanical frequency and photonic mode splitting.
A key parameter to optimize the size of the signal in the longtime homodyne current which encodes the initial phonon number. From \crefeq:Xval, this is directly determined by the steadystate value of . In steady state, we obtain from \crefeq:chi,eq:rhoJss,
(33) 
where is the thermal average of taken at a temperature given by \crefeq:Teff.
It follows from these expressions that to have an appreciable signature of the phonon number in the steadystate value of we need and to be comparable. A finite value of (i.e. manyphoton coupling ) is needed to ensure that the steady state has a transverse component that can be detected by the primary photonic mode, as this mode couples to . A nonzero value of (i.e. mismatch between mechanical frequency and photonic mode splitting) is needed so that the effective temperature describing the steadystate is finite and the spin is not in a completely depolarized state.
To make the above discussion more concrete, consider the simple case where we have one phonon intially, and our effective spin corresponds to a spin . There,
(34) 
fig:chiJ shows the steady state value of for other values of initial phonon number. We find the optimal ratio, for a system with one or a few initial phonons, is
(35) 
While the steadystate homodyne signal is not sensitive to the size of (but rather only to the ratio ), the overall magnitude of the effective field is nonetheless important to our scheme. If , then it follows from \crefeq:Jmaster,eq:effrates that the relaxation of the spin to its steady state is greatly suppressed. In this limit, there is simply very little density of states in the primary mode for transitions that involve changing the energy of our spin. We find that it is generally optimal to have the spin reach the steady state on a time scale much shorter than , requiring . That is, one should avoid the manyphoton optomechanical strong coupling regime.
Finally, to have the simple shelving dynamics depicted in Fig. 1b, the coherent oscillations between and b should be underdamped, leading to the requirement . This condition is easy to meet in experiment given the typical weakness of singlephoton optomechanical coupling strengths. Putting these conditions together, we find that the optimal regime for phonon number measurement requires:
(36) 
An alternative approach to the measurement which results in faster measurement rates but is more complicated to implement experimentally is discussed in \crefsec:Gmeas.
iv.3 Resolving Fock states
We next present numerical results for the measurement output of our scheme based on simulations of the full master equation dynamics, including . We construct an estimator for phonon number from the integrated homodyne current:
(37) 
This estimator has been constructed so that at long times
(38) 
where is the initial phonon number, and can be calculated from the effective temperature above, and goes to for .
The behavior of the measurement output, expressed in terms of the phonon number estimator , are presented in \creffig:QND, again for the ideal case of . As expected from our analysis, different phonon numbers are clearly resolvable after an integration time on the order of . The details of these numerical calculations are given in \crefsec:numerics.
Note that in this ideal situation where there is no spurious dissipation, the resolving power of our measurement continues to increase indefinitely as the integration time is increased. This reflects the shelving nature of the dynamics. We now go on to address the additional limitations placed by unwanted dissipation.
V Spurious dissipation: mechanical and auxiliary mode losses
Our analysis so far has treated the ideal case, where the only dissipation is the necessary coupling between the primary photonic mode and the waveguide (or transmission line) used to collect the output field which serves as the measurement record. This idealized analysis gives a useful picture for understanding the more realistic experimental case where there is also mechanical dissipation (damping rate ) and damping of the auxiliary mode (damping rate ).
All three forms of spurious dissipation (mechanical loss, mechanical heating, auxiliary mode loss) cause the quantity to change, disrupting the shelving physics that is at the heart of our scheme. If spurious dissipation changes , the simple correspondence between the measurement record and the initial mechanical phonon number is lost. Focusing on small initial phonon numbers, the fastest rate at which such spurious processes occur is given by . Thus, one requires that the measurement essentially occur on a timescale much shorter than the inverse of this rate. This translates to the condition
(39) 
In \creffig:gamma we show the results of detailed master equation simulations including all forms of spurious dissipation which show this intuitive understand is correct.
It is interesting to consider two relevant limits of the condition in Eq. (39). In the case where is much faster than the mechanical heating , the condition reduces to . This is the same constraint that limits other proposals for measuring phonon number in twocavity optomechanical systems (15); (18).
given recent progress in constructing twomode optomechanical systems with extremely low auxiliarycavity dissipation (21), it is worth considering the opposite limit, where dominates. In this case, our condition reduces to
(40) 
where we have introduced the singlephoton cooperativity . This regime is illustrated in \creffig:gamma. Numerically, we find the requirement is to distinguish between the vacuum state and a single phonon, if the mechanics are the dominant source of spurious dissipation.
Vi Outlook
We have presented and analyzed a new approach for QND phonon number measurement in quantum optomechanics, one which adapts the successful shelving strategy used in trappedion systems. By considering regimes where the driving of the photonic system is strong (and not perturbative), we are able to obtain measurement rates that are parametrically larger than previously considered approaches. At the level of theory, our approach of mapping a bosonic problem onto a spin system could be useful in more complex optomechanical systems where there are again conserved quantities akin to the excitation number in our system.
Our approach is most suited to systems where the auxiliary (unmeasured) photonic mode has an extremely low dissipation rate; recent experiments on twocavity microwavefrequency optomechanical systems (22) present a promising route to this regime.
Acknowledgements
This research was undertaken in part thanks to funding from NSERC and the Canada Research Chairs program.
Appendix A Analogy with electron shelving
The shelving protocol, pioneered by Hans Dehmelt (23); (24); (25), is a form of QND measurement of the state of twolevel atomic system. Given an ion in its ground state () or a metastable excited state (), the protocol makes use of a shortlived auxiliary state (). By applying a resonant laser drive we induce Rabi oscillations between and ; this is accompanied by the spontaneous decay process . The detection of photoemission in the appropriate wavelength is then confirmation that the atom is not in the excited state . However, throughout the process, the atom cannot be said to be in the original ground state . Instead it is in some hybrid state of and . This process is sketched out in figure \creffig:analogyions.
Notice that the rate of spontaneous emission is low, and distinguishing the signal from the noise requires a prolonged measurement similar to what is discussed in \crefsec:meas. Thus the original shelving protocol also constitutes a weak measurement.
The analogy to our measurement scheme is most direct in the case where we are trying to distinguish an initial mechanical phonon number of or . Here, the system is originally in the state or (where denotes a Fock state with the respective number of phonons and auxiliary photons). The application of the drive generates Rabi oscillations between and , while a decay process allows the relaxation with emission of a primary photon, first into the cavity and then out of it through the dissipation channel . By detecting these emitted photons, we can determine the initial phonon number. This process is sketched out in figure \creffig:analogyom.
Our scheme takes the basic shelving concept further, allowing one to distinguish between multiple possible values of while remaining nondestructive. As discussed in the text, information on the phonon number is encoded in the rate at which emitted photons are produced.
Appendix B Optimized measurement prototocol
As discussed in \crefsec:optpar, phonon number measurement in the steady state is subject to two conflicting constraints: the magnitude of is proportional to while the measurable cavity field is coupled to the transverse components, .
These competing requirements on the effective field can be avoided by arranging parameters so that the system has a long relaxation time, and performing the measurement in the transient period before a steady state is reached. We can then align the spin entirely along the axis of measurement and increase the magnitude of the cavity field. As we have seen in \crefeq:effrates, a long relaxation time is achieved in the regime . We therefore suggest the following protocol:

The system is set up with some large energy mismatch . This ensures a slow relaxation to the steady state once the primary mode drive is turned on, the relevant rates being .

With the driving turned off, , the mechanics are prepared in their initial state. As discussed in \crefsec:model, the auxiliary mode remains at zero population in the absence of driving. This corresponds to the spin pointing along .

We then ramp up the primary mode drive semiadiabatically, until . In the spin language, we increase the component until the effective field is nearly parallel with .

If the driving is increased sufficiently slowly, the spin itself will track the direction of the field and finally point along . We make use of the separation of scales to ramp over a time such that
(41) 
Finally, we perform a homodyne measurement as described in \crefsec:meas. Here we use the estimator
(42)
This procedure is illustrated in \creffig:QNDalt. We see that the required measurement time is about three times shorter for this set of parameters.
Appendix C Details of Numerical Simulations
The numerical calculations presented throughout this work are the results of full master equation treating each of the three bosonic modes in our system:
(43) 
where is given in \crefeq:Heff and
(44) 
is the standard Lindblad superoperator.
The numerical simulations were performed using QuTiP (36). We truncated the Hilbert space to , and to for simulations involving spurious dissipation (truncation is not needed if ). We then solved the eigenvalue problem defined by \crefeq:masternum and calculated at any time as a sum of exponentials.
The initial state in all cases was where
(45) 
To simulate the driving of the primary mode, the parameter is initially ramped from zero to its final value. For \creffig:QND,fig:gamma we ramp linearly over a time span of , while in \creffig:QNDalt we increase it exponentially as explained in the figure caption. Experimentally these behaviors can be achieved using pulseshaping techniques.
The standard deviations plotted are given by
(46) 
They were calculated using the quantum regression theorem (29), taking into account the full dynamics of the system and not just cavity vacuum noise.
References
 M. Aspelmeyer, T. J. Kippenberg, and F. Marquardt. Cavity Optomechanics. Springer, Berlin (2014).
 J. D. Teufel et al. Nature 475 359, (2011).
 J. Chan et al. Nature 478 89, (2011).
 R. W. Peterson et al. Phys. Rev. Lett. 116 063601, (2016).
 D. W. C. Brooks et al. Nature 488 476, (2012).
 A. H. SafaviNaeini et al. Nature 500 185, (2013).
 T. P. Purdy, P.L. Yu, R. W. Peterson, N. S. Kampel, and C. A. Regal. Phys. Rev. X 3 031012, (2013).
 E. E. Wollman et al. Science 349, (2015).
 J.M. Pirkkalainen, E. Damskägg, M. Brandt, F. Massel, and M. Sillanpää. Phys. Rev. Lett. 115 243601, (2015).
 F. Lecocq, J. B. Clark, R. W. Simmonds, J. Aumentado, and J. D. Teufel. Phys. Rev. X 5 041037, (2015).
 D. H. Santamore, A. C. Doherty, and M. C. Cross. Phys. Rev. B 70 144301, (2004).
 I. Martin and W. H. Zurek. Phys. Rev. Lett. 98 120401, (2007).
 J. D. Thompson et al. Nature 452 72, (2008).
 A. M. Jayich et al. New J. Phys. 10 095008, (2008).
 H. Miao, S. Danilishin, T. Corbitt, and Y. Chen. Phys. Rev. Lett. 103 100402, (2009).
 A. A. Gangat, T. M. Stace, and G. J. Milburn. New J. Phys. 13 043024, (2011).
 Y. Yanay, J. C. Sankey, and A. A. Clerk. Phys. Rev. A 93 063809, (2016).
 M. Ludwig, A. H. SafaviNaeini, O. Painter, and F. Marquardt. Phys. Rev. Lett. 109 063601, (2012).
 P. Kómár et al. Phys. Rev. A 87 013839, (2013).
 S. BasiriEsfahani, U. Akram, and G. J. Milburn. New J. Phys. 14 085017, (2012).
 T. K. Paraïso et al. Phys. Rev. X 5 041024, (2015).
 L. D. Tóth, N. R. Bernier, A. Nunnenkamp, A. K. Feofanov, and T. J. Kippenberg. arXiv:1602.05180, (2016).
 H. G. Dehmelt. Adv. At. Mol. Phys. 3 53, (1968).
 W. Nagourney, J. Sandberg, and H. Dehmelt. Phys. Rev. Lett. 56 2797, (1986).
 J. Javanainen. Phys. Rev. A 33 2121, (1986).
 M. Aspelmeyer, T. J. Kippenberg, and F. Marquardt. Rev. Mod. Phys. 86 1391, (2014).
 H. K. Cheung and C. K. Law. Phys. Rev. A 84 023812, (2011).
 A. Xuereb and P. Domokos. New J. Phys. 14 095027, (2012).
 C. W. Gardiner and P. Zoller. Quantum Noise. Springer, Berlin/Heidelberg (2004).
 A. A. Clerk, M. H. Devoret, S. M. Girvin, F. Marquardt, and R. J. Schoelkopf. Rev. Mod. Phys. 82 1155, (2010).
 D. F. Walls and G. J. Milburn. Quantum optics. Springer, Berlin/Heidelberg, 2nd edition (1997).
 J. M. Manley and H. E. Rowe. Proc. Inst. Radio Eng. 44 904, (1956).
 M. I. Dykman, M. Marthaler, and V. Peano. Phys. Rev. A 83 052115, (2011).
 V. Peano and M. Thorwart. Phys. Rev. B 82 155129, (2010).
 M.A. Lemonde and A. A. Clerk. Phys. Rev. A 91 033836, (2015).
 J. R. Johansson, P. D. Nation, and F. Nori. Comput. Phys. Commun. 184 1234, (2013).