Fast topological pumping for the generation of large-scale Greenberger−Horne−Zeilinger states in a superconducting circuit

Abstract

Topological pumping of edge states in the finite lattice with nontrivial topological phases provides a powerful means for robust excitation transfer, requiring extremely slow evolution to follow an adiabatic transfer. Here, we propose fast topological pumping via edge channels to generate large-scale Greenberger−Horne−Zeilinger (GHZ) states in a topological superconducting circuit with a sped-up evolution process. The scheme indicates a conceptual way of designing fast topological pumping related to the instantaneous energy spectrum characteristics rather than relying on the shortcuts to adiabaticity. Based on fast topological pumping, large-scale GHZ states show greater robustness against on-site potential defects, the fluctuation of couplings and losses of the system in comparison with the conventional adiabatic topological pumping. With experimentally feasible qutrit-resonator coupling strengths and moderate decay rates of qutrits and resonators, fast topological pumping drastically improves the scalability of GHZ states with a high fidelity. Our work opens up prospects for the realization of large-scale GHZ states based on fast topological pumping in the superconducting quantum circuit system, which provides potential applications of topological matters in quantum information processing.

Graphical abstract

Keywords

topological pumping / superconducting ciruit / large-scale / Greenberger−Horne−Zeilinger states

1 Introduction

Topological phenomena, rooted in global geometric properties characterized by topological invariants, provide natural protection against perturbation and disorder in systems [1,2]. Topological quantum pumping serves as a transport phenomenon making use of topological properties of a modulated system, including two types: Thouless pumping via topological bulk bands [3,4] and adiabatic pumping via topological edge modes [5-11]. Thouless pumping is a quantized transport directly connected to the Chern number [4]. In the adiabatic pumping, the edge state at one side will be adiabatically coupled to bulk states and then transferred to the other side [5-11]. On account of a topological nature, the transport scheme demonstrates to be robust against the disorder, exploiting topological pumping greatly concerned in applications such as quantum state transfer [8,12,13], quantum information processing [14] and quantum entanglement [15].

A very promising direction towards the realization of an efficient platform to perform fault-tolerant quantum computation combines topological states of matter with superconducting circuits [16-19]. One of the most appealing properties of topological systems is that they host edge states, which are robust to different sources of quantum decoherence. Superconducting quantum circuits with high coherence and straightforward connectivity can be fabricated with well-designed parameters and be accessible to manufacture tens of qubits with various types of couplings [20-25], which enables to implement large-scale quantum simulation [26-30]. Recently, such state-of-the-art techniques enable superconducting quantum circuits to be a promising candidate for investigating topological phenomena theoretically and experimentally [31-47]. Topical examples include analogies to topological quantum walks [39], Zak phases [40], topological phase transitions [41,42], Chern numbers [43,44], topological Uhlmann phase [47] and topological magnon insulator states [46] realizing in the phase space of microwave resonators and superconducting qubit, respectively. Furthermore, the topological robustness of anyonic braiding statistics has also been demonstrated in the superconducting quantum circuit [45].

So far, many topological states and effects have been implemented for quantum entanglement, which is pivotal for quantum communication, quantum information process and quantum metrology [48-50]. It has demonstrated the topological protected entangled N00N state and its transport via edge states in Refs. [51-53]. The adiabatic topological pumping based on edge channels offered a topologically protected way to observe the Hong−Ou−Mandel interference [54] of a photon pair and generate spatial entangled states [9]. Additionally, large-scale Greenberger−Horne−Zeilinger (GHZ) states have been realized by adiabatic topological pumping via edge channels in a superconducting quantum system [15]. Nevertheless, the speed of procedures is extremely slow for a large-scale system, which is intrinsically limited by adiabatic requirements to avoid nonadiabatic transitions between the edge channel and the bulk bands. In order to refrain from the decoherence effects in the system, several proposals have been proposed to accelerate the process of quantum state transfer in topologically non-trivial systems [11-13,55,56]. The fast and robust quantum state transfer can be realized by performing shortcut-to-adiabaticity techniques [13]. A scheme of topological pumping based on adiabatic passage has been proposed by interfacing two dimerized Su−Schrieffer−Heeger (SSH) chains with different topological order [11,56]. In the topological protected quantum state transfer, optimization techniques have been controlled to reach high-fidelity values for short evolution times [12,55]. However, it is also urgent to improve the transfer efficiency in adiabatic topological pumping for generating large-scale entangled states, thereby making sense to enhance the robustness against decoherence effects [57-59] and random fluctuations from quantum system [60-62], and improve the feasibility and scalability of experimental implementations in quantum information processing [63-65].

GHZ state is a typical type of maximally entangled states, which is not only of great interest for fundamental tests of quantum mechanics [66], but also has applications in quantum communications [67,68], error-correction protocols [69], quantum metrology [70] and high-precision spectroscopy [71]. The capability of entangling multiple particles is central to fundamental tests of quantum theory and the prerequisite for quantum information processing. During the past years, theoretical proposals of GHZ states have been investigated in circuit-QED system and other physical platform [15,72-83]. Moreover, experimental realizations of entanglement of three superconducting qubits [84], and GHZ states with three superconducting qubits in circuit-QED system [85], five superconducting qubits via capacitance coupling [86], ten qubits connecting to a bus resonator in a superconducting circuit [87] and eighteen qubits in a 20-qubit superconducting device [25] have been reported. Furthermore, high-dimensional GHZ state with three superconducting transmon qutrits has also generated by the fidelity of 76\mathrm{\%} \pm 1\mathrm{\%} in a superconducting quantum processor [88]. However, the fidelity of GHZ states is limited by the ubiquitous noise and device imperfections based on the superconducting circuit system, indicating that the effectiveness of theoretical proposals and experimental realizations remains to be solved.

In this paper, we present a scheme of fast topological pumping via edge channels with respect to a single-excitation quantum state for generating large-scale GHZ states in a generalized SSH model of superconducting qutrit-resonator circuit. We analytically derive that such a generalized SSH model can be expressed as a two-band structure in momentum space, and its topologically trivial and nontrivial phases can be characterized by the Zak phases equal to 0 and \pi, respectively. The odd-size generalized SSH model is always in the topological non-trivial case, proved by the topologically protected zero-mode edge state in the energy spectrum. The scheme indicates a conceptual way of designing fast topological pumping depending on the instantaneous energy spectrum characteristics for speeding up the generation of large-scale GHZ states, which abides by the adiabatic passage. Furthermore, large-scale GHZ states are robust against on-site potential defects, the fluctuation of couplings and losses of the system which are benefit from the chiral symmetry and topological property of system. Our work may facilitate potential applications of optimization techniques to control topological superconducting circuit system in quantum information processing, due to the following advantages and interests. Firstly, this work focuses on crucial aspects of the energy spectrum to design topological pumping in order to substantially speed up the evolution process, which does not rely on the shortcuts to adiabaticity [89,90]. Secondly, the present scheme shows greater robustness and wider scalability for generating large-scale GHZ states than that in Ref. [15], owing to the fast topological pumping via edge channels. Finally, our method shows more simplicity and feasibility in the generalized SSH model confining to engineer nearest-neighbor couplings in comparison with introducing next-to-nearest-neighbor interactions in Ref. [13].

2 Physical model and engineering of topological pumping

2.1 Generalized SSH model and Zak phase

The schematic of our generalized SSH model is shown in Fig.1(a), which describes a one-dimensional coupled superconducting qutrit-resonator chain and composes of N qutrits and N resonators, respectively. Each flux qutirt holds a three-level structure, forming two ground states |L⟩,|R⟩ and one excited state |e⟩. The frequency of a single-mode resonator B_n is set to match resonantly with the transition frequency of |L⟩\leftrightarrow |e⟩ ( |R⟩\leftrightarrow |e⟩) of the two nearest-neighbor qutrits A_n and A_{n+1}, when n is odd (even). Between two adjacent cells, the resonator B_n is coupled to its front qutrit A_n and the behind qutrit A_{n+1} with alternating coupling strengths J_1 and J_2, respectively. The interaction of the chain can be described by the following interaction-picture Hamiltonian ( \hslash = 1)

(1)\begin{array}{r}{H}_I=\sum _n(J_1|e{⟩}_n⟨j_n|+J_2|e{⟩}_{n+1}⟨j_n|)b_n+\text{H.c.},\end{array}

where j_n=L(R) when n is odd (even), and b_n the annihilation operator of the resonator B_n. Further, the one-dimensional coupled superconducting qutrit-resonator chain is analogous to an SSH model, which can be simplified as

(2)\begin{array}{r}{H}_I=\sum _n{J}_1{b}_n{S}_n^{†}+J_2{b}_n{S}_{n+1}^{†}+\text{H.c.},\end{array}

where S_n^{†}=|e{⟩}_n⟨R|( {\sigma }_n^x{)}^n and {\sigma }_n^x =|L{⟩}_n⟨R|+|R{⟩}_n⟨L|. The operator ( {\sigma }_n^x{)}^n works for two different transitions |L⟩\leftrightarrow |e⟩ and |R⟩\leftrightarrow |e⟩ for qutrit A_n depending on the odevity of n.

For periodic boundary conditions, we can use the Bloch theorem and rewrite H_I as

(3)

with {\mathrm{\Psi }}_{k_n}^{†}=({\psi }_{S,k}^{†},{\psi }_{B,k}^{†})= N^{-1/2}\sum _{n=1}^N{\mathrm{e}}^{\mathrm{i}{a}_{0}nk}(S_n^{†},b_n) and

(4)\begin{array}{r}{\mathcal{H}}^B(k)=J_2\left(\begin{array}{cc}0& \rho (k)\\ {\rho }^{\ast }(k)& 0\end{array}\right),\end{array}

where \rho ( k)=J_1/J_2 +{\mathrm{e}}^{-\mathrm{i}k{a}_0} with k\in [-\pi / a_0,\pi /a_0]. Introducing the vector of Pauli matrices \sigma =({\sigma }_x, {\sigma }_y), the Hamiltonian can be expressed in the form {\mathcal{H}}^B(k)= - J_2\mathit{g}(k)\cdot \sigma with \mathit{g}(k) =(\text{Re}(\rho ),-\text{Im}(\rho )). The eigenvalues can be obtained by diagonalizing {\mathcal{H}}^B(k),

(5)\begin{array}{r}{ϵ}_{k,\pm }=\pm J_2|\mathit{g}(k)|=\pm J_2|\rho (k)|=\pm \sqrt{J_1^2+J_2^2+\cos (k{a}_0)}.\end{array}

The corresponding eigenstates are

(6)\begin{array}{r}{u}_{\pm }(k)=\frac{1}{\sqrt{2}}\left(\begin{array}{c}{\mathrm{e}}^{-\mathrm{i}\varphi (k)}\\ \pm 1\end{array}\right),\end{array}

where the phase \varphi ( k) is given by \cot \varphi (k)=\frac{J_1/J_2}{\sin k{a}_0}+\cot k{a}_0. Analyzing Eqs. (5) and (6), the eigenenergy spectrum of the system is divided into two bands, i.e., the negative eigenenergy band u_{-}( k) and positive eigenenergy band u_{+}(k). This model belongs to the BDI class according to the standard topological classification [91] and possesses two topological distinct phases for and J_1>J_2, where the band gap is closed at the boundaries of the first Brillouin zone ( k=\pm \pi /a_0). The chain is more like a “conductor” when the system lies in the phase transition point at J_1=J_2. Otherwise, the chain behaves as an “insulator”. However, there exists a band gap with width 2 | J_1- J_2| except for the closed points. It is noted that the trajectory of the vector \mathit{g}(k) for all k enclose or does not enclose the origin and |\varphi ( k)|\in [0, 2\pi ] or when or J_1>J_2. This topological behavior of the phase \varphi (k) is bound up with the value of the Zak phase, given by

(7)\begin{array}{r}\mathcal{Z}=\mathrm{i}{\int }_{-\pi /a}^{\pi /a}\mathrm{d}k[u_{-}^{†}(k){\mathrm{\partial }}_k{u}_{-}(k)]=\frac{1}{2}\propto \mathrm{d}k\frac{\mathrm{d}\varphi }{\mathrm{d}k}=\frac{\mathrm{\Delta }\varphi }{2},\end{array}

where \mathrm{\Delta }\varphi is the variation of \varphi ( k) when k varies across the full Brillouin zone. After some algebra, one can obtain \mathcal{Z}=\pi and \mathcal{Z}=0 for the cases of and J_1>J_2, respectively. The winding of the vector \mathit{g}(k) when k varies across the Brillouin zone is shown in Fig.2(a) and (b) for two values of and J_1/J_2>1, respectively. When , the loop encloses the origin and the phase \varphi ( k) can take an any value. As for J_1/J_2>1, the curve \mathit{g}(k) does not enclose the origin and for all k. Thus, topological phenomenon of the phase \varphi (k) is related to the value of Zak phase. The Zak phase \mathcal{Z} is \pi times the winding number of the curve \mathit{g}(k) around the origin and is, therefore, zero if this curve does not enclose the origin while \pi if it does. Accordingly,

(8)

The results show that tuning the ratio J_1/J_2 induces a topological phase transition characterized by the Zak phase.

2.2 Analaysis of fast topological pumping and edge state

As for even- and odd-sized SSH chains, the characteristics of energy spectrum versus the varying J_1/J_2 demonstrate the different distribution of zero energy mode [92]. Under the open boundary condition for the even-sized SSH chain, the two phases are distinguished by the presence and absence of degenerate zero-mode edge states with topologically nontrivial phase and topologically trivial phase J_1/J_2>1, as shown in Fig.2(c). As expected in Fig.2(d), there is only one topologically protected zero-mode edge state in the band gap with regardless to the value of J_1/J_2. Therefore, the odd-sized chain is always in the topological nontrivial phase [93].

In this protocol, the size of superconducting qutrit-resonator chain is restricted to the odd number, the interaction Hamiltonian of which can be written as H_{2N-1}=\sum _{n=1}^{N-1}{J}_1{b}_n S_n^{†}+J_2 b_n{S}_{n+1}^{†}+\text{H.c.}. Accordingly, the topological protected zero-mode edge state of the odd-sized chain with a single excitation are exponentially localized at the boundaries, which can be obtained analytically (see Appendix A for details)

(9)\begin{array}{r}|{\phi }_E⟩=\sum _{n=1}^N{\lambda }^n{\sigma }_n^{+}({\sigma }_n^x{)}^{n-1}\underset{l=1}{\overset{n-1}{⨂}}{\sigma }_l^x|G{⟩}_L,\end{array}

where {\sigma }_n^{+}=|e{⟩}_n⟨R| and |G{⟩}_L=|RLR\cdots {⟩}_{A_N}\otimes |000\cdots {⟩}_{B_{N-1}} denotes a decoupled state of the qutrit-resonator chain with A_1 in |R⟩, B_1 in |0⟩ (zero-photon Fock state), A_2 in |L⟩, B_2 in |0⟩, \cdots and \lambda =- J_1/J_2. The system can perform a topologically protected state transfer on the zero mode from the left edge to the right edge, corresponding to an excitation transfer from one side to the other, when setting initially J_1/J_2=0 while finally J_1/J_2=+\mathrm{\infty }.

It is evident that the above topological protected transfer process needs to satisfy the adiabatic limit [8,12,13]. A sufficient condition for the adiabatic evolution is \mathrm{\Delta }E=|E_m( t)- E_n(t)|\gg |⟨{\dot{\phi }}_m(t) | {\phi }_n(t)⟩|, where E_m, E_n, |{\phi }_m( t)⟩ and |{\phi }_n(t)⟩ are the mth and nth instantaneous eigenenergies and corresponding eigenstates of the Hamiltonian, and the overdot indicates differentiation with respect to time. In general, the eigenenergy difference |E_m(t)-E_n(t )| between |{\phi }_m( t)⟩ and |{\phi }_n(t)⟩ is very small. Therefore, for the sake of satisfying the adiabatic limit, the time evolution of the adiabatic state transfer should be very slow to make the evolution remain in the zero mode (edge state) without exciting other eigenstates (bulk states) due to a trivial \mathrm{\Delta }E. The adiabatic limit condition can be relaxed to realize fast topological pumping by increasing the eigenenergy difference \mathrm{\Delta } E to form a nontrivial band gap in the Hamiltonian spectrum or/and decreasing the derivative of the Hamiltonian.

In order to realize fast topological pumping of edge states by relaxing the adiabatic limit condition, it is crucial to suitably adjust the coupling strengths. Here we specify that the couplings are shaped by an exponential function

(10)\begin{array}{rl}{J}_1& =J_0\left(\frac{1-{\mathrm{e}}^{-\alpha t/T}}{1-{\mathrm{e}}^{-\alpha }}\right),\\ J_2& =J_0\left[\frac{1-{\mathrm{e}}^{-\alpha (T-t)/T}}{1-{\mathrm{e}}^{-\alpha }}\right],\end{array}

where T is the total evolution time and \alpha a free parameter that can be fine-tuned [12]. The forms of J_{1,2} in Eq. (10) satisfy the state transfer conditions J_1/J_2{|}_{t\to 0}=0 and J_1/J_2{|}_{t\to T}=+\mathrm{\infty } well. To illustrate the crucial characteristics of topological pumping for coupling strengths and analyze how the energy spectrum evolves over time, the functions of coupling strengths and the corresponding instantaneous energy spectrum versus time are plotted in the left plane and the right plane of Fig.3, respectively. We take the size L=2N-1= 61 of the chain as an example and select three different values of free parameter \alpha = 1 in (a) and (b), \alpha = 6 in (c) and (d), and \alpha = 11 in (e) and (f), respectively. In the odd-sized chain, there is always one zero-mode state remaining unchanged and is separated from bulk states, shown in the right plane of Fig.3. At the beginning, the energy gap separating the edge state from the bulk states takes its maximum value. Subsequently, the energy gap approaches its minimum, occurring at J_1=J_2, and at the end of evolution time the energy gap regains its maximum value.

Comparing the functions with different values of free parameter \alpha employed in the topological chain, we can notice their qualitative differences. In particular, there is the most significant eigenenergy difference \mathrm{\Delta }E between the bulk state and the edge state for the smallest value of \alpha. The functions with \alpha =1 approach the minimum value of the energy gap with the largest slope, while the functions with \alpha = 11 drive the system more gently in this region, shown in the left planes of Fig.3. Depending on the value of free parameter \alpha, the minimum value of energy gap can be changed, which is the key point to avoid nonadiabatic effects and to speed the state transfer. For the odd-sized chain with length L, the eigenmode energies (besides the zero-energy solution) can be analytically given by

(11)\begin{array}{r}g=|J_1+J_2{\mathrm{e}}^{\mathrm{i}{q}_m}|,\end{array}

where q_m=2m\pi /(L+1) with m=1,2,\cdots ,[L/2], and we define that [z] gives the greatest integer that is less than equal to z [94]. Thus, the energy gap can be analytically determined by \chi =2|g_{[L/2]} |, which is related to the values of J_1 and J_2. Note here that, there always exists a time such that J_1=J_2. In the fast topological pumping with exponential couplings, we consider this time to be t=T/2 in the right plane of Fig.3. In the infinite system, J_1=J_2 corresponds to the closing of the energy gap separating the zero-mode state with the rest of energy states. Thus, the point in the parameter space where the coupling strengths J_1 and J_2 are equal corresponds to the minimum value of energy gap {\chi }_{\text{min}} during the whole evolution process. The minimum values of energy gap are analytically obtained by {\chi }_{\text{min}}^{\alpha =1}=0.0631 J_0, {\chi }_{\text{min}}^{\alpha =6}=0.0965 J_0 and {\chi }_{\text{min}}^{\alpha =11}=0.1009 J_0, which are labeled with tags in the right plane of Fig.3.

Based on time-dependent exponential couplings with different values of the free parameter \alpha, we plot the state distribution of the zero mode versus evolution time in Fig.4(a)−(c). The population of zero-mode state is defined by P(t)= ⟨{\phi }_E|\rho (t) | {\phi }_E⟩, where |{\phi }_E⟩ and \rho ( t) are the target zero-mode state and the time-dependent density matrix of the system by solving Liouville equation \dot{\rho }( t)=- \mathrm{i}[H_{2N-1},\rho ( t)], and t corresponds to the evolution instant. Concretely, the dark blue fringes show that almost no population of zero-mode state is distributed in the lattice sites. However, the bright yellow fringes indicate that population about 0.2 of zero-mode state is emerged in the even lattice sites at the specific evolution moments. The most bright yellow fringes mean that near-unity population of zero-mode state appears in the first and the last lattice sites at the begin and end of evolution. In Fig.4(a)−(c), the zero-mode state is localized at the first site A_1 or the last A_{61} and the distribution of left or right edge state with zero mode is equal to unity when J_{0}t\to 0 or J_{0}t\to +\mathrm{\infty } to realize the topological protected edge state transfer in the odd-sized SSH chain. To be specific, the value of parameter \alpha is related to the total evolution time. It is noted that the shortest evolution time with the size of chain L=2N-1=61 approaches to T=260/J_0 under the condition of \alpha =6. Besides, the distribution of zero-mode state at odd sites are enlarged by increasing the value of parameter \alpha in the time scale, seen from bright fringes in Fig.4(a)−(c). However, the population of zero mode increases exponentially on the odd sites (qutrits) while the zero-mode state exhibits no populations on the even sites (resonators), indicating the virtual photon excitation in the resonators. The reason lies in the interaction property of the system, that is, the formulation of the Hamiltonian, which determines that the distribution of zero-mode state of system on different lattice sites. In the SSH model, according to H_{2N-1}|{\phi }_E⟩=0, one can calculate the analytic form of the zero-mode state of the system (see Appendix A), for which it is exactly one of the most important features of the SSH model that the zero-mode state only distributes on odd-size lattice sites while involves no even-size ones.

As shown in Fig.4(d), by selecting different \alpha and the corresponding minimum energy gap as numerical samples, it can be found that the maximization of minimum energy gap needs a large value of parameter \alpha. A larger value of parameter \alpha ( \alpha \in [1,11 ]) determines J_1 and J_2 to equate at a higher value. As for \alpha \ge 11, the minimum energy gap can be maximized and keeps around 0.1J_0 because equating J_1 and J_2 at values close to the maximum value of couplings can be acquired during the transfer. On the one hand, the minimum energy gap can be actually used to specify a characteristic timescale. It can be assumed that the system is close to the adiabatic following of the zero-mode state, under the condition of the sufficiently longer transfer time by comparing with this timescale. On the other hand, the coupling functions are required with the steeper slope at first and the smaller slope at the middle time of evolution. In view of the above two key characteristics of the coupling functions, the exponential function with the tunable parameter \alpha seems to be one of the most appropriate to increase the speed of the transfer protocol.

3 Fast generation of large-scale GHZ states with great robustness and scalability

3.1 Fast generation of large-scale GHZ states

We focus on the generation of large-scale GHZ states among the N qutrits based on the fast topological pumping via edge states. The implementation of zero-mode edge state transfer by engineering exponential coupling functions via the topological protected edge channel has been shown in Fig.4(a)−(c). Thus, the left and right edge states localized in the odd-sized chain at the beginning and end of evolution time can be expressed, respectively, as

(12)\begin{array}{rl}& |l{⟩}_L=|eLR\cdots m{⟩}_{A_N}\otimes |000\cdots 0{⟩}_{B_{N-1}},\\ & |r{⟩}_L=|LRL\cdots e{⟩}_{A_N}\otimes |000\cdots 0{⟩}_{B_{N-1}},\end{array}

where m=R(L ) when N is odd (even). Initially, we suppose that the state of the chain is |{\mathrm{\Psi }}_0⟩= \left(|G{⟩}_L+|l{⟩}_L\right)/\sqrt{2} meaning that the first qutrit are prepared in a superpostion (|R⟩+ | e⟩)/\sqrt{2}. The evolution via the protected edge channel takes the form of

(13)\begin{array}{rl}& |G{⟩}_L⟶|G{⟩}_L,\\ & |l{⟩}_L⟶-(-1{)}^N|r{⟩}_L,\end{array}

which is irrelevant with the state |G{⟩}_L because of zero-photon in the resonators that cannot excite the ground-state qutrits. We specify that logical states 0 and 1 are carried by the ground levels |L⟩ and |R⟩, respectively, for the qutrit A_n ( 1\le n\le N- 1). As for the last qutrit A_N, we set that |e⟩ is employed by the logical state 1 (0) when N is even (odd). Consequently, a large-scale GHZ state of N qutrits without regard to the zero-photon product state of resonators is obtained and reads as

(14)\begin{array}{r}\frac{1}{\sqrt{2}}(|101\cdots 1(0){⟩}_{A_N}-(-1{)}^N|010\cdots 0(1){⟩}_{A_N}).\end{array}

This is exactly an N-body GHZ state according to a general form of x-qubit GHZ state which can be expressed as (|m_1{m}_2{m}_3\cdots m_x⟩+{\mathrm{e}}^{\mathrm{i}\varphi }|n_1{n}_2{n}_3\cdots n_x⟩)/\sqrt{2} where m_j+n_j= 1 ( m_j,n_j\in 0,1) [95-98].

In the following, we take \alpha =6 as an example and consider L=61 to plot the evolution of population based on the fast topological pumping with exponential couplings and the conventional adiabatic topological pumping with Gauss couplings [15] for the ideal 31-body GHZ state |{\mathrm{\Psi }}_{\text{ideal}}⟩= (|G{⟩}_{61}+|r{⟩}_{61})/\sqrt{2}, the initial state |{\mathrm{\Psi }}_{\text{initial}}⟩= (|G{⟩}_{61}+|l{⟩}_{61})/\sqrt{2}, the left edge state |l{⟩}_{61}=|eLR\cdots R {⟩}_{A_{31}}\otimes |000\cdots 0{⟩}_{B_{30}}, the right edge state |r{⟩}_{61}=|LRL\cdots e{⟩}_{A_{31}}\otimes |000\cdots 0{⟩}_{B_{30}} and the decoupled state |G{⟩}_{61}=|RLR\cdots {⟩}_{A_{31}}\otimes |000\cdots 0{⟩}_{B_{30}}. As expected in Fig.5(a), a complete state transfer from the left edge state |l{⟩}_{61} to the right edge state |r{⟩}_{61} is manifested and the decoupled state |G{⟩}_{61} remains unchanged, resulting in the successful creation of a 31-body GHZ state in both fast topological pumping with exponential couplings and the conventional adiabatic topological pumping with Gauss couplings. However, the populations of states |r{⟩}_{61}, |l{⟩}_{61}, |{\mathrm{\Psi }}_{\text{initial}}⟩ and |{\mathrm{\Psi }}_{\text{ideal}} ⟩ based on the fast topological pumping with exponential couplings tends to 0.5, 0, 0.25 and 0.999, respectively, only at the end of evolution time. In respect to the conventional adiabatic topological pumping, the population of those states can also attain relatively ideal results at t\approx 5000/J_0 and keep unchanged until at the end of evolution time.

In order to verify the generation of the large-scale GHZ state via fast topological pumping in a short evolution time, we make a comparison between this protocol and the conventional adiabatic topological pumping. The quantity that determines how faithfully the evolution has occurred is the fidelity, which can be formulated as F(t) =⟨{\mathrm{\Psi }}_{\text{ideal}}|\rho (t) | {\mathrm{\Psi }}_{\text{ideal}}⟩, where |{\mathrm{\Psi }}_{\text{ideal}}⟩ and \rho (t) are the target state and the time-dependent density matrix of the system by solving Liouville equation \dot{\rho }(t)= - \mathrm{i}[H_{2N-1},\rho ], and t corresponds to the evolution instant. In Fig.5(b), for each protocol, we plot the final fidelity F(T) as a function of the total evolution time T. In the case of reaching the fidelity of 0.999, it can be found that the protocol with exponential couplings based on the fast topological pumping is greatly faster than that with Gaussian couplings, because of the former occurring at T_{\text{Exp}}\approx 390/J_0 as compared to the latter where this happens at T_{\text{Gauss}}\approx 4293/J_0. Fig.5 indicates the successful and fast creation of large-scale GHZ state via the fast topological pumping in this protocol.

3.2 Effect of different values of α

The implementation of fast topological pumping via edge channels for generating large-scale GHZ states in an odd-sized topological chain has been shown above, which takes the free parameter \alpha =6 as an example. Now that the qualitative differences among coupling functions, energy spectrum and the distribution of the zero-energy mode for different values of free parameter \alpha have become apparent in 2.2 here we further examine some quantitative results for generating large-scale GHZ states based on the fast topological pumping.

In the following, Fig.6(a)−(c) show that the final fidelities of 5, 18 and 31-body GHZ state as functions of the total evolution time by setting \alpha =1,6,11,16 and 21, respectively. Apparently, an N-body GHZ state needs a longer total evolution time with the increase of N. In the limit of T\to +\mathrm{\infty }, the fidelity of N-body GHZ state approaches unity for all the \alpha parameters, which indicates a perfect transfer of excitation along the chain. In Fig.6(a), there are slight differences among different values of \alpha for the total evolution time of 5-body GHZ state. However, it can be noticed that the fidelity curve tends to oscillate for T{J}_0\in [ 20,100], owing to the action of resonant processes. As for the 18- and 31-body GHZ states in Fig.6(b) and (c), there are more significant oscillations for \alpha =11 ,16 and 21, respectively.

According to the adiabatic condition, the oscillations during the evolution are strongly associated with the energy gap between the zero-mode energy E_N and the most-closed-to-zero eigenenergy E_{N+1}. In order to study the oscillations appearing in the fidelity of 5-, 18- and 31-body GHZ states, the corresponding instantaneous energy E_{N+1} in cases of \alpha =1,6 ,11,16 and 21 are plotted in the bottom plane of Fig.6. As for the small-scale GHZ state of N=5 in Fig.6(d), the parameter values \alpha = 6,11,16 and 21 have almost same effects on the corresponding instantaneous energy E_{N+1}. When t\in [0.1,0.9]T, the instantaneous energy E_{N+1}\approx 0.6J_0 is always closed to the zero-mode energy so that the slight oscillations appear in the adiabatic evolution process. However, the instantaneous energy E_{N+1}\approx 0.4J_0 in the case of the parameter \alpha =1 approaches to zero-mode energy only at the middle of evolution, closer to zero so that the resonance process is inhibited inconspicuously when compared with cases of other \alpha. In Fig.6(e) for generating 18-body GHZ state, the instantaneous energy E_{N+1} indicates similar evolution in cases of \alpha = 11,16 and 21 to form the most obvious oscillations in Fig.6(b). In respect to the cases of the parameter \alpha =1 and 6, there is few oscillations during the adiabatic evolution owing to the greater and greater separation between the zero-mode energy and the energy E_{N+1} except for the middle of evolution in Fig.6(e). By increasing N in Fig.6(f), the corresponding instantaneous energy E_{N+1} of 31-body GHZ state with different values of \alpha verges on the zero-mode energy more closely in contrast to 5- and 18-body GHZ state. When concentrating on the total evolution time T for different values of \alpha, the decreasing value of energy gap between the zero-mode energy and the energy E_{N+1} indicates longer evolution time to attain high fidelity N-body GHZ states. Nevertheless, the resonant processes can be properly handled to generate the large-scale GHZ state with a high efficiency by tuning the parameter \alpha.

Actually, the above phenomenon can be further understood by characteristics of coupling functions and the corresponding energy spectrum in Fig.3. Smaller values of the parameter \alpha lead to a less steep slope of the coupling function, resulting J_1 and J_2 equating at a smaller value. Simultaneously, the topological edge state at zero energy is separated from bulk states better as the function of evolution time except for t=T/2. Owing to a better separation between the edge state and the bulk states, the resonant processes can be suppressed effectively. Specially, the minimum value of energy gap {\chi }_{\text{min}} is demanded as large as possible, which makes the chain evolve adiabatically along the zero-mode edge state without exciting other eigenstates at t=T/2. Therefore, the parameter \alpha is closely related to two aspects, a smaller value for suppressing the nonadiabatic transition between the edge state and bulk states and a larger value for obtaining a greater {\chi }_{\text{min}}. Typically, a suitable value of parameter \alpha is the key point to generate high fidelity of N-body GHZ state based on the fast topological pumping, which can not only avoid the resonant processes but also shorten the total evolution time.

In the following, we investigate the fidelity of 5-, 18-, 31- and 44-body GHZ states versus the varying \alpha and the total evolution time T in Fig.7(a)−(d), respectively. InFig.7(a), the appropriate range of the parameter \alpha is chosen as \alpha \in [1,6] to generate 5-body GHZ state with the high fidelity and the short total evolution time. The 0.995 and 0.999 fidelity contour lines exhibit strongly oscillations for all values of the parameter \alpha and the total evolution time T, which is similar to Fig.6(a). As for 18- and 31-body GHZ states in Fig.7(c) and (d), respectively, N-body GHZ state with the high fidelity and the short total evolution time demands the parameter \alpha \in [6, 11] approximately. Considering the 0.995 and 0.999 fidelity contour lines of 18-, 31- and 44-body GHZ states, the oscillations are suppressed substantially under the range of \alpha \in [1,6 ]. Based on the fast topological pumping for generation of N-body GHZ state, it is no wonder that the value of parameter \alpha =6 is a trade-off, since the chain is driven strongly to increase the speed but also gently enough to avoid resonant processes.

In order to more intuitively reflect the speed of protocols for the fast topological pumping and the conventional adabatic topological pumping, we select different N and corresponding evolution times of generating GHZ states with the fidelity of 99.9\mathrm{\%} as numerical samples to fit the function of the total evolution time T versus N in Fig.8(a). The fitting functions for Gauss and exponential couplings are J_{0}T=6.9419N^2+2.455N-59.8933 and J_{0}T=0.3691N^2-2.06 N+25.9, respectively. As N increases, the total evolution time of GHZ state shows a quadratic trend with N for both the Gauss and exponential couplings. However, an 80-body GHZ state with 99.9 \mathrm{\%} fidelity requires the total evolution time T_{\text{Exp}}\approx 35.4\text{µ}s but T_{\text{Gauss}}\approx 711\text{µ}s by choosing the value of J_0=2\pi \times 10 MHz [99]. In order to get the information from the shortest time of generating N-body GHZ states in the small scale based on the fast topological pumping with exponential couplings and the conventional adiabatic pumping with Gauss couplings, we numerically calculate the fidelity of ideal N-body GHZ states in the small scale, and plot the time evolution of {\log }_{10}(1-F) with N ranging from 10 to 30 at intervals of 5 in Fig.8(b). The numerical results based on the conventional adiabatic pumping exhibit that a higher fidelity needs a longer evolution time with increasing N. As for the fast topological pumping, the evolution time for attaining 0.999 fidelity of N-body GHZ state has a trivial difference among different values of N in contrast to the conventional adiabatic pumping when by enlarging the range of T\in [0,300]/J_0, owing to the instantaneous energy spectrum characteristics of small-scale systems based on the fast topological pumping. It is evident that the efficiency of evolution based on the fast topological pumping improves more significantly with a larger N than the conventional topological pumping.

3.3 Robustness against disorders and losses in the superconducting qutrit-resonator chain

The one-dimensional superconducting qutrit-resonator chain may possess not only itself defects and perturbation but also dissipative processes, which could result in infidelities for the perfect generation of N-body GHZ states. Three dominant aspects of influence are considered here: (i) unwanted on-site potential defects for qutrits A_n and resonators B_n; (ii) inevitable variation in ideal couplings; (iii) losses of qutrits A_n and resonators B_n with decay rates {\kappa }_n and {\gamma }_n, respectively.

In order to study the robustness of the protocol against the on-site potential defects and the variation of ideal couplings, the chain can be described by

(15)\begin{array}{rl}{H}_d=& \sum _{n=1}^N{\delta }_{A,n}|e{⟩}_n⟨e|+\sum _{n=1}^{N-1}{\delta }_{B,n}{b}_n^{†}{b}_n+\sum _{n=1}^{N-1}(J_{1,n}^{{}^{\prime }}|e{⟩}_n⟨j_n|\\ & +J_{2,n}^{{}^{\prime }}|e{⟩}_{n+1}⟨j_n|)b_n+\text{H.c.},\end{array}

where J_{1(2),n}^{{}^{\prime }}=J_{1(2)}+\text{rand}[-{\delta }_{J_{1,2}},{\delta }_{J_{1,2}}], {\delta }_{A(B),n}=\text{rand}[-{\delta }_{A(B),n},{\delta }_{A(B),n}] and rand [- x,x] denotes a random number in the range of [-x, x]. For convenience, the disorder strengths of on-site potential defects and couplings are chosen as {\delta }_{A(B),n}={\delta }_{J_{1,2},n}=\delta as an example. The relation between the fidelity of 31-body GHZ state and the disorder \delta \in [0, 0.5] is exhibited in Fig.9. It is the disorder that is randomly sampled among 10 000 times, and then the fidelity is taken as an average of the 10 000 results versus the disorder of potential defects for all A-type lattice sites, B-type lattice sites and coupling strengths, respectively. The fidelity shows relatively great robustness for the disorder of coupling strengths and on-site potential defects. By magnifying the disorder for coupling strengths, the fidelity of 31-body GHZ state keeps still above 99.6 \mathrm{\%}. The yellow-squared and purple-crossed lines indicate that the on-site potential defects for all B-type lattice sites and all A-type lattice sites cause the same degree damage to the fidelity of 31-body GHZ state, which can retain 91\mathrm{\%} even if \delta =0.5. The robustness of N-body GHZ state to the disturbance and perturbation benefits not only from the topological protected edge channel but also the improvement of evolutionary efficiency, which will show a more obvious superiority even when increasing N.

Furthermore, we also investigate the relation among the fidelity of 31-body GHZ state, the total evolution time T and the disorder \delta for unexpected couplings, on-site potential defects of A- and B-type lattice sites based on the fast topological pumping in Fig.10(a), (c) and (d), respectively. The damage to fidelity caused by unexpected couplings and on-site potential defects for all A- and B-type lattice sites with \delta /J_0\in [ 0.2,0.5] can be compensated by longer evolution time T. Specially, we show 0.995 and 0.999 fidelity contour lines of 31-body GHZ state in Fig.10(a), (c) and (d), respectively. The robustness against the disorder from unexpected couplings shows a greater result than that from on-site potential defects due to the topological property, which is equivalent to Fig.9. Also, the relation among the fidelity of 31-body GHZ state, the total evolution time T and the disorder \delta for unexpected couplings based on the conventional adiabatic topological pumping is exhibited in Fig.10(b). The conventional adiabatic topological pumping shows the more robustness against the unexpected couplings than the fast topological pumping, learning from the red contour lines in Fig.10(a) and (b). This result comes from the fact that the conventional adiabatic topological pumping needs a longer adiabatic evolution time to offset the damage on the fidelity from the unexpected couplings than the fast topological pumping, indicating that the improvement in speed based on the fast topological pumping is in cost of robustness against disorder from the unexpected couplings.

The robustness against disorders is benefit from the topological property of the system, which can trace back to the chiral symmetry in the system. Concretely, the system is satisfied with \mathrm{\Gamma }{H}_{2N-1}{\mathrm{\Gamma }}^{-1}=- H_{2N-1}, where \mathrm{\Gamma }= \prod _n(S_n^{†}{S}_n-b_n^{†} b_n)= \prod _n{\sigma }_z [93,100]. Meanwhile, this symmetry indicates that the zero-energy mode can exist and results in a symmetric energy spectrum with each positive eigenenergy E accompanied by a negative eigenenergy -E. From an algebraic point of view, the Bloch Hamiltonian of the chain {\mathcal{H}}^B in Eq. (4) can be written as a linear combination of the Pauli matrices {\sigma }_x and {\sigma }_y while is irrelevant to {\sigma }_z. Thus, in the presence of unexpected couplings J_{1,2}^{{}^{\prime }}, the system still obeys the chiral symmetry. Besides, the robustness against the disorder from the unexpected couplings also profits from adiabatic passage of generating large-scale GHZ state based on fast topological pumping. The reason lies in that the exponential couplings are related to the adiabatic evolution time while on-site potential defects for all A- and B-type lattice sites are constant terms. However, the presence of on-site potential defects breaks the chiral symmetry of this system and modifies the energies of the edge state. For a large disorder, these states cannot be distinguished from the bulk. Consequently, the evolution along the zero-mode edge state may be affected gently from the on-site potential defects for qutrits and resonators while hardly from the unexpected couplings.

We now take the effects of losses of qutrits and resonators on the fidelity of GHZ states based on the fast topological pumping. The dynamics of the lossy system affected by the two dominant channels for the losses of qutrits and resonators, respectively, can be governed by the non-Hermitian Liouville equation \dot{\rho }=-\mathrm{i}(H^{{}^{\prime }}\rho -\rho H^{{}^{\prime }†}), where H^{{}^{\prime }}=H_{2N-1}-\mathrm{i}{\kappa }_n \sum _{n=1}^{N-1}{b}_n^{†} b_n/2-\mathrm{i}{\gamma }_n\sum _{n=1}^N|e{⟩}_n⟨e|/2 [58,101]. For convenience, we assume that {\kappa }_n= \kappa, {\gamma }_{q_n}={\gamma }_q and J_0/(2\pi )=10 MHz. Fig.11 shows the final fidelity versus the resonator decay rate \kappa and the qutrit decay rate {\gamma }_q based on the conventional adiabatic topological pumping with Gauss couplings and the fast topological pumping with exponential couplings. The upper and lower surfaces indicate that the qutrit decay brings down the fidelity significantly as {\gamma }_q increases, while the loss of resonators has slight effect on the fidelity, owing to the virtual excitation of photon in the resonator B_n. As for the fast protocol with exponential couplings, the fidelity is improved apparently in contrast to the lower surface for the conventional adiabatic topological pumping with Gauss couplings. When \kappa /(2\pi ) ={\gamma }_q /(2\pi )=3 kHz [102,103], the fidelity of 31-body GHZ state can reach around 25% and 65.25% for Gauss and exponential couplings, respectively. The above results manifest the protocol based on the fast topological pumping enhances robustness largely against losses of the resonators and qutrits.

3.4 Scalability of GHZ states

As stated above, our protocol shows a remarkable improvement against both disorders and losses in the chain. This robustness is not only related to topological protection, as well as the enhancement of evolutionary efficiency, which stems from the short timescales to enable the fast topological pumping according to the crucial characteristics of energy spectrum and zero-mode edge state transfer.

In order to verify more quantitatively the effect of our protocol, one more vital direction is the scalability of GHZ state. In the following, we take into account the losses of qutrits and resonators with different coupling strengths J_0/(2\pi ) =1,10 and 50 MHz. For the fixed and experimentally available decay rates of qutrits and resonators with {\gamma }_q/(2\pi )=\kappa /(2\pi )=1 kHz [104-108], the relation between final fidelities of GHZ states and the scalability of entanglement N are exhibited in Fig.12. The unmarked and marked lines represent the conventional adaibatic topological pumping with Gauss couplings and the fast topological pumping with exponential couplings, respectively. Obviously, the fidelity exhibits decreasing tendency with increasing N for both two protocols. However, a greater value of coupling strength J_0 can be applied to effectively suppress the trend of decline. Furthermore, the advantage of our protocol becomes more pronounced as the scalability of GHZ state is strengthened for the fast topological pumping with exponential couplings. When J_0/(2\pi ) =1 MHz, the fidelity of 30-body GHZ state can stay above 90\mathrm{\%} for exponential couplings while just 35\mathrm{\%} for Gauss couplings. As for the GHZ state with N=150 and the coupling strength J_0/(2\pi ) =50 MHz, the fidelity of the fast topological pumping and the conventional topological pumping can reach 88\mathrm{\%} and 50 \mathrm{\%}, respectively. Thus, by inspecting Fig.12, it is evident that the fast topological pumping outperforms a lot the conventional adaibatic topological pumping in terms of the scalability of entanglement N.

Furthermore, we consider the scalability of N-body GHZ state in a two-dimensional square lattice with the size of L\times L based on the superconducting qutrit-resonator system in Fig.13. The fast topological pumping occurs at the ith line and the (2N-1)th column of SSH chain to generate N-body GHZ states, including three steps:

Step 1 The (2N-1)\text{th} column of SSH chain is initialized in the state |{\mathrm{\Psi }}_{N,L}⟩ =\left(|G{⟩}_{N,L}+|M{⟩}_{N,L}\right)/\sqrt{2}, where |G{⟩}_{N,L}=|RLR\cdots {⟩}_{A_{N,N}}\otimes |000\cdots {⟩}_{B_{N-1,N}} and |M{⟩}_{N,L}= |LRL\cdots LLRLR\cdots m{⟩}_{A_{N,N}}\otimes | 000\cdots 0{⟩}_{B_{N-1,N}}. In particular, it is supposed that the states of the qutrit A_{i-1,N} and the qutrit A_{i,N} are |L{⟩}_{A_{i-1,N}} and |L{⟩}_{A_{i,N}}, respectively. It is noted that such a special initial state can be represented by |{\mathrm{\Psi }}_{N,L}⟩=| {\mathrm{\Psi }}_{N,L}^1⟩\otimes | {\mathrm{\Psi }}_{N,L}^2⟩, where |{\mathrm{\Psi }}_{N,L}^1⟩=(|RL R\cdots RLR{⟩}_{A_{i-1,N}}+|LRL\cdots LRL{⟩}_{A_{i-1,N}})\otimes |000\cdots {⟩}_{B_{i-2,N-1}}/\sqrt{2} is exactly an (i-1 )-body GHZ state that can be attained by fast topological pumping in 1D SSH chain and | {\mathrm{\Psi }}_{N,L}^2⟩= | LRL\cdots m {⟩}_{A_{i,N}}\otimes |000\cdots {⟩}_{B_{i,N-1}}.

Step 2 The topological protected zero-mode state transfer occurs from the qutrit A_{i,1} to the qutrit A_{i,N} along the ith line of SSH chain. Initially, we suppose that the state of the ith line of SSH chain is |l{⟩}_{i,L}, whose evolution process takes form of

(16)| l{⟩}_{i,L}⟶- (-1 {)}^i|r{⟩}_{i,L},

where |l{⟩}_{i,L}=|eLR\cdots m {⟩}_{A_{i,N}}\otimes |000\cdots 0{⟩}_{B_{i,N-1}} and |r{⟩}_{i,L}= |LRL\cdots e {⟩}_{A_{i,N}}\otimes |000\cdots 0{⟩}_{B_{i,N-1}}.

Step 3 Owing to the fast topological pumping from Step 1, the state of the (2N-1)th column of SSH chain is |{\mathrm{\Psi }}_{N,L}^{{}^{\prime }}⟩=(|G{⟩}_{N,L}+|M^{{}^{\prime }}{⟩}_{N,L})/\sqrt{2}, where |M^{{}^{\prime }} {⟩}_{N,L}= |LRL\cdots LeRLR\cdots m{⟩}_{A_{1,N}}\otimes |000\cdots 0{⟩}_{B_{N-1,N}}. Subsequently, the topological protected channel is established from the qutrit A_{i,N} to the qutrit A_{N,N} in the (2N-1)th column of SSH chain, whose evolution process is followed by

(17)\begin{array}{rl}& |G{⟩}_{N,L}⟶|G{⟩}_{N,L},\\ & |M^{{}^{\prime }}{⟩}_{N,L}⟶-(-1{)}^N|r{⟩}_{N,L},\end{array}

where |r{⟩}_{N,L}=|LRL\cdots e {⟩}_{A_{i,N}}\otimes |000\cdots 0{⟩}_{B_{i,N-1}}.

Consequently, we specify that logical states 0 and 1 are carried by the ground levels |L⟩ and |R⟩, respectively, for the qutrit A_{i,N} ( 1\le i\le N- 1) and the zero-photon product state of resonators are discarded, which can also generate N-body GHZ state (14) in two-dimensional square lattice of superconducting qutrit-resonator system.

4 Experiment consideration for superconducting circuit devices

This protocol for generating large-scale GHZ states is applicable to superconducting circuit devices, which is benefit from existing circuit-QED technologies. We can construct a superconducting qutrit-resonator chain to arrange alternately the LC resonators and the flux qutrits in one-dimensional space, whose equivalent circuit of one unit cell is shown in Fig.14. The flux qutrit A_n consists of a superconducting loop interrupted by two large Josephson junctions E_1=E_3 =E_J and one small Josephson junction E_2={\alpha }_2{E}_J. This kind of qutrit based on three junctions enables a reduction of the loop size while retaining a large inductance [109,110]. To create a double-well potential, is used to avoid effects of charge noise. The additional Josephson junction and the additional coupler capacitor parallel with the flux qutrit A_n and the resonator B_n, which acts as a tunable coupler with one junction E_4 controlled by the external flux-bias line (FBL) to modify the coupler properties [111-113]. The resonator B_n is composed by a spiral inductor L_B and a capacitor C_B in analogy with LC harmonic oscillator, which has a single mode. In terms of the capacitor charge Q and the inductor current I, the Hamiltonian of LC oscillator is written as

(18)H_{\text{LC}}={\displaystyle \frac{Q^2}{2C_B}}+{\displaystyle \frac{{\mathrm{\Phi }}_B^2}{2L_B}} ,

where {\mathrm{\Phi }}_B is the flux through the inductor and Q the charge on the capacitor. Based on the standard quantization process of an LC circuit [105], the Hamiltonian of the resonator B_n can be further written as H_{\text{LC}}=\hslash {\omega }_{b_n}{b}_n^{†} b_n in terms of the creation and annihilation operators defined by b_n^{†}=1/\sqrt{2\hslash {\omega }_{b_n}}(Q /\sqrt{C_B}- \mathrm{i}{\mathrm{\Phi }}_B/\sqrt{L_B}) and b_n=1/\sqrt{2\hslash {\omega }_{b_n}}( Q/ \sqrt{C_B}+\mathrm{i}{\mathrm{\Phi }}_B/\sqrt{L_B}), where {\omega }_{b_n}= 1/\sqrt{L_B{C}_B} is the oscillator frequency and can be engineered in a large range of possible values by adjusting the parameters L_B and C_B.

At the same time, the electric dipole interaction between the flux qutrit A_n and the resonator B_n can be realized by using a tunable coupler without considering the charge interaction [114]. A more versatile design, shown in Fig.14, separates the three-junction flux qutrit and the resonator by an additional Josephson junction, while keeping the flux qutrits quantization independent of an additional flux difference along the segment \mathrm{\Delta }\mathrm{\Psi }, which is the source of the coupling. The corresponding Hamiltonian for the three-junction flux qutrit and the additional Josephson junction is [111,115,116]

(19)H_{\text{q}}={\displaystyle \frac{P_p^2}{2M_p}}+{\displaystyle \frac{P_m^2}{2M_m}}+U_{\text{eff}},

where M_p=2C_J [{\mathrm{\Phi }}_0/(2 \pi ){]}^2 and M_m=M_p(1 +2\alpha ) can be considered effective masses, while P_p=- \mathrm{i}\hslash \mathrm{\partial }/\mathrm{\partial }{\varphi }_p and P_m=-\mathrm{i}\hslash \mathrm{\partial }/\mathrm{\partial }{\varphi }_m can be considered effective momenta. Moreover, {\varphi }_p=({\varphi }_1+{\varphi }_3+{\varphi }_2)/2 and {\varphi }_m =({\varphi }_1+{\varphi }_3-{\varphi }_2)/2 are defined by the phase drops {\varphi }_n across the junction E_{j_n} ( n=1,2,3). The effective potential of flux qutrit and the additional Josephson junction is given by U_{\text{eff}} = E_J [−α₂ cos(f₁ − \varphi_p) − 2cos( \varphi_m/2) cos( \varphi_p/2)] + α₄ E_J cos(f₁+ f_2-\mathrm{\Delta }\mathrm{\Psi }+{\varphi }_p), where {\alpha }_4 is the ratio of the added Josephson junction energy E_4 to the larger Josephson junction energy E_J and f_n(n=1 ,2) reduced magnetic flux through the loop n. When f_1\ne 0.5, a flux qutrit holds a \mathrm{\Lambda }-type energy level structure [117,118]. When introducing the capacitive terms from the coupler capacitive C_{\text{coupler}}, a numerical evaluation of the Hamiltonian for f_2=\pi in the qutrit basis reveals that

(20)$$\begin{gathered} H_{\text{q}}\sim {\displaystyle \sum _{k=L,R,e}{\mathrm{\Delta }}_k|k⟩⟨k|+{\alpha }_4{E}_J\mathrm{\Delta }\mathrm{\Psi }{\displaystyle \sum _{m=L,R}{c}_r^1({\alpha }_2,{\alpha }_4,f_1){\sigma }_m,}} \\ \end{gathered}$$

where {\mathrm{\Delta }}_k the frequency of the energy |k⟩ in the qutrit A_n and {\sigma }_m=|m{⟩}_n⟨e|+|e{⟩}_n⟨m | ( m=L,R). c_m^1({\alpha }_2,{\alpha }_4,f_1 ) is liner in the field and has a tunable orientation. To be specific, the explicit expression of the function can be represented by c_r^1({\alpha }_2,{\alpha }_4,f_1)=\sin (f_1+{\varphi }_{+}), where {\varphi }_{+}= {\varphi }_1+{\varphi }_3 and f_1={\varphi }_1+{\varphi }_3-{\varphi }_2-{\varphi }_4. The phases {\varphi }_n across the junction E_n ( n=1,2,3,4) are strongly associated with the values of {\alpha }_2 and {\alpha }_4 [114]. Accordingly, we restrict to the case in which the line forms a single-mode resonator B_n and the phase slip then becomes approximately

(21)\begin{array}{rl}\mathrm{\Delta }\mathrm{\Psi }& =\frac{{\mathrm{\partial }}_{x}u(\mathit{x})\mathrm{\Delta }x}{{\varphi }_0}\sqrt{\frac{\hslash }{{\omega }_{b_n}C}}(b_n+b_n^{†})\\ & =\frac{2\pi {\mathrm{\partial }}_x\mathrm{\Psi }(\mathit{x})\mathrm{\Delta }x}{{\mathrm{\Phi }}_0}(b_n+b_n^{†}),\end{array}

where u(\mathit{x}) is the photon mode eigenfunction in the resonator, \mathrm{\Delta }x the separation between the two qutrit-line intersections and C the total transmission line capacitance. Assuming a flux gradient |{\mathrm{\partial }}_x\mathrm{\Psi }(\mathit{x})|=65\times {10}^{-5}{\mathrm{\Phi }}_0 µm and a qutrit size \mathrm{\Delta }x =5 µm, we reach a coupling J_0=2\times {10}^{-4}{E}_J, which for a typical junction with E_J=250 GHz implies a 50 MHz coupling. By substituting quantized magnetic flux and employing the rotating wave approximation in the qutrit-resonator chain, the quantized interaction Hamiltonian can be obtained in Eq. (1).

Notice that the coupling term is strictly independent of the qutrit Hamiltonian, so it now becomes possible to switch on and off the interaction. The simplest way to tune the coupling strength J_{1,2} is to apply a control magnetic flux to this loop dynamically with \mathrm{\Psi }(\mathit{x})={\int }_s\mathit{B}( \mathit{x},t)\cdot \mathrm{d}\mathit{S}, by adding the FBL to connect with an arbitrary waveforms generator (AWG) by adopting controlled voltage pulses [119], shown in Fig.14. Moreover, the coupling between the flux qutrit and the resonator induced by the capacitive energy of junctions gives a negligible coupling \sim {10}^{-3}\hslash {\omega }_{b_n} [111]. It is an expedient way to add one additional Josephson junction indicating that (i) the tunable couplings between the qutrit and the resonator can switch on and off by the parameter f_1 in contrast to no additional Josephson junction and (ii) the complexity of superconducting circuits and unexpected higher-order couplings can be reduced by adding less number of Josephson junctions. In addition, an appropriate flux difference increasing or decreasing f_1 is applied on the flux qutrit A_1 to unbalance the populations of the two current states and create a superposition state of (|R⟩+ | e⟩)/\sqrt{2} [120]. The resonator stays in the zero-photon Fock state. Finally, the qutrit-resonator chain can be prepared in 1/\sqrt{2}(|RLR\cdots R(L){⟩}_{A_N}+|eLR\cdots R(L) {⟩}_{A_N})\otimes | 000\cdots 0{⟩}_{B_{N-1}}.

In experiment, the above superconducting circuit device can be integrated by a metal chip, which is fabricated on a micron scale and operated in the dilution refrigerator at millikelvin temperatures [119,121-123]. These circuits are driven by currents, voltages and microwave photons that excite the system from one quantum state to another in a controllable manner, which can be used to test fundamental quantum mechanical principles at a macroscopic scale. As shown in Fig.14, the probe microwave signal is continuously sent from a network analyzer and attenuated in the signal input line before arriving at the sample, which is placed in a magnetic shield. The transmitted signal from the sample is amplified and measured by the network analyzer. Furthermore, the superconducting quantum computer on the order of 65-qubit has been realized [124]. The 72-site circuit QED lattice was demonstrated experimentally [30]. Thus, superconducting circuits possessing advantages of flexibility, scalability and tunability [18,125,126], provide an excellent platform for generating large-scale GHZ states with a high fidelity.

5 Conclusion

Summing up, we have proposed the scheme of the fast topological pumping via edge channels with respect to a single excitation to generate large-scale GHZ states in a generalized SSH model of a superconducting qutrit-resonator circuit. We analytically derive that this generalized SSH model can be expressed as a two-band structure in momentum space, and its topologically trivial and nontrivial phase can be characterized by Zak phase equal to 0 and \pi, respectively. There is a nontrivial band gap to protect fast topological pumping via edge channels. The crucial characteristics of exponential couplings for fast topological pumping are that it suitably adapts the slope of the coupling functions based on the value of the instantaneous energy gap, while at the same time it ensures that the minimum value of the energy gap is as greater as possible. We study the effect of the disorder and losses of system for generating large-scale GHZ states, emphasizing the fact that the robustness of large-scale GHZ state are increased. Furthermore, the accessible scalability of entanglement based on the fast topological pumping outperforms a lot the conventional adiabatic topological pumping. Last but not least, the multiple platforms for realization of the fast topological pumping in experimental setups may mitigate additional constraints and possibilities for the control with respect to the other accelerating methods [11-13,55]. The scheme provides a fast topological pumping to generate large-scale GHZ states with high fidelity and robustness, which is expected to make a substantial contribution to speeding up adiabatic protocols with the topological properties of matter in the superconducting circuit device.

6 Derivation of zero-mode edge state

For the odd-sized superconducting qutrit-resonator chain, the translational invariance of the system suggests the following ansatz for an eigenstate of edge state

\begin{array}{}\text{(text{A1})}& |{\phi }_E⟩={\displaystyle \sum _{n=1}^N{\lambda }^n[\gamma {\sigma }_n^{†}({\sigma }_n^x{)}^{n-1}\underset{l=1}{\overset{n-1}{⨂}}{\sigma }_l^x+\eta b_n^{†}\underset{l=1}{\overset{n}{⨂}}{\sigma }_l^x]|G{⟩}_L,}\end{array}

where {\sigma }_n^{+}=|e{⟩}_n⟨R| and |G{⟩}_L=|RLR\cdots {⟩}_{A_N}\otimes |000\cdots {⟩}_{B_{N-1}} denotes a decoupled state of the qutrit-resonator chain by mapping into the parameter coordinate space with A_1 in |R⟩, B_1 in |0⟩ (zero-photon Fock state), A_2 in |L⟩, B_2 in |0⟩, \cdots \lambda is the localized index, \gamma and \eta being the probability amplitude of the gap states. The probability amplitude on site n decays (increases) exponentially with the distance n when ( |\lambda |>1), corresponding to the left (right) edge state, after the wave function is normalized. It is supposed that the eigenenergy of an edge state is E. Through the eigenvalue equation E|{\phi }_E⟩=H_{2N-1}|{\phi }_E⟩, one can obtain

\begin{array}{}\text{(text{A2})}& \begin{array}{rl}& E\sum _{n=1}^N{\lambda }^n[\gamma {\sigma }_n^{†}({\sigma }_n^x{)}^{n-1}\underset{l=1}{\overset{n-1}{⨂}}{\sigma }_l^x+\eta b_n^{†}\underset{l=1}{\overset{n}{⨂}}{\sigma }_l^x]|G{⟩}_L\\ & =J_1{\lambda }^n[\eta {\sigma }_n^{†}({\sigma }_n^x{)}^n\underset{l=1}{\overset{n}{⨂}}{\sigma }_l^x+\gamma b_n^{†}\underset{l=1}{\overset{n}{⨂}}{\sigma }_l^x]|G{⟩}_L\\ & +J_2[\eta {\lambda }^{n-1}{\sigma }_n^{†}({\sigma }_n^x{)}^{n-1}\underset{l=1}{\overset{n-1}{⨂}}{\sigma }_l^x+\gamma {\lambda }^{n+1}{b}_j^{†}\underset{l=1}{\overset{n}{⨂}}{\sigma }_l^x]|G{⟩}_L.\end{array}\end{array}

When \gamma =0 ( \eta =0), according to Eq. (A1) the resonators (qutrits) are occupied by the edge state whose eigenenergy is E=0. In particular, in order to generate large-scale GHZ states of qutrits, we choose \gamma = 1 and \eta =0 to render the qutrit in each unit cell to occupy the left ( ) and right ( \lambda >1) edge states with E=0. Then, the edge state wave function can be derived as

\begin{array}{}\text{(text{A3})}& \begin{array}{r}|{\phi }_E⟩=\sum _{n=1}^N{\lambda }^n{\sigma }_n^{†}({\sigma }_n^x{)}^{n-1}\underset{l=1}{\overset{n-1}{⨂}}{\sigma }_l^x|G{⟩}_L,\end{array}\end{array}

with \lambda =- J_1/J_2.

References

[1] A.Y. Kitaev, Unpaired Majorana fermions in quantum wires, Phys. Uspekhi . 44(103), 131 ( 2001)

[2] C.L. Kane E.J. Mele, Z₂ topological order and the quantum spin Hall effect , Phys. Rev. Lett . 95(14), 146802 ( 2005)

[3] D.J. Thouless. Quantization of particle transport. Phys. Rev. B , 1983, 27( 10): 6083

[4] Q.Niu, D.J. Thouless. Quantised adiabatic charge transport in the presence of substrate disorder and many-body interaction. J. Phys. Math. Gen. , 1984, 17( 12): 2453

[5] Y.E. Kraus, Y.Lahini, Z.Ringel, M.Verbin, O.Zilberberg. Topological states and adiabatic pumping in quasicrystals. Phys. Rev. Lett. , 2012, 109( 10): 106402

[6] M.Verbin, O.Zilberberg, Y.Lahini, Y.E. Kraus, Y.Silberberg. Topological pumping over a photonic Fibonacci quasicrystal. Phys. Rev. B , 2015, 91( 6): 064201

[7] N.Lang H.P. Büchler, Topological networks for quantum commuication between distant qubits, npj Quantum Inf. 3, 47 ( 2017)

[8] F.Mei, G.Chen, L.Tian, S.L. Zhu, S.Jia. Robust quantum state transfer via topological edge states in superconducting qubit chains. Phys. Rev. A , 2018, 98( 1): 012331

[9] J.L. Tambasco, G.Corrielli, R.J. Chapman, A.Crespi, O.Zilberberg, R.Osellame, A.Peruzzo. Quantum interference of topological states of light. Sci. Adv. , 2018, 4( 9): eaat3187

[10] P.Boross, J.K. Asbóth, G.Széchenyi, L.Oroszlány, A.Pályi. Poor man’s topological quantum gate based on the Su−Schrieffer−Heeger model. Phys. Rev. B , 2019, 100( 4): 045414

[11] S.Longhi. Topological pumping of edge states via adiabatic passage. Phys. Rev. B , 2019, 99( 15): 155150

[12] N.E. Palaiodimopoulos, I.Brouzos, F.K. Diakonos, G.Theocharis. Fast and robust quantum state transfer via a topological chain. Phys. Rev. A , 2021, 103( 5): 052409

[13] F.M. D’Angelis, F.A. Pinheiro, D.Guéry-Odelin, S.Longhi, F.Impens. Fast and robust quantum state transfer in a topological Su−Schrieffer−Heeger chain with next-to-nearest-neighbor interactions. Phys. Rev. Res. , 2020, 2( 3): 033475

[14] F.Mei, G.Chen, L.Tian, S.L. Zhu, S.Jia. Topology-dependent quantum dynamics and entanglement-dependent topological pumping in superconducting qubit chains. Phys. Rev. A , 2018, 98( 3): 032323

[15] J.X. Han, J.L. Wu, Y.Wang, Y.Xia, Y.Y. Jiang, J.Song. Large-scale Greenberger−Horne−Zeilinger states through a topologically protected zero-energy mode in a superconducting qutrit-resonator chain. Phys. Rev. A , 2021, 103( 3): 032402

[16] S.Das Sarma, M.Freedman, C.Nayak. Topological quantum computation. Phys. Today , 2006, 59( 7): 32

[17] M.H. Devoret, R.J. Schoelkopf. Superconducting circuits for quantum information: An outlook. Science , 2013, 339( 6124): 1169

[18] J.Q. You F.Nori, Atomic physics and quantum optics using superconducting circuits, Nature 474(189), 8 ( 2005)

[19] J.M. Martinis, Qubit metrdogy for building a fault−tolerant quantum computer, npj Quantum Inf. 1, 15005 ( 2015)

[20] C.Song, K.Xu, H.Li, Y.R. Zhang, X.Zhang, W.Liu, Q.Guo, Z.Wang, W.Ren, J.Hao, H.Feng, H.Fan, D.Zheng, D.W. Wang, H.Wang, S.Y. Zhu. Generation of multicomponent atomic Schrödinger cat states of up to 20 qubits. Science , 2019, 365( 6453): 574

[21] J.Kelly, R.Barends, A.G. Fowler, A.Megrant, E.Jeffrey, T.C. White, D.Sank, J.Y. Mutus, B.Campbell, Y.Chen, Z.Chen, B.Chiaro, A.Dunsworth, I.C. Hoi, C.Neill, P.J. J. O’Malley, C.Quintana, P.Roushan, A.Vainsencher, J.Wenner, A.N. Cleland, J.M. Martinis. State preservation by repetitive error detection in a superconducting quantum circuit. Nature , 2015, 519( 7541): 66

[22] J.S. Otterbach R.Manenti N.Alidoust A.Bestwick M.Block B.Bloom S.Caldwell N.Didier E.S. Fried S.Hong P.Karalekas C.B. Osborn A.Papageorge E.C. Peterson G.Prawiroatmodjo N.Rubin C.A. Ryan D.Scarabelli M.Scheer E.A. Sete P.Sivarajah R.S. Smith N.T. A. Staley W.J. Zeng A.Hudson B.R. Johnson M.Reagor M.P. da Silva C.Rigetti, Unsupervised machine learning on a hybrid quantum computer, arXiv: 1712.05771 ( 2017)

[23] A.Kandala, A.Mezzacapo, K.Temme, M.Takita, M.Brink, J.M. Chow, J.M. Gambetta. Hardware-efficient variational quantum eigensolver for small molecules and quantum magnets. Nature , 2017, 549( 7671): 242

[24] C.Neill P.Roushan K.Kechedzhi S.Boixo S.V. Isakov V.Smelyanskiy A.Megrant B.Chiaro A.Dunsworth K.Arya R.Barends B.Burkett Y.Chen Z.Chen A.Fowler B.Foxen M.Giustina R.Graff E.Jeffrey T.Huang J.Kelly P.Klimov E.Lucero J.Mutus M.Neeley C.Quintana D.Sank A.Vainsencher J.Wenner T.C. White H.Neven J.M. Martinis, A blueprint for demonstrating quantum supremacy with superconducting qubits, Science 360(6385), 195 ( 2018)

[25] K.X. Wei, I.Lauer, S.Srinivasan, N.Sundaresan, D.T. McClure, D.Toyli, D.C. McKay, J.M. Gambetta, S.Sheldon. Verifying multipartite entangled Greenberger−Horne−Zeilinger states via multiple quantum coherences. Phys. Rev. A , 2020, 101( 3): 032343

[26] Y.Salathé, M.Mondal, M.Oppliger, J.Heinsoo, P.Kurpiers, A.Potočnik, A.Mezzacapo, U.Las Heras, L.Lamata, E.Solano, S.Filipp, A.Wallraff. Digital quantum simulation of spin models with circuit quantum electrodynamics. Phys. Rev. X , 2015, 5( 2): 021027

[27] S.Hacohen-Gourgy, V.V. Ramasesh, C.De Grandi, I.Siddiqi, S.M. Girvin. Cooling and autonomous feedback in a Bose−Hubbard chain with attractive interactions. Phys. Rev. Lett. , 2015, 115( 24): 240501

[28] P.J. J. O’Malley, R.Babbush, I.D. Kivlichan, J.Romero, J.R. McClean, R.Barends, J.Kelly, P.Roushan, A.Tranter, N.Ding, B.Campbell, Y.Chen, Z.Chen, B.Chiaro, A.Dunsworth, A.G. Fowler, E.Jeffrey, E.Lucero, A.Megrant, J.Y. Mutus, M.Neeley, C.Neill, C.Quintana, D.Sank, A.Vainsencher, J.Wenner, T.C. White, P.V. Coveney, P.J. Love, H.Neven, A.Aspuru-Guzik, J.M. Martinis. Scalable quantum simulation of molecular energies. Phys. Rev. X , 2016, 6( 3): 031007

[29] Y.P. Zhong, D.Xu, P.Wang, C.Song, Q.J. Guo, W.X. Liu, K.Xu, B.X. Xia, C.Y. Lu, S.Han, J.W. Pan, H.Wang. Emulating anyonic fractional statistical behavior in a superconducting quantum circuit. Phys. Rev. Lett. , 2016, 117( 11): 110501

[30] M.Fitzpatrick, N.M. Sundaresan, A.C. Y. Li, J.Koch, A.A. Houck. Observation of a dissipative phase transition in a one-dimensional circuit QED lattice. Phys. Rev. X , 2017, 7( 1): 011016

[31] J.Q. You, X.F. Shi, X.Hu, F.Nori. Quantum emulation of a spin system with topologically protected ground states using superconducting quantum circuits. Phys. Rev. B , 2010, 81( 1): 014505

[32] J.Koch, A.A. Houck, K.L. Hur, S.M. Girvin. Time−reversal-symmetry breaking in circuit-QED-based photon lattices. Phys. Rev. A , 2010, 82( 4): 043811

[33] M.Hafezi, P.Adhikari, J.M. Taylor. Engineering three-body interaction and Pfaffian states in circuit QED systems. Phys. Rev. B , 2014, 90( 6): 060503

[34] E.Kapit, M.Hafezi, S.H. Simon. Induced self-stabilization in fractional quantum Hall states of light. Phys. Rev. X , 2014, 4( 3): 031039

[35] D.I. Tsomokos, S.Ashhab, F.Nori. Using superconducting qubit circuits to engineer exotic lattice systems. Phys. Rev. A , 2010, 82( 5): 052311

[36] F.Mei, Z.Y. Xue, D.W. Zhang, L.Tian, C.Lee, S.L. Zhu. Witnessing topological Weyl semimetal phase in a minimal circuit-QED lattice. Quantum Sci. Technol. , 2016, 1( 1): 015006

[37] J.Tangpanitanon, V.M. Bastidas, S.Al-Assam, P.Roushan, D.Jaksch, D.G. Angelakis. Topological pumping of photons in nonlinear resonator arrays. Phys. Rev. Lett. , 2016, 117( 21): 213603

[38] T.Goren, K.Plekhanov, F.Appas, K.Le Hur. Topological Zak phase in strongly coupled LC circuits. Phys. Rev. B , 2018, 97( 4): 041106

[39] V.V. Ramasesh, E.Flurin, M.Rudner, I.Siddiqi, N.Y. Yao. Direct probe of topological invariants using Bloch oscillating quantum walks. Phys. Rev. Lett. , 2017, 118( 13): 130501

[40] E.Flurin, V.V. Ramasesh, S.Hacohen-Gourgy, L.S. Martin, N.Y. Yao, I.Siddiqi. Observing topological invariants using quantum walks in superconducting circuits. Phys. Rev. X , 2017, 7( 3): 031023

[41] M.D. Schroer, M.H. Kolodrubetz, W.F. Kindel, M.Sandberg, J.Gao, M.R. Vissers, D.P. Pappas, A.Polkovnikov, K.W. Lehnert. Measuring a topological transition in an artificial spin-1/2 system. Phys. Rev. Lett. , 2014, 113( 5): 050402

[42] P.Roushan C.Neill Y.Chen M.Kolodrubetz C.Quintana N.Leung M.Fang R.Barends B.Campbell Z.Chen B.Chiaro A.Dunsworth E.Jeffrey J.Kelly A.Megrant J.Mutus P.J. J. O’Malley D.Sank A.Vainsencher J.Wenner T.White A.Polkovnikov A.N. Cleland J.M. Martinis, Observation of topological transitions in interacting quantum circuits, Nature 515(7526), 241 ( 2014)

[43] Z.Zhang, T.Wang, L.Xiang, J.Yao, J.Wu, Y.Yin. Measuring the Berry phase in a superconducting phase qubit by a shortcut to adiabaticity. Phys. Rev. A , 2017, 95( 4): 042345

[44] X.Tan, D.W. Zhang, Q.Liu, G.Xue, H.F. Yu, Y.Q. Zhu, H.Yan, S.L. Zhu, Y.Yu. Topological Maxwell metal bands in a superconducting qutrit. Phys. Rev. Lett. , 2018, 120( 13): 130503

[45] C.Song, D.Xu, P.Zhang, J.Wang, Q.Guo, W.Liu, K.Xu, H.Deng, K.Huang, D.Zheng, S.B. Zheng, H.Wang, X.Zhu, C.Y. Lu, J.W. Pan. Demonstration of topological robustness of anyonic braiding statistics with a superconducting quantum circuit. Phys. Rev. Lett. , 2018, 121( 3): 030502

[46] W.Cai, J.Han, F.Mei, Y.Xu, Y.Ma, X.Li, H.Wang, Y.P. Song, Z.Y. Xue, Z.Q. Yin, S.Jia, L.Sun. Observation of topological magnon insulator states in a superconducting circuit. Phys. Rev. Lett. , 2019, 123( 8): 080501

[47] O.Viyuela A.Rivas S.Gasparinetti A.Wallraff A.Wallraff S.Filipp M.A. Martin-Delgado, Observation of topological Uhlmann phases with superconduction qubits, npj Quantum Inform. 4, 10 ( 2018)

[48] R.Horodecki, P.Horodecki, M.Horodecki, K.Horodecki. Quantum entanglement. Rev. Mod. Phys. , 2009, 81( 2): 865

[49] J.W. Pan, Z.B. Chen, C.Y. Lu, H.Weinfurter, A.Zeilinger, M.Żukowski. Multiphoton entanglement and interferometry. Rev. Mod. Phys. , 2012, 84( 2): 777

[50] L.Pezzè, A.Smerzi, M.K. Oberthaler, R.Schmied, P.Treutlein. Quantum metrology with nonclassical states of atomic ensembles. Rev. Mod. Phys. , 2018, 90( 3): 035005

[51] M.C. Rechtsman, Y.Lumer, Y.Plotnik, A.Perez-Leija, A.Szameit, M.Segev. Topological protection of photonic path entanglement. Optica , 2016, 3( 9): 925

[52] A.Blanco-Redondo, B.Bell, D.Oren, B.J. Eggleton, M.Segev. Topological protection of biphoton states. Science , 2018, 362( 6414): 568

[53] M.Wang, C.Doyle, B.Bell, M.J. Collins, E.Magi, B.J. Eggleton, M.Segev, A.Blanco-Redondo. Topologically protected entangled photonic states. Nanophotonics , 2019, 8( 8): 1327

[54] C.K. Hong, Z.Y. Ou, L.Mandel. Measurement of subpicosecond time intervals between two photons by interference. Phys. Rev. Lett. , 1987, 59( 18): 2044

[55] I.Brouzos, I.Kiorpelidis, F.K. Diakonos, G.Theocharis. Fast, robust, and amplified transfer of topological edge modes on a time-varying mechanical chain. Phys. Rev. B , 2020, 102( 17): 174312

[56] Y.X. Shen, L.S. Zeng, Z.G. Geng, D.G. Zhao, Y.G. Peng, X.F. Zhu. Acoustic adiabatic propagation based on topological pumping in a coupled multicavity chain lattice. Phys. Rev. Appl. , 2020, 14( 1): 014043

[57] J.L. Wu, Y.Wang, J.X. Han, Y.K. Feng, S.L. Su, Y.Xia, Y.Jiang, J.Song. One-step implementation of Rydberg-antiblockade SWAP and controlled-SWAP gates with modified robustness. Photon. Res. , 2021, 9( 5): 814

[58] X.R. Huang, Z.X. Ding, C.S. Hu, L.T. Shen, W.Li, H.Wu, S.B. Zheng. Robust Rydberg gate via Landau−Zener control of Förster resonance. Phys. Rev. A , 2018, 98( 5): 052324

[59] Q.Guo, S.B. Zheng, J.Wang, C.Song, P.Zhang, K.Li, W.Liu, H.Deng, K.Huang, D.Zheng, X.Zhu, H.Wang, C.Y. Lu, J.W. Pan. Dephasing-insensitive quantum information storage and processing with superconducting qubits. Phys. Rev. Lett. , 2018, 121( 13): 130501

[60] V.Balachandran, J.Gong. Adiabatic quantum transport in a spin chain with a moving potential. Phys. Rev. A , 2008, 77( 1): 012303

[61] J.Allcock, N.Linden. Quantum communication beyond the localization length in disordered spin chains. Phys. Rev. Lett. , 2009, 102( 11): 110501

[62] J.L. Wu, Y.Wang, J.X. Han, S.L. Su, Y.Xia, Y.Jiang, J.Song. Unselective ground-state blockade of Rydberg atoms for implementing quantum gates. Front. Phys. , 2022, 17( 2): 22501

[63] R.R. Agundez, C.D. Hill, L.C. L. Hollenberg, S.Rogge, M.Blaauboer. Superadiabatic quantum state transfer in spin chains. Phys. Rev. A , 2017, 95( 1): 012317

[64] J.L. Wu, Y.Wang, J.X. Han, Y.Jiang, J.Song, Y.Xia, S.L. Su, W.Li. Systematic-error-tolerant multiqubit holonomic entangling gates. Phys. Rev. Appl. , 2021, 16( 6): 064031

[65] Y.Zhou, D.Y. Lü, W.Y. Zeng. Chiral single-photon switch-assisted quantum logic gate with a nitrogen-vacancy center in a hybrid system. Photon. Res. , 2021, 9( 3): 405

[66] D.M. Greenberger M.A. Horne A.Zeilinger, Bell’s theorem, Quantum Theory, and Conceptions of the Universe, Kluwer Dordrecht, 1989

[67] M.Hillery, V.Bužek, A.Berthiaume. Quantum secret sharing. Phys. Rev. A , 1999, 59( 3): 1829

[68] S.Bose, V.Vedral, P.L. Knight. Multiparticle generalization of entanglement swapping. Phys. Rev. A , 1998, 57( 2): 822

[69] E.Knill. Quantum computing with realistically noisy devices. Nature , 2005, 434( 7029): 39

[70] V.Giovannetti, S.Lloyd, L.Maccone. Quantum-enhanced measurements: Beating the standard quantum limit. Science , 2004, 306( 5700): 1330

[71] D.Leibfried, M.Barrett, T.Schaetz, J.Britton, J.Chiaverini, W.Itano, J.Jost, C.Langer, D.Wineland. Toward Heisenberg-limited spectroscopy with multiparticle entangled states. Science , 2004, 304( 5676): 1476

[72] C.P. Yang, Q.P. Su, S.B. Zheng, F.Nori. Entangling superconducting qubits in a multi-cavity system. New J. Phys. , 2016, 18( 1): 013025

[73] S.Matsuo S.Ashhab T.Fujii F.Nori K.Nagai N.Hatakenaka, Generation of Bell states and Greenberger−Horne−Zeilinger states in superconducting phase qubits, in: Quantum Communication, Measurement and Computing, No. 8, Ed.: O. Hirota et al. , Tokyo: NICT, 2006

[74] L.F. Wei, Y.Liu, F.Nori. Generation and control of Greenberger−Horne−Zeilinger entanglement in superconducting circuits. Phys. Rev. Lett. , 2006, 96( 24): 246803

[75] S.L. Zhu, Z.D. Wang, P.Zanardi. Geometric quantum computation and multiqubit entanglement with superconducting qubits inside a cavity. Phys. Rev. Lett. , 2005, 94( 10): 100502

[76] C.P. Yang, Q.P. Su, S.Han. Generation of Greenberger−Horne−Zeilinger entangled states of photons in multiple cavities via a superconducting qutrit or an atom through resonant interaction. Phys. Rev. A , 2012, 86( 2): 022329

[77] C.P. Yang, Q.P. Su, S.B. Zheng, S.Han. Generating entanglement between microwave photons and qubits in multiple cavities coupled by a superconducting qutrit. Phys. Rev. A , 2013, 87( 2): 022320

[78] S.Aldana, Y.D. Wang, C.Bruder. Greenberger−Horne−Zeilinger generation protocol forN superconducting transmon qubits capacitively coupled to a quantum bus. Phys. Rev. B , 2011, 84( 13): 134519

[79] W.Feng, P.Wang, X.Ding, L.Xu, X.Q. Li. Generating and stabilizing the Greenberger−Horne−Zeilinger state in circuit QED: Joint measurement, Zeno effect, and feedback. Phys. Rev. A , 2011, 83( 4): 042313

[80] J.L. Wu, C.Song, J.Xu, L.Yu, X.Ji, S.Zhang. Adiabatic passage for one-step generation of n-qubit Greenberger–Horne–Zeilinger states of superconducting qubits via quantum Zeno dynamics. Quantum Inform. Process. , 2016, 15( 9): 3663

[81] X.T. Mo, Z.Y. Xue. Single-step multipartite entangled states generation from coupled circuit cavities. Front. Phys. , 2019, 14( 3): 31602

[82] Y.H. Kang, Z.C. Shi, B.H. Huang, J.Song, Y.Xia. Deterministic conversions between Greenberger−Horne−Zeilinger states and W states of spin qubits via Lie-transform-based inverse Hamiltonian engineering. Phys. Rev. A , 2019, 100( 1): 012332

[83] T.Liu, Q.P. Su, Y.Zhang, Y.L. Fang, C.P. Yang. Generation of quantum entangled states of multiple groups of qubits distributed in multiple cavities. Phys. Rev. A , 2020, 101( 1): 012337

[84] J.M. Chow, J.M. Gambetta, E.Magesan, D.W. Abraham, A.W. Cross, B.R. Johnson, N.A. Masluk, C.A. Ryan, J.A. Smolin, S.J. Srinivasan, M.Steffen. Implementing a strand of a scalable fault-tolerant quantum computing fabric. Nat. Commun. , 2014, 5( 1): 4015

[85] L.DiCarlo M.D. Reed L.Sun B.R. Johnson J.M. Chow J.M. Gambetta L.Frunzio M.H. Girvin, S. M. and Devoret, and R. J. Schoelkopf, Preparation and measurement of three-qubit entanglement in a superconduction circuit, Nature 467(7315), 574 ( 2010)

[86] R.Barends, J.Kelly, A.Megrant, A.Veitia, D.Sank, E.Jeffrey, T.C. White, J.Mutus, A.G. Fowler, B.Campbell, Y.Chen, Z.Chen, B.Chiaro, A.Dunsworth, C.Neill, P.O’Malley, P.Roushan, A.Vainsencher, J.Wenner, A.N. Korotkov, A.N. Cleland, J.M. Martinis. Superconducting quantum circuits at the surface code threshold for fault tolerance. Nature , 2014, 508( 7497): 500

[87] C.Song K.Xu W. Liu C.Yang S.B. Zheng H.Deng Q.Xie K.Huang Q.Guo L.Zhang P.Zhang D.Xu D. Zheng X.Zhu H.Wang Y.A. Chen C.Y. Lu S.Han J.W. Pan, 10-qubit entanglement and parallel logic operations with a superconducting circuit, Phys. Rev. Lett . 119(18), 180511 ( 2017)

[88] A.Cervera-Lierta, M.Krenn, A.Aspuru-Guzik, A.Galda. Experimental high-dimensional Greenberger−Horne−Zeilinger entanglement with superconducting transmon qutrits. Phys. Rev. Appl. , 2022, 17( 2): 024062

[89] X.Chen, I.Lizuain, A.Ruschhaupt, D.Guéry-Odelin, J.G. Muga. Shortcut to adiabatic passage in two- and three-level atoms. Phys. Rev. Lett. , 2010, 105( 12): 123003

[90] D.Guéry-Odelin, A.Ruschhaupt, A.Kiely, E.Torrontegui, S.Martínez-Garaot, J.G. Muga. Shortcuts to adiabaticity: Concepts, methods, and applications. Rev. Mod. Phys. , 2019, 91( 4): 045001

[91] A.Altland, M.R. Zirnbauer. Nonstandard symmetry classes in mesoscopic normal-superconducting hybrid structures. Phys. Rev. B , 1997, 55( 2): 1142

[92] S.Ryu, Y.Hatsugai. Topological origin of zero-energy edge states in particle−hole symmetric systems. Phys. Rev. Lett. , 2002, 89( 7): 077002

[93] J.K. Asbóth L.Oroszlány A.Pályi, A short course on topological insulators, Lect. Notes Phys . 919, 85 ( 2016)

[94] A.Coutant, V.Achilleos, O.Richoux, G.Theocharis, V.Pagneux. Robustness of topological corner modes against disorder with application to acoustic networks. Phys. Rev. B , 2020, 102( 21): 214204

[95] D.M. Greenberger, M.A. Horne, A.Shimony, A.Zeilinger. Bell’s theorem without inequalities. Am. J. Phys. , 1990, 58( 12): 1131

[96] R.Mukherjee, H.Xie, F.Mintert. Bayesian optimal control of Greenberger−Horne−Zeilinger states in Rydberg lattices. Phys. Rev. Lett. , 2020, 125( 20): 203603

[97] F.Reiter, D.Reeb, A.S. Sørensen. Scalable dissipative preparation of many-body entanglement. Phys. Rev. Lett. , 2016, 117( 4): 040501

[98] S.B. Zheng. One-step synthesis of multiatom Greenberger−Horne−Zeilinger states. Phys. Rev. Lett. , 2001, 87( 23): 230404

[99] P.Mundada, G.Zhang, T.Hazard, A.Houck. Suppression of qubit crosstalk in a tunable coupling superconducting circuit. Phys. Rev. Appl. , 2019, 12( 5): 054023

[100] W.P. Su, J.R. Schrieffer, A.J. Heeger. Solitons in Polyacetylene. Phys. Rev. Lett. , 1979, 42( 25): 1698

[101] J.Pachos, H.Walther. Quantum computation with trapped ions in an optical cavity. Phys. Rev. Lett. , 2002, 89( 18): 187903

[102] A.P. M. Place, L.V. H. Rodgers, P.Mundada, B.M. Smitham, M.Fitzpatrick, Z.Leng, A.Premkumar, J.Bryon, A.Vrajitoarea, S.Sussman, G.Cheng, T.Madhavan, H.K. Babla, X.H. Le, Y.Gang, B.Jäck, A.Gyenis, N.Yao, R.J. Cava, N.P. de Leon, A.A. Houck. New material platform for superconducting transmon qubits with coherence times exceeding 0.3 milliseconds. Nat. Commun. , 2021, 12( 1): 1779

[103] H.Zhang, S.Chakram, T.Roy, N.Earnest, Y.Lu, Z.Huang, D.K. Weiss, J.Koch, D.I. Schuster. Universal fast-flux control of a coherent, low-frequency qubit. Phys. Rev. X , 2021, 11( 1): 011010

[104] I.M. Pop, M.Ansmann, G.Catelani, R.J. Schoelkopf, L.I. Glazman, M.H. Devoret. Coherent suppression of electromagnetic dissipation due to superconducting quasiparticles. Nature , 2014, 508( 7496): 7496

[105] X.Gu, A.F. Kockum, A.Miranowicz, Y.X. Liu, F.Nori. Microwave photonics with superconducting quantum circuits. Phys. Rep. , 2017, 718– 719,1

[106] M.Reagor, H.Paik, G.Catelani, L.Sun, C.Axline, E.Holland, I.M. Pop, N.A. Masluk, T.Brecht, L.Frunzio, M.H. Devoret, L.Glazman, R.J. Schoelkopf. Reaching 10 ms single photon lifetimes for superconducting aluminum cavities. Appl. Phys. Lett. , 2013, 102( 19): 192604

[107] M.Reagor, W.Pfaff, C.Axline, R.W. Heeres, N.Ofek, K.Sliwa, E.Holland, C.Wang, J.Blumoff, K.Chou, M.J. Hatridge, L.Frunzio, M.H. Devoret, L.Jiang, R.J. Schoelkopf. Quantum memory with millisecond coherence in circuit QED. Phys. Rev. B , 2016, 94( 1): 014506

[108] C.Axline, M.Reagor, R.Heeres, P.Reinhold, C.Wang, K.Shain, W.Pfaff, Y.Chu, L.Frunzio, R.J. Schoelkopf. An architecture for integrating planar and 3D cQED devices. Appl. Phys. Lett. , 2016, 109( 4): 042601

[109] J.E. Mooij, T.P. Orlando, L.Levitov, L.Tian, C.H. van der Wal, S.Lloyd. Josephson persistent-current qubit. Science , 1999, 285( 5430): 1036

[110] C.H. van der Wal, A.C. J. ter Haar, F.K. Wilhelm, R.N. Schouten, C.J. P. M. Harmans, T.P. Orlando, S.Lloyd, J.E. Mooij. Quantum superposition of macroscopic persistent-current states. Science , 2000, 290( 5492): 773

[111] B.Peropadre, P.Forn-Díaz, E.Solano, J.J. García-Ripoll. Switchable ultrastrong coupling in circuit QED. Phys. Rev. Lett. , 2010, 105( 2): 023601

[112] M.S. Allman F.Altomare J.D. Whittaker K.Cicak D.Li A. Sirois J.Strong J.D. Teufel R.W. Simmonds, RF-SQUID-mediated coherent tunable coupling between a superconducting phase qubit and a lumped-element resonator, Phys. Rev. Lett . 104(17), 177004 ( 2010)

[113] M.S. Allman, J.D. Whittaker, M.Castellanos-Beltran, K.Cicak, F.da Silva, M.P. DeFeo, F.Lecocq, A.Sirois, J.D. Teufel, J.Aumentado, R.W. Simmonds. Tunable resonant and nonresonant interactions between a phase qubit and LC resonator. Phys. Rev. Lett. , 2014, 112( 12): 123601

[114] J.Bourassa, J.M. Gambetta, A.A. Abdumalikov, O.Astafiev, Y.Nakamura, A.Blais. Ultrastrong coupling regime of cavity QED with phase-biased flux qubits. Phys. Rev. A , 2009, 80( 3): 032109

[115] T.P. Orlando, J.E. Mooij, L.Tian, C.H. van der Wal, L.S. Levitov, S.Lloyd, J.J. Mazo. Superconducting persistent-current qubit. Phys. Rev. B , 1999, 60( 22): 15398

[116] Y.X. Liu, J.Q. You, L.F. Wei, C.P. Sun, F.Nori. Optical selection rules and phase-dependent adiabatic state control in a superconducting quantum circuit. Phys. Rev. Lett. , 2005, 95( 8): 087001

[117] M.Devoret B.Huard R.Schoelkopf L.F. Cugliandolo, Quantum Machines: Measurement and Control of Engineered Quantum Systems, Oxford University Press, USA, 2014

[118] V.E. Manucharyan, J.Koch, L.I. Glazman, M.H. Devoret. Fluxonium: Single Cooper-pair circuit free of charge offsets. Science , 2009, 326( 5949): 113

[119] L.DiCarlo, J.M. Chow, J.M. Gambetta, L.S. Bishop, B.R. Johnson, D.I. Schuster, J.Majer, A.Blais, L.Frunzio, S.M. Girvin, R.J. Schoelkopf. Demonstration of two-qubit algorithms with a superconducting quantum processor. Nature , 2009, 460( 7252): 240

[120] T.Wang, Z.Zhang, L.Xiang, Z.Jia, P.Duan, Z.Zong, Z.Sun, Z.Dong, J.Wu, Y.Yin, G.Guo. Experimental realization of a fast controlled-Z gate via a shortcut to adiabaticity. Phys. Rev. Appl. , 2019, 11( 3): 034030

[121] A.A. Clerk, K.W. Lehnert, P.Bertet, J.R. Petta, Y.Nakamura. Hybrid quantum systems with circuit quantum electrodynamics. Nat. Phys. , 2020, 16( 3): 257

[122] T.Niemczyk, F.Deppe, H.Huebl, E.P. Menzel, F.Hocke, M.J. Schwarz, J.J. Garcia-Ripoll, D.Zueco, T.Hümmer, E.Solano, A.Marx, R.Gross. Circuit quantum electrodynamics in the ultrastrong-coupling regime. Nat. Phys. , 2010, 6( 10): 772

[123] A.Blais, S.M. Girvin, W.D. Oliver. Quantum information processing and quantum optics with circuit quantum electrodynamics. Nat. Phys. , 2020, 16( 3): 247

[124] G.J. Mooney, G.A. L. White, C.D. Hill, L.C. L. Hollenberg. Whole-device entanglement in a 65-qubit superconducting quantum computer. Adv. Quantum Technol. , 2021, 4( 10): 2100061

[125] A.Blais, R.S. Huang, A.Wallraff, S.M. Girvin, R.J. Schoelkopf. Cavity quantum electrodynamics for superconducting electrical circuits: An architecture for quantum computation. Phys. Rev. A , 2004, 69( 6): 062320

[126] J.Clarke, F.K. Wilhelm. Superconducting quantum bits. Nature , 2008, 453( 7198): 1031

RIGHTS & PERMISSIONS

Higher Education Press