We consider the common spin-1/2 XX model in one dimension with open boundary conditions and a large but finite number of spins. The system is in thermal equilibrium at times t < 0 , and is subject to a weak local perturbation (quantum quench) at t = 0 . Focusing mainly on single-spin perturbations and observables, we show that the system rethermalizes for sufficiently large times t > 0 without invoking any unproven assumptions besides the basic laws of quantum mechanics. Moreover, the time-dependent relaxation behavior is obtained in quantitative detail and is found to exhibit a wealth of interesting features.
We consider the common spin-1/2 XX model in one dimension with open boundary conditions and a large but finite number of spins. The system is in thermal equilibrium at times t < 0 , and is subject to a weak local perturbation (quantum quench) at t = 0 . Focusing mainly on single-spin perturbations and observables, we show that the system rethermalizes for sufficiently large times t > 0 without invoking any unproven assumptions besides the basic laws of quantum mechanics. Moreover, the time-dependent relaxation behavior is obtained in quantitative detail and is found to exhibit a wealth of interesting features.
Absence of thermalization after a global quantum quench is a well-established numerical observation in integrable many-body systems, and can be empirically related to a violation of the eigenstate thermalization hypothesis (ETH) in such models. Still, in many of those examples a weaker version of the conventional ETH (wETH) has been numerically reported or even rigorously proven. In this paper we show analytically and illustrate numerically that the absence of thermalization is already possible after a local quench. A closely related finding is a strong violation of the ETH, meaning that not even the wETH is fulfilled anymore. In our analytical explorations we focus on XX-spin-chain models with open boundary conditions, where the local quench is generated by initiating the system in thermal equilibrium and then suddenly switching on (or slightly changing) a single-spin impurity either at the end or in the center of the chain. Numerically we observe qualitatively similar phenomena also for more general XXZ-models in the case of an end-impurity, but not in the case of a central impurity.
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.
Exactly solvable models provide an opportunity to study different aspects of reduced quantum dynamics in detail. We consider the reduced dynamics of a single spin in finite XX and XY spin 1/2 chains. First we introduce a general expression describing the evolution of the reduced density matrix. This expression proves to be tractable when the combined closed system (i.e. open system plus environment) is integrable. Then we focus on comparing decoherence and thermalization timescales in the XX chain. We find that for a single spin these timescales are comparable, in contrast to what should be expected for a macroscopic body. This indicates that the process of quantum relaxation of a system with few accessible states can not be separated in two distinct stages - decoherence and thermalization. Finally, we turn to finite-size effects in the time evolution of a single spin in the XY chain. We observe three consecutive stages of the evolution: regular evolution, partial revivals, irregular (apparently chaotic) evolution. The duration of the regular stage is proportional to the number of spins in the chain. We observe a "quiet and cold period" in the end of the regular stage, which breaks up abruptly at some threshold time.
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).
We study thermalization in a one-dimensional quantum system consisting of a noninteracting fermionic chain with each site of the chain coupled to an additional bath site. Using a density matrix renormalization group algorithm we investigate the time evolution of the observables in the chain after a quantum quench. For low densities we show that the intermediate time dynamics can be quantitatively described by a system of coupled equations of motion. For higher densities our numerical results show a prethermalization for the local observables at intermediate times and a full thermalization to the grand canonical ensemble at long times. For the case of a weak bath-chain coupling we find, in particular, a Fermi momentum distribution in the chain in equilibrium in spite of the seemingly oversimplified bath in our model.
We study the spin and energy dynamics in one-dimensional spin-21 systems induced by local quantum quenches at finite temperatures using a time-dependent density matrix renormalization group method. System sizes are chosen large enough to ensure that the time-dependent data for the accessible time scales represent the behavior in the thermodynamic limit. As a main result, we observe a ballistic spreading of perturbations of the energy density in the integrable spin-21 XXZ chain for all temperatures and exchange anisotropies, related to the divergent thermal conductivity in this model and the exact conservation of the energy current. In contrast, the spin dynamics is ballistic in the massless phase, but shows a diffusive behavior at high temperatures in the easy-axis phase in the case of a vanishing background spin density. We extract a quantitative estimate for the spin-diffusion constant from the time dependence of the spatial variance of the spin density, which agrees well with values obtained from current-current correlation functions using an Einstein relation. Interestingly, the diffusion constant approaches a constant value deep in the easy-axis regime. As an example for nonintegrable models, we consider two-leg ladders, for which we observe indications of diffusive energy and spin dynamics. The relevance of our results for recent experiments with quantum magnets and bosons in optical lattices is discussed.