Global Quantum Quench Research Articles

Quantum quench and scaling of entanglement entropy

Global quantum quench with a finite quench rate which crosses critical points is known to lead to universal scaling of correlation functions as functions of the quench rate. In this work, we explore scaling properties of the entanglement entropy of a subsystem in a harmonic chain during a mass quench which asymptotes to finite constant values at early and late times and for which the dynamics is exactly solvable. When the initial state is the ground state, we find that for large enough subsystem sizes the entanglement entropy becomes independent of size. This is consistent with Kibble–Zurek scaling for slow quenches, and with recently discussed “fast quench scaling” for quenches fast compared to physical scales, but slow compared to UV cutoff scales.

Small quenches and thermalization

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.

Extended slow dynamical regime close to the many-body localization transition

Many-body localization is characterized by a slow logarithmic growth of the entanglement entropy after a global quantum quench while the local memory of an initial density imbalance remains at infinite time. We investigate how much the proximity of a many-body localized phase can influence the dynamics in the delocalized ergodic regime where thermalization is expected. Using an exact Krylov space technique, the out-of-equilibrium dynamics of the random-field Heisenberg chain is studied up to $L=28$ sites, starting from an initially unentangled high-energy product state. Within most of the delocalized phase, we find a sub-ballistic entanglement growth $S(t)\propto t^{1/z}$ with a disorder-dependent exponent $z\ge1$, in contrast with the pure ballistic growth $z=1$ of clean systems. At the same time, anomalous relaxation is also observed for the spin imbalance ${\cal{I}}(t)\propto t^{-\zeta}$ with a continuously varying disorder-dependent exponent $\zeta$, vanishing at the transition. This provides a clear experimental signature for detecting this non-conventional regime.

Quench dynamics in confined 1 + 1-dimensional systems

We present a framework for investigating the response of conformally invariant confined 1 + 1-dimensional systems to a quantum quench. While conformal invariance is generally destroyed in a global quantum quench, systems that can be described as or mapped to integrable deformations of a CFT may present special instances where a conformal field theory-based analysis could provide useful insight into the non-equilibrium dynamics. We investigate this possibility by considering a quench analogous to that of the quantum Newton’s Cradle experiment (Kinoshita et al 2006 Nature 440 900) and demonstrating qualitative agreement between observables derived in the CFT framework and those of the experimental system. We propose that this agreement may be a feature of the proximity of the experimental system to an integrable deformation of a c = 1 CFT.

Signatures of the many-body localization transition in the dynamics of entanglement and bipartite fluctuations

The many-body localization transition is a dynamical quantum phase transition between a localized and an extended phase. We study this transition in the XXZ model with disordered magnetic field and focus on the time evolution following a global quantum quench. While the dynamics of the bipartite entanglement and spin fluctuations are already known to provide insights into the nature of the many-body localized phases, we discuss the relevance of these quantities in the context of the localization transition. In particular, we observe that near the transition the long time limits of both quantities show behavior similar to divergent thermodynamic fluctuations.

Entanglement tsunami in (1+1)-dimensions

We study the time dependence of the entanglement entropy of disjoint intervals following a global quantum quench in (1+1)-dimensional CFTs at large-$c$ with a sparse spectrum. The result agrees with a holographic calculation but differs from the free field theory answer. In particular, a simple model of free quasiparticle propagation is not adequate for CFTs with a holographic dual. We elaborate on the entanglement tsunami proposal of Liu and Suh and show how it can be used to reproduce the holographic answer.

On the time evolution of holographic n-partite information

We study various scaling behaviors of n-partite information during a process of thermalization after a global quantum quench for n disjoint system consisting of n parallel strips whose widths are much larger than the separation between them. By making use of the holographic description for entanglement entropy we explore holographic description of the n-partite information by which we show that it has a definite sign: it is positive for even n and negative for odd n. This might be thought of as an intrinsic property of a field theory which has gravity dual.

Entanglement negativity after a global quantum quench

We study the time evolution of the logarithmic negativity after a global quantum quench. In a 1 + 1-dimensional conformal invariant field theory, we consider the negativity between two intervals which can be either adjacent or disjoint. We show that the negativity follows the quasi-particle interpretation for the spreading of entanglement. We check and generalise our findings with a systematic analysis of the negativity after a quantum quench in the harmonic chain, highlighting two peculiar lattice effects: the late birth and the sudden death of entanglement.

Quench echo and work statistics in integrable quantum field theories.

We propose a boundary thermodynamic Bethe ansatz calculation technique to obtain the Loschmidt echo and the statistics of the work done when a global quantum quench is performed on an integrable quantum field theory. We derive an analytic expression for the lowest edge of the probability density function and find that it exhibits universal features, in the sense that its scaling form depends only on the statistics of excitations. We perform numerical calculations on the sinh-Gordon model, a deformation of the free boson theory, and we obtain that by turning on the interaction the density function develops fermionic properties. The calculations are facilitated by a previously unnoticed property of the thermodynamic Bethe ansatz construction.

Entanglement contour

In the context of characterizing the structure of quantum entanglement in many-body systems, we introduce the entanglement contour, a tool to identify which real-space degrees of freedom contribute, and how much, to the entanglement of a region A with the rest of the system B. The entanglement contour provides a complementary, more refined approach to characterizing entanglement than just considering the entanglement entropy between A and B, with several concrete advantages. We illustrate this in the context of ground states and quantum quenches in fermionic quadratic systems. For instance, in a quantum critical system in D = 1 spatial dimensions, the entanglement contour allows us to determine the central charge of the underlying conformal field theory from just a single partition of the system into regions A and B (using the entanglement entropy for the same task requires considering several partitions). In D ⩾ 2 dimensions, the entanglement contour can distinguish between gapped and gapless phases that obey the same boundary law for entanglement entropy. During a local or global quantum quench, the time-dependent contour provides a detailed account of the dynamics of entanglement, including propagating entanglement waves, which offers a microscopic explanation of the behavior of the entanglement entropy as a function of time.

Reduced density matrix after a quantum quench

We consider the reduced density matrix (RDM) \rho_A(t) for a finite subsystem A after a global quantum quench in the infinite transverse-field Ising chain. It has been recently shown that the infinite time limit of \rho_A(t) is described by the RDM \rho_{GGE,A} of a generalized Gibbs ensemble. Here we present some details on how to construct this ensemble in terms of local integrals of motion, and show its equivalence to the expression in terms of mode occupation numbers widely used in the literature. We then address the question, how \rho_A(t) approaches \rho_{GGE,A} as a function of time. To that end we introduce a distance on the space of density matrices and show that it approaches zero as a universal power-law t^{-3/2} in time. As the RDM completely determines all local observables within A, this provides information on the relaxation of correlation functions of local operators. We then address the issue, of how well a truncated generalized Gibbs ensemble with a finite number of local higher conservation laws describes a given subsystem at late times. We find that taking into account only local conservation laws with a range at most comparable to the subsystem size provides a good description. However, excluding even a single one of the most local conservation laws in general completely spoils this agreement.