We consider the quantum XY model and study the effects of interacting perturbations on the time evolution of the von Neumann and Rényi entropies of spin blocks after global quenches. We show that the entropies are sensitive to perturbations that break hidden symmetries behind the integrability of the model. At times much larger than the characteristic time of the well-known linear increase of the entropies, we identify a time window characterized by a novel linear growth followed by saturation. The typical time of the phenomenon is inversely proportional to the perturbation strength and the behavior is trigger off by the extinction of an infinite number of local conservation laws following a non-abelian algebra. The universality of the crossover is revealed by a semi-classical picture that captures the leading behavior of the entropies. We check our theoretical predictions against iTEBD simulations. The good agreement between theory and numerics substantiates the method developed in [Bertini and Fagotti, J. Stat. Mech. (2015) P07012] for investigating a pre-relaxation limit in weakly interacting models.
We consider the quantum XY model and study the effects of interacting perturbations on the time evolution of the von Neumann and Rényi entropies of spin blocks after global quenches. We show that the entropies are sensitive to perturbations that break hidden symmetries behind the integrability of the model. At times much larger than the characteristic time of the well-known linear increase of the entropies, we identify a time window characterized by a novel linear growth followed by saturation. The typical time of the phenomenon is inversely proportional to the perturbation strength and the behavior is trigger off by the extinction of an infinite number of local conservation laws following a non-abelian algebra. The universality of the crossover is revealed by a semi-classical picture that captures the leading behavior of the entropies. We check our theoretical predictions against iTEBD simulations. The good agreement between theory and numerics substantiates the method developed in [Bertini and Fagotti, J. Stat. Mech. (2015) P07012] for investigating a pre-relaxation limit in weakly interacting models.
We investigate the crossover of the entanglement entropy toward its thermal value in nearly integrable systems. We employ equations-of-motion techniques to study the entanglement dynamics in a lattice model of weakly interacting spinless fermions after a quantum quench. For weak enough interactions we observe a two-step relaxation of the entanglement entropies of finite subsystems. Initially, the entropies follow a nearly integrable evolution, approaching the value predicted by the generalized Gibbs ensemble (GGE) of the unperturbed model. Then, they start a slow drift toward the thermal stationary value described by a standard Gibbs ensemble (GE). While the initial relaxation to the GGE is independent of the interaction, the slow drift from GGE to GE values happens on timescales proportional to the inverse interaction squared. For asymptotically large times and subsystem sizes the dynamics of the entropies can be predicted using a modified quasiparticle picture that keeps track of the evolution of the fermionic occupations caused by the integrability breaking. This picture gives a quantitative description of the results as long as the integrability-breaking timescale is much larger than the one associated with the (quasi)saturation to the GGE. In the opposite limit, the quasiparticle picture still provides the correct late-time behavior, but it underestimates the initial slope of the entanglement entropy.
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 the time evolution following a quantum quench in spin-1/2 chains. It is well known that local conservation laws constrain the dynamics and, eventually, the stationary behavior of local observables. We show that some widely studied models, like the quantum XY model, possess extra families of local conservation laws in addition to the translation invariant ones. As a consequence, the additional charges must be included in the generalized Gibbs ensemble that describes the stationary properties. The effects go well beyond a simple redefinition of the stationary state. The time evolution of a non-translation invariant state under a (translation invariant) Hamiltonian with a perturbation that weakly breaks the hidden symmetries underlying the extra conservation laws exhibits pre-relaxation. In addition, in the limit of small perturbation, the time evolution following pre-relaxation can be described by means of a time-dependent generalized Gibbs ensemble.
We consider time evolution in models close to integrable points with hidden symmetries that generate infinitely many local conservation laws that do not commute with one another. The system is expected to (locally) relax to a thermal ensemble if integrability is broken, or to a so-called generalised Gibbs ensemble if unbroken. In some circumstances expectation values exhibit quasi-stationary behaviour long before their typical relaxation time. For integrability-breaking perturbations, these are also called pre-thermalisation plateaux, and emerge e.g. in the strong coupling limit of the Bose–Hubbard model. As a result of the hidden symmetries, quasi-stationarity appears also in integrable models, for example in the Ising limit of the XXZ model. We investigate a weak coupling limit, identify a time window in which the effects of the perturbations become significant and solve the time evolution through a mean-field mapping. As an explicit example we study the XYZ spin- chain with additional perturbations that break integrability. One of the most intriguing results of the analysis is the appearance of persistent oscillatory behaviour. To unravel its origin, we study in detail a toy model: the transverse-field Ising chain with an additional nonlocal interaction proportional to the square of the transverse spin per unit length (2013 Phys. Rev. Lett . 111 197203). Despite being nonlocal, this belongs to a class of models that emerge as intermediate steps of the mean-field mapping and shares many dynamical properties with the weakly interacting models under consideration.
We consider the time evolution following a quantum quench in spin-1/2 chains. It is well known that local conservation laws constrain the dynamics and, eventually, the stationary behavior of local observables. We show that some widely studied models, such as the quantum XY model, possess extra families of local conservation laws in addition to the translation invariant ones. As a consequence, the additional charges must be included in the generalized Gibbs ensemble that describes the stationary properties. The effects go well beyond a simple redefinition of the stationary state. The time evolution of a non-translation-invariant state under a (translation-invariant) Hamiltonian with a perturbation that weakly breaks the hidden symmetries underlying the extra conservation laws exhibits pre-relaxation. In addition, in the limit of small perturbation, the time evolution following pre-relaxation can be described by means of a time-dependent generalized Gibbs ensemble.
Quantum Ising model in one dimension is an exactly solvable example of a quantum phase transition. We investigate its behavior during a quench caused by a gradual turning off of the transverse bias field. The system is then driven at a fixed rate characterized by the quench time ${\ensuremath{\tau}}_{Q}$ across the critical point from a paramagnetic to ferromagnetic phase. In agreement with Kibble-Zurek mechanism (which recognizes that evolution is approximately adiabatic far away, but becomes approximately impulse sufficiently near the critical point), quantum state of the system after the transition exhibits a characteristic correlation length $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$ proportional to the square root of the quench time ${\ensuremath{\tau}}_{Q}$: $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}=\sqrt{{\ensuremath{\tau}}_{Q}}$. The inverse of this correlation length is known to determine average density of defects (e.g., kinks) after the transition. In this paper, we show that this same $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$ controls the entropy of entanglement, e.g., entropy of a block of L spins that are entangled with the rest of the system after the transition from the paramagnetic ground state induced by the quench. For large L, this entropy saturates at $\frac{1}{6}\phantom{\rule{0.2em}{0ex}}{\mathrm{log}}_{2}\phantom{\rule{0.2em}{0ex}}\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$, as might have been expected from the Kibble-Zurek mechanism. Close to the critical point, the entropy saturates when the block size $L\ensuremath{\approx}\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$, but---in the subsequent evolution in the ferromagnetic phase---a somewhat larger length scale $l=\sqrt{{\ensuremath{\tau}}_{Q}}\phantom{\rule{0.2em}{0ex}}\mathrm{ln}\phantom{\rule{0.2em}{0ex}}{\ensuremath{\tau}}_{Q}$ develops as a result of a dephasing process that can be regarded as a quantum analog of phase ordering, and the entropy saturates when $L\ensuremath{\approx}l$. We also study the spin-spin correlation using both analytic methods and real time simulations with the Vidal algorithm. We find that at an instant when quench is crossing the critical point, ferromagnetic correlations decay exponentially with the dynamical correlation length $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$, but (as for entropy of entanglement) in the following evolution length scale l gradually develops. The correlation function becomes oscillatory at distances less than this scale. However, both the wavelength and the correlation length of these oscillations are still determined by $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{\ensuremath{\xi}}$. We also derive probability distribution for the number of kinks in a finite spin chain after the transition.