Characterizing decoherence rates of a superconducting qubit by direct microwave scattering

We experimentally investigate a superconducting qubit coupled to the end of an open transmission line, in a regime where the qubit decay rates to the transmission line and to its own environment are comparable. We perform measurements of coherent and incoherent scattering, on- and off-resonant fluorescence, and time-resolved dynamics to determine the decay and decoherence rates of the qubit. In particular, these measurements let us discriminate between non-radiative decay and pure dephasing. We combine and contrast results across all methods and find consistent values for the extracted rates. The results show that the pure dephasing rate is one order of magnitude smaller than the non-radiative decay rate for our qubit. Our results indicate a pathway to benchmark decoherence rates of superconducting qubits in a resonator-free setting.


INTRODUCTION
Superconducting circuits are promising building blocks for implementing quantum computers [1][2][3] . In those devices, the key elements are superconducting artificial atoms made by Josephson junctions which induce a strong and engineerable nonlinearity. Such artificial atoms are also used in the field of superconducting waveguide quantum electrodynamics (waveguide QED) 4,5 , where they interact with a continuum of light modes in a 1D waveguide. In the past decade, many quantum effects from atomic physics and quantum optics have been demonstrated in waveguide QED, e.g., the Mollow triplet 6 , giant cross-Kerr effect 7 , and cooperative effects 5,8,9 . Other recent experiments have shown phenomena which are currently beyond the reach of atomic physics, such as ultra-strong 10,11 and superstrong coupling 12 between light and matter. Waveguide QED is also an enabling quantum technology. One of the key applications is to generate [13][14][15][16][17][18][19] and detect 4,20-26 single photons. It has been proposed to use waveguide QED to create bound states [27][28][29] and implement quantum computers [30][31][32] .
The performance of quantum computers and waveguide-QED devices is often limited by the coherence of the Josephson circuits. For example, the efficiency of producing and detecting single photons, the lifetime of bound states, and the fidelity of logical gates can all be improved by enhancing the coherence. In a waveguide-QED setup, decoherence can be due to decay into the waveguide, pure dephasing, and non-radiative decay into other modes. However, the rates for pure dephasing and nonradiative decay are typically not explored separately. An understanding of which one is dominant will give an insight into the decoherence mechanisms, and thus how the device performance can be improved.
In this work, we probe a superconducting transmon qubit coupled directly to the end of an open transmission line. In previous realizations [6][7][8][33][34][35] , the coupling rates were much larger than intrinsic decoherence mechanisms of the qubit, so the effects of non-radiative decay and pure dephasing were small and could not be well characterized. Here, we investigate a qubit whose radiative decay rate into the transmission line is larger than, yet comparable to, other decoherence mechanisms. This allows us to explore the pure dephasing rate Γ ϕ , the radiative decay rate Γ r from the capacitive coupling to the waveguide, and the nonradiative decay rate Γ n . The total relaxation and decoherence rates are given by Γ 1 = Γ r + Γ n and Γ 2 = Γ 1 /2 + Γ ϕ , respectively. We demonstrate different methods to extract the different rates and find consistent results. Similar to the results in circuit QED [36][37][38][39][40] , our methods can also enable the evaluation of the decoherence of qubits over a broad range of frequencies. They provide a pathway to investigate Josephson junctions or superconducting quantum interference devices (SQUIDs) without any resonator. In addition, we also consider it important to study Γ n and Γ ϕ separately. For instance, this could help to improve the Purcell enhancement factor, Γr Γnþ2Γϕ , in devices such as the one presented in ref. 9 . Moreover, the spontaneous-emission factor β, which is customarily quoted in other waveguide-QED platforms [41][42][43][44] is also related to Γ n , namely, β ¼ Γr ΓrþΓn in our case. The paper is structured as follows: first, we characterize the coherent scattering of the device and obtain the radiative decay rate and the decoherence rate of the qubit as a reference for later measurements. Then, we exploit the fluorescence of the qubit under coherent excitation to find the non-radiative decay rate and the pure dephasing rate. The resonance fluorescence spectrum at strong driving develops into the Mollow triplet 45 , which has been widely used to probe quantum properties in systems based on superconducting qubits such as coherence 6,8 and vacuum squeezing 46 . The resonance spectrum is symmetric around the central peak. However, if pure dephasing exists, the off-resonant spectrum becomes asymmetric, something which has been studied experimentally in quantum dots 47,48 . We take advantage of this fact to extract the pure dephasing rate. We also measure the non-radiative decay rate under a continuous coherent drive, where coherently and incoherently scattered photons provide information about the different decay channels. Afterward, we apply a pulse to the qubit to both obtain the decay rates and find the stability of the qubit frequency and coherence as a function of time. In contrast to other methods, we use the phase information of emitted photons from the qubit to investigate the qubitfrequency stability in superconducting waveguide QED. Finally, we summarize the measured results and compare the advantages and disadvantages of the different methods.
In this paper, we have investigated waveguide QED with a qubit weakly coupled to a transmission line which is an unexplored regime. In this regime, it becomes possible to quantify the nonradiative decay rate of the qubit. We use several different methods to extract the different decay rates, provide a thorough comparison of the different methods, and find consistent results. Among these methods, those based on the off-resonant Mollow triplet and the power loss analysis are original. Moreover, we show that time-domain analysis of decay rates can also be done in waveguide QED.

Device characterization
The device used in our experiment (see Fig. 1a) is a magnetic-fluxtunable Xmon-type transmon qubit 49 , capacitively coupled to the open end of a one-dimensional transmission line with characteristic impedance Z 0 ≃ 50 Ω. The circuit is equivalent to an atom in front of a mirror in 1D space. The device is fabricated from aluminum on a silicon substrate using the same fabrication recipe as in ref. 38 . We denote 0 j i, 1 j i and 2 j i as the ground state, the first and the second excited states of the qubit, respectively. The 0 is the charging energy, e is the elementary charge, C ∑ is the total capacitance of the qubit, and E J (Φ) is the Josephson energy. The Josephson energy can be tuned from its maximum value E J;max by an external magnetic flux Φ using a coil: Figure 1a illustrates the simplified experimental setup for measuring the reflection coefficient of a probe signal coming from a vector network analyzer (VNA), after interacting with the qubit. The probe signal at frequency ω pr and a pump at frequency ω p are combined and attenuated before being fed into the transmission line. Then, the VNA receives the reflected signal to determine the complex reflection coefficient. Figure 1b shows the magnitude of the reflection coefficient, |r|, for a weak probe (with an intensity Ω pr < Γ 2 ) as a function of the external flux Φ. Afterward, all the other measurements are taken at the flux-sweet spot. We use two-tone spectroscopy to determine the anharmonicity of the qubit, α = (ω 12 − ω 01 )/2π, where ω 12 is the frequency of the 1 j i À 2 j i transition. Specifically, we apply a strong pump (with an intensity Ω p ≫ Γ 2 ) at ω 01 to saturate the 0 j i À 1 j i transition and measure the reflection coefficient as a function of probe frequency. The result is shown in Fig. 1c: a dip appears in the reflection at ω pr = ω 12 due to the photon scattering from the 1 j i À 2 j i transition. At a higher pump power, the Autler − Townes splitting is observed 50 . From Fig. 1b, c, we find E C /h ≈ α = 252 MHz, and then, E J;max =h ¼ 16:56 GHz by fitting the data in Fig. 1b to the equation between the qubit frequency and the external flux mentioned previously. In order to obtain the radiative decay and decoherence rates, we perform single-tone spectroscopy with a weak probe (Ω pr ≪ Γ 2 ). Figure 1d, e are the corresponding magnitude and phase response of r, where we obtain Γ r /2π = 227 kHz and Γ 2 /2π = 141 kHz by using the circle fit technique from ref. 51 . The corresponding fitting equation is where the phase ϕ ≈ 0.150 ± 0.004 quantifies the impedance mismatches in the line (for a perfectly matched line, ϕ = 0).

Atomic fluorescence
Even though measuring the reflection coefficient can give the decoherence and radiative decay rates, it cannot distinguish between pure dephasing and non-radiative decay. In order to distinguish them, we study the atomic fluorescence for different pump intensities and frequencies. For this measurement, the VNA is turned off, the pump is used to drive the system, another 50 dB of attenuation is added between the directional coupler and the pump, and the output signal is sampled by a digitizer (see Fig. 1a). When the qubit is pumped, its state evolves at a Rabi frequency Ω.
With a Rabi frequency much larger than the natural linewidth of the qubit (Ω ≫ Γ 2 ), the energy levels of the qubit become dressed, leading to three distinct spectral components known as the Mollow triplet 45 . In particular, the spectrum contains the elastic "Rayleigh" line in the middle in which the scattered wave has the b Single-tone spectroscopy. The magnitude of the reflection coefficient r is measured as a function of the external flux Φ and probe frequency ω pr . The red dashed curve is a fit for the qubit frequency ω 01 . c Two-tone spectroscopy. A strong pump is applied to the 0 j i À 1 j i transition and we sweep the frequency of the weak probe on the 1 j i À 2 j i transition to obtain the corresponding reflection coefficient. d, e show the magnitude and phase response of the qubit at Φ = 0 under weak probing. Red lines are the corresponding fits using the circle fit technique from ref. 51 .
same frequency as the incident wave, with two inelastic sidebands positioned symmetrically on both sides of the central peak.
To measure the power spectrum for the atomic fluorescence, we use a digitizer to sample the amplified signal as V(t) in the time domain. The sample rate and the number of samples are from 2 MHz to 6 MHz and from 500 to 1000, respectively, for the different measured spectra in Figs. 2 and 3. After that, we calculate the power spectrum S meas: i ðωÞ according to the Welch method 52 , after normalization to the system gain G. In order to increase the corresponding signal-to-noise ratio (SNR), a number of averages are required to obtain hS meas: i ðωÞi. For different results shown below, the number of averages varies from 2 × 10 6 to 10 8 .

On-resonant Mollow triplet
When the qubit is driven by a resonant pump (Δ = ω p − ω 01 = 0), as shown in Fig. 2a, the splitting of the spectrum between the sidebands and the central peak increases as a function of the pump power P p . The splitting equals the Rabi frequency and obeys Ω ¼ 2 ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi AΓ r P p =ð_ω 01 Þ p . By fitting the extracted Rabi splitting |Ω| in Fig. 2b, and using Γ r from the previous measurement in the section "Device characterization", we extract a total attenuation A = −145 dB of which about −125 dB attenuation is from attenuators and directional couplers, −7 dB from an Eccosorb filter and the rest is due to cable loss. This allows us to renormalize all applied powers to either the power at the qubit, or the corresponding Rabi frequency. The total gain in the output line of the measurement setup can be calibrated by tuning the qubit away and measuring the power at the output port at room temperature. This results in a total gain G = 115 dB, of which approximately 44 dB gain comes from a high electron mobility transistor (HEMT) amplifier, and the rest is from the room temperature amplifiers and the pre-amplifiers of the digitizer.
The Rabi rate can be made much larger than all the decay rates of the qubit (Ω ≫ Γ 1 , Γ 2 ). Consequently, the overlap in the frequency domain between the sideband emission and the central peak becomes negligible. In Fig. 2c, we use an input Fig. 2 Power spectral density (PSD) of the on-resonant Mollow triplet from the atomic fluorescence. a Resonant fluorescence emission spectrum as a function of the pump power and detuning of the detected radiation, δω 01 = ω − ω 01 . P p is the power from the RF source while P q is the corresponding power on the qubit. Inset: a schematic of the triplet transitions in the dressed-state picture, where the qubit energy levels split by Ω due to strong driving, creating three transitions with frequencies ω 01 − Ω, ω 01 and ω 01 + Ω. b Rabi splitting Ω (dots) vs. drive amplitude, extracted from (a). The black line is the linear fit to obtain the attenuation in the input line which is A = −145 dB. c Power spectral density at −116 dBm power at the sample. The black lines are individual fits to the linewidths of the three peaks, yielding Γ 2 /2π = 141 ± 2 kHz and Γ 1 /2π = 276 ± 5 kHz. The arrows correspond to the transitions in the inset of (a).
Power spectral density of the off-resonant Mollow triplet from the atomic fluorescence. a Off-resonant fluorescence emission spectrum as a function of the pump power at the qubit. The frequency of the pump field is detuned by Δ = ω p − ω 01 = −2π*790 kHz from the qubit frequency. The Rabi splitting is increased with the input power. b Off-resonant fluorescence emission spectrum as a function of the frequency of the pump at P q = −133 dBm. c Off-resonant PSD at Δ/2π = −790 kHz (green dots) with P q = −133 dBm. The solid curve is a fit to Γ 1 /2π = 275 ± 6 kHz and Γ ϕ /2π = 3 ± 3 kHz. d Off-resonant PSD in a second cooldown at Δ/2π = 820 kHz (brown dots), Δ/2π = −830 kHz (green dots). A value of Γ ϕ /2π = 7 kHz gives a good fit to both traces (black).
power to the qubit P q ≈ −116 dBm, equivalent to Ω/2π ≈ 9 MHz. The incoherent part of the corresponding power spectral density (PSD) is given by Γs (see more details in Supplementary section A), where the half width at half maximum of the central peak and the sidebands are Γ 2 and Γ s = (Γ 1 + Γ 2 )/2, respectively. The solid curves in Fig. 2c are fits to 2πS i (ω) which is the PSD expressed in linear frequency. We obtain Γ 2 /2π = 141 ± 2 kHz for the central peak, Γ s,red /2π = 210 ± 3 kHz and Γ s,blue /2π = 206 ± 4 kHz for sidebands. By taking the average of Γ s,red and Γ s,blue , we obtain Γ 1 /2π = 275 ± 7 kHz. We note that the extracted Γ 2 /2π value is fully consistent with the result from the reflection-coefficient measurement in the section "Device characterization". From that measurement, we also know Γ r /2π = 227 ± 1 kHz. Thus, we can now extract both the nonradiative decay rate, Γ n /2π = 48 ± 7 kHz, and the pure dephasing rate, Γ ϕ /2π = 3 ± 4 kHz. By integrating the PSD of each peak in the Mollow triplet we can compare their relative weights. After normalization with hω 01 Γ r , the results are about 0.254, 0.116 and 0.124 for the middle peak, the red and blue sidebands, respectively. According to Eq.
(2), we would expect these numbers to be 0.250, 0.125 and 0.125, respectively, for a fully saturated qubit.

Off-resonant Mollow triplet
We also study the off-resonant Mollow triplet at a variety of pump powers and frequency detunings between the pump and the qubit. In Fig. 3a, the pump power at the qubit is swept from −150 dBm to −130 dBm at detuning Δ/2π = 790 kHz. We find that the PSD is weaker than the on-resonant case, implying that the qubit is less excited. In Fig. 3b, as we sweep the frequency detuning between the pump and the qubit, the spectrum is nearly symmetric over 1 MHz bandwidth, therefore, the extracted pure dephasing rates by the off-resonant fluorescence are insensitive to the frequency detuning Δ. We can either choose a large Δ which will lead to a small excitation of the qubit, or a small Δ which results in an unresolved spectrum between the central peak and sidebands.
Compared to the on-resonant case, the off-resonant Mollow triplet carries additional information in its sideband asymmetry and the approximation used in Eq. (2) is no longer valid. Therefore, the full expression for the PSD must be used, shown in Eq. (16) in Supplementary section A which is an extension of ref. 53 . The fit of the data in Fig. 3c to Eq. (16) yields Γ ϕ /2π = 3 ± 3 kHz (pink solid line). The symmetry of the sidebands around the central peak is due to the relatively small pure dephasing rate. In Fig. 3d, the sample was measured in an earlier cooldown. There, we observed a larger asymmetry for both positive (green dots) and negative detunings (brown dots). In the case of positive detuning, the red sideband is closer to the qubit original frequency than the blue sideband, whereas the blue sideband is closer when the sign of the detuning is changed. We fit the two sets of data simultaneously to obtain Γ ϕ =2π ¼ 7 ± 2 ð ÞkHz, which is slightly larger than the second cooldown. This is likely due to that we used only two isolators in the first cooldown, and four isolators in the second cooldown, leading to less thermal photons from the transmission line in the second case.
The mechanism by which the pure dephasing gives rise to an asymmetry in the Mollow triplet can be explained as follows (for details, see Supplementary section B). Relaxation from the qubit will cause transitions between dressed states jn; ± i (see Fig. 4a) that contain different numbers n of drive photons. As shown in Fig. 4b, these transitions will either be between or within the + and − subspaces. In equilibrium, if the pure dephasing rate is zero, the probabilities P ± for the system to be in these subspaces are given by the detailed-balance condition i.e., the number of emitted photons causing transitions from + to − (the blue sideband) must equal the number of emitted photons causing transitions from − to + (the red sideband). However, the interaction causing pure dephasing has a non-zero matrix element for transitions between n; þ j i and n; À j i, which leads to a modified detailed-balance condition: As Γ ϕ increases, this will push the occupation probabilities towards P + = P − . For off-resonant driving Γ +− ≠ Γ −+ and thus the number of emitted photons in the two sidebands becomes different: Γ + − P + ≠ Γ −+ P − . The larger number of photons will be emitted at the frequency corresponding to the larger of the two transition rates Γ +− and Γ −+ ; from transition-matrix elements, this can be seen to be the frequency that is closer to the qubit frequency. From Fig. 3c, we have Γ 1 /2π = 275 ± 6 kHz and Γ 2 /2π = 140 ± 3 kHz. Again, from the measured radiative decay rate, by subtracting Γ r from Γ 1 , we obtain Γ n /2π = 48 ± 6kHz. Based on the results in this section, the on/off-resonant Mollow spectra allow us to extract the pure dephasing rate and non-radiative rate of a qubit. Specifically, for our qubit in this environment, we find that the non-radiative decay rate is one order of magnitude larger than the pure dephasing rate. Compared to the on-resonant Mollow triplet, the off-resonant Mollow triplet allows us to characterize the qubit decay rates at a lower pump power.
Photon scattering by the qubit To verify the extracted decay rates above, we can also measure the power scattered by the qubit and the dissipated power due to the non-radiative decay channel directly. We normalize all the powers by the single-photon energy hω 01 . The pump is on resonance with the qubit. The output power then consists of a coherent part and an incoherent part, P out = P coh + P incoh , where P coh ¼ Ω 2 4Γr ð1 À Γ1Γr Ω 2 þΓ2Γ1 Þ 2 and P incoh ¼ Γr

(see Supplementary section D).
For our qubit, the pure dephasing rate was verified to be around 3 kHz, i.e., much less than other rates and therefore negligible, so, the expression for the incoherent power can be further simplified to P incoh ' 2Γ r ρ 2 11 , where ρ 11 is the population of the first-excited state of the qubit. In this case, the expression for the dissipated power due to the non-radiative decay is then Experimentally, we use about 4.2 × 10 9 averages to measure all the powers. We denote the measured voltage V and the pump power P in . The subscripts "off" and "on" used in the following contexts mean that the qubit is off/on resonance with the pump, respectively. When the qubit is tuned away, it is off-resonant with the pump; we will have P in ¼ hVi 2 off because of the coherence of the pump. Besides the pump power, the system noise will also make a contribution to the total measured power, P meas: in . Therefore, we have P meas: where P noise is the noise power added by the amplification chain over the measurement bandwith. When instead the qubit is on resonance with the pump, the total measured output power P meas: out ¼ hV 2 i on ¼ P out þ P noise , where P coh ¼ hVi 2 on (Fig. 5, black circles). Therefore, P loss (blue crosses) is obtained by taking P loss ¼ P meas: in À P meas: out ¼ P in À P out with P incoh = P in − P coh − P loss (green circles). Figure 5 shows all the types of measured power as a function of the Rabi frequency. There, we find three interesting regions: (i) At high input power, when Ω > ð1 þ 1 ffiffi 2 p ÞΓ r % 2π Ã 391 kHz, to the right of the dashed line in Fig. 5, the qubit starts to be saturated. The outgoing field is then mainly coherent from the pump itself. By increasing the input power further, the qubit is completely saturated, leading to P in % P coh % Ω 2 4Γr , P incoh % Γr 2 and P loss % Γn 2 . In this case, almost all the incoming photons are reflected by the mirror.
(ii) In the low-power region Ω <

ffiffiffiffiffiffiffiffiffi ffi
Γ1Γ2Γn ΓrÀΓn q % 2π Ã 103 kHz derived from P incoh = P loss , to the left of the dash-dotted line, the scattering process is dominated by the interaction between the qubit and the incoming photons. The incoherent scattering is proportional to ρ 2 11 whereas the power loss depends linearly on the excitation probability. Therefore, the incoherent power can be less than the power loss when ρ 11 < Γn 2Γr % 0:11. Besides the incoherent photons, there is a small coherent scattering by the qubit. Compared to the loss, the coherent power is smaller if the non-radiative decay < Ω < ð1 þ 1 ffiffi 2 p ÞΓ r , both the mirror and the qubit make substantial contributions to the scattering process. The photons reflected by the mirror interfere destructively with those scattered by the qubit, resulting in a suppression of the coherent part of the output field. In particular, the dip around Ω/2π ≈ 160 kHz in the coherent power appears due to the fully destructive interference. In addition, the qubit excitation is not small anymore and the incoherent power is larger than the loss because Γ r > Γ n . We note that this region can be non-existent when either the non-radiative decay or the pure dephasing is sufficiently large.
The results shown here agree well with the values from other sections in the paper, which implies that it is possible to take the Γ r value from the reflection coefficient measurement as a reference for the Mollow triplet in order to separate the non-radiative decay rate and the pure dephasing rate.
Time-resolved dynamics All measurements described previously span time ranges from several hours to tens of hours. It is noteworthy that qubit decay rates extracted by different methods agree relatively well. However, the long duration means that any fluctuations of the decay rates are averaged out. Recently, several groups have characterized such fluctuations in circuit QED, using Rabi pulses, Ramsey interference measurements, and dispersive qubit readout [36][37][38][39]54 . To probe the decoherence of the qubit with a temporal Fig. 5 Normalized powers as a function of Rabi frequency. The input power, P in (red), representing the input photon flux at the qubit, is measured when the qubit is tuned away by the external flux. When the qubit is on resonance with the input signal, we have the coherent power, P coh (black circles), consisting of photons reflected from either the qubit or the end of the transmission line. The qubit can also scatter photons incoherently, P incoh (green circles), due to the decoherence of the qubit. Moreover, the excited qubit has some probability to release a photon to the environment resulting in the power loss, P loss (blue crosses). The solid curves are fits to different types of the scattered powers. The dotted and dashdotted lines separate the qubit response into three interesting regions (see more details in the text). resolution of 7 minutes, we prepare the qubit in a superposition of the ground and the first-excited state and monitor its spontaneous emission into the waveguide by recording both quadratures of the output field with a digitizer as a complex trace in the time domain, i.e., V(t) = I(t) + iQ(t). We measure for 4.27 × 10 5 s (approximately 119 hours) with 975 repetitions, and each trace has 2.30 × 10 6 averages. The averaged trace, 〈V(t)〉, is proportional to the qubit emission. In order to generate a pulse at the qubit frequency, a signal generator at the frequency 250 MHz higher than the qubit frequency is mixed with a pulse with a 250 MHz carrier frequency. This pulse is from an arbitrary waveform generator (AWG) where the bandwidth of the AWG is 1.25 giga samples per second.
In detail, after the 50 ns long π/2-pulse, the qubit superposition state evolves in time τ as 1 ffiffi 2 p ð 0 j i þ e ÀΓ2τÀiδω01τ 1 j iÞ with δω 01 = ω 01 − ω pulse . The emitted field carries information about the qubit operator hσ À i ¼ e ÀΓ2τ e Àiδω01τ , where the amplitude response and the phase response show the decoherence and the qubit frequency shift with τ, respectively. Figure 6a, b shows the magnitude and phase response of a single trace where the decay of the magnitude is fitted to an exponential curve and the phase of the photons emitted from the qubit grows linearly with time due to the free evolution of the qubit, where the slope determines δω 01 /2π = 125 kHz. Figures 6c, d are histograms of Γ 2 and δω 01 for all the repetitions. Both histograms can be fitted to a Gaussian with parameters shown in the figures. In comparison with the decoherence rates extracted from other measurements, we find that the standard deviation here is larger than the previously measured error bar. This shows that the dynamics of the qubit on a short time differs slightly from that over a long measurement time. By taking the average of all the traces in Fig. 6d, we fit to an exponential decay and get an averaged Γ 2 /2π = 145 ± 1 kHz.
To study Γ 1 , we instead send a π-pulse to fully flip the qubit, and then measure the emission from the qubit. The corresponding output power, PðτÞ ¼ ð_ω 01 Γ r =2Þð1 þ hσ z iÞe ÀΓ1τ 55 allows us to determine Γ 1 . The trace is measured with 1.92 × 10 9 averages, shown in Fig. 6e. A fit to an exponential decay with Γ 1 /2π = 273 ± 11 kHz agrees well with the data. Combining these numbers with Γ r from section "Device characterization", we can also calculate Γ n and Γ ϕ from these measurements. The resulting values can are summarized in Table 1. The decay of the qubit state by averaging all the measured traces in (c) to extract the decoherence rate. f A π-pulse is applied to flip the qubit to the excited state with 1.92 × 10 9 averages, where P(τ) is the power emitted by the qubit at time τ after the pulse. By fitting the emitted power (blue circles) to an exponential decay, we extract Γ 1 /2π = 273 ± 11 kHz. Except for the histograms, the error bars are for 95% confidence.

DISCUSSION
We have shown several methods to determine different decay rates of a qubit placed in front of a mirror. In principle, these methods can also be used when the qubit couples to a transmission line without a mirror, except for the scattering method, where the corresponding measurement taken on both the input and output ports is required. The measured rates are consistent between methods within the error bars of two standard deviations corresponding to 95% confidence. The results are summarized in Table 1. The reflection measurement is the baseline to provide the value of Γ r to extract the non-radiative decay rate of the qubit for measurements, except for the scattering measurement. These different methods have advantages and disadvantages that we summarize below: (i) The fastest way to obtain the non-radiative decay rate is to send a strong pump on resonance with the qubit so that the central peak and the sidebands of the Mollow triplet do not overlap. The drawback is that the pump power needed here is much stronger than for the other methods and that may change the rates slightly. (ii) In the second method, we measure the off-resonant Mollow triplet by detuning the pump frequency slightly from the qubit frequency. The sidebands will be asymmetric around the central peak if the pure dephasing rate is non-negligible. In this case, only weak probe power is required. However, the corresponding measurement time is increased by almost a factor of three. (iii) The most accurate way to measure the non-radiative decay rate is to measure the difference between the input and output power, labeled as Scattering in Table 1. Using this method, we can obtain not only the power loss but also the coherent and incoherent power scattered by the qubit. However, the measurement time is much longer. In addition, the attenuation between the sample and the input line, as well as the gain between the detector and the sample have to be calibrated at the beginning, in order to get the absolute power values from the qubit. To simplify the measurement, as we described in "Photon scattering by the qubit", we can see that when the pump saturates the qubit, we get P incoh % Γr 2 and P loss % Γn 2 . Then, the ratio of the non-radiative decay rate to the radiative decay rate can be obtained from the ratio of the lost power to the incoherent power. Knowing the value of Γ r from the reflection measurement, we can obtain the non-radiative decay rate. Therefore, in principle, we do not need to sweep the pump power as was done in Fig. 5. This simple way is labeled as SinglePoint in Table 1. (iv) Finally, pulses can be applied to excite the qubit. Afterward, the exponential decay of the emission and the emitted power trace from the qubit can be recorded to extract the total relaxation rate and the decoherence rate with a much larger measurement bandwidth. The distortion on the scattered photons due to the non-flat frequency response will affect the extracted values of the decay rates. This may be the reason why decay rates from this method are slightly different from those measured by other methods. However, the advantage of this method is that it allows us to study the short-time dynamics of the qubit.
The measurement time for these methods is from 2 to 63 h. The coherent measurement is related to the first moment (amplitude) whereas other methods are related to the second moment (power). In order to estimate the system noise N, we measured the background PSD by turning the drive off (not shown) and comparing the result with the measurement of Fig. 2c. We found N ≈ 49 photons. However, we expect that using a quantum-limited Josephson traveling-wave amplifier 56 would reduce the system noise to about two photons. The corresponding SNR is improved by a factor of 5. Therefore, this would result in a reduction of the measurement time by factors of 25 and 625, respectively, for the coherent measurement and the other methods 57 .
From the measured result, our qubit is T 1 -limited, i.e., the radiative decay dominates the interaction. However, the nonradiative decay rate is one order of magnitude larger than the pure dephasing rate. The corresponding spontaneous-emission factor is β ≈ 85%, which is typically close to 100% when we engineer the radiative decay much larger than the non-radiative decay. Therefore, reducing the non-radiative decay rate will be the next step to improve the intrinsic coherence of our qubit. In addition, it is worthwhile to investigate why the non-radiative decay rate of our qubit is one order of magnitude larger than the qubit coupled to a resonator which was fabricated on the same wafer 38 in the future.
Our methods allow us to analyze all the decay channels in detail. This will be useful to study and engineer decay channels of the qubit, which is the critical element in superconducting circuits. For example, engineering the decay channels can improve the quantum efficiency of generating single photons, which is set by Γ r /2Γ 2 . Also, the fidelity of detecting a single photon can be increased by extending the qubit coherence time. More importantly, compared to circuit QED where a resonator dispersively couples to a qubit, our study provides a straightforward way to investigate superconducting qubits, which are crucial elements in superconducting quantum computers.

DATA AVAILABILITY
The data that supports the findings of this study are available from the corresponding authors upon reasonable request.  (1) 20  23 The first method is a measurement of the reflection coefficient under a weak probe. On-res.MT and Off-res.MT represent the on/off-resonant Mollow triplet from atomic fluorescence, respectively. The energy loss is estimated by calibrating the input power and measuring both the coherent and incoherent power scattered by the atom. Finally, timeresolved measurements of the decay from both a superposition state and the first-excited state were used. BW and T are the measurement bandwidth and the total measurement time, respectively, for each method.