Quantum cooling and squeezing of a levitating nanosphere via time-continuous measurements

With the purpose of controlling the steady state of a dielectric nanosphere levitated within an optical cavity, we study its conditional dynamics under simultaneous sideband cooling and additional time-continuous measurement of either the output cavity mode or the nanosphere's position. We find that the average phonon number, purity and quantum squeezing of the steady-states can all be made more non-classical through the addition of time-continuous measurement. We predict that the continuous monitoring of the system, together with Markovian feedback, allows one to stabilize the dynamics for any value of the laser frequency driving the cavity. By considering state-of-the-art values of the experimental parameters, we prove that one can in principle obtain a non-classical (squeezed) steady-state with an average phonon number $n_{\sf ph}\approx 0.5$.

With the purpose of controlling the steady state of a dielectric nanosphere levitated within an optical cavity, we study its conditional dynamics under simultaneous sideband cooling and additional time-continuous measurement of either the output cavity mode or the nanosphere's position. We find that the average phonon number, purity and quantum squeezing of the steady-states can all be made more non-classical through the addition of time-continuous measurement. We predict that the continuous monitoring of the system, together with Markovian feedback, allows one to stabilize the dynamics for any value of the laser frequency driving the cavity. By considering state-of-the-art values of the experimental parameters, we prove that one can in principle obtain a non-classical (squeezed) steady-state with an average phonon number n ph ≈ 0.5.

I. INTRODUCTION
Bringing physical degrees of freedom to the quantum regime is proving so difficult that quantum control is bound to be a multi-branched endeavour, where techniques developed on different platforms and designed for different aims are blended together. In the context of cooling matter to the quantum ground state, a primary directive of quantum control, various techniques have come to the fore over the last twenty years. Prominent among them in the case of macroscopic mechanical systems is sideband cooling, where the targeted degree of freedom is driven by light on a red sideband of an optical cavity, such that a beam-splitting light-matter interaction is achieved and excitations are drained out of the system, cooling it down. On the other hand, the implementation of efficient non-demolition quantum measurements is another obvious way to extract entropy from a quantum system. This study evaluates the combined performance of sideband laser cooling and continuous quantum measurements (also known as 'monitoring') on a levitating nanosphere, an interesting optomechanical system where both such techniques are applicable to actual experiments.
The research field of quantum opto-mechanics, whose goal is to achieve control at the quantum level of massive mechanical oscillators, has received increasing attention in the last years, both for applicative and fundamental reasons [1]. So far, the main objective pursued by both theorists and experimentalists is the cooling of the oscillator either to its motional ground state or to non-classical states, such as low-number Fock states or squeezed states; several protocols have been proposed in this respect, adapted to different physical settings and using different control strategies [2][3][4][5][6][7][8][9][10][11].
Here, we will focus on a particularly promising setup, where the opto-mechanical system corresponds to a dielectric nanosphere trapped inside an optical cavity [12][13][14][15]. As the nanosphere is levitating, the coupling to the environment is minimised; it is possible then to neglect the thermal background of phonons which is typically one of the most detrimental sources of noise in opto-mechanical systems. Varia- * marco.genoni@ucl.ac.uk tions of this basic paradigm, where the nanosphere is trapped by the cavity field only [13,16], by an optical tweezer within the cavity field [17], or with the help of an electromagnetic trap [15], have been recently proposed. On the theoretical side, a detailed derivation of the master equation for the quantum state corresponding to the nanosphere's motion and to the cavity mode may be found in [18,19], and allows one to study the time behaviour of this opto-mechanical system, as well as its steady-state properties. As a matter of fact, this master equation paves the way to the analysis of protocols combining time-continuous measurement and feedback operations, and it will be the starting point of our study.
We have witnessed constant progress in the understanding of quantum filtering, i.e. of the conditional dynamics of quantum systems subjected to time continuous measurements [20][21][22][23][24][25][26][27][28][29][30][31][32][33][34][35]. In particular, as regards systems described by continuous, canonical degrees of freedom, diffusive dynamics described by multivariate Wiener increments have been characterised in detail, and a general framework is available [36,37]. Such diffusive dynamics correspond to conditional evolutions due to the monitoring of the environment through the class of so called general-dyne detections [21]. These quantum measurements amount to performing homodyne detections on the environment and, possibly, additional ancillary modes which are coupled to the environment itself via Gaussian unitary transformations [38,39]. We recall that the term "homodyne" detection refers to the projective measurement on the eigenbases of canonical position and momentum operatorsx andp. In the case where the overall Hamiltonian is quadratic in the canonical operators, and the system is linearly coupled to the environment, the conditional dynamics due to general-dyne detections preserves the Gaussianity of the quantum state and thus the whole dynamics can be equivalently described by the evolution of first and second statistical moments only. In this case, it is easy to optimise different steady-state properties, such as entanglement, squeezing and purity, over the parameters characterising the detection scheme [40][41][42][43][44][45].
The role of conditional dynamics due to indirect measurements, and of subsequent quantum feedback, for steering the quantum state of a mechanical oscillator towards either its ground state or a certain non-classical state has been already discussed in the literature. In [46], Doherty and Jacobs derive arXiv:1503.05603v1 [quant-ph] 18 Mar 2015 the effective stochastic master equation corresponding to position measurements obtained through adiabatic elimination of a cavity mode continuously monitored through homodyne detection, and discuss the related feedback strategies aimed cooling the motion of the oscillator. In [47], the unravelling corresponding to direct position measurements of an oscillator interacting with a non-zero temperature thermal bath is investigated; in particular, the results obtainable in different regimes, corresponding to different measurement resolutions, are discussed in great detail, and the ensuing stochastic master equation has been seminal in the design of protocols to engineer thermo-mechanical squeezing [11,48]. Also, the usefulness of discrete and repeated measurements on a coupled qubit in a hybrid setup has been investigated in [49], where an effective dynamics able to prepare a squeezed steady-state has been identified. Very recently, Hofer and Hammerer [50] have studied the effect of continuous homodyne detection on the cavity output combined with sideband cooling in a standard optomechanical setup, where the mechanical oscillator interacts with a non-zero temperature thermal bath. First, they consider a single oscillator and discuss the corresponding steady-state average number of phonons; then, they present more complex and sophisticated protocols able, for example, to create entanglement between two distant oscillators.
As already stated above, in this manuscript we focus on the case of a levitating dielectric nanosphere in an optical cavity, described by the master equation derived in [19]. Notice that our treatment is distinct from the existing literature in that the master equation of the levitating nanosphere includes a photon scattering term, and the measurements of both cavity output, via homodyne detection, and oscillator position, through the light scattered by the nanosphere itself, are considered simultaneously. By combining these measurements with sideband cooling and Markovian feedback, we address the possibility of both cooling the oscillator towards its ground state and of generating quantum mechanical squeezing, that is sub-vacuum fluctuations, which is a paradigmatic signature of non-classicality useful for quantum metrology and precision sensing [51]. Finally, we also address in more detail the experimental setup described in [16], where the nanosphere is trapped in a high finesse optical cavity: by considering stateof-the-art values for the experimental parameters and for the measurement efficiencies, we show the possibility to vastly improve the performances of sideband cooling both in terms of steady-state average number of phonons and in terms of generation of squeezed quantum states.
The manuscript is organised as follows: in Sec. II we discuss the master equation of a levitating nanosphere, and introduce the notation and figures of merit that will be discussed in the remainder of the article. In Sec. III we introduce the stochastic master equation describing the time-continuous measurements and then present the results obtainable for different values of their measurement efficiencies. In Sec. IV we discuss the performances of these protocols for a specific experimental setup, while Sec. V concludes the paper with some final remarks.

II. LEVITATING DIELECTRIC NANOSPHERE MASTER EQUATION
We will consider two quantum degrees of freedom; the cavity electromagnetic mode and the mechanical motion of a trapped nanosphere, described respectively by bosonic operators a and b satisfying the commutation relations [a, a † ] = [b, b † ] = 1. We can then define the corresponding position and momentum quadrature operators as By considering the cavity driven by a laser at frequency ω L , the Hamiltonian describing the interaction between the two modes reads where ω m is the mechanical frequency and we have already transformed the Hamiltonian to a frame rotating at the frequency ω L , such that ∆ = ω L − ω c denotes the detuning from the cavity resonance ω c . By considering the open dynamics resulting from the interaction with the environment (i.e. the free electromagnetic modes), one obtains the master equation [19], where The first term is responsible for the unitary dynamics, the second one describes the usual cavity loss, while the third one corresponds to the recoil heating due to photon scattering from the oscillating nanosphere.
By assuming that the system is prepared in a Gaussian state (e.g. in a thermal state) at time t = 0, at every time the dynamics keeps the state Gaussian (see [52] for different reviews on Gaussian states). As a consequence, the whole dynamics can be fully described by means of the first moment vector R and the covariance matrix σ, whose elements are defined as Throughout the article, we will discuss the efficiency of our protocols by considering the effect on the mechanical oscillator properties. Hence, it is useful to introduce the covariance (sub)matrix corresponding to the oscillator quantum state alone, obtained by tracing out the cavity mode, This corresponds to where σ jk are the elements of the global covariance matrix defined in Eq. (4), which in fact correspond to variances and covariances of the quadrature operators x m and p m . All the properties we are interested in can be easily obtained from the matrix σ (m) . We will focus on three figures of merit: quantum squeezing, purity and average number of phonons.
Quantum squeezing can be quantified in a dB scale through the minimum eigenvalue of σ (m) , as dB = 10 log 10 ξ , Negative values of dB corresponds to a quadrature of the mechanical oscillator having sub-vacuum fluctuations, and thus to a non-classical squeezed state. In general this quadrature will be a certain linear combination of position and momentum operators, and the detectability and usefulness of the corresponding squeezing may not be straightforward. This is one of the reasons why, in Sec. IV, we will also focus on the position fluctuations ∆x 2 m only. The usefulness of a certain quantum state in quantum information and quantum communication protocols often strictly depends on its purity, that is on how close a state is to a projector |ψ ψ| on a single Hilbert space vector, rather than a statistical mixture thereof. Such single vector states are also known as pure states. The purity of a quantum state is defined as µ = Tr[ 2 ] and it takes is maximum value µ = 1 if and only if is a pure state. For Gaussian states, the purity can be evaluated through the covariance matrix, and in particular for the single-mode mechanical oscillator we can use the for- Notice that, as we are considering single-mode states, all entropies, including the von Neumann entropy, are moononic functions of the purity, which thus fully characterize the mixedness of the quantum state.
Finally, as we will mainly consider steady-states of the mechanical oscillator having zero first moments (Tr[ x m ] = Tr[ p m ] = 0), the number of phonons can also be evaluated directly from the covariance matrix as Here, tr[·] denotes to the trace of a finite dimensional matrix. Notice that whenever the first moments of the oscillator are not equal to zero, Eq. (7) gives only a lower bound on the actual average number of phonons n ph . In terms of first moment vector and of the global covariance matrix, the master equation for the quantum state in (3) entails the following time evolution: where the matrices A (drift matrix) and D (diffusion matrix) can be easily evaluated from the parameters entering in (3) and in the interaction Hamiltonian (2) [37] and are reported in the Appendix A. In particular, the existence of a steadystate of the dynamics (having zero first moments) can be easily discussed by analyzing the eigenvalues of the drift matrix A: where Re[x] denotes the real part of a complex number x, and α j are the eigenvalues of the drift matrix A. Such a property is often referred to as 'Hurwitz stability' in the control literature. A plot of the stable regions for the master equation (3) as a function of the detuning ∆ and of the coupling constant g is pictured in Fig. 1. As one can see, stability is obtained only for red-detuning (∆ < 0). In particular, in order to reach a steady-state, for larger values of the coupling constant g, one needs also a larger value of the detuning |∆|. In Fig. 2, we  (3), as a function of the laser detuning ∆ and of the oscillator-cavity mode coupling constant g, for κ = ωm/3 and Γ = ωm/10. Notice that, for ∆/ωm = 0, the drift matrix A always has an eigenvalue equal to zero, and thus the system cannot be considered strictly stable.
plot the values of the purity of the oscillator and the average number of phonons obtainable at steady-state as a function of the detuning ∆, where the other parameters are fixed according to the current state of the art in experiments [19]; one can see that for large red-detunings one can cool the oscillator to a state with around n ph = 30 phonons at steady-state.

III. TIME-CONTINUOUS HOMODYNE MEASUREMENT OF CAVITY MODE AND OSCILLATOR POSITION
We now consider the conditional evolution due to generaldyne time-continuous measurement on both the cavity mode and the oscillator. The corresponding stochastic master equation reads where and dw j are uncorrelated Wiener increments, such that dw j dw k = dt δ jk . The term √ η 1 κH[ae iφ ] describes the effect of continuous homodyne on the output cavity mode with efficiency η 1 , where the phase φ can be adjusted by choosing the optical phase of the monitored quadrature operator (e.g. φ = 0 and φ = π/2 correspond respectively to homodyning quadratures x c and p c ) [21,37]. Analogously, the term √ describes the effect of continuous monitoring of the oscillator position, with efficiency η 2 [46,47].
As for the unconditional master equation (3), the dynamics induced by the continuous measurement here considered does not change the Gaussian character of the quantum state; as a consequence we can translate Eq. (11) into equations for the first moment vector and covariance matrix: where dw = (dw 1 , dw 2 ) T and the matrices N, B, A and D can be evaluated starting from the parameters entering the stochastic master equation (11) [37] and are reported in the Appendix A. It is important to observe that the Riccati equation for the covariance matrix is completely deterministic, and yields a steady-state that can be efficiently evaluated numerically. On the other hand, the first moments' evolution is stochastic, i.e. it depends on the outcomes of the continuous measurements. As a consequence, at each time the conditional state is a Gaussian state whose covariances and correlations evolve deterministically according to the dissipative dynamics and the kind of measurement performed, while its first moments evolves randomly in the phase-space, depending on the values of the photocurrents. In the following we will focus on these conditional steady-states only.
Although it is possible to achieve these conditional covariance matrices by pure filtering, i.e. recording the measurement outcomes (photocurrents), in order to remain in the harmonic trap regime, where our treatment applies, it is useful to suppress the drift of the first moments, due to the stochastic evolution, by an active feedback operation. The role of feedback is indeed to use the information contained in the measurement outcomes in order to remove the contribution given by the last term in Eq. (12), which is proportional to the Wiener increment dw. This can always be done by adding a linear feedback term in the Hamiltonian (3) with coupling constants proportional to the photocurrents, i.e.
where r is the vector of quadrature operators introduced in Eq.
(1) and f (t) is an optimized vector of time-dependent coupling constants whose values depends linearly on the continuousmeasurement outcomes [37]. In practice, while for the cavity field this corresponds simply to a linear driving, in the case of a mechanical oscillator it can be obtained by means of a combination of impulses and shifts of the trapping potential (for a more detailed discussion of this issue see [46]). The first important consequence of Eqs. (11) and (13) regards the stability of the opto-mechanical system: we find that, whenever the interaction between the two bosonic modes is on (i.e. for g > 0), if a continuous measurement is performed, i.e. if η 1 > 0 or η 2 > 0, the stochastic master equation satisfies the stability conditions described in [37]. We should remark that this stability condition regards the covariance matrix steady-state, while in principle the first moments could not go to a steady-state value (e.g. to zero). However, as we have just stated above and discussed for example in [46], the information obtained from the measurement can be used to obtain a proper steady-state for the quantum system with zero first moments, as the stochastic drift on the latter may always be canceled by Markovian linear feedback.
In the following, we will concentrate on the the steadystate properties of the harmonic oscilaltor. As anticipated in the previous section, we will analyse the number of phonons, the purity of the state, and at the achievable quantum squeezing, quantified by the minimum eigenvalue of the steady-state covariance matrix. We will consider different measurement strategies: (i) measurement of the cavity mode only (η 1 > 0 and η 2 = 0); (ii) measurement of the oscillator position only (η 1 = 0 and η 2 > 0); (iii) simultaneous measurement of the cavity mode and of the oscillator (η 1 > 0 and η 2 > 0).

A. Time-continuous measurement of the cavity mode
In this subsection we investigate the properties of the steady-state of the oscillator in the case where no measurement is performed directly on the nanosphere, while the output of the cavity is continuously measured. The phase φ of the quadrature which is monitored through homodyne detection is optimized for every set of parameters and for all the figures of merit considered. We notice that the behaviours of these different optimized homodyne phases for squeezing, number of phonons and purity as a function of the detuning parameter are almost identical in all the cases we investigated. In Fig. 3 we plot the steady-state average phonon number, purity and squeezing for a particular (experimentally reasonable [19]) choice of the parameters. Note that very good results are obtained for all the values of the detuning we are considering. This is really important from an experimental point of view as it strongly relaxes the requirement to be sideband resolved in order to cool the nanosphere motion, a condition that in fact is particularly difficult to meet for low frequency oscillators.
As one can notice, if we are interested in cooling the oscillator, both in terms of number of phonons and of purity of the quantum state, the optimal choice is obtained for ∆ ≈ 0, that is near resonance, with the actual optimal point taking a very small negative value. On the other hand, in order to obtain the largest value of quantum squeezing, one may choose a value of the detuning around ∆ = −2.5ω m , obtaining a non-classical state that exhibits around 3dB of squeezing.

B. Adding time-continuous measurement of the oscillator position
As a preliminary analysis, let us consider the effect of continuous monitoring of the oscillator position, whilst the cavity output is left unobserved (η 1 = 0). The results for this case are shown in Fig. 4. We observe that, if we want to minimize the average number of phonons or maximize the purity, the optimal performances are obtained in the case of large detuning |∆| 1, regardless of the red or blue shift of the driving field. This should not come as a surprise, as a large detuning corresponds to decoupling the oscillator from the cavity, and hence directly measuring an isolated degree of freedom. By solving the dynamics for the decoupled mechanical oscillator alone (i.e., for g = 0), one can evaluate analytically the corresponding steady-state. One may prove that its purity simply depends on the measurement efficiency, as µ = √ η 2 , which thus univocally characterize the entropy of the steady-state; on the other hand quantum squeezing and number of phonons do depend also on the ratio between the noise parameter and the mechanical frequency Γ/ω m , and their behaviour is plotted in Fig. 5 for different values of the measurement efficiency. We remind the reader that larger values of Γ correspond to a large amount of scattered light from the nanosphere. On the one hand, this implies a larger amount of incoherent energy acquired by the oscillator due to the recoil heating process; on the other hand, it also corresponds to a large amount of information available for the continuous position measurement, and thus to the possibility to convert such incoherent energy into stead-state quantum squeezing. This explain the different behaviour we observe in Fig. 5 for quantum squeezing and number of phonons. As anticipated above, the results derived here for the decoupled oscillator almost perfectly correspond to the ones reported in Fig. 4 in the case of large detuning. Note that the small values of the number of phonons away from resonance are indeed due to the fact that the oscillator decoherence rate Γ is relatively small with respect to ω m (in the range of 0.1ω m ). For example, for a unit efficiency measure- ment (η 2 = 1), an almost pure quantum state is obtained with around n ph = 0.02 phonons. The state is also squeezed, with a squeezing around 1dB (a value compatible with the number of phonons obtained). Our findings show that, if direct position monitoring with high efficiency were possible, feedback cooling would greatly outperform sideband cooling of the oscillator. However, in practice, a decoupled cavity is not likely to be a favourable condition to work in, as the actual efficiencies of position measurements through scattered light are bound to be severely limited by a number of practical factors (one among all, the geometric impossibility of probing the whole solid angle of scattering). Quite interestingly, at lower values of the measurement efficiency (i.e. η 2 < 0.5) or if our aim is to optimize the squeezing of the steady-state (also with larger values of η 2 ), a combination of sideband cooling and position measurements still yields the best results. In particular, as regards quantum squeezing, the optimal detuning is again around ∆ ≈ −2.5ω m , obtaining a quantum squeezing around 3dB (with only small variations for measurements with efficiency η 2 ∈ [0.5, 1]). Note that, in this section, we are not taking into account the fact that, in principle, varying the detuning will change the number of photons inside the cavity, and thus also the effective coupling constant g (such that larger values of detuning correspond to lower values of g). This effect will be properly taken into account in Sec. IV.
Finally, let us consider the simultaneous monitoring of the cavity output mode (optimizing the phase of the homodyne detection and with unit efficiency: η 1 = 1) and of the oscillator position (with different efficiencies η 2 ). The results for this case are plotted in Fig. 6. The properties of the steady-state are qualitatively similar to the ones we have just discussed, showing the prominent role played by the oscillator position measurement over the other control strategies (i.e. sideband cooling and cavity homodyne measurement). Nevertheless, one observes slightly better results with respect to the unobserved cavity scenario, both in terms of steady-state phonons (in particular if we do not consider the large detuning regime) and in terms of quantum squeezing. For red-detuning with ∆ ≈ −2.5ω m , one obtains the highest value of squeezing of approximately 3.5dB.

IV. RESULTS FOR A NANOPSHERE LEVITATED IN A HIGH FINESSE OPTICAL CAVITY
In this section, we will make specific predictions on what could be achieved by the continuous measurement of a nanosphere levitated by the field of a high finesse optical cavity, as depicted in Fig. 7. The position and dynamics of the nanosphere can be directly measured by collecting the light it scatters, as described in [15], or indirectly monitored through the homodyne monitoring of the light that leaves the optical cavity.
The experimental setup we consider comprises a silica nanosphere of radius r = 200 nm, with mass m = 7.35 · 10 −17 Kg, and a cavity, characterized by a resonance frequency ω c /2π = 2.8 · 10 14 Hz (λ = 1064 nm), length L = 13 mm, finesse F = 400000, and cavity waist w = 60 µm . The corresponding value for the cavity loss parameter appearing in the master equation (3) is κ/(2π) = c/(2FL) = 29 kHz. The cavity average photon number reads As usual in opto-mechanical setups, the cavity-oscillator coupling constant g depends on the cavity average photon number  7. Schematic of the experimental set-up for the quantum control of a levitated dielectric nanosphere within an optical cavity. The nanoparticle can be sideband-cooled by sending in light that is reddetuned from the cavity resonance. The light that leaves the cavity can be continuously monitored to perform generaldyne measurements. The light which is scattered by the particle can be collected to give information upon the position of the nanosphere.
n c = |α| 2 . In particular we adopt the formula where k = 2π/λ, and the volumes of the nanosphere and of the cavity read V s = (4/3)πr 3 and V m = πw 2 L respectively. The remaining parameters, r and ω l denote the dielectric constant and the driving laser frequency. Since the nanosphere is here trapped by the cavity field, the mechanical frequency depends on the average photon number n c too, as per It is important to notice that, because of the dependence of ω m on the average photon number n c , the behaviour of the opto-mechanical coupling constant is modified with respect to standard setups, in that g ∼ n 1/4 c . As far as the oscillator's decoherence rate Γ is concerned, one can follow the results shown in [16]. In particular, one observes that Γ ∼ n 1/2 c , while the corresponding ratio Γ/ω m ≈ 0.15 -evaluated for our specific experimental parametersis fixed for every value of the detuning ∆ [53]. In the case of zero detuning (∆ = 0), we obtain a mechanical frequency ω m0 /2π ≈ 47 kHz and g 0 /2π ≈ 24 kHz.
Like in the previous section, we will consider the steadystate properties as a function of the detuning ∆. As the mechanical frequency ω m decreases for increasing detuning, the corresponding zero-point motion will increase; for this reason besides the quantum squeezing property (which is obtained by considering variances renormalized to the zero-point motion and thus does not take into account the effect of the variation of the frequency ω m ), we will also consider the proper fluctu-ation of the position ∆x 2 m being the fluctuations of the dimensionless position operator x m . Apart from taking into account the different zeropoint fluctuations, ∆X is arguably more interesting from an experimental and practical point of view as it is directly accessible in experiments.
The results are depicted in Figs. 8 and 9. In the first one, we plot the results in the case where the efficiency of the homodyne measurement of the cavity output is maximum (η 1 = 1), the phase of the homodyne is optimized for each figure of merit, and different values of the oscillator's measurement efficiency η 2 are considered. As one should expect, we observe better performances for increasing values of η 2 . In detail, by focusing on the average number of phonons, for small values of η 2 we observe that two local minima of n ph are obtained, one in the red-and one in the blue-detuning regimes. By considering an efficiency over 50%, we observe a local minimum is obtained near resonance, while the optimal performances are obtained in the case of large detuning where, as discussed above, the oscillator's motion decouples from the cavity mode, and thus its steady-state properties become entirely dependent on the position measurement efficiency. If we instead focus on quantum squeezing, in this case too we obtain that for all values of η 2 the optimal working point corresponds to a value of the detuning near ∆ = −ω m0 , that is in the red sideband. Surprisingly, observing the behaviour of the position fluctuations ∆X, its actual minimum is always reached at resonance, while sub-vacuum fluctuations are observed only in the blue sideband, that is for ∆ > 0. One important remark is needed here: as sub-vacuum fluctuations are not obtained for most of the values of the detuning ∆, the amount of quantum squeezing observed in the other plot of Fig. 8 necessarily corresponds to different quadrature operator describing the oscillator; in general, this quadrature will correspond to a certain linear combination of position and momentum of the oscillator, whose usefulness and detectability may not be straightforward.
Similar observations are made also regarding Fig. 9, where we considered state-of-the-art values for the position and homodyne measurement efficiencies, estimated at η 2 = 0.2 and η 1 = 0.9 respectively. In this realistic scenario, we demonstrate that for all the different values of ∆ considered one can achieve a steady-state characterized by n ph < 1, with a minimum reaching n ph = 0.4 and a maximum squeezing of dB ≈ −1.9 obtained in the red sideband.

V. SUMMARY AND CONCLUSIONS
Our study explicitly shows the quantum control possibilities offered, in realistic setups, by the combined simultaneous monitoring of scattered as well as coherent cavity light interacting with a levitating dielectric nanosphere. In particular, it was shown that • time-continuous measurements of either the cavity mode or the oscillator position, accompanied by Markovian feedback, are able to stabilize the nanosphere motion for all the values of the detuning ∆ and of the measurements efficiencies.
• The addition of time-continuous homodyne monitoring of the cavity output plus Markovian feedback greatly improves the performance that one would have obtained with sideband cooling only. For the realistic values of physical parameters considered in our study, while sideband cooling would prepare a phase-insensitive steadystate characterized by n ph 30 phonons on average, the addition of continuous homodyne measurements of the cavity output would prepare a squeezed steady-state, with n ph < 10 phonons. In particular, this is true for a large range of detuning values, which relaxes the requirements of sideband resolution to cool down the oscillator. This is particularly important for low frequency heavy oscillators which cannot be operated in the sideband resolved regime for dispersive coupling or where a dissipative coupling is not available [54]. • In terms of optimizing the purity and the cooling of the oscillator, sideband cooling would cease to be useful if one is able to directly measure the oscillator with very high efficiency. Nevertheless, if we take into account the state-of-the-art values for these measurement efficiencies in an actual experimental setup, the combination of the two control procedures represents still the best choice for experimental realizations, leading in principle to a quantum squeezed steady-state with less than one phonon on average. The performances are further improved if a simultaneous measurement of the cavity mode output is carried out.
• If we consider a state-of-the-art experimental setup where the dielectric nanosphere is trapped by the field of a high-finesse cavity, the proposed measurement protocols are in principle able to prepare a quantum squeezed state with less than n ph = 1 phonons at steady-state.
Over the next few years, it will be possible to perform more and more exhaustive time-continuous measurements on the outputs of interesting micro-and meso-scopic physical systems. It is apparent from our findings that such a possibility will be one of the pathways to reduce the entropic content of such systems, drive them to the quantum regime, and ultimately achieve their full or partial quantum control.

VI. ACKNOWLEDGEMENTS
MGG, JM, PFB and AS acknowledge support from EPSRC through grant EP/K026267/1. In this appendix, we provide the explicit form of the matrices entering the unconditional and conditional evolution equations for first and second moments, corresponding to the Eqs.
As we stated in the manuscript, the master equation (3)) describing the unconditional noisy evolution of the quantum state of the mechanical mode and of the cavity mode, yields the following evolution for the corresponding first moment vector R and for the covariance matrix σ: The drift matrix A and the diffusion matrix D read On the other hand, if we consider simultaneous continuous monitoring of the output cavity mode with homodyne detection, and of the oscillator position, the corresponding stochastic master equation (11) is translated for first moment vector and covariance matrix in: where dw = (dw 1 , dw 2 ) T , A has been defined just above and the other matrices read