Dynamics of a qubit while simultaneously monitoring its relaxation and dephasing Q. Ficheux,1, 2 S. Jezouin,2 Z. Leghtas,3, 2, 4 and B. Huard1, 2, ∗ 1Universit´eLyon, ENS de Lyon, Universit´eClaude Bernard, CNRS, Laboratoire de Physique, F-69342 Lyon, France 2Laboratoire Pierre Aigrain, D´epartement de physique de l'ENS, Ecole´ normale sup´erieure, PSL Research University, Universit´eParis Diderot, Sorbonne Paris Cit´e,Sorbonne Universit´es, UPMC Univ. Paris 06, CNRS, 75005 Paris, France 3Centre Automatique et Syst`emes,Mines ParisTech, PSL Research University, 60 Boulevard Saint-Michel, 75272 Paris Cedex 6, France. 4QUANTIC team, INRIA de Paris, 2 Rue Simone Iff, 75012 Paris, France (Dated: April 26, 2018) Abstract Decoherence originates from the leakage of quantum information into external degrees of freedom. For a qubit the two main decoherence channels are relaxation and dephasing. Here, we report an experiment on a superconducting qubit where we retrieve part of the lost information in both of these channels. We demonstrate that raw averaging the corresponding measurement records provides a full quantum tomography of the qubit state where all three components of the effective spin-1/2 are simultaneously measured. From single realizations of the experiment, it is possible to infer the quantum trajectories followed by the qubit state conditioned on relaxation and/or dephasing channels. The incompatibility between these quantum measurements of the qubit leads to observable consequences in the statistics of quantum states. The high level of controllability of superconducting circuits enables us to explore many regimes from the Zeno effect to underdamped Rabi oscillations depending on the relative strengths of driving, dephasing and relaxation. Introduction fore extending the previously explored case of two in- Decoherence can be understood as the result of mea- compatible measurement outcomes [10, 12] to the case of surement of a system by its environment. For a qubit, three spin directions. By varying the drive amplitudes at the two main sources of decoherence are relaxation by the cavity and qubit transition frequencies, we are able to spontaneous emission and dephasing that can be mod- reach a variety of regimes corresponding to different con- eled by unmonitored readout of coupled quantum sys- figurations for Ω=Γ1 and Γd=Γ1, where Ω is the Rabi fre- tems (Fig. 1a). What becomes of the qubit state if, quency, Γ1 the fixed relaxation rate and Γd the dephasing instead of disregarding the information leaking to the rate. This work hence provides a textbook experimen- environment, we continuously monitor both decoherence tal demonstration of quantum measurement backaction channels? Owing to measurement backaction, the knowl- on a qubit with incompatible and simultaneous measure- edge of the measurement record then leads to a stochas- ments. tic quantum trajectory of the qubit state for each single Results realization of an experiment [1{3]. Recently, diffusive Description of the experiment. Two parallel de- quantum trajectories were observed following the contin- tection setups operate via spatially separated measure- uous homodyne or heterodyne measurements of either a ment lines (see Fig. 1b). The fluorescence heterodyne dephasing channel [4{9] or a relaxation channel [10, 11]. detection setup enables the measurement of both quadra- tures u(t) and v(t) of the spontaneously emitted field out Here we report an experiment in which we have si- of a 3D transmon qubit [13]. The complex amplitude of multaneously monitored the spontaneous emission of a arXiv:1711.01208v2 [quant-ph] 25 Apr 2018 the emitted field is on average proportional to the expec- superconducting qubit by heterodyne measurement (re- tation of the qubit lowering operator σ = (σ − iσ )=2 laxation channel) and the transmitted field through a − x y so that u and v are on average proportional to the ex- dispersively coupled cavity by homodyne measurement pectations of σ and σ . Here, σ = jgi hej + jei hgj, (dephasing channel). We demonstrate that the average x y x σ = i jgi hej − i jei hgj and σ = jei hej − jgi hgj are the outcomes of these two non-projective measurements are y z qubit Pauli operators, where jgi and jei are the ground the three coordinates x, y and z of the Bloch vector. It and excited states. For a single realization, the measure- is remarkable that a full quantum tomography can be ment outcomes read [10] obtained at any time by simply raw averaging measure- p ment outcomes of many realizations of a single experi- u(t)dt = ηf Γ1=2x(t)dt + dWu(t) ment despite the incompatibility of the three observables p ; (1) v(t)dt = ηf Γ1=2y(t)dt + dWv(t) that characterize a qubit state. For single realizations the −1 resulting quantum trajectories show signatures of the in- where Γ1 = (15 µs) is the qubit relaxation rate, compatibility between the measurement channels, there- ηf = 0:14 is the total fluorescence measurement efficiency, 2 excited state (the qubit and cavity are in the dispersive regime and χcq = 5:1 MHz as explained in Supplemen- tary Note 2). The phase of the measured quadrature in the homodyne measurement can then be chosen in such a way that [8, 17] p w(t)dt = 2ηdΓdz(t)dt + dWw(t): (2) Importantly, the measurement induced dephasing rate Γd can be tuned arbitrarily as it is proportional to the drive power at fd. Similarly to the notations above, ηd = 0:34 is the total dispersive measurement efficiency, z(t) = Tr(σzρt) is the last of the three Bloch coordinates (see Fig. 1a) and Ww is another independent Wiener process. Full tomography by direct averaging. As can be seen from Eqs. (1,2), taking a raw average of the FIG. 1: Measurement setup and quantum trajectory result- outcomes (u; v; w) on a large number of realizations of ing from its outputs. a Bloch vector representation of a the experiment directly leads to the Bloch coordinates qubit whose state is described by a density matrix ρt = (x; y; z) of the qubit. In Fig. 2, we show the direct aver- (1 + x(t)σx + y(t)σy + z(t)σz) =2. A quantum trajectory ρt is represented as a blue line. The qubit decoherence can be aging of the three outcomes in two configurations of the modeled as originating from a relaxation channel at a rate input drives: one in the regime of underdamped Rabi os- Γ1 and a dispersive measurement channel at a rate Γd. b A cillations (Fig. 2a) and another in the regime of strong superconducting qubit in a cavity is driven by two microwave dispersive measurement rate, the so-called Zeno regime signals at the weakly coupled input. The one at qubit fre- (Fig. 2b). The raw averaging of (u; v; w), once rescaled by quency fq = 5:353 GHz (orange) induces Rabi oscillations the prefactors in Eqs. (1,2), agrees well with the average of the qubit at frequency Ω. The one at cavity frequency evolution of the qubit, as predicted by the solution of the f = 7:761 GHz (purple) leads to a dispersive measurement d master equation (Eq. 3 below without the last stochastic of the qubit state along σz. A diplexer at the strongly cou- pled output port separates the outgoing signals depending on term). We thus demonstrate that performing a dispersive their frequency. The radiation at fq that is spontaneously measurement and a measurement of fluorescence reveals emitted by the qubit is processed by a Josephson Paramet- information on all three components of a spin-1/2. Such a ric Converter (JPC) [14, 15] so that a following heterodyne direct full tomography cannot be done by measuring two measurement reveals the two quadratures u(t) and v(t) of the records only [10, 12]. Note however that it is possible to fluorescence field [10, 16, 17]. The transmitted signal at fd perform an indirect tomography using a small number of is processed by a doubly pumped Josephson Parametric Am- plifier (JPA) [4, 18] with a pump phase such that a following records and maximum likelihood estimation [19]. A com- homodyne measurement reveals the quadrature w(t) of the parison between our technique and the usual technique field at fd. c Measurement records u (blue), v (red) and w using a qubit rotation followed by a projective measure- (yellow) as a function of time for one realization of the exper- ment is discussed in Supplementary Note 3. iment. These records feed the stochastic master equation (3), Experiments of Figs. 2a and 2b differ by the relative which leads to the trajectory in a rate of the dispersive readout Γd compared to the Rabi frequency Ω. For weak measurement rate Γd; Γ1 < Ω, the Rabi oscillations are underdamped, while they are x(t) = Tr(σxρt) and y(t) = Tr(σyρt) are the qubit Bloch overdamped when Γd Ω; Γ1 owing to the fact that the coordinates corresponding to the density matrix ρt (see Zeno effect prevents any unitary evolution such as Rabi Fig. 1a) and Wu and Wv are two independent stochas- oscillations. For a single realization of the experiment tic Wiener processes describing the measurement noise, though, the trajectory of the qubit state that one can which includes the zero point fluctuations, and such that infer from the measurement records u(t), v(t) and w(t) dW 2 = dt and dW is zero on average. Experimen- can strongly differ from this average behavior. tally, the measurement takes a non infinitesimal time Single quantum trajectories. In order to determine dt = 100 ns, which we chose smaller than the inverse this quantum trajectory, one can use the formalism of the measurement rates and compatible with the detection stochastic master equation [17].