We study analytically the role of initial conditions in nonequilibrium quantum dynamics considering the one-dimensional ferromagnets in the regime of spontaneously broken symmetry. We analyze the expectation value of local operators for the infinite-dimensional space of initial conditions of domain wall type, generally intended as initial conditions spatially interpolating between two different ground states. At large times the unitary time evolution takes place inside a light cone produced by the spatial inhomogeneity of the initial condition. In the innermost part of the light cone the form of the space-time dependence is universal, in the sense that it is specified by data of the equilibrium universality class. The global limit shape in the variable x/t changes with the initial condition. In systems with more than two ground states the tuning of an interaction parameter can induce a transition which is the nonequilibrium quantum analog of the interfacial wetting transition occurring in classical systems at equilibrium. We illustrate the general results through the examples of the Ising, Potts and Ashkin-Teller chains.
We study analytically the role of initial conditions in nonequilibrium quantum dynamics considering the one-dimensional ferromagnets in the regime of spontaneously broken symmetry. We analyze the expectation value of local operators for the infinite-dimensional space of initial conditions of domain wall type, generally intended as initial conditions spatially interpolating between two different ground states. At large times the unitary time evolution takes place inside a light cone produced by the spatial inhomogeneity of the initial condition. In the innermost part of the light cone the form of the space-time dependence is universal, in the sense that it is specified by data of the equilibrium universality class. The global limit shape in the variable x/t changes with the initial condition. In systems with more than two ground states the tuning of an interaction parameter can induce a transition which is the nonequilibrium quantum analog of the interfacial wetting transition occurring in classical systems at equilibrium. We illustrate the general results through the examples of the Ising, Potts and Ashkin-Teller chains.
We consider the unitary time evolution of an interface in the regime of spontaneously broken symmetry of the two-dimensional quantum Ising model. The interface is induced by an initial condition interpolating between the two degenerate ground states in one of the spatial dimensions. The interpolation is left generic in order to investigate the dependence of the late time dynamics on the initial condition. Exploiting the basis of asymptotic quasiparticle states of the bulk theory, the order parameter is analytically determined at large times in the rough phase. The mechanism allowing the breakdown of this phase as the distance from criticality increases emerges from the theory.
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.
The nonequilibrium dynamics of classical random Ising spin chains with nonconserved magnetization are studied using an asymptotically exact real space renormalization group (RSRG). We focus on random field Ising model (RFIM) spin chains with and without a uniform applied field, as well as on Ising spin glass chains in an applied field. For the RFIM we consider a universal regime where the random field and the temperature are both much smaller than the exchange coupling. In this regime, the Imry-Ma length that sets the scale of the equilibrium correlations is large and the coarsening of domains from random initial conditions (e.g., a quench from high temperature) occurs over a wide range of length scales. The two types of domain walls that occur diffuse in opposite random potentials, of the form studied by Sinai, and domain walls annihilate when they meet. Using the RSRG we compute many universal asymptotic properties of both the nonequilibrium dynamics and the equilibrium limit. We find that the configurations of the domain walls converge rapidly toward a set of system-specific time-dependent positions that are independent of the initial conditions. Thus the behavior of this nonequilibrium system is pseudodeterministic at long times because of the broad distributions of barriers that occur on the long length scales involved. Specifically, we obtain the time dependence of the energy, the magnetization, and the distribution of domain sizes (found to be statistically independent). The equilibrium limits agree with known exact results. We obtain the exact scaling form of the two-point equal time correlation function <S0(t)Sx(t)> and the two-time autocorrelations <S0(t('))S0(t)>. We also compute the persistence properties of a single spin, of local magnetization, and of domains. The analogous quantities for the +/-J Ising spin glass in an applied field are obtained from the RFIM via a gauge transformation. In addition to these we compute the two-point two-time correlation function <S0(t)Sx(t)><S0(t('))Sx(t('))> which can in principle be measured by experiments on spin-glass-like systems. The thermal fluctuations are studied and found to be dominated by rare events; in particular all moments of truncated equal time correlations are computed. Physical properties which are typically measured in aging experiments are also studied, focusing on the response to a small magnetic field which is applied after waiting for the system to equilibrate for a time t(w). The nonequilibrium fluctuation-dissipation ratio X(t,t(w)) is computed. We find that for (t-t(w)) approximately t(alpha)(w) with alpha;<1, the ratio equal to its equilibrium value X=1, although time translational invariance does not hold in this regime. For t-t(w) approximately t(w) the ratio exhibits an aging regime with a nontrivial X=X(t/t(w)) not equal to 1, but the behavior is markedly different from mean field theory. Finally the distribution of the total magnetization and of the number of domains is computed for large finite size systems. General issues about convergence toward equilibrium and the possibilities of weakly history-dependent evolution in other random systems are discussed.
We characterize the universal far-from-equilibrium dynamics of the two-dimensional quantum Heisenberg magnet isolated from its environment. For a broad range of initial conditions, we find a long-lived universal prethermal regime characterized by self-similar behavior of spin–spin correlations. We analytically derive the spatial–temporal scaling exponents and find excellent agreement with numerics using phase space methods. The scaling exponents are insensitive to the choice of initial conditions, which include coherent and incoherent spin states with values of total magnetization and energy in a wide range. Compared to previously studied self-similar dynamics in nonequilibrium O( n ) field theories and Bose gases, we find qualitatively distinct scaling behavior originating from the presence of spin modes that remain gapless at long times and are protected by the global SU(2) symmetry. Our predictions, which suggest a distinct nonequilibrium universality class from Bose gases and O( n ) theories, are readily testable in ultracold atoms simulators of Heisenberg magnets.
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].
We study the quantum dynamics resulting from preparing a one-dimensional quantum system in the ground state of initially two decoupled parts which are then joined together (local quench). Specifically we focus on the transverse Ising chain and compute the time-dependence of the magnetization profile, m_l(t), and correlation functions at the critical point, in the ferromagnetically ordered phase and in the paramagnetic phase. At the critical point we find finite size scaling forms for the nonequilibrium magnetization and compare predictions of conformal field theory with our numerical results. In the ferromagnetic phase the magnetization profiles are well matched by our predictions from a quasi-classical calculation.