Arxiv:2004.05826V2 [Quant-Ph] 12 Mar 2021 Hermitian Systems [15–22], Time Modulated Systems [23– Periment Difﬁculties for Large Lattices

Interaction induced non-reciprocal three-level quantum transport

Sai Li,1 Tao Chen,1 Jia Liu,1, 2, ∗ and Zheng-Yuan Xue1, 2, † 1Guangdong Provincial Key Laboratory of Quantum Engineering and Quantum Materials, and School of Physics and Telecommunication Engineering, South China Normal University, Guangzhou 510006, China 2Frontier Research Institute for Physics, South China Normal University, Guangzhou 510006, China (Dated: March 15, 2021) Besides its fundamental importance, non-reciprocity has also found many potential applications in quantum technology. Recently, many quantum systems have been proposed to realize non-reciprocity, but stable non- reciprocal process is still experimentally difficult in general, due to the needed cyclical interactions in artificial systems or operational difficulties in solid state materials. Here, we propose a new kind of interaction induced non-reciprocal operation, based on the conventional STIRAP setup, which removes the experimental difficulty of requiring cyclical interaction, and thus it is directly implementable in various quantum systems. Furthermore, we also illustrate our proposal on a chain of three coupled superconducting transmons, which can lead to a non-reciprocal circulator with high fidelity without a ring coupling configuration as in the previous schemes or implementations. Therefore, our protocol provides a promising way to explore fundamental non-reciprocal quantum physics as well as realize non-reciprocal quantum device.

PACS numbers: 03.65.Ta, 03.65.Aa, 03.67.Lx Keywords: Non-reciprocity, quantum transport, superconducting quantum circuits

I. INTRODUCTION on superconducting circuits, they usually need cyclical interaction among at least three levels of a superconducting qubit Reciprocity, which means that the measured scattering does device [35] or three coupled superconducting qubit devices in not change when the source and the detector are interchanged a two-dimensional configuration [33, 34]. These are experi- [1], is a fundamental phenomenon in both classical and quan- mentally challenging for large scale lattices, as they require tum regimes. Meanwhile, non-reciprocal devices, such as that quantum systems have cyclical transition or at least two- isolators and circulators, are also essential in both classical dimensional configuration for three-level non-reciprocal pro- and quantum information processing. Especially, circulators cess. Therefore, proposals using tunable non-reciprocal pro- can separate opposite signal flows, spanning from classical to cess and its potential applications to achieve non-reciprocity quantum computation and communication systems [2]. Thus, are still highly desired theoretically and experimentally. circulators are vital for the design of full-duplex communica- Here, we propose a general scheme on a three-level quan- tion systems, which can transmit and receive signals through tum system based on the conventional STIRAP setup [39, 40] a same frequency channel, providing the opportunity to en- to realize non-reciprocal operations by time modulation. The hance channel capacity and reduce power consumption [3]. distinct merit of our proposal is that the realization only Therefore, many theoretical and experimental progresses have needs two time-modulation coupling, which removes the ex- been made recently to build non-reciprocal devices in differ- perimental difficulty of requiring cyclical interaction, and ent quantum systems. Specifically, the conventional way of thus it is directly implementable in various quantum sys- realizing non-reciprocity is achieved by adding magnetic field tems, for example, superconducting quantum circuits sys- or using magnetic materials directly [4]. However, the exter- tems [2, 33, 41, 42], nuclear magnetic resonance systems nal magnetic field would affect the transformation and mag- [43], a nitrogen-vacancy center in diamonds [14], trapped ions netic materials is hardly to induce non-reciprocity. [44], hybrid quantum systems [45], and so on. Meanwhile, Recently, realization of non-reciprocity has been proposed the three-level non-reciprocal process can be implemented in in many artificial quantum systems, such as in non-linear a one-dimensional configuration instead of two-dimensional systems [5–7], synthetic magnetism systems [8–14], non- system in previous proposals, which greatly releases the ex- arXiv:2004.05826v2 [quant-ph] 12 Mar 2021 Hermitian systems [15–22], time modulated systems [23– periment difficulties for large lattices. 32, 34, 35], etc. Although these successful methods have It is well known that superconducting quantum circuits sys- been realized in nitrogen-vacancy centers systems [14], cold tem is scalable and controllable, and thus attracts great at- atom systems [22], superconducting circuits [34, 35], and op- tention in many researches. Different from the cold atoms tical systems [36–38], tunable non-reciprocal process is still and optical lattice systems, superconducting circuits possess lacking for quantum manipulation. This is because previous good individual controllability and easy scalability. Follow- schemes rely highly on special material properties, e.g., non- ing that, we illustrate our proposal on a chain of three coupled linear property. In addition, for the experimental implementa- superconducting transmon devices with appropriate parame- tions of non-reciprocity induced by synthetic magnetism, e.g., ters, and achieve a non-reciprocal circulator with high fidelity. As our proposal is based on a one-dimension superconducting lattice, e.g. Refs. [46, 47], and with demonstrated techniques ∗ [email protected] there, thus it can be directly verified. Therefore, our proposal † [email protected] provides a new approach based on time modulation for engi- 2

receive the quantum information from sender |Bi and send in- (a) | A  formation to receiver |Ai. Those quantum processes are very |B important for quantum information transformation and processing. g (t) 1 g 2 (t) III. CONSTRUCTION |M Here, we proceed to introduce details for construction of non-reciprocal scattering matrix. Generally, the evolution pro- (b) (c) cess of time-dependent Hamiltonian with period τ can be ex- pressed as U(τ) = T exp[−i R τ H(t)dt], where T is time- |B S1 0 order operator. When the quantum system satisfies the von- | A R2 Neumann equation ∂I(t)/∂t+i[H(t),I(t)] = 0 with dynam-   ical invariant I(t), the concrete solution of the evolution process with time-dependent Hamiltonian H(t) can be carried out by Lewis-Riesenfeld (LR) invariant method [48]. Here, we R1 |M S2 concentrate on the solution of nonreciprocal transition evolution process with time modulation. With the LR method, the FIG. 1. Schematic diagram. (a) Initial picture for generating non- evolution process with period τ can be exactly derived as reciprocal devices: two subspaces |Ai and |Bi are resonantly cou- X −iθn(τ) pled to subspace |Mi simultaneously with time modulated coeffi- U(τ) = e |µn(τ)ihµn(0)|, (2) cients g1(t) and g2(t). (b) Quantum circulator with one direction n=0,± flow through a period τ. (c) |Mi simultaneously serve as receiver R1 and sender S2 for receiving information from sender S1 and send where in the Hilbert space {|Ai, |Mi, |Bi}, information to receiver R2 through a period τ.  cos γ(t) cos β(t)  |µ0(t)i =  −i sin γ(t)  (3) − cos γ(t) sin β(t) neering non-reciprocal devices, which can find many interesting applications in quantum information processing, including and one-way propagation of quantum information, quantum mea-   surement and readout, and quantum steering. 1 sin γ(t) cos β(t) ± i sin β(t) |µ±(t)i = √  i cos γ(t)  (4) 2 − sin γ(t) sin β(t) ± i cos β(t) II. GENERAL FRAMEWORK are the eigenstates of the invariant [49] Now, we start from a general three-level quantum system   labeled in the Hilbert space {|Ai, |Mi, |Bi}. As shown in µ 0 cos γ sin β −i sin γ I(t) = cos γ sin β 0 cos γ cos β , (5) Fig. 1(a), considering |Ai and |Bi simultaneously coupled to 2   |Mi resonantly. Assuming ¯h = 1 hereafter, the interaction i sin γ cos γ cos β 0 Hamiltonian in the interaction picture can be written as where µ is an arbitrary constant with unit of frequency to H(t) = g1(t)|AihM| + g2(t)|BihM| + H.c., (1) keep I(t) with dimensions of energy, γ(t) and β(t) are auxiliary parameters , which satisfy the von-Neumann equation g (t) where 1,2 are the time-modulation coupling strength. ∂I(t)/∂t + i[H(t),I(t)] = 0, and θ (τ) is the LR phase with H(t) n Based on the Hamiltonian , even if there is no direct cou- θ (τ) = 0 and θ (τ) = −θ (τ), which can be addressed by |Ai |Bi 0 − + pling between bare states and , these two bare states auxiliary parameters γ(t) and β(t). To induce non-reciprocal |Mi are both coupled to state , then, the transition between transition evolution process, we set the boundary conditions |Ai |Bi bare states and can be realized via the middle state as |Mi, i.e. STIRAP [39, 40]. Especially, when the pulse shapes of g1(t) and g2(t) are different from each other, the symme- γ(0) = γ(τ) = 0, β(0) = 0, β(τ) = π/2. (6) try of the system can be broken naturally. Thus, with appropriate designed time-modulation coupling strengths g1,2(t), a After that, the final evolution operator in the Hilbert space quantum circulator with one-direction flow can be achieved {|Ai, |Mi, |Bi} can be determined as through a period evolution with time τ, i.e., |Ai → |Bi→   |Mi → |Ai illustrated in in Fig. 1(b). That means, on the one 0 −i sin θ+(τ) cos θ+(τ) hand, transition |Ai → |Bi is allowed, not vice versa. Mean- U[θ+(τ)] =  0 cos θ+(τ) −i sin θ+(τ)  . (7) while, the process |Bi → |Mi and |Mi → |Ai can be simul- −1 0 0 taneously realized, which means that |Mi can simultaneously 3

(a) H (t) = H (t) + H (t), where H (t) and H (t) are free  |1 T f int f int |1 |1 and interaction Hamiltonian respectively. For the free part, X Hf (t) = {[ωj + f(j(t))]|1ijh1|} + ωM |1iM h1|, (8) |0 |0 |0 j=A,B A M B where (t) = f −1(F˙ (t)) with F (t) = η(t) sin(νt). For the interaction term, X (b) | 100  Hint(t) = [ gj(|0ijh1| + H.c.)] · (|0iM h1| + H.c.). (9)  | 001 j=A,B

Transforming to the rotating frame defined by V = V1 + P V2, where V1 = exp[i( j=A,B ωj|1ijh1| + ωM |1iM h1|)] P | 010  and V2 = exp[i j=A,B Fj(t)|1ijh1|], and the transformed Hamiltonian is FIG. 2. Illustration of our scheme with three transmon devices. (a) dV † H (t) = V †H (t)V + i V Two qubit-frequency driven transmons A and B with the respective t T dt longitudinal field A,B (t) resonantly coupled to the transmon M. X −iωj t iFj (t) (b) Effective resonant coupling architecture in the single-excitation = gj[|0ijh1|e e + H.c.] (10) subspace {|100i, |010i, |001i}. j=A,B −iωM t ⊗ [|0iM h1|e + H.c.], in the single-excitation subspace {|100i, |010i, |001i}, where |ambi ≡ |ai ⊗|mi ⊗|bi labels the product states of three To understand the result clearly, for the case θ+(τ) = π, the A M B final evolution operator U[π] = −|MihM|−|AihB|−|BihA| transmons, after neglecting the high order oscillating terms, represents a normal two-direction transitions. Especially, for the Hamiltonian can be written as another case θ (τ) = 3π/2, the evolution process U[3π/2] = i∆At−iFA(t) + HI (t) =gA|100ih010|e + i|AihM| + i|MihB| − |BihA| shows non-reciprocal transi- i∆B t−iFB (t) tions, that means transition |Ai → |Bi is allowed and transi- gB|001ih010|e + H.c., (11) tion |Bi → |Ai is forbidden for the same process. Obviously, where gj is coupling strength for transmons A, B to M the evolution operator U[3π/2] also means a cyclic transporta- and ∆j = ωj − ωM are the frequency difference. Then, tion. To sum up, the evolution operator U[3π/2] induces a using the Jacobi-Anger identity expansion exp(iFj(t)) = cyclic chiral transportation, which exactly realizes a quantum P∞ m π m=−∞ i Jm[ηj(t)] exp[im(νjt− 2 )], and considering the circulator, from pure time-modulation of the interaction. resonant interaction case ∆j = νj, the effective Hamiltonian with time modulation can be obtained as g0 (t) g0 (t) IV. ILLUSTRATIVE SCHEME WITH TRANSMONS H (t) = A |100ih010| + B |001ih010| + H.c., (12) eff 2 2 0 Here, we propose a scheme on superconducting quantum where gj(t) = 2gjJ1(ηj(t)) are effective time-modulation circuits. For a transmon [51], there are three lowest levels, coupling strength for transmons A, B to M with J1 being the which can be resonantly driven by two microwave fields to Bessel function. We can use the effective Hamiltonian Heff(t) induce the Hamiltonian H(t) in Eq. (1) our proposal. In this to realize quantum circulator in our protocol. case, only nonreciprocal state transfer within a transmon can With the LR invariant method [49], according to the von- Neumann equation ∂I(t)/∂t + i[Heff (t),I(t)] = 0, the form be obtained, and the implementation is straightforward, i.e., 0 letting |0i, |1i, and |2i take the role of |Ai, |Mi, and |Bi. of gj(t) can be given as Furthermore, we consider a more interesting case, that g0 (t) = 2[β˙(t) cot γ(t) sin β(t) +γ ˙ (t) cos β(t)], is, three coupled transmons implementation, with the low- A 0 ˙ est two levels |0i and |1i in superconducting quantum cir- gB(t) = 2[β(t) cot γ(t) cos β(t) − γ˙ (t) sin β(t)]. (13) cuits. As shown in Fig. 2(a), we labeled three transmons with Considering the boundary conditions Eq. (6), the commuta- A, M and B with frequencies ωA,M,B and anharmonicities tion relations [H(0),I(0)] = [H(τ),I(τ)] = 0, and the ex- αA,M,B. Here, we introduce qubit frequency drives f((t)) 0 perimental apparatus restriction, the values gj(t) can be set as [52], which can be determined experimentally by the longi- zeros at time t = 0 and τ, thus, a set of auxiliary parameters −1 ˙ tudinal field (t) = f (F (t)), where F (t) = η(t) sin(νt) γ(t) and β(t) can be selected in a proper form [50] as is intentionally chosen with ν being the frequency of the lon- λ gitudinal field (t), and two qubit-frequency drives f(j(t)) γ(t) = t2(t − τ)2, (j = A, B) are added in transmons A and B respectively (τ/2)4 to induce time-modulation resonant interaction with trans- −10πt7 35πt6 42πt5 35πt4 β(t) = + − + , (14) mon M. Then, the coupled system can be described by τ 7 τ 6 τ 5 2τ 4 4

15 1 1 (a) 12 0.8 (b) 0.8 (c) 9 0.6 0.6 / + 6 0.4 g (t)/g 0.4

A Amax Fidelity 3 0.2 g (t)/g 0.2 B Bmax 0 0 0 1230 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 t/ t/ 1 (d) (e) (f) 0.8 |100 0.6 |010 0.4 |001 Fidelity F 0.2 s 0 0 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 t/ t/ t/

FIG. 3. Numerical performance. (a) The LR phase θ+(τ) with respect to auxiliary parameter λ, where black square represents θ+(τ) = 3π/2. 0 (b) The pulse shapes of asymmetrical time-modulation coupling strength gA,B (t) with θ+(τ) = 3π/2 in a period τ. (c) The ﬁdelity of the quantum circulator U[3π/2] for simultaneously sending and receiving the quantum information in a period τ. (d), (e) and (f) The state populations and ﬁdelity of the quantum circulator U[3π/2] in a period τ with the initial state being |100i, |001i and |010i, respectively.

ˆ ˆ ˆ† ˆ† ˆ ˆ† ˆ where λ is a tunable time-independent auxiliary parameter, L(Ok) = Okρ1Ok − OkOkρ1/2 − ρ1OkOk/2 is the Lind- which directly determines the LR phase θ+(τ) concerned in bladian of the operator Oˆk with Oˆk = |0ikh1| + (|0ikh0| − our proposal shown in Fig. 3(a). Furthermore, the effective |1ikh1|). We first evaluate the quantum circulator U[3π/2] 0 coupling strength gj(t) can be carried out according to Eq. with initial state |ψI1i sequentially prepared on the states (13). Then, we realize the final evolution operator U[θ+(τ)]. |100i, |001i and |010i. Define the state fidelity Fs = Following, we choose appropriate experimental parameters hψT 1|ρ|ψT 1i , where ρ is the density matrix and sequential [53] and show how to realize our protocol to achieve non- target states |ψT 1i are −|001i, i|010i and i|100i. As shown reciprocal operations on superconducting quantum circuits. in Fig. 3 (d), (e) and (f), the sequential state fidelity Fs are Due to the anharmonicity of transmon qubits are relatively 99.08%, 99.25% and 99.28%, thus, the results demonstrate small, thus the second excited state will contribute harmfully the construction of our quantum circulator. In our simulation, for the quantum process. To numerically quantify this effect, we do not include control errors of the external driving field, we set the anharmonicity of three transmons as αA = 2π×220 due to the following. First, they can be well-controlled exper- MHz, αM = 2π × 210 MHz and αB = 2π × 230 MHz. imentally. Second, our protocol is robust against the control Meanwhile, we set the frequency of the longitudinal field errors, as show in Fig. 3 (c) (d), (e) and (f), i.e., the target state νj equal to the corresponding frequency difference ∆j as populations are nearly flat around the final times. νA = ∆A = 2π × 345 MHz and νB = ∆B = 2π × 345 MHz respectively to induce time-modulation resonant inter- Then, we evaluate the ability to simultaneously send and re- action in the single-excitation subspace. Furthermore, we set ceive the quantum states for the quantum circulator U[3π/2]. the decoherence rates of the transmons as ΓA = 2π × 3 kHz, We set transmon M as both sender and receiver, which can ΓM = 2π × 4 kHz and ΓB = 2π × 5 kHz, coupling strength send quantum information to transmon A and simultane- for transmons A, B to M as gA = gB = 2π × 10 MHz and ously receive different quantum information from transmon the quantum evolution period τ = 145 ns. Then, to realize B. To demonstrate this, we set the initial states |ψI2i = the quantum circulator U[3π/2], we modify auxiliary param- cos ϑ|010i + sin ϑ|001i, and the corresponding target states eter λ = 0.4974 to make θ+(τ) = 3π/2 and naturally de- are |ψT 2i = i cos ϑ|100i + i sin ϑ|010i. We define the fi- g0 (t) 1 R 2π termine the time-modulation coupling strength j , whose delity Fm = 2π 0 hψT 2|ρ|ψT 2idϑ with the integration nu- pulse shapes are plotted in Fig. 3(b), which is smooth and merically performed for 1001 input states with ϑ being uni- easily experimentally realized. formly distributed over [0, 2π]. As shown in Fig. 3(c), the ini- We numerically simulate the performance of the quantum tial fidelity is less than 0.2 due to the big deference between circulator U[3π/2] by using Lindblad master equation as initial states and target states; the non-monotonic behavior in the time evolution of the fidelity means that the evolution X ρ˙ = i [ρ, HI (t)] + ΓkL Oˆk , (15) process is a non-reciprocal process; and we get the fidelity k=A,M,B Fm = 99.23%, which shows the power of our scheme can realize an effective quantum transfer station which can receive where ρ is the density matrix of the considered system and and send different quantum information in one direction; the 5 infidelity is caused by decoherence about 0.52% and leakage transmon. Considering the scalability and controllability of error about 0.25%. the superconducting quantum circuits, our scheme provides promising candidates for non-reciprocal quantum information processing and devices in the near future.

V. CONCLUSION ACKNOWLEDGEMENTS In summary, we propose a general scheme based on time modulation to realize non-reciprocal operations. Our proposal This work was supported by the National Natural Sci- can be easily realized in many quantum systems. We illus- ence Foundation of China (Grant No. 11874156 and trate our proposal on superconducting quantum circuits with No. 11904111) and the Project funded by China Postdoctoral two driving transmons simultaneously coupled to the middle Science Foundation (Grant No. 2019M652684).

[1] Deak´ L and Ful¨ op¨ T 2012 Anna. Phys. 327 1050 [28] Fleury R, Khanikaev A and Alu` A 2016 Nat. Commun. 7 11744 [2] Devoret M H and Schoelkopf R J 2013 Science 339 1169 [29] Wang D W, Zhou H T, Guo M J, Zhang J X, Evers J and Zhu S [3] Bharadia D, McMilin E and Katti S 2013 ACM SIGCOMM Y 2013 Phys. Rev. Lett. 110 093901 Comp. Commun. Rev. 43 375 [30] Estep N A, Sounas D L, Soric J and Alu` A 2014 Nat. Phys. 10 [4] Aplet L J and Carson J W 1964 Appl. Opt. 3 544 923 [5] Fan L, Wang J, Varghese L T, Shen H, Niu B, Xuan Y, Weiner [31] Reiskarimian N and Krishnaswamy H 2016 Nat. Commun. 7 A M and Qi M 2012 Science 335 447 11217 [6] Kamal A, Roy A, Clarke J and Devoret M H 2014 Phys. Rev. [32] Sounas D L and Alu` A 2017 Nat. Photonics 11 774 Lett. 113 247003 [33] Gu X, Kockum A F, Miranowicz A, Liu Y X and Nori F 2017 [7] Kamal A and Metelmann A 2017 Phys. Rev. Appl. 7 034031 Phys. Rep. 718 1 [8] Koch J, Houck A A, Hur K L and S M Girvin 2010 Phys. Rev. [34] Roushan P et al., 2017 Nat. Phys. 13 146 A 82 043811 [35] Vepsal¨ ainen¨ A, Danilin S and Paraoanu G S 2019 Sci. Adv. 5 [9] Umucalılar R O and Carusotto I 2011 Phys. Rev. A 84 043804 eaau5999 [10] Habraken S J M, Stannigel K, Lukin M D, Zoller P and P Rabl [36] Shen Z, Zhang Y L, Chen Y, Zou C L, Xiao Y F, Zhou X B, 2012 New J. Phys. 14 115004 Sun F W, Guo G C and Dong C H 2016 Nat. Photonics 10 657 [11] Wang Y P, Wang W, Xue Z Y, Yang W L, Hu Y and Wu Y 2015 [37] Ruesink F, Miri M A, Alu´ A and Verhagen E 2016 Nat. Com- Sci. Rep. 5 8352 mun. 7 13662 [12] Sun F X, Mao D, Dai Y T, Ficek Z, He Q Y and Gong Q H 2017 [38] Xu H, Jiang L Y, Clerk A A and Harris J G E 2019 Nature New J. Phys. 19 123039 (London) 568 65 [13] Seif A, DeGottardi W, Esfarjani K and Hafezi M 2018 Nat. [39] Kumar K S, Vepsal¨ ainen¨ A, Danilin S and Paraoanu G S 2016 Commun. 9 1207 Nat. Commun. 7 10628 [14] Barfuss A, Kolbl¨ J, Thiel L, Teissier J, Kasperczyk M and [40] Xu H K, Song C, Liu W Y, Xue G M, Su F F, Deng H, Tian Y, Maletinsky P 2018 Nat. Phys. 14 1087 Zheng D N, Han S Y, Zhong Y P, Wang H, Liu Y X and Zhao [15] Regensburger A, Bersch C, Miri M A, Onishchukov G, S P 2016 Nat. Commun. 7 11018 Christodoulides D N and Peschel U 2012 Nature (London) 488 [41] You J Q and Nori F 2011 Nature (London) 474 589 167 [42] Xiang Z L, Ashhab S, You J Q and Nori F 2013 Rev. Mod. Phys. [16] Fleury R, Sounas D and Alu` A 2015 Nat. Commun. 6 5905 85 623 [17] Chong Y D, Ge L and Stone A D 2011 Phys. Rev. Lett. 106 [43] Feng G, Xu G and Long G 2013 Phys. Rev. Lett. 110 190501 093902 [44] Leibfried D, Blatt R, Monroe C and Wineland D 2003 Rev. [18] Miao P, Zhang Z, Sun J, Walasik W, Longhi S, Litchinitser N Mod. Phys. 75 281 M and Feng L 2016 Science 353 464 [45] Xiang Z L, Ashhab S, You J Q and Nori F 2013 Rev. Mod. Phys. [19] Peng B, Ozdemir¨ S¸K, Lei F, Moniﬁ F, Gianfreda M, Long G L, 85 623 Fan S, Nori F, Bender C M and Yang L 2014 Nat. Phys. 10 394 [46] Li X, Ma Y, Han J, Chen T, Xu Y, Cai W, Wang H, Song Y P, [20] Chang L, Jiang X, Hua S, Yang C, Wen J, Jiang L, Li G, Wang Xue Z Y, Yin Z Q and Sun L 2018 Phys. Rev. Appl. 10 054009 G and Xiao M 2014 Nat. Photonics 8 524 [47] Cai W, Han J, Mei F, Xu Y, Ma Y, Li X, Wang H, Song Y P, [21] Xu X W, Zhao Y J, Wang H, Chen A X and Liu Y X 2019 arxiv: Xue Z Y, Yin Z Q, Jia S and Sun L 2019 Phys. Rev. Lett. 123 1908.08323 080501 [22] Guo W, Chen T, Xie D Z, Xiao T, Deng T S, Gadway B, Yi W [48] Lewis H R and Riesenfeld W B 1969 J. Math. Phys. 10 1458 and Yang B 2020 Phys. Rev. Lett. 124 070402 [49] Chen X, Lizuain I, Ruschhaupt A, Guery-Odelin´ D and Muga J [23] Fang K, Yu Z and Fan S 2012 Phys. Rev. Lett. 108 153901 G 2010 Phys. Rev. Lett. 105 123003 [24] Hafezi M and Rabl P 2012 Opt. Express 20 7672 [50] Zhou J, Li S, Chen T and Xue Z Y 2019 Ann. Phys. (Berlin) [25] Lira H, Yu Z, Fan S and Lipson M 2012 Phys. Rev. Lett. 109 531 1800402 033901 [51] Koch J, Yu T M, Gambetta J, Houck A A, Schuster D I, Majer [26] Poulton C G, Pant R, Byrnes A, Fan S, Steel M J and Eggleton J, Blais A, Devoret M H, Girvin S M and Schoelkopf R J 2007 B J 2012 Opt. Express 20 21235 Phys. Rev. A 76 042319 [27] Sounas D L, Caloz C and Alu` A 2013 Nat. Commun. 4 2407 6

[52] Chu J, Li D Y, Yang X P, Song S Q, Han Z K, Yang Z, Dong Y [53] Barends R et al., 2014 Nature (London) 508 500 Q, Zheng W, Wang Z M, Yu X M, Lan D, Tan X S and Yu Y 2020 Phys. Rev. Appl. 13 064012