Arxiv:2003.06086V2 [Quant-Ph]

Home , Phase qubit

arXiv:2003.06086v2 [quant-ph] 23 Sep 2020 in qatmgts ymaso unu algorithms quantum of means [ the rota- by qubit gates) encodes emulates of (quantum cycles DQS and tions repeated qubits that through onto evolution fact time systems the the quantum the of of result uni- state This a applicable. is universally versality being of advantage the iistetm vlto foeseicmdlHamil- model specific one [ of tonian evolution time the mimics in tnwhsbe ieyapidt td quan- study to applied widely been has [ chemistry now tum It Hamilto- local nian. finite-dimensional any simulate efficiently ftplgclmte ihapormal unu pro- quantum programmable a with matter topological of parameter [ single-qubit topological space in a reported simulating been digitally has Hamiltonian topology, to regard In [ physics energy the [ of platforms advantages these taking of 11 flexibility materials, and controllability real high in what beyond possible go is to opportunities quan- provide materials, simulators solid-state tum quantum in discovered topological en- first Although were oppor- states matters, devices. exciting new of and for research tunities phases physics topological fascinating [ novel abling effect of Hall ploration quantum the and rset fpoigeoi oooia hss[ phases topological exotic probing of prospects 40 Introduction nsakcnrs oaao unu iuainthat simulation quantum analog to contrast stark In .Adsc ihtnblt ol ral nac the enhance greatly could tunability high such And ]. – iia iuaino oooia atro rgambeQ Programmable on Matter Topological of Simulation Digital 45 .Sc ici-ae iuao a,i principle, in can, simulator circuit-based a Such ]. 65 40 egMei, Feng .Hwvr o opromdgtlsimulation digital perform to how However, ]. , 41 ttsb esrn ui xiaindsrbto tteou the at distribution excitation qubit measuring by states nteeprmns entol bev u lodistinguish also but observe only not both we on b experiments, performed can the is In matter circuit quantum topological topological genera related low-depth current and the coupling on spin-orbital demo circuits and quantum propose protected phase we logically gates, topological quantum of elementary simulation two-qubit and quantum enhanced. digital quan greatly realize digital is to phases of universality topological exotic the simulating Given interest. siderable ,dgtlqatmsmlto DS has (DQS) simulation quantum digital ], icvre ftplgclisltr [ insulators topological of Discoveries . iuaigtetplgclpae fmte nsnhtcqu synthetic in matter of phases topological the Simulating 56 46 4 6 undn rvnilKyLbrtr fQatmEngineerin Quantum of Laboratory Key Provincial Guangdong undn rvnilKyLbrtr fQatmEngineerin Quantum of Laboratory Key Provincial Guangdong – , ,2, 1, 62 nttt fLsrSetocp,Sax nvriy Taiyu University, Shanxi Spectroscopy, Laser of Institute 47 n vngoercpae [ phases geometric even and ] ,mn-oymdl [ models many-body ], ∗ 3 colo cec,Xa ioogUiest,Xa 109 S 710049, Xian University, Jiaotong Xian Science, of School 1 ia Guo, Qihao tt e aoaoyo unu pisadQatmOtc De Optics Quantum and Optics Quantum of Laboratory Key State colo nomto n polcrncSineadEngin and Science Optoelectronic and Information of School colo hsc,NnigUiest,Nnig209,Chi 210093, Nanjing University, Nanjing Physics, of School ot hn omlUiest,Gaghu500,China 510006, Guangzhou University, Normal China South ot hn omlUiest,Gaghu500,China 510006, Guangzhou University, Normal China South 3 PT etrfrQatmPeiinMeasurement, Precision Quantum for Center GPETR aebotdteex- the boosted have ] 2 5 hniUiest,Tiun hni000,China 030006, Shanxi Taiyuan, University, Shanxi rnirRsac nttt o hsc n SPTE, and Physics for Institute Research Frontier olbrtv noainCne fEteeOptics, Extreme of Center Innovation Collaborative ainlLbrtr fSldSaeMicrostructures, State Solid of Laboratory National 3 aFiYu, Ya-Fei 48 Dtd etme 4 2020) 24, September (Dated: – 55 12 4 63 high- ], – inunXiao, Liantuan 39 , 1 64 ]. , 4 2 ]. – ] forrslsaetnbewt os unu circuits, protection. quantum topological noisy of power with the tenable All manifests are which phases. results insulator our SOCs Chern 2D of of topological DQS for anomalous circuit and quantum a design to method unknown. remains cessor h oooia hsso atr[ induce matter to of mechanisms phases prominent the topological of electric the one intrinsic is crystal’s and a field, in electrons of movement the matter [ reported distinguish previously and been observe not we and circuit, 0 the the of output the distribution at excitation qubit the measuring processors. by quantum Strikingly, Rigetti and the IBM on available phases currently topological to (1D) circuit quantum one-dimensional simulate topological low-depth experimen- a we perform Furthermore, measure tally states. to edge used topological and including invariants be matter, topological could also of We hallmarks circuits the to quantum matter. component that of essential phases show the topological quan- various providing two-qubit proto- generate thus and the single- gates, in common tum simulated employing digitally by topo- be col related can and matter (SOC) pro- logical coupling quantum Spin-orbital noisy of generation cessors. current quantum the protected on circuits topologically designing of technique unu icisfrSCadterltdtopological related the and SOC for circuits Quantum unique a demonstrate and develop we Letter, this In ,2 1, u iuaos h rseto digitally of prospect the simulators, tum oee,i ssila pnqeto how question open an still is it However, fmte.Hr,uigcmo single- common using Here, matter. of s nsldsaemtras O rgntsfrom originates SOC materials, solid-state In . ino os unu rcsoswhere processors quantum noisy of tion pto h circuits. the of tput iial iuae.I atclr a particular, In simulated. digitally e h-in Zhu, Shi-Liang B n iet unu processors. quantum Rigetti and IBM h and 0 the nu iuaosi oi fcon- of topic a is simulators antum π srt napoc odsg topo- design to approach an nstrate n hni000,China 030006, Shanxi an, nrytplgcleg tts hc have which states, edge topological energy n unu Materials, Quantum and g n unu Materials, Quantum and g π ani China haanxi, nrytplgcledge topological energy ,6, 5, eering, na atmProcessors uantum † vices, n utn Jia Suotang and 66 .Fnly eapyour apply we Finally, ]. 1 , 2 , 13 , 67 ,2 1, .I our In ]. 2

−iH˜st (a ) operator is expressed as U(t) = e . Note that the Q Q Q Q Q Q 1 2 2x- 1 2x 2N - 1 2N excitation number in the qubit array is conserved, i.e., x,• x,¯ the time evolution only works in the single excitation

(b ) Q + - + - subspace if the initial input state of the circuit is a single 2x Z Y (q ) Y (q ) Z U (q ) = excitation state. The single excitation subspace for the two qubits (Q2x,Q2x+1) is {|gei, |egi}. In this subspace, Q 2x+ 1 the evolution operator can be rewritten as

FIG. 1: (a) Mapping a 1D spinful lattice into qubits in a quan- U(θ) = cos(θ)(|geiQ2x,Q2x+1 hge| + |egiQ2x,Q2x+1 heg|) tum circuit. A pair of qubits (Q2x−1,Q2x) are employed to −i sin(θ)(|gei heg| + |egi hge|), simulate a spin particle at a lattice site (x↑,x↓). (b) Quantum Q2x,Q2x+1 Q2x,Q2x+1 circuit for implementing two-site SOC. where θ = Jst. We find that such a time evolution operation can be precisely programmed through a sequence of work, we shall present a unique approach to implement four two-qubit controlled-NOT (CNOT) gates and four DQS of SOC using programmable quantum circuits. Be- single-qubit rotation gates, i.e., fore illustrating our procedure in detail, we summarize its underlying mechanism: a spinful lattice with N sites is + − U(θ) = CNOTQ ,Q · Z · Y (θ) · CNOTQ ,Q mapped into a qubit array with 2N qubits in a quantum 2x 2x+1 Q2x Q2x 2x+1 2x + − circuit. At each site x, a spin particle is emulated by two · YQ2x (θ) · CNOTQ2x+1,Q2x · ZQ2x · CNOTQ2x,Q2x+1 , qubits (Q2x−1, Q2x), representing spin up and spin down, respectively. By judiciously designing the gate sequences, ± ∓ θ ± ∓ 3π where Y (θ) = e i 2 σy and Z = e i 4 σz . The corre- we can digitally simulate the on-site spin rotation and sponding quantum circuit is shown in Fig. 1(b). One can two-site SOC, which are the basic buliding blocks of our easily demonstrate that this circuit acts only on the single approach. Finally, by means of Floquet sequences in- excitation qubit states {|gei, |egi}, leaving the other two cluding these two circuit-based building blocks, we can states {|ggi, |eei} unchanged. Similarly, the time evolu- simulate all forms of SOC [67], which play essential roles tion governed by the on-site Hamiltonian H˜ (x) can be in driving the phases to be topologically nontrivial. o programmed through the same form of composite quan- As shown in Fig. 1, we map a spinful lattice into a tum gate U acting on two qubits (Q ,Q ). qubit array. The ground and excited states of each qubit 2x−1 2x are labelled as |gi and |ei respectively. When all qubits To digitally simulate the topological quantum matter are prepared in the ground state, the state of the sys- described by an effective Hamiltonian Heff consisting ˜ ˜ tem is corresponding to the vacuum state of the simu- of Ho(x) and Hs(x), one can construct a quantum cir- lated lattice |vaci = |gg ··· ggi. Flipping one of the cuit based on the aforementioned programmable pairwise qubits into the excited state |ei simulates the creation SOC and on-site spin rotations to mimick the evolution of a particle in the lattice. Specifically, exciting the odd operation (even) qubit Q2x−1 (Q2x) creates a spin up (down) par- −iHeff T ticle at the lattice site x. The associated raising oper- UT = e , (3) ators are mapped into hard-core bosonic creation oper- † † ators through ax,↑ = |eiQ2x−1 hg| and ax,↓ = |eiQ2x hg|. where T is a time period. We will illustrate this approach Therefore, the coupling between the qubits Q2x−1 and with typical examples in 1D and 2D topological matter. Q2x simulates the on-site spin rotation The generalization of the method to higher dimensions is † straightforward. Hõ(x)= Joa ax,↓ + H.c., (1) x,↑ Quantum circuits for 1D topological matter. We il- while the coupling between the qubits Q2x and Q2x+1 lustrate the approach by simulating a 1D topological simulates the two-site SOC Floquet phase described by the time evolution operator U ≡ e−iH1T with the effective Hamiltonian H given by ˜ † 1 1 Hs(x)= Jsax,↓ax+1,↑ + H.c.. (2) Eq. (S6) in Supplemental Material (SM) [68]. We find that the U can be decomposed as Here we will show that the dynanmics of the SOC and on- 1 site spin rotations can be exactly programmed on a quan- −iH1T −iHoT/3 −iHsT/3 −iHoT/3 tum circuit consisting of sequences of universal quantum U1 ≡ e = e e e , (4) gates, including single-qubit rotation gates and two-qubit ˜ entangled gates. We now address our protocol to pro- where Ho = x Ho(x) is the on-site spin rotation and ˜ ˜ gramme the time evolution of the two-site SOC Hs with Hs = x Hs(Px) is the 1D SOC. Based on Eq.(4), we can two- and single-qubit quantum gates acting on the qubits constructP a quantum circuit with three steps of parallel (Q2x,Q2x+1). After an evolution time t, the evolution two-qubit composite quantum gates U to simulate the 3

n measuring mean displacements [23, 72], we find that the U(a ) U(a ) winding numbers (ν ,ν ) can be directly detected via ݐ݁݌2 ͳ ܵݐ݁݌ ʹ ܵݐ݁݌2 ͵ 0 πܵ ଵ U(b ) n ′ n ܳ two quantum circuits (U1) and (U1) [68], where the ܽ ଶ U(a ) U(a ) ܳ 2 2 ′ ଷ U(b ) starting point in U1 is the composite gate U(β) instead ܳ α ܳସ of U( 2 ). Suppose the input state for both circuits is U(b ) |ψini = |gg...e...ggi. We measure the mean displace- U(a ) U(a ) ′ ′ N ′ 2 2 ¯ ment P = hψout| x=1 x|eiQ2x he|)|ψouti at the output ଶேିଷ U(b ) ′ ′ n ܳ a a ଶேିଶ U( ) U( ) of the circuit, whereP |ψouti = (U1) |ψini. As shown in ܳ 2 2 ܳଶேିଵ Fig. 2(c), the values of the winding numbers (ν0,νπ) are ܳଶே equal to the oscillation center of the mean displacements (0,1) (0,1) (P¯ = −P¯ − P¯′, P¯ = P¯ − P¯′), which can be derived ܾ ߥ଴ǡ ߥగ ܿ 0 π (1,1) (0,0) (1,1) from the even qubit excitation distribution at the output ʹ ʹ ଴ ܲത of the circuit [68]. గ / (1,0) (1,0) ܲത Experimental implementation on quantum processors. a p Free access IBM [73] and Rigetti [74] quantum processors (0,0) ͳ Ͳ (1,1) (1,1) have recently been successfully applied to study many- (0,1) (0,1) body quantum states [75–78]. Here we perform a topological quantum circuit on both IBM and Rigetti quan- Ͳ െʹ tum processors. Specifically, an eight-qubit quantum cir- െʹ Ͳ ʹ െʹ Ͳ ʹ cuit (U )3 is programmed on ibm 16 melbourne [68]. In ߚΤ ߨ ߚΤ ߨ 1 FIG. 2: (a) Quantum circuits for 1D topological insulator. our quantum circuits, topological edge states, the hall- The basic unit gates consist of three steps of composite two- mark of topological phases, appear at both E = 0 and qubit gates U, realizing DQS of an effective Hamiltonian H1. n is the cycle number of the unit gates. (b) The phase dia- π, with the number of modes determined by ν0 and νπ gram of the quantum circuits varies with rotation angles (α, respectively [68]. β) in the single-qubit quantum gates. (c) By measuring the In SM [68], we analytically derive the wave functions centers of the mean displacements (P¯0, P¯π), one can detect of the edge states associated with the Floquet topological the winding numbers (ν0, νπ) as a function of β for α = 1.3π phases. In particular, the wave functions of the topologi- (the dashed line shown in (b)). cal edge states localized at the left boundary with energy E =0 or π are derived as

N N evolution U1, i.e, 0 x−1 π ˜x−1 |ψLi = λ1 |eiQ2x−1 , |ψLi = λ2 |eiQ2x , N N−1 α xX=1 xX=1 U1 = U( )Q2x−1,Q2x · U(β)Q2x,Q2x+1 2 α β ˜ α β Ox=1 Ox=1 where λ1 = − tan( 4 ) cot( 4 ) and λ2 = cot( 4 ) cot( 4 ). In- N 0(π) α terestingly, the probability of qubit excitations in |ψL i · U( ) − , (5) 2 Q2x 1,Q2x is maximized in the leftmost odd (even) qubit Q1 (Q2). Ox=1 This feature allows us to experimentally observe and dis- where U(η) denotes the quantum gate U(η) act- criminate the 0 and π energy edge states. Qm,Qn 0(π) ing on the qubits Qm and Qn with the parameters To observe the edge states |ψL i, we prepare 0(π) Jo =3α/4T and Js =3β/2T , as illustrated in Fig. 2(a). the initial input state of the circuit at |ψin i = Here the ﬁrst (third) and second steps simulate the time |egggggggi(|geggggggi), which yields the output state evolution of pairwise on-site spin rotations and pairwise 0(π) 3 0(π) |ψout i = (U1) |ψin i. The measured qubit excitation two-site SOCs, respectively. distributions at the output acquired from the IBM quan- The topology of the U1-based quantum circuit is char- tum processor are shown in Fig. 3. When the quantum acterized by the topological winding number [69–71] circuit is tuned into the topological phase with ν0 = 1 (ν = 1) and supports one 0 (π) energy edge state, as i π π ν = dk tr τ U −1∂ U , (6) shown in Figs. 3(e) and 3(g) (Figs. 3(d) and 3(h)), the ǫ Z x z 1ǫ kx 1ǫ 4π −π qubit excitation at the output maximally populates the where U (ǫ = 0, π) is the periodized evolution oper- leftmost odd (even) qubit Q1 (Q2). The reason for this is 1ǫ 0(π) ator deﬁned by U1 in the momentum space, see SM that the initial input state |ψin i has a large overlap with 0(π) [68]. The values of ν0 and νπ are numerically calcu- |ψL i, while the qubit excitation evolves mainly in the lated and the results are plotted in Fig. 2(b). By circuit based on the 0 (π) energy edge state wave packet changing the rotation angles (α, β), we can tune the and always maximally localizes in the leftmost odd (even) quantum circuit into various topological phases charac- qubit Q1 (Q2). In contrast, as shown in Figs. 3(a) and terized by ν0 = 0, 1 and νπ = 0, 1. Interestingly, by 3(c) (Figs. 3(b) and 3(f)), when the quantum circuits are 4

漓漓 (0,0) (0,0) (0,0) 濅澽濇澿濝濨 (0,1) ଴ గ (0,1) ଴ గ ଴ గ ܹ ǡ ܹ ߥ ൌ Ͳ ߥ ൌ Ͳ ߥ ൌ Ͳ ߥ ൌ ͳ澹濬濤濙濦濝濡濙濢濨 ܽ ܿ ʹ (1,1) (0,0) (1,1) ͳ QU (ݕ x, y (1,0) (1,0 ݀ ܿ ܾ ܽ x, y ߨ

D (0,0) (0,0) (0,0) ௘ Qx, y Ȁ ߜ ܲ ͲǤͷ (1,0) (1,0) Ͳ (1,1) (0,0) (1,1) (0,1) (0,1) 漓漓 (0,0) (0,0) (0,0) Ͳ ଴ గ ଴ గ ߥ ൌ ͳ ߥ ൌ Ͳ ߥ ൌ ͳ ߥ ൌ ͳ n ͳ െʹ ʹ U ݔ െʹ Ͳ ݂ ݁ ݄ Qx, y p ܾ U ( ) D 2 ݀ ߛΤ ߨ Qx, y g ௘ U ( )

ܲ U 2 Qx+1, y p ͲǤͷ U ( ) U (d ) D 2 Qx+1, y U Qx, y + 1 p U ( ) U (d ) QD 2 x, y +1 g 濅濩濖濝濨澔澽濢濘濙濬濅濩濖濝濨澔澽濢濘濙濬濅濩濖濝濨澔澽濢濘濙濬濅濩濖濝濨澔澽濢濘濙濬 U ( ) U 2 Qx+1, y + 1 p Ͳ U ( ) ͳ ʹ ͵ Ͷ ͷ ͸ ͹ ͺ ͳ ʹ ͵ Ͷ ͷ ͸ ͹ ͺ ͳ ʹ ͵ Ͷ ͷ ͸ ͹ ͺ ͳ ʹ ͵ Ͷ ͷ ͸ ͹ ͺ D 2 Qx+1, y +1 FIG. 3: Experimental observation of 0 and π energy topological edge states on the IBM quantum processor. The qubit excitation probability distributions are measured at the out- FIG. 4: (a) Simulating a 2D spinful lattice with a quan- put of an eight-qubit topological quantum circuit when the U 0 π input qubit states are (a,c,e,g) |ψini and (b,d,f,h) |ψini. The tum circuit. Each lattice site is encoded by two qubits Qx,y D single-qubit rotation angles are chosen as (a,b) α = π and and Qx,y. (b) A quantum circuit for 2D topological matter. β = 0.9π; (c,d) α = 1.9π and β = 0.4π; (e,f) α = 0.1π and (c) The Floquet winding numbers (W0,Wπ) as a function of β = 0.7π; (g,h) α = π and β = 1.9π. single-qubit rotation angles (γ, δ). (d) Qubit excitation distribution at the output of a circuit with (γ = 1.9π, δ = 0.8π) U and n = 10 after the qubit Qx=1,y=1 (labelled by the star) is excited. tuned into the trivial phases with ν0 = 0 (νπ = 0) and 0(π) there is no overlap with the edge states |ψL i, the qubit excitations are transferred into the bulk of the circuit by Eq. (S66) in SM [68]. We find that the U2 can be and do not maximally populate the leftmost edge qubits. decomposed as This is because in this case the initial state is a superpo- −iHsy T/3 −iHsxT/3 −iHosT/3 sition of different Bloch bulk states and the qubit excita- U2 = e e e , (7) tions transfer in the circuit via bulk state wave packets. ˜ ˜ Therefore, the 0 and π energy edge states can be distin- where Hos = x,y Ho(x, y) with Ho(x, y) = † P guished by observing whether the leftmost odd qubit Q1 Josax,y,↑ax,y,↓ + H.c. and Jos = 3π/4T , Hsx = † or the leftmost even qubit Q2 is maximally excited. The ˜ ˜ x,y Hsx(x, y) with Hsx(x, y) = Jsxax,y,↓ax+1,y,↑ + experimental results acquired from the Rigetti quantum PH.c and Jsx = 3γ/4T , Hsy = H˜sy(x, y) with processor are shown in SM [68] and demonstrate similar x,y ˜ † P results. In Fig. 3, we show the simulation results de- Hsy(x, y) = Jsyax,y,↓ax,y+1,↑ + H.c. and Jsy = 3δ/2T . rived from noisy quantum circuit QISKit. The mismatch As shown in Fig. 4(b), we further find that such a circuit between the experimental and the simulated results is can be constructed as due to the errors caused by the crosstalk in the quan- N−1 N−1 γ tum gates and the imperfections in state readouts [68]. U2 = U(δ)QD ,QU · U( )QD ,QU x,y x,y+1 2 x,y x+1,y However, due to topological protection, edge state lo- Ox=1 Ox=1 calization can still be unambiguously observed in noisy N π quantum circuits. · U( ) U D , (8) 2 Qx,y ,Qx,y Quantum circuits for 2D topological matter. Our Ox=1

approach is universal and can be applied in high- π where the quantum gate U( ) D U in the first step dimensional systems. In the following example, we show 2 Qx,y ,Qx,y programs the time evolution of on-site spin rotation, and how to design quantum circuits protected by 2D topol- γ U the gates U( 2 )QD ,QU and U(δ)QD ,QU in the sec- ogy. As with the 1D case, we exploit two qubits Qx,y x,y x+1,y x,y x,y+1 and QD to simulate a site of a 2D spinful lattice (see ond and third steps program the time evolution of x- and x,y y-direction SOCs, respectively. U Fig. 4(a)), where the raising operators |eiQx,y hg| and † The topology of the U2-based quantum circuit is |ei D hg| are mapped to the creation operators a and Qx,y x,y,↑ characterized by the Floquet winding number Wǫ=0,π † ax,y,↓ respectively. We aim to simulate a 2D topological [69, 79, 80]. In Fig. 4(c), we show the rich topological Floquet phase described by the time evolution operator phase diagram versus single-qubit rotation angles (γ,δ). −iH2T U2 ≡ e with the effective Hamiltonian H2 given There are four 2D Floquet topological phases featured 5 by W0 = 0, 1 and Wπ = 0, 1. Although the Chern number for the Floquet topological phase (W0 = 1, Wπ = 1) is zero, the system supports topological edge states [68], ∗ Electronic address: [email protected] which can be used as topologically protected quantum † Electronic address: [email protected] channels to transfer quantum state [81–84]. We assume [1] X.-L. Qi and S.-C. Zhang, Topological insulators and su- U that the qubit Qx=1,y=1 is initially excited, as shown in perconductors, Rev. Mod. Phys. 83, 1057 (2011). Fig. 4(b). We use noisy quantum circuit QISKit to calcu- [2] M. Z. Hasan and C. L. Kane, Colloquium: topological late qubit excitation distribution at the output [68] and insulators, Rev. Mod. Phys. 82, 3045 (2010). the results are plotted in Fig. 4(d). It turns out that the [3] R. Prange and S. Girvin, The quantum Hall effect, qubit excitation state can still be transferred along the Sringer, 1990. [4] I. Bloch, J. Dalibard, and W. Zwerger, Many-body physics edges via the edge state quantum channel, which reflects with ultracold gases, Rev. Mod. Phys. 80, 885 (2008). the power of quantum circuits with topological protec- [5] M. Lewenstein, et al., Ultracold atomic gases in optical tion. lattices: mimicking condensed matter physics and beyond, Adv. Phys. 56, 243 (2007). In summary, we have theoretically developed and ex- [6] C. Gross and I. Bloch, Quantum simulaitons with ul- perimentally demonstrated an approach to implement tracold atoms in optical lattices, Science 357, 995-1001 DQS of SOCs and the associated topological phases of (2017). matter on current generation of noisy intermediate-scale [7] R. Blatt and C. F. Roos, Quantum simulations with quantum processors. Our approach is generic and can trapped ions, Nat. Phys. 8, 277-284 (2012). be used to simulate a wide variety of complex SOCs, [8] K. Kim, M.-S. Chang, S. Korenblit, R. Islam, E. E. Ed- wards, J. K. Freericks, G.-D. Lin, L.-M. Duan, and C. such as Rashba [85] and Dresselhaus SOC [86], and re- Monroe, Quantum Simulation of Frustrated Ising Spins lated topological phases that are intractable in atomic, with Trapped Ions, Nature 465, 590-593 (2010). photonic, and solid-state systems. By applying the pro- [9] A. A. Houck, H. E. Türeci and J. Koch, On-chip quantum posed method and utilizing current quantum processors, simulation with superconducting circuits, Nat. Phys. 12, we can, in principle, build and explore the topological 292-299 (2012). matters we are interested in, such as high-dimensional [10] F. Arute et al., Quantum supremacy using a pro- 574 and high-order topological matters that are currently grammable superconducting processor, Nature , 505- 510 (2019). highly pursued [1, 2, 87]. When multiple qubit excita- [11] A. Aspuru-Guzik and P. Walther, Photonic quantum sim- tions are put into topological quantum circuits, inter- ulators, Nat. Phys. 12, 285-291 (2012). acting topological states [88] also can be investigated on [12] N. Goldman, J. C. Budich, and P. Zoller, Topological noisy intermediate-scale quantum processors, which are quantum matter with ultracold gases in optical lattices, beyond the computation capability of current classical Nat. Phys. 12, 639 (2016). computers and expected to show quantum advantages. [13] D.-W. Zhang, Y.-Q. Zhu, Y. X. Zhao, H. Yan, and S.-L. In addition, it would be interesting to investigate topo- Zhu, Topological quantum matter with cold atoms, Adv. Phys. 67, 253 (2019). logical quantum information processing tasks via topo- [14] N. R. Cooper, J. Dalibard, and I. B. Spielman, Topologi- logical quantum circuits [89–91]. cal bands for ultracold atoms, Rev. Mod. Phys. 91, 015005 (2019). Acknowledgements. F. M. thanks Heng Shen for [15] E. J. Meier, F. Alex An, A. Dauphin, M. Maffei, P. invaluable help in improving the presentation of our Massignan, T. L. Hughes, B. Gadway, Observation of manuscript and thanks Ying Hu and Dong-Ling Deng the topological Anderson insulator in disordered atomic for helpful discussions. This work was supported by wires, Science 362, 929 (2018). the National Key Research and Development Program [16] H. Cai, J. Liu, J. Wu, Y. He, S.Y. Zhu, J.X. of China (2017YFA0304203), National Natural Science Zhang, and Da-Wei Wang, Experimental Observation Foundation of China (NSFC) (11604392), Changjiang of Momentum-Space Chiral Edge Currents in Room- Temperature Atoms, Phys. Rev. Lett. 122, 023601 (2019). Scholars and Innovative Research Team in University of [17] D. Xie, T.S. Deng, T. Xiao, W. Gou, W. Yi, and B. Yan, Ministry of Education of China (PCSIRT)(IRT 17R70), Topological quantum walks in momentum space with a Fund for Shanxi 1331 Project Key Subjects Construc- Bose-Einstein condensate, arXiv:1906.12016. tion, and 111 Project (D18001). Y.F. Yu was sup- [18] T. Kitagawa, M. A. Broome, A. Fedrizzi, M. S. Rud- ported by National Natural Science Foundation of China ner, E. Berg, I. Kassal, A. Aspuru-Guzik, E. Demler, A. (Nos.61775062). S. L. Z. was supported by the Key- G. White, Observation of topologically protected bound 3 Area Research and Development Program of Guang- states in photonic quantum walks, Nat. Commun. , 882 (2012). Dong Province (Grant No. 2019B030330001), the Na- [19] J. S. Xu, K. Sun, Y.J. Han, C.F. Li, J. K. Pachos, G. C. tional Key Research and Development Program of China Guo, Simulating the exchange of Majorana zero modes (Grant No. 2016YFA0301800), the National Natural Sci- with a photonic system, Nat. Commun. 7, 13194 (2016). ence Foundation of China (Grant No. 91636218), and the [20] L. Xiao, X. Zhan, Z. H. Bian, K. K. Wang, X. Zhang, Key Project of Science and Technology of Guangzhou X. P. Wang, J. Li, K. Mochizuki, D. Kim, N. Kawakami, (Grant No. 201804020055). W. Yi, H. Obuse, B. C. Sanders, P. Xue, Observation of 6

topological edge states in parity-time-symmetric quantum Polkovnikov, and K.W. Lehnert, Measuring a Topolog- walks, Nat. Phys. 13, 1117 (2017). ical Transition in an Artificial Spin-1/2 System, Phys. [21] F. Cardano, M. Maffei, F. Massa, B. Piccirillo, C. d. Rev. Lett. 113, 050402 (2014). Lisio, G. D. Filippis, V. Cataudella, E. Santamato, L. [34] P. Roushan, C. Neill, Y. Chen, M. Kolodrubetz et al., Marrucci, Statistical moments of quantum-walk dynamics Observation of topological transitions in interacting quan- reveal topological quantum transitions, Nat. Commun. 7, tum circuits, Nature (London) 515, 241 (2014). 11439 (2016). [35] Y. P. Zhong, D. Xu, P. Wang, C. Song, Q.J. Guo, W.X. [22] X. Zhan, L. Xiao, Z. Bian, K. Wang, X. Qiu, B. C. Liu, K. Xu, B.X. Xia, C.-Y. Lu, Siyuan Han, Jian-Wei Sanders, W. Yi, and P. Xue, Detecting Topological In- Pan, and H. Wang, Emulating Anyonic Fractional Sta- variants in Nonunitary Discrete-Time Quantum Walks, tistical Behavior in a Superconducting Quantum Circuit, Phys. Rev. Lett. 119, 130501 (2017). Phys. Rev. Lett. 117, 110501 (2016). [23] F. Cardano, A. D’Errico, A. Dauphin, M. Maffei, B. Pic- [36] T. Wang, Z. Zhang, L. Xiang, Z. Gong, J. Wu, and cirillo, C. d. Lisio, G. D. Filippis, V. Cataudella, E. San- Y. Yin, Simulating a topological transition in a super- tamato, L. Marrucci, M. Lewenstein, P. Massignan, De- conducting phase qubit by fast adiabatic trajectories, Sci. tection of Zak phases and topological invariants in a chi- China Phys. Mech. Astron. 61, 047411 (2018). ral quantum walk of twisted photons, Nat. Commun. 8, [37] X. Tan, D.-W. Zhang, Q. Liu, G. Xue, H.-F. Yu, Y.-Q. 15516 (2017). Zhu, H. Yan, S.-L. Zhu, and Y. Yu, Topological Maxwell [24] Y. Wang, Y.H. Lu, F. Mei, J. Gao, Z.M. Li, H. Tang, S.L. metal bands in a superconducting qutrit, Phys. Rev. Lett. Zhu, S. Jia, X.M. Jin, Direct Observation of Topology 120, 130503 (2018). from Single-photon Dynamics on a Photonic Chip, Phys. [38] C. Song, D. Xu, P. Zhang, J. Wang, Q. Guo, W. Liu, Rev. Lett. 122, 193903 (2019). K. Xu, H. Deng, K. Huang, D. Zheng, Shi-Biao Zheng, [25] J.S. Xu, K. Sun, J. K. Pachos, Y.J. Han, C. F. Li, and H. Wang, X.B. Zhu, C.Y. Lu, and J.W Pan, Demonstra- G. C. Guo, Photonic implementation of Majorana-based tion of topological robustness of anyonic braiding statis- Berry phases, Sci. Adv. 4, 6533 (2018). tics with a superconducting quantum circuit, Phys. Rev. [26] X.Y. Xu, Q. Wang, W. Pan, K. Sun, J.S. Xu, G. Chen, Lett. 121, 030502 (2018). J.S. Tang, M. Gong, Y.J. Han, C.F. Li, G.C. Guo, Mea- [39] X. Tan, Y.X. Zhao, Q. Liu, G. Xue, H.F. Yu, Z. D. Wang, suring the winding number in a large-scale chiral quan- and Y. Yu, Simulation and manipulation of tunable Weyl- tum walk, Phys. Rev. Lett. 120, 260501 (2018). semimetal bands using superconducting quantum circuits, [27] C. Chen, X. Ding, J. Qin, Y. He, Y. Luo, M. Chen, Phys. Rev. Lett. 122, 010501 (2019). C. Liu, X. Wang, W. Zhang, H. Li, L. You, Z. Wang, [40] S. Lloyd, Universal quantum simulators, Science 273, D. Wang, B. C. Sanders, C.Y. Lu, and J.W. Pan, Ob- 1073-1078 (1996). servation of topologically protected edge states in a pho- [41] I. Buluta, F. Nori, Quantum Simulators, Science 326, tonic two-dimensional quantum walk, Phys. Rev. Lett. 108-111 (2009). 121, 100502 (2018). [42] B. P. Lanyon, C. Hempel, D. Nigg, M. Müller, R. Ger- [28] B. Wang, T. Chen, and X. Zhang, Experimental observa- ritsma, F. Zähringer, P. Schindler, J. T. Barreiro, M. tion of topologically protected bound states with vanishing Rambach, G. Kirchmair, M. Hennrich, P. Zoller, R. Blatt, chern numbers in a two-dimensional quantum walk, Phys. C. F. Roos, Universal digital quantum simulation with Rev. Lett. 121, 100501 (2018). trapped ions, Science 334, 57 (2011). [29] W. Cai, J. Han, Feng Mei, Y. Xu, Y. Ma, X. Li, H. Wang, [43] M. Müller, S. Diehla, G. Pupillo, P. Zoller, Engineered Y. Song, Z.Y. Xue, Z.Q. Yin, S. Jia, Luyan Sun, Obser- Open Systems and Quantum Simulations with Atoms vation of topological magnon insulator states in a super- and Ions, Advances In Atomic, Molecular, and Optical conducting circuit, Phys. Rev. Lett. 123, 080501 (2019). Physics 61, 1 (2012). [30] A. D. King, J. Carrasquilla, I. Ozfidan, J. Raymond, E. [44] I. M. Georgescu, S. Ashhab, F. Nori, Quantum simula- Andriyash, A. Berkley, M. Reis, T. M. Lanting, R. Har- tion, Rev. Mod. Phys. 86, 153 (2014). ris, G. Poulin-Lamarre, A. Yu. Smirnov, C. Rich, F. Al- [45] L. Lamata, A. Parra-Rodriguez, M. Sanz, E. Solano, tomare, P. Bunyk, J. Whittaker, L. Swenson, E. Hoskin- Digital-analog quantum simulations with superconducting son, Y. Sato, M. Volkmann, E. Ladizinsky, M. Johnson, circuits, Advances in Physics: X 3, 1457981 (2018). J. Hilton, M. H. Amin, Observation of topological phe- [46] S. McArdle, S. Endo, A. Aspuru-Guzik, S. Ben- nomena in a programmable lattice of 1,800 qubits, Nature jamin, X. Yuan, Quantum computational chemistry, (London) 560, 456 (2018). arXiv:1808.10402. [31] E. Flurin, V. V. Ramasesh, S. Hacohen-Gourgy, L. S. [47] Y. Cao, J. Romero, J. P. Olson, M. Degroote, P. D. Martin, N. Y. Yao, and I. Siddiqi, Observing Topologi- Johnson, M. Kieferová, I. D. Kivlichan, T. Menke, B. cal Invariants Using Quantum Walks in Superconducting Peropadre, N. P. D. Sawaya, S. Sim, L. Veis, A. Aspuru- Circuits, Phys. Rev. X 7, 031023 (2017). Guzik, Quantum Chemistry in the Age of Quantum Com- [32] P. Roushan, C. Neill, A. Megrant, Y. Chen, R. Bab- puting, arXiv:1812.09976. bush, R. Barends, B. Campbell, Z. Chen, B. Chiaro, A. [48] H. Weimer, M. Müller, I. Lesanovsky, P. Zoller and H. Dunsworth, A. Fowler, E. Jeffrey, J. Kelly, E. Lucero, Peter Büchler, A Rydberg quantum simulator, Nat. Phys. J. Mutus, P. J.J. O’Malley, M. Neeley, C. Quintana, D. 6, 382 (2010). Sank, A. Vainsencher, J. Wenner, T. White, E. Kapit, H. [49] A. Mezzacapo, J. Casanova, L. Lamata, E. Solano, Digi- Neven, J. Martinis, Chiral ground-state currents of inter- tal Quantum Simulation of the Holstein Model in Trapped acting photons in a synthetic magnetic field, Nat. Phys. Ions, Phys. Rev. Lett. 109, 200501 (2012). 13, 146 (2017). [50] J. Casanova, A. Mezzacapo, L. Lamata, E. Solano, Quan- [33] M. D. Schroer, M. H. Kolodrubetz, W. F. Kindel, M. tum Simulation of Interacting Fermion Lattice Models in Sandberg, J. Gao, M. R. Vissers, D. P. Pappas, A. Trapped Ions, Phys. Rev. Lett. 108, 190502 (2012). 7

[51] U. L. Heras, A. Mezzacapo, L. Lamata, S. Filipp, A. tum phase transition in topological insulators, Phys. Rev. Wallraff, E. Solano, Digital Quantum Simulation of Spin B 89, 045432 (2014). Systems in Superconducting Circuits, Phys. Rev. Lett. [66] Q. Cheng, Y. Pan, H. Wang, C. Zhang, D. Yu, A. Gover, 112, 200501 (2014). H. Zhang, T. Li, L. Zhou, S. Zhu, Observation of anoma- [52] R. Barends, L. Lamata, J. Kelly, L. Garc´ıa-Alvarez,´ A. lous π modes in photonic Floquet engineering, Phys. Rev. G. Fowler, A. Megrant, E. Jeffrey, T. C. White, D. Sank, Lett. 122, 173901 (2019). J. Y. Mutus, B. Campbell, Yu Chen, Z. Chen, B. Chiaro, [67] V. Galitski and I. B. Spielman, Spin-orbit coupling in A. Dunsworth, I.-C. Hoi, C. Neill, P. J. J. O’Malley, quantum gases, Nature 494, 49-54 (2013). C. Quintana, P. Roushan, A. Vainsencher, J. Wenner, [68] See Supplemental Material. E. Solano, J. M. Martinis, Digital quantum simulation [69] M. S. Rudner, N. H. Lindner, E. Berg, and M. Levin, of fermionic models with a superconducting circuit, Nat. Anomalous edge states and the bulk-edge correspondence Commun. 6, 7654 (2015). for periodically driven two-dimensional systems, Phys. [53] Y. Salathé, M. Mondal, M. Oppliger, J. Heinsoo, P. Rev. X 3, 031005 (2013). Kurpiers, A. Potoˇcnik, A. Mezzacapo, U. Las Heras, L. [70] M. Fruchart, Complex classes of periodically driven topo- Lamata, E. Solano, S. Filipp, A. Wallraff, Digital quan- logical lattice systems, Phys. Rev. B 93, 115429 (2016) tum simulation of spin models with circuit quantum elec- [71] J. K. Asboth, B. Tarasinski, P. Delplace, Chiral sym- trodynamics, Phys. Rev. X 5, 021027 (2015). metry and bulk–boundary correspondence in periodi- [54] R. Barends, A. Shabani, L. Lamata, J. Kelly, A. Mezza- cally driven one-dimensional systems, Phys. Rev. B 90, capo, U. Las Heras, R. Babbush, A. G. Fowler, B. Camp- 125143 (2014). bell, Yu Chen, Z. Chen, B. Chiaro, A. Dunsworth, E. [72] S. Longhi, Probing one-dimensional topological phases Jeffrey, E. Lucero, A. Megrant, J. Y. Mutus, M. Neeley, in waveguide lattices with broken chiral symmetry, Opt. C. Neill, P. J. J. O’Malley, C. Quintana, P. Roushan, D. Lett. 43, 4639 (2018). Sank, A. Vainsencher, J. Wenner, T. C. White, E. Solano, [73] https://raw.githubusercontent.com/Qiskit/qiskit/master H. Neven, J. M. Martinis, Digitized adiabatic quantum /Qiskit.bib. computing with a superconducting circuit, Nature (Lon- [74] R. S. Smith, M. J. Curtis, W. J. Zeng, A Practical Quan- don) 534, 222 (2016). tum Instruction Set Architecture, arXiv:1608.03355. [55] A. Smith, M. S. Kim, F. Pollmann, and J. Knolle, Simu- [75] A. Chiesa, F. Tacchino, M. Grossi, P. Santini, I. Taver- lating quantum many-body dynamics on a current digital nelli, D. Gerace, S. Carretta, Quantum hardware simu- quantum computer, arXiv:1906.06343. lating four-dimensional inelastic neutron scattering, Nat. [56] S. P. Jordan, K. S. M. Lee, J. Preskill, Quantum al- Phys. 15, 455 (2019). gorithms for quantum field theories, Science 336, 1130 [76] M. Motta, C. Sun, A. T. K. Tan, M. J. O’Rourke, E. Ye, (2012). A. J. Minnich, F. G. S. L. Brandão, G. Kin-Lic Chan, [57] L. Garc´ıa-Alvarez,´ J. Casanova, A. Mezzacapo, I. L. Determining eigenstates and thermal states on a quantum Egusquiza, L. Lamata, G. Romero, E. Solano, Fermion- computer using quantum imaginary time evolution, Nat. Fermion Scattering in Quantum Field Theory with Super- Phys. (2019) doi:10.1038/s41567-019-0704-4. conducting Circuits, Phys. Rev. Lett. 114, 070502 (2015). [77] K. Choo, C. W. von Keyserlingk, N. Regnault, and T. [58] A. Mezzacapo, E. Rico, C. Sab´ın, I. L. Egusquiza, L. Neupert, Measurement of the Entanglement Spectrum of Lamata, E. Solano, Non-Abelian SU(2) Lattice Gauge a Symmetry-Protected Topological State Using the IBM 121 Theories in Superconducting Circuits, Phys. Rev. Lett. Quantum Computer, Phys. Rev. , 086808 (2018). 115, 240502 (2015). [78] A. Smith, B. Jobst, A. G. Green, and F. Pollmann, Cross- [59] E. Zohar, A. Farace, B. Reznik, J. I. Cirac, Digital ing a topological phase transition with a quantum com- quantum simulation of Z2 lattice gauge theories with dy- puter, arXiv:1910.05351v1. namical fermionic matter, Phys. Rev. Lett. 118, 070501 [79] T. Kitagawa, E. Berg, M. Rudner, E. Demler, Topological (2017). characterization of periodically driven quantum systems, 82 [60] L. Garc´ıa-Alvarez,´ I. L. Egusquiza, L. Lamata, A. del Phys. Rev. B , 235114 (2010). [80] M. D. Reichl and E. J. Mueller, Floquet edge states with Campo, J. Sonner, E. Solano, Digital quantum simula- 89 tion of minimal AdS/CFT, Phys. Rev. Lett. 119, 040501 ultracold atoms, Phys. Rev. A , 063628 (2014). (2017). [81] N. Y. Yao, C. R. Laumann, A. V. Gorshkov, H.Weimer, [61] E. A. Martinez, C. A. Muschik, P. Schindler, D. Nigg, A. L. Jiang, J. I. Cirac, P. Zoller, and M. D. Lukin, Topo- Erhard, M. Heyl, P. Hauke, M. Dalmonte, T. Monz, P. logically protected quantum state transfer in a chiral spin Zoller, R. Blatt, Real-time dynamics of lattice gauge the- liquid, Nat. Commun. 4, 1585 (2013). ories with a few-qubit quantum computer, Nature (Lon- [82] C. Dlaska, B. Vermersch, and P. Zoller, Robust quantum don) 534, 516 (2016). state transfer via topologically protected edge channels in [62] C. Kokail, C. Maier, R. van. Bijnen, T. Brydges, M. K. dipolar arrays, Quantum Sci. Technol. 2, 015001 (2017). Joshi, P. Jurcevic, C. A. Muschik, P. Silvi, R. Blatt, [83] Feng Mei, G. Chen, L. Tian, S.L. Zhu, and Suotang Jia, C. F. Roos, P. Zoller, Self-verifying variational quantum Robust quantum state transfer via topological edge states simulation of lattice models, Nature (London) 569, 355 in superconducting qubit chains, Phys. Rev. A 98, 012331 (2019). (2018). [63] B. Murta, G. Catarina, and J. F. Rossier, Phys. Rev. A [84] S. Tan, R. Weda Bomantara, J. Gong, High-fidelity and 101, 020302(R) (2020). long-distance entangled-state transfer with Floquet topo- [64] S. Tamiya and Y.O. Nakagawa, arXiv:2003.01706v1. logical edge modes, arXiv:1909.03646. [65] C. Ju, C. Lei, X. Xu, D. Culcer, Z. Zhang, and J. Du, [85] Y. A. Bychkov and E. I. Rashba, Oscillatory effects and NV-center-based digital quantum simulation of a quan- the magnetic susceptibility of carriers in inversion layers. J. Phys. Chem. 17, 6039C6045 (1984). 8

[86] G. Dresselhaus, Spin-orbit coupling eﬀects in zinc blende information processing in 1D wire networks, Nat. Phys. structures. Phys. Rev. 100, 580C586 (1955). 7, 412 (2011). [87] W. A Benalcazar, B. A.Bernevig, and T. L. Hughes, [90] C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. Quantized electric multipole insulators, Science 357, 61, D. Sarma, Non-Abelian anyons and topological quantum (2017). computation Rev. Mod. Phys. 80, 1083 (2008). [88] X. Chen, Z.-C. Gu, Z.-X. Liu, and X.-G. Wen, Symme- [91] A Stern, N.H. Lindner, Topological quantum try protected topological orders in interacting bosonic sys- computation-from basic concepts to ﬁrst experiments, tems, Science 338, 1604 (2012). Science 339, 1179 (2013). [89] J. Alicea, Y. Oreg, G. Refael, F. von Oppen, M. P. A. Fisher, Non-Abelian statistics and topological quantum