Correlations of zero-entropy critical states in the XXZ model:
integrability and Luttinger theory far from the ground state
R. Vlijm, I.S. Eliëns, J.-S. Caux
Institute for Theoretical Physics, University of Amsterdam,
Science Park 904, 1098 XH Amsterdam, The Netherlands
July 8, 2019
Pumping a finite energy density into a quantum system typically leads to ‘melted’ states characterized by exponentially-decaying correlations, as is the case for finite-temperature equilibrium situations. An important exception to this rule are states which, while being at high energy, maintain a low entropy. Such states can interestingly still display features of quantum criticality, especially in one dimension. Here, we consider high-energy states in anisotropic Heisenberg quantum spin chains obtained by splitting the ground state’s magnon Fermi sea into separate pieces. Using methods based on integrability, we provide a detailed study of static and dynamical spin-spin correlations. These carry distinctive signatures of the Fermi sea splittings, which would be observable in eventual experimental realizations. Going further, we employ a multi-component Tomonaga-Luttinger model in order to predict the asymptotics of static correlations. For this effective field theory, we fix all universal exponents from energetics, and all non-universal correlation prefactors using finite-size scaling of matrix elements. The correlations obtained directly from integrability and those emerging from the Luttinger field theory description are shown to be in extremely good correspondence, as expected, for the large distance asymptotics, but surprisingly also for the short distance behavior. Finally, we discuss the description of dynamical correlations from a mobile impurity model, and clarify the relation of the effective field theory parameters to the Bethe Ansatz solution.
The ground states and low-lying excitations of one-dimensional many-body quantum systems often display interesting features associated to the Tomonaga-Luttinger liquid universality class [1, 2], notable examples being one-dimensional quantum gases and spin chains. Bosonization techniques [3, 4] then provide a route for the description of the long-range asymptotics of basic correlations. One-dimensional quantum physics is also famous for possessing a class of integrable theories for which all eigenstates can be exactly obtained by Bethe Ansatz . Prototypical examples of Bethe Ansatz solvable models are the Lieb-Liniger model [6, 7] of delta-interacting bosons and the Heisenberg spin chain [8, 9]. For these and other integrable models, the algebraic Bethe Ansatz  and more specifically Slavnov’s theorem [11, 12] provides a method to compute matrix elements of local operators between two Bethe states. Summations over relevant matrix elements of particle-hole like excitations can then be implemented, leading to efficient quantitative evaluations of dynamical correlation functions and expectation values of local operators.
The correspondence between results from integrability on the one hand, and Tomonaga-Luttinger theory on the other, are well-known for ground states of the above-mentioned models. The bosonized Tomonaga-Luttinger liquid descriptions are characterized by universal Luttinger parameters related to the compressibility and sound velocity, which can be fitted from the energy levels of the model at finite system size computed by alternative computational methods . For example, the compressibility is readily fitted by the energy levels calculated from integrability upon addition or removal of particles from the Fermi sea. Moreover, analysis of the system size scaling of matrix elements of Umklapp states [14, 15, 16, 17] calculated by algebraic Bethe Ansatz methods allows for the determination of non-universal exponents and prefactors appearing in correlation functions of the bosonized theory.
It is interesting to ask whether these correspondences between exact results and field-theory predictions are strictly limited to the vicinity of the ground state, or whether other regions of Hilbert space can similarly be ‘captured’ both by integrability and by an appropriate field theory. From integrability, things look promising: since the Bethe Ansatz does not care whether the wave functions one writes are near or far from the ground state (these remain exact, irrespective of which energy they have), correlations on such high-energy states can in principle be obtained by extensions of existing ground state-based summation methods. On the other hand, the applicability of Tomonaga-Luttinger theory relies on the linearity of the dispersion relations of the effective excitations around the given state, which cannot be generic for high-energy states. Here, we concentrate on a class of states for which this bosonization procedure can be applied. The setup is easy to visualize: starting from the ground state’s Fermi sea configuration, we give a finite opposite momentum to two groups of particles, in fact splitting the Fermi sea in two seas , yielding what we coin a ‘Moses’ state which, while having a macroscopic energy above the ground state, still possesses a zero entropy density (and thus a potential for displaying critical power-law behavior in its correlations. The corresponding multi-component Tomonaga-Luttinger liquid theory and its correlations were recently obtained for the Lieb-Liniger model [19, 20]. The entanglement entropy of these states is nontrivial  and consistent with predictions from the effective conformally invariant field theory.
The physical motivation to study such states with a double Fermi sea with opposite momentum kicks originates from experiments where a Bragg pulse is applied to a gas of interacting bosons to create an initial state with counterpropagating particles, leading to a quantum realization of the famous Newton’s cradle . It has recently been shown that late-time correlations in the quantum Newton’s cradle experiment are better described by a finite entropy version of these states, which can be constructed in the Lieb-Liniger model by theoretically mimicking the application of a Bragg pulse on the ground state and for which the peaks are smoothened and show the characteristic ghost-shaped momentum distribution function at late times . The zero-entropy nature of the states we consider here translates into a quasi-condensate momentum distribution with sharp peaks at finite momenta similar to the ones observed in cold atoms after domain-wall melting of a one-dimensional Mott insulator . At present it is not clear whether a connection exists between the steady state for this quantum quench and the split-Fermi-sea states we consider. Finally, another motivation is related to spin ladders, where states with a Fermi sea consisting of distinct pockets can appear as the ground state [25, 26]. The comparison of their work with ours puts into sharp focus which characteristics should be attributed to the out-of-equilibrium nature of the state and which should not.
The aim of the present article is to extend the approach for the Bose gas with double Fermi seas elaborated in Ref.  to the anisotropic Heisenberg spin chain. We compute the dynamical structure factor for this class of excited zero-entropy states from integrability. The static spin-spin correlations are subsequently studied from both the matrix element summation approach from algebraic Bethe Ansatz and the multi-component Tomonaga-Luttinger predictions supported by parameter fitting from integrability, which show surprisingly good correspondence for both the long-range asymptotics and short distances. Going further, we also extend recent advances in the computation of time-dependent correlations by means of mobile impurity models and apply these extensions here. The correspondence with the effective field theory becomes more delicate in this case, but if the separation between seas is large enough the method still gives adequate results.
This article is structured in the following way. Section 2 describes the setup from Bethe Ansatz of the zero-entropy critical states consisting of a state with two Fermi seas, while the dynamical structure factor of these states are evaluated using algebraic Bethe Ansatz matrix element summations in section 3. Section 4 gives the multi-component Tomonaga-Luttinger liquid approach to the real space correlations, which are compared to the algebraic Bethe Ansatz results in section 5. Finally, a description of dynamical correlations by a mobile impurity model is given in section 6.
2 Zero-entropy critical states in the XXZ model
with periodic boundary conditions . We fix such that the ground state in zero field is antiferromagnetic, and furthermore restrict our analysis to the quantum critical cases with and . The Tomonaga-Luttinger theory is applicable to these regimes, allowing for a comparison of both approaches.
The commutation of the total spin operator along the -axis with Hamiltonian (1) assures a splitting of the Hilbert space in sectors of fixed magnetization . In such a fixed- subsector, the Bethe Ansatz wave functions [5, 9] are constructed from the fully polarized vacuum state as plane waves of magnons
for which the amplitudes are set by the scattering phases. Each eigenstate of Hamiltonian (1) is specified by a unique, non-coinciding set of rapidities satisfying Bethe equations, which are derived from the scattering phases upon the interchange of two magnons and the imposition of periodic boundary conditions. In logarithmic form, the Bethe equations are given by
where for and , for . The momentum in Eq. (2) is expressible as .
The Bethe quantum numbers are (half-odd) integers for (even) odd, and form the starting point of the construction of Bethe states at finite system size. The set of rapidities can be obtained by solving the logarithmic Bethe equations numerically by an iterative procedure for a given set of non-coinciding quantum numbers . With the set of rapidities of a Bethe state at hand, properties of the state such as its energy can be calculated directly,
with for and for , while the momentum of a Bethe state is expressed directly in terms of its quantum numbers as
In general, rapidities can take on complex values, while the full set solving Bethe equations must remain self-conjugate, leading to an arrangement of the rapidities in terms of string solutions. String solutions will not be considered for the initial state consisting of two separate Fermi seas of real rapidities. However, string solutions can occur in the intermediate states of the matrix element summations described in section 3.
The sets of quantum numbers specifying the Bethe states are bounded by limiting quantum numbers , which are derived by taking the limit of one of the rapidities to infinity and computing the associated quantum number from the Bethe equations (3). The maximum allowed quantum number to still give a finite valued rapidity is then . For a Bethe state consisting of finite real rapidities, one obtains and , meaning there are available possibilities to distribute quantum numbers.
The ground state in both zero and finite magnetic field consists of real rapidities , where the magnetization sector (and therefore also the number of rapidities) is fixed by the magnetic field, with the field strength acting as a Lagrange multiplier. The quantum numbers of the ground state are
At zero field (), the ground state is the only state with real, finite rapidities, implying there is no room for magnon-like excitations in the quantum numbers. In order to apply a shift in the quantum numbers of the Fermi sea, one therefore has to resort to finite field, since there are now many more possibilities for the quantum numbers available beyond the occupation of the Fermi sea alone.
We define the state with separated Fermi seas at magnetization by applying a shift to the ground state quantum numbers as (we take to be even)
The effect of shifting the quantum numbers, separating the Fermi sea, is illustrated in Fig. 1. We like to call such a state a Moses state.
Four Fermi points are identified, labeled by for the left and right sea, and for the left or right Fermi point respectively. Furthermore, Fermi momenta are defined by , while signs for each left or right Fermi point in a sea are defined as , . With these notations, generalization to separated Fermi seas is straightforward.
3 Dynamical structure factor
The dynamical structure factor constitutes an important connection between inelastic neutron scattering experiments and theory for quantum spin systems. It is directly related to the differential cross section from scattering experiments on the one hand, while it is computable theoretically from both analytical and (exact) numerical methods. The dynamical structure factor is defined as the Fourier transform of the spin-spin correlation
where the label distinguishes the longitudinal and transversal structure factors respectively. Transverse structure factors have been computed using field theoretical methods in Ref. .
For spin chains at zero temperature, the expectation value in the expression for the spin-spin correlation is evaluated with respect to the ground state of Hamiltonian (1). In the context of computations for the state with split Fermi sea, the reference state is taken to be the state defined from the quantum numbers given by Eq. (7).
The dynamical structure factor can be evaluated by inserting a resolution of the identity in Eq. (8), such that the correlator turns into a sum over matrix elements of the split Fermi sea state and its excitations,
The intermediate states are composed of excitations on the state with double Fermi seas. An adaptive scanning procedure through the most relevant intermediate states is applied using the ABACUS algorithm  in order to evaluate the dynamical structure factor. The rapidities of the intermediate states can be obtained by solving Bethe equations (3) with corresponding quantum numbers by an iterative numerical procedure. Subsequently, determinant expressions in terms of the rapidities [30, 31, 32, 33] based on Slavnov’s formula [11, 12] are employed to evaluate the matrix elements . Besides intermediate states consisting of real rapidities, the states can also contain string solutions. In this case, the Bethe-Gaudin-Takahashi equations  in terms of string centers are solved, while their matrix elements can be computed by using reduced determinant expressions .
By integrating over energy and summing over all momenta for Eq. (8), sum rules for the total intensities of the structure factors can be derived, yielding a quantitative measure on the completeness of the truncated matrix element summations. The sum rules are given by
where the magnetization is . Sum rule saturations for all ABACUS computation results of dynamical structure factors used throughout this article are listed in Tab. 1. The table displays significantly lower sum rule saturations for the structure factor computations as opposed to the and structure factors. An intuitive explanation for the computational difficulties of is provided in the Bethe Ansatz language, where the operator removes a rapidity from the state on which it acts, while the operator must add a rapidity to the state. For the latter, it might however be the case that all available quantum number slots around momentum are already filled, such that much more extensive reorganizations of the state are necessary. Thus, computing requires summing over a much more extensive set of intermediate states than for . Given limited computational resources, the saturation of the former will thus be markedly lower than those of the latter.
The longitudinal dynamical structure factors (Eq. (8) with ) are shown in Fig. 2 for various values of the anisotropy, containing plots for the ground state, as well as for states with a double Fermi sea with varying momentum shifts in the quantum numbers. The transverse dynamical structure factors (Eq. (8) with and ) have been plotted in Figs. 3 and 4, respectively.
In order to be able to visualize the delta functions in energy appearing in Eq. (9), Gaussian smoothening has been applied to the data as , the width being of the order of .
The boundaries of the spectra shown in Fig. 2 can be explained by tracking the energies of particle-hole excitations on top of the Moses sea for which either the particle or hole is created at one of the Fermi points. The boundaries of the spectra are equivalent to Fig. 2 in Ref. , for the particle-hole excitations for the double Fermi sea in the Lieb-Liniger model.
In particular, for increasing momentum split , the energy of the reference state increases with respect to the ground state, opening up the possibility to populate branches of the spectrum at negative energy. Moreover, the incommensurate points (at zero energy) start moving in momentum. Broad continua of the spectrum remain with sharply defined thresholds, such that the Tomonaga-Luttinger liquid paradigm retains its validity.
4 Multi-component Tomonaga-Luttinger model
The Hamiltonian of the XXZ model can be mapped exactly to spinless fermions on a lattice by the Jordan-Wigner transformation
with . The Hamiltonian then reads (neglecting chemical potential terms)
The starting point for an effective description of correlations for the split-Fermi-sea states under consideration is the XX point , corresponding to free fermions. The Hamiltonian is then diagonal in momentum space and our state is characterized by four Fermi points .
In constructing a multi-component Tomonaga-Luttinger liquid [35, 36, 37] we follow Ref.  and introduce a species of chiral fermions for each of the Fermi points . In the continuum limit, , we expand the Jordan-Wigner fermion in terms of the chiral fermions as
This is essentially a mode expansion where each of the chiral fermion species encodes the modes close to which determine the correlations, so that the fields can be considered to vary slowly on the scale of the lattice spacing. When , the are non-interacting but have nonlinear dispersion. When we switch on interactions but linearize the dispersion relation and neglect backscattering-like terms, the effective Hamiltonian can be written as
where is the density of species , are effective -ology-like interaction parameters, are the ‘Fermi velocities’ corresponding to the bare, cosine dispersion of the XX model. Note that we define the velocity by taking the derivative of the dispersion to the right, also at left Fermi points. The combination would be positive in equilibrium and corresponding to the Fermi velocity, however, it may be negative in our out-of-equilibrium context. Normal ordering is left implicit.
The validity of Hamiltonian (15) beyond weak interactions to describe correlations cannot be justified by renormalization group arguments in the usual sense since we are describing a high-energy state. Still, the approximations made rely on the idea that we keep the most important terms for the long range physics determined by the modes that are ‘close’ to the Fermi points , which is done by keeping operators of scaling dimension . We will adhere to the equilibrium terminology and call these marginal while terms of higher scaling dimension are called irrelevant.
Note that by considering the momentum of particle or hole excitations constructed by creating a hole or particle at the Fermi points in quantum number space, it is clear that the do not change when we vary the interaction parameter , as the total momentum is completely determined by the quantum number configuration. We thus observe that a kind of generalized Luttinger’s theorem fixes the independent of interactions.
We bosonize the chiral fermions according to
(where ). The Hamiltonian then becomes quadratic in terms of the bosonic fields and can be diagonalized by a Bogoliubov transformation:
This results in the diagonal form of the effective Hamiltonian
where the effective velocities are now related to the dressed dispersion of the XXZ model with . While the interaction parameters cannot be reliably obtained, the Bogoliubov parameters —which also determine the exponents of physical correlation functions—are related to finite-size energy contributions when we extend the filled quantum number blocks by particles at Fermi point . The correction is then to order given by
Here is the number of added (or removed when ) particles corresponding to chiral species and is the energy associated to Fermi point . Eq. (19) gives a relation between the and the finite-size energy differences upon addition or removal of a particle at the Fermi points . Thanks to the properties of the matrix , this relation can in fact be inverted and leads to a way to determine and directly from these finite-size corrections. The generalizes the universal compressibility parameter used in equilibrium situations. In the case of a symmetric combination the derivation may be simplified as detailed in appendix A.
The Bogoliubov parameters have a beautiful interpretation in terms of the phase shifts of the modes at Fermi point upon addition of a particle at Fermi point . This can be argued upon refermionization of the effective Tomonaga-Luttinger theory and can be made precise in terms of the shift function describing the change of the rapidities when the system is excited. In the thermodynamic limit the shift function is determined by the integral equation
where . The relation to the Bogoliubov parameters is then
which can be shown by comparing the finite size corrections to the energy. A derivation of this relation will be presented elsewhere  (also see appendix B). In equilibrium it is known  and the shift function plays an important role in going beyond the Luttinger liquid approximation in computing dynamic correlation functions [40, 41].
Physical correlations generally translate into (products of) two-point functions of vertex operators in bosonized language, which in our conventions are evaluated according to
Here, is measured in units of the lattice spacing.
Asymptotes of spin correlation functions are now obtained by applying the Jordan-Wigner transformation, taking the continuum limit and using the bosonization identities in order to obtain a correlator of the bosonic fields . The Bogoliubov transformation then expresses this in terms of the free fields for which correlation functions are easily evaluated, leading to an expression involving the .
For example, the prediction for the real space longitudinal correlation from the multi-component Tomonaga-Luttinger model is
Here, non-universal prefactors are included for the fluctuating terms and behave as . These do not follow from the Tomonaga-Luttinger construction but can be obtained from the finite-size scaling of matrix element detailed in the next section.
The real space transverse correlation is given from the multi-component Luttinger theory as
where and again non-universal prefactors, now denoted , are included. The exponents are
In the next section we discuss how all parameters are obtained from numerical evaluation of the spectrum and matrix elements and how these predictions compare to real-space correlations obtained from the ABACUS data.
5 Real space correlations
The real space spin-spin correlations of the double Fermi sea state are directly obtained from the ABACUS dynamical structure factor data from section 3 by applying an inverse Fourier transform,
All computations have been carried out at system size at half-magnetization , , for various values of anisotropy and momentum shifts in the Fermi seas. The sum rule saturations of all ABACUS data used in this section are listed in Tab. 1. The results for the static real space correlations () are plotted as data points in Figs. LABEL:fig:mcmpLLSzz-LABEL:fig:mcmpLLSzzshort (longitudinal, ) and Figs. LABEL:fig:mcmpLLSmp-LABEL:fig:mcmpLLSmpshort (transverse, ) respectively. The multi-component Tomonaga-Luttinger model predictions for the correlations with fitted parameters from integrability are incorporated in the figures as well. The predictions for real-space correlations from the transverse structure factor with differ only at from . Although the sum-rule saturation for this correlation (see Tab. 1) is considerably less, we have checked that the fit of the real-space correlation is comparable to the case .
The longitudinal real space correlation from Eq. (24) requires the determination of three classes of parameters. First, the parameters can be deduced from the behavior of the energy upon removal or addition of particles to all four Fermi points. We therefore consider the second derivative of the finite size corrections to the Tomonaga-Luttinger model (Eq. (19)) and define a matrix as
By considering all possible combinations of adding and removing particles to the four Fermi points, the second derivative with respect to energy can be calculated from the energy levels obtained from Bethe Ansatz (Eq. (4)). Subsequently, the eigenvectors of the matrix yield all the parameters .
The remaining non-universal prefactors and exponents can be obtained using the system size scaling behavior of the Umklapp matrix elements, from the relation
where is defined as the Umklapp state by removing a particle at the Fermi point labeled by and placing it back at the Fermi point labeled by . By scaling the system size , the parameters and are directly obtained by a linear fit to the logarithm of Eq. (30) for all combinations of Umklapp states. This procedure is repeated for all values of anisotropy and momentum shifts in the quantum numbers considered here. The values of the prefactors and exponents are plotted in Fig. 5 as function of anisotropy for a fixed value of the momentum shift to the Fermi seas.
In order to compare the multi-component Tomonaga-Luttinger model predictions to finite size results, a conformal transformation
is applied to the scaling behavior of Eq. (24). The resulting expressions for the longitudinal real space correlations, along with the parameters obtained by the procedure described above, are plotted in Figs. LABEL:fig:mcmpLLSzz and LABEL:fig:mcmpLLSzzshort. Fig. LABEL:fig:mcmpLLSzz shows the correspondence of the multi-component Tomonaga-Luttinger model to the matrix element summations obtained from ABACUS at , , for a momentum shift in the Fermi seas of . For all distances but the very smallest, both approaches show good agreement. Fig. LABEL:fig:mcmpLLSzzshort displays a comparison at very short distances for different momentum shifts and anisotropy, still showing a large agreement in both approaches, in a regime where the Tomonaga-Luttinger model is not a priori expected to give bonafide predictions.
Similar to the determination of the non-universal prefactors and exponents for the longitudinal case, the scaling relation
allows for the determination of the parameters and for the purpose of Eq. (26). The Umklapp state is defined by the removal of a particle at the Fermi point labeled by , while dictates the direction of the shift in the quantum numbers due to change in the parity after changing the number of particles. Again, the prefactors and exponents are obtained by fitting the finite size scaling behavior of Eq. (32) and their values are plotted in Fig. 5 as function of anisotropy for a fixed value of the momentum shift to the Fermi seas.
The transverse real space correlations from Eq. (26) are plotted for several values of the anisotropy in Fig. LABEL:fig:mcmpLLSmp for a fixed value of the momentum shift and in Fig. LABEL:fig:mcmpLLSmpshort for short distances and three separate values of the momentum shift, respectively. The non-universal prefactors and exponents are obtained by the aforementioned method, while the conformal transformation to finite size from Eq. (31) has been applied as well. Both figures show again perfect agreement for large distances, while the agreement is also good for short distances with respect to system size. The smallest momentum shift ( in Fig. LABEL:fig:mcmpLLSmpshort) shows the worst agreement at very short distances.
Finally, all previous procedures have been applied to a state where an asymmetric momentum shift is employed to separate the Fermi seas. Fig. 6 shows the corresponding longitudinal and transverse dynamical structure factors and real space correlations obtained from ABACUS at , , and the multi-component Tomonaga-Luttinger model. The parameters for the latter have been obtained from the fitting procedure described in this section, applied to system size scaling behavior of Umklapp matrix elements on the asymmetric Fermi points. Once again, the real space correlations display agreement for both the asymptotics as well as for short distances for both approaches.
6 Time dependence of correlations
The fact that the zero-entropy states we are considering are far from equilibrium is not visible when we restrict attention to the static correlations, which would be the same if this state was obtained as the ground state of a different Hamiltonian. In order to make the out-of-equilibrium nature apparent we have to probe the energies of ‘excitations’, i.e. modifications of the Moses sea by creating additional particles and holes, which may now have both positive and negative energy differences with respect to the reference state. A physically meaningful way to probe the energy landscape is by computing time-dependent correlations which can in principle be related to observable quantities. These are already encoded in the dynamical structure factors computed with ABACUS and can be obtained by Fourier transformation.
Recent years have witnessed a revolutionary increase in understanding of dynamical correlations in critical one-dimensional systems from the perspective of both effective field theory methods and integrability [42, 43, 44, 45, 46, 47, 48, 49, 40, 50, 51, 52, 53, 54, 17, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67]. The threshold behavior of many dynamical correlations in energy and momentum space can be understood in terms of specific configurations of particle and hole excitations. These lead to a scattering phase shift of the modes close to the Fermi energy which is identified as the cause of the characteristic power-law singularities by means of Anderson’s orthogonality catastrophe. This threshold behavior, which also determines the asymptotic behavior of the correlations in real space and time, is then described by an effective model in which the high energy particle or hole excitation is treated as a mobile impurity interacting with the low-energy modes.
We generalize this mobile impurity approach to the present out-of-equilibrium context by extending our multi-component Tomonaga-Luttinger model to include the appropriate impurity configurations and interactions. To be specific, we will compute the spin autocorrelation
where denotes the Jordan-Wigner fermion and we used translational invariance. By imagining to obtain as a Fourier transform in -space taking the -integral first, one can argue that as a function of singular behavior stems from the ‘Fermi points’ and points where the edge of support has a tangent with vanishing velocity. This identifies the important impurity configurations for this function as corresponding to a particle or hole with vanishing velocity, i.e. either at the bottom or the top of the band. Let us assume that the Moses state leaves the corresponding quantum numbers unoccupied, which is valid for a symmetric configuration with an even number of seas. This means that there are only high-energy particle impurities. We will use an index to label the particle at the bottom/top of the band respectively. The mobile impurity model becomes
Note that the last term is just a density-density interaction between the impurity modes and the chiral fermions parametrized by the coupling constants , while the first two terms correspond to the impurity dispersions and the TL modes. In the small limit all parameters in can be obtained from in Eq. (13) as in appendix C. In general they can be obtained from integrability.
The impurity modes are decoupled from the Tomonaga-Luttinger model up to irrelevant operators by the unitary transformation
The effect is that correlators of the TL model are still computed in the same way, but the impurity operator obtains an extra vertex operator in terms of the bosonic modes
For the computation of the autocorrelation , the crucial observation is now that the asymptotic behavior, determined by the behavior around a few singular points of the longitudinal structure factor, is well captured by certain correlations computable using the Hamiltonian . In marked contrast to the equilibrium case, the TL model does not account for zero-energy states only, and therefore even the contributions to Eq. (33) that do not involve the impurity will display the energy difference of the Fermi points leading to fluctuating terms
where is the energy associated to Fermi point . The TL contributions sum up to an expression similar to the static correlation in Eq. (23). The impurity contributions are of the form
Using the impurity correlator
we find the result
The prefactor can in principle be obtained from finite-size scaling of matrix elements similar to but with the Umklapp state replaced by the appropriate impurity state . We have checked this expression for the autocorrelation against Fourier transformed ABACUS data for small values of , and find that it converges to the exact result on quite short time scales for a configuration when the Fermi points are well away from the band top and bottom at (Fig. LABEL:fig:autoplot), but the correspondence for short to moderate times becomes noticeably worse when we decrease the separation between the two seas. This could be a finite size effect since the number of states in between the impurity mode and the Fermi edges becomes small, but rather we believe that the correlation is not well-captured by the impurity model in that case as a clear separation in sub-bands becomes questionable.
When the current mobile-impurity approach works well, this tells us that the time-dependent correlation is determined by the modes very close to the Fermi points which correspond to particle-hole excitations involving only the quantum numbers close to the edges of the two Fermi seas. The role of the spectrum at the Fermi points and of the impurity is two-fold: (i) The energy differences determine the frequencies of fluctuations. (ii) The Fermi velocities and impurity mass change the prefactor of the separate terms. Note that the decay of the correlation on the other hand is determined by energy independent data, namely by the appropriate phase shifts and Anderson’s orthogonality catastrophe, very similar to the equilibrium case.
We have considered high energy zero-entropy states for the anisotropic Heisenberg spin chain defined by a double-Fermi sea quantum-number configuration. Our focus was on dynamical correlations computed by summing over relevant matrix elements of particle-hole excitations at finite system size, the matrix elements being given by algebraic Bethe Ansatz. Correlations in real space and time are obtained from these by numerical Fourier transformation.
Zero-entropy Bethe states provide an interesting class of eigenstates of integrable models which, while far from equilibrium, share many features with the ground state. In particular we have shown that, when in the critical regime of the anisotropic Heisenberg spin chain, these states display critical correlations characterized by fluctuations with power-law decay. Starting from the Bethe Ansatz solution, we have constructed an effective field theory in terms of a multi-component Tomonaga-Luttinger model capturing these correlations with great accuracy for large, but also surprisingly short, distances in real space. Similar to the ground state case, fluctuations correspond to the differences of the generalized Fermi momenta , the power-law exponents are related to a Bogoliubov transformation diagonalizing the effective Hamiltonian, while non-universal prefactors determine the relative amplitudes of the fluctuating terms.
The Bogoliubov parameters can be related to finite-size energy differences upon adding or removing particles at the Fermi points, and the non-universal correlation prefactors to the finite-size scaling of matrix elements with Umklapp excitations. We have used this to obtain all field theory parameters from integrability yielding parameter-free fits for the static correlations. A surprising fact is that, while static correlations do not know about the Hamiltonian for time evolution, the energies computed to obtain the do feature the Fermi velocities from the spectrum of the actual XXZ Hamiltonian, which does not determine the statistical ensemble in the present case. This implies that the relation of the spectrum and the Bogoliubov parameters is universally valid irrespective of the energy function one uses.
The parameters turn out to correspond to the phase shifts of the modes at one of the generalized Fermi points upon creating an excitation at Fermi point . At zero temperature, the Luttinger parameter can similarly be regarded as a parametrization of the phase shift at the left and right Fermi points upon creating single-particle excitations with momentum , and it has recently been realized that this implies a universal description of dynamical correlations that go beyond the linear-spectrum approximation made in the Tomonaga-Luttinger model . The applicability of the multi-component Tomonaga-Luttinger model suggests that, similar to the ground state case, a universal description of dynamical correlations in a domain close to the Fermi points may be obtained. The description of domains near the edge of support is likely captured by a mobile impurity model. As a first step of this generalized use of what is known as non-linear Luttinger liquid theory, we have compared the longitudinal autocorrelation function to predictions obtained from a mobile impurity model at small . The correlation seems to converge to the mobile impurity result in quite reasonable times at least when the two Fermi seas are well-separated.
In conclusion, we have presented results on dynamical correlations of zero-entropy states in the anisotropic Heisenberg chain. The distinctive features, which may serve for identification in experiment, can be understood by adapting the familiar ground state reasoning. By making the appropriate adjustments to equilibrium techniques based on the Tomonaga-Luttinger model many aspects can furthermore be understood in great detail and with quantitative agreement once a handful of parameters is fixed from integrability. Although we have focused on states which have vanishing entropy density, we may consider thermal-like dressings to the split seas. In the Tomonaga-Luttinger description, finite temperatures are treated by a simple functional replacement for the fundamental correlators. On the side of integrability, recent work has shown that finite temperature correlators are also numerically accessible, at least for the Lieb-Liniger model . How the correspondence between integrability results and the field theory works out in split-sea configurations at finite temperature remains to be investigated.
We thank Yuri van Nieuwkerk for useful comments. The authors acknowledge support from the Foundation for Fundamental Research on Matter (FOM) and from the Netherlands Organization for Scientific Research (NWO). This work forms part of the activities of the Delta-Institute for Theoretical Physics (D-ITP).
The work of R. P. V. and J.-S. C. was supported by NWO VICI grant number 680-47-605. The work of I. S. E. and J.-S. C. was supported by FOM program 128, ‘The singular physics of 1D electrons’.
Appendix A Simplification for symmetric seas
In the case of a symmetric configuration of seas, we have and and in general that the system is symmetric under simultaneously exchanging and . It is then convenient to combine the fields at Fermi points at opposite momenta and define
with inverse transformation
The Hamiltonian for the multi-component TL model can then be written as
with and .
In order to respect the canonical commutation relations and diagonalize we introduce new fields according to
where the are related to the Bogoliubov parameters as
by virtue of the symmetry . The effective Hamiltonian takes the familiar diagonal form
The matrix is straightforwardly obtained from calculations in finite size and finite particle number from corrections to the energy upon creating particle-number or current excitations
in terms of the numbers of particles added at Fermi point .
Appendix B Moses states and generalized TBA
Finite-size corrections to the spectrum for the ground state are usually discussed by taking the thermodynamic limit of the Bethe equations and the energy, keeping finite size corrections using the Euler-Maclaurin formula . The class of Moses states can be considered as the zero temperature limit of a generalized Gibbs ensemble for which a generalized thermodynamic Bethe Ansatz and a treatment of finite-size corrections similar to the ground state exists [70, 71]. However, in our case we still use the energy dispersion obtained from the true XXZ Hamiltonian determining the time-evolution rather than the statistical ensemble (which we define by hand in the micro-canonical sense). Since we also obtain the critical exponents using the true XXZ Hamiltonian and not a GGE Hamiltonian, the fact that the results coincide requires clarification. The standard derivation of the finite size corrections of the energy makes use of the property that the corresponding dressed energy vanishes at the Fermi points, but in our case if we measure energies by using the true XXZ Hamiltonian.
Upon changing the particle number at one of the Fermi points the extremal rapidities experience a shift of order . The change in the Fermi points and rapidities can be shown to satisfy
The finite-size corrections to the energy due to adding or removing particles from the extremities of the Fermi-sea-like blocks is then
where is the energy of a particle (hole) with rapidity on top of the Moses state. In contrast to the equilibrium case where energy is measured or a GGE dressed energy satisfying , the funtion cannot be defined through an integral equation but has an additional contribution stemming from the finite energy at the Fermi points. Details will be presented elsewhere.
Appendix C Perturbative expressions for effective-field-theory parameters
A convenient way to obtain perturbative expressions for the exponents of asymptotes of correlation functions is to derive the effective Hamiltonian to first order in the interaction. Starting from the Hamiltonian of spinless fermions on the lattice we use the mode expansion
and bosonize the chiral fermions. The kinetic term leads to the velocities for the chiral fermions and the impurity parameters in the noninteracting limit and . The interaction term
renormalizes these values. By plugging in the bosonization identities, normal ordering and neglecting irrelevant terms we get the first order in expressions
Next, we focus on the terms from the interaction involving the impurity. This leads to a density-density interaction of the impurity modes with the chiral fermions with parameters
There is also a first order correction to the impurity energy appearing from Eq. (54) from the terms proportional to (after normal ordering), which is
This corresponds to the Hartree-Fock correction
with , which corresponds to the interaction potential in Eq. (54): .
The non-universal prefactors can also be obtained perturbatively using the methods discussed in Ref. , but we have not done this calculation.
-  F. D. M. Haldane, ‘Luttinger liquid theory’ of one-dimensional quantum fluids. I. Properties of the Luttinger model and their extension to the general 1D interacting spinless Fermi gas, J. Phys C: Sol. St. Phys. 14(19), 2585 (1981), doi:10.1088/0022-3719/14/19/010.
-  F. D. M. Haldane, Demonstration of the “Luttinger liquid” character of Bethe-ansatz-soluble models of 1-D quantum fluids, Phys. Lett. A 81(2-3), 153 (1981), doi:10.1016/0375-9601(81)90049-9.
-  A. O. Gogolin, A. A. Nersesyan and A. M. Tsvelik, Bosonization and Strongly Correlated Systems, Cambridge University Press (1998).
-  J. von Delft and H. Schoeller, Bosonization for Beginners - Refermionization for Experts, Ann. Phys. (Leipzig) 7(4), 225 (1998), doi:10.1002/(SICI)1521-3889(199811)7:4<225::AID-ANDP225>3.0.CO;2-L.
-  H. A. Bethe, Zur Theorie der Metalle. i. Eigenwerte und Eigenfunktionen der linearen Atomkette, Zeit. für Physik 71, 205 (1931), doi:10.1007/BF01341708.
-  E. H. Lieb and W. Liniger, Exact Analysis of an Interacting Bose Gas. I. The General Solution and the Ground State, Phys. Rev. 130(4), 1605 (1963), doi:10.1103/PhysRev.130.1605.
-  E. H. Lieb, Exact Analysis of an Interacting Bose Gas. II. The Excitation Spectrum, Phys. Rev. 130(4), 1616 (1963), doi:10.1103/PhysRev.130.1616.
-  W. Heisenberg, Zur Theorie des Ferromagnetismus, Z. Phys. 49, 619 (1928), doi:10.1007/BF01328601.
-  R. Orbach, Linear antiferromagnetic chain with anisotropic coupling, Phys. Rev. 112(2), 309 (1958), doi:10.1103/PhysRev.112.309.
-  V. E. Korepin, N. M. Bogoliubov and A. G. Izergin, Quantum Inverse Scattering Method and Correlation Functions, Cambridge Univ. Press (1993).
-  N. A. Slavnov, Calculation of scalar products of wave functions and form factors in the framework of the algebraic Bethe Ansatz, Theor. Math. Phys. 79, 502 (1989), doi:10.1007/BF01016531.
-  N. A. Slavnov, Nonequal-time current correlation function in a one-dimensional bose gas, Theor. Math. Phys. 82, 273 (1990), doi:10.1007/BF01029221.
-  T. Giamarchi, Quantum Physics in One Dimension, Oxford University Press (2004).
-  N. Kitanine, K. K. Kozlowski, J. M. Maillet, N. A. Slavnov and V. Terras, On the thermodynamic limit of form factors in the massless XXZ Heisenberg chain, J. Math. Phys. 50(9), 095209 (2009), doi:10.1063/1.3136683.
-  N. Kitanine, K. K. Kozlowski, J. M. Maillet, N. A. Slavnov and V. Terras, Algebraic Bethe ansatz approach to the asymptotic behavior of correlation functions, J. Stat. Mech.: Th. Exp. 2009(04), P04003 (2009), doi:10.1088/1742-5468/2009/04/P04003.
-  A. Shashi, L. I. Glazman, J.-S. Caux and A. Imambekov, Nonuniversal prefactors in the correlation functions of one-dimensional quantum liquids, Phys. Rev. B 84(4), 045408 (2011), doi:10.1103/PhysRevB.84.045408.
-  N. Kitanine, K. K. Kozlowski, J. M. Maillet, N. A. Slavnov and V. Terras, A form factor approach to the asymptotic behavior of correlation functions in critical models, J. Stat. Mech.: Th. Exp. 2011(12), P12010 (2011), doi:10.1088/1742-5468/2011/12/P12010.
-  M. de Llano and J. P. Vary, Generalized Fermi sea for plane-wave Hartree-Fock theory: One dimensional model calculation, Phys. Rev. C 19, 1083 (1979), doi:10.1103/PhysRevC.19.1083.
-  M. Kormos, Y.-Z. Chou and A. Imambekov, Exact Three-Body Local Correlations for Excited States of the 1D Bose Gas, Phys. Rev. Lett. 107, 230405 (2011), doi:10.1103/PhysRevLett.107.230405.
-  T. Fokkema, I. S. Eliëns and J.-S. Caux, Split Fermi seas in one-dimensional Bose fluids, Phys. Rev. A 89, 033637 (2014), doi:10.1103/PhysRevA.89.033637.
-  V. Alba, M. Fagotti and P. Calabrese, Entanglement entropy of excited states, Journal of Statistical Mechanics: Theory and Experiment 2009(10), P10020 (2009), doi:10.1088/1742-5468/2009/10/P10020.
-  T. Kinoshita, T. Wenger and D. S. Weiss, A quantum Newton’s cradle, Nature 440, 900 (2006), doi:10.1038/nature04693.
-  R. van den Berg, B. Wouters, S. Eliëns, J. De Nardis, R. M. Konik and J.-S. Caux, Separation of time scales in a quantum newton’s cradle, Phys. Rev. Lett. 116, 225302 (2016), doi:10.1103/PhysRevLett.116.225302.
-  L. Vidmar, J. P. Ronzheimer, M. Schreiber, S. Braun, S. S. Hodgman, S. Langer, F. Heidrich-Meisner, I. Bloch and U. Schneider, Dynamical quasicondensation of hard-core bosons at finite momenta, Phys. Rev. Lett. 115, 175301 (2015), doi:10.1103/PhysRevLett.115.175301.
-  H. Frahm and C. Rödenbeck, A generalized spin ladder in a magnetic field, Eur. Phys. J. B 10(3), 409 (1999), doi:10.1007/s100510050870.
-  A. Zvyagin, A. Klümper and J. Zittartz, Integrable correlated electron model with next-nearest-neighbour interactions, Eur. Phys. J. B 19(1), 25 (2001), doi:10.1007/s100510170347.
-  C. N. Yang and C. P. Yang, One-Dimensional Chain of Anisotropic Spin-Spin Interactions. I. Proof of Bethe’s Hypothesis for Ground State in a Finite System, Phys. Rev. 150(1), 321 (1966), doi:10.1103/PhysRev.150.321.
-  H. Karimi and I. Affleck, Transverse spectral functions and dzyaloshinskii-moriya interactions in spin chains, Phys. Rev. B 84, 174420 (2011), doi:10.1103/PhysRevB.84.174420.
-  J.-S. Caux, Correlation functions of integrable models: a description of the ABACUS algorithm, J. Math. Phys. 50(9), 095214 (2009), doi:10.1063/1.3216474.
-  N. Kitanine, J. M. Maillet and V. Terras, Form factors of the XXZ Heisenberg finite chain, Nucl. Phys. B 554(3), 647 (1999), doi:10.1016/S0550-3213(99)00295-3.
-  M. Gaudin, B. M. McCoy and T. T. Wu, Normalization sum for the Bethe’s hypothesis wave functions of the Heisenberg-Ising chain, Phys. Rev. D 23(2), 417 (1981), doi:10.1103/PhysRevD.23.417.
-  V. E. Korepin, Calculation of norms of Bethe wave functions, Commun. Math. Phys. 86, 391 (1982), doi:10.1007/BF01212176.
-  J.-S. Caux, R. Hagemans and J. M. Maillet, Computation of dynamical correlation functions of Heisenberg chains: the gapless anisotropic regime, J. Stat. Mech.: Th. Exp. 2005(09), P09003 (2005), doi:10.1088/1742-5468/2005/09/P09003.
-  M. Takahashi and M. Suzuki, One-Dimensional Anisotropic Heisenberg Model at Finite Temperatures, Prog. Theor. Phys. 48(6), 2187 (1972), doi:10.1143/PTP.48.2187.
-  K. Penc and J. Sólyom, One-dimensional Hubbard model in a magnetic field and the multicomponent Tomonaga-Luttinger model, Phys. Rev. B 47, 6273 (1993), doi:10.1103/PhysRevB.47.6273.
-  J. Voit, One-dimensional Fermi liquids, Rep. Prog. Phys. 58(9), 977 (1995), doi:10.1088/0034-4885/58/9/002.
-  M. Pletyukhov and V. Gritsev, Persistent currents in a multicomponent Tomonaga-Luttinger liquid: Application to a mesoscopic semiconductor ring with spin-orbit interaction, Phys. Rev. B 70, 165316 (2004), doi:10.1103/PhysRevB.70.165316.
-  I. S. Eliëns and J.-S. Caux, General finite-size effects for zero-entropy states in one-dimensional quantum integrable models, arXiv (2016), 1608.00378.
-  V. Korepin and N. Slavnov, The new identity for the scattering matrix of exactly solvable models, Eur. Phys. Jour. B 5, 555 (1998), doi:10.1007/s100510050477.
-  V. V. Cheianov and M. Pustilnik, Threshold Singularities in the Dynamic Response of Gapless Integrable Models, Phys. Rev. Lett. 100(12), 126403 (2008), doi:10.1103/PhysRevLett.100.126403.
-  A. Imambekov and L. I. Glazman, Exact Exponents of Edge Singularities in Dynamic Correlation Functions of 1D Bose Gas, Physical Review Letters 100(20), 206805 (2008), doi:10.1103/PhysRevLett.100.206805.
-  J. M. P. Carmelo and P. D. Sacramento, Finite-energy landau liquid theory for the one-dimensional hubbard model: Pseudoparticle energy bands and degree of localization/delocalization, Phys. Rev. B 68, 085104 (2003), doi:10.1103/PhysRevB.68.085104.
-  A. V. Rozhkov, Fermionic quasiparticle representation of tomonaga-luttinger hamiltonian, Eur. Phys. J. B 47(2), 193 (2005), doi:10.1140/epjb/e2005-00312-3.
-  R. G. Pereira, J. Sirker, J.-S. Caux, R. Hagemans, J. M. Maillet, S. R. White and I. Affleck, Dynamical spin structure factor for the anisotropic spin- heisenberg chain, Phys. Rev. Lett. 96, 257202 (2006), doi:10.1103/PhysRevLett.96.257202.
-  M. Pustilnik, M. Khodas, A. Kamenev and L. I. Glazman, Dynamic response of one-dimensional interacting fermions, Phys. Rev. Lett. 96, 196405 (2006), doi:10.1103/PhysRevLett.96.196405.
-  M. Khodas, M. Pustilnik, A. Kamenev and L. I. Glazman, Dynamics of excitations in a one-dimensional bose liquid, Phys. Rev. Lett. 99, 110405 (2007), doi:10.1103/PhysRevLett.99.110405.
-  M. Khodas, M. Pustilnik, A. Kamenev and L. I. Glazman, Fermi-luttinger liquid: Spectral function of interacting one-dimensional fermions, Phys. Rev. B 76(15), 155402 (2007), doi:10.1103/PhysRevB.76.155402.
-  A. Imambekov and L. I. Glazman, Universal Theory of Nonlinear Luttinger Liquids, Science 323(5911), 228 (2009), doi:10.1126/science.1165403.
-  R. G. Pereira, S. R. White and I. Affleck, Exact edge singularities and dynamical correlations in spin- chains, Phys. Rev. Lett. 100, 027206 (2008), doi:10.1103/PhysRevLett.100.027206.
-  R. G. Pereira, S. R. White and I. Affleck, Spectral function of spinless fermions on a one-dimensional lattice, Phys. Rev. B 79, 165113 (2009), doi:10.1103/PhysRevB.79.165113.
-  A. Imambekov and L. I. Glazman, Phenomenology of one-dimensional quantum liquids beyond the low-energy limit, Phys. Rev. Lett. 102, 126405 (2009), doi:10.1103/PhysRevLett.102.126405.
-  F. H. L. Essler, Threshold singularities in the one-dimensional hubbard model, Phys. Rev. B 81, 205120 (2010), doi:10.1103/PhysRevB.81.205120.
-  T. L. Schmidt, A. Imambekov and L. I. Glazman, Spin-charge separation in one-dimensional fermion systems beyond luttinger liquid theory, Phys. Rev. B 82, 245104 (2010), doi:10.1103/PhysRevB.82.245104.
-  T. L. Schmidt, A. Imambekov and L. I. Glazman, Fate of 1d spin-charge separation away from fermi points, Phys. Rev. Lett. 104, 116403 (2010), doi:10.1103/PhysRevLett.104.116403.
-  K. K. Kozlowski and V. Terras, Long-time and large-distance asymptotic behavior of the current–current correlators in the non-linear schrödinger model, Journal of Statistical Mechanics: Theory and Experiment 2011(09), P09013 (2011), doi:10.1088/1742-5468/2011/09/P09013.
-  R. G. Pereira, K. Penc, S. R. White, P. D. Sacramento and J. M. P. Carmelo, Charge dynamics in half-filled hubbard chains with finite on-site interaction, Phys. Rev. B 85, 165132 (2012), doi:10.1103/PhysRevB.85.165132.
-  A. Imambekov, T. L. Schmidt and L. I. Glazman, One-dimensional quantum liquids: Beyond the luttinger liquid paradigm, Rev. Mod. Phys. 84, 1253 (2012), doi:10.1103/RevModPhys.84.1253.
-  R. G. Pereira, Long time correlations of nonlinear luttinger liquids, Int. J. Mod. Phys. B 26(22), 1244008 (2012), doi:10.1142/S0217979212440080.
-  L. Seabra, F. H. L. Essler, F. Pollmann, I. Schneider and T. Veness, Real-time dynamics in the one-dimensional hubbard model, Phys. Rev. B 90, 245127 (2014), doi:10.1103/PhysRevB.90.245127.
-  T. Price and A. Lamacraft, Fine structure of the phonon in one dimension from quantum hydrodynamics, Phys. Rev. B 90, 241415 (2014), doi:10.1103/PhysRevB.90.241415.
-  S. Gangadharaiah, T. L. Schmidt and D. Loss, Structure factor of interacting one-dimensional helical systems, Phys. Rev. B 89, 035131 (2014), doi:10.1103/PhysRevB.89.035131.
-  C. Karrasch, R. G. Pereira and J. Sirker, Low temperature dynamics of nonlinear luttinger liquids, New J. Phys. 17(10), 103003 (2015).
-  F. H. L. Essler, R. G. Pereira and I. Schneider, Spin-charge-separated quasiparticles in one-dimensional quantum fluids, Phys. Rev. B 91, 245150 (2015), doi:10.1103/PhysRevB.91.245150.
-  K. K. Kozlowski and J. M. Maillet, Microscopic approach to a class of 1d quantum critical models, Journal of Physics A: Mathematical and Theoretical 48, 484004 (2015).
-  K. K. Kozlowski, Large-distance and long-time asymptotic behavior of the reduced density matrix in the non-linear schrödinger model, In Annales Henri Poincaré, vol. 16, pp. 437–534. Springer, doi:10.1007/s00023-014-0327-3 (2015).
-  T. Veness and F. H. L. Essler, Mobile impurity approach to the optical conductivity in the hubbard chain, Phys. Rev. B 93, 205101 (2016), doi:10.1103/PhysRevB.93.205101.
-  I. S. Eliëns, F. B. Ramos, J. C. Xavier and R. G. Pereira, Boundary versus bulk behavior of time-dependent correlation functions in one-dimensional quantum systems, Phys. Rev. B 93, 195129 (2016), doi:10.1103/PhysRevB.93.195129.
-  A. Shashi, M. Panfil, J.-S. Caux and A. Imambekov, Exact prefactors in static and dynamic correlation functions of one-dimensional quantum integrable models: Applications to the calogero-sutherland, lieb-liniger, and models, Phys. Rev. B 85, 155136 (2012), doi:10.1103/PhysRevB.85.155136.
-  M. Panfil and J.-S. Caux, Finite-temperature correlations in the Lieb-Liniger one-dimensional Bose gas, Phys. Rev. A 89, 033605 (2014), doi:10.1103/PhysRevA.89.033605.
-  J. Mossel and J.-S. Caux, Generalized TBA and generalized Gibbs, J. Phys. A: Math. Theor. 45(25), 255001 (2012), doi:10.1088/1751-8113/45/25/255001.
-  E. Eriksson and V. Korepin, Finite-size effects from higher conservation laws for the one-dimensional Bose gas, J. Phys. A: Math. Theor. 46(23), 235002 (2013), doi:10.1088/1751-8113/46/23/235002.