This thesis reports the results obtained during my PhD research in the field of out of equilibrium quantum many-body systems. Chapter 1 consists in a brief introduction of the field and the introduction of concept that are useful for the following chapters. In Chapter 2 the statistics of the work as a tool for characterizing the dynamics of many-body quantum systems is introduced its general features discussed. Then, such a statistics is computed for generic time-dependent protocols (both global and local) in the quantum Ising chain and in the Gaussian field theory, showing, in particular, that in its low-energy part there are features that are independent of the details of the specific chosen protocol. Chapter 3 is devoted to the study of the dynamical phase transition in the $O(N)$ quantum vector model in the $N \rightarrow \infty$ limit, whose critical properties in generic dimensions are characterized. Moreover, a strong connection between such a transition and the statistics of excitations produced in a double quench as a function of the waiting time is showed. The chapter ends by studying the fate of the dynamical transition and the its critical properties when a ramp of finite duration $\tau$ is applied to the system instead of a sudden quench. In particular, we will show that when $\tau \rightarrow \infty$ the critical point tends to the equilibrium critical point (at zero temperature) in a power-law fashion and that for every finite $\tau$ the critical properties are always the same (and different from the equilibrium critical properties). Finally in Chapter 4 we will discuss the emergence of a non adiabatic behavior in the dynamics of the order parameter for a low dimensional quantum system driven within a gapped phase by considering in detail the case of a quantum Ising chain subject to a linear variation in time of the transverse field, showing that, no matter how slowly the ramp is performed, such a change leads eventually to the disruption of the order.
This thesis reports the results obtained during my PhD research in the field of out of equilibrium quantum many-body systems. Chapter 1 consists in a brief introduction of the field and the introduction of concept that are useful for the following chapters. In Chapter 2 the statistics of the work as a tool for characterizing the dynamics of many-body quantum systems is introduced its general features discussed. Then, such a statistics is computed for generic time-dependent protocols (both global and local) in the quantum Ising chain and in the Gaussian field theory, showing, in particular, that in its low-energy part there are features that are independent of the details of the specific chosen protocol. Chapter 3 is devoted to the study of the dynamical phase transition in the O(N) quantum vector model in the N→∞ limit, whose critical properties in generic dimensions are characterized. Moreover, a strong connection between such a transition and the statistics of excitations produced in a double quench as a function of the waiting time is showed. The chapter ends by studying the fate of the dynamical transition and the its critical properties when a ramp of finite duration τ is applied to the system instead of a sudden quench. In particular, we will show that when τ→∞ the critical point tends to the equilibrium critical point (at zero temperature) in a power-law fashion and that for every finite τ the critical properties are always the same (and different from the equilibrium critical properties). Finally in Chapter 4 we will discuss the emergence of a non adiabatic behavior in the dynamics of the order parameter for a low dimensional quantum system driven within a gapped phase by considering in detail the case of a quantum Ising chain subject to a linear variation in time of the transverse field, showing that, no matter how slowly the ramp is performed, such a change leads eventually to the disruption of the order.
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.
The theory of statistical mechanics provides a powerful conceptual framework within which the relevant (macroscopic) features of systems at equilibrium can be described. As there is currently no equivalent capable of encompassing the much richer class of non-equilibrium phenomena, research in this direction proceeds mainly on an instance-by-instance basis. The aim of this Thesis is to describe in some detail three such attempts, which involve different dynamical aspects of classical and quantum systems. As summarised below, each of the last three Chapters of this document delves into one of these different topics, while Chapter 2 provides a brief introduction on the study of non-equilibrium dynamics. In Chapter 3 we investigate the purely relaxational dynamics of classical critical ferromagnetic systems in the proximity of surfaces, paying particular attention to the effects that the latter induce on the early stages of the evolution following an abrupt change in the temperature of the sample. When the latter ends close enough to the critical value which separates the paramagnetic from the ferromagnetic phase, it effectively introduces a temporal boundary which can be treated as if it were a surface. Within this picture, we highlight the emergence of novel effects near the effective edge formed by the intersection of the two spatial and temporal boundaries. Our findings are apparently in disagreement with previous predictions which were based on the assumption that the presence of such an edge would not affect the scaling behaviour of observables; in order to explain this discrepancy, we propose an alternative for the original power-counting argument which, at least, correctly predicts the emergence of novel field-theoretical divergences in our one-loop calculations. We show that said singularities are associated with the scaling at the edge. Moreover, by encoding our findings in a boundary renormalisation group framework, we argue that the new predicted behaviour represents a universal feature associated to the short-distance expansion of the order parameter of the transition near the edge; we also calculate explicitly its anomalous dimension at the first-order in a dimensional expansion. As a qualitative feature, this anomalous dimension depends on the type of phase transition occurring at the surface. We exploit this fact in order to provide numerical support to our predictions via Monte Carlo simulations of the dynamical behaviour of a three-dimensional Ising model. The main results reported in Chap. 3 have appeared in Ref. [1]. In Chapter 4 we revisit the Euclidean mapping to imaginary times which has been recently proposed [2, 3] as an alternative for approaching the problem of quantum dynamics following a quench. This is expected to allow one to reformulate the original problem as a static one confined in a film geometry. We show that this interpretation actually holds only if the initial state of the dynamics is pure. Statistical mixtures, instead, intertwine the effects due to the two boundaries, which therefore cannot be regarded as being independent. We emphasize that, although the aforementioned reinterpretation as a confined static problem fails, one is still able, in principle, to write down and solve the corresponding equations. We also discuss in some detail the relation between this approach and the real-time field-theoretical one which makes use of the two-time Keldysh contour. For this purpose, we study the analytical structure of relevant observables — such as correlation functions — in the complex plane of times, identifying a subdivision of this domain into several sectors which depend on the ordering of the imaginary parts of the involved time coordinates. Within each of these subdomains, the analytic continuation to the real axis provides in principle a different result. This feature allows one to reconstruct from the Euclidean formalism all possible non-time-ordered functions, which in particular include all those which can be calculated via the Keldysh two-time formalism. Moreover, we give a prescription on how to retrieve response functions, discussing some simple examples and rationalising some recent numerical data obtained for one of these observables in a one-dimensional quantum Ising chain [4]. We also highlight the emergence of a light-cone effect fairly similar to the one previously found for correlation functions [2], which therefore provides further confirmation to the fact that information travels across the system in the form of the entanglement of quasi-particles produced by the quenching procedure. We have reported part of this analysis in Ref. [5]. Chapter 5 presents part of our recent work on effective relaxation in quantum systems following a quench and on the observed prethermalisation. We analyse the effects caused by the introduction of a long-range integrability-breaking interaction in the early stages of the dynamics of an otherwise integrable quantum spin chain following a quench in the magnetic field. By employing a suitable transformation, we redefine the theory in terms of a fully-connected model of hard-core bosons, which allows us to exploit the (generically) low density of excitations for rendering our model exactly solvable (in a numerical sense, i.e., by numerically diagonalising an exact matrix). We verify that, indeed, as long as the parameters of the quench are not too close to the critical point, the low-density approximation captures the dynamical features of the elementary operators, highlighting the appearance of marked plateaux in their dynamics, which we reinterpret as the emergence of a prethermal regime in the original model. As expected, the latter behaviour is reflected also on extensive observables which can be constructed as appropriate combinations of the mode populations. For these quantities, the typical approach to the quasi-stationary value is algebraic with exponent a ≈ 3, independently of the size of the system, the strength of the interaction and the amplitude of the magnetic field (as long as it is kept far from the critical point). The plateaux mentioned above last until a recurrence time — which can be approximately identified with tR ≈ N/2 for single modes and t′R ≈ N/4 for extensive quantities — after which quantum oscillations due to the finite size of the chain reappear. Our procedure allows us to shed some light over prethermal features without having to considerably limit the size of the system, which we can choose to be quite large, as we discuss in Ref. [6].
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.
in particular, the order parameter profile and, also comparing with the results of recent numerical simulations, we will draw conclusions about the mass of the topological quantum particle.After these applications to classical spin systems in which the field theoretical time is taken imaginary, we will finally consider a situation in which it is kept real, namely that of non-equilibrium quantum dynamics in isolated systems.In this case the most interesting question is about the fate of the non-equilibrium evolution at large times, and it is this limit that makes low-energy excitations dominant.We will formulate the theory of quantum quenches -both homogeneous and inhomogeneous -in generic space dimensionality and will show under which conditions one-point functions of local operators (e.g. the order parameter) can exhibit oscillations that remain undamped at large times.We will also more generally study the dependence of the dynamics on the initial conditions and will exhibit the mechanisms allowing universal properties to emerge at late times.The thesis is mainly based on the results obtained in the following research articles:
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.
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.