Active-feedback quantum control of an integrated low-frequency mechanical resonator

Preparing a massive mechanical resonator in a state with quantum limited motional energy provides a promising platform for studying fundamental physics with macroscopic systems and allows to realize a variety of applications, including precise sensing. While several demonstrations of such ground-state cooled systems have been achieved, in particular in sideband-resolved cavity optomechanics, for many systems overcoming the heating from the thermal bath remains a major challenge. In contrast, optomechanical systems in the sideband-unresolved limit are much easier to realize due to the relaxed requirements on their optical properties, and the possibility to use a feedback control schemes to reduce the motional energy. The achievable thermal occupation is ultimately limited by the correlation between the measurement precision and the back-action from the measurement. Here, we demonstrate measurement-based feedback cooling on a fully integrated optomechanical device fabricated using a pick-and-place method, operating in the deep sideband-unresolved limit. With the large optomechanical interaction and a low thermal decoherence rate, we achieve a minimal average phonon occupation of 0.76 when pre-cooled with liquid helium and 3.5 with liquid nitrogen. Significant sideband asymmetry for both bath temperatures verifies the quantum character of the mechanical motion. Our method and device are ideally suited for sensing applications directly operating at the quantum limit, greatly simplifying the operation of an optomechanical system in this regime.


INTRODUCTION
Observing and utilizing quantum effects of a macroscopic mechanical resonator is of significant interest in physics.It offers great opportunities in understanding fundamental physics, such as quantum mechanics with massive objects [1,2], and in quantum applications with mechanical resonators, including quantum metrology [3] and tasks in quantum communications [4,5].However, the mechanical motion is typically overwhelmed by excess classical noise due to its contact with its surrounding thermal bath [6,7], making it challenging to use in quantum applications.Preparing and initializing the mechanical resonator close to its motional ground state reduces the classical noise, thus is a key pre-requisite for observing quantum behavior of the mechanical system.Several seminal demonstrations of reducing the phonon occupancy below 1 have been achieved over the past years, often in the sideband-resolved limit where the cavity linewidth is smaller than the mechanical frequency [8][9][10][11][12][13][14][15][16].These achievements have paved the way for the experimental observation of mechanical quantum behavior [17][18][19][20].In contrast, operating the optomechancial system in the sideband-unresolved regime has received increasing attention recently, in particular for sensing applications [2,[21][22][23][24][25][26][27][28][29].In such a sensor, the force to be measured is coupled to the displacement of the mechanical resonator, which is then read-out by the optical field.* s.groeblacher@tudelft.nl A large bandwidth of the optical cavity allows obtaining the displacement information faithfully, without any significant filtering by the cavity.Furthermore, in applications where a low frequency mechanical resonator is required or unavoidable due to its large size [7,22,26,[29][30][31][32][33][34][35], working in the sideband-unresolved regime significantly relaxes the stringent requirement on the optical cavity.For most experiments, low-mass, high frequency mechanical resonators [8,15] or Millikelvin bath temperatures achieved with dilution refrigerators [9,12] are typically used.The low temperature and high mechanical frequency reduces the initial phonon number n th = k B T /(ℏΩ M ) [21], where k B is the Boltzmann constant, ℏ is the reduced Planck constant, T is the bath temperature, and Ω M is the mechanical resonance (angular) frequency.In contrast, performing cooling on a more macroscopic system where the resonance frequency is lower, or at a higher temperature, is a more demanding task due to the larger initial n th .
Several experimental demonstrations with a variety of trapped or bulk systems have recently reached phonon numbers n < 1 [13,16,36,37] starting from large bath temperatures and with low frequency resonators.Here, we demonstrate feedback cooling [2,13,16,[38][39][40][41] of a fully integrated optomechanical resonator with a mechanical resonance frequency of only 1 MHz, an effective mass of 16 pg and an in-plane dimension of 0.5 mm, which is part of a fully integrated optomechanical system and only moderately pre-cooled in a continuousflow cryostat.Such a system is ideally suited for sensing applications due to its compact size and easy-to-use experimental setup.We measure the displacement of the mechanical resonator using a homodyne measurement, which is processed and sent to a controller that reduces the mechanical motional energy through active feedback control.By operating at sufficiently large feedback gain we observe sideband-asymmetry using an outof-loop heterodyne measurement.This imbalance in scattering rates between the Stokes and anti-Stokes processes is a hallmark of optomechanical interaction in the quantum regime [41][42][43][44][45].The asymmetry allows us to independently calibrate the absolute energy measurement of the mechanical mode [16,36,43], which is in good agreement with the inferred phonon number obtained from the calibrated homodyne measurement.This then allows us to perform feedback cooling to reach a minimum phonon number of 0.76 ± 0.16 starting from a mechanical mode temperature of 18 K, corresponding to an initial phonon number of 3.6 × 10 5 .We further demonstrate cooling from 77 K, where liquid nitrogen is used in the cryostat, reaching a minimum phonon number of 3.45 ± 0.15 and sideband-asymmetry is also observed.

Integrated high-Q, high-g0 optomechanical device
Our integrated device consists of a "soft-clamped" mechanical resonator inspired by a fractal structure [34] and an optical cavity formed by a photonic crystal.They are assembled using a pick-and-place method.Details of the device structure and the assembly method are discussed in [46].The mechanical resonator is fabricated from a 50nm-thick high-stress silicon nitride layer.The mechanical structure and the simulated mechanical mode are shown in Figure 1(a).The structure has a high aspect ratio, with dimensions on the 100's µm scale in the in-plane direction.Its fundamental mode oscillates at 1.1 MHz with an effective mass of 16 pg.The mechanical motion couples to the evanescent field of the photonic crystal cavity, which is made from a separate silicon nitride layer, and is placed above the center of the mechanical structure.The designed optical resonance wavelength is around 1550 nm.These two structures are separated by a small gap of 130 nm.We measure a mechanical quality factor Q M = 1.8 × 10 7 at room temperature, which increases to around 5.1 × 10 7 at 18 K (see Figure 1(b)).Due to the small gap, a large optomechanical coupling g 0 can be achieved, which is crucial for efficient feedback cooling.As shown in Figure 1(c, d), for the assembled device, we measure the reflectivity and the mechanical frequency shift at different laser frequency.The optical cavity has a linewdith (FWHM, full width at half maximum) of κ/(2π) = 8.8 GHz, including an external coupling of κ e /(2π) = 6.9 GHz, putting the device firmly in the sideband-unresolved cavity limit κ ≫ Ω M .The overcoupled optical cavity is required for a high total detection efficiency η det .The change of the measured mechanical frequency at different detuning arises due to the optical spring effect [21], from which we extract g 0 /2π = 224 ± 4 kHz.With the device, we also achieve a large single-photon cooperativity, C 0 = 4g 2 0 /(κΓ M ) ≈ 10 3 at 18 K, manifesting the high single photon interaction rate compared to the energy damping of the system [40].
As our device is in the sideband-unresolved regime, a measurement-based feedback protocol [39] can be used to reduce the energy in its mechanical motion.The displacement of the mechanical mode is continuously measured and sent to a controller.The controller then processes the information in real-time, and it applies a feedback force onto the mechanical resonator.In optomechanics, this force can be introduced by modulating the optical input power [13,47].For an optical cavity with large bandwidth, the modulation of the input power is converted to a change of the radiation force [21], without significant filtering effect in the frequency range of interest.In this scheme, realizing a fast and precise measurement is crucial.In a quantum-noise limited measurement with coherent light as the input, a measurement rate is defined to characterize the time-scale at which the zero-point fluctuation can be distinguished from shot noise, Γ meas = 4η det g 2 /κ, where g = √ n c g 0 is the multi-photon optomechanical coupling rate with a cavity photon number n c [40,47].In order to achieve ground state cooling, the measurement rate should be comparable to the decoherence rate, Γ dec = n tot Γ M [40].
Here, n tot = n th + n ba is the total effective bath phonon number, including the excess phonon number due to the FIG. 2. Feedback cooling at 18 K.(a) Experimental setup.EOM: electro-optic modulator for phase (blue) and intensity (yellow) modulation.φ: fiber stretcher for phase locking.In the full setup two sets of measurements are performed simultaneously.The orange lines show the optical paths for the feedback cooling, where the laser is on resonance with the optical cavity, and is modulated by a phase EOM to generate a phase calibration tone.A homodyne scheme is used to measure the phase quadrature of the reflected light.The measured signal is sent to a controller, which controls the light intensity through the intensity EOM.A heterodyne detection scheme (red laser) is used to measure the sideband-asymmetry.A phase EOM, driven by a sawtooth signal, is used to shift the optical frequency by Ω het .(b) Power spectral density (PSD) of the homodyne measurement, converted to the mechanical displacement quanta (XM), at various feedback gains (g fb ).The gray dashed line shows the shot noise level, and other dashed curves are the fits to the data.(c) Power spectral density of the heterodyne signal, normalized to shot-noise, showing two asymmetric sidebands.Integrating the area gives the corresponding phonon number, shown in (d).The left panel is obtained by fitting the spectrum and integrating the fit, while the data on the right results from directly integrating the measured spectra.Triangles show the sideband power, crosses show the magnitude of the asymmetry, and the dashed line are for the expected asymmetry.The reduction of the asymmetry on the right panel at high feedback gain is due to the finite integration range, which does not fully capture the broadened mechanical spectrum.(e) Phonon number obtained by homodyne measurement (blue), integrating the fit of the sideband-asymmetry spectrum (orange), and the direct integration of the sideband-asymmetry spectrum (green).More details on the data processing can be found in the Methods section.Error bars in all panels represent standard deviations.
back-action noise n ba = n c C 0 [40].It originates from the quantum fluctuation of the cavity field, which disturbs the mechanical motion through optomechanical coupling.
Comparing the two rates, This ratio should approach 1 for ground state cooling.The first term compares the thermal decoherence rate to the measurement rate and shows that a high C 0 is beneficial, as provided by our device.It reduces the need for a large cavity photon number, which would generate heat and hence raise the bath temperature due to photon absorption [48][49][50].Furthermore, as indicated by the second term due to the back-action noise, there is an ultimate lower bound for the ratio.Any loss in the detection process reduces the precision of the measurement as less information is obtained and therefore a high detection efficiency η det is critical.

Feedback cooling and sideband asymmetry with LHe
The experimental setup is shown in Figure 2(a).The device is placed in a continuous-flow cryostat, which can be cooled either by liquid helium or liquid nitrogen.A thermometer underneath the sample stage measures the cold-fingers temperature, which reaches 6 K when using liquid helium.The mechanical frequency slightly reduces to 1.045 MHz at low temperature, which corresponds to a thermal phonon occupation of around 10 5 .A laser, on resonance with the optomechanical cavity, is used for the feedback cooling.A phase EOM, directly after the laser and driven by a sinusoidal signal, generates a phase modulation tone of known amplitude.Both the phase modulation signal and the mechanical motion experience the same transduction due to the optical cavity [51].This allows us to calibrate the mechanical displacement using the independently obtained g 0 .A balanced homodyne measurement [52] is set up to measure the phase quadrature of the light, which contains displacement information of the mechanical resonator [21].The measured signal is then fed to a controller (RedPitaya 125-14), whose output is connected to an intensity EOM and modulating the input light intensity.The controller can realize a complicated transfer function, which is required for a system with multiple mechanical modes and with significant delay [47].In this scheme, parasitic detection of the intensity modulation occurs.We have verified that the effect is small, introducing an error of only ∼10 −6 to the feedback force (see Methods section).Another laser, red-detuned from cavity resonance, is used to perform an out-of-loop heterodyne measurement [16,41].The light is scattered by the mechanical resonator into two sidebands with lower and higher frequency, whose magnitude ratio is given by n/(n + 1), where n is the average phonon number of the mechanical mode [9,16,41,45].The extra quanta in the ratio, which is due to the quantum fluctuation of the mechanical motion [41,43,44], provides an absolute energy scale in our system and can be obtained by measuring the spectrum of the heterodyne signal [16,36].In our setup, the heterodyne probing light, with both sidebands, are selected by a filter cavity (MicronOptics FFP-TF2, κ/2π ≈ 35 MHz).It is then combined with the heterodyne local oscillator whose frequency is shifted by Ω het /(2π) = 2.81 MHz with a serrodyne scheme [19,53].Since our system is in a regime where Ω M /κ ∼ 10 −4 ≪ 1, correction of the two sidebands due to the optical cavity [9,16] is not required.Classical amplitude noise might generate spurious sideband asymmetry in the heterodyne detection [41].In our setup, the heterodyne probe beam merges with the cooling beam only after the intensity EOM, thus the intensity EOM does not introduce additional classical noise and information to the heterodyne probe beam, and the heterodyne detection is unaffected by our cooling modulation.Furthermore, its detuning with respect to the cooling laser (≳1 GHz) is much larger than the detection bandwidth of our photodetector, which eliminates any interference from our feedback control onto the heterodyne probe.
The spectra of the mechanical displacement, obtained from calibrated homodyne measurements, is shown in Figure 2(b).We keep a fixed input power and change the feedback strength by changing the electrical gain of the feedback controller.The feedback introduces an extra damping channel to the mechanical resonator.When the feedback strength is increased the extra damping becomes stronger removing energy from the mechanical motion, manifesting in a broadened mechanical peak with reduced amplitude.At large gain, it enters a noise squashing regime [13] where the measured mechanical signal cancels with the measurement noise.We fit the mechanical spectrum to obtain the average phonon occupancy, as shown in Figure 2(e).The phonon number initially reduces at small gain, but eventually increases again due to the noise being fed into the mechanical resonator.By fitting the spectra we obtain a minimum phonon occupancy of 1.06 ± 0.06.The heterodyne measurement is shown in Figure 2(c), where the power spectral density is normalized to the shot noise level.We clearly see a difference in the amplitude of the two sidebands.We fit the two sidebands to extract the corresponding energy.The result is plotted in the left panel in Figure 2(d), where all the values are normalized by the average energy difference of all curves.The energy difference remains stable over different feedback strength.At high gain, the mechanical peak is broad and its amplitude is small, which makes fitting no longer possible.We also directly integrate the measured heterodyne spectrum over a frequency range between 1.035 MHz and 1.05 MHz, around the mechanical central frequency.We then subtract the noise floor obtained from the fit.The result is shown in the right panel of 2(d).Since the mechanical peak broadens at larger gain, the energy inside the integration range reduces and thus the energy difference of the two sideband reduces at larger gain.It nevertheless agrees with the theoretical expectation from the broadening of the peak.We compare the phonon number calibrated from the sideband asymmetry to the phonon number obtained by fitting the homodyne measurement in Figure 2(e), and find both methods to be in good agreement.

Ground-state cooling with improved detection efficiency
As indicated by Eq. ( 1), achieving a low Γ dec /Γ meas requires a high detection efficiency η det .Performing two sets of measurements simultaneously requires additional optical components to separate the measurement results, which introduces losses.Furthermore, an out-of-loop heterodyne measurement requires sending additional light into the optical cavity, which leads to larger quantum back-action noise.This is equivalent to a lower η det as the light for the heterodyne measurement does not contribute to the feedback cooling.Still, the quantum fluctuation of the light perturbs the mechanical resonator, leading to an increase of the motional energy.As we have confirmed the agreement for the extracted phonon number between the homodyne measurement and the sideband asymmetry measurement, which shows the validity of our calibration of the in-loop homodyne measurement, A single, slightly red-detuned laser is used to perform the measurement.(b) A measured spectrum, converted into the displacement quadrature quanta, of the mechanical resonator, at different feedback strength.The fits (dashed lines) are used to extract the system parameters and the phonon occupancy.The extracted effective bath temperature (c) is higher than the expected value, indicating excess decoherence (e.g., higher bath temperature).(d) Shows the extracted phonon number, where the gray region shows an occupancy below 1. Error bars in all panels represent standard deviations.
we performed a feedback cooling measurement without the sideband asymmetry measurement setup.Extra optical components are removed, as shown in Figure 3(a), to increase the detection efficiency.We fix the cavity photon number to n c ∼ 350.With a detection efficiency of η det ≈ 0.49, we have Γ dec /Γ meas ≈ 2.7.We measure and fit the spectrum for different feedback strengths (see Figure 3(b)), and extract the phonon number, shown in Figure 3(d).In these measurements a minimum phonon number of 0.76±0.16 is achieved.We note that the fitted effective bath temperature, with the contribution from the back-action noise n ba subtracted, corresponds to a thermal temperature of 18.3 K (cf. Figure 3(c)), which is significantly higher than the readout from the thermometer.Such a discrepancy is indeed quite common due to heating of the mechanical structure and non-ideal thermalization of our chip [8].For comparison, a temperature of 6 K would allow a minimum phonon number of 0.45.

Feedback cooling with LN2
Cooling a mechanical resonator to close to its motional ground state and revealing quantum effects at a higher and more accessible temperature is technologically of utmost importance but also significantly more challenging.
Here we demonstrate feedback cooling starting at 77 K, where the bath is pre-cooled by flowing liquid nitrogen.We otherwise use the same setup as shown in Figure 2(a).
The spectrum in the homodyne measurement is shown in Figure 4(a).At higher temperature, the thermal decoherence rate is higher and thus a stronger measurement is needed to achieve a low phonon occupancy.However, still within the regime Γ meas ≪ Γ dec , the measurement noise floor increases.This is a direct result of the mechanical motion of the low Q M mode of the photonic crystal [46], which shows a strong 1/f feature [6,55] and affects the measurement of the low frequency mechanical mode.In addition, the initial mechanical quality factor is lower at this elevated temperature with Q M = 4.1 × 10 7 , limiting the minimum achievable phonon number.Nevertheless, sideband-asymmetry (see Figure 4(b)), as a unique feature of quantum physics [41][42][43][44]56], is still observed at this higher temperature.We extract and compare the phonon occupancy obtained by fitting the calibrated homodyne spectrum, the heterodyne spectrum and by directly integrating the area of the heterodyne spectrum.The results are presented in Figure 4(c), showing consistency among the methods.The minimal phonon number of 3.45 ± 0.15 is determined from the fitting of the homodyne measurement.

DISCUSSION
We have fabricated an integrated optomechanical device, achieving large optomechanical coupling (g 0 /(2π) = 224 kHz) and a high mechanical quality factor (Q M = 5.1 × 10 7 ) for the fundamental out-of-plane mechanical mode at 1 MHz.We pre-cool the bath with liquid helium to an effective mode temperature of around 18 K, where the corresponding thermal phonon occupation is around 3.6 × 10 5 , and perform measurementbased feedback cooling to reduce the motional energy.Using sideband-asymmetry, we verify our measurement scheme, which agrees well with the result obtained from the calibrated in-loop homodyne measurement, confirming the validity of our calibration.As this double measurement scheme reduces the detection efficiency we obtain an occupation of 1.06 ± 0.06, close to the ground state, compared to 0.76 ± 0.16, when the setup's efficiency is improved.We further demonstrate that even starting from 77 K, which is only a factor of 4 from room temperature, our technique allows us to reach a regime with significant asymmetry in the two motional sidebands, underlining the quantum nature of the created state.Small improvements in the mechanical quality factor will allow to directly reach the ground state even from ambient temperatures.The use of our fully integrated optomechanical device, deep in the sidebandunresolved regime, and our simple experimental setup with the ability to cool its fundamental mode into its motional ground state is highly relevant for real-world applications and will allow for using this technology as an easy-to-use quantum technology.This could include,   4. Feedback cooling from an initial bath temperature of 77 K, pre-cooled with liquid nitrogen.(a) Power spectral density of the mechanical displacement for various feedback gains.(b) Heterodyne spectrum signal, normalized to the shot noise, for various feedback gains.The asymmetry in the two sidebands at large gain is clearly visible, where the Stokes scattering (red trace) has a higher power spectral density.Extra peaks are present due to the mixing higher order modes in the detection [54].(c) Average phonon number extracted from the homodyne (blue) and from the heterodyne measurement (orange: integrating the fit of the heterodyne spectrum; green: integrating the heterodyne spectrum directly).

Experimental setup
We use a Santec TSL-510 laser as our cooling laser, and a New Focus TLB-6728 as the heterodyne probe laser.In both homodyne and heterodyne setups, the path lengths of the local oscillator and the signal arm are matched.We confirm this by scanning the laser over a broad wavelength range and observing the interference fringes.An additional free-space path is added to the homodyne local oscillator to fine-tune the optical path length.
In our experiment, a phase modulation tone is generated to calibrate the mechanical displacement.[51] The phase modulation δϕ PM = ϕ 0 cos Ω PM t is equivalent to a frequency modulation, δω L = dϕPM dt = −ϕ 0 Ω PM sin Ω PM t.Here, ϕ 0 is the phase modulation depth, which is much smaller than 2π.Note that the optical cavity transduces both the laser frequency modulation and the mechanical displacement to the phase quadrature.By comparing the homodyne signal, the laser frequency modulation due to the phase EOM is compared to the cavity frequency modulation due to the mechanical resonator, and the mechanical displacement is calibrated.Due to the large optical linewidth of the optomechanical cavity, the resulting transduction to the phase quadrature is very small around 1 MHz.It is then susceptible to the residual amplitude modulation from the phase EOM and the residual amplitude detection in the homodyne measurement scheme.We therefore increase the frequency of the phase modulation tone to 60 MHz.While it is still much smaller than the cavity linewdith, the corresponding frequency modulation is larger, and therefore the transduction is larger.It allows us to unambiguously measure the phase modulation.In the detection, the frequency responses of electronics are not flat across the large frequency range.We normalize all the measured values to the optical shot noise, which is measured by blocking the signal arm.Since the shot noise is flat in spectrum, the measured spectra in different frequencies are calibrated and comparable.

Data processing
In the experiment, we measure the spectrum of the photodetector voltage.In order to account for any filtering effect from the photodetector and the electronics, we normalize all the spectra to the spectrum corresponding to the shot noise.The shot noise spectrum is measured by blocking the light entering the optomechanical cavity and leaving the other experimental parts unchanged.
For the cooling and the homodyne measurement, it is convenient to model our system in a semi-classical way [13,50].The displacement of the mechanical resonator in the frequency domain and normalized to the zero-point motion is [68] where the susceptibility F (ω) = F N (ω) + F fb (ω) is the input "force", consisting of input noise noise and the feedback control.The input noise has multiple origins.The thermal noise, in the high temperature regime (k B T /(ℏΩ M ) ∼ 10 5 ≫ 1 in our experiment), the spectrum is flat [68] S F th (ω) = 2Γ M (2n th + 1).
It is therefore not possible to distinguish F ba and F th directly in the experiment.We introduce an effective bath temperature to include the effect of the backaction noise, T eff = T th + T ba , and where n bath = k B T eff /(ℏΩ M ).Since our system is in the sideband-unresolved limit, the optical cavity does not have a filtering effect in the control of the mechanical resonator.The feedback force is directly proportional to the output of the electronic controller with a transfer function H fb , Here, H fb (ω) = Hfb (ω)e iωτ fb includes signal transmission and processing delay τ fb , and Hfb (ω) is the transfer function without time delay.For a sideband-unresolved system, The readout of the actual displacement of the mechanical resonator is not filtered by the optical cavity, giving the first term.The second term is the imprecision of the measurement, including the quantum fluctuation of the measured optical field (vacuum noise) [13,68] and any possible classical measurement noise.We also consider the third term, which is due to the imperfection of our experiment.In our experiment, the feedback controller modulates the intensity of the input light, and we measure the phase quadrature of the light.The modulation and the readout belong to the same optical mode, but they are on different quadratures.However, experimentally, it is challenging to keep the modulation and the measurement perfectly orthogonal.A small error results in measuring a small amount (quantified by a free parameter ϵ fb ) of the modulation signal.Experimentally, before performing the cooling, we generate an amplitude modulation tone from the feedback controller, and we finetune the homodyne phase locking such that the detected power of the modulation tone is minimized.This guarantees an "almost" perfect detection, where the detected quadrature is orthogonal to the modulation quadrature.In the fit (see below) we obtain ϵ fb ∼ 10 −6 .
Experimentally, the power spectral density (normalized to the shot noise level) is measured.Combining the above equations, By using S FN , S X imp M , g fb , ϵ fb , and τ fb as fitting parameters to fit the experimental curve, we can extract their values.The spectrum of the actual displacement can be inferred (9) The effective average phonon number is the integral of the inferred mechanical spectrum [39] n = 1 2 which is evaluated numerically.Our heterodyne measurement is performed by shifting the frequency of the local oscillator beam by Ω het .It is then possible to access the two sidebands resulted from the optomechanical interaction [16,36,41,42].For our out-of-loop heterodyne detection, we fit the spectrum of the two sidebands by where l, r are for the left and right sidebands.k j represents the magnitude of the sideband, and n j is for the noise floor at the vicinity of the sideband.The spectral frequency of the two sidebands are shifted due to the frequency difference on the local oscillator, ω l = Ω Het − ω, and ω r = ω − Ω Het and are fitting parameters.The susceptibility is approximated to have the same form as the susceptibility of a bare mechanical resonator, where the effective resonance frequency ΩM and the effective mechanical dissipation rate ΓM is introduced due to the effect of the feedback cooling, and they are left as fitting parameters.From the fitting, it is then possible to obtain the phonon occupation number n = 1 2 since the energy difference between the two sidebands corresponds to 1 phonon, and the sum corresponds to 2n+1 phonons [16,36,41,42].We further verify that the power difference between the two sidebands is a constant by integrating Equation (11) (see Figure 2(d, left)).Alternatively, it is possible to integrate the measured and normalized heterodyne spectrum directly.We perform the integration by summing all the power spectrum near the mechanical peaks, with the noise floor n j subtracted.Let the summed results for the two sidebands be s l and s r , the phonon number is then given by n = 1 2 Due to the finite integration range and the broadening of the mechanical spectrum, the power difference between the two sidebands is no longer a constant.We compare it to the theoretical value, obtained by integrating the inferred mechanical spectrum from the homodyne measurement and within the corresponding frequency range.They show good agreement (cf. Figure 2(d, right)).

Feedback filter design
In our device, the targeted high-Q M mode is the fundamental mode.Other mechanical modes are far away, in particular compared to other mechanical structures based on phononic crystals [13,47].However, high-order out-of-plane modes still exhibit relatively large optomechanical coupling in our structure.In the presence of the signal transmission and processing delay, which is comparable to the oscillation period of the the mechanical modes, an extra phase lag is therefore introduced.We tune the filter to have a phase response such that other mechanical modes are not heated, as described in detail in [47].At different temperature, the mechanical frequencies shift slightly, and we tune the filter separately.The total feedback delay is tuned to 640 ns (with LHe cooling) and 680 ns (with LN 2 cooling).

Laser noise for sideband-asymmetry detection
The phase and amplitude noise of the heterodyne detection laser affects the observed sideband asymmetry [41].We verify that they do not have significant effect on the heterodyne detection.
Classical laser amplitude noise might enhance the asymmetry [41].We measure the laser amplitude noise by sending light to a variable coupler (beamsplitter) and then a balanced photodetector, as shown in Figure S1(a).We then take the spectrum of the received signal.When the splitting ratio of the variable coupler is set to 50:50, the classical noise is canceled, and the spectrum of the measured signal shows the amplitude shot noise, and additional electronic noise.When the variable coupler is set to 100:0, the measured spectrum consists of classical laser noise, shot noise, and electronic noise.The electronic noise is measured by blocking the laser.After subtracting the electronic noise, we can estimate the ratio of the classical amplitude noise to the shot noise.For our typical measurement power (∼ 1.5 µW), excess classical noise, normalized to shot noise, is about 2 × 10 −3 .
The measured spectrum and the classical amplitude noise are shown in Figure S1(b,c).
The impact of the classical phase noise to the laser, whose PSD is C Y when normalized to the shot noise, is approximately 4(∆ω/κ 2 )C Y [41].The additional noise is added to both sidebands [41], reducing the sideband asymmetery.The laser phase noise is measured by sending the laser through a filter cavity.The laser is on resonance with the filter cavity, and we do a phase measurement through a homodyne setup (Figure S1(d)).The measured spectrum is compared to a phase calibration tone (1.05 MHz) generated by a phase EOM at the common path of the homodyne setup, and it allows us to calibrate the phase noise PSD C θ .We then calculate the normalized phase quadrature spectrum by using the input power to the filter cavity [41], C Y = 2nC θ , where n is the input photon flux.The classical phase noise around the mechanical frequency C Y ≈ 430, and the correction due to the classical laser phase noise is 4(∆ω/κ 2 )C Y ≈ 0.05.The measured spectrum and the calculated classical phase noise are shown in Figure S1(e,f).

FIG. 1 .
FIG. 1. Device design and characterization.(a) Mechanical structure and the simulated displacement field of the fundamental mode.At the center (blue box and inset), a photonic crystal cavity (purple) is placed above the mechanics.The spacers (yellow) define the distance between the photonic crystal and the mechanical structure.(b) Ringdown measurement of the fundamental mechanical mode at 18 K.The fit (orange dashed line) allows us to extract a QM = 5.1 × 10 7 .(c) Cavity reflection and (d) mechanical frequency shift due to the optical spring effect at different laser detuning and at room temperature.A fit to a simple model (orange dashed curves) is used to extract g0.

FIG. 3 .
FIG. 3.Feedback cooling without sidebandasymmetry measurement.(a) Feedback cooling setup.A single, slightly red-detuned laser is used to perform the measurement.(b) A measured spectrum, converted into the displacement quadrature quanta, of the mechanical resonator, at different feedback strength.The fits (dashed lines) are used to extract the system parameters and the phonon occupancy.The extracted effective bath temperature (c) is higher than the expected value, indicating excess decoherence (e.g., higher bath temperature).(d) Shows the extracted phonon number, where the gray region shows an occupancy below 1. Error bars in all panels represent standard deviations.

TABLE I .
Experimental parameters.