Circuit quantum electrodynamics : beyond the linear dispersive regime

Transcription

1 Circuit quantum electrodynamics : beyond the linear dispersive regime 1 Jay Gambetta 2 Alexandre Blais 1 1 Département de Physique et Regroupement Québécois sur les matériaux de pointe, 2 Institute for Quantum Computing and Department of Physics and Astronomy, University of Waterloo June 23 th, 2008 Boissonneault, Gambetta and Blais, Phys. Rev. A (R) (2008)

2 1 Introduction Atom and cavity Cavity QED Charge qubit and coplanar resonator Circuit QED The linear dispersive limit Circuit VS cavity QED 2 Beyond linear dispersive Understanding the dispersive transformation The dispersive limit Dissipation in the system Dissipation in the transformed basis 3 Results Reduction of the SNR Measurement induced heat bath The case of the transmon 4 Conclusion Conclusion

3 Atom and cavity... Energy ω 01 = ω a Two-levels system Hamiltonian H = ωa σz σz = «

4 Atom and cavity x y z Energy... L Two-levels system Hamiltonian H = ωa σz σz = 2 ω 01 = ω a «Cavity Hamiltonian H = X k ω k a k a k + 1 «2

5 Atom and cavity x y z Energy... L Two-levels system Hamiltonian H = ωa σz σz = 2 ω 01 = ω a «Cavity Hamiltonian H = X k ω k a k a k + 1 «2

6 Atom and cavity ω r = ω 1 ω 2 Energy Response κ = ω r/q... Input frequency, rf Two-levels system Hamiltonian H = ωa σz σz = 2 ω 01 = ω a «Cavity Hamiltonian H = X k Single-mode : ω k a k a k + 1 «2 H = ω ra a

7 Cavity QED g

8 Cavity QED g Atom-cavity interaction H I = D E g(a + a)σ x g(a σ + aσ + ) g(z) = d 0 r ω V ǫ 0 sinkz

9 Cavity QED g Atom-cavity interaction H I = D E g(a + a)σ x g(a σ + aσ + ) Jaynes-Cummings Hamiltonian g(z) = d 0 r ω V ǫ 0 sinkz H = ωa 2 σz + ωra a + g(a σ + aσ + ) Jaynes and Cummings, Proc. IEEE (1963) Raimond, Brune and Haroche, Rev. Mod. Phys (2001) Mabuchi and Doherty, Science (2002)

10 Charge qubit and coplanar resonator Classical Hamiltonian Vg C g E J C J Vg C g n E J C J E C = H = 4E C (n n g) 2 E J cos δ E J = I 0Φ 0 2π e 2 CgVg, ng = 2(C g + C J ) 2e

11 Charge qubit and coplanar resonator Classical Hamiltonian Vg C g E J C J Vg C g n E J C J E C = H = 4E C (n n g) 2 E J cos δ E J = I 0Φ 0 2π e 2 CgVg, ng = 2(C g + C J ) 2e Quantum Hamiltonian H = X n 4E C (n n g) 2 n n X n E J ( n n h.c.) 2 Restricting to n g [0,1] : H = ω aσ z/2 Shnirman, Schön and Hermon, Phys. Rev. Lett (1997) Bouchiat et al., Physica Scripta T (1998) Nakamura, Pashkin and Tsai, Nature (London) (1999)

12 Charge qubit and coplanar resonator Classical Hamiltonian Vg C g E J C J Vg C g n E J C J E C = H = 4E C (n n g) 2 E J cos δ E J = I 0Φ 0 2π e 2 CgVg, ng = 2(C g + C J ) 2e Energie [Arb. Units] EJ/4EC=0.1 EJ Gate charge, n g = C g V g /2e Quantum Hamiltonian H = X n X n 4E C (n n g) 2 n n E J ( n n h.c.) 2 Restricting to n g [0,1] : H = ω aσ z/2 Shnirman, Schön and Hermon, Phys. Rev. Lett (1997) Bouchiat et al., Physica Scripta T (1998) Nakamura, Pashkin and Tsai, Nature (London) (1999)

13 Charge qubit and coplanar resonator

14 Charge qubit and coplanar resonator L r C r Classical Hamiltonian H = Φ2 2L r CrV 2 ω r = s 1 L rc r

15 Charge qubit and coplanar resonator L r C r Classical Hamiltonian H = Φ2 2L r CrV 2 ω r = s 1 L rc r Quantum Hamiltonian r ωr V = (a + a), 2C r H = ω r a a + 1 «2 r ωr Φ = i (a a) 2L r Quantum Fluctuations in Electrical Circuits, M. H. Devoret, Les Houches Session LXIII, Quantum Fluctuations p (1995).

16 Circuit QED Measurement output C g E J CJ Qubit control and readout ~ 10 GHz Atom: super conducting charge qubit Cavity: super conducting 1D transmission line resonator

17 Circuit QED Measurement output C g E J CJ Qubit control and readout ~ 10 GHz Atom: super conducting charge qubit Cavity: super conducting 1D transmission line resonator Parameters g : Qubit-cavity interaction ω a : Qubit frequency ω r : Resonator frequency = ω a ω r : Detuning H = ωa 2 σz + ωra a + g(a σ + aσ + ) Blais et al., Phys. Rev. A (2004) Wallraff et al., Nature (2004) Wallraff et al., Phys. Rev. Lett (2005) Leek et al., Science (2007) Schuster et al., Nature (2007) Houck et al., Nature (2007) Majer et al., Nature (2007)

18 Circuit QED C g E J CJ Parameters g : Qubit-cavity interaction ω a : Qubit frequency ω r : Resonator frequency = ω a ω r : Detuning H = ωa 2 σz + ωra a + g(a σ + aσ + ) Blais et al., Phys. Rev. A (2004) Wallraff et al., Nature (2004) Wallraff et al., Phys. Rev. Lett (2005) Leek et al., Science (2007) Schuster et al., Nature (2007) Houck et al., Nature (2007) Majer et al., Nature (2007)

19 The linear dispersive limit Jaynes-Cummings H = ω ra a+ω a σ z 2 +g(a σ +aσ + ) Small parameter λ = g/ Measurement output Qubit control and readout ~ 10 GHz Atom: super conducting charge qubit Cavity: super conducting 1D transmission line resonator

20 The linear dispersive limit Jaynes-Cummings H = ω ra a+ω a σ z 2 +g(a σ +aσ + ) Small parameter λ = g/ Measurement output Linear dispersive H D = (ω a + χ ) σz 2 + (ωr + χσz )a a Lamb shift (χ = gλ = g 2 / ) Stark shift or cavity pull Valid if n n crit., where n crit. = 1/4λ 2. Qubit control and readout ~ 10 GHz Atom: super conducting charge qubit Cavity: super conducting 1D transmission line resonator 2CP T ransmi ssio n (arb. units) κ Δ ~ 2π 1GHz 2χ κ Phase -2χ κ ωr CP ωr + CP ω a

21 The linear dispersive limit Jaynes-Cummings H = ω ra a+ω a σ z 2 +g(a σ +aσ + ) Small parameter λ = g/ Measurement output Linear dispersive H D = (ω a + χ ) σz 2 + (ωr + χσz )a a Lamb shift (χ = gλ = g 2 / ) Stark shift or cavity pull Valid if n n crit., where n crit. = 1/4λ 2. Qubit control and readout ~ 10 GHz Atom: super conducting charge qubit Cavity: super conducting 1D transmission line resonator 2CP T ransmi ssio n (arb. units) κ Δ ~ 2π 1GHz 2χ κ Phase -2χ κ ωr CP ωr + CP ω a

22 The linear dispersive limit Jaynes-Cummings H = ω ra a+ω a σ z 2 +g(a σ +aσ + ) Small parameter λ = g/ Measurement output Linear dispersive H D = (ω a + χ ) σz 2 + (ωr + χσz )a a Lamb shift (χ = gλ = g 2 / ) Stark shift or cavity pull Qubit control and readout ~ 10 GHz Atom: super conducting charge qubit Cavity: super conducting 1D transmission line resonator Valid if n n crit., where n crit. = 1/4λ 2. Rabi π-pulse Wallraff et al., Phys. Rev. Lett (2005) 2CP T ransmi ssio n (arb. units) κ ωr CP Δ ~ 2π 1GHz ωr + CP ω a 2χ κ Phase -2χ κ Averaged times. SNR for single-shot is 0.1.

23 Circuit VS cavity QED Symbol Optical cavity Microwave cavity Circuit ω r/2π or ω a/2π 350 THz 51 GHz 10 GHz g/π 220 MHz 47 khz 100 MHz g/ω r Hood et al., Science (2000) Raimond, Brune and Haroche, Rev. Mod. Phys (2001) Blais et al., Phys. Rev. A (2004)

24 Circuit VS cavity QED Motivation Symbol Optical cavity Microwave cavity Circuit ω r/2π or ω a/2π 350 THz 51 GHz 10 GHz g/π 220 MHz 47 khz 100 MHz g/ω r Hood et al., Science (2000) Raimond, Brune and Haroche, Rev. Mod. Phys (2001) Blais et al., Phys. Rev. A (2004) Circuit QED is harder than cavity QED on the dispersive limit (n crit. is smaller)

25 Circuit VS cavity QED Motivation Symbol Optical cavity Microwave cavity Circuit ω r/2π or ω a/2π 350 THz 51 GHz 10 GHz g/π 220 MHz 47 khz 100 MHz g/ω r Hood et al., Science (2000) Raimond, Brune and Haroche, Rev. Mod. Phys (2001) Blais et al., Phys. Rev. A (2004) Circuit QED is harder than cavity QED on the dispersive limit (n crit. is smaller) The SNR is low, we want to measure harder... how does higher order terms affect measurement?

26 Circuit VS cavity QED Motivation Symbol Optical cavity Microwave cavity Circuit ω r/2π or ω a/2π 350 THz 51 GHz 10 GHz g/π 220 MHz 47 khz 100 MHz g/ω r Hood et al., Science (2000) Raimond, Brune and Haroche, Rev. Mod. Phys (2001) Blais et al., Phys. Rev. A (2004) Circuit QED is harder than cavity QED on the dispersive limit (n crit. is smaller) The SNR is low, we want to measure harder... how does higher order terms affect measurement? Must consider higher order corrections in perturbation theory

27 Understanding the dispersive transformation J-C : block diagonal H = ω ra a + ω aσ z/2 + g(a σ + aσ +)

28 Understanding the dispersive transformation J-C : block diagonal H = ω ra a + ω aσ z/2 + g(a σ + aσ +) 1x1 block : H 0 = ω ai/2 2x2 blocks : H n = 2 σn z + g nσ n x Total Hamiltonian H = H 0 H 1 H 2 H

29 Understanding the dispersive transformation J-C : block diagonal H = ω ra a + ω aσ z/2 + g(a σ + aσ +) 1x1 block : H 0 = ω ai/2 Z θ n H n 2x2 blocks : H n = 2 σn z + g nσ n x Total Hamiltonian H = H 0 H 1 H 2 H X θ n = arctan(2g n/ )

30 Understanding the dispersive transformation J-C : block diagonal H = ω ra a + ω aσ z/2 + g(a σ + aσ +) 1x1 block : H 0 = ω ai/2 Z θ n H n 2x2 blocks : H n = 2 σn z + g nσ n x Total Hamiltonian H = H 0 H 1 H 2 H Diagonalization Rotation around Y axis In all subspaces E n E 0 = { g,0 } E n = { g, n, e, n 1 } { g n, e n } X θ n = arctan(2g n/ ) Z Z D θ n D H n D X X

31 Understanding the dispersive transformation J-C : block diagonal H = ω ra a + ω aσ z/2 + g(a σ + aσ +) 1x1 block : H 0 = ω ai/2 Z θ n H n 2x2 blocks : H n = 2 σn z + g nσ n x Total Hamiltonian H = H 0 H 1 H 2 H Diagonalization Rotation around Y axis In all subspaces E n E 0 = { g,0 } E n = { g, n, e, n 1 } { g n, e n } The qubit is now part photon and vice-versa X θ n = arctan(2g n/ ) Z Z D θ n D H n X X D

32 Understanding the dispersive transformation J-C : block diagonal H = ω ra a + ω aσ z/2 + g(a σ + aσ +) 1x1 block : H 0 = ω ai/2 Z θ n H n 2x2 blocks : H n = 2 σn z + g nσ n x Total Hamiltonian H = H 0 H 1 H 2 H Diagonalization Rotation around Y axis In all subspaces E n E 0 = { g,0 } E n = { g, n, e, n 1 } { g n, e n } The qubit is now part photon and vice-versa X θ n = arctan(2g n/ ) 2g n/ Z D Z θ n D H n X X D

33 The dispersive limit Dispersive limit Jaynes-Cummings hamiltonian H = ω ra a + ω a σ z 2 + g(a σ + aσ + ) Exact transformation : D Small parameter λ = g/ 4λ 2 n 1

34 The dispersive limit Dispersive limit Jaynes-Cummings hamiltonian H = ω ra a + ω a σ z 2 + g(a σ + aσ + ) Exact transformation : D Small parameter λ = g/ 4λ 2 n 1 Result at order λ H D = (ω a + χ ) σz 2 +(ωr + χσz )a a Lamb shift (χ = gλ = g 2 / ) Stark shift or cavity pull T ransmi ssio n (arb. units) κ ωr CP 2CP Δ ~ 2π 1GHz ωr + CP ω a 2χ κ Phase -2χ κ

35 The dispersive limit Dispersive limit Jaynes-Cummings hamiltonian H = ω ra a + ω a σ z 2 + g(a σ + aσ + ) Exact transformation : D Small parameter λ = g/ 4λ 2 n 1 Result at order λ 2 H D = (ω a + χ ) σz 2 + [ωr + (χ ζa a)σ z ]a a χ = χ(1 λ 2 ), ζ = λ 2 χ The cavity pull decrease : CP = χ ζ a a Result at order λ H D = (ω a + χ ) σz 2 +(ωr + χσz )a a Lamb shift (χ = gλ = g 2 / ) Stark shift or cavity pull T ransmi ssio n (arb. units) κ 2CP Δ ~ 2π 1GHz 2χ κ Phase -2χ κ ωr CP ωr + CP ω a

36 Dissipation in the system κ γ ϕ γ 1 Model for dissipation Coupling to a bath H κ = R p 0 gκ(ω)[b κ(ω) + b κ(ω)][a + a]dω H γ = R p 0 gγ(ω)[b γ(ω) + b γ(ω)]σ xdω Parameters κ : Rate of photon loss γ 1 : Transverse decay rate

37 Dissipation in the system κ γ ϕ Energie [Arb. Units] EJ/4EC=0.1 δn g γ Gate charge, n g = C g V g /2e Model for dissipation Coupling to a bath H κ = R p 0 gκ(ω)[b κ(ω) + b κ(ω)][a + a]dω H γ = R p 0 gγ(ω)[b γ(ω) + b γ(ω)]σ xdω Parameters κ : Rate of photon loss γ 1 : Transverse decay rate

38 Dissipation in the system κ γ ϕ Energie [Arb. Units] EJ/4EC=0.1 δn g γ Gate charge, n g = C g V g /2e Model for dissipation Coupling to a bath H κ = R p 0 gκ(ω)[b κ(ω) + b κ(ω)][a + a]dω H γ = R p 0 gγ(ω)[b γ(ω) + b γ(ω)]σ xdω Dephasing H ϕ = ηf(t)σ z Parameters κ : Rate of photon loss γ 1 : Transverse decay rate γ ϕ : Pure dephasing rate

39 Dissipation in the transformed basis Z H n θ n γ 1 γ ϕ X

40 Dissipation in the transformed basis Z θ n γ 1 γ ϕ H n X Transformation of system-bath hamiltonian a D a + λσ + O `λ2 D σ σ + λaσ z + O `λ2 D σ z σz 2λ(a σ + aσ + ) + O `λ2 Z D Z D H n θ n γ ϕeff γ γ ϕeff γ γ, X X D

41 Dissipation in the transformed basis Z D Z Z θ n γ 1 H γ ϕ D n H n X Transformation of system-bath hamiltonian a D a + λσ + O `λ2 D σ σ + λaσ z + O `λ2 D σ z σz 2λ(a σ + aσ + ) + O `λ2 Method Transform the system-bath hamiltonian Trace out heat bath and cavity degrees of freedom (Gambetta et al., Phys. Rev. A (2008)) θ n γ ϕeff γ γ ϕeff γ γ, X X D

42 Dissipation in the transformed basis Z D Z Z θ n γ 1 H γ ϕ D n H n X Transformation of system-bath hamiltonian a D a + λσ + O `λ2 D σ σ + λaσ z + O `λ2 D σ z σz 2λ(a σ + aσ + ) + O `λ2 Method Transform the system-bath hamiltonian Trace out heat bath and cavity degrees of freedom (Gambetta et al., Phys. Rev. A (2008)) New rates (assuming white noises) θ n γ ϕeff γ γ γ, X X D γ = γ 1 ˆ1 2λ 2 ` n γκ + 2λ 2 γ ϕ n γ = 2λ 2 γ ϕ n γ κ = λ 2 κ γ ϕeff

43 Reduction of the SNR Parameters for the SNR Number of measurement photons SNR Num. phot.

44 Reduction of the SNR Parameters for the SNR Number of measurement photons Output rate : κ SNR κ Num. phot.

45 Reduction of the SNR Parameters for the SNR Number of measurement photons Output rate : κ Fraction of photons detected : η κ η Num. phot. SNR

46 Reduction of the SNR Parameters for the SNR Number of measurement photons Output rate : κ Fraction of photons detected : η Info per photon : cavity pull κ η Num. phot. Info per phot. SNR

47 Reduction of the SNR Parameters for the SNR Number of measurement photons Output rate : κ Fraction of photons detected : η Info per photon : cavity pull Mixing rate : γ +γ = γ 1 ˆ1 2λ 2 ` n γκ+4λ 2 γ ϕ n SNR κ η Num. phot. Info per phot. Mixing rate

48 Reduction of the SNR Parameters for the SNR Number of measurement photons Output rate : κ Fraction of photons detected : η Info per photon : cavity pull Mixing rate : γ +γ = γ 1 ˆ1 2λ 2 ` n γκ+4λ 2 γ ϕ n SNR κ η Num. phot. Info per phot. Mixing rate SNR Cooper-Pair Box Linear With corrections Conclusion SNR levels off with non-linear effects! Explains low experimental SNR Applies to all dispersive homodyne measurement n/n crit. /2π = 1.7 GHz, g/2π = 170 MHz κ/2π = 34 MHz, γ 1 /2π = 0.1 MHz γ ϕ = 0.1 MHz, η = 1/80 n crit. = 1/4λ 2 = 25

49 Reduction of the SNR Parameters for the SNR Number of measurement photons Output rate : κ Fraction of photons detected : η Info per photon : cavity pull Mixing rate : γ +γ = γ 1 ˆ1 2λ 2 ` n γκ+4λ 2 γ ϕ n SNR Conclusion κ η Num. phot. Info per phot. Mixing rate SNR levels off with non-linear effects! Explains low experimental SNR Applies to all dispersive homodyne measurement SNR Cooper-Pair Box Linear With corrections Transmon n/n crit. /2π = 1.7 GHz, g/2π = 170 MHz κ/2π = 34 MHz, γ 1 /2π = 0.1 MHz γ ϕ = 0.1 MHz, η = 1/80 n crit. = 1/4λ 2 = 25

50 Measurement induced heat bath Mixing rates Downward rate : γ ( n) Upward rate : γ ( n) Heat bath with temperature T( n) = ( ω r/k B )/log(1 + 1/ n) σ z Time [µs]

51 Measurement induced heat bath Mixing rates Downward rate : γ ( n) Upward rate : γ ( n) Heat bath with temperature T( n) = ( ω r/k B )/log(1 + 1/ n) σ z σ z Time [µs] Increasing power

52 Measurement induced heat bath Mixing rates Downward rate : γ ( n) Upward rate : γ ( n) Heat bath with temperature T( n) = ( ω r/k B )/log(1 + 1/ n) σ z σ z σ z Time [µs] Increasing power

53 The case of the transmon C g Vg E J C J Energie [Arb. Units] EJ/4EC=0.1 EJ Gate charge, n g = C gv g/2e Koch et al., Phys. Rev. A (2007)

54 The case of the transmon C g C g Vg E J C J Vg C S Energie [Arb. Units] EJ/4EC=0.1 EJ Gate charge, n g = C gv g/2e Koch et al., Phys. Rev. A (2007)

55 The case of the transmon C g C g C g Vg E J C J Vg C S Vg C S Energie [Arb. Units] EJ/4EC=0.1 EJ Gate charge, n g = C gv g/2e Koch et al., Phys. Rev. A (2007)

56 The case of the transmon C g C g C g Vg E J C J Vg C S Vg C S Energie [Arb. Units] EJ/4EC=0.1 EJ Gate charge, n g = C gv g/2e Koch et al., Phys. Rev. A (2007)

57 Conclusion Main results Simple model describing the physics of measurement Side-effect of measuring harder : you heat your qubit (even with photons that can t be directly absorbed) Side-effect of measuring harder : each photon you add carries less information than the previous one Measuring harder bigger SNR Coming soon The transmon (3 level system) (Koch et al., Phys. Rev. A (2007)) Taking advantage of the non-linearity Comparison with experiments More information : Boissonneault, Gambetta and Blais, Phys. Rev. A (R) (2008) FQRNT