The steady state after a quantum quench from the Néel state to the anisotropic Heisenberg model for spin chains is investigated. Two methods that aim to describe the postquench non-thermal equilibrium, the generalized Gibbs ensemble and the quench action approach, are discussed and contrasted. Using the recent implementation of the quench action approach for this Néel-to-XXZ quench, we obtain an exact description of the steady state in terms of Bethe root densities, for which we give explicit analytical expressions.Furthermore, by developing a systematic small-quench expansion around the antiferromagnetic Ising limit, we analytically investigate the differences between the predictions of the two methods in terms of densities and postquench equilibrium expectation values of local physical observables. Finally, we discuss the details of the quench action solution for the quench to the isotropic Heisenberg spin chain. For this case we validate the underlying assumptions of the quench action approach by studying the large-system-size behavior of the overlaps between Bethe states and the Néel state.

The dynamical structure factor of the Babujan-Takhtajan antiferromagnetic spin-1 chain is computed numerically at zero temperature and zero magnetic field, using the higher-spin generalization of an Algebraic Bethe Ansatz-based method previously used for spin-1/2 integrable chains. This method, which consists in the explicit construction of eigenstates and the summation of the Lehmann representation of the correlator, is here particularly challenging to implement in view of the presence of strongly deviated string solutions to the Bethe equations. We show that a careful treatment of these deviations makes it possible to obtain perfect saturation of sum rules for small system sizes, and extremely good saturation for large system sizes where the dynamical structure factor is computed by including all two-spinon and four-spinon contributions. The real-space spin-spin correlation, obtained by Fourier transforming our results, displays asymptotics fitting predictions from conformal field theory.

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. SciPost Physics SubmissionA Simplification for symmetric seas 23 B Moses states and generalized TBA 24 C Perturbative expressions for effective-field-theory parameters 24References 25

The quantum scattering of magnon bound states in the anisotropic Heisenberg spin chain is shown to display features similar to the scattering of solitons in classical exactly solvable models. Localized colliding Gaussian wave packets of bound magnons are constructed from string solutions of the Bethe equations and subsequently evolved in time, relying on an algebraic Bethe ansatz based framework for the computation of local expectation values in real space-time. The local magnetization profile shows the trajectories of colliding wave packets of bound magnons, which obtain a spatial displacement upon scattering. Analytic predictions on the displacements for various values of anisotropy and string lengths are derived from scattering theory and Bethe ansatz phase shifts, matching time evolution fits on the displacements. The time evolved block decimation (TEBD) algorithm allows for the study of scattering displacements from spin-block states, showing similar scattering displacement features.Comment: 15 pages, 7 figures. (v2: citations added