Quantum nondemolition measurement of mechanical motion quanta

The fields of optomechanics and electromechanics have facilitated numerous advances in the areas of precision measurement and sensing, ultimately driving the studies of mechanical systems into the quantum regime. To date, however, the quantization of the mechanical motion and the associated quantum jumps between phonon states remains elusive. For optomechanical systems, the coupling to the environment was shown to make the detection of the mechanical mode occupation difficult, typically requiring the single-photon strong-coupling regime. Here, we propose and analyse an electromechanical setup, which allows us to overcome this limitation and resolve the energy levels of a mechanical oscillator. We found that the heating of the membrane, caused by the interaction with the environment and unwanted couplings, can be suppressed for carefully designed electromechanical systems. The results suggest that phonon number measurement is within reach for modern electromechanical setups.

E nergy quantization is one of the hallmarks of quantum mechanics. First theorized for light by Einstein and Planck, it was found to be ubiquitous in nature and represents a cornerstone of modern physics. It has been observed in various microscopic systems starting from nuclei, atoms and molecules, to larger mesoscopic condensed matter systems such as superconductors 1 . For macroscopic systems, however, the observation of energy quantization is hindered by the smallness of the Planck constant. Thus, although being a milestone of contemporary physics, up to date the discrete energy spectrum of mechanical resonators has never been seen directly.
Extreme progress in studying mechanical systems has been achieved in experiments exploiting radiation pressure. This is the core of optomechanics 2 , where photons and phonons of the optical and mechanical subsystems interact with each other. A similar type of coupling can be realized in the microwave domain with electrical circuits, leading to the field of electromechanics [3][4][5][6][7][8] . The numerous advances of optomechanics and electromechanics include ground state cooling 4,5,[9][10][11] , ultra precise sensing [12][13][14][15] , generation of squeezed light and mechanical states 7,8,[16][17][18] , back action cancellation 19,20 and detection of gravitational waves 21 . In all of these systems, however, the operation in the single-photon/ phonon regime is challenging due to the small value of the bare coupling 3,22 . Instead, experiments exploit an enhanced linearized effective coupling induced by a large driving field. This severely limits the nature of the interactions 23 and possible quantum effects. In particular, it precludes the observation of the energy quantization in mechanical resonators.
Quantization of mechanical energy can be observed by a quantum nondemolition (QND) measurement 24,25 of an oscillator's phonon number operatorn b . Here, QND means that the interaction, which couples the mechanical system with the measurement apparatus, does not affect the observable we are interested in. This is achieved if the total Hamiltonian commutes withn b , and the influence of the environment is minimized.
Considering the electromechanical setups in Fig. 1, we show that QND detection is feasible for a capacitor in which one of the electrodes is a light micromechanical oscillator. By choosing an antisymmetric mode for the oscillator, the interaction between the electrical and mechanical subsystems is quadratic in the displacement. Along with the suppression of the linear coupling, this ensures the QND nature of the measurement, as originally proposed in refs. 26,27 for an optomechanical system. In that system, however, it was shown in refs. 28,29 that the combination of unwanted losses and the coupling to an orthogonal electromagnetic mode spoils the interaction, unless strong single-photon coupling is achieved. Here, we show that for the considered electromechanical setup the equivalent orthogonal mode can have dramatically different properties, allowing for the phonon QND detection. We derive general conditions under which the QND measurement is possible, and characterize its experimental signatures. As compared to most approaches to phonon QND measurement 26,27,[30][31][32] , our procedure does not impose stringent requirements on the single-photon optomechanical coupling, but relies on the ratio of the involved coupling constants. This makes our approach attractive even for systems where the interaction is limited, for example, due to stray capacitances in the setup. For a measurement of the square displacement, a similar advantage was identified in ref. 31 .

Results
Proceeding. We first study an RLC circuit with one capacitor plate being an oscillating membrane, without assuming the symmetry discussed above (Fig. 1b). The mechanical motion of the plate shifts the resonance frequency of the circuit, while the electric potential exerts a force on the membrane. In order to perform a QND measurement of the phonon number, we require this interaction to be proportional ton b . We therefore Taylor expand the inverse of the capacitance to second order in the displacement, 1=CðxÞ ≃ C À1 0 +g 1 ðb þb y Þ +g 2 ðb þb y Þ 2 =2, where we replaced the positionx with the creationb y and annihilationb operators of the mechanical motion, andg 1;2 denote linear and quadratic coupling constants. Within the rotating wave approximation,g 2 ðb þb y Þ 2 =2 ≃g 2nb , leading to the desired QND interaction, while theg 1 term adds unwanted heating that spoils the phonon measurement.
The main aim of this work is to identify conditions under which the QND measurement is feasible, despite the presence of heating. We initially consider the simple circuit in Fig. 1b, and assume the incoming signalV in to be in a coherent state resonant with the circuit. The quadratic interaction then shifts the electrical resonance frequency proportionally to the phonon numberg 2nb . For smallg 2 , this shift leads to a phase change of the outgoing signalV out that can be determined by homodyne measurement. Different phononic states will thus lead to distinct outcomes V M , as shown in Fig. 2. The distance d between output signals for differentn b and the standard deviation σ of the noise define the signal-to-noise ratio D = d/σ (see Fig. 2), which needs to be maximized.
In order to have a successful QND measurement, the phonon numbern b must be conserved. If the mechanical state jumps during a measurement, the outcome V M ends up between the desired peaks. This leads to a reduced contrast, as illustrated by the distribution in the background of Fig. 2. The probability forn b to change is generally state-dependent, in the sense that higher Fock states are more likely to jump. A state-independent characterization of this heating is given by the average phonons Δn b added to the ground state during the measurement time T. The jump probability for any state can then be derived from Δn b using standard results for harmonic oscillators (for details see Supplementary Note 3 available in Supplementary Material online).
Both D and Δn b are proportional to the incoming intensity. We therefore characterize a setup by the parameter λ = D 2 /Δn b , where λ ) 1 is required for successful QND detection. For the RLC circuit in Fig. 1b, we find below that where g 1 ¼g 1 C 0 ω s , g 2 ¼g 2 C 0 ω s and n e is the thermal occupation of R 0 and Z out (assumed equal, R 0 = Z out ). Here, ω m and ω s = (C 0 L 0 ) −1/2 ≫ ω m are the mechanical and electrical frequencies, respectively, and γ t = Z out /L 0 corresponds to the output coupling rate. A result similar to Eq. (1) is derived in ref. 33 .
Despite progresses in reaching the resolved sideband regime ω m ) γ t in both optomechanical and electromechanical systems, g 2 is generally much smaller than g 1 , implying λ ( 1 in Eq. (1). To circumvent this problem, we use the second fundamental mode of the membrane in the capacitor, as depicted in Fig. 1a. The first-order coefficientg 1 of the 1=CðxÞ expansion then vanishes, leavingg 2 to be the largest contribution to the electromechanical coupling. In this situation λ seemingly grows indefinitely, the induced heating disappears and the QND measurement of the phonon number is easily realized. In practice, however, two effects will limit the achievable value of λ. First, inaccuracies in the nanofabrication can cause misalignments and, consequently, a residual linear coupling. Second, the oscillation of the membrane induces a charge redistribution in the capacitor to maintain it at an equipotential. The associated antisymmetric electrical mode introduces an effective linear coupling, and a similar heating mechanism as the one identified in ref. 28 for the optomechanical setup of refs 26,27 . In these papers, the quadratic interaction results from a hybridization of two modes linearly coupled to the mechanical position, and the QND detection was found to be impossible unless the single-photon coupling g 1 exceeded the intrinsic cavity damping. In our case, the QND interaction arises directly from the Taylor expansion of the capacitance. Hence, there is no constraint tying the second-order coupling g 2 to the properties of the symmetric and antisymmetric electrical modes, which can have vastly different resonance frequencies and dampings. This inhibits the mechanical heating and ultimately allows for the QND detection of the phonon number. We model the charge redistribution in the capacitor by parasitic inductances (L) and resistances (R) in the equivalent circuit of Fig. 1c. Each of the two arms containing R and L represents one half of the capacitor, with opposite dependence on the membrane position, CðxÞ and CðÀxÞ.
Single-arm RLC circuit. In the following, we derive Eq. (1) for the RLC circuit in Fig. 1b. The methods sketched here will then be generalised for the double-arm circuit in Fig. 1c. Using the standard approach 34 , we write the circuit Hamiltonian asĤðxÞ =Φ 2 = 2L 0 ½ þQ 2 = 2CðxÞ ½ , where the conjugate variableŝ Q andΦ are the charge and magnetic flux, respectively. We can expandĤðxÞ in the mechanical positionx /b þb y , in order to obtain the circuit HamiltonianĤ e =Ĥx ¼ 0 ð Þ and the coupling HamiltonianĤ em = g 1 ω s L 0Q 2 ðb þb y Þ=2 + g 2 ω s L 0Q 2 ðn b þbb=2 þb yby =2Þ. The total HamiltonianĤ tot = H e þĤ em þĤ m is therefore the sum of the circuit, interaction and the mechanical HamiltonianĤ m = hω mb yb .
Next, we describe the environmental effects corresponding to decay and heating of the modes. Associating each resistor R i with its own Johnson-Nyquist noiseV R i , we find the equations of motion of the composite system : : where γ r = R 0 /L 0 , γ b is the intrinsic mechanical damping rate with associated noiseF b and is the amplitude of the zero-point motion for a membrane of mass m. From now on, we consider optimally loaded setups with γ r = γ t . Equations (2)-(4) fully characterize the dynamics of the system, and represent the starting point for our detailed analysis.
a Sketch of a capacitor with an oscillating plate, here represented by a graphene membrane. We consider an antisymmetric (2, 1) mechanical mode. b RLC oscillator formed by the inductance L 0 , resistance R 0 , and position-dependent capacitance CðxÞ. The circuit is driven by the input voltageV in through a transmission line of impedance Z out .V out is the reflected signal. c Model circuit for an RLC system where the capacitor has the same form as in a. The membrane has a vanishing linear coupling to the symmetric electrical mode used for probing the system. The antisymmetric mode, residing in the small loop containing parasitic inductances L and resistances R, describes the redistribution of charge on the capacitor The feedback of the membrane's motion on the electrical circuit is described by Eq. (3). Driving the system at the electrical resonance frequency ω s , the terms proportional to g 1 ðb þb y Þ and g 2 ðbb þb yby Þ give rise to sidebands at frequencies ω s ± ω m and ω s ± 2ω m , respectively, whereas g 2nb induces a phonon-dependent frequency shift of the microwave cavity. Since homodyne detection is only sensitive to signals at the measured frequency, the sidebands are removed in the outcome V M , which is defined as the phase quadrature ofV out ¼V in À γ tΦ . This allows us to neglect oscillating terms in the calculation of V M (the linear term also leads to mechanically induced damping, but this is typically negligible compared to γ t ). The only contribution to V M is therefore the phonon-dependent frequency shift, which allows us to resolve the mechanical state. On the contrary, the electrically induced mechanical heating only involves the sidebands ω s ± ω m and ω s ± 2ω m , being unaffected by the term g 2nb in the Hamiltonian. For the RLC circuit in Fig. 1b, the heating is dominated by the linear term, since g 1 ) g 2 , and we shall neglect g 2 for the calculation of Δn b .
Below, we quantify the heating of the membrane and the phonon-dependent LC frequency shift. We first assume that the mechanical state does not jump during the measurement. Then, the equations of motion of the two subsystems decouple and we , where the number of photons α j j 2 sent into the circuit within the measurement time T sets the measurement strength. As discussed above, Δn b is the average phonon number at the end of the measurement Δn b =n b ðTÞ h i, with the mechanics initially in its ground state. For T much shorter than the mechanical lifetime γ À1 b , Δn b can be linearized to find the rate at which the membrane heats up. For the RLC circuit in Fig. 1b The parameter λ given in Eq. (1) is then found as the ratio λ = D 2 /Δn b . For details see Supplementary Note 1 available in Supplementary Material online.
Double-arm circuit. With the overall linear coupling vanishing, the parameter λ will be limited by fabrication imperfections and coupling to the antisymmetric mode. To model these phenomena, we consider the circuit in Fig. 1c, where the antisymmetric mode resides inside the small loop containing the two capacitors, and the symmetric one probes the system. We derive g 1 and g 2 from the expansion of each of the two capacitors: 1=C ±x ð Þ≃ C À1 0 ±g 1 ðb þb y Þ þg 2nb , so that in the absence of fabrication imperfections the total capacitor C tot = Cx ð Þ þ C Àx ð Þ is not linearly coupled to the symmetric mode. The coefficients g 1 and g 2 are related to their tilde counterparts in the same way as before, and the parameters D 2 and Δn b are evaluated in a similar fashion as we did for the RLC circuit. Since we quantify two sources of heating, it is convenient to write λ = ðλ À1 b þ λ À1 p Þ À1 , where λ b takes into account heating from charge redistribution, and λ p describes the influence of fabrication imperfections. With the details presented in Supplementary Note 2 (for details, see Supplementary Material available online) and Methods, we find where ω s = [C 0 (L + 2L 0 )] −1/2 is the frequency of the symmetric mode, γ t = [2Z out ]/[L + 2L 0 ] is the decay to the transmission line and g r = 2C 0 x 0 ω s ∂ x C À1 tot ðxÞ is the residual linear coupling induced by fabrication imperfection. We use the same notation introduced for the RLC circuit to allow a direct comparison. Equations (5) and (6) express the gain of our approach to QND detection. First, Eq. (6) quantifies the advantage of symmetry: λ dramatically improves compared to Eq. (1) by having a small residual linear coupling g r ( g 1 . Second, Eq. (5) is multiplied by the factor (ω s /ω m ) 2 with respect to Eq. (1). For microwave readout of a megahertz oscillator, this factor can be substantial. Furthermore, the mechnical oscillator is now only susceptible to the noise associated with charge redistribution on the capacitor, and not to the resistance in the inductor. This gives an additional improvement if R < Z out . Fig. 2 Sketch of the experimental outcome. Distribution of outcomes V M for two different phonon numbers: n b = 0 (first peak to the left) and n b = 1 (last peak to the right). For a given value of n b , repeated measurements are Gaussian distributed with a variance σ 2 / 1 þ 2 n e of the outgoing signalV out , consisting of vacuum and thermal noise. The distance d between the two peaks depends on the circuit parameters and the number of incident photons, and identifies the signal-to-noise ratio D = d/σ. Ideally, for each shot of the measurement, the mechanics is either in its ground or first excited state. However, for Δn b > 0 there will be events where the mechanical state jumps, resulting in outcomes V M in between the peaks relative to n b = 0 and n b = 1 (smaller peaks in the figure). This leads to the smeared distribution shown in the back. The visibility of the QND measurement is quantified by the values at the peaks and valleys, as indicated by I 0 , I 1 and I R (see Eq. (9)). The figure is for illustration only, and is not to scale To describe a realistic situation, we numerically simulate the case in which the parasitic resistances R, inductances L and the two bare capacitances C 0 differ from each other. In Fig. 3, we test the system with these asymmetries and the physical parameters given below. In the left plot, the role of a residual linear coupling g r is investigated. In the right one, we consider unbalanced resistances R ± δR, inductances L ± δL and capacitances C 0 ± δC. The results show that our analytical predictions accurately describe a system with non-zero g r and δC. Furthermore, the numerical points confirm that δR and δL enter as higher order perturbations. In fact, we generally find that Eqs. (5) and (6) are accurate for relatively large perturbations (up to 25%).
Inspired by recent experiments 15,35-38 , we estimate the value of λ, which can be reached in state-of-the-art setups. We consider a rectangular monolayer graphene membrane of length 1 μm and width 0.3 μm, with a mechanical frequency of ω m = (2π)80 MHz and a quality factor Q = 10 6 . It is suspended d 0 = 10 nm above a conducting plate, forming the capacitor (see sketch in Fig. 1a). Assuming that the membrane is clamped to the substrate along its boundaries, we identify the ratio of the coupling coefficients for each capacitor Cð ±xÞ in Fig. 1c to be g 2 /g 1 = π 2 x 0 = 8d 0 ð Þ 39 . Considering that for these geometries stray capacitances C s are typically preponderant with respect to C 0 , we take g 1 ≃ (2π)7 kHz and g 2 ≃ (2π)1 Hz, corresponding to C s ≃ 100C 0 . For comparison, a value of C s = 50 fF is obtained in ref. 35 , for a graphene membrane about two and a half times the size considered here. This stray capacitance would be 376 times C 0 ≃ 13 fF. Assuming a reduction of C s due to the smaller dimensions, we take C s = 100C 0 .
With an electrical reservoir at zero temperature n e ≃ 0 (valid for milliKelvin experiments), an electrical frequency ω s = (2π)7 GHz and decay rate γ t = (2π)150 kHz, we get λ b = 105 × Z out /R and λ p = 0.014 × (g 1 /g r ) 2 . Since the graphene coupling can be tuned via electric fields 40-42 , we assume g 1 /g r~1 00, which fixes λ between 60 (R = Z out ) and 122 (R = Z out /10), mostly restricted by λ p . This limit is well above the threshold for having a good visibility of the phonon number states (see below), and can be further improved by either increasing the sideband resolution ω m /γ t , the electrical frequency ω s or by reducing the size of the membrane. In Fig. 4b, we show the linear coupling g 1 as a function of the stray capacitance. For small values of C s , we reach the strong-coupling regime, where g 1 ≥ γ t . In the realistic scenario described above, where C s ) C 0 , our scheme still allows for phonon QND measurement even for g 1 ; g 2 ( γ t . This is in contrast to the typical optomechanical approach, where the quadratic interaction results from a hybridization of two optical modes, and strong coupling g 1 > γ t is required 28 . Regardless of how much C s reduces the coupling constants, it is in principle always possible to compensate by using stronger power. For details see the Supplementary Note 4 available in Supplementary Material online).
Measurement. We now evaluate how well a given value of λ allows for the QND detection of the phonon number. To this end, we consider a situation where the system is continuosly probed and measured. The output is then turned into discrete results by averaging over a suitable time T, and a histogram is constructed from the measured values V M . We assume that the heating of the continuous QND probing is in equilibrium with the mechanical damping and the associated reservoir. In this case, one also needs to consider the thermal bath of the membrane. In addition to Δn b determined above, the total heating out of the ground state is thus n m T. This additional term leads to a redefinition of the parameter λ to and the equilibrium average mechanical occupation, resulting from both the mechanical reservoir and the QND probe, becomes The phonon QND measurement is then characterized by λ′, which is desirable to have as close as possible to its maximum λ. This can be achieved by choosing a sufficiently strong probing power and a short measurement time T, such that the mechanical heating can be neglected. This leads to a large N eff , which does Given λ′, we now want to optimize all remaining parameters of the system, to be able to discern the ground and first excited states with the largest contrast. We simulate the mechanical system with the quantum-jump method, and pick Gaussian distributed random values for the electrical vacuum and thermal noise. From this, we make the histogram of the resulting output voltages V M presented in Fig. 5a, where the induced heating Δn b is optimized numerically. For the optimization we consider the visibility where I 0 and I 1 are the heights of the peaks corresponding to n b = 0 and n b = 1 phonons, while I R is the lowest height in between I 0 and I 1 (see Fig. 2). Additionally, we make an analytical model where we allow for one jump during each measurement period. We can extract the asymptotic behaviour of the visibility reflecting the compromise between the contributions to I R from the noise ∝ exp(−D 2 /8) and from the jumps during the measurements ∝ Δn b . The results of simulations and model are shown in Fig. 5a. The blue points are the numerical optimization, which are in good agreement with the analytical result (red, dotted line). Notice that for small values of λ′, the optimal Δn b is sufficiently high to allow multiple jumps during the measurement time T, leading to minor discrepancies. The black, solid line is Eq. (10), and the shadowed region corresponds to the predicted values of λ for the parameters introduced above. Qualitatively, clear signatures of the mechanical energy quantization are present for λ′ ≳ 40, where the visibility exceeds 20%.
For the experimental parameters considered above, the maximum attainable value of λ′ is λ = 122 (for R = Z out /10), and is achieved with a strong probe such that N eff ) n m . The incident power and the measurement time T provide a handle to optimize the performance for given experimental conditions. Qualitatively, a short value of T minimizes the effects of the mechanical heating, and makes λ′ ≃ λ. On the other hand, the required power to reach such a regime can be troublesome 43 , and we may need to integrate for too long time to have sufficient statistics (since N eff ) 1). This last problem can be solved by adding an electrical cooling, red-detuned by ω m ) γ t from the QND probe. This cooling would not affect the parameter λ′, since it does not heat up the system, but only reduces N eff . The visibility ξ thus remains almost unaltered (see Eq. (10) and Fig. 5b), but the probability to find the membrane in low excited states is increased, reducing the experimental time.
As an example, assume that the heating from the electrical feedback and the mechanical bath are equal, such that λ′ = λ/2 = 61. Considering a cryogenic temperature of 14 mK 37 , the average mechanical occupation is n m ≃ 3, implying N eff ¼ 6. The optimal Δn b is then 0.3, and can be obtained with a driving power of 16 nW and a measurement time of 0.1 ms for a mechanical quality factor Q = 10 6 and a stray capacitance C s = 100C 0 . For other values of Q and C s , the driving power can be varied to fulfil the constraint N eff ¼ 2 n m , as shown in Fig. 4a. The incident field is rather intense, which may cause additional heating to the system. In the set up of ref. 43 , such additional heating has been observed above an intracavity photon number of 10 8 . For comparison, in Fig. 4a we show the intracavity photon numberα j j 2 =γ t for our system, whereα j j 2 ¼ α j j 2 =T is the photon flux. Depending on the parameters, we see thatα j j 2 =γ t will be similar or higher than 10 8 for C s ≳ 100C 0 . These devices cannot, however, be compared directly. Nevertheless, since ref. 43 indicates that the source of this heating is electrical, we believe that it would be strongly suppressed for the QND measurement considered here. Since the linear coupling is almost cancelled by symmetry, the resulting heating rate is likely reduced by a factor (g r /g 1 ) 2 ≃ 10 −4 . In absence of this suppression, conducting our experiment in a pulsed regime may substantially reduce other heating mechanisms.

Discussion
We have revisited the challenge of performing a phonon QND measurement. Employing symmetry to inhibit the linear coupling, the detrimental heating is suppressed while retaining the desired quadratic coupling. Contrary to the generally studied optomechanical case 28 , the residual coupling to the antisymmetric mode is strongly suppressed by its higher frequency and reduced resistance. A particularly attractive feature of the current approach is that it is only sensitive to the ratio g 2 /g 1 , and not to their absolute values. Stray capacitances, which reduce the electromechanical couplings, can thus be compensated using stronger input fields.
These attractive features put QND detection within reach of presently available technology. A successful realization of a QND detection will not only represent a demonstration of genuine nonclassical behaviour of mechanical systems, but also extend the interactions available in electro/optomechanics to non-Gaussian operations 44 . This will considerably expand the realm of effects that can be studied with these systems, and facilitate their application for quantum information processing 23 . a Average intracavity photonsα j j 2 =γ t required for the QND measurement, as a function of the relative value of the stray capacitance C s /C 0 . The three lines correspond to different values of the mechanical quality factor, as indicated in the legend. We assume Δn b = 0.3 and equal contributions from the mechanical and electrically induced reservoirs n m ¼ N eff =2 ¼ 3. As a reference, the grey dashed lines indicate the associated powers of the probe. b Linear coupling g 1 as a function of C s /C 0 . For both figures, the shadowed region indicates the strong coupling g 1 ≥ γ t , where QND detection is feasible with other approaches 26,31,32 As an outlook, it is desirable to extend this work to the optomechanical case, where mechanical systems with a similar quadratic coupling have recently been studied 33,45,46 , but conditions to have a successful phonon QND measurement have not been yet determined. The electromechanical systems considered here can be described with Kirchoff's laws, which give rigorous results within a well-defined model. The physical mechanisms behind the heating are identified to be the Johnson-Nyquist noises associated to the resistors, and fabrication imperfections. For comparison, the exact description of dissipation in a multimode optomechanical system may be more involved. Nevertheless, the results presented here could be useful for guiding the intuition towards QND detection in the optical regime. As a further extension, it would be interesting to investigate the effect of squeezing. By reducing the vacuum noise, squeezing can lead to a direct improvement in λ, thus reducing the physical requirements for the QND detection.

Methods
The double-arm circuit. The Hamiltonian for the system in Fig. 1c where subscripts "a" and "s" indicate the asymmetric and the symmetric electrical fields, respectively. From Eq. (11) and using Kirchoff's laws, it is possible to determine the equations of motions, including noises and decays. The normalized distance D 2 = d 2 /σ 2 is obtained assuming the phonon number to be constant within T-that is: setting g 1 = 0-so that the asymmetric and symmetric fields decouple. Looking at the phase quadrature of the reflected signal V out ¼V in À γ tΦs , we determine d. The noise σ is the sum of vacuum noise from the input coherent field, and the Johnson-Nyquist noises of the resistors. As discussed above, the heating Δn b ¼n b ðTÞ h ihas two contributions: asymmetries leading to a non-vanishing linear coupling g r , and the charge redistribution. The first is found by assuming R ( R 0 and L ( L 0 , such that the circuit in Fig. 1c is equivalent to the one in Fig. 1b, for which we already know Δn b . The contribution from charge redistribution is determined from the Hamiltonian in Eq. (11) neglecting the quadratic interaction, which does not alter the phonon number. The strongly driven symmetric electrical field is then substituted with its steady state, obtained assuming a constant photon flux. The time evolution ofn b h i is finally found by looking at the equations of motion for the asymmetric field and the mechanical creation/annihilation operators. With the amplitude of the symmetric mode replaced by its steady state, these equations are now linear in the annihilation (creation) operatorsb ðb y Þ and can be solved by standard optomechanics techniques.
Asymmetric circuit. To obtain Fig. 3, we analyse the system in the presence of asymmetries. First, we derive the generalization of the Hamiltonian in Eq. (11) with unequal rest capacitors, resistors, inductors and linear couplings. Differently from above, we linearize the symmetric/asymmetric electrical fields around their mean values (Q a=s ! hQ a=s i þδQ a=s andφ a=s ! hφ a=s i þδϕ a=s ), and the mechanical creation/annihilation operators ðb ðyÞ ! hb ðyÞ i þδb ðyÞ Þ. Here, besides the usual oscillatory behaviour of the mechanical operators hb ðyÞ i, the amplitude is generally time dependent 47 . This can be understood by looking at Eq. (11); since both the electrical fields have now non-zero average, the three body interaction /Q aQs ðb þb y Þ is equivalent to a force directly driving the mechanical system. Once solutions for the averages are found, it is possible to determine the variations, and finally the time evolution of the phonon number.
Optimization of the visibility. To obtain Fig. 5 we rely on both an analytical and a numerical optimization of the visibility ξ. To determine the red, dotted curve, we assume that the initial mechanical state is thermal, such that the occupations of the Fock states can be found. Given Δn b , the probability to jump once either up or down during the measurement time T is a Poissonian process. The probability distribution function for the outcomes V M can then be obtained and maximized, by varying Δn b . The histograms and the blue points are derived with Monte Carlo simulations, where the time evolution of single mechanical trajectories are replicated with the stochastic wave-function method 48 . Importantly, every measurement interval of duration T has been discretized, to allow for multiple jumps. The parameter Δn b is then varied to find the best visibility ξ.

Data availability
All material related to this work can be found at https://sid.erda.dk/share_redirect/ eUaGoI8JbN.