Dynamic correlation and response functions of classical and quantum systems in thermal equilibrium are connected by fluctuation-dissipation theorems, which allow an alternative definition of their (unique) temperature. Motivated by this fundamental property, we revisit the issue of thermalization of closed many-body quantum systems long after a sudden quench, focusing on the nonequilibrium dynamics of the Ising chain in a critical transverse field. We show the emergence of distinct observable-dependent effective temperatures, which rule out Gibbs thermalization in a strict sense but might still have a thermodynamic meaning.
Dynamic correlation and response functions of classical and quantum systems in thermal equilibrium are connected by fluctuation-dissipation theorems, which allow an alternative definition of their (unique) temperature. Motivated by this fundamental property, we revisit the issue of thermalization of closed many-body quantum systems long after a sudden quench, focusing on the nonequilibrium dynamics of the Ising chain in a critical transverse field. We show the emergence of distinct observable-dependent effective temperatures, which rule out Gibbs thermalization in a strict sense but might still have a thermodynamic meaning.
Abstract Fluctuation-dissipation relations, i.e., the relation between two-time correlation and linear response functions, were successfully used to search for signs of equilibration and to identify effective temperatures in the non-equilibrium behavior of a number of macroscopic classical and quantum systems in contact with thermal baths. Among the most relevant cases in which the effective temperatures thus defined were shown to have a thermodynamic meaning one finds the stationary dynamics of driven super-cooled liquids and vortex glasses, and the relaxation of glasses. Whether and under which conditions an effective thermal behavior can be found in quantum isolated many-body systems after a global quench is a question of current interest. We propose to study the possible emergence of thermal behavior long after the quench by studying fluctuation-dissipation relations in which (possibly time- or frequency-dependent) parameters replace the equilibrium temperature. If thermalization within the Gibbs ensemble eventually occurs these parameters should be constant and equal for all pairs of observables in ‘partial’ or ‘mutual’ equilibrium. We analyze these relations in the paradigmatic quantum system, i.e., the quantum Ising chain, in the stationary regime after a quench of the transverse field. The lack of thermalization to a Gibbs ensemble becomes apparent within this approach.
In this Thesis we discussed several issues related to the dynamics of isolated quantum many-body systems. The brief overview of Chapter 1 highlighted a series of open fundamental theoretical questions which are crucial for the investigation of nonequilibrium phenomena. Among them, we mainly focused on the emergence of thermal behaviors and on the characterization of dynamical phase transitions and of the associated universal dynamics. To this end, we analyzed some paradigmatic models, such as the Fermi-Hubbard model, the quantum Ising chain, and the O(N) model. -Chapter 2: Absence of thermalization in a Fermi liquid. In this Chapter we report two results that question the common expectation of thermalization in a Landau-Fermi liquid after an interaction quench. We first consider the perturbative expansion in the interaction strength of the long-wavelength structure factor S(q) and show that it does not satisfy the hypothesis that steady-state averages correspond to thermal ones. In particular, S(q) has an analytical component ∼ const. + O(q^2) for q → 0, compatible with thermalization, but it retains also a nonanalytic term ∼ |q| characteristic of a Fermi liquid at zero-temperature. In real space, this nonanalyticity turns in a power-law decay of the density-density correlations, in contrast with the exponential decay associated to thermalization. We next analyze the low-density case, where one can obtain results valid at any order in interaction but at leading in the density, and find that in the steady state the momentum distribution jump at the Fermi surface is strictly finite, though smaller than at equilibrium. -Chapter 3: Nonadiabatic stationary behavior in a driven Ising chain. In this Chapter we discuss the emergence of a nonadiabatic behavior in the dynamics of the order parameter in a low-dimensional system subject to a linear ramp of one of its parameters within a gapped phase, which should be the most favorable situation for an adiabatic evolution. We study in details the dynamics of the spontaneous magnetization m_x(t) in a quantum Ising chain after a linear variation in time of the transverse field within the ordered phase. In particular, focusing on the value of the order parameter at the end of the ramp m_x(τ), we show that the smaller the switching rate of the transverse field is the closer m_x(τ) gets to its ground state value. Nonetheless, this apparently small correction to adiabaticity eventually leads to the disruption of the order exponentially fast in the subsequent time evolution, no matter how slowly the ramp is performed. -Chapter 4: Aging and coarsening in isolated quantum systems after a quench. This Chapter analyzes the nonequilibrium dynamics of an isolated quantum system after a sudden quench to the dynamical critical point, where the emergence of scaling and universal exponents is expected, due to the absence of time and length scales. We explore these features for a bosonic interacting scalar field theory with O(N) symmetry in the large-N limit, where the model is exactly solvable and the exponents and scaling forms of the relevant two-time correlation functions can be analytically calculated. Moreover, we provide evidence of the emergence of a dynamical scaling behavior also for quenches below the dynamical critical point, associated with coarsening. We find that the latter case is characterized by the same scaling functions as those describing the critical case, yet with different exponents. -Chapter 5: Dynamical transitions and statistics of excitations. In this Chapter we study the dynamics of the O(N) model (introduced in the previous Chapter) resulting from a different protocol: a linear ramp of one of its parameters. We find that the presence of a dynamical phase transition, as well as its critical properties, are robust against the change of the protocol. We show that a characterization based on the critical dimensions and exponents would suggest that the dynamical phase transition is analogous to the equilibrium thermal one. However, its nonequilibrium nature becomes evident analyzing the statistics of excitations produced in the ramp process. In particular, the critical properties are encoded in the fluctuations in the number of excitations, which display qualitatively different behaviors depending on the ramp being performed above, at, or below the dynamical critical point. These behaviors bear no dependence on the duration of the protocol.
We explore the robustness of universal dynamical scaling behavior in a quantum system near criticality with respect to initialization in a large class of states with finite energy. By focusing on a homogeneous XY quantum spin chain in a transverse field, we characterize the nonequilibrium response under adiabatic and sudden quench processes originating from a pure as well as a mixed excited initial state, and involving either a regular quantum critical or a multicritical point. We find that the critical exponents of the ground-state quantum phase transition can be encoded in the dynamical scaling exponents despite the finite energy of the initial state. In particular, we identify conditions on the initial distribution of quasiparticle excitation that ensure that Kibble-Zurek scaling persists. The emergence of effective thermal equilibrium behavior following a sudden quench toward criticality is also investigated, with focus on the long-time expectation value of the quasiparticle number operator. Despite the integrability of the XY model, this observable is found to behave thermally in quenches to a regular quantum critical point, provided that the system is initially prepared at sufficiently high temperature. However, a similar thermalization behavior fails to occur in quenches toward a multicritical point. We argue that the observed lack of thermalization originates in this case in the asymmetry of the impulse region that is also responsible for anomalous multicritical dynamical scaling.
We study the nonequilibrium dynamics of the quantum Ising model following an abrupt quench of the transverse field. We focus on the on-site autocorrelation function of the order parameter, and extract the phase-coherence time tau(Q)(phi) from its asymptotic behavior. We show that the initial state determines tau(Q)(phi) only through an effective temperature set by its energy and the final Hamiltonian. Moreover, we observe that the dependence of tau(Q)(phi) on the effective temperature fairly agrees with that obtained in thermal equilibrium as a function of the equilibrium temperature.
We study the statistics of the work done, fluctuation relations, and irreversible entropy production in a quantum many-body system subject to the sudden quench of a control parameter. By treating the quench as a thermodynamic transformation we show that the emergence of irreversibility in the nonequilibrium dynamics of closed many-body quantum systems can be accurately characterized. We demonstrate our ideas by considering a transverse quantum Ising model that is taken out of equilibrium by an instantaneous change of the transverse field.
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.