This thesis investigates the dynamical properties of equilibrium and nonequilibrium systems, both quantum and classical, under the guise of a thermodynamic formalism. Large deviation functions associated with the generating functions of time-integrated observables play the role of dynamical free energies and thus determine the trajectory phase structure of a system. The 1d Glauber-Ising chain is studied using the time-integrated energy as the dynamical order parameter and a whole curve of second order trajectory transitions are uncovered in the complex counting field plane. The leading dynamical Lee-Yang zeros of the associated generating function are extracted directly from the time dependent high order cumulants. Resolving the cumulants into constituent contributions the motion of each contribution’s leading Lee-Yang zeros pair allows one to infer the positions of the trajectory transition points. Contrastingly if one uses the full cumulants only the positions of those closest to the origin, in the limit of low temperatures, can be inferred. Motivated by homodyne detection schemes this thermodynamic approach to trajectories is extended to the quadrature trajectories of light emitted from open quantum systems. Using this dynamical observable the trajectory phases of a simple “blinking” 3-level system, two weakly coupled 2-level systems and the micromaser are studied. The trajectory phases of this observable are found to either carry as much information as the photon emission trajectories or in some cases capture extra dynamically features of the system (the second example). Finally, the statistics of the time-integrated longitudinal and transverse magnetization in the 1d transverse field quantum Ising model are explored. In both cases no large deviation function exists but the generating functions are still calculable. From the singularities of these generating functions new transition lines emerge. These were shown to be linked to: (a) the survival probability of an associated open system, (b) PT-symmetry, (c) the temporal scaling of the cumulants and (d) the topology of an associated set of states.
This thesis investigates the dynamical properties of equilibrium and nonequilibrium systems, both quantum and classical, under the guise of a thermodynamic formalism. Large deviation functions associated with the generating functions of time-integrated observables play the role of dynamical free energies and thus determine the trajectory phase structure of a system. The 1d Glauber-Ising chain is studied using the time-integrated energy as the dynamical order parameter and a whole curve of second order trajectory transitions are uncovered in the complex counting field plane. The leading dynamical Lee-Yang zeros of the associated generating function are extracted directly from the time dependent high order cumulants. Resolving the cumulants into constituent contributions the motion of each contribution’s leading Lee-Yang zeros pair allows one to infer the positions of the trajectory transition points. Contrastingly if one uses the full cumulants only the positions of those closest to the origin, in the limit of low temperatures, can be inferred. Motivated by homodyne detection schemes this thermodynamic approach to trajectories is extended to the quadrature trajectories of light emitted from open quantum systems. Using this dynamical observable the trajectory phases of a simple “blinking” 3-level system, two weakly coupled 2-level systems and the micromaser are studied. The trajectory phases of this observable are found to either carry as much information as the photon emission trajectories or in some cases capture extra dynamically features of the system (the second example). Finally, the statistics of the time-integrated longitudinal and transverse magnetization in the 1d transverse field quantum Ising model are explored. In both cases no large deviation function exists but the generating functions are still calculable. From the singularities of these generating functions new transition lines emerge. These were shown to be linked to: (a) the survival probability of an associated open system, (b) PT-symmetry, (c) the temporal scaling of the cumulants and (d) the topology of an associated set of states.
We use high-order cumulants to investigate the Lee-Yang zeros of generating functions of dynamical observables in open quantum systems. At long times the generating functions take on a large-deviation form with singularities of the associated cumulant generating functions-or dynamical free energies-signifying phase transitions in the ensemble of dynamical trajectories. We consider a driven three-level system as well as the dissipative Ising model. Both systems exhibit dynamical intermittency in the statistics of quantum jumps. From the short-time behavior of the dynamical Lee-Yang zeros, we identify critical values of the counting field which we attribute to the observed intermittency and dynamical phase coexistence. Furthermore, for the dissipative Ising model we construct a trajectory phase diagram and estimate the value of the transverse field where the stationary state changes from being ferromagnetic (inactive) to paramagnetic (active).
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.
The large-deviation method can be used to study the measurement trajectories of open quantum systems. For optical arrangements this formalism allows to describe the long time properties of the (nonequilibrium) photon counting statistics in the context of a (equilibrium) thermodynamic approach defined in terms of dynamical phases and transitions between them in the trajectory space [J. P. Garrahan and I. Lesanovsky, Phys. Rev. Lett. 104, 160601 (2010)]. In this paper, we study the thermodynamic approach for fluorescent systems coupled to complex reservoirs that induce stochastic fluctuations in their dynamical parameters. In a fast modulation limit the thermodynamics corresponds to that of a Markovian two-level system. In a slow modulation limit, the thermodynamic properties are equivalent to those of a finite system that in an infinite-size limit is characterized by a first-order transition. The dynamical phases correspond to different intensity regimes, while the size of the system is measured by the transition rate of the bath fluctuations. As a function of a dimensionless intensive variable, the first and second derivatives of the thermodynamic potential develop an abrupt change and a narrow peak, respectively. Their scaling properties are consistent with a double-Gaussian probability distribution of the associated extensive variable.
The dynamical behavior of many-body systems is often richer than what can be anticipated from their static properties. Here we show that in closed quantum systems this becomes evident by considering the statistics of time-integrated observables. In particular, the analytic properties of their generating functions, as estimated by full counting statistics (FCS), allow one to identify FCS phases, i.e., phases with specific fluctuation properties of time-integrated observables, and to locate transitions between these phases. We discuss in detail the case of the quantum Ising chain in a transverse field. We show that this model displays a continuum of full counting statistics transitions, of which the static transition is just an end point. These singularities are not a consequence of particular choices of initial conditions or other external nonequilibrium protocols such as quenches in coupling constants. They can be probed generically through quantum jump statistics of an associated open problem, and for the case of the quantum Ising chain we outline a possible experimental realization of this scheme by digital quantum simulation with cold ions.
We apply the large-deviation method to study trajectories in dissipative quantum systems. We show that in the long time limit the statistics of quantum jumps can be understood from thermodynamic arguments in terms of dynamical phases and transitions between them in trajectory space. We illustrate our approach with three simple examples: a driven 2-level system where we find a particular scale invariance point in the ensemble of trajectories of emitted photons; a blinking 3-level system, where we argue that intermittency in the photon count is related to a crossover between distinct dynamical phases; and a micromaser, where we find an actual first-order phase transition in the ensemble of trajectories.
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.