The confinement of elementary excitations induces distinctive features in the non-equilibrium quench dynamics. One of the most remarkable is the suppression of entanglement entropy, which in several instances turns out to oscillate rather than grow indefinitely. While the qualitative physical origin of this behavior is clear, till now no quantitative understanding away from the field theory limit was available. Here we investigate this problem in the weak quench limit, when mesons are excited at rest, hindering entropy growth and exhibiting persistent oscillations. We provide analytical predictions of the entire entanglement dynamics based on a Gaussian approximation of the many-body state, which captures numerical data with great accuracy and is further simplified to a semiclassical quasiparticle picture in the regime of weak confinement. Our methods are valid in general and we apply explicitly to two prototypical models: the Ising chain in a tilted field and the experimentally relevant long-range Ising model.
The confinement of elementary excitations induces distinctive features in the non-equilibrium quench dynamics. One of the most remarkable is the suppression of entanglement entropy, which in several instances turns out to oscillate rather than grow indefinitely. While the qualitative physical origin of this behavior is clear, till now no quantitative understanding away from the field theory limit was available. Here we investigate this problem in the weak quench limit, when mesons are excited at rest, hindering entropy growth and exhibiting persistent oscillations. We provide analytical predictions of the entire entanglement dynamics based on a Gaussian approximation of the many-body state, which captures numerical data with great accuracy and is further simplified to a semiclassical quasiparticle picture in the regime of weak confinement. Our methods are valid in general and we apply explicitly to two prototypical models: the Ising chain in a tilted field and the experimentally relevant long-range Ising model.
The nonequilibrium dynamics of disordered many-body quantum systems after a quantum quench unveils important insights about the competition between interactions and disorder, yielding, in particular, an interesting perspective toward the understanding of many-body localization. Still, the experimentally relevant effect of bond randomness in long-range interacting spin chains on their dynamical properties have so far not been investigated. In this Letter, we examine the entanglement entropy growth after a global quench in a quantum spin chain with randomly placed spins and long-range tunable interactions decaying with distance with power $\ensuremath{\alpha}$. Using a dynamical version of the strong disorder renormalization group we find for $\ensuremath{\alpha}>{\ensuremath{\alpha}}_{c}$ that the entanglement entropy grows logarithmically with time and becomes smaller with larger $\ensuremath{\alpha}$ as $S(t)={S}_{p}ln(t)/(2\ensuremath{\alpha})$. Here, ${S}_{p}=2ln2\ensuremath{-}1$. We present results of numerical exact diagonalization calculations for system sizes up to $N\ensuremath{\sim}16$ spins, in good agreement with the analytical results for sufficiently large $\ensuremath{\alpha}>{\ensuremath{\alpha}}_{c}\ensuremath{\approx}1.8$. For $\ensuremath{\alpha}<{\ensuremath{\alpha}}_{c}$, we find that the entanglement entropy grows as a power law with time, $S(t)\ensuremath{\sim}{t}^{\ensuremath{\gamma}(\ensuremath{\alpha})}$ with $0<\ensuremath{\gamma}(\ensuremath{\alpha})<1$ a decaying function of the interaction exponent $\ensuremath{\alpha}$.
The non-equilibrium dynamics of disordered many-body quantum systems after a global quantum quench unveils important insights about the competition between interactions and disorder, yielding in particular an insightful perspective on many body localization (MBL). Still, the experimentally relevant effect of bond randomness in long-range interacting spin chains on the quantum quench dynamics have so far not been investigated. In this letter, we examine the entanglement entropy growth after a global quench in a quantum spin chain with randomly placed spins and long-range tunable interactions decaying with distance with power α. Using a dynamical version of the strong disorder renormalization group (SDRG) we find for α>αc that the entanglement entropy grows logarithmically with time and becomes smaller with larger α as S(t)=Spln(t)/(2α). Here, Sp=2ln2−1. We use numerical exact diagonalization (ED) simulations to verify our results for system sizes up to $ N\sim 16$ spins, yielding good agreement for sufficiently large α>αc≈1.8. For α<αc, we find that the entanglement entropy grows as a power-law with time, S(t)∼tγ(α) with 0<γ(α)<1 a decaying function of the interaction exponent α.
The task of exploring and understanding various aspects of far-from-equilibrium dynamics of closed and generic quantum many-body systems has received a thrust of attention in recent years, driven partly by remarkable advances in ultracold experimental technologies. In this work, for the paradigmatic Ising spin chain with transverse and longitudinal fields and partly motivated by the practice of site-resolved control in contemporary ultracold experiments, we present numerical observations of several fine-grained (small-subsystem level) features of far-from-equilibrium dynamics from a quantum informational point of view, induced by quantum quenches far across the Ising critical point between states deep inside the para- and ferro-magnetic regimes. Rather featureless dynamics is seen for ferromagnetic to paramagnetic quenches, but paramagnetic to ferromagnetic quenches exhibit rich behaviour, including recurrences of an approximately Page-like dynamics of entanglement entropies of one- and two-spin subsystems, periodic but short-lived occurrences of approximately 1−uniform states, a series of sudden deaths and revivals of entanglement between two spins in the system's bulk, non-analytic cusps in single-copy entanglement entropy for three-spin and bigger subsystems, insufficient mixedness and a series of scrambling-unscrambling of local mutual information between neighboring spins. Moreover, essentially indistinguishable dynamics is seen at very early times between the integrable limit (zero longitudinal field) and non-integrable cases, with the former eventually showing signatures of better mixing and faster approach to equilibration than the latter. These features are expected to hold for quench dynamics across Ising quantum critical points in more complicated systems.
The time evolution of the entanglement entropy in non-equilibrium quantum systems provides crucial information about the structure of the time-dependent state. For quantum quench protocols, by combining a quasiparticle picture for the entanglement spreading with the exact knowledge of the stationary state provided by Bethe ansatz, it is possible to obtain an exact and analytic description of the evolution of the entanglement entropy. Here we discuss the application of these ideas to several integrable models. First we show that for non-interacting systems, both bosonic and fermionic, the exact time-dependence of the entanglement entropy can be derived by elementary techniques and without solving the dynamics. We then provide exact results for interacting spin chains that are carefully tested against numerical simulations. Finally, we apply this method to integrable one-dimensional Bose gases (Lieb-Liniger model) both in the attractive and repulsive regimes. We highlight a peculiar behaviour of the entanglement entropy due to the absence of a maximum velocity of excitations.
We consider non-equilibrium dynamics after quantum quenches in the mixed-field three-state Potts quantum chain in the ferromagnetic regime. Compared to the analogous setting for the Ising spin chain, the Potts model has a much richer phenomenology, which originates partly from baryonic excitations in the spectrum and partly from the various possible relative alignments of the initial magnetisation and the longitudinal field. We obtain the excitation spectrum by combining semiclassical approximation and exact diagonalisation, and we use the results to explain the various dynamical behaviours we observe. Besides recovering dynamical confinement, as well as Wannier-Stark localisation due to Bloch oscillations similar to the Ising chain, a novel feature is the presence of baryonic excitations in the quench spectroscopy. In addition, when the initial magnetisation and the longitudinal field are misaligned, both confinement and Bloch oscillations only result in partial localisation, with some correlations retaining an unsuppressed light-cone behaviour together with a corresponding growth of entanglement entropy.
The transverse field Ising model of N fully connected spins provides a mathematical accessible, yet non-trivial, system to investigate the non-thermal dynamics after a sudden quantum quench. The focus in this thesis is three-fold. First, in a 1/N expansion the dynamics of the mean magnetization within its variance is obtained analytically in the large N limit. The variance constitutes a leading order correction to the mean field limit at infinite N. While mean field theories are one of the most accessible approximations to deal with the quantum complexity, its validity in time out of equilibrium has to be checked on a case by case basis. By studying the dynamics of the variance, we see that the mean field approximation can break down at surprisingly early times scaling with the square root of N. The underlying mechanism in the fully connected Ising model is identified as a dephasing effect between nearby effective orbits. The dynamics of the magnetization variance discriminates four qualitatively different regions in the dynamical phase diagram that are indistinguishable by just looking at the mean magnetization. Second, bipartite entanglement between disjoint groups of the spins in the fully connected Ising model is studied. In the large N limit the entanglement Hamiltonian is shown to be a harmonic oscillator. A quantitative relation between the angular frequency of the oscillator and a squeezing measure of the spin state is derived. Explicit expressions for the Renyi entanglement entropies are given. In contrast to the variance, which vanishes in the thermodynamic limit, the entanglement entropy saturate to an N-independent positive value. The dynamics of the entanglement Hamiltonian and the entanglement entropies is discussed. Linear growth as a function of time is linked to the unstable fixed point and the homoclinic orbit of a pitchfork bifurcation, while logarithmic increase is explained due to a dephasing mechanism of nearby effective orbits. A generalization to other models that permit a semiclassical effective description (on a possibly higher dimensional phase space) is given. As a byproduct the connection between entangled Gaussian states and the symplectic eigenvalues of its reduced covariance is derived. An upper bound on the von Neumann entanglement entropy of any bipartite pure state is given in terms of the symplectic properties of its reduced covariance. Third, the existence of dynamical phase transitions (DPTs) in the fully connected Ising model is argued by two different approaches. DPTs are defined as non-analyticities in the time-dependence of the return probability rate function in the thermodynamic limit. Like entanglement entropy, the behavior of the return rate function is an observable-independent quantity used to characterize non-equilibrium states. A numerically based connection between the DPT and the analytic property of the infinite time averaged order parameter as a function of quench strength is discussed. The concept of Feynman history states is used to devise the notion of history entanglement. An information theoretic interpretation of the history entanglement entropy as the precision of a quantum clock follows from Holevo’s bound. The dependence of this entanglement entropy on the observation time constitutes an operator-independent characterization of the quantum dynamics in closed systems. For observation times comparable to the Heisenberg time, the history entanglement becomes sensitive to the statistics of the energy gaps, and thus distinguishes integrable and non-integrable dynamics. This is demonstrated for two examples, a random matrix ensemble, and a one-dimensional hardcore boson model with integrability breaking next nearest neighbor interaction.