## Abstract

We study a cavity-QED setup consisting of a two-level system coupled to a single cavity mode with two-photon relaxation. The system dynamics is modeled via a Lindblad master equation consisting of the Rabi Hamiltonian and a two-photon dissipator. We show that an even-photon relaxation preserves the Z2 symmetry of the Rabi model, and provide a framework to study the corresponding non-Hermitian dynamics in the number-parity basis. We discuss the role of different terms in the two-photon dissipator and show how one can extend existing results for the closed Rabi spectrum to the open case. Furthermore, we characterize the role of the Z2 symmetry in the excitation-relaxation dynamics of the system as a function of light-matter coupling. Importantly, we observe that initial states with even-odd parity manifest qualitatively distinct transient and steady state behaviors, contrary to the Hermitian dynamics that is only sensitive to whether or not the initial state is parity invariant. Moreover, the parity-sensitive dynamical behavior is not a creature of ultrastrong coupling and is present even at weak coupling values.

Original language | English (US) |
---|---|

Article number | 043601 |

Journal | Physical review letters |

Volume | 122 |

Issue number | 4 |

DOIs | |

State | Published - 2019 |

## All Science Journal Classification (ASJC) codes

- Physics and Astronomy(all)

## Access to Document

## Fingerprint Dive into the research topics of 'Quantum Rabi Model with Two-Photon Relaxation'. Together they form a unique fingerprint.

## Cite this

- APA
- Author
- BIBTEX
- Harvard
- Standard
- RIS
- Vancouver

*Physical review letters*,

*122*(4), [043601]. https://doi.org/10.1103/PhysRevLett.122.043601

**Quantum Rabi Model with Two-Photon Relaxation**. In: Physical review letters. 2019 ; Vol. 122, No. 4.

}

*Physical review letters*, vol. 122, no. 4, 043601. https://doi.org/10.1103/PhysRevLett.122.043601

**Quantum Rabi Model with Two-Photon Relaxation.** / Malekakhlagh, Moein; Rodriguez, Alejandro W.

Research output: Contribution to journal › Article › peer-review

TY - JOUR

T1 - Quantum Rabi Model with Two-Photon Relaxation

AU - Malekakhlagh, Moein

AU - Rodriguez, Alejandro W.

N1 - Funding Information: We appreciate helpful discussions with Pengning Chao. This work was supported by the National Science Foundation under Grants No. DMR-1454836, No. DMR 1420541, and Award No. EFMA-1640986. Funding Information: Here, we study the dissipative dynamics of the system and discuss the role of Z 2 symmetry. For concreteness, we consider the situation in which the cavity is initially prepared with an even or odd number of photons, and describe the ensuing dynamics of the cavity photon and qubit population as a function of both time and g . In particular, we consider two scenarios of starting with two [ ρ ^ ( 0 ) = | 2 , g ⟩ ⟨ 2 , g | ] or three [ ρ ^ ( 0 ) = | 3 , g ⟩ ⟨ 3 , g | ] initial cavity photons and a qubit in the ground state, as representatives of the plus or minus parity subspaces. Because of pair exchange of photons with the environment, we intuitively expect states with even or odd initial cavity photons to exhibit different transient and steady state behavior. First, consider the simplest case of g = 0 . This choice of parameter decouples the qubit and hence corresponds to the problem of a single cavity mode with two-photon relaxation, which has been studied in detail using multiple methods [36–40] . In this case, initial states having even (odd) numbers of cavity photons end up with zero (one) cavity photons in the steady state [55] . Next, we move on to characterize the interplay of two-photon relaxation and the qubit for g ≠ 0 . Here, closed form analytical solutions of the evolution operator at arbitrary g seem intractable, and instead we employ numerical integration of the Lindblad Eq. (6a) . The time evolution of the cavity or qubit excitations as a function of g is studied in Figs. 3 and 4 for the cases of two and three initial cavity photons, correspondingly. In both cases, it is generally observed that as g is increased, more complex beatings between various normal modes emerge. Such beatings can be approximately understood from the mapping of the initial cavity state to the corresponding eigenmodes of the open Rabi model. This shows which modes are more active at a given value of g in each parity subspace [Figs. 3(c) and 4(c) ]. For example, for the case of ρ ^ ( 0 ) = | 2 , g ⟩ ⟨ 2 , g | , the initial probability is shared between states | 1 g , + ⟩ and | 2 g , + ⟩ up to intermediate values of g ( 0 < g ≲ 0.5 ν c ), beyond which | 1 g , + ⟩ and | 3 g , + ⟩ dominate. The corresponding frequency and decay rate of the modes can be obtained from Figs. 2(a)–2(b) . 3 10.1103/PhysRevLett.122.043601.f3 FIG. 3. Excitation-relaxation dynamics of the system of Fig. 2 when the system is prepared with two cavity photons and the qubit is in the ground state, i.e., ρ ^ ( 0 ) = | 2 , g ⟩ ⟨ 2 , g | = | 2 , + ⟩ ⟨ 2 , + | as a function of light-matter coupling g . (a) Cavity photon number, (b) qubit excitation number, and (c) mapping of the bare state | 2 , g ⟩ to the eigenmodes in the even ( + ) parity subspace. For convenience, we omit the parity index in the x axis. (d) Steady state populations. Model parameters are the same as in Fig. 2 . The time axis in (a) and (b) is normalized to half of the cavity round-trip time T c ≡ π / ν c . The two-photon relaxation time reads T κ 2 ≡ 1 / κ c 2 = 40 T c / π . The cavity mode Hilbert space cutoff is chosen as N c = 9 . 4 10.1103/PhysRevLett.122.043601.f4 FIG. 4. Excitation-relaxation dynamics when the system is prepared with three cavity photons and the qubit is in the ground state, i.e., ρ ^ ( 0 ) = | 3 , g ⟩ ⟨ 3 , g | = | 3 , - ⟩ ⟨ 3 , - | , as a function of light-matter coupling g . The figure follows the same format as Fig. 3 , except that the bare state | 3 , g ⟩ is instead mapped to eigenmodes in the odd ( − ) parity subspace. Other parameters are the same as in Fig. 3 . Despite this generic similarity, it is observed that due to the nontrivial interplay of light-matter coupling and two-photon relaxation, the two cases under consideration have different transient and steady state characteristics. For the case of two initial cavity photons, we observe that the system reaches steady state on a timescale that is more or less given by the two-photon relaxation rate κ c 2 [Figs. 3(a)–3(b) ]. On the other hand, in the case of three initial cavity photons, the transient dynamics has more features. Generally, at small g , the dynamics can be described as follows [Figs. 4(a)–4(b) ]: First, a fast depletion of the initial three cavity photons into one photon, with timescale roughly determined by κ c 2 . This can be seen by the sharp transition of the cavity excitation number from 3 to approximately 1 [red to blue in Fig. 4(a) ]. Second, a slower depletion of the remaining cavity photon after a large number of Rabi exchanges between the qubit and the cavity, with timescale roughly determined by the decay rate of state | 1 g , - ⟩ . Essentially, since two-photon relaxation only allows pairs of exchange with the environment, the quantum state | 1 g , - ⟩ acts like a dark state at g = 0 (i.e., | 1 , g ⟩ ). As g is increased, the decay rate of this state is barely modified up until g / ν c ≈ 0.5 [See Fig. 2(b) ], consistent with the observed long-lived excitations in the qubit or cavity dynamics [Figs. 4(a)–4(b) ]. Steady state excitations have also been studied as a function of g in [Figs. 3(d)–4(d) ]. In the case of two initial photons, we observe that the steady state populations of the cavity and qubit increase nonmonotonically with increasing g , exhibiting a local maximum close to g ≈ κ c 2 . The case of three initial photons is more complicated. For small g < κ c 2 , one observes fast relaxation of two photons, while the remaining photon energy is transferred to the qubit at steady state. At intermediate values of g , the excitation is shared between the cavity and the qubit while at very large g , the qubit excitation saturates and the cavity photon population increases linearly [Fig. 4(d) ]. The overall increase observed in the steady state populations arises from the fact that the coupling in Eq. (6b) appears effectively as an incoherent drive on the cavity. Lastly, we note that the steady state quantities obtained from the Linbdlad formalism will become less accurate at large values of g , as one needs to account for the renormalization of the dissipator arising from the underlying system-bath formalism [56] , resulting in a Bloch-Redfield master equation [57] . Using Rayleigh-Schrödinger perturbation theory, however, one can show that disspator renormalizations are higher order in g compared to the ones for the Hamiltonian. This leaves a window, at intermediate values of coupling, where the use of bare dissipators is still justified. We appreciate helpful discussions with Pengning Chao. This work was supported by the National Science Foundation under Grants No. DMR-1454836, No. DMR 1420541, and Award No. EFMA-1640986. [1] 1 I. I. Rabi , Phys. Rev. 51 , 652 ( 1937 ). PHRVAO 0031-899X 10.1103/PhysRev.51.652 [2] 2 M. Brune , F. Schmidt-Kaler , A. Maali , J. Dreyer , E. Hagley , J. M. Raimond , and S. Haroche , Phys. Rev. Lett. 76 , 1800 ( 1996 ). PRLTAO 0031-9007 10.1103/PhysRevLett.76.1800 [3] 3 J. M. Raimond , M. Brune , and S. Haroche , Rev. Mod. Phys. 73 , 565 ( 2001 ). RMPHAT 0034-6861 10.1103/RevModPhys.73.565 [4] 4 H. Mabuchi and A. Doherty , Science 298 , 1372 ( 2002 ). SCIEAS 0036-8075 10.1126/science.1078446 [5] 5 H. Walther , B. T. Varcoe , B.-G. Englert , and T. Becker , Rep. Prog. Phys. 69 , 1325 ( 2006 ). RPPHAG 0034-4885 10.1088/0034-4885/69/5/R02 [6] 6 A. Wallraff , D. I. Schuster , A. Blais , L. Frunzio , R.-S. Huang , J. Majer , S. Kumar , S. M. Girvin , and R. J. Schoelkopf , Nature (London) 431 , 162 ( 2004 ). NATUAS 0028-0836 10.1038/nature02851 [7] 7 A. Blais , R.-S. Huang , A. Wallraff , S. M. Girvin , and R. J. Schoelkopf , Phys. Rev. A 69 , 062320 ( 2004 ). PLRAAN 1050-2947 10.1103/PhysRevA.69.062320 [8] 8 I. Chiorescu , P. Bertet , K. Semba , Y. Nakamura , C. Harmans , and J. Mooij , Nature (London) 431 , 159 ( 2004 ). NATUAS 0028-0836 10.1038/nature02831 [9] 9 D. I. Schuster , A. A. Houck , J. A. Schreier , A. Wallraff , J. M. Gambetta , A. Blais , L. Frunzio , J. Majer , B. Johnson , M. H. Devoret , S. M. Girvin , and R. J. Schoelkopf , Nature (London) 445 , 515 ( 2007 ). NATUAS 0028-0836 10.1038/nature05461 [10] 10 J. Clarke and F. K. Wilhelm , Nature (London) 453 , 1031 ( 2008 ). NATUAS 0028-0836 10.1038/nature07128 [11] 11 M. Hofheinz , H. Wang , M. Ansmann , R. C. Bialczak , E. Lucero , M. Neeley , A. O’connell , D. Sank , J. Wenner , J. M. Martinis , Nature (London) 459 , 546 ( 2009 ). NATUAS 0028-0836 10.1038/nature08005 [12] 12 E. K. Irish and K. Schwab , Phys. Rev. B 68 , 155311 ( 2003 ). PRBMDO 0163-1829 10.1103/PhysRevB.68.155311 [13] 13 A. N. Cleland and M. R. Geller , Phys. Rev. Lett. 93 , 070501 ( 2004 ). PRLTAO 0031-9007 10.1103/PhysRevLett.93.070501 [14] 14 K. C. Schwab and M. L. Roukes , Phys. Today 58 , No. 7 , 36 ( 2005 ). PHTOAD 0031-9228 10.1063/1.2012461 [15] 15 M. LaHaye , J. Suh , P. Echternach , K. C. Schwab , and M. L. Roukes , Nature (London) 459 , 960 ( 2009 ). NATUAS 0028-0836 10.1038/nature08093 [16] 16 K. Hennessy , A. Badolato , M. Winger , D. Gerace , M. Atatüre , S. Gulde , S. Fält , E. L. Hu , and A. Imamoğlu , Nature (London) 445 , 896 ( 2007 ). NATUAS 0028-0836 10.1038/nature05586 [17] 17 D. Leibfried , R. Blatt , C. Monroe , and D. Wineland , Rev. Mod. Phys. 75 , 281 ( 2003 ). RMPHAT 0034-6861 10.1103/RevModPhys.75.281 [18] 18 J. Pedernales , I. Lizuain , S. Felicetti , G. Romero , L. Lamata , and E. Solano , Sci. Rep. 5 , 15472 ( 2015 ). SRCEC3 2045-2322 10.1038/srep15472 [19] 19 E. T. Jaynes and F. W. Cummings , Proc. IEEE 51 , 89 ( 1963 ). IEEPAD 0018-9219 10.1109/PROC.1963.1664 [20] 20 T. Niemczyk , F. Deppe , H. Huebl , E. Menzel , F. Hocke , M. Schwarz , J. Garcia-Ripoll , D. Zueco , T. Hümmer , E. Solano , Nat. Phys. 6 , 772 ( 2010 ). NPAHAX 1745-2473 10.1038/nphys1730 [21] 21 P. Forn-Díaz , J. Lisenfeld , D. Marcos , J. J. García-Ripoll , E. Solano , C. J. P. M. Harmans , and J. E. Mooij , Phys. Rev. Lett. 105 , 237001 ( 2010 ). PRLTAO 0031-9007 10.1103/PhysRevLett.105.237001 [22] 22 F. Yoshihara , T. Fuse , S. Ashhab , K. Kakuyanagi , S. Saito , and K. Semba , Nat. Phys. 13 , 44 ( 2017 ). NPAHAX 1745-2473 10.1038/nphys3906 [23] 23 E. Irish , Phys. Rev. Lett. 99 , 173601 ( 2007 ). PRLTAO 0031-9007 10.1103/PhysRevLett.99.173601 [24] 24 V. V. Albert , G. D. Scholes , and P. Brumer , Phys. Rev. A 84 , 042110 ( 2011 ). PLRAAN 1050-2947 10.1103/PhysRevA.84.042110 [25] 25 D. Braak , Phys. Rev. Lett. 107 , 100401 ( 2011 ). PRLTAO 0031-9007 10.1103/PhysRevLett.107.100401 [26] 26 Q.-H. Chen , C. Wang , S. He , T. Liu , and K.-L. Wang , Phys. Rev. A 86 , 023822 ( 2012 ). PLRAAN 1050-2947 10.1103/PhysRevA.86.023822 [27] 27 V. Bargmann , Commun. Pure Appl. Math. 14 , 187 ( 1961 ). CPMAMV 0010-3640 10.1002/cpa.3160140303 [28] 28 F. A. Wolf , M. Kollar , and D. Braak , Phys. Rev. A 85 , 053817 ( 2012 ). PLRAAN 1050-2947 10.1103/PhysRevA.85.053817 [29] 29 H. Zhong , Q. Xie , M. T. Batchelor , and C. Lee , J. Phys. A 46 , 415302 ( 2013 ). JPAMB5 1751-8113 10.1088/1751-8113/46/41/415302 [30] 30 F. A. Wolf , F. Vallone , G. Romero , M. Kollar , E. Solano , and D. Braak , Phys. Rev. A 87 , 023835 ( 2013 ). PLRAAN 1050-2947 10.1103/PhysRevA.87.023835 [31] 31 D. Z. Rossatto , C. J. Villas-Bôas , M. Sanz , and E. Solano , Phys. Rev. A 96 , 013849 ( 2017 ). PLRAAN 2469-9926 10.1103/PhysRevA.96.013849 [32] 32 Q. Xie , H. Zhong , M. T. Batchelor , and C. Lee , J. Phys. A 50 , 113001 ( 2017 ). JPAMB5 1751-8113 10.1088/1751-8121/aa5a65 [33] 33 H. Simaan and R. Loudon , J. Phys. A 8 , 539 ( 1975 ). JPHAC5 0305-4470 10.1088/0305-4470/8/4/016 [34] 34 H. Simaan and R. Loudon , J. Phys. A 11 , 435 ( 1978 ). JPHAC5 0305-4470 10.1088/0305-4470/11/2/018 [35] 35 H. Voigt , A. Bandilla , and H.-H. Ritze , Z. Phys. B 36 , 295 ( 1980 ). ZPBBDJ 0340-224X 10.1007/BF01325293 [36] 36 J. Gea-Banacloche , Phys. Rev. Lett. 62 , 1603 ( 1989 ). PRLTAO 0031-9007 10.1103/PhysRevLett.62.1603 [37] 37 L. Gilles and P. L. Knight , Phys. Rev. A 48 , 1582 ( 1993 ). PLRAAN 1050-2947 10.1103/PhysRevA.48.1582 [38] 38 A. Klimov and J. Romero , J. Opt. B 5 , S316 ( 2003 ). JOBOFD 1464-4266 10.1088/1464-4266/5/3/363 [39] 39 A. Voje , A. Croy , and A. Isacsson , New J. Phys. 15 , 053041 ( 2013 ). NJOPFM 1367-2630 10.1088/1367-2630/15/5/053041 [40] 40 V. V. Albert and L. Jiang , Phys. Rev. A 89 , 022118 ( 2014 ). PLRAAN 1050-2947 10.1103/PhysRevA.89.022118 [41] 41 Z. Leghtas , S. Touzard , I. M. Pop , A. Kou , B. Vlastakis , A. Petrenko , K. M. Sliwa , A. Narla , S. Shankar , M. J. Hatridge , Science 347 , 853 ( 2015 ). SCIEAS 0036-8075 10.1126/science.aaa2085 [42] 42 M. Wolinsky and H. J. Carmichael , Phys. Rev. Lett. 60 , 1836 ( 1988 ). PRLTAO 0031-9007 10.1103/PhysRevLett.60.1836 [43] 43 See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevLett.122.043601 for further discussion on physical realization of two-photon relaxation, Z 2 symmetry of the Rabi model, and comparison between the phenomenological and full spectrum, which includes the following additional Refs. [41,42,44–48]. [44] 44 J. Casanova , G. Romero , I. Lizuain , J. J. García-Ripoll , and E. Solano , Phys. Rev. Lett. 105 , 263603 ( 2010 ). PRLTAO 0031-9007 10.1103/PhysRevLett.105.263603 [45] 45 X. Yi and S. Yu , J. Opt. B 3 , 372 ( 2001 ). JOBOFD 1464-4266 10.1088/1464-4266/3/6/304 [46] 46 S. Felicetti , D. Z. Rossatto , E. Rico , E. Solano , and P. Forn-Díaz , Phys. Rev. A 97 , 013851 ( 2018 ). PLRAAN 2469-9926 10.1103/PhysRevA.97.013851 [47] 47 T. Prosen , J. Stat. Mech. ( 2010 ) P07020 . JSMTC6 1742-5468 10.1088/1742-5468/2010/07/P07020 [48] 48 T. G. Kolda and B. W. Bader , SIAM Rev. 51 , 455 ( 2009 ). SIREAD 0036-1445 10.1137/07070111X [49] 49 M. Mirrahimi , Z. Leghtas , V. V. Albert , S. Touzard , R. J. Schoelkopf , L. Jiang , and M. H. Devoret , New J. Phys. 16 , 045014 ( 2014 ). NJOPFM 1367-2630 10.1088/1367-2630/16/4/045014 [50] 50 L. Sun , A. Petrenko , Z. Leghtas , B. Vlastakis , G. Kirchmair , K. Sliwa , A. Narla , M. Hatridge , S. Shankar , J. Blumoff , Nature (London) 511 , 444 ( 2014 ). NATUAS 0028-0836 10.1038/nature13436 [51] 51 N. Ofek , A. Petrenko , R. Heeres , P. Reinhold , Z. Leghtas , B. Vlastakis , Y. Liu , L. Frunzio , S. Girvin , L. Jiang , Nature (London) 536 , 441 ( 2016 ). NATUAS 0028-0836 10.1038/nature18949 [52] 52 C. Wang , Y. Y. Gao , P. Reinhold , R. W. Heeres , N. Ofek , K. Chou , C. Axline , M. Reagor , J. Blumoff , K. Sliwa , Science 352 , 1087 ( 2016 ). SCIEAS 0036-8075 10.1126/science.aaf2941 [53] 53 G. Lindblad , Commun. Math. Phys. 48 , 119 ( 1976 ). CMPHAY 0010-3616 10.1007/BF01608499 [54] 54 V. Gorini , A. Kossakowski , and E. C. G. Sudarshan , J. Math. Phys. (N.Y.) 17 , 821 ( 1976 ). JMAPAQ 0022-2488 10.1063/1.522979 [55] Note that states with mixed parity can always be expressed as linear combinations of two independent problems. [56] 56 F. Beaudoin , J. M. Gambetta , and A. Blais , Phys. Rev. A 84 , 043832 ( 2011 ). PLRAAN 1050-2947 10.1103/PhysRevA.84.043832 [57] 57 A. Redfield , Adv. Magn. Opt. Reson. 1 , 1 ( 1965 ). 10.1016/B978-1-4832-3114-3.50007-6

PY - 2019

Y1 - 2019

N2 - We study a cavity-QED setup consisting of a two-level system coupled to a single cavity mode with two-photon relaxation. The system dynamics is modeled via a Lindblad master equation consisting of the Rabi Hamiltonian and a two-photon dissipator. We show that an even-photon relaxation preserves the Z2 symmetry of the Rabi model, and provide a framework to study the corresponding non-Hermitian dynamics in the number-parity basis. We discuss the role of different terms in the two-photon dissipator and show how one can extend existing results for the closed Rabi spectrum to the open case. Furthermore, we characterize the role of the Z2 symmetry in the excitation-relaxation dynamics of the system as a function of light-matter coupling. Importantly, we observe that initial states with even-odd parity manifest qualitatively distinct transient and steady state behaviors, contrary to the Hermitian dynamics that is only sensitive to whether or not the initial state is parity invariant. Moreover, the parity-sensitive dynamical behavior is not a creature of ultrastrong coupling and is present even at weak coupling values.

AB - We study a cavity-QED setup consisting of a two-level system coupled to a single cavity mode with two-photon relaxation. The system dynamics is modeled via a Lindblad master equation consisting of the Rabi Hamiltonian and a two-photon dissipator. We show that an even-photon relaxation preserves the Z2 symmetry of the Rabi model, and provide a framework to study the corresponding non-Hermitian dynamics in the number-parity basis. We discuss the role of different terms in the two-photon dissipator and show how one can extend existing results for the closed Rabi spectrum to the open case. Furthermore, we characterize the role of the Z2 symmetry in the excitation-relaxation dynamics of the system as a function of light-matter coupling. Importantly, we observe that initial states with even-odd parity manifest qualitatively distinct transient and steady state behaviors, contrary to the Hermitian dynamics that is only sensitive to whether or not the initial state is parity invariant. Moreover, the parity-sensitive dynamical behavior is not a creature of ultrastrong coupling and is present even at weak coupling values.

UR - http://www.scopus.com/inward/record.url?scp=85061003405&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85061003405&partnerID=8YFLogxK

U2 - 10.1103/PhysRevLett.122.043601

DO - 10.1103/PhysRevLett.122.043601

M3 - Article

C2 - 30768294

AN - SCOPUS:85061003405

VL - 122

JO - Physical Review Letters

JF - Physical Review Letters

SN - 0031-9007

IS - 4

M1 - 043601

ER -