In integrable many-particle systems, it is widely believed that the stationary state reached at late times after a quantum quench can be described by a generalized Gibbs ensemble (GGE) constructed from their extensive number of conserved charges. A crucial issue is then to identify a complete set of these charges, enabling the GGE to provide exact steady-state predictions. Here we solve this long-standing problem for the case of the spin-1=2 Heisenberg chain by explicitly constructing a GGE which uniquely fixes the macrostate describing the stationary behavior after a general quantum quench. A crucial ingredient in our method, which readily generalizes to other integrable models, are recently discovered quasilocal charges. As a test, we reproduce the exact postquench steady state of the Néel quench problem obtained previously by means of the Quench Action method. Introduction.-Understanding and describing the equilibration of isolated many-particle systems is one of the main current challenges of quantum physics. The presence of higher conserved charges (above the Hamiltonian) is linked to the absence of full relaxation to a thermalized state; the conjectured appropriate framework to characterize the steady-state properties in such a situation is the generalized Gibbs ensemble (GGE) [1], in which all available charges are ascribed an individual "chemical potential" set by the initial conditions, and the steady state is the maximal entropy state fulfilling all the constraints associated to the conserved charges . The basic idea underlying the GGE is as follows. Let H ≃ H ð1Þ be the Hamiltonian of an integrable model, and fH ðnÞ g a set of conserved charges fulfilling ½H ðnÞ ; H ðmÞ ¼ 0. The situation we are interested in is that of a quantum quench, where we initially prepare our system in the ground state jΨð0Þi of a local Hamiltonian H 0 and then consider unitary time evolution with respect to our integrable Hamiltonian

We review recent progress in understanding the notion of locality in integrable quantum lattice systems. The central concept are the so-called quasilocal conserved quantities, which go beyond the standard perception of locality. Two systematic procedures to rigorously construct families of quasilocal conserved operators based on quantum transfer matrices are outlined, specializing on anisotropic Heisenberg XXZ spin-1/2 chain. Quasilocal conserved operators stem from two distinct classes of representations of the auxiliary space algebra, comprised of unitary (compact) representations, which can be naturally linked to the fusion algebra and quasiparticle content of the model, and non-unitary (non-compact) representations giving rise to charges, manifestly orthogonal to the unitary ones. Various condensed matter applications in which quasilocal conservation laws play an essential role are presented, with special emphasis on their implications for anomalous transport properties (finite Drude weight) and relaxation to non-thermal steady states in the quantum quench scenario.

For fundamental integrable quantum chains with deformed symmetries we outline a general procedure for defining a continuous family of quasilocal operators whose time derivative is supported near the two boundary sites only. The program is implemented for a spin 1/2 XXZ chain, resulting in improved rigorous estimates for the high temperature spin Drude weight.

Nonergodic dynamical systems display anomalous transport properties. Prominent examples are integrable quantum systems, whose exceptional properties are diverging dc conductivities. In this Letter, we explain the microscopic origin of ideal conductivity by resorting to the thermodynamic particle content of a system. Using group-theoretic arguments we rigorously resolve the long-standing controversy regarding the nature of spin and charge Drude weights in the absence of chemical potentials. In addition, by employing a hydrodynamic description, we devise an efficient computational method to calculate exact Drude weights from the stationary currents generated in an inhomogeneous quench from bipartitioned initial states. We exemplify the method on the anisotropic Heisenberg model at finite temperatures for the entire range of anisotropies, accessing regimes that are out of reach with other approaches. Quite remarkably, spin Drude weight and asymptotic spin current rates reveal a completely discontinuous (fractal) dependence on the anisotropy parameter.

We present an identification of the spectra of local conserved operators of integrable quantum lattice models and the density distributions of their thermodynamic particle content. This is derived explicitly for the Heisenberg XXZ spin chain. As an application we discuss a quantum quench scenario, in both the gapped and critical regimes. We outline an exact technique which allows for an efficient implementation on periodic matrix product states. In addition, for certain simple product states we obtain closed-form expressions for the density distributions in terms of solutions to Hirota difference equations. Remarkably, no reference to a maximal entropy principle is invoked.ArXiv ePrint: 1512.04454

We present a general formulation of Floquet states of periodically time-dependent open Markovian quasi-free fermionic many-body systems in terms of a discrete Lyapunov equation. Illustrating the technique, we analyze periodically kicked XY spin 1/2 chain which is coupled to a pair of Lindblad reservoirs at its ends. A complex phase diagram is reported with re-entrant phases of long range and exponentially decaying spin-spin correlations as some of the system's parameters are varied. The structure of phase diagram is reproduced in terms of counting non-trivial stationary points of Floquet quasi-particle dispersion relation. Introduction. Understanding and controlling dynamics of many-body quantum systems when they are open to the environment and driven far from equilibrium is an exciting and important topic of current research in theoretical [1, 2] and experimental quantum physics [3]. In particular, since it has been recently realized that certain emergent phenomena, such as quantum phase transitions and long range order -previously known only in equilibrium zero-temperature quantum states [4] -can appear also in far from equilibrium steady states of quantum Liouville evolution [2,[5][6][7]. In investigating dynamical and critical many-body phenomena, quasi-free (quadratic) quantum systems play an important role as they are amenable to analytical treatment (see e.g. [8]), so many effects can be analyzed exactly or in great detail. For example, quantum phase transitions in nonequilibrium steady states have been observed either in quasi-free [5,9], or strongly interacting [6], or even dissipative [7,10] quantum systems in one dimension.

We outline a general formalism of hydrodynamics for quantum systems with multiple particle species which undergo completely elastic scattering. In the thermodynamic limit, the complete kinematic data of the problem consists of the particle content, the dispersion relations, and a universal dressing transformation which accounts for interparticle interactions. We consider quantum integrable models and we focus on the one-dimensional fermionic Hubbard model. By linearizing hydrodynamic equations, we provide exact closed-form expressions for Drude weights, generalized static charge susceptibilities and charge-current correlators valid on hydrodynamic scale, represented as integral kernels operating diagonally in the space of mode numbers of thermodynamic excitations. We find that, on hydrodynamic scales, Drude weights manifestly display Onsager reciprocal relations even for generic (i.e. non-canonical) equilibrium states, and establish a generalized detailed balance condition for a general quantum integrable model. We present the first exact analytic expressions for the general Drude weights in the Hubbard model, and explain how to reconcile different approaches for computing Drude weights from the previous literature.

We identify a class of one-dimensional spin and fermionic lattice models which display diverging spin and charge diffusion constants, including several paradigmatic models of exactly solvable strongly correlated many-body dynamics such as the isotropic Heisenberg spin chains, the Fermi-Hubbard model, and the t-J model at the integrable point. Using the hydrodynamic transport theory, we derive an analytic lower bound on the spin and charge diffusion constants by calculating the curvature of the corresponding Drude weights at half filling, and demonstrate that for certain lattice models with isotropic interactions some of the Noether charges exhibit super-diffusive transport at finite temperature and half filling.