Long-time large-distance asymptotics of the transverse correlation functions of the XX chain in the spacelike regime

We derive an explicit expression for the leading term in the long-time, large-distance asymptotic expansion of a transverse dynamical two-point function of the XX chain in the spacelike regime. This expression is valid for all nonzero finite temperatures and for all magnetic fields below the saturation threshold. It is obtained here by means of a straightforward term-by-term analysis of a thermal form factor series, derived in previous work, and demonstrates the usefulness of the latter.


Introduction
The XX chain is a spin chain with Hamiltonian [15] where the σ α j , α = x, y, z, are Pauli matrices acting on site j ∈ {1, . . ., L} of an L-site periodic lattice, σ α 0 = σ α L .The parameters J > 0 and h > 0 denote the strengths of the spin-spin interaction and of the applied magnetic field.We shall restrict the magnetic field to values below the saturation threshold, 0 < h < 4J.
In our recent work [9] we have derived a novel form factor series for the transverse dynamical correlation function of the XX chain in equilibrium with a heat bath at temperature T .It measures the space-time evolution of a local perturbation relating two points at distance m and temporal separation t.
Our series originates from a form factor expansion related to the quantum transfer matrix [7].
It can be resummed into a 'Fredholm determinant representation' consisting of a prefactor times a Fredholm determinant of an integrable integral operator [11].The latter is different from the Fredholm determinant representation derived by Colomo et al. in [5].
For Fredholm determinants and resolvent kernels of integrable integral operators a general method [6] is available that allows one to analyse their asymptotic dependence on parameters.Starting with the Fredholm determinant representation obtained in [5] the authors of [12] applied this 'nonlinear steepest-descent method' to the late-time, largedistance analysis of (2) at a fixed ratio α = m/(4Jt).They found exponential decay of the form where C, ν and ξ depend on T , h and α.The functional dependence differs according to whether α > 1 or α < 1.The former regime, in which the spatial distance in units of 4J is larger than the temporal separation, is called 'spacelike', while the latter is referred to as 'the timelike regime'.
In [12] the authors considered magnetic fields below the saturation threshold, 0 < h < 4J.They obtained explicit expressions for ν and ξ in both, space-and timelike regimes.Later the 'constant term' C was obtained for h > 4J in [13].Although the nonlinear steepest descent method would allow one to calculate C for 0 < h < 4J as well, it seems that nobody has ever attempted to do so.This may be partially attributed to the cumbersome nature of the required calculations.
In this work we reconsider the late-time, large-distance asymptotic analysis of the two-point function σ − 1 σ + m+1 (t) T in the spacelike regime.It turns out that the novel thermal form factor series derived in [9] allows us to obtain the asymptotics, including the constant term C, by a rather elementary term-by-term analysis of the series that avoids the use of any Riemann-Hilbert problem.
On the other hand, our thermal form factor series can be resummed into a Fredholm determinant representation as well.As we shall see below this Fredholm determinant representation is rather different from the one of Its et al. [12] in that the term that provides the leading late-time, large-distance asymptotics in the spacelike regime appears to be pulled out as a pre-factor.Our finding strikingly resembles in structure the Borodin-Okounkov, Geronimo-Case formula [2,3,8] for a Toeplitz determinant generated by a symbol satisfying the hypotheses of the Szegö theorem.
We should point out that the late-time, large-distance asymptotics considered here do not commute with the low and high-temperature asymptotics.At any finite temperature the asymptotic decay of the transverse two-point functions is exponential and given by (3).If, however, the temperature is send to zero first, the correlation functions will vary algebraically [14].We shall consider this limit for the more general XXZ chain in subsequent work.If we send the temperature to infinity first, then the behaviour of the correlation functions in 'time-direction' becomes Gaussian [4,16].We have recently analysed the latter situation in full generality in [11], which is one of two companion papers of this work.In the other one [10] we evaluate the two-point function numerically, for a wide range of temperature and space-time separations, directly from the novel Fredholm determinant representation.

Thermal form factor series representation
The starting point of our analysis will be a thermal form factor series for the transversal two-point function (2) derived in [9].The series is a series of multiple integrals which is most compactly expressed in terms of certain functions characteristic of the XX chain.These are in first place the momentum p and the energy of the single-particle excitations of the Hamiltonian expressed in terms of the rapidity variable, Here we choose the principal branch of the logarithm in the definition of the momentum function p(λ), cutting the complex plane from −iπ/2 to zero modulo iπ.Because of the πi-periodicity of the momentum, shared by all other functions in our form factor series, we may think of these functions as being defined on a cylinder of circumference π, which is equivalent to restricting their values to the 'fundamental strip' It is easy to see that has precisely two roots in S.These roots are called the Fermi rapidities.The value of the momentum function evaluated at the left Fermi rapidity is the Fermi momentum.
Using the Fermi rapidities we can represent the energy function as Energy and momentum functions and p are real on the lines x ± iπ/4, x ∈ R, where they take the values Figure 1: Sketch of the hole and particle contours C h and C p .The Fermi rapidity The one-particle energy determines the function Most of the functions occurring in the form-factor series below are defined as integrals over two simple closed contours C h and C p , involving p, , z and some hyperbolic functions.The 'hole contour' C h and the 'particle contour' C p are sketched in Fig. 1.They are defined in such a way that C h encloses all roots of e − (x)/T − 1 located inside the strip −π/4 < Im x < π/4 ('the holes') as well as the left Fermi rapidity λ − F , whereas C p encloses the roots of e − (x)/T − 1 inside the strip π/4 < Im x < 3π/4 ('the particles') as well as the right Fermi rapidity λ + F .Given these contours we define the Cauchy transforms for all x ∈ S \ C h , and for all x ∈ S \ C p .For fixed . Since the integrands in ( 11), (12) are rapidly decaying for λ → ±∞, we may deform the contours and conclude that Another function needed below is the square of a generalized Cauchy determinant, After these preparations we can now recall the form factor series derived in [9].Using the above notation and performing several more or less obvious simplifications it can be written as where The contour C h is tightly enclosed by C h .
3 Asymptotics in the spacelike regime Theorem.In the spacelike regime m > 4Jt the form factor series ( 15) is absolutely convergent and determines the late-time, large-distance asymptotics of the transverse dynamical correlation function of the XX chain as where In preparation of the proof we introduce the short-hand notations and the function with real and imaginary parts u(λ) = Re g(λ) and v(λ) = Im g(λ).Then the 'wave factors' in (15) take the form e ±i(mp(λ)−t (λ)) = e ∓iht±τ g(λ) .
We will be interested in the asymptotic behaviour of (15) for large positive τ and fixed α > 1.As we shall see below, it is determined by the poles of the integrands at λ ± F .The saddle points contribute only to the subleading corrections.This becomes clear when we consider the function g close to the lines R ± iπ/4 and on the lines Re λ = ±R for R > 0 large enough.
(ii) Define the oriented contours where R, δ > 0. Then R and δ can be chosen in such a way that u(λ) < 0 for all λ ∈ C h,sd , u(λ) > 0 for all λ ∈ C p,sd and all hole roots are inside C h,sd , while all particle roots are inside C p,sd .
Since α > 1, there is a unique ϕ > 0 such that α = cth(2ϕ).Using this parameterization we find for any λ = x + iy ∈ S that Here the first term in the square bracket is unbounded from above for x → ±∞, implying that the only roots of ∂ y u(λ) in S are at y = 0, π/2 if |x| is large enough.Taking into account (26) we see that, if the latter is the case, then It follows that u(λ) < 0 of the lines ±R + i(−π/4, π/4), while u(λ) > 0 on ±R + i(π/4, 3π/4), if R > 0 large enough.The statement about the location of the particle and hole roots follows by straightforward inspection of the integrands in (15).
Proof of the Theorem.The function D {x j } n j=1 , {y k } n−1 k=1 is symmetric separately in all x j and y k .It satisfies if x j = x k or y j = y k for all j = k.Setting V h (x) = Φ p (x)e i(mp(x)−t (x)) πi e (x)/T − 1 , V p (x) = e −i(mp(y)−t (y)) πiΦ h (y) 1 − e − (y)/T (31) and using the above lemma we therefore obtain Here C h,sd and C p,sd are the contours introduced in ( 22).Notice that we consider res dx V h (x), x = λ − F as a functional acting on functions f holomorphic in a disc D ε (λ − F ) of sufficiently small radius ε centered about λ − F as and similarly for res dy V p (y), y = λ + F .In particular, Equation ( 32) implies that where the four series S (m, t) can be written as follows. with In order to show the convergence of the series and to estimate their asymptotic behaviour, we have to establish bounds on the individual terms.We start with the functions D {x j } n j=1 , {y k } n k=1 and recall the Hadamard bound for the determinant of an n × n matrix det j,k=1,...,n Since the contours C h,sd and C p,sd are finite and disjoint, we can use (38) to estimate where Likewise we set As follows from the above lemma, there exist κ, c > 0 such that With this we obtain a bound on every individual term in the series S 1 , for some constant C 1 > 0. This implies absolute convergence of the series S 1 and shows that, asymptotically for large τ , the series behaves like In a similar way one obtains for constants C j > 0, j = 2, 3, 4. It follows that the series S j , j = 2, 3, 4, converge absolutely and behave asymptotically as Inserting ( 44), ( 46) into (35) and recalling the explicit form (16) of F(m) we have arrived at the statement of the theorem.
The theorem fixes the constant term of the asymptotics in the spacelike regime that remained undetermined in [12].Note that the function (λ − F ) can be easily calculated explicitly, For the other factors composing the constant C(T, h) we did not find any further simplification so far.

Discussion
For the interpretation of our result we would like to recall a Fredholm determinant representation of the transversal two-point function (2) that was obtained in [11], where it was used for the asymptotic analysis of the correlation function in the high-temperature limit.
Referring to [11] we define the functions ϕ(x, y) = e y−x sh(y − x) and Using these functions we define two integral operators V and P acting on functions on the contour C p , Then (cf.[11]) the transversal correlation functions of the XX chain admit the Fredholm determinant representation Comparing with the asymptotic behaviour of the correlation function in the spacelike regime m > 4Jt we see that meaning that the Fredholm determinant collects the higher-order corrections to the main asymptotics.This is the analogy with the Borodin-Okounkov-Geronimo-Case formula [3,8] mentioned in the introduction.
On the level of the Fredholm determinant representation it is easiest to compare our result with that of Its, Izergin, Korepin and Slavnov [12].For this purpose we rewrite their integral operators acting on functions on the the unit circle as integral operators acting on functions on . This is achieved by employing the map z → e ip(λ) to the Fredholm determinant representation in [12].Then where W is an integrable operator with kernel and Q is a one-dimensional projector acting as Comparing ( 51) and (53) we see that in (53) the late-time, large-distance asymptotics is entirely inside the Fredholm determinants and therefore harder to analyse.The fact that the late-time, large-distance asymptotic behaviour of the transverse dynamical correlation functions of the XX chain, including the constant term, can be obtained directly from the series representation (15) raises a number of interesting questions.
(i) Is a similar analysis possible for the XXZ quantum spin chain?Unlike the XX chain treated in this work no Fredholm determinant representation for its two-point function is expected to exist, but a thermal form-factor series similar to (15) is still available [9].As the structure of the saddle-point equations is very similar, there seems to be a good chance that the answer will turn out to be positive.
(ii) What can be done in the timelike regime?Here all terms in the series (15) contribute to the late-time, large-distance asymptotics.A further resummation would be necessary.Can we devise a method to find such a generalization?We would like to close with two remarks.First, in our recent work [10] we have compared the asymptotic formula of our theorem with a numerical evaluation based on the Fredholm determinant representation (51).As should be clear from the fact that the corrections are exponentially small for large m and t the asymptotic formula turns out to be very efficient.For an example see Fig. 2. Second, the constant term C(t, h), equation (18), does not depend on α.For this reason it should agree with the constant obtained by Barouch and McCoy [1] in form of infinite double products in their analysis of the static correlation functions (see equations (3.17)-(3.19) of their paper).We have numerical evidence that this is indeed the case.