We provide numerical evidence that after a local quench in an isolated infinite quantum spin chain, the quantum state locally relaxes to the ground state of the post-quenched Hamiltonian, i.e. dissipates. This is a consequence of the unitary quantum dynamics. A mechanism similar to the eigenstate thermalization hypothesis is shown to be responsible for the dissipation observed. We also demonstrate that integrability obstructs dissipation. The numerical simulations are done directly in the thermodynamic limit with a time-evolution algorithm based on matrix product states. The area law of entanglement entropy is observed to hold after the local quench. As a result, the simulations can be performed for long times with small bond dimensions. Various local quenches on the Ising chain and the three-state Potts chain are studied.
We provide numerical evidence that after a local quench in an isolated infinite quantum spin chain, the quantum state locally relaxes to the ground state of the post-quenched Hamiltonian, i.e. dissipates. This is a consequence of the unitary quantum dynamics. A mechanism similar to the eigenstate thermalization hypothesis is shown to be responsible for the dissipation observed. We also demonstrate that integrability obstructs dissipation. The numerical simulations are done directly in the thermodynamic limit with a time-evolution algorithm based on matrix product states. The area law of entanglement entropy is observed to hold after the local quench. As a result, the simulations can be performed for long times with small bond dimensions. Various local quenches on the Ising chain and the three-state Potts chain are studied.
Abstract We review the dynamics after quantum quenches in integrable quantum spin chains. We give a pedagogical introduction to relaxation in isolated quantum systems, and discuss the description of the steady state by (generalized) Gibbs ensembles. We then turn to general features in the time evolution of local observables after the quench, using a simple model of free fermions as an example. In the second part we present an overview of recent progress in describing quench dynamics in two key paradigms for quantum integrable models, the transverse field Ising chain and the anisotropic spin-1/2 Heisenberg chain.
We investigate the eigenstate thermalization hypothesis (ETH) in integrable models, focusing on the $\text{spin}\ensuremath{-}\frac{1}{2}$ isotropic Heisenberg (XXX) chain. We provide numerical evidence that the ETH holds for typical eigenstates (weak ETH scenario). Specifically, using a numerical implementation of state-of-the-art Bethe ansatz results, we study the finite-size scaling of the eigenstate-to-eigenstate fluctuations of the reduced density matrix. We find that fluctuations are normally distributed, and their standard deviation decays in the thermodynamic limit as ${L}^{\ensuremath{-}1/2}$, with L the size of the chain. This is in contrast with the exponential decay that is found in generic nonintegrable systems. Based on our results, it is natural to expect that this scenario holds in other integrable spin models and for typical local observables. Finally, we investigate the entanglement properties of the excited states of the XXX chain. We numerically verify that typical midspectrum eigenstates exhibit extensive entanglement entropy (i.e., volume-law scaling).
Quantum quench is a non-equilibrium process where the Hamiltonian is suddenly changed during the quantum evolution. The change can be made by spatially local perturbations (local quench) or globally switching to a completely different Hamiltonian (global quench). This thesis investigates the post-quench non-equilibrium dynamics with an emphasis on the time dependence of the quantum entanglement. We inspect the scaling of entanglement entropy (EE) to learn how correlation and entanglement built up in a quench. We begin with two local quench examples. In Chap. 2, we apply a local operator to the groundstate of the quantum Lifshitz model and monitor the change of the EE. We find that the entanglement grows according to the dynamical exponent z = 2 and then saturates to the scaling dimension of the perturbing operator -- a value representing its strength. In Chap. 3, we study the evolution after connecting two different one-dimensional critical chains at their ends. The Loschmidt echo which measures the similarity between the evolved state and the initial one decays with a power law, whose exponent is the scaling dimension of the defect (junction). Among other conclusions, we see that the local quench dynamics contain universal information of the (critical) theory. In the global quench scenario, the change of the Hamiltonian affects all parts of the system. In this thesis, we focus on the global chaotic quench driven by generic non-integrable Hamiltonians. In Chap. 4, we propose to use the operator entanglement entropy of the unitary operator as a probe. Its fast linear entanglement production is sharply contrasted to the slow logarithmic spreading of the many-body localized system. The entanglement saturation suggests that the evolution operator in the long time can be modeled by a random unitary matrix. In Chap. 5, we construct a random tensor network which consists of random unitary matrices connected locally to model chaotic evolution with local interactions. We find that the entanglement dynamics is mapped to the statistical mechanics of interacting random walks. This appealing emergent picture allows us to understand the universal linear growth as well as the fluctuations of entanglement in a chaotic quench.
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.
We study dissipative phase transitions in a system of two coupled fully-connected quantum Ising models interacting with an environment. The dynamics is governed by a Lindblad master equation combining coherent unitary evolution and incoherent dissipative processes, where the unitary part is described within a self-consistent mean-field framework effectively acting on the local Hilbert space of two coupled spins at each site. We analyze two fundamentally different classes of dissipators. In the first case, the jump operators are defined in the instantaneous eigenbasis of the mean-field Hamiltonian and satisfy a detailed-balance condition. In this setting, the relaxation dynamics depends strongly on the quench protocol: a parametric quench of the Hamiltonian leads to conventional relaxation, whereas a temperature quench gives rise to a dynamical phase transition characterized by nonanalytic behavior in time. Yet, in both cases, the system relaxes toward a steady state determined solely by the post-quench parameters and the bath temperature, which closely resembles a thermal Gibbs state of the mean-field Hamiltonian. As a result, the dissipative phase transition occurs at a critical point consistent with the corresponding equilibrium transition. In contrast, when the dissipators are realized via local spin raising and lowering operators, the steady state is genuinely nonequilibrium, leading to a significantly richer phase diagram. In particular, for sufficiently strong system-bath coupling, we observe a reentrant phase featuring a symmetry-broken region bounded by two continuous dissipative phase transitions. Our results evidence how the structure of dissipative processes controls the emergence of equilibrium-like versus genuinely nonequilibrium critical behavior in open quantum systems.
Entanglement and entropy are key concepts standing at the foundations of quantum and statistical mechanics. Recently, the study of quantum quenches revealed that these concepts are intricately intertwined. Although the unitary time evolution ensuing from a pure state maintains the system at zero entropy, local properties at long times are captured by a statistical ensemble with nonzero thermodynamic entropy, which is the entanglement accumulated during the dynamics. Therefore, understanding the entanglement evolution unveils how thermodynamics emerges in isolated systems. Alas, an exact computation of the entanglement dynamics was available so far only for noninteracting systems, whereas it was deemed unfeasible for interacting ones. Here, we show that the standard quasiparticle picture of the entanglement evolution, complemented with integrability-based knowledge of the steady state and its excitations, leads to a complete understanding of the entanglement dynamics in the space-time scaling limit. We thoroughly check our result for the paradigmatic Heisenberg chain.