arxiv: v2 [cond-mat.str-el] 8 Aug 2017

Transcription

1 Spinon decay in the spin-/ Heisenberg chain with weak next nearest neighbour exchange Stefan Groha The Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford, OX 3NP, United Kingdom arxiv: v [cond-mat.str-el] 8 Aug 07 Fabian H.L. Essler The Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford, OX 3NP, United Kingdom Abstract. Integrable models support elementary excitations with infinite lifetimes. In the spin-/ Heisenberg chain these are known as spinons. We consider the stability of spinons when a weak integrability breaking perturbation is added to the Heisenberg chain in a magnetic field. We focus on the case where the perturbation is a next nearest neighbour exchange interaction. We calculate the spinon decay rate in leading order in perturbation theory using methods of integrability and identify the dominant decay channels. The decay rate is found to be small, which indicates that spinons remain well-defined excitations even though integrability is broken.

2 Spinon Decay. Introduction Integrable many-particle quantum systems are special in that they support stable elementary excitations. These typically are related in very complicated ways to the basic degrees of freedom. For example, in the Heisenberg antiferromagnet the elementary excitations are interacting spin-/ objects called spinons [, ]. Crucially, these elementary excitations are protected from decay into multi-particle excitations by the existence of local integrals of motion, even in cases where decay is kinematically allowed. In such situations application of integrability breaking perturbations has the immediate effect of inducing particle decay, and an important question is how large the corresponding decay rates are. If they are small, the elementary excitations of the integrable model will remain a good basis for describing the physics of the perturbed model. Such questions have been investigated in some detail for integrable quantum field theories [3, 4, 5, 6]. The case of integrable lattice models is considerably harder, and to the best of our knowledge has not been investigated so far. The added difficulty compared to field theory cases is that the description of the ground and excited states is more complicated see below). The question of what effects weak integrability breaking perturbations have on the excitation spectrum of lattice models is also of importance in so-called mobile impurity approaches to the calculation of threshold singularity exponents in lattice models [7]. As pointed out in Ref. [8] in the context of the Hubbard model, there exist different formulations of mobile impurity models [9], which correspond to different choices of bases of elementary excitations. One may argue that for integrable models the integrable basis of elementary excitations ought to be the preferred choice. An obvious question is then whether this remains the case even if integrability is weakly broken. This is intimately related to how large the decay rate of the excitations is once a perturbation is applied. For the Hubbard model the available integrable model technology [0] does not currently permit to answer this question. In this work we therefore consider the simpler case of the spin-/ Heisenberg XXZ chain of length L in a magnetic field h HJ,, h) = J L S x j Sj+ x + S y j Sy j+ + Sz j Sj+) L z h Sj z. ) j= j= Here J < 0 and S α j are spin operators with commutation relations [S α i, S β j ] = iδ ijε αβγ S γ i. ) The spectrum of ) is gapless for and h < J + ) [, ]. The model is integrable and elementary excitations over the ground state carry S z quantum number / and are known as spinons [,,, ]. A simple way of perturbing the model away from the integrable point is by introducing a next nearest neighbour interaction δh = κ j S z j S z j+. 3) This interaction destroys integrability, but still commutes with the total spin operator along the z-axis S z = j Sz j. Hence the z-component of the total spin remains a good quantum number. In presence of the perturbation spinons cease to be exact elementary excitations and we expect them to acquire a finite life-time. Using Fermi s

3 Spinon Decay 3 golden rule for small perturbations the decay rate can be expressed in the form Γ = π f Mi f) ρ f E i ) δ pf,p i 4) where E f and p f E i and p i ) are the energy and momentum of the final initial) state [3], ρ f E) the density of states of the final state and the matrix element M is given by M = f δh i. 5) We are interested in the case where the initial state is an exact one-spinon eigenstate of ), while the final state is any exact eigenstate of the unperturbed system. Most of our analysis will focus on the isotropic Heisenberg model at = and h > 0. Other values of can be treated in the same way and h < 0 by considering the spin overturned sector. The outline of this paper is as follows. Section presents a brief summary of the Bethe Ansatz solution of the Heisenberg model. We then describe the excited states that contribute to the decay rate in section 3. We then use the Algebraic Bethe Ansatz to obtain explicit expressions for the matrix elements describing the spinon decay, cf. 4. In section 5 we then numerically determine the contributions of various decay channels to the decay rate. We end with a discussion of our results in section 6.. Bethe Ansatz solution of the XXX-chain.. Coordinate Bethe Ansatz Eigenstates of the XXZ Hamiltonian ) can be constructed by means of the Bethe Ansatz [4] for any value of the anisotropy. As S z commutes with the XXZ Hamiltonian and with the perturbation δh, it is convenient to work in a sector with a fixed number of down-spins N with respect to the ferromagnetic state 0 := Sz = L = L j, 6) which will be used as a reference state in the following. Energy eigenstates with N down-spins take the form N = L j,...,j N = j= aj,..., j N ) S j a 0, 7) where j < j < < j N L. The wave functions have Bethe Ansatz form [,, 5, 6] aj,..., j N ) = ) N ) [P ] A P exp i k P a j a, P S N a= A P = ) e ik Pa +k P ) b + e ik Pa. 8) a<b a=

4 Spinon Decay 4 The energy of the state with wave function 8) is given by E = J N ) L cos k a ) h N. 9) a=.. Bethe equation for the XXX model Imposing periodic boundary conditions on the wave functions 8) leads to quantization conditions for the wave numbers k a known as Bethe Ansatz equations e ikal = a b [ eika e ika+ik b ] e ik b e ik a+ik. 0) b From here on we set =. It is convenient to introduce rapidity variables λ a defined by e ika = λ a i/ λ a + i/. ) In terms of the rapidity variables the Bethe Ansatz equations read ) L λa i/ = λ a + i/ b= b a λ a λ b i λ a λ b + i, a =,..., N. ) A standard way of analyzing ) is by employing the string hypothesis. This assumes that all solutions of ) are composed of strings of the form λ n,j α = λ n α + i n + j) + δn,j α, j =,..., n. 3) Here δα n,j are deviations from ideal strings and are assumed to be exponentially small in system size. Let us now consider a solution to ) that contains M n strings of length n with corresponding string centres λ n α this implies that n M nn = N). Substituting 3) into ) and neglecting the deviations we obtain a set of coupled equations for the set {λ n α}. Taking logarithms we arrive at ) λ n Lθ α = πiα n + θ nm λ n α λ m β ). 4) n m,β) n,α) Here Iα n are integer or half-odd integers numbers arising from taking logarithms), θx) = arctanx), and { θ x ) n n + j= θ nm x) = θ ) x j for m = n θ ) ) ) x n m + θ x n m θ x n+m + θ x ) n+m for m n. 5) Equations 4) are called Takahashi s equations. They relate the solutions of the BAE to a set of integer of half-odd integer numbers, which therefore can be considered as quantum numbers of our problem. The permitted ranges of the Iα n are [] [ ] Iα n L minm, n) δ n,m ) M m. 6) m=

5 Spinon Decay 5 Energy and momentum of solutions to 4) are given by M m E = πjan λ m β ) + mh ) hl, P = π + where we have defined m= β= M m m= β= πi m β L 7) a n x) = n π x + n/). 8) All solutions to Takahashi s equations correspond to highest weight states of the spin SU) algebra [] S + {λ n α} = 0. 9) A complete set of energy eigenstates is then obtained by acting with the spin lowering operator on these highest weight states S ) m {λ n α }, m = 0,,..., L n nm n. 0) 3. Low lying excitations and spectrum In order to have access to single-spinon excitations we need to consider odd chain lengths L. For even values of L the lowest excitations involve at least two spinons []. 3.. One particle and one hole excitations For odd L with an odd number N of down spins there are two degenerate lowest energy states. They are obtained by considering real solutions -strings) to the Bethe Ansatz equations and choosing either I α = N + α, α =,..., N, ) or I α = N + + α, α =,..., N, ) The corresponding configurations of half-odd integers for N = look as follows: I α I α The energy density eh) of these two states in the thermodynamic limit L, N, N = n = fixed 3) L

6 Spinon Decay 6 can be expressed in terms of the solution of a linear integral equation for the root density ρ λ), cf. Ref. [] B ρ λ) = a λ) dη a λ η) ρ η). 4) B Here the integration boundary B is determined by the density of down spins n through B The energy per site is then given by where eh) = B B B dλ ρ λ) = n. 5) dλ ρ λ)ε 0) λ), 6) ε 0) λ) = Jπa λ) + h. 7) The two states above are particular limits of one-parameter particle-like and holelike excitations. The particle excitation corresponds to I α configurations of the form whereas the hole-like excitation is obtained by promoting one half-odd integer I α in the ground state configuration to the Fermi edge that has one fewer half-odd integer: I p I α I h I α Both types of excitations involve a single parameter: I p for the particle excitation and I h for the hole excitation. For asymptotically large system sizes L the energies and momenta of these excitations are given by [] E p = Le + ε λ p ) + o), λ p > B, 8) λ p P p = π + π dλ ρ λ) + OL ), 9) 0 E h = Le ελ h ) + o), λ h < B, 30) P h = π π λ h 0 dλ ρ λ) + OL ), 3) where the dressed energy ɛ λ) is a solution to the linear integral equation ε λ) = ε 0) λ) B dµ a µ λ) ε µ). 3) π B The rapidities λ p and λ h are continuous parameters above and below the Fermiedge respectively. They are related to the parameters I p and I h through Takahashi s equations 4). The excitation energy for a one-spinon excitation can now be extracted by simply subtracting the extensive part of the energy which equals the ground state energy

7 Spinon Decay 7 E sp Le h =0 h = h = P sp /π Figure : Energy-momentum dispersion of the one-parameter excitations for different values of the magnetic field h. per site of the Heisenberg chain), which allows us to extract the spinon energy and momentum ɛ s λ) = ε λ), p s λ) = π + π sgn λ B) λ 0 dµ ρ µ). 33) The corresponding dispersion relation is plotted for several values of magnetic field h in Fig. where the value for the magnetic field is fixed by imposing ε B) = 0. We note that by construction the spinon dispersion is identical to the one extracted from the two-spinon excitation of the Heisenberg model with even chain lengths L, apart from a shift in momentum by π. 3.. Excitations involving several particles and/or holes As δh commutes with S z, the decay of the single-particle hole) excitation described above can only involve excited states with the same S z quantum number. These are obtained in the following ways: i) One can consider solutions of Takahashi s equations only involving -strings. These will involve additional particle-hole excitations on top of the -spinon excitation constructed above. ii) One can consider solutions of Takahashi s equations involving n-strings with n. As a result of the magnetic field these excitations have a gap. iii) One can consider excitations of the form 0) that are not SU) highest-weight states. These again have a gap for h > 0 because [S, HJ,, h)] = hs 34) As we are dealing with an interacting theory, this leaves us with an infinite number of possible decay channels, i.e. even to first order in perturbation theory in κ, a single

8 Spinon Decay 8 spinon can decay into excitations involving 3, 5, 7,... particles. As in one dimension the accessible phase space shrinks with the number of particles involved [7], it is reasonable to assume that the dominant decay channels will involve excitations with low numbers of particles. In the following we will focus on excitations involving 3 particles. We have considered a class of five-particle excitations where we excite two particle- and hole-type excitations in addition to the one-spinon excitation, and found the corresponding decay rate to be smaller see section 5) pph-excitation This excitation involves only -strings and corresponds to configurations of the half-odd) integers I α looking as follows I p I h I p I α States of this kind can be thought of as a sub-class of 3-spinon excitations that involves two particles and one hole, which are parametrized by I h and I p, respectively or equivalently by the corresponding rapidities λ h, λ p, ). Energy and momentum of this excitation are given by E pph = Le + ɛ s λ p ) + ɛ sλ p ) + ɛ sλ h ) + o), 35) P pph = p s λ p ) + p sλ p ) + p sλ h ) + OL ), 36) where ɛ s λ) and p s λ) are defined in 33). The excitation energy is obtained subtracting the ground state energy, and the corresponding continuum of 3-spinon excited states is shown in Fig. a. The grey shading reflects the density of excitations at given values of energy and momentum. Darker regions correspond to higher densities. The intensity of the shading is obtained by considering large but finite L and varying I h and I p, over all allowed values for a given excitation, and calculating approximate values of λ h, λ p, by solving the equation p s λ h ) = πih L, p sλ p j ) = πip j L, j =,. 37) The corresponding approximate excitation energy is then obtained by substituting these values into 36). Each set {I h, I p, } provides one point in the P pph -E pph -plane and the collection of all these points generates a shading that reflects the density of states. We see that for momenta π + m) < p sλ) < π 3 m) decay of the -spinon excitation is kinematically forbidden, while it is allowed for some values in the regions p s λ) > π 3 m) and p sλ) < π + m) phh-excitation This excitation involves only -strings and corresponds to configurations of the half-odd) integers I α looking as follows I h I h I p I α

9 Spinon Decay 9 E pph Le P pph /π 5 4 a) E phh Le P phh /π.5 b) E ppp Le 3 E hhh Le P ppp /π c) P hhh /π d) Figure : Excitation continua grey) for a) two particle one hole b) one particle two holes c) three particle and d) three holes at magnetization m = 3/0. The -spinon dispersion is shown in red. The shading of the continuum reflects the density of states see main text). Decay of the single spinon is kinematically allowed in part of the Brillouin zone. States of this kind are a sub-class of 3-spinon excitations that involves one particle and two holes, which are parametrized by I p and I h, or equivalently by the corresponding rapidities λ p, λ h,. Energy and momentum of this excitation are given by E phh = Le + ɛ s λ p ) + ɛ s λ h ) + ɛ s λ h ) + o), 38) P phh = p s λ p ) + p s λ h ) + p s λ h ) + OL ), 39) where ɛ s λ) and p s λ) are defined in 33). The excitation energy is again obtained by subtracting the ground state energy and is shown as a function of the total momentum in Fig. b ppp-excitations This excitation involves only -strings and corresponds to configurations of the half-odd) integers I α looking as follows I p I p I p 3 I α

10 Spinon Decay 0 States of this kind are a sub-class of 3-spinon excitations that involves three particles. Energy and momentum of this excitation are E ppp = Le + P ppp = 3 j= ɛ s λ p j ) + o), 40) 3 p s λ p j ) + OL ), 4) j= where ɛ s λ) and p s λ) are defined in 33). The excitation energy is again obtained by subtracting the ground state energy and is shown as a function of the total momentum in Fig. c hhh-excitations This excitation involves only -strings and corresponds to configurations of the half-odd) integers I α looking as follows I h I h I h 3 I α States of this kind are a sub-class of 3-spinon excitations that involves three holes. Energy and momentum of this excitation are E hhh = Le + P hhh = 3 ɛ s λ h j ) + o), 4) j= 3 p s λ h j ) + OL ), 43) j= where ɛ s λ) and p s λ) are defined in 33). The excitation energy is again obtained by subtracting the ground state energy and is shown as a function of the total momentum in Fig. d Excitations involving a single -string We now turn to the simplest excitation involving a single -string. This corresponds to solutions of 4) with M = N, M = and configurations of the half-odd integers I α, I of the kind I p I α We note that the permitted values for I have range I [L N]. 44) The excitation is parametrized by the two half-odd integers I p, I or equivalently the corresponding rapidities λ p, λ s. Energy and momentum of this excitation are given by E sp = Le + ɛ s λ p ) + ε λ s ), λ p > B, 45) P sp = p s λ p ) + p λ s ), 46) I

11 Spinon Decay where ε and p λ) are given by [] ε λ) = h + dµ a λ µ) ε µ), λ ) p s λ) = θ µ >B B B dµ θ λ µ) ρ µ). 47) Excitation continua that encompass the two-particle continuum 46) are obtained by adding particle-hole excitations, e.g. E s3ph = Le + P s3ph = 3 ɛ s λ p j ) + j= 3 p s λ p j ) + j= ɛ s λ h k) + ε λ s ), k= λ h j < B < λ p k, 48) p s λ h k) + p λ s ). 49) k= The continuum 49) is shown in Fig. 3 for several magnetizations. We see that the single spinon excitation cannot decay into the -string excitation for kinematic reasons. If we keep on adding particle-hole excitation at small magnetic fields decay of the - spinon excitation will eventually become kinematically allowed. However, the decay rate is then expected to be negligible on the basis of aforementioned phase-space arguments, cf. Ref. [7] Excitations involving longer strings Excitations involving longer strings have larger gaps at finite magnetic fields []. We expect contributions from decay channels involving such excitations to be small for the same reasons we put forward in the - string case above Excitations that are not highest weight states As mentioned above, excitations which are not highest weight states have gaps that are proportional to the magnetic field h. Nevertheless, decay of a single spinon into excitations that are not highest weight states will generally be allowed at sufficiently small h. As an example let us consider highest-weight states with M = N, M n = 0. The lowest energy states in this sector correspond to integers I α = N + α, α =,..., N, 50) or I α = N + α, α =,..., N, 5) I α I α In complete analogy to our discussion above, these can be viewed as particular limits of a -spinon excitation. Acting with the spin lowering operator gives a -parameter

12 Spinon Decay E sppphh Le E sppphh Le P sppphh /π 4 3 a) E sppphh Le E sppphh Le P sppphh /π b) P sppphh /π c) P sppphh /π d) Figure 3: One two-string and 3ph excitation continuum grey) at magnetization a) m = 0.0, b) m = 0.03, c) m = 0.05 and d) m = 0.07, the corresponding magnetic fields are h = 0.09, h = 0.6, h = 0.4 and h = 0.58 respectively). The -spinon dispersion is shown in red. The shading of the continuum reflects the density of states see main text). Decay of the single spinon is kinematically not allowed. excited state with a dispersion that equals the -spinon dispersion shifted upwards in energy by h. Hence decay of our -spinon excitation into this particular descendant state is kinematically not allowed. However, if we add an additional particle-hole pair the situation changes. Let us consider configurations of integers such as I p I h I p I α States of this kind can be thought of as a sub-class of 3-spinons. Energy and momentum of the excitation obtained by acting with the spin-lowering operator on this state are given by E pph desc = Le + ɛ sλ p ) + ɛ sλ p ) + ɛ sλ h ) + h + o), 5) P pph desc = π + p sλ p ) + p sλ p ) + p sλ h ) + OL ). 53) Inspection of Fig. a shows that decay of the -spinon excitation into this continuum is kinematically allowed at sufficiently weak fields. However, as shown in Appendix A this decay is strongly suppressed for large system sizes L.

13 Spinon Decay 3 4. Algebraic Bethe Ansatz and Derivation of the Matrix Element 4.. Algebraic Bethe Ansatz In order to determine decay rates we require matrix elements of the perturbing operator between energy eigenstates. These can be obtained using the Algebraic Bethe Ansatz []. In the following we will first consider the XXZ case with anisotropy parameter = cos γ and only later specialize to the isotropic limit =. A key object is the monodromy matrix ) Aλ) Bλ) T λ) =, 54) Cλ) Dλ) where A, B, C, D are operators acting on the Hilbert space of the chain and λ is known as spectral parameter. The mondromy matrix fulfils the Yang-Baxter equation Rλ µ)t λ) ) T µ)) = T µ))t λ) )Rλ µ), 55) where the R-matrix has the form Rλ, µ) = 0 bλ, µ) cλ, µ) 0 0 cλ, µ) bλ, µ) 0. 56) Here we have defined bλ) = sinhλ) sinhλ + iγ) cλ) = sinhiγ) sinhλ + iγ) 57) 58) The Yang-Baxter algebra determines intertwining relations for the operators A, B, C, D. Eigenstates of ) can be constructed as λ,..., λ N = Bλ i ) 0, 59) where the set of rapidities {λ i } i {,...,N} are solutions to the Bethe equations The reference state 0 satisfies i= aµ j ) dµ j ) = bµ k µ j ) bµ j µ k ) k j, j =,..., N. 60) Aλ) 0 = 0 Dλ) 0 = dλ) 0 Cλ) 0 = 0. 6) The functions aµ) and dµ) are given by aµ) = dµ) = bµ i γ ) ) L. 6)

14 Spinon Decay 4 The isotropic limit = corresponds to taking γ 0, while rescaling the spectral parameters µ j = γλ j, λ j fixed. 63) This recovers the Bethe Ansatz equations ) from 60). The global spin lowering operator is obtained as [] i lim λ Bλ) = λ S. 64) 4.. Determinant Formulas for Matrix Elements in the XXZ chain The Algebraic Bethe Ansatz provides a convenient setting for calculating scalar products as well as the norm of Bethe states [8, 9, 0]. Matrix elements can be analyzed by utilizing the expression of local spin operators σj z in terms of the operators A, B, C, D, cf. Ref. []. With the help of these relations matrix elements of spin operators Sj α between eigenstates of the XXZ Hamiltonian were derived in Ref. [], and general operators were considered in Ref. []. Explicit expressions for the operator Sj zsz j+ were obtained in Ref. [3]. Following the derivation of Ref. [3] we obtain see Appendix B for details) λ,..., λ N Sj z Sj+ z µ,..., µ M = j Le ip {λ}+ip {µ} δ P{λ},P {µ} δ M,N { φiγ) n ) A n det G ) n + B n ) + det G ) n + B n ) k φλ k + iγ/) φµ k + iγ/)φµ k iγ/) 3 a>b φλ b λ a ) b>a φµ b µ a ) ) det G )) det G ))) k φλ k + iγ/) 3 φµ k iγ/) 3 a>b φλ b λ a ) a<b φµ b µ a ) A nm [det G nm + B nm ) det G nm )]. 65) n= m n Here P {λ} is the total momentum of the state parametrized by the rapidities {λ j } P {λ} = [ N λj + iγ π + ln λ j iγ j= )], 66)

15 Spinon Decay 5 and φiγ) H ab = aµ b ) φλ k µ b + iγ) dµ b ) ) φλ k µ b iγ), 67) φλ a µ b ) k a k a G ) n ) ab = G ) n ) ab = { Hab b n φiγ)φλ a) b = n, 68) φλ a iγ ) φλ a+ iγ ) { Hab b n x= φiγ) x φλ a x)φλ a x+iγ)) iγ B ) n ) ab = δ bn )dµ b ) i n φµ i µ b + iγ)φµ b + iγ/) φ iγ) φiγ) φiγ) φµ d iγ/)φµ d + iγ/) + d n b = n, 69) d m,n ) φ µ d iγ/) φµ d iγ/) φ λ b + iγ/) φλ b + iγ/) + φiγ), 70) φµ m + iγ/)φµ m iγ/) b B n ) ) ab = δ bn ) dµ b) φµ i µ b iγ)φµ b + iγ)φiγ). 7) φλ a + iγ/)φλ a iγ/) i n G nm ) ab = H ab φiγ)φλ a) φλ a iγ ) φλ a+ iγ ) φiγ) x φλ a x)φλ a x+iγ)) x= iγ B nm ) ab = δ bm ) δ bn )dµ b ) i m,n b m, n b = m b = n φµ i µ b iγ)φµ b + iγ/) φiγ) φλ a iγ/)φλ a + iγ/) 7). 73) Finally, the function φλ) is given by { λ = φλ) = sinhλ) <. 74) In the isotropic case = of interest to us the matrix element in the rescaled rapidity variables 63) is obtained by simply setting γ = in the above expressions. 5. Decay rates We are now in a position to compute the rates of decay of the one-spinon excitation into the various multi-particle excitations considered above. In practice the calculation is carried out in a large, finite volume L. The energy eigenstates are of the form 59) and involve N rapidity variables, which constitute a solution to the Bethe Ansatz equations. We will denote the states corresponding to the one-spinon and the twoparticle one-hole continuum by N; J, N; I p, Ip, Ih, 75)

16 Spinon Decay 6 Our notations for the respective energies 7) are EJ), E pph I p, Ip, Ih ). 76) Here J and I p,, Ih denote the half-odd integers corresponding to the spinon and the particles/hole respectively. Other excitations are labelled analogously. The decay rate is then given by L Γ sp pph p J ) κ γ sp pph p J ) = πκ N; I p, Ip, Ih Sj z Sj+ N; z J I p,ip,ih j= δ J,I p +Ip Ih δepph I p, Ip, Ih ) EJ)), 77) where we have used 7) to simplify the momentum conservation delta function. The momentum p J of the initial spinon excitation is given by { πj p J = L if J > N, πj L if J N. 78) We regularize the delta function expressing energy conservation by δ η x) = πη e x η ), 79) where lim η 0 δ η x) = δx). For very small η, but still with a sufficient number of final states in the regime where δ η x) is large, we expect the result to be close to the answer in the thermodynamic limit. For 77) to be finite in the thermodynamic limit, the matrix elements should scale as L. As shown in Fig. 4, the decay in L is very slightly faster than L and is compatible with the functional form M pph L ) = a + b/l L c, 80) where c is a very small exponent. In the range of lattice lengths accessible to us, equally good fits can be obtained by replacing L c by ln L) c in 80). The situation is analogous to that for the dynamical structure factor [4, 5, 6, 7, 8, 9, 30, 3]. For the latter it was shown that in order to obtain finite results in the thermodynamic limit, an infinite summation over states that contain additional particle-hole pairs located at the Fermi points ±B was required. On the other hand, the result obtained by working at a fixed value of L 000 and not carrying out this summation was found to give an excellent approximation to the thermodynamic limit. We expect the decay rate to behave in an analogous way. In the following we determine the contributions of the 3-particle excitations described in section 3 to the decay rate for finite system sizes in the range L We have verified that taking into account states with one additional particle-hole excitation gives only small corrections. We now fix L and then compute 77) for several values of the broadening η. The decay rates into the excitations considered in section 3 are shown in Figs 5a, 5b, 5c and 5d. We see that the dominant decay channel for a one spinon excitation is decay into a 3-spinon excitation of pph type. We have argued above excitations involving higher numbers of particles should give smaller contributions to the decay rate. In order to check this assumption we have calculated the decay rate into a 5-spinon excitation

17 Spinon Decay Fit Numerical data L M L Figure 4: Scaling of matrix elements with system size L. The initial state contains one spinon with momentum p = 5.5, while the final state contains three high-energy excitations with momenta p p = 4.5, p p = 4.65, p h = The fit is to the functional form 80) with a = 0.4, b = 5.6, c = We have also considered additional excitations around the Fermi points in the final state as well as different momenta and found similar behaviour. of type ppphh, which we expect to provide the largest contribution among the 5- spinon excitations. The result is shown in Fig. 6b. As expected the contribution is small. Moreover, it is mostly due to umklapp-type terms in the pph-channel, meaning particle-hole type excitations around the Fermi sea on top of the pph-type excitations cf. Fig. 6a). It is clear from Fig. 5 that all other 3-spinon decay channels can be neglected compared to the pph one. Moreover, the decay rate coefficient is of order unity, which means that the decay rate itself is small and proportional to the square of the strength of the integrability breaking perturbation. 5.. Extrapolation to η = 0 The results of the previous section are for finite values of the system size and require a small, finite regularization parameter η. We will consider the extrapolation of these results to the thermodynamic limit and η = 0. As we have seen above, the matrix elements of the perturbing operator scale as L up to corrections that decay algebraically with very small exponent or logarithmically, cf. Fig. 4. We expect that in order to take the thermodynamic limit, one would have to sum over an infinite number of particle-hole excitations at the Fermi points, in analogy with available results for the spin-spin correlation function [7, 8, 9, 30, 3]. Our situation is more complicated as we need to consider excited states with several elementary excitations at finite energies and summing over an infinite number of particle-hole excitations on top of these is beyond the scope of our work. However, we note that the main source of finite-size effects in our calculation is the necessity to have a sufficiently large value of the broadening η. This is required in order to obtain a good approximation to density of final states ρ f E i ). This imposes a restriction η > η 0 L) = const/l. Importantly,

18 Spinon Decay 8 γsp pph γsp ppp η = η = 0.00 η = P sp /π a) η = 0.0 η = η = P sp /π c) γsp phh γsp hhh η = 0.05 η = 0.0 η = P sp /π η = 0.05 η = 0.05 η = 0.0 b) P sp /π d) Figure 5: Coefficients of the decay rate for magnetization m = 3/0 and several values of the broadening η for a) pph-type processes for L = 65, b) phh-type processes for L = 435, c) ppp-type processes for L = 455, d) hhh-type processes for L = 675. γsp pph pph without umklapp pph with umklapp P sp /π a) γsp ppphh η = 0.05 η = 0.05 η = P sp /π b) Figure 6: Coefficients of the decay rate for magnetization m = 3/0 for a) pphtype processes with and without umklapp terms for η = 0.00 at L = 45 and b) ppphh-type processes for L = 55 and several values of the broadening η. η 0 L) tends to zero much faster than M L. This allows us to extrapolate our results to η = 0 as follows. We construct a smooth interpolation function M cont pph p, zp, zp ) for the matrix element multiplied by L, and then turn the sums over Bethe Ansatz

19 Spinon Decay 9 half-odd) integers into integrals using the Euler-Maclaurin sum formula. Taking the limit η 0 results in γsp pphp) extra = dz p dz p 4π δɛ sz p ) + ɛ sz p ) + ɛ sz p + zp p) εp)) D D M cont pph p, z p, zp ), 8) where D is the domain where the one-spinon exciation exists the interval [π/5, π π/5] for magnetization m = 3/0). We stress again that we do not claim that this integral is the exact form one would get after summing all particle-hole pairs in the thermodynamic limit, but that the result of such a calculation is expected to be numerically very close to what is obtained here. One of the integrals can be carried out using the delta-function, which gives γsp pphp) extra = dz p 4π D M cont pph p, zp, z) ε z) + ε z p, 8) + z p) where z is the solution of the equation ɛ s z p ) + ɛ sz) + ɛ s z p + z p) εp) = 0. Carrying out the remaining integral numerically leads to the result shown in Fig γsp pph P sp /π Figure 7: Decay rate to the pph-channel for η 0 at magnetization m = 3/0. The wiggles are due to the finite size effects of the interpolated matrix element at L = Density of kinematically allowed states and finiteness in the thermodynamic limit A simpler quantity of interest is the density of final states to which transitions from the -spinon excitation are kinematically allowed. For free fermions this density of states exhibits a van Hove singularity that leads to logarithmic divergence at the threshold [33]. In the thermodynamic limit the pph channel density of kinematically relevant states is given by ρ pph p) = dz p 4π D ε z) + ε z p + z p) 83)

20 Spinon Decay 0 where z is the same as in 8). Analogous expressions hold in the other channels. Results for the various possible types of 3-spinon final states are shown in Fig. 8. We see that densities of states are finite and do not display the kind of singularity encountered for free fermions. DOSsp pph DOSsp ppp P sp /π a) P sp /π c) DOSsp phh DOSsp hhh P sp /π b) P sp /π d) Figure 8: Density of states for the a) pph-channel b) phh-channel, c) ppp-channel and d) hhh-channel in the thermodynamic limit at magnetization m = 3/0. 6. Conclusions We have considered decay rates of the elementary spinon excitation in the spin-/ Heisenberg XXX model in a magnetic field perturbed by a weak integrability breaking interaction κ j Sz j Sz j+. We have argued that the leading contribution arises from three spinon decay and have determined the corresponding rate. The latter is found to be small, indicating that spinons remain long-lived excitations in the non-integrable theory. Decay of elementary string excitations can be analyzed in an analogous fashion. This would be particularly interesting to do in the attractive regime < 0 of the XXZ chain in a field, where they play an important role in the dynamics. 7. Acknowledgments We are grateful to J.-S. Caux, L. Glazman and R. Pereira for helpful discussions. This work was supported by the EPSRC under grant EP/N0930X/ and by the Clarendon

21 Spinon Decay Scholarship fund SG). Appendix A. Matrix elements and suppression for non-highest weight states We want to consider the normed matrix element of δh between a highest weight state and a non-highest weight state {λ} S + ) k j Sz j Sz j+ {µ} {λ} S+ S ) k {λ} {µ} {µ} A.) where {λ} and {µ} are highest weight Bethe ansatz states. We see immediately from the commutation relation [S +, S z j ] = S + j A.) and from the relation for highest weight states S + {µ} = 0 that for k 3 the matrix element is exactly 0. Furthermore inserting the cyclic shift operator cf.[]) exp i ˆP ) and using e i ˆP S a j e i ˆP = S a j+ a = z, +, A.3) and the fact that the highest weight states are eigenstates of the shift operator with eigenvalue exp ip {λ} ), where P{λ} is the momentum of the highest weight state, we obtain {λ} e i ˆP e i ˆP S + ) k j S z j S z j+ {µ} = e ip {λ} P {µ} ) {λ} S + ) k j S z j S z j+ {µ} A.4) and therefore we see that the momenta have to conincide. Using the commutation relation A.) we obtain for the normed matrix element for k = and k = : {λ} S + j Sz j Sz j+ {µ} {µ} {µ} {λ} S+ S + S S {λ} = {λ} j S+ j Sz j+ + Sz j S+ j+ {µ} L N ) {µ} {µ} {λ} {λ} A.5) {λ} S + ) j Sz j Sz j+ {µ} {µ} {µ} {λ} S+ S + S S {λ} = {λ} j S+ j S+ j+ {µ} L N ) L N ) {µ} {µ} {λ} {λ} A.6) We can now check numerically for solutions of the Bethe equation with same momenta using similar determinant expression as in cf. Appendix B) for these matrix elements, that due to the normalization factor the matrix element is suppressed for large L at finite magnetic field.

22 Spinon Decay Appendix B. Calculation of the next-nearest neighbor spin operator matrix element We want to calculate the matrix element {λ} σ z j σ z j+ {µ} B.) with {µ}, {µ} Bethe states and {µ}, {λ} satisfying the Bethe equations ). We do the calculation for all. To obtain the formula for = the general functions a and d have to be replaced for the ones mentioned above and γ has to be set to corresponding to the rescaling of {λ} with γ and taking the limit γ 0. The σ operators are given in terms of the Bethe-Ansatz operators A, B, C, D, as obtained in []: j σj z = A + D)ξ i )Dξ j ) i= L k=j+ A + D)ξ k ) + B.) where ξ i is an inhomogeneity parameter, introduced at every site in the chain for technical reasons. We will set ξ i iγ/ in the end, but will keep them for the calculation. We note that now dλ) is defined as: L dλ) = bλ ξ l ) l= With this we can write the matrix element as: B.3) {λ} σj z σj+ z {µ} = {λ} σj z {µ} + {λ} σj+ z {µ} {λ} {µ} j L + 4 {λ} A + D)ξ i )Dξ j )A + D)ξ j+ )Dξ j+ ) A + D)ξ k ) {µ} i= k=j+3 B.4) The maxtrix elements and the overlap in the first line are known [, 0]. However as we are interested in j σz j σz j+ and as {λ} and {µ} are orthogonal and eigenstates of σ z, we only need to calculate the expression in the second line. From the Yang-Baxter algebra one can derive the commutation relations between the operators A, B, C, D and from this one gets []): Aµ) Bλ j ) 0 = aµ) b λ j µ) Bλ j ) 0 j= j= N n= j= aλ n ) cλ n µ) bλ n µ) j= j= b λ j λ n )Bµ) Bλ j ) 0 j n Dµ) Bλ j ) 0 = dµ) b µ λ j ) Bλ j ) 0 N n= dλ n ) cµ λ n) bµ λ n ) j= j n j n b λ n λ j )Bµ) Bλ j ) 0 j n B.5) B.6)

23 Spinon Decay 3 We furthermore know that j A + D)ξ i ) {λ} = exp ) ijp {λ} {λ} B.7) i= where {λ} is a Bethe state and P {λ} is the total momentum of {λ} and L A + D)ξ i ) {λ} =. B.8) i= We now need to calculate the matrix element: D = {λ} Dξ j )A + D)ξ j+ )Dξ j+ ) {µ} B.9) Using the commutation relations for A, B, C, D cf. []) we obtain: D = I + II B.0) where I is dependent on simple matrix elements where two rapidities are replaced with inhomogeneities and II is dependent on a matrix element with three insertions of inhomogeneities: I = N a= dµ a )cξ j+ µ a ) bξ j+ µ a ) N b a i a dµ b )cξ j µ b ) bξ j µ b )bµ b ξ j+ ) + cξ j+ ξ j+ ) bξ j+ ξ j+ ) i a,b i a { bµ a µ i ) bξ j+ ξ j+ ) i a,b bµ i ξ j+ ) i a bµ i ξ j+ ) bµ b µ i ) {λ} Bξ j)bξ j+ ) N dµ b )cξ j µ b ) bξ j µ b ) bµ b ξ j+ ) } Bµ j ) 0 b a bµ b µ i ) {λ} Bξ j)bξ j+ ) j a,b N dµ a )cξ j+ µ a ) II = bξ a= j+ µ a ) bµ a µ i ) i a { N [ cµ b ξ j+ ) bµ b ξ j+ ) bµ i µ b ) bξ j+ µ b ) b a N c a,b i a,b cξ j+ µ b ) + dµ b ) bξ j+ µ b )bµ b ξ j+ ) dµ c )cξ j µ c ) bξ j µ c )bµ c ξ j+ )bµ c ξ j+ ) {λ} Bξ j )Bξ j+ )Bξ j+ ) i a,b,c Bµ i ) 0 i a,b,c } i a,b j a,b ] bµ b µ i ) bµ c µ i ) Bµ j ) 0 B.)

24 Spinon Decay 4 Now we can use Slavnov s formula [0] for the overlap of two states. One of the states has to be a Bethe state, the other state can be parametrized by an arbitrary set of rapidities. Let {λ} be solutions of the Bethe equations ) and {µ} arbitrary, then one gets: det H{λ}, {µ})) {λ} {µ} = j>k φλ k λ j ) j<k φµ k µ j ) where H is a matrix defined as φiγ) H ab = aµ b ) φλ k µ b + iγ) dµ b ) φλ k µ b iγ) φλ a µ b ) k a k a with φx) = x and γ set to in the = scaling limit. We will now treat I and II seperately. Appendix B.. Part II B.) B.3) For the part II, the limit of the ξ j, ξ j+ and ξ j+ can be taken seperately for the matrix element and the prefactor. Taking the limit ξ i iγ/ for the prefactor amounts to replacing the ξ i with iγ/. For the matrix element we obtain using Slavnov s determinant formula: {λ} Bξ j )Bξ j+ )Bξ j+ ) Bµ i ) 0 = i k,m,n det H{λ}, {µ i k,m,n, ξ j, ξ j+, ξ j+ }) a>b φλ a λ b ) a<b φµ b µ a )) µk ξ j,µ m ξ j+,µ n ξ j+ B.4) where φµ b µ a ) µk ξ j,µ m ξ j+,µ n ξ j+ = φξ j ξ j+ )φξ j ξ j+ )φξ j+ ξ j+ ) φµ k µ m )φµ k µ n )φµ m µ n ) a<b c k,m,n φξ j µ c )φξ j+ µ c )φξ j+ µ c ) φµ k µ c )φµ m µ c )φµ n µ c ) φµ b µ a ) a<b B.5) and in the determinant we have to replace µ k ξ j, µ m ξ j+ and µ n ξ j+. Therefore the important part when taking the limits is: lim ξ j,ξ j+ iγ/ det H{λ}, {µ i k,m,n, ξ j, iγ/, ξ j+ }) φξ j iγ/)φξ j ξ j+ )φiγ/ ξ j+ ) B.6) where the limit ξ j+ iγ/ is already taken. Let us now take the limit ξ j+ iγ/. We see that both numerator and denominator go to zero here. Therefore using the rule of l Hospital we obtain: lim ξ j,ξ j+ iγ/ lim ξ j iγ/ det H{λ}, {µ i k,m,n, ξ j, iγ/, ξ j+ }) φξ j iγ/)φξ j ξ j+ )φiγ/ ξ j+ ) = φξ j iγ/)φξ j iγ/) lim x iγ/ d dx det H{λ}, {µ i k,m,n, ξ j, iγ/, x})) B.7)

25 Spinon Decay 5 Analogous to [34] we can now use a Laplace expansion of the determinant for the column that is dependent on x and evaluate the derivative and limit: lim x iγ/ d dx det H{λ}, {µ i k,m,n, ξ j, iγ/, x})) = d = lim ) n+i) C i x) minor ni H{λ}, {µ i k,m,n, ξ j, iγ/, x})) x iγ/ dx = lim x iγ/ i i ) n+i) d dx C ix) ) minor ni H{λ}, {µ i k,m,n, ξ j, iγ/, x})) B.8) where the minor is not dependent on x and C i x) is the ith element of the column n. Therefore we get: lim x iγ/ d dx det H{λ}, {µ i k,m,n, ξ j, iγ/, x})) = det H {λ}, {µ}; m, n) ) µk ξj B.9) with H ) ab {λ}, {µ}; m, n) = H ab φiγ) φλ a iγ )φλa+ iγ ) φiγ)φλ a) φλ a iγ ) φλ a+ iγ ) b m, n b = m b = n B.0) Repeating this step for the limit ξ j iγ/ using the rule of l Hospital twice we finally obtain: M ) ijk lim {λ} Bξ j)bξ j+ )Bξ j+ ) Bµ i ) 0 ξ j,ξ j+,ξ j+ iγ/ where = i i k,m,n φ λ i + iγ ) 3 φµ m µ n )φµ k µ m )φµ k µ n ) c k,m,n φµ c µ k )φµ c µ m )φµ c µ n ) φµ c iγ )3 det H ) a>b φλ b λ a ) a<b φµ b µ a ) H ) ab {λ}, {µ}; k, m, n) = H ab φiγ) φλ a iγ )φλa+ iγ ) φiγ)φλ a) φλ a iγ ) φλ a+ iγ ) φiγ) x φλ a x)φλ a x+iγ)) x= iγ b m, n, k b = m b = n b = k B.) B.)

26 Spinon Decay 6 With this we obtain after some algebra: II = k φλ k + iγ/) 3 φµ k iγ/) 3 N a>b φλ b λ a ) a<b φµ dµ n )φµ n iγ/) φµ i µ n iγ) b µ a ) n= i [ φµm 3iγ/) φµ i µ m + iγ)φµ m iγ/) φµ m µ n iγ) φµ ] m + iγ/) φµ m µ m + iγ) m n i k m,n dµ k ) i φµ k + iγ/) φµ i µ k iγ) φµ m µ k iγ)φµ n µ k iγ) det H {λ}, {µ}; k, m, n)) B.3) Using a Lemma from Laplace s determinant formula [3]) we finally obtain: II = k φλ k + iγ/) 3 φµ k iγ/) 3 a>b φλ b λ a ) a<b φµ b µ a ) N A nm [det G nm + B nm ) det G nm )] n= m n B.4) with G nm ) ab = H ab φiγ)φλ a) φλ a iγ ) φλ a+ iγ ) φiγ) x φλ a x)φλ a x+iγ)) x= iγ b m, n b = m b = n B.5) B nm ) ab = δ bm ) δ bn )dµ b ) φµ i µ b iγ)φµ b + iγ/) φiγ) φλ a iγ/)φλ a + iγ/) B.6) i m,n Appendix B.. Part I First we can take the limit ξ j, ξ j+ iγ/. This again just amounts to replacing the ξ j and ξ j+ with iγ/ in the prefactors and doing a similar analysis for the matrix element depending on both ξ j and ξ j+ as for the part II. However performing the limit ξ j+ iγ/ is a bit more involved, as the limit can not be taken independently for prefactor and matrix elements. This is due to terms ξ j+ iγ/ appearing in the prefactors. Doing a consistent series expansion in ξ j+ iγ/ again utilizing the Laplace determinant expansion and then taking the limit ξ j+ iγ/ we obtain after

27 Spinon Decay 7 some calculation: I = φiγ) { N dµ n ) n= i k φλ k + iγ/) φµ k + iγ/)φµ k iγ/) 3 a>b φλ b λ a ) b>a φµ b µ a ) φµ i µ n iγ) φµ n iγ/) φµ n + iγ/) dµ m ) φµ i µ m iγ)φµ m + iγ/) m n i n φ iγ) φiγ) φiγ) φµ d iγ/)φµ d + iγ/) + d n φ λ b + iγ/) φλ b + iγ/) + d m,n φiγ) φµ m + iγ/)φµ m iγ/) φ µ d iγ/) φµ d iγ/) ) det H 3 {λ}, {µ}; m, n)) + b dµ m )φµ m + iγ/) φµ i µ m iγ) det H 4{λ}, {µ}; m, n)) m n i n B.7) with H 3 ) ab = H 4 ) ab = H ab φiγ) φλ a+iγ/)φλ a iγ/) φiγ)φλ a) φλ a iγ ) φλ a+ iγ ) H ab φiγ) φλ a+iγ/)φλ a iγ/) φiγ) x b m, n b = m b = n φλ a x)φλ a x+iγ)) x= iγ b m, n b = m b = n B.8) B.9) using the Lemma from Laplace s determinant formula again, we finally obtain: I = φiγ) n k φλ k + iγ/) φµ k + iγ/)φµ k iγ/) 3 a>b φλ b λ a ) b>a φµ b µ a ) A n det G ) n + B n ) ) + det G ) n + B n ) det G )) det G )) ) ) B.30)

28 Spinon Decay 8 where G ) n ) ab = G ) n ) ab = { Hab b n φiγ)φλ a) b = n φλ a iγ ) φλ a+ iγ { ) Hab b n x x= φiγ) φλ a x)φλ a x+iγ)) iγ b = n B ) n ) ab = δ bn )dµ b ) i n φµ i µ b + iγ)φµ b + iγ/) φ iγ) φiγ) φiγ) φµ d iγ/)φµ d + iγ/) + d n b φ λ b + iγ/) φλ b + iγ/) + B n ) ) ab = δ bn ) d µ b d m,n φiγ) φµ m + iγ/)φµ m iγ/) φ µ d iγ/) φµ d iγ/) ) φiγ) φµ i µ b iγ)φµ b + iγ) φλ a + iγ/)φλ a iγ/) i n B.3) B.3) B.33) B.34) Appendix B.3. Total matrix element We can now put together the total matrix element. We obtain {λ} Sj z Sj+ z {µ} = e ip {λ}+ip {µ} e ijp {λ} P {µ} ) D j j B.35) = Le ip {λ}+ip {µ} δ P{λ},P {µ} D B.36) [] L. D. Faddeev, L. A. Takhtajan, What is the spin of a spin wave?, Lett. A ) [] L.D. Faddeev and L. Takhtajan, Spectrum and scattering of excitations in the one-dimensional isotropic Heisenberg model, J. Sov. Math. 4, 4 984). [3] G. Delfino, G. Mussardo, P. Simonetti, Non-integrable quantum field theories as perturbations of certain integrable models, Nucl. Phys. B473, ). [4] G. Delfino, G. Mussardo, Non-integrable aspects of the multi-frequency sine-gordon model, Nucl. Phys. B56, ). [5] G. Delfino, P. Grinza, G. Mussardo, Decay of particles above threshold in the Ising field theory with magnetic field, Nucl. Phys. B737, 9 006). [6] B. Pozsgay and G. Takacs, Characterization of resonances using finite size effects, Nucl. Phys. B 748, ). [7] A. Imambekov, T. L. Schmidt, and L. I. Glazman, One-dimensional quantum liquids: Beyond the Luttinger liquid paradigm, Rev. Mod. Phys. 84, 53 0); R. G. Pereira, S. R. White, and I. Affleck, Spectral function of spinless fermions on a one-dimensional lattice, Phys. Rev. B79, ) and references therein. [8] F.H.L. Essler, R.G. Pereira and I. Schneider, Spin-charge-separated quasiparticles in onedimensional quantum fluids, Phys. Rev. B9, ). [9] T.L. Schmidt, A. Imambekov, and L. I. Glazman, Spin-charge separation in one-dimensional fermion systems beyond Luttinger liquid theory, Phys. Rev. B ); F.H.L. Essler, Threshold singularities in the one-dimensional Hubbard model, Phys. Rev. B 8, ). [0] F. H. L. Essler, H. Frahm, F. Göhmann, A. Klümper, and V. E. Korepin, The One-Dimensional Hubbard Model, Cambridge University Press, Cambridge 005). [] M. Takahashi, Thermodynamics of One-Dimensional Solvable Models, Cambridge University Press 999)

29 Spinon Decay 9 [] V. Korepin, N. Bogoliubov, and A. Izergin, Quantum Inverse Scattering Method and Correlation Functions, Cambridge Monographs on Mathematical Physics, Cambridge University Press 997) [3] M. Peskin, D. Schroeder, An introduction to quantum field theory, Westview Press Reading Mass.) 995) [4] H. Bethe, Zur Theorie der Metalle. I. Eigenwerte und Eigenfunktionen der linearen Atomkette, Zeitschrift für Physik, ) [5] M. Gaudin, J.S. Caux, The Bethe Wavefunction, Cambridge University Press 04) [6] R. Orbach, Linear Antiferromagnetic Chain with Anisotropic Coupling, Phys. Rev., ) [7] G. Mussardo, Statistical field theory an introduction to exactly solved models in statistical physics, Oxford University Press 00) [8] M. Gaudin, B. M. McCoy, T. T. Wu, Normalization sum for the Bethe s hypothesis wave functions of the Heisenberg-Ising chain, Phys. Rev. D 3, 47 98) [9] V. E. Korepin, Calculation of norms of Bethe wave functions, Comm. Math. Phys. 86, 39 98) [0] N. A. Slavnov, Calculation of scalar products of wave functions and form factors in the framework of the algebraic Bethe ansatz, Teor. Mat. Fiz. 79, 3 989) [] N. Kitanine, J. M. Maillet, V. Terras, Form factors of the XXZ Heisenberg spin-/ finite chain, Nucl. Phys. B 554, ) [] N. Kitanine, J. M. Maillet, V. Terras, Correlation functions of the XXZ Heisenberg spin-/ chain in a magnetic field, Nucl. Phys. B 567, ) [3] A. Klauser, J. Mossel, J. S. Caux, Adjacent spin operator dynamical structure factor of the S = / Heisenberg chain, J. Stat. Mech. P030 0) [4] J. S. Caux, Correlation functions of integrable models: A description of the ABACUS algorithm, Journal of Mathematical Physics 50, ) [5] R. G. Pereira, J. Sirker, J. S. Caux, R. Hagemans, J. M. Maillet, S. R. White, I. Affleck, Dynamical structure factor at small q for the XXZ spin-/ chain, J. Stat. Mech. P ) [6] J. S. Caux, R. Hagemans, J. M. Maillet, Computation of dynamical correlation functions of Heisenberg chains: the gapless anisotropic regime, J. Stat. Mech. P ) [7] A. Shashi, L. I. Glazman, J. S. Caux, A. Imambekov, Nonuniversal prefactors in the correlation functions of one-dimensional quantum liquids, Phys. Rev. B 84, ) [8] A. Shashi, M. Panfil, J. S. Caux, A. Imambekov, Exact prefactors in static and dynamic correlation functions of one-dimensional quantum integrable models: Applications to the Calogero-Sutherland, Lieb-Liniger, and XXZ models, Phys. Rev. B 85, ) [9] N. Kitanine, K. K. Kozlowski, J. M. Maillet, N. A. Slavnov, V. Terras, On the thermodynamic limit of form factors in the massless XXZ Heisenberg chain, J. Math. Phys ) [30] N. Kitanine, K. K. Kozlowski, J. M. Maillet, N. A. Slavnov, V. Terras, Thermodynamic limit of particle-hole form factors in the masless XXZ Heisenberg chain, J. Stat. Mech. 05 P0508 0) [3] N. Kitanine, K. K. Kozlowski, J. M. Maillet, N. A. Slavnov, V. Terras, Form factor approach to dynamical correlation functions in critical models, J. Stat. Mech., P0900 0) [3] N. Kitanine, K. K. Kozlowski, J. M. Maillet, N. A. Slavnov, V. Terras, A form factor approach to the asymptotic behavior of correlation functions in critical models, J. Stat. Mech. : Th. and Exp., P00 0) [33] R. G. Pereira, S. R. White, I. Affleck, Spectral function of spinless fermions on a one-dimensional lattice, Phys. Rev. B ) [34] R. L. Hagemans, Dynamics of Heisenberg Spin Chains, PhD thesis 007), retrieved from UvA- DARE database, uvapub:5466