Quench action and Renyi entropies in integrable systems
Abstract
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 coldatom experiments, here we consider the generalisation to Rényi entropies. We develop a new technique to calculate the diagonal Rényi entropy in the quench action formalism. In the spirit of the replica treatment for the entanglement entropy, the diagonal Rényi entropies are generalised free energies evaluated over a thermodynamic macrostate which depends on the Rényi index and, in particular, it is not the same describing the von Neumann entropy. The technical reason for this, maybe 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ényi entropies encode information about different regions of the spectrum of the postquench Hamiltonian. Our approach provides a very simple proof of the longstanding issue that, for integrable systems, the diagonal entropy is half of the thermodynamic one and it allows us to generalise this result to the case of arbitrary Rényi entropy.
Introduction.–An intriguing question for an outofequilibrium isolated system evolving from a pure state is to understand microscopically the thermodynamic entropy of the emergent statistical ensemble. Indeed, it has been undoubtedly established that for large times the local properties of a quenched system are correctly captured by a statistical ensemble which is either the Gibbs or generalised Gibbs depending on whether the system is generic (2); (3); (4); (5); (6); (7); (8) or integrable (9); (10); (11); (12); (13); (14); (15); (16); (17); (18); (19); (20); (21); (22); (23); (24); (25); (26); (27); (28); (29); (30); (31); (32); (33); (34); (35); (36) respectively. However, since the entire system is always in a zeroentropy pure state because of the unitary time evolution, it is not obvious whether the entropy of the statistical ensemble has any physical meaning. This apparent paradox has been nowadays clarified and it is widely accepted that the thermodynamic entropy is (37) the stationary value of the entanglement of a large subsystem embedded in an infinite system (38); (40); (39); (41); (42); (43); (44); (45); (46); (47). This entanglement is quantified by the von Neumann entropy of the reduced density matrix of the subsystem . Very recently, the equivalence between stationary entanglement entropy and thermodynamic entropy has been exploited to provide analytic exact predictions for the entire time dependence of the entanglement entropy in integrable models (46) as carefully tested against numerical simulations in the XXZ spinchain (46).
An alternative and natural definition of entropy after a quench is the diagonal entropy (48), which is the von Neumann entropy of the diagonal ensemble
(1) 
Here the states are the eigenstates of the postquench Hamiltonian and the coefficients are the overlaps with the initial state . The great advantage of the diagonal entropy is that can be computed from the initial state, without solving the complicated quench dynamics. A lot of effort has been devoted to understand the relation between the diagonal entropy and the stationary subsystem entanglement entropy (which, as already stressed, is the entropy of the stationary ensemble). It has been suggested that for generic systems relaxing to a Gibbs ensemble the diagonal entropy is the same as the thermodynamic one (48), as tested also in numerical simulations (49); (40). Contrarily, for quenches to integrable models it has been found numerically and analytically that the diagonal entropy is half of the thermodynamic one (40); (50); (51); (52); (41); (42); (53). A solid physical argument to explain this factor has been provided by Gurarie (50), but so far there was no proof, neither a complete understanding of its generality.
For a generic density matrix , the Rényi entropies
(2) 
represent a more general class of entropies generalising the von Neumann one. Here is the index of the Rényi entropy and it is an arbitrary positive real number. In the limit , (2) gives back the von Neumann entropy. Although Rényi entropies for generic do not have some important properties of the von Neumann one (the most physical one being the strong subadditivity), they are attracting a lot of attention, especially in connection with entanglement. Indeed, Rényi entanglement entropies (i.e. (2) with ) are entanglement monotones (54), i.e. quantities that do not increase under local operations and classical communication which is the technical requirement for any good measure of entanglement. Furthermore, for integer , they are the core of the replica approach to the entanglement entropy (55). While the replica method was introduced mainly as a theoretical trick (55), it became a useful tool for the experimental measurement of entanglement: Rényi entanglement entropies (for ) have been measured experimentally with cold atoms, both in equilibrium (56) and after a quantum quench (44).
The main goal of this paper is to initiate the investigation of Rényi entropies in interacting integrable models also, but not only, to provide predictions testable in experiments. To this aim, we adapt the Thermodynamic Bethe Ansatz (TBA) techniques for quantum quenches (overlap TBA or Quench Action method (57); (58)) and we derive analytically the diagonal Rényi entropies (i.e. (2) with ). In order to make this paper as readable as possible, we report only the general framework for calculating the diagonal entropies and we discuss in details the main physical consequences of our findings. We postpone a detailed analytical and numerical technical analysis of the TBA solutions for specific integrable models to forthcoming publications (59).
The Thermodynamic Bethe ansatz (TBA).–In Bethe ansatz solvable models, the eigenstates are in correspondence with a set of pseudomomenta , usually called rapidities. These rapidities satisfy a set of nonlinear coupled equation known as Bethe equations. In the thermodynamic limit (60) the rapidities form a continuum and the eigenstates are characterised by the densities of particles and of holes (holes are unoccupied Bethe modes that in an interacting model are not simply related to ). With bold symbols and we denoted the full set of all stable Bethe quasiparticles, i.e. , where depends on the considered model (e.g. there is only one species for the repulsive LiebLiniger model, but infinite for the gapped XXZ chain). In a given eigenstate, the densities and are related by the continuum limit of the Bethe equation. A set of densities identifies a thermodynamic macrostate that corresponds to an exponentially large number of microscopic eigenstates. The total number of microstates with the same and is , with the thermodynamic YangYang entropy of the macrostate (61)
(3) 
where the total density is . For systems in thermal tequilibrium, is the thermal entropy (62).
In the thermodynamic limit, a sum over eigenstates is replaced by a functional integral over the rapidity densities as
(4) 
where the factor counts the exponentially large number of microscopic eigenstates leading to the same densities.
The quench action.–The out of equilibrium dynamics of integrable systems can be cast in the TBA language with the quench action method (57); (58). Although this approach can be used to study the time evolution (as done in a few cases (63); (64); (65)), it mainly provides a calculable and manageable representation of the stationary state (or equivalently of the diagonal ensemble), as nowadays explicitly worked out for many integrable models (66); (67); (68); (69); (70); (71); (72); (73); (77); (74); (75); (76). The stationary value of a local observable , i.e. its value in the diagonal ensemble, is
(5) 
According to the recipe (4), in the thermodynamic limit this becomes
(6) 
where we defined the thermodynamic limit of the logarithm of the overlaps as (60)
(7) 
This functional integral can be evaluated in the limit using the saddle point approximation, since both the and are generally estensive. One has to minimise the functional appearing in the exponential in (6), i.e.
(8) 
which has been termed quench action. The minimisation is carried out solving
(9) 
under the constraint that the thermodynamic Bethe equations hold. Notice that appearing in (6) and (8) counts only those microstates with a non zerooverlap with the initial state (i.e. those contributing to the sum (5) before taking the thermodynamic limit). This is only a small fraction of the total ones: in all quenches solved so far only parity invariant Bethe states (i.e. those containing only pairs of rapidities with opposite sign, i.e., such that ) have nonzero overlap (78); (79); (77); (80); (81); (82); (83); (84); (85); (86). This fact can be effectively taken into account either by performing the integral defining (3) only on positive ’s or multiplying by . We will adopt the latter choice.
Finally for the expectation values of observables (6), one obtains
(10) 
under the reasonable assumption that the considered operator does not shift the saddle point (as obvious for local ones). As a particular case, the identity operator in (5) provides the normalisation of the initial state , implying that at the saddle point the quench action vanishes, i.e.
(11) 
Overlap TBA for the diagonal entropies.– So far we have collected all the needed ingredients to study the diagonal entropy with the quench action. Taking the thermodynamic limit (4) of the definition (2) with in (1), we simply have
(12) 
where we introduced the factor in front of to explicitly take into account that only parityinvariant eigenstates have nonzero overlap with most initial states, as discussed above. The most important aspect of (12) is that the Rényi index appears in the exponential term and so it shifts the saddle points from the one in (9).
The modified quench action which should be considered in this case explicitly depends on the Rényi index as
(13) 
and the quench action saddlepoint equation is
(14) 
The Rényi diagonal entropies are then the saddle point expectation of this quench action
(15) 
can be thought of as a generalised free energy, with playing the role of the Hamiltonian. Using the property that as in (11), it is convenient to rewrite (15) as
(16) 
which is a form that closely resembles the replica approach to the entanglement entropy (55).
The set of coupled equations (14) can be solved, at least numerically, by standard methods. This analysis will be presented elsewhere (59), while here we only discuss some general properties of the diagonal Rényi entropies following from (14) and (15).
Proof that the diagonal entropy is half of the YangYang entropy.– We now consider the limit
(17) 
The normalisation (11), implying , ensures the limit to be finite. We then need to expand at the first order in the saddle point quench action as
(18) 
We have
(19) 
in which the first term in the rhs is zero because of the saddlepoint condition (14) while the second term is simply (cf. (13)). Hence we have
(20) 
i.e.
(21) 
where in the last equality we used the normalisation conditions (11). This proves that the diagonal entropy is half of the thermodynamic one. The sole assumption to derive this general result is that only parity invariant states have nonzero overlaps. While this is a technical assumption, it is exactly on the same line of thoughts of the physical interpretation of Gurarie (50). If for some quenches (still unknown in integrable models) all states would have nonzero overlaps, then the above proof would lead to the equality between the diagonal and the thermodynamic entropy. If one imagine of generalising the quench action approach to nonintegrable model (as speculated in (58)) then the above argument would immediately lead to the equality of the two entropies, as numerically observed (49).
Indeed, our formalism can be used in general to show that Rényi diagonal entropies are half of the thermodynamic Rényi entropies for arbitrary . Anyhow, the derivation is slightly more technical and so we report it in the appendix.
Free fermionic systems.– For translational invariant systems of free fermions (or models mappable to), the diagonal Rényi entropies after a quench can be written in terms of the momentum occupation numbers as
(22) 
which simply tells that each mode is populated with probability and empty with probability . Again, assuming parity invariance, the integral is only on positive . In TBA language, the rapidities are the real momenta ; there is only one root density and the hole density is .
Thus, Eq. (22) implies that, for free fermions, the Rényi entropies for arbitrary can be written in terms of the root density at without the need of generalised TBA equations. It is instructive to recover (22) from the generalised TBA saddle point (15) and show that, in this particular case, it is possible to write in terms of . In the following we will work with to get rid of factors ’s and often omit the dependence.
The YangYang entropy for free fermions is
(23) 
because . (For lattice free fermions, the integral runs between and , but this does not influence any of the following results.) Exploiting the parity invariance, the driving term can be always written as
(24) 
for some calculable depending on the initial state (e.g. for the quench from zero to infinite interaction in the LiebLiner gas (78) with being the density). The solution of the generalised quench action saddle point equation (14) is straightforwardly worked out
(25) 
The mode occupation is , i.e. .
In terms of the generalised saddlepoint free energy is, after simple algebra,
(26) 
One can use the relation between and to write as function of as
(27) 
The Renyi entropy is straightforwardly obtained from (16)
(28) 
which is the same as (22). Notice that to show this equality in a simple and general way we had to add to the saddlepoint generalised action the term which is zero, but written in terms of an integral of a nonzero integrand.
The same derivation above also holds for free bosons with very minor modifications.
A numerical example for an interacting model.– It is far beyond the aim of this paper to present detailed calculations of the Rényi diagonal entropies for interacting integrable models which require a large amount of technicalities. However, it is instructive to show at least one example of the comparison between the quench action approach and the numerical computation for finite systems. We focus on the quench in the spin XXZ chain starting from the Néel state which has been studied intensively in the literature (69); (71); (87). The exact diagonalisation data for the diagonal entropies are reported in Fig. 1 showing for and . Finitesize effects are visible for both and . The star symbols denote the diagonal entropy densities in the thermodynamic limit, which are obtained using (15). The dashdotted lines are linear fits to with fixed to the TBA prediction. The data are clearly compatible with the quench action prediction for all values of although large oscillations in affect the numerical data.
Conclusions.– In this paper we developed a Bethe ansatz framework to calculate the diagonal Rényi entropies after a quantum quench in integrable models. The diagonal Rényi entropies are generalised free energies evaluated over a thermodynamic macrostate which depends on the Rényi index and it is not the same describing the von Neumann entropy and the local observables. We checked that this prescription reproduces known results for free models and we tested numerically its validity for interacting integrable models. An intriguing consequence is that the steady state contains information about different regions of the spectrum of the chain, which can be accessed by varying . This is similar, in spirit, to the observation (88) that a single eigenstate of a generic (nonintegrable) Hamiltonian at finite energy density contains information about the full spectrum of the Hamiltonian. A byproduct of our approach is a very simple proof of the long standing issue that the diagonal entropy is half of the thermodynamic one in the steady state. This is a direct consequence of the parity invariance of the Bethe eigenstates with nonzero overlap with the initial state.
This Letter opens several new questions. First, it would be important to understand whether the Rényi entropies can be rewritten in terms of the saddle point root densities describing local observables. This does not look neither simple nor reasonable from our approach, but, physically, the thermodynamic Rényi entropy is the same as the entanglement one for a large subsystem embedded in an infinite system. Thus it is, in principle, expected to be describable by the same saddlepoint of the von Neumann entropy. Having the Rényi entropies written in terms of the standard root densities would allow to reconstruct the entire time dependence of the entanglement Rényi entropies exploiting the quasiparticle picture (38), similarly to what done for the von Neumann entropy (46).
Another issue is the role played by the parity invariance of the relevant postquench eigenstates. Up to now this has been seen only as a technical remark (see however (89)), but we showed that it affects the ratio between diagonal and thermodynamic entropies. Is it possible to design a solvable quench where nonparity invariant states have non vanishing overlaps so that the ratio between the entropies equals one? (the arguments in (89) suggests this to be unlikely). Alternatively, are there quenches in which only a small subclass of parity invariant eigenstates have nonzero overlap making the ratio between the two entropies even smaller than ?
Acknowledgments.– We are grateful to Bruno Bertini for suggesting the proof reported in the appendix. PC acknowledges support from the ERC under the Starting Grant No 279391 EDEQS and VA the European Union’s Horizon 2020 under the Marie SklodowoskaCurie grant agreement No 702612 OEMBS.
Note added.– After the completion of this work, one of us and coworkers found a set of initial states for freelike theories in which the eigenstates with nonzero overlap are different from parity invariant states (90). The analysis of the corresponding entropies is not trivial and still under examination.
Appendix A Rényi diagonal entropies are half of thermodynamic ones
In this appendix we prove that Rényi diagonal entropies are half of thermodynamic one for arbitrary index , generalizing the result for presented in the main text. However, this proof turn out to be more cumbersome and technical than the one for the von Neumann entropy, mainly because the thermodynamic Rényi entropies must be derived from the GGE since there is no analogous of the YangYang entropy. In particular, the proof requires the manipulation of the explicit TBA equations that we avoided in the main text.
In TBA language, the saddle point equation (9) is a set of coupled integral equations for the root and hole densities and for each string of length . These equations are written in a more compact form in terms of the ratios , assuming the standard form (62); (58)
(29) 
where the symbol denotes the convolution and we introduced the driving term from
(30) 
In Eq. (29) the functions encode the interaction in the thermodynamic limit and are simply read from the Bethe equations of the specific model (see e.g. (62)). Consequently the saddle point equation for general in Eq. (14) are simply
(31) 
The generalized Gibbs ensemble describing local observables after a quench to an integrable model has the form
(32) 
where the operators ’s form a complete set of local and quasilocal integrals of motion and are fixed by the condition that the charges assume the same value in the initial state and in the GGE, i.e. .
The expectation values of local operator in the GGE are obtained by generalized TBA (91). Analogously to (6) for the diagonal ensemble, we write the expectation value of a local operator as a path integral
(33) 
where we used the fact that Bethe states are eigenstates of the conserved charges with eigenvalues , i.e.
(34) 
The resulting saddle point equations are easily read from the above path integral (91)
(35) 
where we introduced the functions as
(36) 
These functions are the representation of the local and quasilocal integrals of motion in terms of the rapidities and are generically known but can be cumbersome. For example in the case of the repulsive LiebLiniger model where there are no strings (i.e. there is not the sum over above), these functions are simply , i.e. .
At this point, since the expectation values of all the local observables in the diagonal ensemble and in the GGE must be equal, these must have the same saddle point root distribution which implies Eqs. (29) and (35) to be equal. This equality leads to the identification
(37) 
for where is defined. This identification implies
(38) 
where we used the fact that the integral in (30) is on the positive axis, while the one in (36) is on the entire real line. Notice that those functions which are odd in do not contribute to the GGE because the considered initial state is even and the corresponding ’s are vanishing. This consideration ensures the validity of (38).
We are finally ready to obtain thermodynamic Rényi entropies. These are given by
(39) 
Thus, the saddle point equation are the same as in (35) with the replacement for all . Because of the identification (37), these generalised TBA for the GGE equations are the same as those for the diagonal ensemble (31). This equivalence could have been expected, but it is far from trivial because Rényi entropies are global and not local quantities (and indeed we will soon find a factor 2 between the two). At this point, the GGE Rényi entropy is read off from the saddlepoint of the above path integral, in analogy to the diagonal one (15), as
(40) 
which using (38) is exactly twice the diagonal Rényi entropy (15).
References
References
 J. von Neumann, Beweis des Ergodensatzes und des HTheorems, Z Phys. 57, 30 (1929).
 R. V. Jensen and R. Shankar, Statistical Behavior in Deterministic Quantum Systems with Few Degrees of Freedom, Phys. Rev. Lett. 54, 1879 (1985).
 J. M. Deutsch, Quantum statistical mechanics in a closed system, Phys. Rev. A 43, 2046 (1991)..
 M. Srednicki, Chaos and quantum thermalization, Phys. Rev. E 50, 888 (1994).
 M. Rigol, V. Dunjko, and M. Olshanii, Thermalization and its mechanism for generic isolated quantum systems, Nature 452, 854 (2008).
 M. Rigol and M. Srednicki, Alternatives to Eigenstate Thermalization, Phys. Rev. Lett. 108, 110601 (2012)..
 L. D’Alessio, Y. Kafri, A. Polkovnikov, and M. Rigol, From Quantum Chaos and Eigenstate Thermalization to Statistical Mechanics and Thermodynamics, Adv. Phys. 65, 239 (2016).
 M. Rigol, V. Dunjko, V. Yurovsky, and M. Olshanii, Relaxation in a Completely Integrable ManyBody Quantum System: An Ab Initio Study of the Dynamics of the Highly Excited States of 1D Lattice HardCore Bosons. Phys. Rev. Lett. 98, 050405 (2007).
 M. A. Cazalilla, Effect of Suddenly Turning on Interactions in the Luttinger Model, Phys. Rev. Lett. 97, 156403 (2006).
 T. Barthel and U. Schollwöck, Dephasing and the Steady State in Quantum ManyParticle Systems. Phys. Rev. Lett. 100, 100601 (2008).

M. Cramer, C. M. Dawson, J. Eisert, and T. J. Osborne,
Exact Relaxation in a Class of Nonequilibrium Quantum Lattice Systems,
Phys. Rev. Lett. 100, 030602 (2008);
M. Cramer and J. Eisert, A quantum central limit theorem for nonequilibrium systems: exact local relaxation of correlated states, New J. Phys. 12, 055020 (2010).  M. Fagotti and P. Calabrese, Evolution of entanglement entropy following a quantum quench: Analytic results for the XY chain in a transverse magnetic field. Phys. Rev. A 78, 010306 (2008).
 M. A. Cazalilla, A. Iucci, and M.C. Chung, Thermalization and quantum correlations in exactly solvable models, Phys. Rev. E 85, 011133 (2012).

P. Calabrese, F. H. L. Essler, and M. Fagotti, Quantum Quench in the TransverseField Ising Chain,
Phys. Rev. Lett. 106, 227203 (2011);
P. Calabrese, F. H. L. Essler, and M. Fagotti, Quantum quench in the transverse field Ising chain: I. Time evolution of order parameter correlators, J. Stat. Mech. (2012) P07016;
P. Calabrese, F. H. L. Essler, and M. Fagotti, Quantum quenches in the transverse field Ising chain: II. Stationary state properties, J. Stat. Mech. (2012) P07022. 
D. Fioretto and G. Mussardo, Quantum quenches in integrable field theories,
New J. Phys. 12, 055015 (2010);
S. Sotiriadis, D. Fioretto, and G. Mussardo, ZamolodchikovFaddeev algebra and quantum quenches in integrable field theories, J. Stat. Mech. (2012) P02017. 
M. Collura, S. Sotiriadis, and P. Calabrese, Equilibration of a TonksGirardeau Gas Following a Trap Release,
Phys. Rev. Lett. 110, 245301 (2013);
M. Collura, S. Sotiriadis, and P. Calabrese, Quench dynamics of a TonksGirardeau gas released from a harmonic trap, J. Stat. Mech. (2013) P09025.  M. Fagotti and F. H. L. Essler, Stationary behaviour of observables after a quantum quench in the spin Heisenberg chain, J. Stat. Mech. (2013) P07012.
 B. Pozsgay, The generalized Gibbs ensemble for Heisenberg spin chains, J. Stat. Mech. (2013) P07003.
 M. Fagotti and F. H. L. Essler, Reduced Density Matrix after a Quantum Quench, Phys. Rev. B 87, 245107 (2013).
 M. Kormos, M. Collura, and P. Calabrese, Analytic results for a quantum quench from free to hardcore onedimensional bosons, Phys. Rev. A 89, 013609 (2014).
 S. Sotiriadis and P. Calabrese, Validity of the GGE for quantum quenches from interacting to noninteracting models, J. Stat. Mech. (2014) P07024.
 M. Fagotti, M. Collura, F. H. L. Essler, and P. Calabrese, Relaxation after quantum quenches in the spin12 Heisenberg XXZ chain, Phys. Rev. B 89, 125101 (2014).

E. Ilieveski, J. De Nardis, B. Wouters, J.S. Caux, F. H. L. Essler, and T. Prosen,
Complete Generalized Gibbs Ensembles in an Interacting Theory,
Phys. Rev. Lett. 115, 157201 (2015);
E. Ilievski, E. Quinn, J. D. Nardis, and M. Brockmann, Stringcharge duality in integrable lattice models, J. Stat. Mech. (2016) 063101.  V. Alba, Simulating the Generalized Gibbs Ensemble (GGE): a Hilbert space Monte Carlo approach. arXiv:1507.06994.
 T. Langen, S. Erne, R. Geiger, B. Rauer, T. Schweigler, M. Kuhnert, W. Rohringer, I. E. Mazets, T. Gasenzer, and J. Schmiedmayer, Experimental observation of a generalized Gibbs ensemble, Science 348, 207 (2015).

F. H. L. Essler, G. Mussardo, and M. Panfil, Generalized Gibbs ensembles for quantum field theories,
Phys. Rev. A 91, 051602 (2015);
F. H. L. Essler, G. Mussardo, and M. Panfil, On Truncated Generalized Gibbs Ensembles in the Ising Field Theory, J. Stat. Mech. (2017) 013103.  J. Cardy, Quantum quenches to a critical point in one dimension: some further results, J. Stat. Mech. (2016) 023103.
 S. Sotiriadis, Memorypreserving equilibration after a quantum quench in a 1d critical model, Phys. Rev. A 94, 031605 (2016).
 A. Bastianello and S. Sotiriadis, Quasi locality of the GGE in interactingtofree quenches in relativistic field theories, J. Stat. Mech. (2017) 023105..
 E. Vernier and A. Cortés Cubero, Quasilocal charges and the complete GGE for field theories with non diagonal scattering, J. Stat. Mech. (2017) 23101..
 B. Pozsgay, E. Vernier, and M. A. Werner, On Generalized Gibbs Ensembles with an infinite set of conserved charges, arXiv:1703.09516 (2017).
 L. Vidmar and M. Rigol, Generalized Gibbs ensemble in integrable lattice models, J. Stat. Mech. (2016) 064007.
 C. Gogolin and J. Eisert, Equilibration, thermalisation, and the emergence of statistical mechanics in closed quantum systems, Rep. Prog. Phys. 79, 056001 (2016).
 P. Calabrese, F. H. L. Essler, and G. Mussardo, Introduction to “Quantum Integrability in Out of Equilibrium Systems”, J. Stat. Mech. (2016) P064001.
 F. H. L. Essler and M. Fagotti, Quench dynamics and relaxation in isolated integrable quantum spin chains, J. Stat. Mech. (2016) 064002.
 Hereafter, when stating that two entropies are the same, we mean that they have the same extensive behaviour (i.e. the same density), but they can have different subleading terms (and in most cases they do).
 P. Calabrese and J. Cardy, Evolution of Entanglement Entropy in OneDimensional Systems, J. Stat. Mech. (2005) P04010.
 J. M. Deutsch, H. Li, and A. Sharma, Microscopic origin of thermodynamic entropy in isolated systems, Phys. Rev. E 87, 042135 (2013).
 L. F. Santos, A. Polkovnikov, and M. Rigol, Weak and strong typicality in quantum systems, Phys. Rev. E 86, 010102 (2012).
 M. Collura, M. Kormos, and P. Calabrese, Stationary entropies following an interaction quench in Bose gas, J. Stat. Mech. P01009 (2014).

M. Kormos, L. Bucciantini, and P. Calabrese, Stationary entropies after a quench from excited states in the Ising chain,
EPL 107, 40002 (2014);
L. Bucciantini, M. Kormos, and P. Calabrese, Quantum quenches from excited states in the Ising chain, J. Phys. A 47, 175002 (2014).  W. Beugeling, A. Andreanov, and M. Haque, Global characteristics of all eigenstates of local manybody Hamiltonians: participation ratio and entanglement entropy, J. Stat. Mech. (2015) P02002.
 A. M. Kaufman, M. E. Tai, A. Lukin, M. Rispoli, R. Schittko, P. M. Preiss, and M. Greiner, Quantum thermalization through entanglement in an isolated manybody system, Science 353, 794 (2016).
 P. Calabrese and J. Cardy, Quantum quenches in 1+1 dimensional conformal field theories, J. Stat. Mech. (2016) 064003.
 V. Alba and P. Calabrese, Entanglement and thermodynamics after a quantum quench in integrable systems, PNAS 114, 7947 (2017).
 J. Dubail, Entanglement scaling of operators: a conformal field theory approach, with a glimpse of simulability of longtime dynamics in 1+1d, J. Phys. A 50, 234001 (2017).
 A. Polkovnikov, Microscopic diagonal entropy and its connection to basic thermodynamic relations, Ann. Phys. 326, 486 (2011).
 L. F. Santos, A. Polkovnikov, and M. Rigol, Entropy of isolated quantum systems after a quench, Phys. Rev. Lett. 107, 040601 (2011).
 V. Gurarie, Global large time dynamics and the generalized Gibbs ensemble, J. Stat. Mech. (2013) P02014.
 M. Fagotti, Finitesize corrections vs. relaxation after a sudden quench, Phys. Rev. B 87, 165106 (2013).
 B. Dóra, Escort distribution function of work done and diagonal entropies in quenched Luttinger liquids, Phys. Rev. B 90, 245132 (2014).
 L. Piroli, E. Vernier, P. Calabrese, and M. Rigol, Correlations and diagonal entropies after quantum quenches in chains, Phys. Rev. B 95, 054308 (2017).
 G. Vidal, Entanglement monotones, J. Mod. Opt. 47, 355 (2000).

P. Calabrese and J. Cardy, Entanglement entropy and quantum field theory,
J. Stat. Mech. P06002 (2004);
P. Calabrese and J. Cardy, Entanglement entropy and conformal field theory, J. Phys. A 42, 504005 (2009).  R. Islam, R. Ma, P. M. Preiss, M. E. Tai, A. Lukin, M. Rispoli, and M. Greiner, Measuring entanglement entropy in a quantum manybody system, Nature 528, 77 (2015).
 J.S. Caux and F. H. L. Essler, Time evolution of local observables after quenching to an integrable model, Phys. Rev. Lett. 110, 257203 (2013).
 J.S. Caux, The Quench Action, J. Stat. Mech. (2016) 064006.
 V. Alba and P. Calabrese, in preparation.
 The thermodynamic limit is taken with the number of particles and the length going to infinity at fixed density .
 C. N. Yang and C. P. Yang, Thermodynamics of a One Dimensional System of Bosons with Repulsive Delta Function Interaction, J. Math. Phys. 10, 1115 (1969).
 M. Takahashi, Thermodynamics of onedimensional solvable models, Cambridge University Press, Cambridge, 1999.
 J. De Nardis and J.S. Caux, Analytical expression for a postquench time evolution of the onebody density matrix of onedimensional hardcore bosons, J. Stat. Mech. (2014), P12012.
 J. De Nardis, L. Piroli, and J.S. Caux, Relaxation dynamics of local observables in integrable systems, J. Phys. A 48, 43FT01 (2015).
 L. Piroli and P. Calabrese, Exact dynamics following an interaction quench in a onedimensional anyonic gas, Phys. Rev. A 96, 023611 (2017).
 B. Bertini, D. Schuricht, and F. H. L. Essler, Quantum quench in the sineGordon model, J. Stat. Mech. (2014) P10035.
 B. Bertini, L. Piroli, and P. Calabrese, Quantum quenches in the sinhGordon model: steady state and onepoint correlation functions, J. Stat. Mech. (2016) 063102.
 B. Pozsgay, The dynamical free energy and the Loschmidt echo for a class of quantum quenches in the Heisenberg spin chain, J. Stat. Mech. (2013) P10028.

B. Wouters, J. De Nardis, M. Brockmann, D. Fioretto, M. Rigol, and J.S. Caux,
Quenching the Anisotropic Heisenberg Chain: Exact Solution and Generalized Gibbs Ensemble Predictions,
Phys. Rev. Lett. 113, 117202 (2014);
M. Brockmann, B. Wouters, D. Fioretto, J. D. Nardis, R. Vlijm, and J.S. Caux, Quench action approach for releasing the Néel state into the spin1/2 XXZ chain, J. Stat. Mech. (2014) P12009. 
B. Pozsgay, M. Mestyán, M. A. Werner, M. Kormos, G. Zaránd, and G. Takács,
Correlations after Quantum Quenches in the Spin Chain: Failure of the Generalized Gibbs Ensemble,
Phys. Rev. Lett. 113, 117203 (2014);
M. Mestyán, B. Pozsgay, G. Takács, and M. A. Werner, Quenching the XXZ spin chain: quench action approach versus generalized Gibbs ensemble, J. Stat. Mech. (2015) P04001.  V. Alba, and P. Calabrese, The quench action approach in finite integrable spin chains, J. Stat. Mech. (2016), 043105.
 L. Piroli, B. Pozsgay, and E. Vernier, From the quantum transfer matrix to the quench action: the Loschmidt echo in XXZ Heisenberg spin chains, J. Stat. Mech. (2017) 23106.
 M. Mestyan, B. Bertini, L. Piroli, and P. Calabrese, Exact solution for the quench dynamics of a nested integrable system, J. Stat. Mech. (2017) 083103.

L. Piroli, P. Calabrese, and F. H. L. Essler,
Multiparticle BoundState Formation following a Quantum Quench to the OneDimensional Bose Gas with Attractive Interactions,
Phys. Rev. Lett. 116, 070408 (2016);
L. Piroli, P. Calabrese, and F. H. L. Essler, Quantum quenches to the attractive onedimensional Bose gas: exact results, SciPost Phys. 1, 001 (2016).  L. Bucciantini, Stationary State After a Quench to the LiebLiniger from Rotating BECs, J Stat Phys 164, 621 (2016).

J. De Nardis, M. Panfil, A. Gambassi, L. F. Cugliandolo, R. Konik, and L. Foini,
Probing nonthermal density fluctuations in the onedimensional Bose gas,
arXiv:1704.06649;
J. De Nardis and M. Panfil, Exact correlations in the LiebLiniger model and detailed balance outofequilibrium, SciPost Phys. 1, 015 (2016).  J. De Nardis, B. Wouters, M. Brockmann, and J.S. Caux, Solution for an interaction quench in the LiebLiniger Bose gas, Phys. Rev. A 89, 033601 (2014).
 V. Gritsev, T. Rostunov, and E. Demler, Exact methods in the analysis of the nonequilibrium dynamics of integrable models: application to the study of correlation functions for nonequilibrium 1D Bose gas, J. Stat. Mech. (2010) P05012.

P. Le Doussal and P. Calabrese,
The KPZ equation with flat initial condition and the directed polymer with one free end
J. Stat. Mech. (2012) P06001;
P. Calabrese and P. Le Doussal, Interaction quench in a LiebLiniger model and the KPZ equation with flat initial conditions J. Stat. Mech. (2014) P05004.  B. Pozsgay, Overlaps between eigenstates of the XXZ spin1/2 chain and a class of simple product states, J. Stat. Mech. (2014) P06011.

M. Brockmann, J. D. Nardis, B. Wouters, and J.S. Caux, A Gaudinlike determinant for overlaps of Néel and XXZ Bethe states,
J. Phys. A 47, 145003 (2014);
M. Brockmann, Overlaps of qraised Néel states with XXZ Bethe states and their relation to the LiebLiniger Bose gas, J. Stat. Mech. (2014) P05006;
M. Brockmann, J. De Nardis, B. Wouters, and J.S. Caux, NéelXXZ state overlaps: odd particle numbers and LiebLiniger scaling limit, J. Phys. A 47, 345003 (2014).  L. Piroli and P. Calabrese, Recursive formulas for the overlaps between Bethe states and product states in XXZ Heisenberg chains, J. Phys. A 47, 385003 (2014).

M. de Leeuw, C. Kristjansen, and K. Zarembo,
Onepoint Functions in Defect CFT and Integrability,
JHEP 08 (2015) 098;
I. BuhlMortensen, M. de Leeuw, C. Kristjansen, and K. Zarembo, Onepoint Functions in AdS/dCFT from Matrix Product States, JHEP 02 (2016) 052;
O. Foda and K. Zarembo, Overlaps of partial Néel states and Bethe states, J. Stat. Mech. (2016) 23107.  M. de Leeuw, C. Kristjansen, and S. Mori, AdS/dCFT onepoint functions of the SU(3) sector, Phys. Lett. B 763, 197 (2016).

D. X. Horváth, S. Sotiriadis, and G. Takács, Initial states in integrable quantum field theory quenches from an integral equation hierarchy,
Nucl. Phys. B 902, 508 (2016);
D. X. Horváth and G. Takács, Overlaps after quantum quenches in the sineGordon model, Phys. Lett. B 771, 539 (2017).  P. P. Mazza, J.M. Stéphan, E. Canovi, V. Alba, M. Brockmann, and M. Haque, Overlap distributions for quantum quenches in the anisotropic Heisenberg chain, J. Stat. Mech. (2016) 013104.
 L. Piroli, E. Vernier, and P. Calabrese, Exact steady states for quantum quenches in integrable Heisenberg spin chains, Phys. Rev. B 94, 054313 (2016).
 J. R. Garrison and T. Grover, Does a single eigenstate encode the full Hamiltonian?, arXiv:1503.00729.
 G. Delfino, Quantum quenches with integrable prequench dynamics, J. Phys. A 47 (2014) 402001.
 B. Bertini, E. Tartaglia, and P. Calabrese, Quantum Quench in the Infinitely Repulsive Hubbard Model: The Stationary State, arXiv:1707.01073.
 J. Mossel and J.S. Caux, Generalized TBA and generalized Gibbs, J. Phys. A 45, 255001 (2012).