We study the dynamics of relaxation and thermalization in an exactly solvable model of a quantum particle interacting with a thermal bath of oscillators with the goal of understanding the effects of off-shell processes. The focus is to compare the exact evolution of the distribution function with different approximations to the relaxational dynamics: Boltzmann, non-Markovian and Markovian quantum kinetics. The time evolution of the distribution function is evaluated exactly using different approaches where each method provides different insights. There are two different cases that are studied in detail: (i) no stable particle states below threshold of the bath and a quasiparticle resonance above it and (ii) a stable discrete exact ``particle'' state below threshold. The exact solution for the evolution allows us to investigate the concept of the formation time of a quasiparticle and to study the difference between the relaxation of the distribution of particles and quasiparticles. For the case of quasiparticles in the continuum (resonances) the exact quasiparticle distribution asymptotically tends to a statistical equilibrium distribution that differs from a simple Bose-Einstein form as a result of off-shell processes such as the strength of the quasiparticle poles, the width of the unstable particle and proximity to thresholds. In case (ii), the distribution of particles does not thermalize with the bath. We study the kinetics of thermalization and relaxation by deriving a non-Markovian quantum kinetic equation which resums the perturbative series and includes off-shell effects. A Markovian approximation that includes off-shell contributions and the usual Boltzmann equation are obtained from the quantum kinetic equation in the limit of a wide separation of time scales upon different coarse-graining assumptions. The relaxational dynamics predicted by the non-Markovian, Markovian and Boltzmann approximations are compared to the exact result of the model. The Boltzmann approach is seen to fail in the case of wide resonances and when threshold and renormalization effects are important. Implications for thermalization in field theory models are discussed.
We study the dynamics of relaxation and thermalization in an exactly solvable model of a particle interacting with a harmonic oscillator bath. Our goal is to understand the effects of non-Markovian processes on the relaxational dynamics and to compare the exact evolution of the distribution function with approximate Markovian and non-Markovian quantum kinetics. There are two different cases that are studied in detail: (i) a quasiparticle (resonance) when the renormalized frequency of the particle is above the frequency threshold of the bath and (ii) a stable renormalized "particle" state below this threshold. The time evolution of the occupation number for the particle is evaluated exactly using different approaches that yield to complementary insights. The exact solution allows us to investigate the concept of the formation time of a quasiparticle and to study the difference between the relaxation of the distribution of bare particles and that of quasiparticles. For the case of quasiparticles, the exact occupation number asymptotically tends to a statistical equilibrium distribution that differs from a simple Bose-Einstein form as a result of off-shell processes whereas in the stable particle case, the distribution of particles does not thermalize with the bath. We derive a non-Markovian quantum kinetic equation which resums the perturbative series and includes off-shell effects. A Markovian approximation that includes off-shell contributions and the usual Boltzmann equation (energy conserving) are obtained from the quantum kinetic equation in the limit of wide separation of time scales upon different coarse-graining assumptions. The relaxational dynamics predicted by the non-Markovian, Markovian, and Boltzmann approximations are compared to the exact result. The Boltzmann approach is seen to fail in the case of wide resonances and when threshold and renormalization effects are important.
The quantum time evolution of ${\ensuremath{\varphi}}^{4}$-field theory for a spatially homogeneous system in 2+1 space-time dimensions is investigated numerically for out-of-equilibrium initial conditions on the basis of the Kadanoff-Baym equations including the tadpole and sunset self-energies. Whereas the tadpole self-energy yields a dynamical mass, the sunset self-energy is responsible for dissipation and an equilibration of the system. In particular we address the dynamics of the spectral (``off-shell'') distributions of the excited quantum modes and the different phases in the approach to equilibrium described by Kubo-Martin-Schwinger relations for thermal equilibrium states. The investigation explicitly demonstrates that the only translation invariant solutions representing the stationary fixed points of the coupled equation of motions are those of full thermal equilibrium. They agree with those extracted from the time integration of the Kadanoff-Baym equations for $\stackrel{\ensuremath{\rightarrow}}{t}\ensuremath{\infty}.$ Furthermore, a detailed comparison of the full quantum dynamics to more approximate and simple schemes such as that of a standard kinetic (on-shell) Boltzmann equation is performed. Our analysis shows that the consistent inclusion of the dynamical spectral function has a significant impact on relaxation phenomena. The different time scales that are involved in the dynamical quantum evolution towards a complete thermalized state are discussed in detail. We find that far off-shell $1\ensuremath{\leftrightarrow}3$ processes are responsible for chemical equilibration, which is missed in the Boltzmann limit. Finally, we briefly address the case of (bare) massless fields. For sufficiently large couplings $\ensuremath{\lambda}$ we observe the onset of Bose condensation, where our scheme within symmetric ${\ensuremath{\varphi}}^{4}$ theory breaks down.
A new approach to the dynamics of relaxation and kinetics of thermalization in a scalar field theory is presented that incorporates the relevant time scales through the resummation of hard thermal loops. An alternative derivation of the kinetic equations for the "quasiparticle" distribution functions is obtained that allows a clear understanding of the different "coarse-graining" approximations usually involved in a kinetic description. This method leads to a systematic perturbative expansion to obtain the kinetic equations including hard thermal loop resummation and to an improvement including renormalization, off-shell effects, and contributions that change chemical equilibrium on short time scales. As a by-product of these methods we establish the equivalence between the relaxation time scale in the linearized equation of motion of the quasiparticles and the thermalization time scale of the quasiparticle distribution function in the "relaxation time approximation" including hard thermal loop effects. Hard thermal loop resummation dramatically modifies the scattering rate for long wavelength modes as compared to the usual (semi)classical estimate. Relaxation and kinetics are studied both in the unbroken and broken symmetry phases of the theory. The broken symmetry phase also provides the setting to obtain the contribution to the kinetic equations from processes that involve decay of a heavy scalar into light scalar particles in the medium.
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