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.
Conserved-charge densities are very special observables in quantum many-body systems as, by construction, they encode information about the dynamics. Therefore, their evolution is expected to be of much simpler interpretation than that of generic observables and to return universal information on the state of the system at any given time. Here, we study the dynamics of the fluctuations of conserved U ( 1 ) charges in systems that are prepared in charge-asymmetric initial states. We characterize the charge fluctuations in a given subsystem using the full-counting statistics of the truncated charge and the quantum entanglement between the subsystem and the rest resolved to the symmetry sectors of the charge. We show that, even though the initial states considered are homogeneous in space, the charge fluctuations generate an effective inhomogeneity due to the charge-asymmetric nature of the initial states. We use this observation to map the problem into that of charge fluctuations on inhomogeneous, charge-symmetric states and treat it using a recently developed space-time duality approach. Specializing the treatment to interacting integrable systems we combine the space-time duality approach with generalized hydrodynamics to find explicit predictions. Published by the American Physical Society 2024
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.
Owing to its probabilistic nature, a measurement process in quantum mechanics produces a distribution of possible outcomes. This distribution-or its Fourier transform known as full counting statistics (FCS)-contains much more information than say the mean value of the measured observable, and accessing it is sometimes the only way to obtain relevant information about the system. In fact, the FCS is the limit of an even more general family of observables-the charged moments-that characterize how quantum entanglement is split in different symmetry sectors in the presence of a global symmetry. Here we consider the evolution of the FCS and of the charged moments of a U(1) charge truncated to a finite region after a global quantum quench. For large scales these quantities take a simple large-deviation form, showing two different regimes as functions of time: while for times much larger than the size of the region they approach a stationary value set by the local equilibrium state, for times shorter than region size they show a nontrivial dependence on time. We show that, whenever the initial state is also U(1) symmetric, the leading order in time of FCS and charged moments in the out-of-equilibrium regime can be determined by means of a space-time duality. Namely, it coincides with the stationary value in the system where the roles of time and space are exchanged. We use this observation to find some general properties of FCS and charged moments out of equilibrium, and to derive an exact expression for these quantities in interacting integrable models. We test this expression against exact results in the Rule 54 quantum cellular automaton and exact numerics in the XXZ spin-1/2 chain.
We study the entanglement dynamics generated by quantum quenches in the quantum cellular automaton Rule 54. We consider the evolution from a recently introduced class of solvable initial states. States in this class relax (locally) to a one-parameter family of Gibbs states and the thermalisation dynamics of local observables can be characterised exactly by means of an evolution in space. Here we show that the latter approach also gives access to the entanglement dynamics and derive exact formulas describing the asymptotic linear growth of all Rényi entropies in the thermodynamic limit and their eventual saturation for finite subsystems. While in the case of von Neumann entropy we recover exactly the predictions of the quasiparticle picture, we find no physically meaningful quasiparticle description for other Rényi entropies. Our results apply to both homogeneous and inhomogeneous quenches.
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.
Re acute accent nyi entropies are conceptually valuable and experimentally relevant generalizations of the celebrated von Neumann entanglement entropy. After a quantum quench in a clean quantum many-body system they generically display a universal linear growth in time followed by saturation. While a finite subsystem is essentially at local equilibrium when the entanglement saturates, it is genuinely out of equilibrium in the growth phase. In particular, the slope of the growth carries vital information on the nature of the system's dynamics, and its characterization is a key objective of current research. Here we show that the slope of Re acute accent nyi entropies can be determined by means of a spacetime duality transformation. In essence, we argue that the slope coincides with the stationary density of entropy of the model obtained by exchanging the roles of space and time. Therefore, very surprisingly, the slope of the entanglement is expressed as an equilibrium quantity. We use this observation to find an explicit exact formula for the slope of Re acute accent nyi entropies in all integrable models treatable by thermodynamic Bethe ansatz and evolving from integrable initial states. Interestingly, this formula can be understood in terms of a quasiparticle picture only in the von Neumann limit.