The emergence of statistical ensembles from the unitary evolution of isolated quantum systems is one of the most fundamental conundrums of Statistical Physics. One of the main avenues for its investigation is the out-of-equilibrium dynamics in many-body 1D quantum systems, and in particular the study of thermalization in generic models, or relaxation to a generalized Gibbs ensemble in integrable models. We consider quenches in the XXZ model. Our analysis relies on the thermodynamic and algebraic Bethe ansatz, generalized hydrodynamics, and the quasiparticle picture. Until recently, in integrable quenches, low-entangled initial states have been the primary choice. However, highly entangled initial states are attracting growing interest. We define a family of states with increasing entanglement while preserving macroscopic homogeneity. They are inspired by the rainbow state and we call them 2n-rainbow states owing to their 2n-site periodicity. A relatively recent conjecture states that, in the XXX model, random states exhibit local stationary states approaching infinite temperature as their unit cell grows. Given the similarity between such random states and the 2n-rainbow states, the central question addressed in our work is: does the 2n-rainbow state, as n grows, thermalize to infinite temperature in the XXZ model? We provide analytical and numerical evidence supporting this scenario. Moreover, we obtain some nice results concerning the entanglement dynamics.
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.
In recent years the out-of-equilibrium dynamics of many-body quantum systems has attracted growing interest, both for its fundamental implications in the understanding of quantum thermalization and for potential applications in quantum information and quantum computation. In particular, quantum quenches provide a natural setting to investigate the mechanisms of correlation propagation and entanglement growth far from equilibrium. Interest in these phenomena stems not only from the possibility of testing theoretical predictions coming from emergent hydrodynamic descriptions, but also from understanding whether, and to what extent, such descriptions remain valid in the presence of spatial gradients or interfaces. Two paradigmatic frameworks in this context are the Generalized Hydrodynamics (GHD) and quasi-particle picture. The former is an effective theory devised to describe the long-time evolution after a quench in the thermodynamic limit. Within this perspective, inhomogeneous quenches are a crucial testbed for extending hydrodynamic theories to more realistic situations where, for example, scattering at an interface can occur. This is the main goal of the present thesis. The reference physical model is the one-dimensional Kitaev chain, an integrable spinless p-wave superconductor that exhibits both a topological and a trivial phase. The bilinear structure of the Hamiltonian in the fermionic creation and annihilation operators admits an exact rewriting in terms of Majorana fermions, in which the Hamiltonian remains quadratic. This representation allows a compact description of the ground-state properties via the Majorana correlation matrix (Gamma) at time t=0: the entries of Gamma completely encode the quantum information of the Gaussian ground state. The time evolution of Gamma is studied in the Heisenberg picture, where unitary dynamics correspond to an orthogonal rotation of Gamma implemented by an orthogonal matrix determined by the Kitaev Hamiltonian in the Majorana basis. Consequently, the entanglement entropy dynamics following a quench, both homogeneous and inhomogeneous, can be simulated through simple matrix products, making the approach numerically efficient even for sizable systems. The evolution code for Gamma was validated on a homogeneous quench with periodic boundary conditions: subsystem-size (l) curves exhibit the expected finite-size scaling, (Δ \propto 1/l), and numerical growth rates in the linear regime (t < v/2l), extracted by linear fits at fixed l, match the quasi-particle predictions for various values of the new chemical potential. The rate discrepancy behaves asymptotically as (Δ_{rate}} \propto 1/t_{\min}\) for large (t_{\min}\), up to oscillations from sub-extensive terms. The physical problem of interest is then an inhomogeneous quench of the chemical potential on the infinite Kitaev chain ground state with open boundaries, (μ_0 to μ_L, μ_R), producing a central interface. Initial and final Hamiltonians are free in the quasiparticle sense (ballistic propagation), but the presence of two distinct energy bands implies interface scattering. Simulations varying only μ_R (within the same trivial phase) show a nonmonotonic dependence of the entanglement growth rate, with a maximum at a particular μ_R. To explain this, an analytic ansatz for the stationary Majorana wavefunction was constructed: plane waves multiplied by a Majorana spinor, with bulk spinors obtained from the homogeneous PBC problem via the time-independent Schrödinger equation. Imposing matching at the interface yields four nontrivial equations for two unknowns (overdetermined in general). In the special cases (Δ ̸= t) two equations become trivial, the system is determined, and one obtains an explicit transmission coefficient. Combining this result with known quasi-particle formulas for two semi-infinite chains yields a hydrodynamic prediction for entanglement growth with interface scattering. Energy conservation across the interface selects the momenta that can be transmitted and thus contribute to inter-region correlations. Numerical tests for several total lengths (L) confirm the finite-size scaling (Δ_{\mathrm{e.e.}} \propto 1/L\), and show excellent agreement between numerical and theoretical growth rates for both quenches restricted to the right half and quenches applied to both halves. For the right-only case the rate discrepancy again follows (Δ_{rate} \propto 1/t_{min}\) as
tmin→∞. Possible future directions include: the study of the general inhomogeneous quench with (Δ ̸= t) via a direct extension of the analytic ansatz; the analysis of quenches starting from the topological ground state; the investigation of universal effects near the quantum critical point (μ=2t); the introduction of multiple interfaces or regions with spatially varying parameters to explore interference phenomena between quasiparticle fronts (e.g. periodic alternation (μ_1,μ_2,μ_1,μ_2,....) and, finally, the extension to interacting models where the Gaussian formalism is no longer applicable.
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.
A bstract The dynamics of quantum entanglement plays a central role in explaining the emergence of thermal equilibrium in isolated many-body systems. However, entanglement is notoriously hard to measure, and can in fact be “forged”: recent works have introduced a notion of pseudoentanglement describing ensembles of many-body states that, while only weakly entangled, cannot be efficiently distinguished from states with much higher entanglement, such as random states in the Hilbert space. This prompts the question: how much entanglement is truly necessary to achieve thermal equilibrium in quantum systems? In this work we address this question by introducing random circuit models of quantum dynamics that, at late times, equilibrate to pseudoentangled ensembles — a phenomenon we name ensemble pseudothermalization . These models replicate all the efficiently observable predictions of thermal equilibrium, while generating only a small (and tunable) amount of entanglement, thus deviating from the “maximum-entropy principle” that underpins thermodynamics. We examine (i) how a pseudoentangled ensemble on a small subsystem spreads to the whole system as a function of time, and (ii) how a pseudoentangled ensemble can be generated from an initial product state. We map the above problems onto a family of classical Markov chains on subsets of the computational basis. The mixing times of such Markov chains are related to the time scales at which the states produced from the dynamics become indistinguishable from Haar-random states at the level of each statistical moment, or number of copies. Based on a combination of rigorous bounds and conjectures supported by numerics, we argue that each Markov chain’s relaxation time and mixing time have different asymptotic behavior in the limit of large system size. This is a necessary condition for a cutoff phenomenon : an abrupt dynamical transition to equilibrium. We thus conjecture that our random circuits give rise to asymptotically sharp distinguishability transitions.
Understanding how complex entanglement structures emerge is a central problem in quantum many-body physics. Recent work by Zhang et al. [Phys. Rev. Lett. 137, 020404 (2026)] has considered structured initial states prepared by evolving a product state under a chaotic Hamiltonian for a finite time before quenching to the target Hamiltonian. In this setup, total entanglement entropy growth in many-body localized systems exhibits two distinct regimes, first increasing and then decreasing as the initial entanglement is tuned. In this work, we identify the physical origin of this behavior by analyzing the dynamics of both the Rényi entanglement entropy and the Wehrl-Rényi entropy in the random-field XXZ model, the latter of which characterizes multipartite entanglement. We show that a similar non-monotonic dependence on the initial entanglement also appears in the net growth of the Wehrl-Rényi entropy for product states polarized along the
z-direction. The first regime is governed by a finite magnetization associated with local integrals of motion, while the second reflects inter-site correlations. In contrast, for product states in the
x/y-direction, the entanglement growth exhibits a monotonic decay. Our results provide a more fine-grained picture of how distinct initial-state properties shape entanglement dynamics in many-body localized 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.