Entropy is a fundamental concept in equilibrium statistical mechanics, yet its origin in the nonequilibrium dynamics of isolated quantum systems is not fully understood. A strong consensus is emerging around the idea that the stationary thermodynamic entropy is the von Neumann entanglement entropy of a large subsystem embedded in an infinite system. Also motivated by cold-atom experiments, here we consider the generalization to R\'enyi entropies. We develop a new technique to calculate the diagonal R\'enyi entropy in the quench action formalism. In the spirit of the replica treatment for the entanglement entropy, the diagonal R\'enyi entropies are generalized free energies evaluated over a thermodynamic macrostate which depends on the R\'enyi index and, in particular, is not the same state describing von Neumann entropy. The technical reason for this perhaps surprising result is that the evaluation of the moments of the diagonal density matrix shifts the saddle point of the quench action. An interesting consequence is that different R\'enyi entropies encode information about different regions of the spectrum of the postquench Hamiltonian. Our approach provides a very simple proof of the long-standing issue that, for integrable systems, the diagonal entropy is half of the thermodynamic one and it allows us to generalize this result to the case of arbitrary R\'enyi entropy.
Entropy is a fundamental concept in equilibrium statistical mechanics, yet its origin in the nonequilibrium dynamics of isolated quantum systems is not fully understood. A strong consensus is emerging around the idea that the stationary thermodynamic entropy is the von Neumann entanglement entropy of a large subsystem embedded in an infinite system. Also motivated by cold-atom experiments, here we consider the generalization to R'enyi entropies. We develop a new technique to calculate the diagonal R'enyi entropy in the quench action formalism. In the spirit of the replica treatment for the entanglement entropy, the diagonal R'enyi entropies are generalized free energies evaluated over a thermodynamic macrostate which depends on the R'enyi index and, in particular, is not the same state describing von Neumann entropy. The technical reason for this perhaps surprising result is that the evaluation of the moments of the diagonal density matrix shifts the saddle point of the quench action. An interesting consequence is that different R'enyi entropies encode information about different regions of the spectrum of the postquench Hamiltonian. Our approach provides a very simple proof of the long-standing issue that, for integrable systems, the diagonal entropy is half of the thermodynamic one and it allows us to generalize this result to the case of arbitrary R'enyi entropy.
Abstract We study the Rényi entropies in the spin- 1 / 2 anisotropic Heisenberg chain after a quantum quench starting from the Néel state. The quench action method allows us to obtain the stationary Rényi entropies for arbitrary values of the index α as generalised free energies evaluated over a calculable thermodynamic macrostate depending on α . We work out this macrostate for several values of α and of the anisotropy Δ by solving the thermodynamic Bethe ansatz equations. By varying α different regions of the Hamiltonian spectrum are accessed. The two extremes are α → ∞ for which the thermodynamic macrostate is either the ground state or a low-lying excited state (depending on Δ) and α = 0 when the macrostate is the infinite temperature state. The Rényi entropies are easily obtained from the macrostate as function of α and a few interesting limits are analytically characterised. We provide robust numerical evidence to confirm our results using exact diagonalisation and a stochastic numerical implementation of Bethe ansatz. Finally, using tDMRG we calculate the time evolution of the Rényi entanglement entropies. For large subsystems and for any α , their density turns out to be compatible with that of the thermodynamic Rényi entropies.
We study the behaviour of Renyi entropies in a generic thermodynamic macrostate of an integrable model. In the standard quench action approach to quench dynamics, the Renyi entropies may be derived from the overlaps of the initial state with Bethe eigenstates. These overlaps fix the driving term in the thermodynamic Bethe ansatz (TBA) formalism. We show that this driving term can be also reconstructed starting from the macrostate's particle densities. We then compute explicitly the stationary Renyi entropies after the quench from the dimer and the tilted Neel state in XXZ spin chains. For the former state we employ the overlap TBA approach, while for the latter we reconstruct the driving terms from the macrostate. We discuss in full detail the limits that can be analytically handled and we use numerical simulations to check our results against the large time limit of the entanglement entropies.
Abstract We study the behaviour of Rényi entropies in a generic thermodynamic macrostate of an integrable model. In the standard quench action approach to quench dynamics, the Rényi entropies may be derived from the overlaps of the initial state with Bethe eigenstates. These overlaps fix the driving term in the thermodynamic Bethe ansatz (TBA) formalism. We show that this driving term can be also reconstructed starting from the macrostate’s particle densities. We then compute explicitly the stationary Rényi entropies after the quench from the dimer and the tilted Néel state in XXZ spin chains. For the former state we employ the overlap TBA approach, while for the latter we reconstruct the driving terms from the macrostate. We discuss in full detail the limits that can be analytically handled and we use numerical simulations to check our results against the large time limit of the entanglement entropies.
Abstract We develop a general approach to compute the symmetry-resolved Rényi and von Neumann entanglement entropies (SREE) of thermodynamic macrostates in interacting integrable systems. Our method is based on a combination of the thermodynamic Bethe ansatz and the Gärtner–Ellis theorem from large deviation theory. We derive an explicit simple formula for the von Neumann SREE, which we show to coincide with the thermodynamic Yang–Yang entropy of an effective macrostate determined by the charge sector. Focusing on the XXZ Heisenberg spin chain, we test our result against iTEBD calculations for thermal states, finding good agreement. As an application, we provide analytic predictions for the asymptotic value of the SREE following a quantum quench.
Using arguments built on ergodicity, we derive an analytical expression for the Renyi entanglement entropies which, we conjecture, applies to the finite-energy density eigenstates of chaotic many-body Hamiltonians. The expression is a universal function of the density of states and is valid even when the subsystem is a finite fraction of the total system-a regime in which the reduced density matrix is not thermal. We find that in the thermodynamic limit, only the von Neumann entropy density is independent of the subsystem to the total system ratio V_{A}/V, while the Renyi entropy densities depend nonlinearly on V_{A}/V. Surprisingly, Renyi entropies S_{n} for n>1 are convex functions of the subsystem size, with a volume law coefficient that depends on V_{A}/V, and exceeds that of a thermal mixed state at the same energy density. We provide two different arguments to support our results: the first one relies on a many-body version of Berry's formula for chaotic quantum-mechanical systems, and is closely related to the eigenstate thermalization hypothesis. The second argument relies on the assumption that for a fixed energy in a subsystem, all states in its complement allowed by the energy conservation are equally likely. We perform an exact diagonalization study on quantum spin-chain Hamiltonians to test our analytical predictions.
In this Thesis I present the results of four research projects I carried out during my Ph.D. training at SISSA. The broad topic of my research was the mechanism of relaxation in isolated, integrable quantum systems. Roughly speaking, integrable quantum systems are one-dimensional lattice models or field theories that feature an extensively large number of conservation laws beyond energy and momentum conservation, which greatly influence their dynamical behaviour. Integrable systems have enjoyed a considerable interest in the last decade, mainly due to the development of ultracold atomic techniques. These techniques made it possible to implement experimentally a wealth of one-dimensional quantum models whose significance had been purely theoretical before. The theoretical interest in integrable systems stems in part from the fact that in many aspects they are exactly solvable. Exact analytical formulas are available for their spectrum and their equilibrium thermodynamic properties, creating an incentive to study analytically their non-equilibrium behaviour as well. Exact results for the dynamics of integrable systems may serve as a benchmark for confirming theories in statistical mechanics or testing approximative methods. To this day, computing analytically the entire time evolution of observables in a genuinely interacting integrable system is prohibitively difficult in general. In this Thesis, I concentrate on the description of the system at long times after it is brought out of equilibrium. Homogeneous, translationally invariant integrable systems are expected to relax to a stationary state, in which local observables are described by a generalization of the Gibbs ensemble that takes into consideration all the conservation laws beyond the energy. Two of the projects presented here concern this relaxed state. In Chapter 2, a method is presented for computing the Rényi entropies of the Generalized Gibbs Ensemble describing the relaxed state in the integrable Heisen- berg XXZ spin chain. Rényi entropies provide information about the relaxed state and its entanglement properties, going beyond the ordinary von Neumann thermo- dynamic entropy. The method presented there is very general; in principle it works in any integrable model and any translationally invariant thermodynamic macrostate. An alternative to the GGE in describing relaxed states after a quantum quench is the Quench Action approach. This approach is based on the explicit analytical knowledge of the overlaps between the initial state and the Hamiltonian eigenstates. In Chapter 3, this approach is generalized to a spin-1 integrable spin chain, the SU(3) invariant Lai–Sutherland model. The evolution of entanglement is then studied using the quasiparticle picture, which features ballistically propagating classical particle pairs carrying entanglement. If the requirement of translational invariance of the initial state is relaxed, various transport phenomena arise in integrable systems. The evolution of slowly varying, smoothly distributed states is described by a generalization of hydrodynamics, which takes into consideration the infinitely many local conservation laws. At the scale of non-diffusive, ballistic transport, the main feature of this generalized hydrodynamics (GHD) is the existence of continuity equations for the occupation number of each mode. Further two projects are presented in this Thesis which are concerned with GHD.I n Chapter 4, a numerical method is presented for computing the evolution of entanglement entropies within GHD. This method combines the quasiparticle picture of entanglement with a molecular dynamics simulation of GHD, the so-called fleagas. Numerical evidence is shown that this method reproduces previously known analytical results in GHD. Its versatility is demonstrated by computing another entanglement-related quantity, the mutual information, for which no analytical results are available. Finally, in Chapter 5, the GHD framework is generalized to a system of spin-1/2 fermions with repulsive contact interaction, the Yang–Gaudin model. It is shown that at low temperatures, the dynamics of the system is described by the superposition of two uncoupled, inhomogeneous conformal field theories. This fact can be interpreted as a separation of the spin and the charge degrees of freedom of the system.