Converting quasiclassical states into arbitrary Fock state superpositions in a superconducting cavity
1 1 1 1 1 2, 3 1 1, 4, † 1, ‡ W. Wang, L. Hu, Y. Xu, K. Liu, Y. Ma, Shi-Biao Zheng, ∗ R. Vijay, Y. P. Song, L.-M. Duan, and L. Sun 1Center for Quantum Information, Institute for Interdisciplinary Information Sciences, Tsinghua University, Beijing 100084, China 2Department of Physics, Fuzhou University, Fuzhou 350116, China 3Department of Condensed Matter Physics and Materials Science, Tata Institute of Fundamental Research, Homi Bhabha Road, Mumbai 400005, India 4Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, USA We propose and experimentally demonstrate a new method to generate arbitrary Fock-state superpositions in a superconducting quantum circuit, where a qubit is dispersively coupled to a microwave cavity mode without the need of fine-frequency tuning. Here the qubit is used to conditionally modulate the probability amplitudes of the Fock state components of a coherent state to those of the desired superposition state, instead of pumping photons one by one into the cavity as in previous schemes. Our method does not require the adjustment of the qubit frequency during the cavity state preparation, and is more robust to noise and accumulation of experimental errors compared to previous ones. Using the method, we experimentally generate high-fidelity phase eigenstates under various Hilbert-space dimensions and squeezed states, which are useful for quantum walk and high- precision measurements.
The creation and manipulation of complex superposition ber of the largest Fock state component in the superposition states of a harmonic oscillator (e.g., a light field) is an impor- by 1. After N steps, the oscillator evolves from the ground tant subject in quantum optics [1], and has various applica- state to any superposition of the first N + 1 number states. tions in quantum technologies. In contrast to nonlinear quan- Such a technique was first demonstrated on the harmonic mo- tum systems (e.g., atoms), a harmonic oscillator possesses an tion of a trapped ion trap [16], and then implemented in a infinite-dimensional Hilbert space where the energy is linear superconducting resonator [17, 18] using a qubit with fast fre- in the quantum number. Due to this feature, a harmonic os- quency tunability to coherently pump photons from external cillator has extremely rich nonclassical effects and serves as microwave drive pulses into the resonator. More recently, an ideal candidate for investigation of fundamental quantum it was shown that 2N selective number-dependent arbitrary theory. In addition to fundamental interest, it plays impor- phase gates, together with 2N + 1 displacement operations, tant roles in many quantum computation protocols. On one can be used to construct such superposition states [19]. Fol- hand, it acts as a data bus to link different qubits, as in the lowing this proposal, a one-photon Fock state was created in ones with trapped ions [2], cavity quantum electrodynamics a circuit QED system [20]. With a similar setup, the Gradi- (QED) [3], and circuit QED [4–6]. It can also be used as a ent Ascent Pulse Engineering method [21, 22] was recently quantum memory to store quantum bits (qubits), offering the used for creating the six-photon state and four-component cat possibility to construct a quantum processor with fewer re- states [23]. In the context of cavity QED, Fock states and cat sources [7] and to realize quantum error correction [8–10]. states of a cavity were conditionally generated with Rydberg On the other hand, harmonic oscillators are the basic building atoms dispersively coupled to the cavity mode [24]. block for implementation of quantum computation based on Here we propose and experimentally demonstrate a dis- continuous variables [11, 12]. In addition, the nonclassical ef- tinct method to synthesize any superposition of Fock states fects of harmonic oscillators have applications in other quan- for the field stored in a cavity that is dispersively coupled to tum technologies, e.g., when a quantum state exhibits squeez- a qubit. Our method relies on modulating the probability am- ing, i.e., the reduction of the fluctuation in one certain quadra- plitudes associated with the superposed Fock states in an ini- ture below the vacuum level, it can be used in high-precision tial coherent state, a superposition of photon number states measurement [13]. Many of these applications rely on the with a Poissonian distribution, to those of the desired target ability to generate an arbitrary superposition state of a har- state. Coherent states are quasiclassical states characterized arXiv:1703.03316v1 [quant-ph] 9 Mar 2017 monic oscillator, however, its linear energy spectrum prevents by a complex amplitude [24], and can be generated by using one from doing so by using classical drives only; the appli- a classical drive. The photon-number-dependent shift in the cation of a classical drive produces a displacement in phase qubit frequency, arising from the dispersive qubit-cavity cou- space without enhancing the nonclassical behaviors. pling, enables one to use carefully tailored microwave drive Previous methods employed a step-by-step algorithm [14, signals to individually control the qubit transition amplitude 15], where in each step the quantum state of a two-level sys- associated with each Fock state. The subsequent measure- tem is manipulated in a controllable manner and then trans- ment of the qubit transition conditionally projects the cavity ferred to the harmonic oscillator, increasing the quantum num- to a superposition of Fock states, with the probability ampli- tude of each state component being proportional to the associ- ated qubit transition amplitude, controllable by the amplitude and phase of the corresponding drive signal. Unlike the previ- ous methods, it is unnecessary to individually pump photons ∗Electronic address: [email protected] †Electronic address: [email protected] into the cavity one by one, and the modulating drive signals ‡Electronic address: [email protected] can be simultaneously applied and incorporated into a single 2 pulse. Furthermore, during the preparation process the qubit- 훿0, Ω0, 휙0 … cavity detuning needs not to be adjusted, which is important N+1 for quantum state engineering in circuit QED systems with 훿푁, Ω푁, 휙푁 three-dimensional cavities [20, 25–27], where it is difficult to Qubit adjust the qubit frequency. We implement the experiment in a superconducting circuit, where a three-dimensional cavity is Cavity coupled to a transmon qubit. We analyze the produced states 퐷(훼) by measuring the Wigner functions, showing good agreement Initialization Post-selection Tomography with the desired ones. to |g , 0 on |푒 The system under investigation consists of a qubit disper- sively coupled to a cavity and driven by a classical microwave FIG. 1: (Color online) Experimental sequence. The experiment containing N + 1 frequency components. In the interaction starts with the application of a microwave pulse to the cavity, cor- picture, the Hamiltonian for the total system is responding to a phase-space displacement operation D(α), which transforms the cavity from the ground state 0 to the coherent state α . Then the qubit, initially in the ground| statei g , is driven by a N | i | i † i(δnt+φn) classical microwave pulse comprising N + 1 frequency components, HI = χqca a e e + ∑ Ωne e g + h.c. , (1) − | ih | "n=0 | ih | # with the nth component resonant with the qubit transition conditional on the cavity’s photonic field being in the Fock state with n photons. where e and g denote the excited and ground states of the This corresponds to simultaneous implementation of N + 1 condi- | †i | i tional qubit rotations on the Bloch sphere: N R ( ). Subse- qubit, a and a are the creation and annihilation operators of ∏n=0 en n βn quent measurement of the qubit in the excited state e| collapses the the photonic field in the cavity, χqc is the qubit’s frequency | i shift induced by per photon in the cavity due to the dispersive cavity to a superposition of N + 1 Fock states, whose amplitudes are qubit-cavity coupling, δ is the detuning between the qubit controllable by the parameters of the corresponding microwave com- n ponents. A final Wigner tomography is used to examine the produced and the nth driving component with a Rabi frequency Ωn and states. a phase φn, and h.c. denotes the Hermitian conjugate. The qubit is initially in the ground state g , and the cavity in a | i coherent state the maximum photon number in the desired target state. The detection of the qubit in the state e projects the cavity to the ∞ | i α = ∑ cn n , (2) superposition state | i n=0 | i N i(φn+nχqcτ) 2 n ψc = N cne sinβn n , (5) where cn = exp( α /2)α /√n! is the probability ampli- ∑ | i n=0 | i tude for having −n |photons.| The photon-number-dependent shift of the qubit transition frequency enables individually 1/2 N 2 2 − where N = ∑ cn sin βn is a normalization factor. modulating the probability amplitudes of the superposed Fock n=0 | | states in α with a carefully tailored microwave drive ad- The final cavityh state ψc is equivalenti to the desired state | i | i dressing the qubit. Our aim is to convert such a quasiclassical ψd when the parameters βn and φn are appropriately chosen | i i(nχ τ+φ ) state to a target quantum state so that dn = N cn sinβne qc n . To maximize the success probability P = 1/N 2, we numerically optimize the param- N eters α and β , with the optimal values depending upon the ψ = d n , (3) n d ∑ n desired| state.| | i n=0 | i Our experiment is implemented with a superconducting cir- with the desired complex amplitudes dn for the corresponding cuit, where one transmon qubit is dispersively coupled to two Fock state components. three-dimensional cavities [25–29]. The qubit, with a transi- Under the condition δ = nχ and Ω χ , the nth com- tion frequency ωq/2π = 5.345 GHz, has an energy relaxation n qc n qc ponent resonantly drives the qubit transition when the cavity time T1 = 13.0 µs and a pure dephasing time Tφ = 13.8 µs. is in the Fock state n , but does not affect the qubit state if the Cavity 1, with a frequency ωs/2π = 8.230 GHz and a photon | i cavity is in any other Fock state due to large detunings. The lifetime τs = 80 µs, is used for storage of the photonic state. dynamics governed by HI enables individually controlling the Cavity 2, with a frequency ωr/2π = 7.290 GHz and a photon transition in each subspace g,n , e,n (the notation de- lifetime τr = 44 ns, is connected to a Josephson parametric notes qubit,cavity ); there is{ no| interactioni | i} between different amplifier [30–33], allowing for a high-fidelity (nearly-unity) subspaces.| After ani interaction time τ the qubit-cavity system and quantum non-demolition single-shot readout of the qubit evolves to state. The average readout photon number in the readout cav- ity is about five and the readout time is 240 ns. For simplic- ∞ i(φn+nχqcτ) ity, the term “cavity” refers to the cavity 1. The dispersive ψqc(τ) = ∑ cn cosβn g,n ie sinβn e,n , n=0 | i − | i frequency shift is χqc/2π = 1.44 MHz. The experimental h i(4) apparatus and readout properties− can be found in the Supple- τ where βn = 0 Ωndt is the pulse area associated with the nth mentary Material and Ref. [29]. frequency component.´ We take βn = 0 for n > N, where N is The experimental sequence is outlined in Fig. 1. The initial- 3 ization of the system to g,0 is realized by the post-selection TABLE I: Amplitudes α of the initial coherent states, Wigner func- on the qubit’s ground state| i and the subsequent cavity par- tion reduction factors R, success probabilities P, and fidelities F of ity measurement. The initial coherent state of the cavity is the produced truncated phase states with θN,k = 0 (Fig. 2). The un- achieved by application of a microwave pulse, which pro- certainties are determined by bootstrapping on the measured Wigner duces a displacement operation D(α), transforming the cavity functions [38]. Bootstrapping allows random sampling with replace- from the ground state 0 to α . This is followed by driving ment at each pixel point of the measured Wigner function to obtain a | i | i the qubit, initially in the ground state g , with a pulse with pool of new Wigner functions that are used to extract R and F for a N + 1 frequency components resonant| withi the correspond- statistical analysis. ing qubit-cavity transition frequencies, simultaneously imple- N, k 5, 0 6, 0 7, 0 N menting N + 1 conditional qubit rotations: ∏n=0 Ren n(βn), | α 1.63 1.74 1.85 where Re n(βn) represents the qubit rotation on the Bloch | | n| R 0.84 0.01 0.82 0.01 0.80 0.01 sphere by an angle βn around the en direction conditional on ± ± ± P 0.37 0.31 0.23 the cavity containing n photons, with en being on the equa- F 0.97 0.01 0.96 0.01 0.96 0.01 torial plane with the azimuth angle φn + nχqcτ. In order to ± ± ± have enough selectivity, each frequency component in the qubit pulse simply has a broad Gaussian envelop truncated TABLE II: R, P, F, and (∆X(θ/2))2 for the produced squeezed h i to 4σ = 1.44 µs (σ f = 0.44 MHz). Post-selecting on the states with r = 0.8 (Fig. 3). All have an initial coherent state ampli- excited qubit state, we obtain the desired cavity state. Af- tude α = 1.45. ter preparation of the desired cavity state, the qubit is dis- | | θ 0 π/2 π entangled with the cavity and can be used to measure the R 0.82 0.01 0.84 0.01 0.82 0.01 Wigner function of the cavity for examination of the produced ± ± ± state. Following a technique devised by Lutterbach and Davi- P 0.15 0.14 0.14 dovich [34] and demonstrated in cavity QED [35] and circuit F 0.96 0.01 0.94 0.01 0.96 0.01 ± ± ± QED [26, 27, 29], the Wigner quasiprobability distribution (∆X(θ/2))2 0.067 0.001 0.070 0.001 0.070 0.001 W(β) is measured by a combination of the cavity’s displace- h i ± ± ± ment operation D( β) and qubit Ramsey-type measurement, the qubit’s excited state e , its fidelity is actually insensitive where a conditional− cavity π phase shift C(π) is sandwiched to the decoherence of the| qubiti (see Supplementary Material) in between two unconditional qubit rotations R (π/2). y as demonstrated by the measured results in Table I. As an example, we prepared the truncated phase states [36] In Fig. 2, we also display the real and imaginary parts of 1 N the density matrices in the Hilbert space obtained from the inθN,k θN,k = e n (6) normalized measured Wigner functions (the third and fourth √N + 1 ∑ | i n=0 + rows, respectively). As expected, the first N 1 photon num- with N = 5 to 7, where θN,k = 2kπ/(N + 1). Apart from ber states are almost equally populated, while the others re- the fundamental interest, such states are useful for imple- main unpopulated. The quantum coherence of the superposed mentation of quantum walk [37]. The Wigner functions of Fock states is manifested in the off-diagonal elements ρm,n the ideal states and the produced states with θ = 0 are with m = n. The non-vanishing imaginary parts of the off- N,k 6 shown in the first and second rows of Fig. 2, respectively. diagonal elements are mainly due to the imperfect setting We note that, mainly due to the decoherence of the system of the relative phases of different driving microwave compo- (see Supplementary Material), the contrast of the measured nents, and they can be minimized further by more careful opti- Wigner function is reduced by a factor R: W (β) = RW (β), mization. The value of the matrix element ρm,n (0 m,n N) m r ≤ ≤ where Wm(β) and Wr(β) are the measured and real Wigner slightly decreases as m and n increases since the decay rate is functions of the produced state [38]. The reduction factor proportional to (m + n)/2. The results further show that the 2 is R = Wm(β)d β. To make the results physically mean- produced states are in good agreement with the ideal ones, ingful, we´ normalize the measured Wigner functions, taking demonstrating high-fidelity generation of the desired states. Wm(β)/R as the real Wigner functions of the produced states. As another example to demonstrate the ability to generate The shapes of the Wigner functions of the produced states any superposition state, we produce the squeezed state agree very well with those of the ideal ones, implying that 1 ∞ ( eiθ tanhr)n (2n)! the quasiclassical states have been converted into the desired ξ = − 2n , (7) √ ∑ n quantum states with high accuracy. The amplitudes α of the | i coshr n=0 n!2 p | i initial coherent states, Wigner function reduction factors R, success probabilities P, and fidelities F between the produced with the squeeze parameter ξ = reiθ . For such a state the † iθ/2 iθ/2 states and desired ones are detailed in Table I. Here the fidelity fluctuation of the quadrature X(θ/2) = (a e + ae− )/2 2r is defined as F = θN,k ρp θN,k , where ρp is the density is reduced by a factor e− compared to the vacuum level, operator of the produced state reconstructed by least-square so that it can be used for high-precision measurement [13]. regression using a maximum likelihood estimation [38, 39]. The squeezed state can be well approximated by truncating The infidelity is mainly due to the decoherence of the system the expansion into a superposition of the first N + 1 Fock and imperfection of the conditional qubit rotations. However, states, with N depending on r. We approximately produce since the desired cavity state is obtained by post-selection on the squeezed states for ξ = 0.8, 0.8i, and 0.8 with a cutoff − 4
3 3 3 (a) (b) (c) 2 2 2 W(β) 1 1 1
β) β) β) 0.6 0 0 0 Im( -1 Im( -1 Im( -1 0.4 -2 -2 -2 0.2 3 3 3 0.0 2 2 2 1 1 1 -0.2 β) 0 β) 0 β) 0 -0.4 Im( Im( Im( -1 -1 -1 -0.6 -2 -2 -2
-2 -1 0 1 2 3 -2 -1 0 1 2 3 -2 -1 0 1 2 3 β) β) β) Re( Re( Re( 0.2 (d) Re(ρ) (e) Re(ρ) (f) Re(ρ) 0.2 0.2 0.2 0.1 0.1 0.1 0.15 0 0 0 0 0 0 10 10 10 0.1 5 5 5 5 5 5 10 0 10 0 10 0 Im(ρ) m n Im(ρ) m n Im(ρ) m n 0.2 0.2 0.2 0.05 0.1 0.1 0.1 0 0 0 0 0 0 0 10 10 10 5 5 5 5 5 5 m 10 0 n m 10 0 n m 10 0 n −0.05
FIG. 2: (Color online) Wigner tomography and density matrix reconstruction of the truncated phase states with θN,k = 0. The Wigner functions of the ideal and produced truncated phase states, with N = 5 to 7, are displayed in the first and second rows, respectively. Due to 2 the readout errors, the measured Wigner functions Wm(β) are reduced by a factor R compared to the real ones, where R = Wm(β)d β. We take Wm(β)/R as the real Wigner functions of the produced states. The real and imaginary parts of the density matrices, obtained´ from the real Wigner functions, are shown in the third and fourth rows. 3 2 (a) 3 3 3 (a) (b) (c) 1 2 (a) 2 2 β) 0 1 1 1 W(β) β) β) β) Im( -1 0 0 0 0.6 Im( Im( -1 -1 Im( -1 -2 0.4 3 -2 -2 -2 0.2 2 3 3 3 1 2 2 2 0.0 β) 0 1 1 1 -0.2 β) β) β) Im( -1 0 0 0 -0.4 Im( -1 Im( -1 Im( -1 -2 -0.6 -2 -2 -2 -2 -1 0 1 2 3 Re(β) -2 -1 0 1 2 3 -2 -1 0 1 2 3 -2 -1 0 1 2 3 Re(β) Re(β) Re(β)
FIG. 3: (Color online) Wigner tomography of the squeezed states. Reconstruction of the squeezed states defined by Eq. 7 for ξ = 0.8 (a), 0.8i (b), and 0.8 (c). We take the cutoff in the Fock-state expansion to be N = 8. The calculated and measured Wigner functions are displayed in the upper and− lower rows, respectively.
of N = 8. The good agreement in the calculated and measured initial quasiclassical states have been converted into quantum Wigner functions (Fig. 3) indicates the full controllability of states with strong squeezing. both the amplitude and phase associated with each Fock state component. Table II shows the corresponding R, P, F, and We note that, in the effective Hamiltonian of Eq. 1, we have (∆X(θ/2))2 [obtained from the Gaussian fit of the probabil- ignored the cross-Kerr terms between the cavity and the read- h i † † † ( / ) out cavity, χcra aar ar, where ar and ar are the creation ity distribution of X θ 2 calculated by integrating the mea- − sured Wigner function over X(θ/2 + π/2)] in the produced and annihilation operators of the readout cavity. This effect states. The measured (∆X(θ/2))2 1/4, showing that the can either be corrected by properly setting the initial phase h i of α or the qubit driving phases. We have also ignored | i 5 higher order terms, among which the dominant terms are ble in quantum walk, which allows exponential speedup over the cavity’s self-Kerr effect Ka†2a2/2 and the correction to classical computation for certain problems [40] and more im- − 0 †2 2 the dispersive qubit-cavity coupling χqca a e e /2. Each portantly, can be used as a primitive for universal quantum of these terms causes additional photon-number-dependent| ih | computation [41]. phase shifts; however, the resulting effect is small, as evi- denced by the high-fidelity of our experimental result. In summary, we have proposed and experimentally demon- strated a new single-step scheme to produce arbitrary super- Acknowledgments position of Fock states for a cavity. We show that a coher- ent state can be conditionally converted to any desired su- We acknowledge Zhangqi Yin, Andrei Petrenko, Brian perposition state by tailoring the probability amplitudes as- Vlastakis, and Haidong Yuan for helpful discussions, and sociated with the superposed Fock state components. This thank Tanay Roy and Madhavi Chand for the help with the is achieved based on the experimentally convenient disper- parametric amplifier measurements. This work was sup- sive coupling between the cavity and the qubit. The success ported by the Ministry of Science and the Ministry of Ed- probability of the scheme depends on the desired state. The ucation of China through its grant to Tsinghua University, achieved high fidelity of the experimental states prepared with the National Natural Science Foundation of China under this method opens up interesting applications of such non- Grant Nos. 11674060 and 11474177, the Major State Ba- classical states for implementation of quantum information sic Research Development Program of China under Grant and precision measurements. Of particular interest are the No.2012CB921601, and the 1000 Youth Fellowship program demonstrated truncated phase states that are directly applica- in China.
[1] V. Buzekˇ and P. L. Knight, Progress in Optics 34, 1 (1995). E. Lucero, M. Neeley, A. D. O’Connell, D. Sank, J. Wenner, [2] J. I. Cirac and P. Zoller, Phys. Rev. Lett. 74, 4091 (1995). J. M. Martinis, et al., Nature 459, 546 (2009). [3] J. M. Raimond, M. Brune, and S. Haroche, Rev. Mod. Phys. 73, [19] S. Krastanov, V. V. Albert, C. Shen, C.-L. Zou, R. W. Heeres, 565 (2001). B. Vlastakis, R. J. Schoelkopf, and L. Jiang, Phys. Rev. A 92, [4] M. A. Sillanpa¨a,¨ J. I. Park, and R. W. Simmonds, Nature 449, 040303 (2015). 438 (2007). [20] R. W. Heeres, B. Vlastakis, E. Holland, S. Krastanov, V. V. Al- [5] J. Majer, J. M. Chow, J. M. Gambetta, J. Koch, B. R. Johnson, bert, L. Frunzio, L. Jiang, and R. J. Schoelkopf, Phys. Rev. Lett. J. A. Schreier, L. Frunzio, D. I. Schuster, A. A. Houck, A. Wall- 115, 137002 (2015). raff, et al., Nature 449, 443 (2007). [21] N. Khaneja, T. Reiss, C. Kehlet, T. Schulte-Herbruggen,¨ and [6] L. DiCarlo, J. M. Chow, J. M. Gambetta, L. S. Bishop, B. R. S. J. Glaser, Journal of Magnetic Resonance 172, 296 (2005). Johnson, D. I. Schuster, J. Majer, A. Blais, L. Frunzio, S. M. [22] P. de Fouquieres, S. Schirmer, S. Glaser, and I. Kuprov, Journal Girvin, et al., Nature 460, 240 (2009). of Magnetic Resonance 212, 412 (2011). [7] M. F. Santos, Phys. Rev. Lett. 95, 010504 (2005). [23] R. W. Heeres, P. Reinhold, N. Ofek, L. Frunzio, L. Jiang, M. H. [8] Z. Leghtas, G. Kirchmair, B. Vlastakis, R. J. Schoelkopf, M. H. Devoret, and R. J. Schoelkopf, arXiv:1608.02430 (2016). Devoret, and M. Mirrahimi, Phys. Rev. Lett. 111, 120501 [24] S. Deleglise,´ I. Dotsenko, C. Sayrin, J. Bernu, M. Brune, J. M. (2013). Raimond, and S. Haroche, Nature 455, 510 (2008). [9] N. Ofek, A. Petrenko, R. Heeres, P. Reinhold, Z. Leghtas, [25] H. Paik, D. I. Schuster, L. S. Bishop, G. Kirchmair, G. Cate- B. Vlastakis, Y. Liu, L. Frunzio, S. M. Girvin, L. Jiang, et al., lani, A. P. Sears, B. R. Johnson, M. J. Reagor, L. Frunzio, L. I. Nature 536, 441 (2016). Glazman, et al., Phys. Rev. Lett. 107, 240501 (2011). [10] M. H. Michael, M. Silveri, R. T. Brierley, V. V. Albert, [26] B. Vlastakis, G. Kirchmair, Z. Leghtas, S. E. Nigg, L. Frun- J. Salmilehto, L. Jiang, and S. M. Girvin, Phys. Rev. X 6, zio, S. M. Girvin, M. Mirrahimi, M. H. Devoret, and R. J. 031006 (2016). Schoelkopf, Science 342, 607 (2013). [11] S. Lloyd and S. L. Braunstein, Phys. Rev. Lett. 82, 1784 (1999). [27] L. Sun, A. Petrenko, Z. Leghtas, B. Vlastakis, G. Kirchmair, [12] D. Gottesman, A. Kitaev, and J. Preskill, Phys. Rev. A 64, K. M. Sliwa, A. Narla, M. Hatridge, S. Shankar, J. Blumoff, 012310 (2001). et al., Nature 511, 444 (2014). [13] C. M. Caves, K. S. Thorne, R. W. P. Drever, V. D. Sandberg, [28] A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, R.-S. Huang, and M. Zimmermann, Rev. Mod. Phys. 52, 341 (1980). J. Majer, S. Kumar, S. M. Girvin, and R. J. Schoelkopf, Nature [14] K. Vogel, V. M. Akulin, and W. P. Schleich, Phys. Rev. Lett. 71, 431, 162 (2004). 1816 (1993). [29] K. Liu, Y. Xu, W. Wang, Z. Shi-Biao, R. Tanay, K. Suman, [15] C. K. Law and J. H. Eberly, Phys. Rev. Lett. 76, 1055 (1996). C. Madhavi, A. Ranadive, R. Vijay, Y. P. Song, et al., [16] A. Ben-Kish, B. DeMarco, V. Meyer, M. Rowe, J. Britton, arXiv:1608.04908 (2016). W. M. Itano, B. M. Jelenkovic,´ C. Langer, D. Leibfried, [30] M. Hatridge, R. Vijay, D. H. Slichter, J. Clarke, and I. Siddiqi, T. Rosenband, et al., Phys. Rev. Lett. 90, 037902 (2003). Phys. Rev. B 83, 134501 (2011). [17] M. Hofheinz, E. M. Weig, M. Ansmann, R. C. Bialczak, [31] T. Roy, S. Kundu, M. Chand, V. A. M., A. Ranadive, N. Nehra, E. Lucero, M. Neeley, A. D. O’Connell, H. Wang, J. M. Marti- M. P. Patankar, J. Aumentado, A. A. Clerk, and R. Vijay, Appl. nis, and A. N. Cleland, Nature 454, 310 (2008). Phys. Lett. 107, 262601 (2015). [18] M. Hofheinz, H. Wang, M. Ansmann, R. C. Bialczak, [32] A. Kamal, A. Marblestone, and M. H. Devoret, Phys. Rev. B 6
79, 184301 (2009). Phys. Rev. A 67, 042305 (2003). [33] K. W. Murch, S. J. Weber, C. Macklin, and I. Siddiqi, Nature [38] B. Vlastakis, A. Petrenko, N. Ofek, L. Sun, Z. Leghtas, 502, 211 (2013). K. Sliwa, Y. Liu, M. Hatridge, J. Blumoff, L. Frunzio, et al., [34] L. G. Lutterbach and L. Davidovich, Phys. Rev. Lett. 78, 2547 Nat. Commun. 6 (2015). (1997). [39] J. A. Smolin, J. M. Gambetta, and G. Smith, Phys. Rev. Lett. [35] P. Bertet, A. Auffeves, P. Maioli, S. Osnaghi, T. Meunier, 108, 070502 (2012). M. Brune, J. M. Raimond, and S. Haroche, Phys. Rev. Lett. [40] A. M. Childs, R. Cleve, E. Deotto, E. Farhi, S. Gutmann, and 89, 200402 (2002). D. A. Spielman, Proc. 35th ACM Symposium on Theory of [36] D. T. Pegg and S. M. Barnett, Journal of Modern Optics 44, 225 Computing pp. 59–68 (2003). (1997). [41] A. M. Childs, Phys. Rev. Lett. 102, 180501 (2009). [37] B. C. Sanders, S. D. Bartlett, B. Tregenna, and P. L. Knight, Supplementary Material for “Converting quasiclassical states into arbitrary Fock state superpositions in a superconducting cavity”
W. Wang,1 L. Hu,1 Y. Xu,1 K. Liu,1 Y. Ma,1 Shi-Biao Zheng,2 R. Vijay,3 Y. P. Song,1 L.-M. Duan,1, 4 and L. Sun1 1Center for Quantum Information, Institute for Interdisciplinary Information Sciences, Tsinghua University, Beijing 100084, China 2Department of Physics, Fuzhou University, Fuzhou 350116, China 3Department of Condensed Matter Physics and Materials Science, Tata Institute of Fundamental Research, Homi Bhabha Road, Mumbai 400005, India 4Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, USA
I. DEVICE PARAMETERS the histogram of the qubit readout. The histogram is clearly bimodal and well separated. A threshold Vth = 0 is chosen Our experimental device consists of two waveguide cavities to digitize the readout signal to +1 and 1 for the ground state g and the excited state e , respectively.− The qubit strongly coupled to a transmon qubit which is lithographically | i | i patterned on a sapphire chip, as shown in Fig. 1. One cavity shows an unexpected excited state population of 3.9% in the called storage cavity is used for the preparation, manipula- steady state, presumably due to stray infrared photons or other tion, and storage of the photonic state, with a transition fre- background noise leaking into the cavity. Together with the measured T1 = 13.0 µs, we can derive 1/Γ 330 µs and quency ωs/2π = 8.230 GHz and a decay rate κs/2π = 2 kHz. ↑ ≈ The other cavity noted as the readout cavity with a transi- 1/Γ 13.5 µs. We use measurement-based post-selection to purify↓ ≈ the qubit to g state. After purification, the subsequent tion frequency ωr/2π = 7.290 GHz and a decay rate κr/2π = | i 3.62 MHz, is used to detect the state of the qubit. The trans- qubit measurement shows that the probability of the qubit be- ing populated in g is as high as 0.999 (dashed histogram mon qubit has a transition frequency ωq/2π = 5.345 GHz in Fig. 2a), demonstrating| i the high quantum non-demolition with an anharmonicity αq/2π = (ωge ωe f )/2π = 246 MHz, − property of the qubit readout. If we post-select the qubit at an energy relaxation time T1 = 13.0 µs and a pure dephasing e state, the subsequent measurement finds the qubit remains time Tφ = 13.8 µs. The strong couplings between the qubit | i and the cavities result in state-dependent frequency shifts in e with a fidelity of 0.953 (Fig. 2b), with the errors dom- inantly| i from the T process during the readout time (240 ns) χqs/2π = 1.44 MHz (between qubit and storage cavity) and 1 − and during the waiting time after the initialization measure- χqr/2π = 4.71 MHz (between qubit and readout cavity). For simplicity,− we refer to the storage cavity as “the cavity” ment (250 ns). These readout properties are measured with henceforth. the cavity left in vacuum.
Figure 2c shows the property of a π pulse (σ = 4 ns) which II. READOUT PROPERTIES is used to reset the qubit after generating the desired state by post-selection of the qubit at e state. Figure 2d shows In order to have a high-fidelity single-shot readout, besides the associated measurement pulses.| i The numbers outside the the high-electron-mobility-transistor amplifier at 4K, we have parenthesis correspond to the cavity in the vacuum state, with a Josephson parametric amplifier (JPA) [1, 2] connected to the the loss of fidelity mainly due to qubit T1 process during the output of the readout cavity as the first stage of amplification. measurement time and the waiting time between the two con- The JPA works in a pulsed double-pumped mode [3, 4] to min- secutive measurements. The numbers inside the parenthesis imize pump leakage to the readout cavity. The schematic of correspond to the cavity in a coherent state with an average the measurement setup can be found in Ref. 5. photon number ofn ¯ = 3.5, with about 2% lower fidelities due The characterizations of the readout properties of the qubit to the photon number occupations in the cavity. arXiv:1703.03316v1 [quant-ph] 9 Mar 2017 and cavity photon state are shown in Fig. 2. Figure 2a shows
The Wigner quasiprobability distribution W(β) is mea- sured by a combination of the cavity’s displacement operation Readout rr44 ns f 7.29 GHz D( β) and a parity measurement, following a technique de- Qubit vised− by Lutterbach and Davidovich [6] and demonstrated in X fq 5.345 GHz cavity QED [7] and circuit QED [5, 8, 9]. The parity mea- Storage 80 μs f 8.23 GHz ss surement is achieved in a qubit Ramsey-type measurement, where a conditional cavity π phase shift C(π) is sandwiched Y in between two unconditional qubit rotations Rπ/2. The par- FIG. 1: Experimental device. A cross section and diagram of the ity readout property is shown in Fig. 3. The fidelity loss is Y device, which contains two 3D Al cavities and a sapphire chip with dominantly from the qubit decoherence between the two Rπ/2 a transmon qubit located in a trench. The transmon qubit is strongly rotations. T1 and T2 processes together contribute about 3.8% coupled to both cavities. error. 2
(a)
0.999 푌 g g 푅휋 (b) 0.001 (d) 4휎 = 16 ns Qubit 250 ns 0.047 푒 푒 M2 0.953 M1 Readout 240 ns 240 ns 0.014 (0.035) (c) g g 푌 푌 0.986 (0.965) 푅휋 D() 푅휋 4휎 = 16 ns Qubit 250 ns 0.948 (0.927) 푒 푒 M1 M2 0.052 (0.073) Readout 240 ns 240 ns
FIG. 2: Readout properties. (a) Bimodal and well-separated histogram of the qubit readout. A threshold Vth = 0 is chosen to digitize the readout signal. The solid line corresponds to the case without the qubit purification measurement, showing the qubit has a 3.9% probability of being populated in e state, while the dashed line represents the result after the qubit has been purified. The inset is the measurement pulse | i sequence. (b) Basic qubit readout matrix with the cavity left in vacuum. The readout errors are mainly due to the T1 process. (c) Readout property of a π pulse with a Gaussian envelope truncated to 4σ = 16 ns when the cavity is in the vacuum state. Numbers in parenthesis correspond to the case withn ¯ = 3.5 in the cavity. The loss of fidelity mainly due to qubit T1 process during the measurement time and the waiting time between the two consecutive measurements. There is about an extra 2% fidelity loss for the case withn ¯ = 3.5 photons in the cavity. (d) The pulse sequence used for the measurement in (c).
III. MEASUREMENT PULSE SEQUENCE ble qubit leakage to higher qubit levels, we also apply the so- called “derivative removal by adiabatic gate” technique [10] = , The pulse sequence for the experiment is shown in Fig. 4. for pulses with σ 2 4, and 6 ns. The whole experiment can be divided into three main parts: 1) initialization of the system to g,0 by post-selecting re- sults of the qubit measurement and| thei cavity parity measure- IV. CALIBRATION OF DRIVING AMPLITUDE ment; 2) creation of the desired arbitrary Fock-state superpo- sitions by first displacing the cavity to α state, then rotating The calibration of the qubit driving amplitude is essential | i the qubit simultaneously and independently in the subspaces in creating a high-fidelity Fock-state superposition with our g,n , e,n , and finally post-selecting the excited state of protocol. This calibration is performed by a typical Rabi os- {| i | i} the qubit; 3) Wigner tomography measurement following a π cillation with a selective driving pulse (σ = 360 ns) at each pulse reset of the qubit to the ground state. The rotation of resonant transition frequency of the qubit corresponding to a qubit Re n(βn) (already defined in the main text) in each sub- specific photon number n in the cavity (inset of Fig. 5). The n| space g,n , e,n is appropriately chosen depending on the qubit decoherence during the long rotation only reduces the {| i | i} probability amplitude of Fock state n in the desired superpo- Rabi signal contrast without affecting the calibration of the | i sition state. driving amplitude with a deviation typically less than 1%. The classical microwave pulse to displace the cavity state The calibration of the driving amplitude to create a desired has a narrow square envelop with a width of 20 ns. The am- coherent state in the cavity is also crucial. This calibration plitude has been calibrated carefully as shown in the section of is performed by measuring the Poisson distribution of photon “calibration of driving amplitude”. All qubit drive pulses have numbers in the cavity (Fig. 5). We measure the probability a Gaussian envelope truncated to 2σ. To eliminate the possi- of the first seven Fock states N = 0,1,2,...6 through a selec- ± 3
D() (a) 20 ns
Storage
푌 푌 푅휋/2 푅±휋/2 337 ns Qubit Initialization 4휎 = 8 ns 240 ns
Readout P2 P3 P4 P5 P6
P1
(b) (c) 푌 푌 푌 푌 Cavity parity readout 푅휋/2 , π/χ , 푅휋/2 Cavity parity readout 푅휋/2 , π/χ , 푅−휋/2 Prepared state Prepared state 0.958 0.967 훼 + −훼 Even +1 훼 + −훼 Even +1 0.042 0.033
0.043 0.050 훼 − −훼 Odd -1 훼 − −훼 Odd -1 0.957 0.950
FIG. 3: Parity readout properties. (a) Protocol to measure parity readout fidelity. An initialization creates the g,0 state, followed by a | i cavity displacement D(α = 1) and six consecutive parity measurements. Post-selection of the first five consecutive identical parity results would give a photon state parity with good confidence, and the sixth parity measurement is used for the final parity fidelity measurement. The pulse sequence for each parity measurement has been shown in P1. (b) Parity readout property for given even and odd parity states for two Y protocols ( R π/2 as the second qubit pulse) used in the tomography process. To minimize the system error in the tomography, we take both protocols and± average the measurement results together.
Initialization to g,0 State generation Tomography
D D
Storage
20 ns 20 ns 푌 푅휋 푅푌 푅푌 푌,푛 푅푌 푅푌 휋/2 −휋/2 푅휋 휋/2 ±휋/2 Qubit 320 ns 337 ns 240 ns 4휎 = 24 ns 240 ns 4휎 = 1440 ns 240 ns 4휎 = 16 ns 4휎 = 8 ns 240 ns 250 ns 250 ns 250 ns Readout
M1 M2 M3 M4
FIG. 4: Schematic of the pulse sequence in the experiments. The whole experiment can be divided into three main parts: 1. initialization of the system to g,0 by post-selecting results of the qubit measurement and the cavity parity measurement; 2. creation of the desired arbitrary | i Fock-state superpositions by first displacing the cavity to α state, then rotating the qubit simultaneously and independently in the subspaces | i g,n , e,n , and finally post-selecting the excited state of the qubit; 3. Wigner tomography measurement following a π pulse reset of the {qubit| i to| thei ground} state. tive π pulse (σ = 360 ns) as a function of the displacement the coherent state in the cavity. pulse amplitude but with a fixed width of 20 ns. A nearly perfect global fit to the measurement results with a Poisson distribution not only gives the required calibration of the cav- ity displacement amplitude, but also indicates good control of 4
1.0
1.0
0.8 0.8
0.6 pulse Pg 0.4 0.6 0.2
0.0 -40 -20 0 20 40 Probability 0.4 Qubit drive amplitude (DAC)
0.2
0.0 0 1000 2000 3000 4000 5000 6000 Displacement pulse amplitude (DAC)
FIG. 5: Poisson distribution of photon numbers in the cavity. Dots are data for the probability of the first seven Fock states N = 0,1,2,...6 as a function of the displacement pulse amplitude but with a fixed width 20 ns. The probability measurement is performed with a selective π pulse on the qubit at each resonant transition frequency corresponding to a specific photon number n in the cavity. There is no normalization of the measurement data and the loss of probability dominantly comes from the T1 process in the long duration of the π pulse. Lines are from a global 2N α 2 fit with a Poisson distribution P( α ,N) = A α e−| | /N!, where A is a scaling factor accounting for the probability loss. The excellent | i agreement indicates good control of the coherent| | state in the cavity. The displacement amplitude α is proportional to the measurement pulse amplitude. In this case, DAC level=1343 is calibrated to correspond to α = 1. Inset: Rabi oscillation to calibrate the qubit drive amplitude (σ = 360 ns). The loss of contrast comes from the qubit relaxation in the long pulse duration.
TABLE I: Error budgets for the final Wigner tomography used to examine the produced states. These error sources in total account for the large reduction of R. source of errors duration (µs) coherence time (µs) contribution to reduction of R from 1 qubit T between M3 and π pulse 0.370 13.0 5.7% 1 ∼ qubit T1 between two π/2 pulses in tomography 0.337 13.0 2.6% qubit Tφ between two π/2 pulses in tomography 0.337 13.8 4.9% inaccuracy of π pulse 0.016 NA 2% 4% − cavity decay between two π/2 pulses in tomography 0.337 80/n¯ 1% ∼
V. ANALYSIS OF ERRORS duced state is insensitive to the qubit decay as this process cannot mix the cavity state correlated with the qubit’s ground state into that with the qubit’s excited state; instead, this de- Our selective pulse is relatively slow (about 1.44 µs) due cay mainly leads to slight changes of the coefficients of the to the small dispersive frequency shift, however, to generate superposed Fock states, to which the fidelity of the produced the truncated phase states and the squeezed states presented in state is robust. A numerical simulation shows this T1 decay re- our main text our scheme is not sensitive to the qubit energy duces the fidelity by only 1-3%. On the other hand, the qubit relaxation T1, qubit dephasing Tφ , and cavity photon energy dephasing corresponds to a random qubit frequency fluctua- relaxation. First of all, we apply microwave drives to qubit tion, which just results in slight changes of the probability only associated with Fock states present in the target state and amplitudes associated with the Fock state components of the only select the qubit’s excited state. The probability of select- produced state. The fidelity of the produced state is insensi- ing unwanted Fock states in the final state (orthogonal to the tive to such slight changes and thus robust against the qubit desired state) can be neglected because the probability of Γ dephasing. This is confirmed by numerical simulation, which process ( g to e ) is only about 1.44µs/330µs =0.4%. ↑ | i | i shows that the fidelity F of the produced cavity state is even Secondly, random T1 decay from e to g during the more insensitive to the qubit dephasing than to the qubit en- photon-number-dependent qubit rotation| i obviously| i decreases ergy relaxation. the probability of the qubit being excited; however, the fi- Thirdly, we consider the effect from the cavity photon de- delity of the produced state is insensitive to the qubit energy cay. This effect increases with the photon number (as shown decay. This has a simple explanation: The purity of the pro- in Fig. 2 of the main text, the matrix element ρm,n has a smaller 5 value for larger m and n), however, numerical simulations two comes from the fact that parity is from -1 to 1. 2) Qubit T1 show that with our long photon lifetime (80 µs), photon de- process during the parity measurement causes a reduction of cay actually contributes less to the fidelity loss than the qubit 0.337µs/13µs=2.6%. 3) Qubit dephasing Tφ during the parity decoherence for the generated states in the main text. measurement causes a reduction of 0.337µs/13.8µs 2=4.9%. Therefore, after considering all the above contributions to 4) Due to the inaccuracy of the π pulse (calibrated× when the fidelity loss, the overall fidelity could still be as high as the cavity is in the vacuum state), the qubit has a probabil- 96%. ity of 1-2% to remain in the exited state e for large N due | i In contrast, both T1 and Tφ process cause a much more sig- to the dispersive frequency shift (see Fig. 2c), leading to a nificant reduction of R in the following Wigner tomography reduction of 2-4%. 5) We note that the cavity photon de- to reconstruct the cavity state. The error budgets are shown in cay only contributes to the reduction of R during the disper- Table I. Explicitly, 1) the qubit T1 process between the final sive qubit-cavity interaction sandwiched in between the final detection of e to generate the target state and the tomog- two Ramsey pulses. This contribution is roughly estimated as raphy causes| ai reduction of 0.370µs/13µs 2=5.7%. Here n¯ 0.337µs/80µs 1%. With all of the effects being included, × × ∼ 0.370 µs includes half of the measurement time of M3 to esti- the resulting Wigner function reduction factor is about 0.83, mate the decoherence during the measurement, and the factor in agreement with the measured results.
[1] M. Hatridge, R. Vijay, D. H. Slichter, J. Clarke, and I. Siddiqi, (1997). Phys. Rev. B 83, 134501 (2011). [7] P. Bertet, A. Auffeves, P. Maioli, S. Osnaghi, T. Meunier, [2] T. Roy, S. Kundu, M. Chand, V. A. M., A. Ranadive, N. Nehra, M. Brune, J. M. Raimond, and S. Haroche, Phys. Rev. Lett. M. P. Patankar, J. Aumentado, A. A. Clerk, and R. Vijay, Appl. 89, 200402 (2002). Phys. Lett. 107, 262601 (2015). [8] B. Vlastakis, G. Kirchmair, Z. Leghtas, S. E. Nigg, L. Frun- [3] A. Kamal, A. Marblestone, and M. H. Devoret, Phys. Rev. B zio, S. M. Girvin, M. Mirrahimi, M. H. Devoret, and R. J. 79, 184301 (2009). Schoelkopf, Science 342, 607 (2013). [4] K. W. Murch, S. J. Weber, C. Macklin, and I. Siddiqi, Nature [9] L. Sun, A. Petrenko, Z. Leghtas, B. Vlastakis, G. Kirchmair, 502, 211 (2013). K. M. Sliwa, A. Narla, M. Hatridge, S. Shankar, J. Blumoff, [5] K. Liu, Y. Xu, W. Wang, Z. Shi-Biao, R. Tanay, K. Suman, et al., Nature 511, 444 (2014). C. Madhavi, A. Ranadive, R. Vijay, Y. P. Song, et al., [10] F. Motzoi, J. M. Gambetta, P. Rebentrost, and F. K. Wilhelm, arXiv:1608.04908 (2016). Phys. Rev. Lett. 103, 110501 (2009). [6] L. G. Lutterbach and L. Davidovich, Phys. Rev. Lett. 78, 2547