We consider quenches of a quantum system that is prepared in a canonical equilibrium state of one Hamiltonian and then evolves unitarily in time under a different Hamiltonian. Technically, our main result is a systematic expansion of the pre- and postquench canonical ensembles in the quench strength. We first demonstrate how this can be used to predict the system's temperature after the quench from equilibrium properties at the prequench temperature. For a thermalizing postquench system, it furthermore allows us to calculate equilibrium observable expectation values. Finally, in the presence of additional conserved quantities besides the Hamiltonian, we obtain a hierarchy of necessary conditions for thermalization towards the (postquench) canonical ensemble. At first order, these thermalization conditions have a nice geometric interpretation in operator space with the canonical covariance as a semi-inner product: The quench operator (difference between post- and prequench Hamiltonians) and the conserved quantity must be orthogonal in the orthogonal complement of the postquench Hamiltonian. We illustrate the results numerically for a variety of setups involving integrable and nonintegrable models.
We consider quenches of a quantum system that is prepared in a canonical equilibrium state of one Hamiltonian and then evolves unitarily in time under a different Hamiltonian. Technically, our main result is a systematic expansion of the pre- and postquench canonical ensembles in the quench strength. We first demonstrate how this can be used to predict the system's temperature after the quench from equilibrium properties at the prequench temperature. For a thermalizing postquench system, it furthermore allows us to calculate equilibrium observable expectation values. Finally, in the presence of additional conserved quantities besides the Hamiltonian, we obtain a hierarchy of necessary conditions for thermalization towards the (postquench) canonical ensemble. At first order, these thermalization conditions have a nice geometric interpretation in operator space with the canonical covariance as a semi-inner product: The quench operator (difference between post- and prequench Hamiltonians) and the conserved quantity must be orthogonal in the orthogonal complement of the postquench Hamiltonian. We illustrate the results numerically for a variety of setups involving integrable and nonintegrable models.
Abstract The observable long-time behavior of an isolated many-body system after a quantum quench is considered, i.e. an eigenstate (or an equilibrium ensemble) of some pre-quench Hamiltonian H serves as initial condition which then evolves in time according to some post-quench Hamiltonian H p . Absence of thermalization is analytically demonstrated for a large class of quite common pre- and post-quench spin Hamiltonians. The main requirement is that the pre-quench Hamiltonian must exhibit a Z 2 (spin-flip) symmetry, which would be spontaneously broken in the thermodynamic limit, though we actually focus on finite (but large) systems. On the other hand, the post-quench Hamiltonian must violate the Z 2 symmetry, but for the rest may be non-integrable and may obey the eigenstate thermalization hypothesis for (sums of) few-body observables.
We study the expectation values of observables and correlation functions at long times after a global quantum quench. Our focus is on metallic (`gapless') fermionic many-body models and small quenches. The system is prepared in an eigenstate of an initial Hamiltonian, and the time evolution is performed with a final Hamiltonian which differs from the initial one in the value of one global parameter. We first derive general relations between time-averaged expectation values of observables as well as correlation functions and those obtained in an eigenstate of the final Hamiltonian. Our results are valid to linear and quadratic order in the quench parameter g and generalize prior insights in several essential ways. This allows us to develop a phenomenology for the thermalization of local quantities up to a given order in g. Our phenomenology is put to a test in several case studies of one-dimensional models representative of four distinct classes of Hamiltonians: quadratic ones, effectively quadratic ones, those characterized by an extensive set of (quasi-) local integrals of motion, and those for which no such set is known (and believed to be nonexistent). We show that for each of these models, all observables and correlation functions thermalize to linear order in g. The more local a given quantity, the longer the linear behavior prevails when increasing g. Typical local correlation functions and observables for which the term O(g) vanishes thermalize even to order g^2. Our results show that lowest order thermalization of local observables is an ubiquitous phenomenon even in models with extensive sets of integrals of motion.
We study properties of isolated integrable quantum systems after a sudden quench starting from thermal states. We show that, even if the system is initially in thermal equilibrium at finite temperature, the diagonal entropy after a quench remains a fraction of the entropy in the generalized ensembles introduced to describe integrable systems after relaxation. The latter is also, in general, different from the entropy in thermal equilibrium. Furthermore, we examine the difference in the distributions of conserved quantities in the thermal versus the generalized ensembles after a quench and show that they are also, in general, different from each other. This explains why these systems fail to thermalize. A finite-size scaling analysis is presented for each quantity, which allows us to make predictions for thermodynamically large lattice sizes.
We study the expectation values of observables and correlation functions at long times after a global quantum quench. Our focus is on metallic (``gapless'') fermionic many-body models and small quenches. The system is prepared in an eigenstate of an initial Hamiltonian, and the time evolution is performed with a final Hamiltonian which differs from the initial one in the value of one global parameter. We first derive general relations between time-averaged expectation values of observables as well as correlation functions and those obtained in an eigenstate of the final Hamiltonian. Our results are valid to linear and quadratic order in the quench parameter g and generalize prior insights in several essential ways. This allows us to develop a phenomenology for the thermalization of local quantities up to a given order in g. Our phenomenology is put to a test in several case studies of one-dimensional models representative of four distinct classes of Hamiltonians: quadratic ones, effectively quadratic ones, those characterized by an extensive set of (quasi-) local integrals of motion, and those for which no such set is known (and believed to be nonexistent). We show that for each of these models, all observables and correlation functions thermalize to linear order in g. The more local a given quantity, the longer the linear behavior prevails when increasing g. Typical local correlation functions and observables for which the term O(g) vanishes thermalize even to order g2. Our results show that lowest-order thermalization of local observables is an ubiquitous phenomenon even in models with extensive sets of integrals of motion.
We identify and study classes of initial states in integrable quantum systems that, after the relaxation dynamics following a sudden quench, lead to near-thermal expectation values of few-body observables. In the systems considered here, those states are found to be insulating ground states of lattice hard-core boson Hamiltonians. We show that, as a suitable parameter in the initial Hamiltonian is changed, those states become closer to Fock states (products of single site states) as the outcome of the relaxation dynamics becomes closer to the thermal prediction. At the same time, the energy density approaches a Gaussian. Furthermore, the entropy associated with the generalized canonical and generalized grand-canonical ensembles, introduced to describe observables in integrable systems after relaxation, approaches that of the conventional canonical and grand-canonical ensembles. We argue that those classes of initial states are special because a control parameter allows one to tune the distribution of conserved quantities to approach the one in thermal equilibrium. This helps in understanding the approach of all the quantities studied to their thermal expectation values. However, a finite-size scaling analysis shows that this behavior should not be confused with thermalization as understood for nonintegrable systems.
Nonequilibrium time evolution in isolated many-body quantum systems generally results in thermalization. However, the relaxation process can be very slow, and quasistationary nonthermal plateaux are often observed at intermediate times. The paradigmatic example is a quantum quench in an integrable model with weak integrability breaking; for a long time, the state cannot escape the constraints imposed by the approximate integrability. We unveil a new mechanism of prethermalization, based on the presence of a symmetry of the prequench Hamiltonian, which is spontaneously broken at zero temperature and is explicitly broken by the postquench Hamiltonian. The typical time scale of the phenomenon is proportional to the thermal correlation length of the initial state, which diverges as the temperature is lowered. We show that the prethermal quasistationary state can be approximated by a mixed state that violates cluster decomposition property. We consider two examples: the transverse-field Ising chain, where the full-time evolution is computed analytically, and the (nonintegrable) anisotropic next-nearest-neighbor Ising model, which is investigated numerically.