Abstract We study the dynamics of a classical disordered macroscopic model completely isolated from the environment reproducing, in a classical setting, the ‘quantum quench’ protocol. We show that, depending on the pre and post quench parameters, the system approaches equilibrium, succeeding to act as a bath on itself, or remains out of equilibrium, in two different ways. In the latter one, the system stays confined in a metastable state in which it undergoes stationary dynamics characterised by a single temperature. In the other, the system ages and its dynamics are characterised by two temperatures associated with observations made at short and long time differences (high and low frequencies). The parameter dependence of the asymptotic states is rationalised in terms of a dynamic phase diagram with one equilibrium and two out of equilibrium phases. Aspects of pre-thermalisation are observed and discussed. Similarities and differences with the dynamics of the dissipative model are also explained.
Abstract We study the dynamics of a classical disordered macroscopic model completely isolated from the environment reproducing, in a classical setting, the ‘quantum quench’ protocol. We show that, depending on the pre and post quench parameters, the system approaches equilibrium, succeeding to act as a bath on itself, or remains out of equilibrium, in two different ways. In the latter one, the system stays confined in a metastable state in which it undergoes stationary dynamics characterised by a single temperature. In the other, the system ages and its dynamics are characterised by two temperatures associated with observations made at short and long time differences (high and low frequencies). The parameter dependence of the asymptotic states is rationalised in terms of a dynamic phase diagram with one equilibrium and two out of equilibrium phases. Aspects of pre-thermalisation are observed and discussed. Similarities and differences with the dynamics of the dissipative model are also explained.
Abstract We study the Hamiltonian dynamics of the spherical spin model with fully-connected two-body random interactions. In the statistical physics framework, the potential energy is of the so-called p = 2 kind, closely linked to the scalar field theory. Most importantly for our setting, the energy conserving dynamics are equivalent to the ones of the Neumann integrable model. We take initial conditions from the Boltzmann equilibrium measure at a temperature that can be above or below the static phase transition, typical of a disordered (paramagnetic) or of an ordered (disguised ferromagnetic) equilibrium phase. We subsequently evolve the configurations with Newton dynamics dictated by a different Hamiltonian, obtained from an instantaneous global rescaling of the elements in the interaction random matrix. In the limit of infinitely many degrees of freedom, , we identify three dynamical phases depending on the parameters that characterise the initial state and the final Hamiltonian. We next set the analysis of the system with finite number of degrees of freedom in terms of N non-linearly coupled modes. We argue that in the limit the modes decouple at long times. We evaluate the mode temperatures and we relate them to the frequency-dependent effective temperature measured with the fluctuation-dissipation relation in the frequency domain, similarly to what was recently proposed for quantum integrable cases. Finally, we analyse the N − 1 integrals of motion, notably, their scaling with N , and we use them to show that the system is out of equilibrium in all phases, even for parameters that show an apparent Gibbs–Boltzmann behaviour of the global observables. We elaborate on the role played by these constants of motion after the quench and we briefly discuss the possible description of the asymptotic dynamics in terms of a generalised Gibbs ensemble.
We study the Hamiltonian dynamics of the spherical spin model with fully-connected two-body random interactions. In the statistical physics framework, the potential energy is of the so-called p = 2 kind, closely linked to the scalar field theory. Most importantly for our setting, the energy conserving dynamics are equivalent to the ones of the Neumann integrable model. We take initial conditions from the Boltzmann equilibrium measure at a temperature that can be above or below the static phase transition, typical of a disordered (paramagnetic) or of an ordered (disguised ferromagnetic) equilibrium phase. We subsequently evolve the configurations with Newton dynamics dictated by a different Hamiltonian, obtained from an instantaneous global rescaling of the elements in the interaction random matrix. In the limit of infinitely many degrees of freedom, , we identify three dynamical phases depending on the parameters that characterise the initial state and the final Hamiltonian. We next set the analysis of the system with finite number of degrees of freedom in terms of N non-linearly coupled modes. We argue that in the limit the modes decouple at long times. We evaluate the mode temperatures and we relate them to the frequency-dependent effective temperature measured with the fluctuation-dissipation relation in the frequency domain, similarly to what was recently proposed for quantum integrable cases. Finally, we analyse the N - 1 integrals of motion, notably, their scaling with N, and we use them to show that the system is out of equilibrium in all phases, even for parameters that show an apparent Gibbs-Boltzmann behaviour of the global observables. We elaborate on the role played by these constants of motion after the quench and we briefly discuss the possible description of the asymptotic dynamics in terms of a generalised Gibbs ensemble.
Abstract The collective and purely relaxational dynamics of quantum many-body systems after a quench at temperature T = 0, from a disordered state to various phases is studied through the exact solution of the quantum Langevin equation of the spherical and the O ( n )-model in the limit n → ∞. The stationary state of the quantum dynamics is shown to be a non-equilibrium state. The quantum spherical and the quantum O ( n )-model for n → ∞ are in the same dynamical universality class. The long-time behaviour of single-time and two-time correlation and response functions is analysed and the universal exponents which characterise quantum coarsening and quantum ageing are derived. The importance of the non-Markovian long-time memory of the quantum noise is elucidated by comparing it with an effective Markovian noise having the same scaling behaviour and with the case of non-equilibrium classical dynamics.
It is believed that an isolated and far-from-equilibrium quantum many-body system should try to attain equilibrium via a mechanism whereby any given subsystem acts as an open quantum system that is coupled to an environment, which is the complementary part of the full system, and undergoes a complicated equilibration process such that all the subsystems in the long-time limit attain equilibrium states compatible with the global equilibrium state. This picture begs the question whether the dynamics of any given subsystem is Markovian (monotonic loss of information and memory) or non-Markovian. In this work, by numerically probing the dynamical behaviour of the quantum distances between temporally-separated states of small subsystems, we reveal the telltale signatures of (non-)Markovianity of the dynamics of subsystems of an isolated quantum spin system brought in the far-from-equilibrium regime, exemplified with the mixed-field Ising spin chain quenched between parameter regimes deep inside its magnetically ordered and disordered regimes. Additionally, remarkably systematic behaviour is seen in a measure of classical distances between the quantum states of the considered subsystems. These features strongly depend on the direction of quenching in the parameter space, with paramagnetic-to-ferromagnetic quenches offering considerably stronger signatures of subsystem non-Markovianity, for which we offer heuristic arguments.
Among the many non-equilibrium processes encountered in nature we deal with two different but related aspects. One is the non-equilibrium relaxation process that is at the origin of \'aging phenomena••, and the other one is a non-equilibrium phase transition, called ••dynamic phase transition••. One of the main purposes of our research is to explore more realistic situations than studied previously. Indeed, in the study of aging phenomena certain kinds of disorder effects are considered, and we introduce the ••surface•• as a spatial boundary to the system undergoing the dynamic phase transition. In order to observe these processes as clearly as possible, we study in both cases simple spin systems. Using Monte Carlo simulations we first investigate aging in three-dimensional Ising spin glasses as well as in two-dimensional Ising models with disorder quenched to low temperatures. The time-dependent dynamical correlation length L(t) is determined numerically and the scaling behavior of various two-time quantities as a function of L(t)/L(s) is discussed where t and s are two different times. For disordered Ising models deviations of L(t) from algebraic growth law show up. The generalized scaling forms as a function of L(t)/L(s) reveal a generic simple aging scenario for Ising spin glasses as well as for disordered Ising ferromagnets. We also study the local critical phenomena at a dynamic phase transition by means of numerical simulations of kinetic Ising models with surfaces subjected to a periodic oscillating field. We examine layer-dependent quantities, such as the period-averaged magnetization per layer Q(z) and the layer susceptibility ¥ö(z), and determine local critical exponents through finite size scaling. Both for two and three dimensions, we find that the values of the surface exponents differ from those of the equilibrium critical surface. It is revealed that the surface phase diagram of the non-equilibrium system is not identical to that of the equilibrium system in three dimensions.
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.