This thesis addresses different aspects of the dynamics of isolated quantum many-body systems far from equilibrium, ranging from transient phenomena at short times to questions about relaxation and equilibration at late times. The first part deals with nonequilibrium critical phenomena. Dynamical quantum phase transitions are characterized in a two-dimensional spin system by studying the distribution of Fisher zeros in the complex time plane. Moreover, a new topologically driven transition in the nonequilibrium steady state of quenched Chern insulators is investigated, resulting in the identification of universal nonanalytic behavior of the Hall conductance as a function of the quench parameter. The second part of the thesis comprises a study of transient dynamics of the transverse-field Ising model in two and three dimensions based on a perturbative derivation of classical networks to encode the time-evolved wave function; this construction is also applicable to represent the state as an artificial neural network. In the last part of the thesis a possible definition of the irreversibility of quantum many-body dynamics is introduced, which is based on the echo dynamics of observables under imperfect effective time reversal. Echo dynamics is studied in generic and integrable quantum systems as well as in the semiclassical limit. The relation to out-of-time-order correlators is explored, which were recently introduced to probe the quantum butterfly effect.
This thesis addresses different aspects of the dynamics of isolated quantum many-body systems far from equilibrium, ranging from transient phenomena at short times to questions about relaxation and equilibration at late times. The first part deals with nonequilibrium critical phenomena. Dynamical quantum phase transitions are characterized in a two-dimensional spin system by studying the distribution of Fisher zeros in the complex time plane. Moreover, a new topologically driven transition in the nonequilibrium steady state of quenched Chern insulators is investigated, resulting in the identification of universal nonanalytic behavior of the Hall conductance as a function of the quench parameter. The second part of the thesis comprises a study of transient dynamics of the transverse-field Ising model in two and three dimensions based on a perturbative derivation of classical networks to encode the time-evolved wave function; this construction is also applicable to represent the state as an artificial neural network. In the last part of the thesis a possible definition of the irreversibility of quantum many-body dynamics is introduced, which is based on the echo dynamics of observables under imperfect effective time reversal. Echo dynamics is studied in generic and integrable quantum systems as well as in the semiclassical limit. The relation to out-of-time-order correlators is explored, which were recently introduced to probe the quantum butterfly effect.
The work presented in this thesis is divided into two main parts: strong zero modes in quantum spin chains with open boundaries and large deviations in classical spin systems. The results concerning the strong zero modes are presented in chapters 2 and 3, while chapter 4 contains the results about large deviations in systems of many, but independent, classical degrees of freedom. These two distinct classes of problems in spin chains are both about aspects of slow relaxation to equilibrium in many-body systems. Nowadays slow dynamics and non-ergodicity in many-body quantum systems are receiving much attention due to their notable technological potential for storing and processing, for very long times, the quantum information encoded in boundary degrees of freedom. In the first part of the work presented in this thesis we mainly focus on quantum spin chains with open boundaries, which can display coherence times for edge spins, that diverge with the system size, as a consequence of the emergence, in the thermodynamic limit, of conserved local quantities, called strong zero modes. First, we discuss the fate of these coherence times when dynamics becomes dissipative. Understanding the behavior of these long correlation times in open quantum systems, undergoing Markovian (memoryless) dynamics, would allow for the possibility of exploiting the protected information encoded in the boundary spins for applications in quantum devices in realistic non-equilibrium settings. Second, we study one possible simple generalization of the transverse field Ising model to a class of ladder spin systems with plaquette interactions in the presence of a transverse field. Also for these systems one can explicitly construct the strong zero mode operators and understand the behavior of the consequential long coherence times for the edge degrees of freedom under isolated (unitary) dynamics. This is possible because the models we consider, while appearing complex superficially, can be brought to simple forms by exploiting an extensive number of conserved quantities they possess. The results presented in chapters 2 and 3 are closely related, even if they focus on different models and investigate different aspects of their dynamics. Nevertheless, the underlying concept they are based on is the same one, that is, an eigenstate phase transition for the strong zero modes that has interesting consequences for the full energy spectrum. In the second part of the thesis we focus on a different type of phase transitions. Here we address the issue of dynamical large deviation phase transitions from a disordered phase to an ordered one, in classical systems which consist in a collection of independent Ising spins evolving stochastically. Surprisingly, by reweighing trajectories accordingly to a certain type of time-integrated observable, even a system of independently evolving spins can display some critical behavior, at the fluctuation level, in trajectory space in the large size and long time limit. By making use of the tilted generator techniques we study the statistics of two particular time-integrated observables. Their associated tilted generators have an exact Kramers-Wannier duality and they look like, up to a sign and an additive constant, some specific well known quantum Hamiltonians. In this sense, the dynamical large deviation phase transitions we observe in these classical stochastic systems translate in ground state quantum phase transitions taking place at the self dual point of the tilted generators we study.
In this Thesis we discussed several issues related to the dynamics of isolated quantum many-body systems. The brief overview of Chapter 1 highlighted a series of open fundamental theoretical questions which are crucial for the investigation of nonequilibrium phenomena. Among them, we mainly focused on the emergence of thermal behaviors and on the characterization of dynamical phase transitions and of the associated universal dynamics. To this end, we analyzed some paradigmatic models, such as the Fermi-Hubbard model, the quantum Ising chain, and the O(N) model. -Chapter 2: Absence of thermalization in a Fermi liquid. In this Chapter we report two results that question the common expectation of thermalization in a Landau-Fermi liquid after an interaction quench. We first consider the perturbative expansion in the interaction strength of the long-wavelength structure factor S(q) and show that it does not satisfy the hypothesis that steady-state averages correspond to thermal ones. In particular, S(q) has an analytical component ∼ const. + O(q^2) for q → 0, compatible with thermalization, but it retains also a nonanalytic term ∼ |q| characteristic of a Fermi liquid at zero-temperature. In real space, this nonanalyticity turns in a power-law decay of the density-density correlations, in contrast with the exponential decay associated to thermalization. We next analyze the low-density case, where one can obtain results valid at any order in interaction but at leading in the density, and find that in the steady state the momentum distribution jump at the Fermi surface is strictly finite, though smaller than at equilibrium. -Chapter 3: Nonadiabatic stationary behavior in a driven Ising chain. In this Chapter we discuss the emergence of a nonadiabatic behavior in the dynamics of the order parameter in a low-dimensional system subject to a linear ramp of one of its parameters within a gapped phase, which should be the most favorable situation for an adiabatic evolution. We study in details the dynamics of the spontaneous magnetization m_x(t) in a quantum Ising chain after a linear variation in time of the transverse field within the ordered phase. In particular, focusing on the value of the order parameter at the end of the ramp m_x(τ), we show that the smaller the switching rate of the transverse field is the closer m_x(τ) gets to its ground state value. Nonetheless, this apparently small correction to adiabaticity eventually leads to the disruption of the order exponentially fast in the subsequent time evolution, no matter how slowly the ramp is performed. -Chapter 4: Aging and coarsening in isolated quantum systems after a quench. This Chapter analyzes the nonequilibrium dynamics of an isolated quantum system after a sudden quench to the dynamical critical point, where the emergence of scaling and universal exponents is expected, due to the absence of time and length scales. We explore these features for a bosonic interacting scalar field theory with O(N) symmetry in the large-N limit, where the model is exactly solvable and the exponents and scaling forms of the relevant two-time correlation functions can be analytically calculated. Moreover, we provide evidence of the emergence of a dynamical scaling behavior also for quenches below the dynamical critical point, associated with coarsening. We find that the latter case is characterized by the same scaling functions as those describing the critical case, yet with different exponents. -Chapter 5: Dynamical transitions and statistics of excitations. In this Chapter we study the dynamics of the O(N) model (introduced in the previous Chapter) resulting from a different protocol: a linear ramp of one of its parameters. We find that the presence of a dynamical phase transition, as well as its critical properties, are robust against the change of the protocol. We show that a characterization based on the critical dimensions and exponents would suggest that the dynamical phase transition is analogous to the equilibrium thermal one. However, its nonequilibrium nature becomes evident analyzing the statistics of excitations produced in the ramp process. In particular, the critical properties are encoded in the fluctuations in the number of excitations, which display qualitatively different behaviors depending on the ramp being performed above, at, or below the dynamical critical point. These behaviors bear no dependence on the duration of the protocol.
Understanding how quantum matter behaves when driven out of equilibrium is one of the key focuses in quantum physics. Thanks to impressive progress in the control and precision achieved in quantum synthetic matter over the past decades, the nonequilibrium quantum many-body physics has become one of the most active research areas today, especially after the experimental realization of Bose-Einstein condensates and optical lattices, which allows us to directly observe and study nonequilibrium quantum matter with great accuracy and controllability. In this dissertation, I explore the rich landscape of nonequilibrium quantum many-body physics and how quantum phase transitions, both symmetry-breaking and topological, can be extended to the nonequilibrium setting. In the first part of the dissertation, I focus on spinor Bose-Einstein condensates as an isolated quantum many-body system, and reveal their various dynamical behaviors, including quantum collapse and revivals, thermalization and nonthermal equilibration with no revival even though the system has finite degrees of freedom. In contrast to typical integrable systems, which usually do not thermalize, we find that spinor condensates have a parameter range in which the system thermalizes via the Eigenstate Thermalization Hypothesis (ETH). We show that this observation is linked to the presence of rare nonthermal states whose fraction vanishes with system size, and contributes to the notion of thermalization via weak ETH. Next, I explore a dynamical process that is complementary to thermalization in isolated quantum systems: information scrambling, which could be probed via out-of-time-order correlators (OTOC). I propose a nonintegrable, disordered and quasi-1D spin model, the ladder-XX model, for a feasible detection of information scrambling in a cold atom simulator. This chapter poses a fundamental question: "What are the signatures of quantum phases and phase transitions in isolated interacting systems driven out-of-equilibrium?" I study the ladder-XX model in both clean and disordered potentials, and characterize different nonequilibrium phases, i.e., ergodic and many-body localized, of the model based on the decay properties of OTOCs. Emergent light cone shows sublinear behavior, while the butterfly cones drastically differ from the light cone by demonstrating superlinear spread of information with a velocity that is bounded by the light cone velocity. In the second part of the dissertation, I continue to search for answers to the question posed above, however this time with a particular focus on symmetry-breaking and topological quantum phase transitions. I pin down a universal mechanism underlying the relation between information scrambling at any temperature and quantum phases at low temperatures. Our method points to key ingredients to dynamically detect long-range order in gapped phases through OTOCs for symmetry-breaking quantum phase transitions and Z2 topological order associated with Majorana zero modes localized at the edges. Our results pave the way to an intriguing observation that phases of quantum matter could protect the information from scrambling and thermalization, even when the system is interacting and nonintegrable. Finally, I explore and propose utilizing short-time transient temporal regimes and single-site probes to detect the phases and phase transitions in quantum matter. These studies reveal a dynamical crossover and a dynamical phase transition, respectively for periodic and open-boundary chains. In both cases, a nonequilibrium scaling law appears in the vicinity of the crossover/transition with associated exponents that differ from the analytical predictions for long times. Feasibility of detecting such dynamical criticality in experimental systems are discussed.
The study of quantum dynamics is an increasingly important aspect of atomic, molecular and optical physics, as well as condensed-matter physics and other disciplines. It is essential in order to understand non-equilibrium phenomena which do not fit within the framework of equilibrium thermodynamics, such as many-body localization, dynamical quantum phase transitions, and non-equilibrium steady states in driven-dissipative systems. However, the exact calculation of dynamical observables is extremely difficult in general for several reasons, which makes it challenging for theoretical predictions and analyses to keep pace with rapidly advancing experiments. In this thesis I use and develop several computational methods to investigate the dynamics of quantum many-body systems, with the two goals of better understanding the evolution of local observables and correlation functions in spin and Fermion systems as well as creating and testing more effective numerical methods for solving quantum dynamics problems. I explore quantum dynamics by using Hamiltonian parameter quenches in several commonly-used models in AMO physics and elsewhere: the Fermi-Hubbard model, the XXZ model, and the transverse-Ising model. The first represents a system of short-range interacting Fermions, while the latter two represent spin systems with different spin-spin couplings and magnetic field terms. All of these models can be realized in current experiments. First, I derive an exact solution for an extreme interaction quench in the Fermi-Hubbard model, from a strongly-interacting initial state to a noninteracting final Hamiltonian. From this result I show that there are nontrivial transient connected correlations, despite the fact that the system has very high temperature and a noninteracting Hamiltonian. Those features would suppress such correlations in equilibrium, but the dynamical correlations are still significant. Then, I develop and examine two numerical methods for solving dynamics of the spin expectation values and correlation functions in the transverse-Ising and XXZ models. Both methods rely on approximating a thermodynamically large system with a much smaller and more computationally tractable cluster of lattice sites, or a series of clusters. The first method is a dynamical numerical linked cluster expansion (NLCE). NLCEs have already been used successfully on systems in thermodynamic equilibrium, and I extended the technique to directly compute the dynamics of an observable after a quantum quench. I review the general NLCE method, describe my dynamical variant of it, and compute the quench dynamics of spin and correlation observables in the transverse-Ising and XXZ models. The second method, which I term the adaptive-boundary cluster method, relies on recently refined bounds on information propagation in quantum systems to choose a cluster of lattice sites that minimizes the finite size error bound of a dynamical computation. I briefly review these bounds, and describe a new algorithm which uses them to construct a cluster whose boundary is chosen according to a related error estimate. Such an adaptive-boundary cluster has the property that there exists a computable upper bound on the finite size error of observables measured at the cluster center, and this error bound is nearly minimal when compared to that of other clusters of the same size. I then compute the quench dynamics of spin and correlation observables in the transverse-Ising model using the adaptive boundary clusters. I apply this method to both uniform and strongly disordered systems. To evaluate the effectiveness of both dynamical NLCE and adaptive-boundary clusters, I compare the results obtained using them to results from a standard and general purpose numerical method: exact solvers applied a rectangular cluster with periodic boundary conditions. Both of the new methods usually converge faster and are more accurate than the standard method for comparable cluster sizes. This thesis illustrates that choosing numerical methods which take advantage of the geometric and entanglement structure of the initial state and Hamiltonian can substantially enhance our ability to study quantum dynamics, by reducing the computational cost of numerically exact solutions. I conclude with a discussion of various avenues for future research to build on this work.
Quantum phase transitions occur strictly at zero temperature and are accessed by varying a physical parameter such as a magnetic field or chemical potential. While these transitions themselves are experimentally inaccessible, the existence of a zero temperature quantum critical point has remarkable effects on the dynamical properties of interacting quantum many-body systems at finite temperatures. In this thesis, we study dynamics in four strongly interacting many-body systems near quantum critical points in the context of NMR relaxation, quantum quenches, many-body quantum chaos, and charge diffusion. In chapter 1, we begin by revisiting the issue of NMR relaxation in the quantum disordered phase of Ising spin chains. We find that the imaginary part of the spin susceptibility at low frequences is dominated by rare quantum processes which leads to an NMR relaxation rate that scales as 1/T1∼exp(−T2Δ) near the quantum critical point rather than the previously conjectured 1/T1∼exp(−TΔ). These results resolve certain discrepancies between the energy scales measured with different experimental probes in the quantum disordered paramagnetic phase of the Ising chain system CoNb2O6. In chapter 2, we study the non-equilibrium dynamics of Sachdev-Ye-Kitaev (SYK) models with q body interactions after a quantum quench, considering a quench protocol which drives the system from an initial state with long-lived quasiparticle excitations to a final state lacking any long-lived quasiparticle excitations. For q=4, we find that the timescale for the system to locally thermalize τeq∼kBTℏ and for q→∞ we find that the system thermalizes instantaneously. In chapter 3, we study the real time dynamics of quantum electrodynamics in 2+1 dimensions with Nf flavors of fermions. We focus on thermalization and the onset of many-body quantum chaos and compute the quantum Lyapunov exponent in a Nf1 expansion. We briefly comment on chaos, locality, and gauge invariance. Finally, in Chapter 4 we consider whether the butterfly velocity - a speed at which quantum information propagates - may provide a fundamental bound on diffusion constants in dirty incoherent metals by computing the charge diffusion constant and the butterfly velocity in charge-neutral holographic matter with long wavelength ``hydrodynamic" disorder. In this limit, we find that the butterfly velocity does not set a sharp lower bound for the charge diffusion constant.
Thermalization in quantum many-body systems is a central problem of quantum mechanics. Despite the unitary nature of quantum evolution, isolated many-body systems often appear to relax toward states that are locally described by statistical ensembles. However, phenomena such as integrability and quantum many-body scars demonstrate that thermalization is not universal and can fail in unexpected ways, leading to nonergodic dynamics and long-lived nonequilibrium behavior. In this thesis, the problem is addressed from multiple complementary perspectives, including thermalization, integrability breaking, constrained dynamics, and quantum many-body scars, together with the development of new tensor-network methods for nonequilibrium quantum systems and quantum metrological applications. First, we study energy-filtering protocols for the preparation of thermal states, analyzing the possibility of simulating them classically through the behavior of entanglement and correlations. The effects of weak integrability breaking are then investigated in models relevant to Rydberg-atom experiments, leading to the identification of novel scattering mechanisms and the introduction of subsystem variance as a sensitive probe of weak integrability breaking. A part of the thesis is devoted to quantum many-body scars and their consequences for relaxation and transport phenomena. Their dynamical signatures are studied in both closed and open quantum systems, revealing new mechanisms that sustain long-lived nonequilibrium behavior. The thermalization processes arising during the quasi-adiabatic preparation of the ground state in the trapped Bose-Hubbard model are also examined, showing how transport constraints can stabilize low-entropy states under experimentally realistic conditions. Finally, a novel tensor-network approach based on Tree Tensor Operators is introduced for the efficient computation of the Quantum Fisher Information of generic mixed states. Applicable to both closed and open quantum systems, the method provides a promising numerical tool for quantum sensing and metrological applications.