We investigate the relaxation dynamics in an isolated quantum many-body system. The stationary value of an observable after relaxation is a topic of researches in the field of quantum thermalization, since thermalization is a relaxation phenomena where this stationary value coincides with the equilibrium value. Therefore, computing the stationary value in quantum many-body systems is regarded as an important problem. We, however, prove that the stationary value in quantum many-body systems is incomputable. More precisely, we show that whether the stationary value is in the vicinity of a given value or not is an undecidable problem. Our undecidable result is still satisfied when we restrict our system to a one-dimensional shift-invariant system with nearest-neighbor interaction, our initial state to a product state of a state on a single site, and our observable to a shift-sum of a one-body observable. This result clearly shows that there is no general theorem or procedure to decide the presence or absence of thermalization in a given quantum many-body system.
Thermalization is the process through which a physical system evolves toward a state of thermal equilibrium. Determining whether or not a physical system will thermalize from an initial state has been a key question in condensed matter physics. Closely related questions are determining whether observables in these systems relax to stationary values, and what those values are. Using tools from computational complexity theory, we demonstrate that given a Hamiltonian on a finite-sized system, determining whether or not it thermalizes or relaxes to a given stationary value is computationally intractable, even for a quantum computer. In particular, we show that the problem of determining whether an observable of a finite-sized quantum system relaxes to a given value is PSPACE-complete, and so no efficient algorithm for determining the value is expected to exist. Further, we show the existence of Hamiltonians for which the problem of determining whether the system thermalizes to the Gibbs expectation value is PSPACE-complete. We also show that the related problem of determining whether the system thermalizes to the microcanonical expectation value is contained in PSPACE and is PSPACE-hard under quantum polynomial time reductions. In light of recent results demonstrating undecidability of thermalization in the thermodynamic limit, our work shows that the intractability of the problem is due to inherent difficulties in many-body physics rather than particularities of infinite systems.
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.
The lack of knowledge about the detailed many-particle motion on the microscopic scale is a key issue in any theoretical description of a macroscopic experiment. For systems at or close to thermal equilibrium, statistical mechanics provides a very successful general framework to cope with this problem. However, far from equilibrium, only very few quantitative and comparably universal results are known. Here a quantum mechanical prediction of this type is derived and verified against various experimental and numerical data from the literature. It quantitatively describes the entire temporal relaxation towards thermal equilibrium for a large class (in a mathematically precisely defined sense) of closed many-body systems, whose initial state may be arbitrarily far from equilibrium.
Thermalization, the evolution of an interacting many-body system towards a thermal Gibbs ensemble after initialization in an arbitrary non-equilibrium state, is currently a phenomenon of great interest, both in theory and experiment. As the time evolution of a quantum system is unitary, the proposed mechanism of thermalization in quantum many-body systems corresponds to the so-called eigenstate thermalization hypothesis (ETH) and the typicality of eigenstates. Although this formally solves the contradiction of thermalizing but unitary dynamics in a closed quantum many-body system, it does neither make any statement on the dynamical process of thermalization itself nor in which way the coupling of the system to an environment can hinder or modify the relaxation dynamics. In this thesis, we address both the question whether or not a quantum system driven away from equilibrium is able to relax to a thermal state, which fulfills detailed balance, and if one can identify universal behavior in the non-equilibrium relaxation dynamics. As a first realization of driven quantum systems out of equilibrium, we investigate a system of Ising spins, interacting with the quantized radiation field in an optical cavity. For multiple cavity modes, this system forms a highly entangled and frustrated state with infinite correlation times, known as a quantum spin glass. In the presence of drive and dissipation, introduced by coupling the intra-cavity radiation field to the photon vacuum outside the cavity via lossy mirrors, the quantum glass state is modified in a universal manner. For frequencies below the photon loss rate, the dissipation takes over and the system shows the universal behavior of a dissipative spin glass, with a characteristic spectral density
mathcalA(omega)simsqrtomega. On the other hand, for frequencies above the loss rate, the system retains the universal behavior of a zero temperature, quantum spin glass. Remarkably, at the glass transition, the two subsystems of spins and photons thermalize to a joint effective temperature, even in the presence of photon loss. This thermalization is a consequence of the strong spin-photon interactions, which favor detailed balance in the system and detain photons from escaping the cavity. In the thermalized system, the features of the spin glass are mirrored onto the photon degrees of freedom, leading to an emergent photon glass phase. Exploiting the inherent photon loss of the cavity, we make predictions of possible measurements on the escaping photons, which contain detailed information of the state inside the cavity and allow for a precise, non-destructive measurement of the glass state. As a further set of non-equilibrium systems, we consider one-dimensional quantum fluids driven out of equilibrium, whose universal low energy theory is formed by the so-called Luttinger Liquid description, which, due to its large degree of universality, is of intense theoretical and experimental interest. A set of recent experiments in research groups in Vienna, Innsbruck and Munich have probed the non-equilibrium time-evolution of one-dimensional quantum fluids for different experimental realizations and are pushing into a time regime, where thermalization is expected. From a theoretical point of view, one-dimensional quantum fluids are particular interesting, as Luttinger Liquids are integrable and therefore, due to an infinite number of constants of motion, do not thermalize. The leading order correction to the quadratic theory is irrelevant in the sense of the renormalization group and does therefore not modify static correlation functions, however, it breaks integrability and will therefore, even if irrelevant, induce a completely different non-equilibrium dynamics as the quadratic Luttinger theory alone. In this thesis, we derive for the first time a kinetic equation for interacting Luttinger Liquids, which describes the time evolution of the excitation densities for arbitrary initial states. The resonant character of the interaction makes a straightforward derivation of the kinetic equation, using Fermi\'s golden rule, impossible and we have to develop non-perturbative techniques in the Keldysh framework. We derive a closed expression for the time evolution of the excitation densities in terms of self-energies and vertex corrections. Close to equilibrium, the kinetic equation describes the exponential decay of excitations, with a decay rate
sigmaR=mboxImSigmaR, determined by the self-energy at equilibrium. However, for long times
tau, it also reveals the presence of dynamical slow modes, which are the consequence of exactly energy conserving dynamics and lead to an algebraic decay
simtau−etasubD with
etaD=0.58. The presence of these dynamical slow modes is not contained in the equilibrium Matsubara formalism, while they emerge naturally in the non-equilibrium formalism developed in this thesis. In order to initialize a one-dimensional quantum fluid out of equilibrium, we consider an interaction quench in a model of interacting, dispersive fermions in Chap.~\\ref{sec:Quench}. In this scenario, the fermionic interaction is suddenly changed at time
t=0, such that for
t>0 the system is not in an eigenstate and therefore undergoes a non-trivial time evolution. For the quadratic theory, the stationary state in the limit
trightarrowinfty is a non-thermal, or prethermal, state, described by a generalized Gibbs ensemble (GGE). The GGE takes into account for the conservation of all integrals of motion, formed by the eigenmodes of the Hamiltonian. On the other hand, in the presence of non-linearities, the final state for
trightarrowinfty is a thermal state with a finite temperature
T>0. . The spatio-temporal, dynamical thermalization process can be decomposed into three regimes: A prequench regime on the largest distances, which is determined by the initial state, a prethermal plateau for intermediate distances, which is determined by the metastable fixed point of the quadratic theory and a thermal region on the shortest distances. The latter spreads sub-ballistically
simtalpha in space with
0<alpha<1 depending on the quench. Until complete thermalization (i.e. for times
t<infty), the thermal region contains more energy than the prethermal and prequench region, which is expressed in a larger temperature
Tt>Tinfty, decreasing towards its final value
Tinfty. As the system has achieved local detailed balance in the thermalized region, energy transport to the non-thermal region can only be performed by the macroscopic dynamical slow modes and the decay of the temperature
Tt−Tinftysimt−etasubD again witnesses the presence of these slow modes. The very slow spreading of thermalization is consistent with recent experiments performed in Vienna, which observe a metastable, prethermal state after a quench and only observe the onset of thermalization on much larger time scales. As an immediate indication of thermalization, we determine the time evolution of the fermionic momentum distribution after a quench from non-interacting to interacting fermions. For this quench scenario, the step in the Fermi distribution at the Fermi momentum
ksubF decays to zero algebraically in the absence of a non-linearity but as a stretched exponential (the exponent being proportional to the non-linearity) in the presence of a finite non-linearity. This can serve as a proof for the presence or absence of the non-linearity even on time-scales for which thermalization can not yet be observed. Finally, we consider a bosonic quantum fluid, which is driven away from equilibrium by permanent heating. The origin of the heating is atomic spontaneous emission of laser photons, which are used to create a coherent lattice potential in optical lattice experiments. This process preserves the system\'s
U(1)-invariance, i.e. conserves the global particle number, and the corresponding long-wavelength description is a heated, interacting Luttinger Liquid, for which phonon modes are continuously populated with a momentum dependent rate
partialtnqsimgamma∣q∣. In the dynamics, we identify a quasi-thermal regime for large momenta, featuring an increasing time-dependent effective temperature. In this regime, due to fast phonon-phonon scattering, detailed balance has been achieved and is expressed by a time-local, increasing temperature. The thermal region emerges locally and spreads in space sub-ballistically according to
xtsimt4/5. For larger distances, the system is described by an non-equilibrium phonon distribution
nqsim∣q∣, which leads to a new, non-equilibrium behavior of large distance observables. For instance, the phonon decay rate scales universally as
gammaqsim∣q∣5/3, with a new non-equilibrium exponent
eta=5/3, which differs from equilibrium. This new, universal behavior is guaranteed by the
U(1) invariant dynamics of the system and is insensitive to further subleading perturbations. The non-equilibrium long-distance behavior can be determined experimentally by measuring the static and dynamic structure factor, both of which clearly indicate the exponents for phonon decay,
eta=5/3 and for the spreading of thermalization
etaT=4/5. Remarkably, even in the presence of this strong external drive, the interactions and their aim to achieve detailed balance are strong enough to establish a locally emerging and spatially spreading thermal region. The physical setups in this thesis do not only reveal interesting and new dynamical features in the out-of-equilibrium time evolution of interacting systems, but they also strongly underline the high degree of universality of thermalization for the classes of models studied here. May it be a system of coupled spins and photons, where the photons are pulled away from a thermal state by Markovian photon d
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.