Quick viewing(Text Mode)

Arxiv:1602.03217V2 [Quant-Ph] 18 Jul 2017 Egxzcan Efidta Ti Osbet En a Define to Possible Is It That find We Heisen- Quantum Chain

Arxiv:1602.03217V2 [Quant-Ph] 18 Jul 2017 Egxzcan Efidta Ti Osbet En a Define to Possible Is It That find We Heisen- Quantum Chain

arXiv:1602.03217v2 [quant-ph] 18 Jul 2017 egXZcan efidta ti osbet en a define to possible is it that find We Heisen- quantum chain. a XXZ in berg bound-states topological features topological interactions. between inter-particle interplay and the invariant clarify topological a to define counter- to and easy interacting not is weakly it correlations, the of quite parts that is states from topological different interacting strongly of ization ist ohtheory both difficul- general, to a much ties In pose in states topological emerge interacting system. strongly can multi-magnon clar- bound-states interacting to strongly topological importance great whether of single- ify is it of excitations, and feature topological theory both topological in studied the experiment been has years, excitations few magnon past the ei esnegcanapa sqaiatceexcita- as ferromag- so-called a appear the of tions, chain state Heisenberg collec- ground well-known the netic typ- The around the excitations of spins. sys- one tive interacting is of magnetic chain models quantum Heisenberg ical quantum of The properties tems. excitations low-energy collective into and insights fundamental provide oia hsshsbe salse completely established been topo- has noninteracting phases of logical understanding theoretical now, esnegcan opiigbsncaosi optical in quantum the bosonic lattices in comprising observed chains been Heisenberg interaction have bound-states the their when enough appear strong states, is do scattering bound-states continuous to magnon addition in interaction, a nuaos h tde ftplgclmteshave years matters recent topological in of interests extensive studies attracted the insulators, cal oooia annbudsae nproial modulate periodically in bound-states magnon Topological nti ril,w useanwdrcint xlr the explore to direction new a pursue we article, this In h td fmgosi unu pnmdl can models quantum in magnons of study The nteohrhn,snetedsoeyo topologi- of discovery the since hand, other the On 20,21 4,5 npriua,det h toginter-particle strong the to due particular, In . . 11–17 nihsit h hrceiainadudrtnigo st of topological understanding pursue and to characterization prospect topol the new of a into existence opens insights only the not show study and Our distrib bound-states density two-magnon t spectra, characterize the energy definitio to calculate the us We to enables excitations. key co which the under states, is quantum chain particle symmetry Heisenberg cotranslational quantum inter the strongly the that a in explore we bound-state Here, magnon lattices. optical Heise in quantum atoms in observed eleme experimentally of been states have bound magnets Recently, experiment. and theory both togyitrcigtplgclsae nmulti-particl in states topological interacting Strongly eodteseai o single-magnon for scenario the Beyond . .INTRODUCTION I. 1–3 ihuQin Xizhou 18 eety annecttosand excitations magnon Recently, . n experiment and u a-e nvriy(unzo aps,Gaghu5102 Guangzhou Campus), (Guangzhou University Yat-Sen Sun u a-e nvriy(hhiCmu) hhi598,China 519082, Zhuhai Campus), (Zhuhai University Yat-Sen Sun 1 1 2 ini eerhCne colo hsc n Astronomy, and of School & Center Research TianQin u oteferromagnetic the to Due . e aoaoyo polcrncMtrasadTechnologi and Materials Optoelectronic of Laboratory Key 1 egMei Feng , 19 h character- The . 1 ogunKe Yongguan , 6–8 pto Up . Dtd uy1,2017) 19, July (Dated: 9,10 In . 1 , 2 pncan opiigutaodatoms quantum ultracold possible via comprising is predictions chains spin our it examine Moreover, experimentally de- to theory. many-body using perturbation we by generate bound-states, Hamiltonian effective magnon an to these derive insights of give origin To topological the bound-states. magnon the topological explore of and mod- chain we XXZ periodically a Heisenberg invariant, in ulated topological excitations two-magnon our the cotrans- With study the upon symmetry. based defined lational directly is invariant ical enn oooia nain o neatn ytm is systems interacting method Green-function for the invariant topological defining clivratcnb sdt hrceietetopological bound-states. the characterize multi-magnon to of used feature be can topolog- invariant our ical that demonstrate We system. multi-particle a presents symmetry system tional the if invariant topological KNivratfritrcigssesi yuigthe using condition by to boundary is twisted systems only so-called interacting of applies strongly for generalization thus invariant zone, well-known and TKNN A Brillouin states, systems. first Bloch noninteracting the the on over depends states Bloch Chern of first the number) (i.e. Thouless-Kohmoto- invariant famous (TKNN) in The Nightingale-Nijs symmetry system. intrinsic considered the the with associates directly nition on-tts nodrt lrf h oooia rgnof origin topological the magnon clarify of to order spectrum In two- butterfly-like bound-states. the reveal give and Chern defined We numbers period- the calculate a spectrum, chain. excitation in magnon Heisenberg excitations two-magnon modulated In study ically we symmetry. this IV, proper- a sym- with Sec. define cotranslational basic associates and the invariant systems its introduce topological multi-particle we discuss interacting II, in III, Sec. and metry Sec. In In model our ties. motivation. describe and we background the troduces iZhang Li , oprdwt te oooia nains u defi- our invariants, topological other with Compared hsatcei raie sflos hsscinin- section This follows. as organized is article This unu ytm oegetcalne to challenges great pose systems quantum e tos orltosadCennmesof numbers Chern and correlations utions, ogyitrcigtplgclstates. topological interacting rongly cigtplgclsaecle topological called state topological acting etplgclfaue fmulti-magnon of features topological he tr pnwvs(ann)i quantum in (magnons) waves spin ntary br hiscmrsn lrcl Bose ultracold comprising chains nberg 22 fatplgclivratfrmulti- for invariant topological a of n annbudsae,btas gives also but bound-states, magnon hc sa nerlo h er curvature Berry the of integral an is which , gclpoetdeg bound-states. edge protected ogical 1 , 2 rnltoa ymty efind We symmetry. translational n hoogLee Chaohong and , olcietasainlivrac fa of invariance translational collective : esnegXZchains XXZ Heisenberg d 5 China 75, es, 24,25 and ieety u topolog- our Differently, . 1 , 2 , 23 ∗ nte a to way Another . 4,5 rions or cotransla- 26–28 . 2 magnon bound-states in our model, we derive an effec- III. TOPOLOGICAL INVARIANT UNDER tive single-particle Hamiltonian by using the degenerate COTRANSLATIONAL SYMMETRY perturbation theory in Sec. V. Finally, we discuss the experimental possibilities and summarize our results in The quantum Heisenberg chain (1) actually respects Sec. VI. a cotranslational symmetry associated with the cotrans- lational operator. Under periodic boundary condition (BC), the center-of-mass (c.o.m) quasi-momentum asso- ciated with the cotranslational operator is a good quan- tum number. Based on the above c.o.m quasi-momentum II. MODEL and the periodic parameter of the modulation field, we successfully construct an artificial torus and then define a Chern topological invariant. We consider a Heisenberg XXZ chain described by the In order to define the cotranslational symmetry, we following Hamiltonian, introduce the N-magnon wavefunction. For a Lt- lattice system under the periodic BC, the N-magnon Hilbert space can be spanned by the basis, B(N) = Hˆ = − J Sˆ+Sˆ− + Sˆ−Sˆ+ + ∆SˆzSˆz ˆ+ ˆ+ ˆ+ 0 0 l l+1 l l+1 l l+1 |l1l2 ...lN i = Sl1 Sl2 ... SlN | i , in which | i = l |↓↓ ... ↓i and 1 ≤ l < l < ... < l ≤ X h  i  1 2 N ˆz ˆz Lt. An arbitrary N-magnon state can be expanded + BlSl + B0 Sl , (1) as, |Ψi = l

4

FIG. 1. (Color online) Two-magnon excitation spectra. (a) Spectrum for periodic BC (Lt = 99). (b) Magnon-bound-state spectrum for periodic BC (Lt = 99). (c) Spectrum for open BC (Lt = 100). (d) Magnon-bound-state spectrum for open BC (Lt = 100). (e) Magnon density distribution nl = hψ| nˆl |ψi corresponding to the eigenstates of the red pentagrams (A - edge state localized at right, B - bulk state, C - edge state localized at left) in the spectrum (d). (f) Two-magnon correlations Γqr corresponding to (e) (here only the region with non-zero correlations is presented). (g) Two-magnon minor diagonal correlations Γl,l+1 of (f). The other parameters are chosen as ∆/J = 100, λ/J = 0.04, β = 1/3 and δ = π/6.

Chern numbers indicate that the magnon bound-states in these three subbands are nontrivial topological states. Meanwhile, one can find that the Chern number distribu- tions of these subbands resemble the ones of the Harper- Hofstadter model35,36, which has been used to describe the integer quantum Hall effect22,37,38. This feature is quite remarkable because our system is of strong inter- particle interactions and the Harper-Hofstadter model is interaction-free. The identical topological invariants of the two models means that they share the same topo- logical origins. This point is further supported by the butterfly-like spectrum for the two-magnon states, see Fig. 2. Interestingly, the butterfly-like spectrum also re- sembles the one for the Harper-Hofstadter model. FIG. 2. The butterfly-like spectrum for the two-magnon In addition to the nontrivial Chern numbers, another bound-states. We diagonalize the Hamiltonian (1) with hallmark of topological phases is the appearance of edge ∆/J = 100, Lt = 100 and different values of δ: (a) δ = 0 and − states in the system under open BC. Due to the open (b) δ = π/4. The values of λ are set as λ/J = J[∆ cos (πβ)] 1 BC, the c.o.m quasi-momentum k is no longer a good for β 6= 1/2, and λ = 0 for β = 1/2. The energies E/J are quantum number, see Fig. 1 (c). Besides the extended added with a constant ∆/J + 2J/∆. bound-states, the edge bound-states appear in the energy gaps, see the pentagrams A and C in Fig. 1 (d). To il- lustrate these bound-states are edge states, we plot their to 2π, the two-magnon spectrum varies continuously and density distribution in Fig. 1 (e). The density distribu- periodically with a period 2π/q. In Fig. 3, we show how tion clearly shows that these bound-states do localize on the two-magnon spectrum changes with δ. According to the left or right. In Fig. 1 (f), we calculate the two- the bulk-edge correspondence, the topological bulk state ˆ+ ˆ+ ˆ− ˆ− magnon spin correlations Γqr = hψ| Sq Sr Sr Sq |ψi for under periodic BC corresponds to the edge state under localized edge and extended bulk bound-states. The left open BC. We do find one edge state traversing in the and right edge bound-states manifest as the correlation energy gap for the system under open BC. peaks on lower left and upper right quarters along the two It is important to clarify the role of interaction in minor diagonal lines, respectively. In contrast, the bulk the emerging of topological magnon bound-states. Un- bound-states extensively distribute along the whole mi- der strong interaction, the energy spectrum will be sep- nor diagonal lines. In Fig. 1 (g), we show the two-magnon arated into two bands: the upper band for continuum correlations along the minor diagonal line Γl,l+1. More- states and the lower band for bound-states, see Fig. 1 over, when the modulation parameter δ changes from 0 (a). Our calculation shows that the topological states ly- 5

hopping term and the modulation term,

ˆ ˆ+ ˆ− ˆ− ˆ+ ˆz H1 = −J Sl Sl+1 + Sl Sl+1 + BlSl , (11) l l X  X as a perturbation to the sum of the interaction term and the uniform field term,

ˆ ˆz ˆz ˆz H0 = −∆ Sl Sl+1 + B0 Sl . (12) Xl Xl Within a many-body perturbation theory up to second order (see Appendix B), the effective Hamiltonian reads as

L L J 2 t t Hˆ = − (ˆb† ˆb + H.c.)+ µ ˆb† ˆb , (13) eff ∆ m m+1 m m m m=1 m=1 X X ′ ′ ′ where µm = λ cos(2πβm + δ ) with the parameters λ = 2λ cos(πβ) and δ′ = δ + πβ. In our derivation, we have chosen the periodic BC and omitted a constant energy shift. In the effective Hamiltonian, the two-magnon bound- states are described as a single-particle states, which ˆ† are associated with the creation operator bm. Obvi- ˆ† 0 ˆ+ ˆ+ ously, bm | i is equivalent to SmSm+1 |↓↓ ... ↓i. This ˆ† FIG. 3. (Color online) The two-magnon spectrum versus the means that, bm simultaneously flips two neighbouring periodic parameter δ. (a-b) periodic BC (Lt = 99); (c-d) spins: one on the m-th lattice site and the other one open BC (Lt = 100). (b) and (d) are the enlarged magnon- on the (m + 1)-th lattice site. It is easy to find that bound-state spectra in (a) and (c), respectively. The other the form of the effective Hamiltonian is same as the parameters are chosen as ∆/J = 100, λ/J = 0.04, and β = Harper-Hofstadter model35,36. Thus, the origin of topo- 1/3. logical magnon bound-states can be understood from the integer quantum Hall states in the Harper-Hofstadter model. That is, the topological magnon bound-state can ing in the lower band are bound-states. This means that be looked as an analogue to integer quantum Hall state. the topological states in our interacting system are indeed Therefore, the corresponding topological features could interaction-assisted. Moreover, the interaction-assisted be explained by the Chern topological invariant and the topological states are not limited to two-magnon exci- edge states (see Appendix B). The topological features tations, but also could be generalized to multi-magnon given by the effective Hamiltonian are well consistent excitations, which provide extra robustness for these with the numerical ones without using the perturbation topological bound-states. In contrast, most topologi- theory, which are presented in the previous section. cal states found in will be Since the effective Hamiltonian is derived in the con- destroyed when inter-particle interaction has been in- text of the second-order perturbation theory, it is unclear troduced. Thus, our study adds a new member for whether the perturbation treatment retains the topolog- interaction-enabled topological states and also sheds light ical nature of the original system. Through introducing on understanding interacting topological states. an auxiliary Hamiltonian, which is a combination of the effective Hamiltonian and the original Hamiltonian, we demonstrate the topological equivalence between the ef- V. EFFECTIVE SINGLE-PARTICLE fective system and the original system (see Appendix C). HAMILTONIAN

To give an effective picture for the topological magnon VI. SUMMARY AND DISCUSSIONS bound-states and reveal the resemblance between our model and the Harper-Hofstadter model35,36, we derive In summary, we have explored the existence of topolog- an effective Hamiltonian by employing the many-body ical magnon bound-states in strongly interacting quan- perturbation theory. Under strong anisotropy (|J/∆| ≪ tum spin chains under cotranslational symmetry. We 1 and |λ/∆| ≪ 1), which corresponds to strong inter- find that the cotranslational symmetry allows us to define actions between magnons, one can treat the sum of the a topological invariant for the interacting multi-magnon 6 system. Through analyzing their topological features, is very similar to the use of c.o.m momentum for calcu- we demonstrate that topological magnon bound-states lating our defined Chern number. In Ref.44, an effective do emerge and they are interaction-enabled. By em- Hamiltonian for two-hole states in strongly interacting ploying the many-body degenerate perturbation theory, systems are given, which connects the two-hole model to we derive an effective model for the topological magnon a noninteracting single-hole model. bound-states which resembles the interaction-free inte- ger quantum Hall model, thus the topological magnon bound-states could be understood as quantum Hall states ACKNOWLEDGMENTS of bounded magnons. Our study not only provides a new prospect to explore topological magnon excitations, but We thank T. Graß for providing us their recent work28 also sheds lights on understanding strongly interacting on the experimental realization of one-dimensional spin topological states. systems with synthetic magnetic fluxes. This work is Based on the current state-of-art in well-controlled supported by the National Basic Research Program of many-body quantum systems, it is possible to realize China (Grant No. 2012CB821305) and the National Nat- our quantum spin chain with ultracold atoms. In re- ural Science Foundation of China (Grant No. 11374375, cent years, ultracold atoms have provided a leading plat- 11574405). form for simulate various many-body quantum models in condensed matter physics and other fields. The Heisen- berg XXZ chain has been experimentally realized with Appendix A: Calculation of the Chern number two-component bosonic atoms in a one-dimensional opti- 4,5 cal lattice , in which both single-magnon and magnon- Below we give the procedure of calculating the Chern bound-states have been experimentally observed. To en- number via our definition (8). Firstly, the 2D parameter 4,5 gineer the desired on-site modulation on the systems , space is discretized as {kα =2πα/Lt,δj =2πj/Nδ} with one can apply an additional state-dependent optical lat- α = {1, 2,...,Lt/q} and j = {1, 2,...,Nδ}, thus the 39 tice , the Hamiltonian reads as Chern number is given as ˆ ˆ† ˆ Uσσ 1 HBH = − tσb bkσ + nˆjσ(ˆnjσ − 1)ˆ C = S (F ), (A1) jσ 2 n 2π α,j n hj,ki,σ=↑,↓ jσ α,j X X X + U↑↓ nˆj↑nˆj↓ + Vjσ nˆjσ, (14) with j jσ X X kα+1 δj+1 where the state-dependent optical lattice potential gives Sα,j (Fn)= dk dδ Fn(k,δ), (A2) kα δj Vj↑ = −Vj↓ = λ cos(2πβj + δ). Using the experimen- Z Z tal techniques in4,5, one can prepare an atomic Mott in- being the integral of the Berry curvature over the rect- sulator which obeys the spin chain model (1) with the angular region, Rα,j = [kα, kα+1] × [δj ,δj+1]. Accord- parameters ing to Stokes’ formula, the integral Sα,j (Fn) is equal to the line integral of the Berry connection, (Ak, Aδ) = 2 2 2 2 2t↑t↓ 4t↑ 4t↓ 2(t↑ + t↓) i hφ | (∂ , ∂ ) |φ i, along the boundary of R (de- J = , ∆= + − , (15) k,δ k δ k,δ α,j U↑↓ U↑↑ U↓↓ U↑↓ noted as Bα,j ), i.e., the Berry along Bα,j. The Berry phase can be given as the limit of the dis- ˆ+ 45 where the pseudospin operators can be defined as Sl = crete geometric phase as Nδ → ∞ and Lt → ∞ . If the ˆ† ˆ ˆ− ˆ† ˆ ˆz integral region Rα,j is sufficiently small, the Berry phase bl↑bl↓, Sl = bl↓bl↑ and Sl = (ˆnl↑ − nˆl↓)/2. The interac- tion strength between magnons can be tuned by Feshbach along Bα,j can be approximately replaced by the discrete resonance40,41. geometric phase along the four vertices of the rectangular 46 It is also possible to realize our quantum spin chain region Rα,j , that is, and observe the topological magnon bound-states with ultracold trapped ions. In last years, the rapidly devel- Sα,j (Fn)= −Im ln U12U23U34U41 . (A3) oping techniques in engineering ultracold trapped ions h i ′ ′ ′  have become sufficiently mature to simulate the mag- Here, Ujj = hφj | φj i / |hφj | φj i|, |φ1i = φkα,δj , |φ i = φ , |φ i = φ , |φ i = φ , netic phases and the quantum dynamics in Heisenberg 2 kα+1,δj 3 kα+1,δj+1 4 kα,δj+1 spin chains26–28. and the principal branch of the natural logarithm is spec-

Note added: During the finalization of the present ver- ified in (−π, π]. sion (see Ref.42 for our previous version), we became aware of the recent studies on two-particle topological states in an interacting bosonic43 or fermionic44 SSH Appendix B: Derivation of the effective Hamiltonian chain. In Ref.43, the c.o.m momentum of the two in- teracting particle is used as a good quantum number for Under strong anisotropy (|J/∆| ≪ 1 and |λ/∆| ≪ 1), calculating the two-particle generalized Zak phase, which the two magnons can form a magnon bound-state and 7

ˆz ′ thus can be recognized as a single quasiparticle. In this Since hEl1l2 | Sl′ |Gm i = 0, the second-order perturba- section, we will use the 2nd order perturbation theory for tion hˆ2 = Pˆ0Hˆ1SˆHˆ1Pˆ0 reads as many-body degenerate quantum systems47 to derive an 2 effective Hamiltonian for the two-magnon bound-states. ˆ J ˆ h2 =− |Gmi hGm| Tl |El1l2 i As |J/∆| ≪ 1 and |λ/∆| ≪ 1, separating the Hamil- ∆ ′ ′ mm ll l1l2 tonian (1) as Hˆ = Hˆ0 + Hˆ1, one can treat X  ˆ ′ ′ ′ × hEl1l2 | Tl |Gm i hGm | , (B8) Hˆ = −J Sˆ+Sˆ− + Sˆ−Sˆ+ + B Sˆz. (B1) 1 l l+1 l l+1 l l ˆ ˆ+ ˆ− ˆ− ˆ+  l l where, Tl = Sl Sl+1 + Sl Sl+1, the summation indices ′ ′ X  X {m,m ,l,l } take values from {1, 2,...,Lt} and {l1,l2} as a perturbation to is summed over all states of El1l2 6= E0. Introducing the notations ˆ ˆz ˆz ˆz H0 = −∆ Sl Sl+1 + B0 Sl . (B2) l l T lm = hG | Tˆ |E i , (B9) X X l1l2 m l l1l2 For clarity, we first assume the system (1) obeying the periodic BC. G′ = T lm |G i , (B10) To implement the perturbation treatment, we give l1l2 l1l2 m lm the projection operator onto the subspace involved by X ′ ′ the magnon bound-states and the projection operator ′ l m ∗ we have G = ′ ′ hG ′ | T and onto the orthogonal component of the involved subspace. l1l2 l m m l1l2 The unperturbed term Hˆ0 has only two eigenvalues: (i) 2 P ′ ′ 1 1 J lm l m ∗ E = −∆( L − 1) − B ( L − 2) for the L -fold degen- ˆ ′ 0 4 t 0 2 t t h2 = − |Gmi Tl1l2 Tl1l2 hGm | ∆ ′ ′ erated ground-states {|Gli = |l,l +1i : 1 ≤ l ≤ Lt} and mm ll l1l2 1 1 X (ii) El1l2 = −∆( 4 Lt − 2) − B0( 2 Lt − 2) for excited eigen- J 2 states {|E i = |l ,l i : l +1 < l , (l ,l ) 6= (1,L )}. = − G′ G′ . (B11) l1l2 1 2 1 2 1 2 t ∆ l1l2 l1l2 The two-magnon bound-states only involve the subspace l1l2 X U0 spanned by Lt independent ground states |Gli. The projection operator onto U reads as Using the commutation relations of the spin operators, 0 after some algebra, we obtain

Pˆ0 = |Gli hGl| . (B3) ′ Gl l = δl1,l2−2(|Gl1 i + |Gl1+1i) l 1 2 X +δ (|G i + |G i). (B12) l1−2,l2 l2 l2+1 Introducing V0 as the orthogonal complement of U0, the projection operator onto V0 reads as Inserting Eq. (B12) into Eq. (B11), we get 1 2 ˆ ˆ J S = |El1l2 i hEl1l2 | . (B4) h2 = − δl1,l2−2(|Gl1 i + |Gl1+1i)(hGl1 | + hGl1+1|) E0 − El1l2 ∆ El1l2 6=E0 l1l2 X X  +δl1−2,l2 (|Gl2 i + |Gl2+1i)(hGl2 | + hGl2+1|) . (B13) Using the notations Pˆ0 and Sˆ and including the pertur- bations up to second order47, the effective Hamiltonian Since δl1,l2−2 = 1 for (l1,l2) = (m,m + 2) with m ∈ is given as, {1, 2,...,Lt − 2} and δl1−2,l2 = 1 for (l1,l2)=(1,Lt − 1) or (2,Lt), we have Hˆeff = hˆ0 + hˆ1 + hˆ2 = E0Pˆ0 + Pˆ0Hˆ1Pˆ0 + Pˆ0Hˆ1SˆHˆ1Pˆ0. (B5) L J 2 t One can benefit from the fact that the transverse cou- hˆ = − (|G i +|G i)(hG |+hG |). (B14) ˆ 2 ∆ m m+1 m m+1 pling term only contributes to h2 while the modulation m=1 X term only contributes to hˆ1. Using Eqs. (B1) and (B3), the first-order perturbation From Eqs. (B5), (B7), and (B14), one can obtain the reads as effective Hamiltonian 2 ˆ ˆ ˆ ˆ ˆ ′ ′ J h1 = P0H1P0 = |Gmi hGm| H1 |Gm i hGm | . (B6) Hˆeff = − (|Gmi hGm+1| + |Gm+1i hGm|) ′ ∆ mm m X X 1 eff ˆz ′ ′ + (Bm + Bm+1 + V ) |Gmi hGm| . (B15) Since hGm| Sl |Gm i = δmm (δl,m + δl,m+1 − 2 ), we have 0 m ˆ ′ ′ 1 hGm| H1 |Gm i = δmm (Bm + Bm+1 − 2 l Bl) and so X that with the constant P L t 1 1 1 2J 2 hˆ = B + B − 1 B |G i hG | . (B7) V eff =2B − L B − B − ∆ L − 1 − . 1 m m+1 2 l m m 0 0 2 t 0 2 l 4 t ∆ m=1 ! X Xl Xl   8

In order to capture the single-quasiparticle nature of the magnon bound-states, we introduce creation oper- ˆ† ators bm for the magnon bound-states, which flips the spin on the m-th lattice site and the spin on the (m +1)- th lattice site from the vacuum states |0i = |↓↓ ... ↓i. ˆ† ˆ+ ˆ+ MBS ˆ† 0 That is, bm ⇔ SmSm+1 and nm =1 = bm | i ⇔ |n =1,n =1i = Sˆ+Sˆ+ |0i. Therefore, from m m+1 m m+1 Eq. (B15), the magnon bound-states obey the effective single-particle Hamiltonian,

L L J 2 t t Hˆ = − (ˆb† ˆb + H.c.)+ µ ˆb† ˆb , (B16) eff ∆ m m+1 m m m m=1 m=1 X X eff ˆ where µm = Bm+Bm+1+V0 . Since Heff commutes with ˆ ˆ† ˆ ˆ ˆ eff ˆ Nb = m bmbm, i.e. [Heff , Nb] = 0, the term V0 Nb only causes an energy shift. Therefore, we can ignore eff P V0 in µm without changing the spectrum structure of the effective Hamiltonian (B16). FIG. 4. (Color online) The bound-state spectrum for the To calculate the single-particle spectrum of the effec- systems under the periodic BC. The black circles denote the tive Hamiltonian (B16), assuming the wavefunction as eigenenergies for the effective Hamiltonian (B16) and the red dots are the Lt lowest eigenenergies for the original Hamilto- ˆ† 0 ˆ |Ψi = m ψ(m)bm | i, the eigenequation Heff |Ψi = nian (1). The parameters are chosen as λ/J = 0.04, β = 1/3, E |Ψi reads as δ = π/6, Lt = 99, and different values of ∆/J: (a) ∆/J = 1, P (b) ∆/J = 3, (c) ∆/J = 5, (d)∆/J = 100. The eigenenergies eff eff Heff ψ(m)= −J ψ(m +1)+ ψ(m − 1) + µm ψ(m), E/J are added with a constant ∆/J + 2J/∆. = Eψ(m),  (B17) eff 2 eff where, J = J /∆ > 0 and µm = Bm + Bm+1 = with q independent functions φ (m). Then the single- ′ ′ ′ n,k λ cos(2πβm + δ ) with the parameters λ = 2λ cos(πβ) particle spectrum of the effective model is given by diag- ′ and δ = δ + πβ. Obviously, the above eigenequation onalizing the q × q matrix for the effective Hamiltonian. 29 [Eq. (B17)] reproduces the Harper equation with the In Fig. 4, we show the spectrum of the effective Hamil- parameters determined by the anisotropy ∆ and the fre- tonian (B16) and the magnon bound-state bands of the quency β. Therefore, our effective model can be mapped original Hamiltonian (1). Under sufficiently strong in- onto the 2D Hofstadter model that exhibits the well- teractions, the effective model (B16) well describes the known butterfly spectrum. Our calculations do confirm magnon bound-states in the original system (1). the appearance of the butterfly-like spectrum in the two- Under the periodic BC, the Bloch momentum k for magnon bound-state band. the effective Hamiltonian (B16) is a good quantum num- Moreover, the cotranslation operator T (2) commutes (2) ber and (k,δ) makes up a 2D parameter space to define with the effective Hamiltonian, [T ,Heff ] = 0. The ac- the Chern number. We calculate the Chern numbers for tion of T (2) on the states of the effective model reads (2) eff eff the whole spectrum of the effective Hamiltonian using as T ψ(m) = ψ(m + q). Since µm+q = µm , the the same parameters given in Tab. I. The results are (2) (2) quasi-momentum k corresponding to T [T ψ(m) = consistent with the ones in Tab. I for the original sys- ikq e ψ(m)] is the Bloch momentum of the effective model. tem (1). This means that the topological features of the Hence, the quasi-momentum of the cotranslation opera- two-magnon bound-states in the original Hamiltonian (1) tor is a generalization of the Bloch momentum in non- can be described by the effective single-particle Hamilto- interacting 1D systems with translational symmetry to nian (B16). our interacting 1D systems with cotranslational symme- So far, we only consider the system under the periodic try, and the Chern number (8) is an analog to the one 29 BC. However, the BCs may have strong effects on the introduced in the noninteracting 1D systems . spectrum and the emergence of the edge states. Below we Denoting ψn,k = ψn,k(m) as the common eigenstates of consider the system under the open BC and its effective (2) (2) Heff and T , that is Heff ψn,k = Enψn,k and T ψn,k = single-particle Hamiltonian. eikqψ , we have ψ (m)= eikmφ (m) with φ (m + n,k n,k n,k n,k Under the open BC, the dimension of the Hilbert space q) = φn,k(m) due to the Bloch’s theorem. In terms of spanned by the ground states |Gli = |l,l +1i is Lt − 1 φn,k, the eigenequation (B17) can be expressed as rather than Lt for the system under the periodic BC. eff ik −ik The Hilbert space for the magnon bound-states is then Eφn,k(m)= −J e φn,k(m +1)+ e φn,k(m − 1) U1 = {|Gli : 1 ≤ l ≤ Lt − 1}, and the corresponding eff +µm φ n,k(m) (B18) projection operator onto U1 is still Eq. (B3) but with 9

FIG. 6. (Color online) The energy gaps versus the deforming parameter η. We consider the system (C1) under the periodic BC. The parameters are chosen as ∆/J = 100, λ/J = 0.04, β = 1/3, and Lt = 99. FIG. 5. (Color online) The bound-state spectrum for the sys- tems under the open BC. The black circles denote the eigenen- ergies for the effective Hamiltonian (B21) and the red dots are turbation theory, one can not guarantee that such a per- the Lt lowest eigenenergies for the original Hamiltonian (1). turbation treatment retains the topological nature of the The parameters are given as λ/J = 0.04, β = 1/3, δ = π/6, original model. Below, through introducing an auxiliary Lt = 46, and different values of ∆/J: (a) ∆/J = 1, (b) ∆/J = 3, (c) ∆/J = 5, (d) ∆/J = 100. The energies E/J model, we demonstrate the topological equivalence be- are added with a constant ∆/J + 2J/∆. tween the original model and the effective model numer- ically. According to the topological band theory6,7, the gap the summation index l from 1 to Lt − 1. Following the in a topological system is associated with an index that similar procedure of deriving the effective Hamiltonian characterizes the topological nature of the bands below under the periodic BC, we obtain the effective Hamilto- this gap. In fact, such an index is a topological invari- ˆ O ˆO ˆO ˆO nian Heff = h0 + h1 + h2 with ant and it equals to the sum of the Chern numbers of the bands below this gap. Two gapped systems belong L −1 t to the same topological state if they can be continuously hˆO = B + B − 1 B |G i hG | , 1 m m+1 2 l m m deformed from one into the other without gap closing. ! m=1 l And vice versa, when a system with a topological index X X (B19) is continuously deformed into another system with dif- and ferent topological index (i.e. the two systems are not L −2 J 2 t topologically equivalent), the bands will invert and the hˆO = − (|G i+|G i)(hG |+hG |). (B20) 2 ∆ m m+1 m m+1 gap closes during the deforming process. Correspond- m=1 X ingly, a so-called topological occurs in Introducing the creation operators ˆb† (in which the index the gap closing process. m To show the topological equivalence between the orig- m = {1, 2, ··· ,Lt − 1}), the effective Hamiltonian reads as inal model and the effective model, we introduce an aux- iliary model, L −2 L −1 J 2 t t Hˆ O = − (ˆb† ˆb + H.c.)+ µ ˆb† ˆb Hˆ A(η)= Hˆ + η(Hˆ − Hˆ ), (C1) eff ∆ m m+1 m m m eff m=1 m=1 X X ˆ J 2 which is a combination of the original model H described ˆ†ˆ ˆ† ˆ by the Hamiltonian (1) and the effective model Hˆ de- + (b1b1 + bLt−1bLt−1), (B21) eff ∆ scribed by the Hamiltonian (B16). Here, η is the deform- eff with µm = Bm+Bm+1+V0 . The validity of the effective ing parameter. To discuss the two-magnon excitations, Hamiltonian under the open BC is shown in Fig. 5. Hˆ and Hˆeff are both restricted in the two-magnon Hilbert space. Therefore Hˆeff can be explicitly rewritten in the form of Appendix C: Topological equivalence between the J 2 effective Hamiltonian and the original Hamiltonian Hˆ =− (Sˆ+Sˆ+ Sˆ− Sˆ− + H.c.) eff ∆ m m+1 m+2 m+1 m X Since the mapping from the original model to the ef- ˆ+ ˆ+ ˆ− ˆ− + µmSmSm+1Sm+1Sm, (C2) fective model is in the context of the second-order per- m X 10

eff with µm = Bm + Bm+1 + V0 . To explore whether the effective model H and the orig- A A inal model Hˆeff are topological equivalent, we calculate Obviously, Hˆ = Hˆ for η = 0 and Hˆ = Hˆeff for η = 1. Thus, if the parameter η is continuously changed the energy gap versus the deforming parameter η. As an from 0 to 1, the system Hˆ A continuously deforms from example, for β = 1/3, there are three subbands in the magnon bound-state band and therefore there are two Hˆ into Hˆeff . The energy gap between the j-th subband and the (j + 1)-th subband for the magnon bound-state energy gaps δE12 and δE23. As shown in Fig. 6, the en- ergy gaps are nonzero for all η, that is the energy gaps band of Hˆ A can be given as keep open during the whole deforming process. This in- dicates that the original model and the effective model δEj,j+1 = min Ej+1(k,δ) − Ej (k,δ) . (C3) are topological equivalent. k,δ  

∗ Corresponding author. Email: [email protected] 19 F. Grusdt, N. Y. Yao, D. Abanin, M. Fleischhauer, and E. 1 H. Bethe, Zur Theorie der Metalle, Z. Phys. 71, 205 (1931). Demler, Interferometric measurements of many-body topo- 2 M. Wortis, Bound States of Two Spin Waves in the Heisen- logical invariants using mobile impurities, Nat. Commun. berg Ferromagnet, Phys. Rev. 132, 85 (1963). 7, 11994 (2016). 3 J. Hanus, Bound States in the Heisenberg Ferromagnet, 20 X. Chen, Z.-C. Gu, Z.-X. Liu, and X.-G. Wen, Symmetry- Phys. Rev. Lett. 11, 336 (1963). Protected Topological Orders in Interacting Bosonic Sys- 4 T. Fukuhara, A. Kantian, M. Endres, M. Cheneau, P. tems, Science 338, 1604 (2012). Schauß, S. Hild, D. Bellem, U. Schollw¨ock, T. Giamarchi, 21 C. Wang, A. C. Potter, and T. Senthil, Classification of C. Gross, I. Bloch, and S. Kuhr, Quantum dynamics of a Interacting Electronic Topological Insulators in Three Di- mobile spin impurity, Nat. Phys. 9, 235 (2013). mensions, Science 343, 629 (2014). 5 T. Fukuhara, P. Schauß, M. Endres, S. Hild, M. Cheneau, 22 D. J. Thouless, M. Kohmoto, M. P. Nightingale, and I. Bloch, and C. Gross, Microscopic observation of magnon M. den Nijs, Quantized Hall Conductance in a Two- bound states and their dynamics, Nature 502, 76 (2013). Dimensional Periodic Potential, Phys. Rev. Lett. 49, 405 6 M. Z. Hasan and C. L. Kane, Colloquium: Topological (1982). insulators, Rev. Mod. Phys. 82, 3045 (2010). 23 Q. Niu, D. J. Thouless, and Y.-S. Wu, Quantized Hall 7 X.-L. Qi and S.-C. Zhang, Topological insulators and su- conductance as a topological invariant, Phys. Rev. B 31, perconductors, Rev. Mod. Phys. 83, 1057 (2011). 3372 (1985). 8 J. E. Moore, The birth of topological insulators, Nature 24 Z. Wang, X.-L. Qi, and S.-C. Zhang, Topological Order 464, 194 (2010). Parameters for Interacting Topological Insulators, Phys. 9 A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Lud- Rev. Lett. 105, 256803 (2010). wig, Classification of topological insulators and supercon- 25 Z. Wang and S.-C. Zhang, Simplified Topological Invari- ductors in three spatial dimensions, Phys. Rev. B 78, ants for Interacting Insulators, Phys. Rev. X 2, 031008 195125 (2008). (2012). 10 A. Kitaev, Periodic table for topological insulators and su- 26 P. Richerme, Z.-X. Gong, A. Lee, C. Senko, J. Smith, M. perconductors, AIP Conf. Proc. 1134, 22 (2009). Foss-Feig, S. Michalakis, A. V. Gorshkov, and C. Monroe, 11 H. Katsura, N. Nagaosa, and P. A. Lee, Theory of the Non-local propagation of correlations in quantum systems Thermal Hall Effect in Quantum Magnets, Phys. Rev. with long-range interactions, Nature 511, 198 (2014). Lett. 104, 066403 (2010). 27 P. Jurcevic, B. P. Lanyon, P. Hauke, C. Hempel, P. Zoller, 12 Y. Onose, T. Ideue, H. Katsura, Y. Shiomi, N. Nagaosa, R. Blatt, and C. F. Roos, Quasiparticle engineering and and Y. Tokura, Observation of the Magnon Hall Effect, entanglement propagation in a quantum many-body sys- Science 329, 297 (2010). tem, Nature 511, 202 (2014). 13 L. Zhang, J. Ren, J.-S. Wang, and B. Li, Topological 28 T. Graß, C. Muschik, A. Celi, R. W. Chhajlany, and M. magnon in insulating ferromagnet, Phys. Rev. Lewenstein, Synthetic magnetic fluxes and topological or- B 87, 144101 (2013). der in one-dimensional spin systems, Phys. Rev. A 91, 14 M. Pereiro, D. Yudin, J. Chico, C. Etz, O. Eriksson, and 063612 (2015). A. Bergman, Topological excitations in a kagome magnet, 29 L.-J. Lang, X. Cai, and S. Chen, Edge States and Topo- Nat. Commun. 5, 4815 (2014). logical Phases in One-Dimensional Optical Superlattices, 15 M. Punk, D. Chowdhury, and S. Sachdev, Topological ex- Phys. Rev. Lett. 108, 220401 (2012). citations and the dynamic structure factor of spin 30 Pl. Kannappan, in Functional Equations and Inequalities on the kagome lattice, Nat. Phys. 10, 289 (2014). with Applications, Springer Monographs in Mathematics 16 M. Hirschberger, R. Chisnell, Y. S. Lee, and N. P. Ong, (Springer-Verlag US, 2009). Thermal Hall Effect of Spin Excitations in a Kagome Mag- 31 F. Mei, S.-L. Zhu, Z.-M. Zhang, C. H. Oh, and N. Gold- net, Phys. Rev. Lett. 115, 106603 (2015). man, Simulating Z2 topological insulators with cold atoms 17 R. Chisnell, J. S. Helton, D. E. Freedman, D. K. Singh, in a one-dimensional optical lattice, Phys. Rev. A 85, R. I. Bewley, D. G. Nocera, and Y. S. Lee, Topological 013638 (2012). Magnon Bands in a Kagome Lattice Ferromagnet, Phys. 32 S.-L. Zhu, Z.-D. Wang, Y.-H. Chan, and L.-M. Duan, Rev. Lett. 115, 147201 (2015). Topological Bose-Mott Insulators in a One-Dimensional 18 J. Maciejko and G. A. Fiete, Fractionalized topological in- Optical Superlattice, Phys. Rev. Lett. 110, 075303 (2013). sulators, Nat. Phys. 11, 385 (2015). 33 F. Mei, J.-B. You, D.-W. Zhang, X. C. Yang, R. Fazio, S.- 11

L. Zhu, and L. C. Kwek, and particle Precision Measurement of Atomic Scattering Properties, pumping in a one-dimensional shaken optical lattice, Phys. Phys. Rev. Lett. 92, 160406 (2004). Rev. A 90, 063638 (2014). 41 C. Gross, T. Zibold, E. Nicklas, J. Est`eve, and M. K. 34 X. Qin, Y. Ke, X. Guan, Z. Li, N. Andrei, and C. Lee, Oberthaler, Nonlinear interferometer surpasses clas- Statistics-dependent quantum co-walking of two particles sical precision limit, Nature 464, 1165 (2010). in one-dimensional lattices with nearest-neighbor interac- 42 X. Qin, F. Mei, Y. Ke, L. Zhang, and C. Lee, Topolog- tions, Phys. Rev. A 90, 062301 (2014). ical magnon bound-states in quantum Heisenberg chains, 35 P. G. Harper, Single Band Motion of Conduction arXiv: 1602.03217. in a Uniform Magnetic Field, Proc. Phys. Soc. A 68, 874 43 M. Di Liberto, A. Recati, I. Carusotto, and C. Menotti, (1955). Two-body physics in the Su-Schrieffer-Heeger model, Phys. 36 D. R. Hofstadter, Energy levels and wave functions of Rev. A 94, 062704 (2016). Bloch electrons in rational and irrational magnetic fields, 44 A. M. Marques and R. G. Dias, Multihole edge states in Phys. Rev. B 14, 2239 (1976). Su-Schrieffer-Heeger chains with interactions, Phys. Rev. 37 Y. Hatsugai, Edge states in the integer quantum Hall effect B 95, 115443 (2017). and the Riemann surface of the Bloch function, Phys. Rev. 45 R. Resta, Manifestations of Berry’s phase in and B 48, 11851 (1993). condensed matter, J. Phys.: Condens. Matter 12, R107 38 Y. Hatsugai, Chern number and edge states in the integer (2000). quantum Hall effect, Phys. Rev. Lett. 71, 3697 (1993). 46 T. Fukui, Y. Hatsugai, and H. Suzuki, Chern Numbers in 39 O. Mandel, M. Greiner, A. Widera, T. Rom, T. W. H¨ansch, Discretized Brillouin Zone: Efficient Method of Comput- and I. Bloch, Coherent Transport of Neutral Atoms in ing (Spin) Hall Conductances, J. Phys. Soc. Jpn. 74, 1674 Spin-Dependent Optical Lattice Potentials, Phys. Rev. (2005). Lett. 91, 010407 (2003). 47 M. Takahashi, Half-filled Hubbard model at low tempera- 40 A. Widera, O. Mandel, M. Greiner, S. Kreim, T. W. ture, J. Phys. C: State Phys. 10, 1289 (1977). H¨ansch, and I. Bloch, Entanglement Interferometry for

Top View