Topological magnon bound-states in periodically modulated Heisenberg XXZ chains Xizhou Qin1, Feng Mei1, Yongguan Ke1,2, Li Zhang1,2, and Chaohong Lee1,2,∗ 1TianQin Research Center & School of Physics and Astronomy, Sun Yat-Sen University (Zhuhai Campus), Zhuhai 519082, China and 2Key Laboratory of Optoelectronic Materials and Technologies, Sun Yat-Sen University (Guangzhou Campus), Guangzhou 510275, China (Dated: July 19, 2017) Strongly interacting topological states in multi-particle quantum systems pose great challenges to both theory and experiment. Recently, bound states of elementary spin waves (magnons) in quantum magnets have been experimentally observed in quantum Heisenberg chains comprising ultracold Bose atoms in optical lattices. Here, we explore a strongly interacting topological state called topological magnon bound-state in the quantum Heisenberg chain under cotranslational symmetry. We find that the cotranslational symmetry is the key to the definition of a topological invariant for multi- particle quantum states, which enables us to characterize the topological features of multi-magnon excitations. We calculate energy spectra, density distributions, correlations and Chern numbers of the two-magnon bound-states and show the existence of topological protected edge bound-states. Our study not only opens a new prospect to pursue topological magnon bound-states, but also gives insights into the characterization and understanding of strongly interacting topological states. I. INTRODUCTION topological invariant if the system presents a cotransla- tional symmetry: collective translational invariance of a The study of magnons in quantum spin models can multi-particle system. We demonstrate that our topolog- provide fundamental insights into collective excitations ical invariant can be used to characterize the topological and low-energy properties of quantum magnetic sys- feature of multi-magnon bound-states. tems. The quantum Heisenberg chain is one of the typ- Compared with other topological invariants, our defi- ical models of interacting spins. The well-known collec- nition directly associates with the intrinsic symmetry in tive excitations around the ground state of a ferromag- the considered system. The famous Thouless-Kohmoto- netic Heisenberg chain appear as quasiparticle excita- Nightingale-Nijs (TKNN) invariant (i.e. the first Chern tions, the so-called magnons1. Due to the ferromagnetic number)22, which is an integral of the Berry curvature interaction, in addition to continuous scattering states, of Bloch states over the first Brillouin zone, strongly magnon bound-states do appear when the interaction depends on the Bloch states, and thus applies only to is strong enough1–3. Recently, magnon excitations and noninteracting systems. A well-known generalization of their bound-states have been observed in the quantum TKNN invariant for interacting systems is by using the Heisenberg chains comprising bosonic atoms in optical so-called twisted boundary condition23. Another way to lattices4,5. defining topological invariant for interacting systems is 24,25 On the other hand, since the discovery of topologi- the Green-function method . Differently, our topolog- cal insulators, the studies of topological matters have ical invariant is directly defined based upon the cotrans- attracted extensive interests in recent years6–8. Up to lational symmetry. With our topological invariant, we now, theoretical understanding of noninteracting topo- study the two-magnon excitations in a periodically mod- logical phases has been established completely9,10. In ulated Heisenberg XXZ chain and explore the emergence the past few years, the topological feature of single- of topological magnon bound-states. To give insights to magnon excitations has been studied in both theory and the topological origin of these magnon bound-states, we experiment11–17. Beyond the scenario for single-magnon derive an effective Hamiltonian by using many-body de- topological excitations, it is of great importance to clar- generate perturbation theory. Moreover, it is possible ify whether topological bound-states can emerge in a to experimentally examine our predictions via quantum 4,5 26–28 strongly interacting multi-magnon system. In general, spin chains comprising ultracold atoms or ions . arXiv:1602.03217v2 [quant-ph] 18 Jul 2017 strongly interacting topological states pose much difficul- This article is organized as follows. This section in- ties to both theory18 and experiment19. The character- troduces the background and motivation. In Sec. II, ization of strongly interacting topological states is quite we describe our model and discuss its basic proper- different from that of the weakly interacting counter- ties. In Sec. III, we introduce the cotranslational sym- parts20,21. In particular, due to the strong inter-particle metry in interacting multi-particle systems and define a correlations, it is not easy to define a topological invariant topological invariant associates with this symmetry. In and to clarify the interplay between topological features Sec. IV, we study two-magnon excitations in a period- and inter-particle interactions. ically modulated Heisenberg chain. We give the two- In this article, we pursue a new direction to explore the magnon excitation spectrum, calculate the defined Chern topological magnon bound-states in a quantum Heisen- numbers and reveal butterfly-like spectrum of magnon berg XXZ chain. We find that it is possible to define a bound-states. In order to clarify the topological origin of 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 <l <...<l ψ(l1,l2,...,lN ) |l1l2 ...lN i, Xl Xl 1 2 N where the N-magnon wavefunction ψ(l1,l2,...,lN ) = 0 ˆ− ˆ− ˆ−P h | SlN ... Sl2 Sl1 |Ψi. According to Hamiltonian (1), the 1 ˆx ˆy ˆz ˆ± ˆx ˆy N-magnon wavefunction obeys, with the spin- 2 operators (Sl , Sl , Sl ) and Sl = Sl ±iSl for the l-th lattice site. For simplicity, we use the units Hψ(l ,l ,...,l ) of the reduced Planck constant ~ = 1. Here, J and ∆ 1 2 N N are the transverse and longitudinal spin-exchange cou- =−J ψ(l ,...,l ,l + α, l ,...,l ) plings, respectively. In addition to the uniform field 1 j−1 j j+1 N j=1 α=±1 B0, the magnetic fields include periodic modulations X X Bl = λ cos(2πβl + δ) characterized by the modulation +Vl1l2...lN ψ(l1,l2,...,lN ), (2) amplitude λ, the periodic parameter δ and the rational with number β = p/q (in which p and q are coprime integers). N−1 N Similar to the modulated single-particle systems29, the V = −∆ δ −∆δ δ + B , (3) periodic parameter δ may act as an additional dimension l1l2...lN lj +1,lj+1 l1,1 lN Lt lj j=1 j=1 for constructing the torus and defining the topological X X invariant. However, differently, the single-particle quasi- where a constant energy shift has been removed. momenta for interacting multi-particle systems are no In the N-magnon sector, the cotranslational operator longer good quantum numbers. This difference imposes Tq(τ) can be defined as a great challenge to define topological invariants for in- teracting multi-particle systems. Tq(τ)ψ(l1,...,lN )= ψ(l1 + τq,...,lN + τq), (4) For convenience, we concentrate on studying magnon where τ is an integer. Here we assume the lattice size Lt excitations over the fully magnetized state |↓↓ ... ↓i with is commensurate with the periodic modulation (i.e. Lt = ′ ′ all spins pointing downwards. For a positive and suffi- qs with s an integer). As Tq(τ)Tq (τ ) = Tq(τ )Tq(τ) = ′ −1 ciently large B0, the ground state is the fully magne- Tq(τ + τ ) and [Tq(τ)] = Tq(−τ), an Abelian group tized state |↓↓ ... ↓i. The excited states over |↓↓ ... ↓i can be constructed. We call this Abelian group as the can be obtained by flipping spins. Regarding the ground cotranslation group. Since Vl1+q,l2+q,...,lN +q = Vl1l2...lN , state as a vacuum state, the excited states can be viewed one can easily demonstrate that Tq(τ)Hψ(l1,...,lN ) = as a gas of quasi-particles called magnons. Thus the HTq(τ)ψ(l1,...,lN ) holds for an arbitrary N-magnon transverse and longitudinal spin-exchange couplings re- wavefunction ψ(l1,...,lN ). Therefore the Hamiltonian spectively represent the hopping and interaction between commutes with all cotranslation operators, that is, it is magnons, and the magnetic fields (B0,Bl) describe the invariant under the cotranslation group on-site energy of magnons. Since the Hamiltonian (1) [T (τ)]−1HT (τ)= H, (5) respects the U(1) symmetry under rotation transforma- q q ˆ± ±iθ ˆ± which represents a cotranslational symmetry. tions, Sl → e Sl , the total spin along z-direction Sˆz = Sˆz is conserved.