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.
Abstract The large-scale properties of homogeneous states after quantum quenches in integrable systems have been successfully described by a semiclassical picture of moving quasiparticles. Here we consider the generalisation for the entanglement evolution after an inhomogeneous quench in noninteracting systems in the framework of generalised hydrodynamics. We focus on the protocol where two semi-infinite halves are initially prepared in different states and then joined together, showing that a proper generalisation of the quasiparticle picture leads to exact quantitative predictions. If the system is initially prepared in a quasistationary state, we find that the entanglement entropy is additive and it can be computed by means of generalised hydrodynamics. Conversely, additivity is lost when the initial state is not quasistationary; yet the entanglement entropy in the large-scale limit can be exactly predicted in the quasiparticle picture, provided that the initial state is low entangled.
Abstract We investigate the non-equilibrium dynamics of the symmetry-resolved Rényi entropies in a one-dimensional gas of non-interacting spinless fermions by means of quantum generalised hydrodynamics, which recently allowed to obtain very accurate results for the total entanglement in inhomogeneous quench settings. Although our discussion is valid for any quench setting accessible with quantum generalised hydrodynamics, we focus on the case of a quantum gas initially prepared in a bipartite fashion and subsequently let evolve unitarily with a hopping Hamiltonian. For this system, we characterise the symmetry-resolved Rényi entropies as function of time t and of the entangling position x along the inhomogeneous profile. We observe an asymptotic logarithmic growth of the charged moments at half system and an asymptotic restoration of equipartition of entropy among symmetry sectors with deviations which are proportional to the square of the inverse of the total entropy.
Abstract We review the recent advances on exact results for dynamical correlation functions at large scales and related transport coefficients in interacting integrable models. We discuss Drude weights, conductivity and diffusion constants, as well as linear and nonlinear response on top of equilibrium and non-equilibrium states. We consider the problems from the complementary perspectives of the general hydrodynamic theory of many-body systems, including hydrodynamic projections, and form-factor expansions in integrable models, and show how they provide a comprehensive and consistent set of exact methods to extract large scale behaviours. Finally, we overview various applications in integrable spin chains and field theories.
Abstract. We study the problem of a quantum quench in which the initial state is the ground state of an inhomogeneous hamiltonian, in two different models, conformal field theory and ordinary free field theory, which are known to exhibit thermalisation of finite regions in the homogeneous case. We derive general expressions for the evolution of the energy flow and correlation functions, as well as the entanglement entropy in the conformal case. Comparison of the results of the two approaches in the regime of their common validity shows agreement up to a point further discussed. Unlike the thermal analogue, the evolution in our problem is non-diffusive and can be physically interpreted using an intuitive picture of quasiparticles emitted from the initial time hypersurface and propagating semiclassically.
We consider a molecular dynamics method, the so-called flea gas for computing the evolution of entanglement after inhomogeneous quantum quenches in an integrable quantum system. In such systems the evolution of local observables is described at large space-time scales by the Generalized Hydrodynamics approach, which is based on the presence of stable, ballistically propagating quasiparticles. Recently it was shown that the GHD approach can be joined with the quasiparticle picture of entanglement evolution, providing results for entanglement growth after inhomogeneous quenches. Here we apply the flea gas simulation of GHD to obtain numerical results for entanglement growth. We implement the flea gas dynamics for the gapped anisotropic Heisenberg XXZ spin chain, considering quenches from globally homogeneous and piecewise homogeneous initial states. While the flea gas method applied to the XXZ chain is not exact even in the scaling limit (in contrast to the Lieb-Liniger model), it yields a very good approximation of analytical results for entanglement growth in the cases considered. Furthermore, we obtain the full-time dynamics of the mutual information after quenches from inhomogeneous settings, for which no analytical results are available.