The eigenstate thermalization hypothesis provides one picture of thermalization in a quantum system by looking at individual eigenstates. However, it is also important to consider how local observables reach equilibrium values dynamically. Quench protocol is one of the settings to study such questions. A recent numerical study [Ba\~nuls, Cirac, and Hastings, Phys. Rev. Lett. 106, 050405 (2007)] of a nonintegrable quantum Ising model with longitudinal field under such a quench setting found different behaviors for different initial quantum states. One particular case called the ``weak-thermalization'' regime showed apparently persistent oscillations of some observables. Here we provide an explanation of such oscillations. We note that the corresponding initial state has low energy density relative to the ground state of the model. We then use perturbation theory near the ground state and identify the oscillation frequency as essentially a quasiparticle gap. With this quasiparticle picture, we can then address the long-time behavior of the oscillations. Upon making additional approximations which intuitively should only make thermalization weaker, we argue that the oscillations nevertheless decay in the long-time limit. As part of our arguments, we also consider a quench from a BEC to a hard-core boson model in one dimension. We find that the expectation value of a single-boson creation operator oscillates but decays exponentially in time, while a pair-boson creation operator has oscillations with a ${t}^{\ensuremath{-}3/2}$ decay in time. We also study dependence of the decay time on the density of bosons in the low-density regime and use this to estimate decay time for oscillations in the original spin model.
The eigenstate thermalization hypothesis provides one picture of thermalization in a quantum system by looking at individual eigenstates. However, it is also important to consider how local observables reach equilibrium values dynamically. Quench protocol is one of the settings to study such questions. A recent numerical study [Ba~nuls, Cirac, and Hastings, Phys. Rev. Lett. 106, 050405 (2007)] of a nonintegrable quantum Ising model with longitudinal field under such a quench setting found different behaviors for different initial quantum states. One particular case called the ``weak-thermalization'' regime showed apparently persistent oscillations of some observables. Here we provide an explanation of such oscillations. We note that the corresponding initial state has low energy density relative to the ground state of the model. We then use perturbation theory near the ground state and identify the oscillation frequency as essentially a quasiparticle gap. With this quasiparticle picture, we can then address the long-time behavior of the oscillations. Upon making additional approximations which intuitively should only make thermalization weaker, we argue that the oscillations nevertheless decay in the long-time limit. As part of our arguments, we also consider a quench from a BEC to a hard-core boson model in one dimension. We find that the expectation value of a single-boson creation operator oscillates but decays exponentially in time, while a pair-boson creation operator has oscillations with a t\ensuremath−3/2 decay in time. We also study dependence of the decay time on the density of bosons in the low-density regime and use this to estimate decay time for oscillations in the original spin model.
We study the dynamics of the quantum Ising chain following a zero-temperature quench of the transverse field strength. Focusing on the behavior of two-point spin correlation functions, we show that the correlators of the order parameter display an effective asymptotic thermal behavior; i.e., they decay exponentially to zero with a phase coherence rate and a correlation length dictated by the equilibrium law with an effective temperature set by the energy of the initial state. On the contrary, the two-point correlation functions of the transverse magnetization or the density-of-kinks operator decay as a power law and do not exhibit thermal behavior. We argue that the different behavior is linked to the locality of the corresponding operator with respect to the quasiparticles of the model: nonlocal operators, such as the order parameter, behave thermally, while local ones do not. We study which features of the two-point correlators are a consequence of the integrability of the model by analyzing their robustness with respect to a sufficiently strong integrability-breaking term.
The results presented in this PhD Thesis aim at understanding thermalization dynamics and, strictly speaking, pre-thermalization in both open and closed quantum spin chains driven out-of-equilibrium. In order to make analytical progress, we consider perturbations of the Quantum Ising Chain (QIC), one of the simplest integrable models, whose non-equilibrium dynamics has been studied in great details in the last few years. In Chapter 2 we introduce for the non-technical reader the main recent ideas behind quantum thermalization and relaxation in integrable and non-integrable quantum many body systems, while Chapter 3 contains an introduction to the current understanding of pre-thermalization. In Chapter 4 we start our study considering a QIC perturbed by a time dependent noise in the direction of the transverse magnetic field [6] [7]. We generalize the concept of pre-thermalization to a noisy quantum many body system driven away from equilibrium simultaneously by a quantum quench of the transverse field and by the noise, highlighting the mechanism behind pre-thermalization and the subsequent relaxation towards the thermal state in a wide class of physically relevant observables. Our main achievement is to show that pre-thermalization occurs quite generally up to time scales where the noise is dominated by the non-equilibrium dynamics induced by the quench of the transverse field, while for longer times pre-thermal plateaux disappear and the noise leaves different signatures in different observables. In Chapter 5 we consider the more challenging problem of pre-thermalization in a closed non-integrable quantum spin chain [8]. We show using a mapping of a weakly non-integrable spin chain into a hard core boson model, that pre-thermalization is successfully captured by a low-density description of the theory, which results in a integrable, bosonic model, that constitutes an excellent description for the out-of-equilibrium dynamics up to intermediate time scales. On the other hand, in order to investigate the late time dynamics we resort to a diagrammatic approach, which serves for a two-fold scope. First, we focus on the effect of the leading inelastic processes, realizing that they constitute only a sub-leading correction to the pre-thermal plateau, found within the effective low-energy description; and secondly, we show from a Quantum Boltzmann equation perspective how inelastic processes are the fundamental ingredient in order to have thermalization in the late time dynamics. The main technical tool used to derive the results in Chapter 4 and 5 is the Keldysh contour technique which is reviewed in the Appendix. This PhD Thesis is based on the following works: 1) J. Marino, A. Silva, Relaxation, pre-thermalization and diffusion in a noisy Quantum Ising Chain, Phys. Rev. B 86, 060408 (Rapid Communications) (2012). 2) M. Marcuzzi, J. Marino, A. Gambassi, A. Silva, Pre-thermalization in a non-integrable quantum spin chain after a quench, arXiv:1307.3738 (2013) (submitted to Phys. Rev. Lett.). 3) J. Marino, A. Silva, Non-Equilibrium Dynamics of a Noisy Quantum Ising Chain: statistics of the work and prethermalization after a sudden quench of the transverse field, arXiv:1309.7595 (2013) (submitted to Phys. Rev. B).
Quantum chaos and the eigenstate thermalization hypothesis are based on the assumption of the validity of random matrix theory description on the spectrum and eigenstates. They provide the foundation and descriptions for the typical dynamics and thermalization in generic closed quantum systems. In this thesis, we investigate situations where the systems show atypical dynamics or anomalous thermalization, conflicting with the usual expectations from quantum chaos and eigenstate thermalization hypothesis. We first examine weak thermalization in a nonintegrable spin chain. The system shows long-lived strong oscillations and relaxes to the thermal equilibrium weakly. We identify the dynamics describable by quasiparticles and recognize the oscillation frequency to be the quasiparticle mass gap. We also estimate the damping time for the oscillations. Next, we study prethermalization, a phenomenon where a system relaxes to an intermediate almost-equilibrium stage before reaching the true thermal equilibrium. We study a nonintegrable spin chain in the strong coupling limit, where an almost-conserved quantity emerges and gives rise to the prethermalization. We also study a newly proposed diagnostic for quantum chaos: out-of-time-ordered correlators. Contrasting to the chaotic systems, we inspect their behaviors in various noninteracting integrable models. Finally, we dig into the quantum many-body scar states in the PXP model which describes a Rydberg atom chain. These special states do not satisfy the random matrix theory description nor the eigenstate thermalization hypothesis, therefore defying quantum chaos.
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.
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.
After a quench in a quantum many-body system, expectation values tend to relax towards long-time averages. However, temporal fluctuations remain in the long-time limit, and it is crucial to study the suppression of these fluctuations with increasing system size. The particularly important case of nonintegrable models has been addressed so far only by numerics and conjectures based on analytical bounds. In this work, we are able to derive analytical predictions for the temporal fluctuations in a nonintegrable model (the transverse Ising chain with extra terms). Our results are based on identifying a dynamical regime of "many-particle dephasing," where quasiparticles do not yet relax but fluctuations are nonetheless suppressed exponentially by weak integrability breaking.