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 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 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 study the dynamics and thermalization of strongly correlated fermions in finite one-dimensional lattices after a quantum quench. Our calculations are performed using exact diagonalization. We focus on one- and two-body observables such as the momentum distribution function [n(k)] and the density-density structure factor [N(k)], respectively, and study the effects of approaching an integrable point. We show that while the relaxation dynamics and thermalization of N(k) for fermions is very similar to the one of hardcore bosons, the behavior of n(k) is distinctively different. The latter observable exhibits a slower relaxation dynamics in fermionic systems. We identify the origin of this behavior, which is related to the off-diagonal matrix elements of n(k) in the basis of the eigenstates of the Hamiltonian. More generally, we find that thermalization occurs far away from integrability and that it breaks down as one approaches the integrable point.
We study the unitary dynamics and the thermalization properties of free-fermion-like Hamiltonians after a sudden quantum quench, extending the results of S. Ziraldo et al. [Phys. Rev. Lett. 109, 247205 (2012)]. With analytical and numerical arguments, we show that the existence of a stationary state and its description with a generalized Gibbs ensemble (GGE) depend crucially on the observable considered (local versus extensive) and on the localization properties of the final Hamiltonian. We present results on two one-dimensional (1D) models, the disordered 1D fermionic chain with long-range hopping and the disordered Ising/XY spin chain. We analytically prove that, while time averages of one-body operators are perfectly reproduced by GGE (even for finite-size systems, if time integrals are extended beyond revivals), time averages of many-body operators might show clear deviations from the GGE prediction when disorder-induced localization of the eigenstates is at play.
We consider dynamic (non-equal-time) correlation functions of local observables after a quantum quench. We show that, in the absence of long-range interactions in the final Hamiltonian, the dynamics is determined by the same ensemble that describes static (equal-time) correlations. For many integrable models, static correlation functions of local observables after a quantum quench relax to stationary values, which are described by a generalized Gibbs ensemble. The same generalized Gibbs ensemble then determines dynamic correlation functions, and the basic form of the fluctuation dissipation theorem holds, although the absorption and emission spectra are not simply related as in the thermal case. For quenches in the transverse field Ising chain, we derive explicit expressions for the time evolution of dynamic order parameter correlators after a quench.
Motivated by recent experiments, we study the relaxation dynamics and thermalization in the one-dimensional Bose-Hubbard model induced by a global interaction quench. Specifically, we start from an initial state that has exactly one boson per site and is the ground state of a system with infinitely strong repulsive interactions at unit filling. Using exact diagonalization and the density-matrix renormalization-group method, we compute the time dependence of such observables as the multiple occupancy and the momentum distribution function. Typically, the relaxation to stationary values occurs over just a few tunneling times. The stationary values are identical to the so-called diagonal ensemble on the system sizes accessible to our numerical methods and we further observe that the microcanonical ensemble describes the time averages of many observables reasonably well for small and intermediate interaction strength. The expectation values of observables in the canonical ensemble agree quantitatively with the time averages obtained from the quench at small interaction strengths, and qualitatively provide a good description even in parameter regimes where the microcanonical ensemble is not applicable due to finite-size effects. We discuss our numerical results in the framework of the eigenstate thermalization hypothesis. Moreover, we also observe that the diagonal and the canonical ensembles are practically identical for our initial conditions already on the level of their respective energy distributions for small interaction strengths. Finally, we discuss implications of our results for the interpretation of a recent sudden expansion experiment [Phys. Rev. Lett. 110, 205301 (2013)], in which the same interaction quench was realized.
We consider the question of thermalization for isolated quantum systems after a sudden parameter change, a so-called quantum quench. In particular we investigate the pre-requisites for thermalization focusing on the statistical properties of the time-averaged density matrix and of the expectation values of observables in the final eigenstates. We find that eigenstates, which are rare compared to the typical ones sampled by the micro-canonical distribution, are responsible for the absence of thermalization of some infinite integrable models and play an important role for some finite size non-integrable systems, such as the Bose-Hubbard model. We stress the importance of finite size effects for the thermalization of isolated quantum systems and propose a condition, relying on the Kullback-Leibler entropy, to obtain thermalization.