A microwave scattering spectral method to detect the nanomechanical vibrations embedded in a superconducting qubit

Nanomechanical resonators (NMRs), as the quantum mechanical sensing probers, have played the important roles for various high-precision quantum measurements. Differing from the previous emission spectral probes (i.e., the NMR modified the atomic emission), in this paper we propose an alternative approach, i.e., by probing the scattering spectra of the quantum mechanical prober coupled to the driving microwaves, to characterize the physical features of the NMR embedded in a rf-SQUID based superconducting qubit. It is shown that, from the observed specifical frequency points in the spectra, i.e., either the dips or the peaks, the vibrational features (i.e., they are classical vibration or quantum mechanical one) and the physical parameters (typically such as the vibrational frequency and displacements) of the NMR can be determined effectively. The proposal is feasible with the current technique and should be useful to design the desired NMRs for various quantum metrological applications.


Introduction
In recent years, nanomechanical resonators (NMRs) have been recently highlighted for various precise measurements, as their mechanical vibrations can reach very high frequency such as up to the GigaHertz (see, e.g., 1 ).This makes them be directly utilized as the electronic devices for radio communications and various precise measurements such as for sensing masses, weak forces, and charge, etc. 2,3 .In fact, the ability of the resonator to detect the physical quantities is closely related to its resonant frequency.For example, the mass-loaded sensitivity can be written as s = δ f n /δ m = f n /2m [4][5][6][7] , if the resonator with the eigenfrequency f n and mass m is added by a mass δ m, which leads to the frequency shift δ f n .This implies that the higher frequency of the resonator corresponds to the stronger ability to detect the smaller masses, and also the vibrational frequency of the NMR should be precisely calibrated beforehand.
Physically, a classical oscillator has a well-defined amplitude of motion; while, for the quantum oscillator, its displacement z is related to the vibrational quantum state with the quantum fluctuations of the momentum and displacement.In fact, by various cooling techniques, the vibrations of the NMRs can be cooled to the approaching quantum ground state, starting from a thermal state 8 .As a consequence, the NMR with the GigaHertz vibrational frequency can be served as the coherent quantum devices to test the fundamental principles in quantum mechanics (e.g., the Heisenberg uncertainty relation, quantum superposition, and macroscopic quantum, etc.) and also generate various hybrid quantum systems (e.g., coupling it to the superconducting circuits 9,10 , semiconducting quantum dots 11,12 , and NV centers 13,14 , etc. 15 ) for implementing the desired quantum metrologies and quantum information processings [16][17][18] .Therefore, characterizing the physical features (i.e., the vibration is classical or quantum mechanical) and further measuring the relevant physical parameters (typically such as the vibrational frequency and displacement) of the vibration of the NMR are particularly important.
Basically, compared with the calibrations of the physical parameters of the classical mechanical vibration reviewed below, the precise measurement of the vibrational displacement of a quantized high frequency harmonic resonator (HO) is still an open problem.This is because that, besides the various background noises, the internal quantum fluctuations of the vibration plays a key role.For example, the amplitude A 0 of a quantized HO at vibrational ground state |0 is directly determined by the quantum fluctuation of the measured vibrational displacement, i.e., A 0 = ρ z / √ 2 with ρ z = 0|ẑ 2 |0 − 0|ẑ|0 2 .Also, the sensitivity of the frequency measurement of a quantized HO is also limited by the uncertainty relation between the quantized vibrational energy and the lifetime of the operated energy stationary state.Up to our knowledge, a few methods have been demonstrated to detect the displacement of the quantized vibration of the NMR at low temperature.For example, the vibrational frequency of the quantized NMR, which is embedded in a superconducting transmission line resonator (STLR), can be measured by using a microwave interferometer configuration 19 .The motion of the quantized vibrations of the NMR can be measured and control by probing the optical-mechanical effects in the cavity optomechanical system 20 .Typically, in Ref. 21 we proposed an effective approach to probe the tiny motion of the NMR by coupling it to a half-wavelength STLR, mediated by a SQUID-based qubit.In that configuration, the quantized mechanical vibration of the NMR modifies the energy structure of the qubit and thus its spontaneous mission spectrum, which can be indirectly detected by the spectral measurement of the STLR.Interestingly, the vibrational features, i.e., the vibration is quantum mechanical or classical, can be identified by observing the modifications in the spectrum spontaneously emitted from the SQUID-based qubit.Alternatively, in the present work we propose an active approach (rather than the passive one in Ref. 21 ) to probe the mechanical vibrations of the NMR by measuring the transmission spectra of the travelling microwave along a one-dimensional transmission line.The detected NMR in the present configuration is embedded in a rf-SQUID-based superconducting qubit.As a consequence, the vibrational frequency and displacement of the NMR can be estimated by observing a few specific frequency points in the transmitted spectra of driven travelling microwave scattered by the qubit.Furthermore, we show that the proposal works also for the alternative measurement configuration, i.e., the scatter of the microwave is replaced by a quarter-wavelength STLR (with a sufficiently-high quality factor), which is inductively coupled to the NMR via the qubit.Due to the use of the STLR, the electromagnetically induced transparency-like effects in microwave band are modified and thus the relevant parameter estimations could be more conveniently achieved with the more observable data.Importantly, a recent experiment 22 had demonstrated the resolution detection of the energy levels of a NMR by using the scattering spectral measurements of the travelling microwaves.
We propose a spectral approach, by probing the transmitted and phase shift spectra of the travelling wave scattered directly by the qubit-NMR device, to estimate the physical parameters of the NMR embedded in a rf-SQUID qubit.Then we treat the problem with a more complicated one, i.e., the STLR-qubit-NMR system, and demonstrate the corresponding spectral measurements of the NMRs.One can see that, with these spectra the desired physical parameters can be more easily estimated by using the more observable data, due to the microwave electromagnetically induced transparency-like effects.Finally, in order to the completeness, we review how the physical parameters of a classical mechanical resonator were measured and give he derivations of the relevant Hamiltonians.

Measuring the frequency and amplitude with a mircowave driven qubit-NMR system
The NMR considered here is generated by the mechanical vibration of the part of a flux-biased rf-SQUID loop, which generates a qubit encoded by the two lowest energy eigenstates of the loop.Another magnetic field B 0 can be applied along the loop plane to provide a restoring force for generating a mechanical vibration of the part of the loop, i.e., the NMR oscillates along the direction perpendicular to the loop plane.The significantly weak vibration of the NMR along the direction parallel to the loop plane 23,24 can be omitted for the simplicity.Certainly, the other configurations can also be utilized to realize the qubit-NMR couplings, see, e.g., in Refs. 25,26 The qubit can served as a probe to measure the physical features of the coupled NMR by using the scattering measurements of the travelling microwave.Without loss of the generality, let us consider the specifical configuration shown schematically in Fig. 1.The Hamiltonian of the system can be expressed as (h = 1) where describes the quantized traveling microwave (with the group speed v g ) transporting along the feedline 27 .ĉL/R (x) and ĉ † L/R (x) are the x-dependent annihilation and creation operators of the left/right-moving microwave photons, respectively.Next, describes the interaction (with the strength V 1 ) between the microwave photons transporting along the feedline and the qubit 27 .σ † and σ− are the Pauli operators of the qubit.As the wavelength of the microwave is significantly longer than the scale of the rf-SQUID loop, the interaction between them can be treated as a δ -function, taking place at x = 0. Thirdly, the Hamiltonian Ĥs , describing the NMR and its coupling to the qubit, takes the forms depending on the specifical features of the vibration of the NMR 21 .Theoretically, the microwave scattering features of the qubit-NMR system can be calculated by solving the time-dependent Schrödinger equation: i∂ |ψ(t) /∂t = Ĥ1 |ψ(t) , for the usual elastic scattering, the problem becomes to solve the stationary Schrödinger equation: by letting |ψ(t) = e −iωt |ψ .Here, ω = v g k (k being the wave vector of the applied microwave transporting along the feedline.

Measuring the eigenfrequency of the qubit
First, if the NMR is absent, i.e., the magnetic field B 0 is not applied, then Ĥs in Eq. ( 1) is nothing but the Hamiltonian of the rf-SUID-based qubit: with ω 0 being the energy of the qubit's excited state |1 .In this case, the generic solution to the equation ( 4) can be expressed as where |φ 0 = |0, 0 represents that the electromagnetic field in feedline is at the vacuum and the qubit is at the ground state.Instituting Eqs.(2-3) and (5-6) into Eq.( 4), we get: For the sake of the convenient calculation, we assume that 28 where t/r is the transmission/reflection amplitude of the travelling-wave photons.Instituting Eq. ( 8) into Eq.( 7), the transmission amplitude t can be solved as Here, the relevant parameters are typically set as: ω 0 = 2.1 × 10 9 Hz, γ c = 3.3 × 10 7 Hz, and δ ω = 6.6 × 10 7 Hz.
which is dependent of the frequency of the driven microwave.Consequently, the transmitted spectrum of the microwave photons scattered by the rf-SQUID qubit, without the NMR vibration, can be calculated as Correspondingly, the phase shift spectrum of the transmitted microwave reads It is seen from Fig. 2(a) that, the transmitted tip is at the microwave-qubit resonant point, i.e., the driving microwave with the frequency ω = ω 0 is completely reflected.Therefore, by observing the frequency at the dip in the transmitted spectrum of the travelling microwave scattered by the qubit, the eigenfrequency ω 0 of the qubit can be determined.Certainly, due to the coupling dissipation γ c , the transmitted dip is not the δ -function.Instead, it shows a Lorentz tip shape with the full width at half minimum (FWHM): δ ω 0 = 2γ c .Experimentally, such a FWHM could be served as the uncertainty of the observed dip.In this model, we omitted the internal dissipations of the qubit and treated γ c is the total dissipation of the qubit.On the other hand, Fig. 2(b) shows that, if ω = ω 0 the phase of the reflected microwave is shifted a π-phase, which is independent of the dissipation of the qubit.

Measuring the vibrational frequency of a quantum mechanical NMR
As shown schematically in Fig. 1, if the magnetic field B 0 is applied, then the vibration of the MNR is excited.Furthermore, let us assume that the vibration of the embedded NMR has been cooled into the quantum regime, i.e., the NMR is treated as a quantum mechanical oscillator (called as the QNMR afterwards) and described by the bosonic operators b and b † .By a long but direct derivation, the Hamiltonian Ĥs in Eq. ( 1) can be effectively expressed as with g Q being the qubit-QNMR coupling strength (See Methods), B 0 is the applied magnetic field in Fig. 1, I p is the amplitude of the circular supercurrent along the rf-SQUID qubit loop, l is the length of the QNMR.Accordingly,the generic solution of Eq. ( 4) can be expressed as Here, |φ Q = |0, 0, 0 refers to the scattering ground state, i.e., the electromagnetic field in feedline is at the vacuum, the NMR is at the vibration ground state, and the qubit is prepared at its ground state |0 .Instituting Eqs.(2-3) and (12-13) into Eq.( 4), one can easily proof that the coefficients in Eq. ( 13) are determined by With the same method used in the above subsections, we get the transmitted spectrum and also the phase shift spectrum: of the travelling microwave in the feedline, respectively.Fig. 3 shows the calculated transmitted-and phase shift spectra for the case wherein the vibration of the NMR embedded in the qubit is quantum mechanical.
Interestingly, Eqs. ( 15) and ( 16) imply also that, if the frequency ω of the travelling microwave is equivalent to ω b , i.e., the vibrational frequency of the QNMR, then a frequency point with |t Q (ω)| 2 = 1 can be observed between two dips in the transmitted spectrum.This indicates that, the frequency ω b of the QNMR could be directly determined by observing the frequency point of the travelling microwaves without any reflection, i.e., the frequency of the microwave is completely transmitted without any phase shift.To measure the vibrational displacement of the QNMR, we need to detect the qubit-QNMR coupling strength g Q .This can be achieved as follows.First, one can see that two transmitted dips centered respectively at ω + and ω − are observed in the spectrum shown in Fig. 3(a).They are determined by solving the equation: T Q (ω) = 0, and read: As a consequence, the qubit-QNMR coupling strength between the qubit and QNMR can be calculated as whose estimated accuracy depends on those of the measured frequency points ω b and ω ± .From the simulated spectra shown in Fig. 3(a), the frequency uncertainty of the two dips are observed as: δ ω − = 1.9 × 10 7 Hz for ω + = 2.162 × 10 9 Hz, and δ ω + = 4.8 × 10 7 Hz for ω − = 1.938 × 10 9 Hz.While, from Eqs. ( 15) and ( 16), the observed eigenfrequency ω b of the NMR could be threaten as the precise value.

Measuring the amplitude (i.e., phonon number) of a quantum mechanical NMR
Next, with Eq. ( 70), we know that either the amplitude of the applied magnetic field or the mass of the NMR can be determined as: Then, by using the qubit-QNMR coupling strength g Q measured above, the mean displacement of the vibrational QNMR can be determined as which is obviously dependent of the quantum state of the system, not only the vibrational quantum state of the QNMR.Therefore, the mean displacement of the QNMR, including the influence from the quantum fluctuation, can be determined.Certainly, such a measurement is a quantum demolition one, as the energy exchange takes place frequently between the qubit and the QNMR.However, following Ref. 22, the quantum nondemolition measurement of the quantized vibration of the NMR could be implemented with the present configuration.This can be achieved by adjusting the eigenfrequency ω 0 of the qubit to let the qubit-QNMR work in the dispersive regime, i.e., g Q /∆ ≪ 1, with ∆ = ω 0 − ω r .under this condition, Ĥs in Eq. ( 1) can be reduced as instead Ĥq−QNMR in Eq. ( 12).Correspondingly, the generic solution to the Schrödinger equation Eq. ( 4) can be written as with |φ ′ Q = |0, 0 being the scattered ground state and |ψ QNMR the vibrational quantized state of the QNMR, which will be nondemolition during the scattering spectral measurements.The coefficients in the above wave function are determined by where n = ψ NMR | b † b|ψ NMR is the average phonon number of the vibrational QNMR.Solving Eq. ( 22) similarly, we get the relevant transmitted-and phase shift spectra: and respectively.It is seen schematically from Fig. 4 that, the dips in the transmitted spectrum and the phase shifts, near the completely reflected frequency points, are really related to the average number of phonon n of the QNMR.The centre frequency of the transmission spectrum.The shifted frequency of the dip is dependent on the phonon number of the QNMR.This implies that the average phonon number of the QNMR could be measured by probing how the frequency of the dip is shifted.Again, the FWHM ∆ω = 2V 2 1 /v g of the observed dip limits the measurement accuracies of the phonon numbers of the QNMR.Experimentally, to differentiate the frequencies corresponding to the nearest two dips induced respectively by the phonon states |n and |n + 1 , the condition: should be satisfied.

Measuring the vibrational amplitude of a classical NMR
Physically, due to the unavoidable dissipation, the quantum feature of the NMR would be lost and the quantum vibration of the NMR becomes the classical one.For the classical NMR (called as the CNMR later) with the amplitude A C , the Hamiltonian Ĥs in Eq. ( 1) reads which is simply obtained by replacing the q-number Bosonic operators b and b † in Eq. ( 12) as the c-number quantities.Here, g C = B 0 lI p A C is the coupling strength between the qubit and the CNMR and A C is the amplitude of the CNMR.In the rotating frame defined by the transformation the Hamiltonian in Eq. ( 26) can be rewritten as Here, σ † and σ− are the Pauli operators of the qubit.The Hamiltonian ( 27) can be easily diagonalized as with σz = | 1 1|− | 0 0| = (ω 0 + ω b ) σz /(2 ω0 )+ g C σx / ω0 .Obviously, the embedded CNMR just modifies the eigenfrequency of the qubit without the NMR.Therefore, replacing just the ω 0 , in the Eqs. (10)and (11), by ω0 , the transmitted and phase shift spectra of the travelling microwave scattered by the present qubit-CNMR system can be easily expressed as respectively.This indicates that, the spectral shape is the same as the situation that the microwave scattered by the qubit without the NMR, but the resonance point is shifted from ω 0 into ω0 .Typically, one can see that, differing from the two dips in Fig. 3(a) for the qubit-QNMR scattering, Fig. 5(a) shows a single dip in the transmitted spectrum of the travelling wave microwave.Therefore, the vibrational feature of the NMR, i.e., either the quantum mechanical or classical vibration, could be calibrated by observing the transmitted spectrum of the travelling microwave scattered by the qubit-NMR; if the observed spectrum has two dips, then the vibration of the NMR is quantum mechanical, while the CNMR refers to one dip spectrum.
Given the vibration frequency ω b of the NMR and also the eigenfrequency ω 0 had been measured, the qubit-CNMR coupling strength can be calculated.Then, the vibrational amplitude of the CNMR can be determined as: With the set parameters, we get g C = 9.05759 × 10 7 Hz, we refer to the following parameters in 6 and assume the applied magnetic field to be B 0 = 5mT, yielding the vibrational amplitude of the CNMR: A C ≃ 2nm.
In the above configuration, wherein the travelling microwave is scattered by the rf-SQUID loop qubit, whose coherence might be easily broken by the driving microwave.Alternatively, in the following we consider another configuration, wherein a transmission line resonator (TLR) is introduced to isolate the qubit from the driving microwave and let the later be scattered by the TLR.We shows that, with such an TLR-qubit-NMR configuration, the physical parameters of the NMR can also be measured.
Improving the accuracies with a driven TLR-qubit-NMR system Following Ref. 22 , we now consider the configuration shown in Fig. 6, wherein the traveling-wave is scattered by the STLR, rather than the qubit.The NMR is still embedded in the rf-SQUID-based qubit.The Hamiltonian of this system can be expressed as with Ĥr = ω r â † r âr (33) describing the fundamental-mode standing-wave photons (with the frequency ω r ) in the STLR, â † r and âr represent the generation and annihilation operators of the photons in TLR respectively.The interaction between the travelling wave photons transporting along the feedline and the standing-wave photons in the resonator reads where V 2 is the coupling strength between the travelling wave photons in the feedline and the standing wave photons in the STLR 29 .The coupling between the photons in resonator and the qubit reads with g rq being the STLR-qubit coupling strength.Similarly, the transmitted spectra of the travelling waves scattered by the present STLR-qubit-NMR system can be calculated by solving the stationary Schrödinger equation: for various vibrational forms of the NMR, i.e., the different forms of the Ĥs shown in Eqs.(5,12,26).Noted that, the spectrum of the travelling wave scattered by a single quarter-wavelength STLR had been calculated exactly in Ref. 29 and verified by a series of experimental measurements 30 .Certainly, if the travelling-wave microwave transporting along the feedline is scattered only by the STLR, a single dip centred at ω = ω r could be observed, with the FWHMs γ r being described by the quality factor of the STLR.Next, we will investigate how the spectra are modified by the scatterings of the STLR-qubit-NMR system.

The qubit eigenfrequency measurement
First, if the NMR is absent, i.e., Ĥs in Eq. ( 32) reduces to Ĥq in Eq. ( 5), then the generic solution of the Schrödinger equation (36), with Hamiltonian (32), can be written as Here, | φ0 = |0, 0, 0 is the ground state of the system with the qubit being at the ground state |0 and the electromagnetic fields in feedline and in TLR at electromagnetic vacuum.The coefficients in the above generic wave function are determined by

9/19
and can be analytically solved as: As a consequence, the spectra of the transmitted and phase shifted of the travelling microwave can be calculated as respectively.Here, the relevant paramters are set as: ω r = ω b = 2 × 10 9 Hz, v g = 3 × 10 8 m/s, ω 0 = 2.1 × 10 9 Hz, V 2 = 10 8 Hz, g rq = 10 8 Hz.The red and yellow lines refer to without and with the CNMR, respectively.
It is seen from the transmitted spectrum shown in Fig. 7 that, due to the coupling of the qubit, the single dip centered at ω = ω 0 in the spectrum of the travelling microwave scattered only by the STLR 29 is splitted as the two dips centered at ω = ω ′ ± , with This is the typical vacuum Rabi splitting phenomenon.Interestingly, in this case, Fig. 7 and also Eq. ( 40) indicate that the eigenfrequency ω 0 of the qubit can be determined by observing the frequency point, at which the input microwave is completely transmitted.The observed phenomenon in the spectra, i.e., the frequency at the original completely reflected dip (scattered by a single STLR) is changed as the completely transmitted point of the microwave scattered by the STLRqubit system.This could be called as the qubit-induced Electromagnetically-induced-transparency (EIT)-like effect of the electromagnetic waves, wherein the qubit is served as the control field to modify the energy structure of the scatter (i.e., the STLR here).Again, by observing the ω ′ ± , the STLR-qubit coupling strength can be determined as: Typically, if the qubit is resonance with the TLR, then distance between the two dips is ω + − ω − = 2g rq .

10/19
Vibrational frequency measurement of the QNMR Now, let us consider the situation that the vibration of the embedded NMR is quantum mechanical.In this case, the Hamiltonian in Eq. ( 32) reads Ĥq−QNMR shown in Eq. ( 12).The generic solution to the corresponding Schrödinger equation is expressed as Here, | φQ = |0, 0, 0 b , 0 represents the ground state of the present TLR-qubit-NMR system, which means that the electromagnetic fields in the feedline and the TLR are both in vacuum, the quantized vibration of the NMR is cooled to the vibrational ground state |0 b and the qubit is prepared at the ground state |0 .The coefficients in Eq. ( 36) are determined by They can be analytically solved as: with Again, the transmitted and phase shifted spectra of the travelling microwave are calculated as respectively.One can see from Fig. 8 that, the quantized vibration of the NMR significantly changes the transmitted and phase shifted spectra of the travelling microwave, typically inducing two EIT-like transparent windows, wherein the travelling microwaves are completely transmitted for At these frequency points, the phase shifts of the travelling microwave are zero.These behaviors are significantly different from the cases, wherein the vibration of the NMR is absent, and provide more data to measure the physical parameters of the NMR.For example, with the eigenfrequency ω 0 (the observed completely transmitted frequency ω = ω 0 in Fig. 7) and the two completely transmitted frequency points ω ′′ − and ω ′′ + (shown in Fig. 8), Consequently, the vibrational frequency of the NMR can be estimated as Furthermore, with the observed frequencies: ω 0 and ω ′′ ± , shown in Figs.7 and 8, the qubit-QNMR coupling strength could be calculated as: similarly, the vibrational displacement and also the phonon number of the QNMR can be measured.

Vibrational amplitude measurement of the CNMR
If the qubit is embedded by a CNMR, then Ĥs is taken as Ĥq in Eq. (28).Consequently, the spectra of the transmitted and phase shifted of the travelling microwave can be calculated as respectively.Above, In Fig. 9 we shows specifically the transmitted and phase shifted spectra of the travelling microwave scatted by the STLR coupled to the qubit embedded by the classical NMR.As mentioned in Fig. 3, due to the coupling of the qubit, the transport feature of the travelling microwave scattered by the STLR shows the electromagnetic-induced-transparent (EIT)-like behavior; a transparent window centered at ω = ω r is generated between the two dips (centered at ω = ω ± ) in the transmitted spectra.In the present case, due to the existence of the classical vibration of the NMR embedded in the qubit, the original EIT-like phenomena (without the NMR) is modified, i.e., the completely transmitted frequency point is shifted into ω = ω0 , although the width of the transparent window is unchanged.This means that, the modified EIT-like effect can be served as the evidence of the existence of the classical vibration of the NMR.Fig. 10 shows the comparison between classical vibration (a) and quantum vibration (b) of the system before and after adding STLR.In general, with the addition of STLR, the trough width of transmission spectrum decreases significantly under the same parameters, which is very important for the accuracy of actual measurement.Here, the relevant paramters are set as: ω r = ω b = 2 × 10 9 Hz, v g = 3 × 10 8 m/s, ω 0 = 2.1 × 10 9 Hz, V 2 = 10 8 Hz, g rq = 10 8 Hz.The red and yellow lines refer to without and with the CNMR, respectively.

Discussion
In conclusion, we proposed a spectral measurement method to detect the vibrational and displacement of the NMR, embedded in the rf-SQUID-based qubit.By observing certain specifical frequency points in the measured transmitted and phase shift spectra of the travelling microwave scattered by either the qubit-NMR system or the STLR-qubit-NMR one, we showed that the vibrational frequency and vibrational displacement of the NMR can be measured effectively.By observing the completely transmitted frequency points to determine the qubit-NMR coupling.Interestingly, the proposal provides a quantifiable way to identify the vibrational feature of the NMR, i.e., is the vibration classical or quantum mechanical.In fact, this is the key problem of the NMR being used to implement various precise measurements and quantum information processings.
In the present work, we only consider the coupling dissipation and omitted all the internal dissipations of the devices in the system, typically including the qubit, STLR, and the NMR.By simply introducing the relevant non-Hermitian term in the Hamiltonian to phenologically describe the relevant dissipations, we argued that the method proposed here should still be effective.

Methods
The spectral method for a classical harmonic oscillator Historically, the stationary spectral method has been widely applied to calibrate the vibrational frequency ω b and displacement z of the classical harmonic resonator (HO).First, by applying a driving force F(t), the equation of the motion of the HO with the internal dissipation γ can be expressed as 31 : with m, ω b , and γ being the mass, frequency and the dissipative coefficient of the HO, respectively.If the applied force is periodic, i.e., F(t) = a cos(ω d t) (with the amplitude a and frequency ω d ), the stationary solution (for t ≫ 1/γ) of the dynamical equation (1) reads: , with the amplitude-and phase spectra: As a consequence, by observing the amplitude-and phase spectra under the different frequency driving, the frequency ω b of the HO can be determined.Specifically, from the amplitude spectrum, ω b can be estimated as the peak value frequency with the accuracy δ ω b = γ (i.e, full width at half maximum), which means that the higher precision of the estimation can be 13/19  Hz, v g = 3 × 10 8 m/s, ω 0 = 2.1 × 10 9 Hz, V 2 = 10 8 Hz, g rq = 10 8 Hz, g C = g Q = 1 × 10 8 Hz.
obtained for the lower dissipation.While, from the phase spectrum we have φ (ω d ) = π/2 for ω d = ω b , which is independent of the dissipation of the HO.Next, two categories: (i) coupling it directly to the sensor 32,33 , and (ii) using the remote radio-wave or optical interferometries [34][35][36] , are usually applied to detect the vibrational displacement of the HO.However, various unavoidable noises, typically such a the stochastic force ξ (t), limit the sensitivity of these displacement detections.This is because that, in the noise background a dissipative HO is described by the Langiven equation 37 : with ξ (t) = 0, but ξ (t)ξ (t ′ ) = 0. Simply, for the white noise is the Fourier transform of the correlation function of the stochastic force reads: S ξ (ω) = ∞ −∞ dte iωt ξ (t)ξ (0) = 2mγk B T .Consequently, the spectrum of the detected displacement is obtained as 38 with the spectral density: Obviously, the reachable sensitivity of the displacement measurement is related to the environment temperature T , dissipation parameter γ and also the measurement bandwidth.
Physically, although the influence of the noises can be reduced by developing various techniques, typically such as the resonance force microscopy techniques 39 and the fluid viscosity 40 , the accuracies of the parameter measurements are very limited as the used probers behave still the classical motions.Basically, the possible improvements should be achieved by using the quantum mechanical probers 22 .In the following, we discuss how to implement such an improvement by using the qubit and quantum oscillator as the probers.

The derivations of Ĥs for various cases
In this section, we provide the derivations of Ĥs for various cases, in detail.

Hamiltonian of the rf-SQUID based qubit
For a flux-biased rf-SQUID loop, the Lagrangian can be expressed as where Φ and L are the total flux and inductance of the rf-SQUID loop, respectively.C J and I c are the capacitance and critical current of the Josephson junction, respectively.Φ e is the biased magnetic flux and Φ 0 = h/(2e) the flux quanta 41 .Defining the canonical momenta we have the classical Hamiltonian of the flux-based rf-SQUD loop, Formally, this Hamiltonian is equivalent to that for describing the motion of a particle with the "mass" m = C J moving in the potential: Noted that the supercurrent in the loop is determined by 42 .This means that a clockwise supercurrent (i.e., I < 0) is generated along the loop if Φ e > Φ, while for Φ e < Φ an anti-clockwise supercurrent (i.e., I > 0) is generated.After the usual canonical quantization, the classical Hamiltonian H0 becomes the quantized Hamiltonian: Specifically, for Φ e /Φ 0 = 0.5 the potential shows a symmetric double wells around the point Φ = Φ 0 .With the typical parameters: πI c L J /Φ 0 = 1,C J = 1.7 × 10 −14 F, and L J = 6 × 10 −9 H, the eigenvalue problem of the Hamiltonian Ĥ0 can be numerically solved and the lowest two eigenvalues are: E 0 = 2.7025 × 10 −23 J and E 1 = 2.7225 × 10 −23 J, with the transition frequency between them being ω 0 = (E 0 − E 0 )/h = 2.1GHz.Fig. 11 shows the corresponding eigenfunctions, marked respec- tively by |0 and |1 , of the corresponding eigenvalues E 0 and E 1 .Physically, a rf-SQUID based qubit can be encoded by the lowest energy eigenstates |0 and |1 of the Hamiltonian Ĥ0 .Its free evolution can be described by the Hamiltonian shown in Eq. ( 5).In fact, limited in the subspace of the qubit, we have This is because that, for the typical parameters: 0|H 0 |1 = 1|H 0 |0 = 8.02 × 10 −31 J, is much less than 0| Ĥ0 |0 = E 0 = 2.7025 × 10 −23 J and 1| Ĥ0 |1 = 2.725 × 10 −23 J. Fig. 11 shows the distributions of the eigenfunction of the qubit states |0 and |1 , respectively.One can easily seen from Fig. 12 that, the states are respectively localized in the symmetric double wells.Therefore, they refer to the clockwise and anti-clockwise supercurrent states 43,44 , respectively.Here, we provide the derivation of the Hamiltonians Ĥq−QNMR in Eq. ( 12) and Ĥq−CNMR in Eq. ( 26), respectively.As shown in Fig. 1, the magnetic field B 0 can be applied to excite the vibrational of the NMR along the z-direction.The Lagrangian of the rf-SQUID loop embedded by a NMR can be written as with m, ω b and z being the mass, vibrational frequency and displacement of the NMR, respeectively.B 0 is applied along the y-direction to provide a restoring force for exciting the mechanical vibration of the NMR.The corresponding Hamiltonian reads: under the usual rotating-wave approximation (RWA).Therefore, the Hamiltonian Ĥs in Eq. ( 1) becomes that in Eq. ( 12).Obviously, if the vibration of the NMR is classical, then the displacement of the CNMR can be represented as a c-number (instead the q-number): z = A C cos(ω b t) with A C being the vibrational amplitude.As a consequence, under the RWA the interaction between the qubit and the CNMR can be expressed as: with g C = B 0 lI p A C being the qubit-CNMR coupling strength.Replacing the ĤI in Eq. ( 69) by Ĥ′ I and removing the pure c-number term, Hamiltonian in Eq. (69) becomes that marked as Eq. ( 26), which describes nothing but the dynamics of the qubit-CNMR system.

The Hamiltonian of the STLR-qubit-NMR system
The Hamiltonian of the standing wave photons in the quarter-wavelength superconducting transmission line resonator (STLR) had been in previous work 29 , wherein the mcirowave current near the grounded point in the STLR reads Here, ω r is the frequency of the STLR.C r and L r are the total capacitance and inductance of the resonator, respectively.The quantized Hamiltonian of the rf-SQUID loop, without the NMR, can be expressed as Here, Ĥ0 = ω 0 |1 1| being Hamiltonian describing the rf-SQUID based qubit, and describing the STLR-qubit interaction, under the RWA approximation.Here, g rq = πMI p /(2L r √ ω r C r ) is the STLR-qubit coupling strength with M rq the mutual inductance between them.

Figure 1 .
Figure 1.Travelling microwave scattered by a qubit-NMR device for calibrating the physical parameters of the NMR embedded in a rf-SQUID loop.Here, the black part represents the Josephson junction, the red arrow perpendicular to the loop indicates the direction of the biased external flux Φ e .The black arrow perpendicular to NMR indicates the direction of another external magnetic field B 0 , which excites the mechanical vibration of the NMR (green) along the z direction.
gviven.Physically, by solving the Heisengberg equations: d b/dt = i[ b, Ĥ1 ] and d b † /dt = i[ b † , Ĥ1 ], one can get the time-evolution of the Bosonic operators; b(t) and b † (t).

Figure 6 .
Figure 6.Transmitted measurements of the travelling wave scattered by a STLR-qubit-NMR system.Here, the quarter-wavelength STLR is capacitively to the feedline and inductively coupled to the qubit embedded by the NMR.

Figure 12 .
Figure 12.The wave functions of the clockwise-(blue line) and anticlockwise (red line) supercurrent states, localized respectively in the left-and right wells of the symmetric potential.The parameters are the same as those used in Fig. 11