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 \(t_{\min} \to \infty\). 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.
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 (t_{\min} \to \infty). 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.
There are many ways to characterize quantum systems, for example, by measuring system-specific observables or probing their response to external perturbations. A complementary approach is based on quantum information measures. Independent of the microscopic details, it offers a universal perspective on quantum systems and enables a unified characterization of diverse quantum systems. Central, in this regard, is the determination of the quantum entanglement within a quantum system. Recently, in the interest of studying genuine nonequilibrium properties of a quantum system, the dynamics of the quantum entanglement have been established as a powerful approach to gain deep insights into the structure of quantum matter. In this Thesis, we investigate the entanglement dynamics in one-dimensional non-interacting fermionic chains following a quantum quench, i.e. after a sudden change of the system parameters. We propose a new and generic setting for the analysis of entanglement spreading in complex inhomogeneous systems, which we call a hybrid quench-probe setup. The key feature of this new approach is a separation between the quenched part of the system and a static probe region, where the quench induced entanglement spreading is measured. This separation allows a spatial and energy-resolved analysis of the entanglement dynamics. We demonstrate the capability of our protocol by identifying fractionally quantized entanglement jumps of log(2)/2 as a signature of fractionalized topological quasi-particles, such as isolated Majorana modes, that appear at the edges of an open-boundary Kitaev chain in the topological phase. In addition, we observe in this setup an interesting nonequilibrium effect: a static and a propagating contribution to the entanglement dynamics. To understand the origin of this effect, a local perspective on the entanglement structure is necessary. For this purpose, we employ the information lattice, which offers a local scale-resolved decomposition of the correlations present in a quantum state. In nonequilibrium systems with unitary dynamics, this framework enables a hydrodynamic description of the information flow through local densities and currents. Using the information lattice, we analyze three quantum quenches in one-dimensional fermionic systems and identify spatiotemporal features of the local information flow, such as information interfaces and long-range transport from topological edge modes. This analysis enables us to derive an analytical description of the nonequilibrium effect observed in the quenched Kitaev chain.
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.
We investigate entanglement dynamics in bipartite systems governed by inhomogeneous Hamiltonians of the form H=HL+HR, where HL/R acts only on the left or right region and is homogeneous within each region. Focusing on the XX chain and the transverse-field Ising chain, we derive analytical formulas for the entanglement entropy between the two regions in the hydrodynamic limit of long times. In this regime, fermions incident on the interface undergo scattering, generating entanglement between reflected and transmitted modes. The resulting quasiparticle picture is controlled by the transmission coefficient, which we obtain analytically by solving the stationary lattice Schrödinger equation. Due to the bounded dispersion, strong inhomogeneity suppresses both transport and entanglement growth. We benchmark our analytical predictions against numerical simulations in paradigmatic setups. Finally, we extend the analysis to the interacting XXZ chain using tDMRG. The numerical data show qualitative agreement with the quadratic case: entanglement growth remains suppressed in the strongly inhomogeneous limit. Notably, however, entanglement continues to increase even when transport is suppressed, at least at intermediate times.
Quantum systems display counterintuitive features, escaping our everyday-life intuition, whose origin is related to entanglement. This phenomenon arises when different parts of the same systems display correlations, which cannot be understood via classical mechanics, and they are dubbed as quantum. To characterize the correlations above, it is necessary to find good entanglement measures which are able to capture quantitatively the quantum features of a given system. This motivates the investigation and the computation of quantities such as the von Neumann entropy of regions, which provides information about how a system is entangled. Our focus is quantum one-dimensional systems close to their critical point, where quantum correlations are enhanced, and a characterization via Quantum Field Theory (QFT) can be provided. In particular, we aim to understand the role of symmetries and localized defects on the entanglement properties. Those aspects display an intimate connection. First, via replica trick, the entropy is expressed via a replica permutation symmetry. Then, single defects break symmetries that are present in the corresponding homogeneous systems. Due to these reasons, the characterization of symmetry actions on subsystems can be viewed as a unifying framework for quantifying entanglement. We split this thesis into two parts. The first part is devoted to the symmetry-resolved entanglement. We first compute it in massive integrable field theories, giving analytical predictions for free fermions/bosons and the 3-state Potts model. For instance, we characterize the universal corrections to area law via form factor bootstrap. We also study excited states of field theories, and we characterize both low-energy states of Conformal Field Theory (CFT) and quasi-particle states. The second part deals with quantum correlations across defects. We provide a CFT framework to compute Rényi entropy and negativity in the ground state of conformal junctions, giving analytical predictions for free fermions and bosons. We find logarithmic growth of the entanglement measures above, and we relate the universal prefactor to the boundary conditions, encoded in a boundary scattering matrix. Then, we focus on the entanglement dynamics, and we study the time evolution of correlations for two paradigmatic classes of quench protocols. For a global quench, we find a linear growth of entropy across the defect whose slope is computed both in CFT, for thermal states, and for free fermions, via a modification of the quasi-particle picture. For a domain wall melting protocol, also known as inhomogeneous quench, we find that, while local observables are easily understood via hydrodynamics and local generalized Gibbs ensembles, the same is not true for the entropy. In contrast, by employing the quasi-particle picture, we provide precise predictions for the dynamics of entanglement, revealing the influence of long-range correlations between points that are specular with respect to the defect.
Recent experimental advances on ultracold atomic gases and trapped ions have made it possible to simulate exactly solvable quantum systems of interacting particles. In particular, the feasibility of making rapid changes, so-called quantum quenches, to such set-ups has allowed experimentalists to probe non-equilibrium phenomena in closed interacting quantum systems. This, in turn, has spurred a considerable theoretical interest in quantum many-body systems out of equilibrium. In this thesis, we study non-equilibrium properties of quantum many-body systems in the framework of exactly solvable quantum field theory in one spatial dimension. Specific systems include interacting fermions described by the Luttinger model and effective descriptions of spin chains using conformal field theory (CFT). Special emphasis is placed on heat and charge transport, studied from the point of view of quench dynamics, and, in particular, the effects of breaking conformal symmetries on transport properties. Examples include the Luttinger model with non-local interactions, breaking Lorentz and scale invariance, and inhomogeneous CFT, which generalizes standard CFT in that the usual propagation velocity v is replaced by a function v(x) that depends smoothly on the position x, breaking translation invariance. The quench dynamics studied here is for quantum quenches between, in general, different smooth inhomogeneous systems. An example of this is the so-called smooth-profile protocol, in which the initial state is defined by, e.g., smooth inhomogeneous profiles of inverse temperature and chemical potential, and the time evolution is governed by a homogeneous Hamiltonian. Using this protocol, we compute exact analytical results for the full time evolution of the systems mentioned above. In particular, we derive finite-time results that are universal in the sense that the same relations between the non-equilibrium dynamics and the initial profiles hold for any unitary CFT. These results also make clear that heat and charge transport in standard CFT are purely ballistic. Finally, we propose and study an inhomogeneous CFT model with v(x) given by a random function. We argue that this model naturally emerges as an effective description of one-dimensional quantum many-body systems with certain static random impurities. Using tools from wave propagation in random media, we show that such impurities lead to normal and anomalous diffusive contributions to heat transport on top of the ballistic one known from standard CFT.
Characterising the universal aspects of non-equilibrium quantum many-body dynamics is one of the key goals of this century's physics research. Progress, however, is hindered by the lack of general theoretical frameworks for studying interacting quantum matter far from equilibrium. A recent breakthrough has been the realization that several key non-equilibrium quantities, such as the rate of growth of entanglement or the fluctuations of conserved charges within finite subsystems, can be related to equilibrium properties through a space-time duality that effectively exchanges the roles of space and time. This observation effectively enables the study of non-equilibrium phenomena using tools and concepts borrowed from equilibrium statistical mechanics and thermodynamics. A first proof of principle of this framework, dubbed space-time duality approach (SDA), was provided by interacting integrable systems, where thermodynamic properties can often be characterized exactly, while dynamical quantities typically remain beyond analytical reach. Subsequent developments, however, revealed that the SDA suffered from an intrinsic ambiguity, restricting its applicability to homogeneous quenches and to charge fluctuations arising from symmetric initial states. Here we resolve this ambiguity from first principles and derive closed-form predictions for entanglement growth and charge fluctuations after general quantum quenches. We benchmark our results against the exact analytical solution of the Rule 54 quantum cellular automaton and extensive TEBD simulations of the XXZ chain. Moreover we show that, when specialised to the entanglement entropy, our framework naturally reproduces the predictions of the quasiparticle picture.