Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Thermal relaxation for particle systems in interaction with several bosonic heat reservoirs

Matthias Mück

Published Jan 1, 2005• 7 citations
Publisher

Abstract

The present thesis is concerned with the study of a quantum physical system composed of a small particle system (such as a spin chain) and several quantized massless boson fields (as photon gasses or phonon fields) at positive temperature. The setup serves as a simplified model for matter in interaction with thermal "radiation" from different sources. Hereby, questions concerning the dynamical and thermodynamic properties of particle-boson configurations far from thermal equilibrium are in the center of interest. We study a specific situation where the particle system is brought in contact with the boson systems (occasionally referred to as heat reservoirs) where the reservoirs are prepared close to thermal equilibrium states, each at a different temperature. We analyze the interacting time evolution of such an initial configuration and we show thermal relaxation of the system into a stationary state, i.e., we prove the existence of a time invariant state which is the unique limit state of the considered initial configurations evolving in time. As long as the reservoirs have been prepared at different temperatures, this stationary state features thermodynamic characteristics as stationary energy fluxes and a positive entropy production rate which distinguishes it from being a thermal equilibrium at any temperature. Therefore, we refer to it as non-equilibrium stationary state or simply NESS. The physical setup is phrased mathematically in the language of C*-algebras. The thesis gives an extended review of the application of operator algebraic theories to quantum statistical mechanics and introduces in detail the mathematical objects to describe matter in interaction with radiation. The C*-theory is adapted to the concrete setup. The algebraic description of the system is lifted into a Hilbert space framework. The appropriate Hilbert space representation is given by a bosonic Fock space over a suitable L2-space. The first part of the present work is concluded by the derivation of a spectral theory which connects the dynamical and thermodynamic features with spectral properties of a suitable generator, say K, of the time evolution in this Hilbert space setting. That way, the question about thermal relaxation becomes a spectral problem. The operator K is of Pauli-Fierz type. The spectral analysis of the generator K follows. This task is the core part of the work and it employs various kinds of functional analytic techniques. The operator K results from a perturbation of an operator L0 which describes the non-interacting particle-boson system. All spectral considerations are done in a perturbative regime, i.e., we assume that the strength of the coupling is sufficiently small. The extraction of dynamical features of the system from properties of K requires, in particular, the knowledge about the spectrum of K in the nearest vicinity of eigenvalues of the unperturbed operator L0. Since convergent Neumann series expansions only qualify to study the perturbed spectrum in the neighborhood of the unperturbed one on a scale of order of the coupling strength we need to apply a more refined tool, the Feshbach map. This technique allows the analysis of the spectrum on a smaller scale by transferring the analysis to a spectral subspace. The need of spectral information on arbitrary scales requires an iteration of the Feshbach map. This procedure leads to an operator-theoretic renormalization group. The reader is introduced to the Feshbach technique and the renormalization procedure based on it is discussed in full detail. Further, it is explained how the spectral information is extracted from the renormalization group flow. The present dissertation is an extension of two kinds of a recent research contribution by Jakšić and Pillet to a similar physical setup. Firstly, we consider the more delicate situation of bosonic heat reservoirs instead of fermionic ones, and secondly, the system can be studied uniformly for small reservoir temperatures. The adaption of the Feshbach map-based renormalization procedure by Bach, Chen, Fröhlich, and Sigal to concrete spectral problems in quantum statistical mechanics is a further novelty of this work.

DOI
10.25358/openscience-1884

Related papers

Interacting Quantum Field Theory at Positive Temperature: States and Properties

Stefano Galanda

In this thesis interacting relativistic and non relativistic Quantum Field Theories (QFT) are studied, investigating their equilibrium configurations and their long time evolution in the framework of perturbative Algebraic Quantum Field Theory, studying also convergence in a specific model. In the first part, concerning the relativistic case, we study, in formal perturbation theory, interacting KMS states generalising results by Fredenhagen and Lindner to fermionic field theories. As an application their stability is investigated together with the characterisation of the Debye screening effect in plasmas via the explicit solution of the Semiclassical Maxwell equations. Connected to these results, we also initiate a systematic study of secular growths appearing in the truncated perturbative series and affecting the validity of perturbation theory in its long time predicting capability. We restrict to theories on Minkowski spacetime and provide, for interactions depending on time via a switching process, a sufficient condition preventing secular growths in the truncated perturbative series. In the second part, we focus on a model of non relativistic quantum field theory. We construct equilibrium states at positive temperature in the presence of a condensation phase for a gas of non relativistic Bose particles on an infinite space interacting through a localised two body interaction. The convergence is proven, both for the interacting KMS state and for the entropy of the interacting KMS state relative to the free KMS state. The strategy of our proof combines ideas of Algebraic QFT, Constructive QFT and methods of statistical mechanics in the operator algebraic setting. Limits where the localisation of the two-body interaction is removed are eventually discussed, in combination with limits of the other physical parameters of the theory such as temperature, strength of the interaction and chemical potential. Finally, in the appendix, we present results on the Araki-Uhlmann relative entropy in fermionic quantum field theories, as well as an entropy-area relation, obtained via semiclassical methods, for the cosmological horizon in the static patch of de Sitter spacetime sourced by a coherent excitation of a free quantum scalar field. This completes the collection of original results obtained during the PhD.

•Dec 18, 2025•0 citations

Dynamics of Open Quantum Systems I, Oscillation and Decay

Marco Merkli

We develop a framework to analyze the dynamics of a finite-dimensional quantum system S in contact with a reservoir R . The full, interacting SR dynamics is unitary. The reservoir has a stationary state but otherwise dissipative dynamics. We identify a main part of the full dynamics, which approximates it for small values of the SR coupling constant, uniformly for all times t≥0 . The main part consists of explicit oscillating and decaying parts. We show that the reduced system evolution is Markovian for all times. The technical novelty is a detailed analysis of the link between the dynamics and the spectral properties of the generator of the SR dynamics, based on Mourre theory. We allow for SR interactions with little regularity, meaning that the decay of the reservoir correlation function only needs to be polynomial in time, improving on the previously required exponential decay.In this work we distill the structural and technical ingredients causing the characteristic features of oscillation and decay of the SR dynamics. In the companion paper [27] we apply the formalism to the concrete case of an N -level system linearly coupled to a spatially infinitely extended thermal bath of non-interacting Bosons.

arXiv:2105.00015•Quantum•Jan 3, 2022•19 citations

Heat operator approach to quantum stochastic thermodynamics in the strong-coupling regime

Sheikh Parvez Mandal, Mahasweta Pandit, Khalak Mahadeviya +2 more

Heat exchanged between an open quantum system and its environment exhibits fluctuations that carry crucial signatures of the underlying dynamics. Within the well-established two-point measurement scheme, we identify a “heat operator,” whose moments with respect to the vacuum state of a thermofield-doubled Hilbert space correspond to the stochastic moments of the heat exchanged with a bath. This recasts heat statistics as a unitary time evolution problem, which we solve by combining chain-mapped reservoirs with tensor-network propagation. In a multibath setup, all total and bath-resolved heat moments then follow from a single pure state evolution. We employ this approach to compute transient and steady-state heat fluctuations in Ohmic spin-boson models in and out of equilibrium, accessing the challenging low-temperature and long memory time regimes of the environment. In the nonequilibrium case, we show a crossover in the Fano factor from super-Poissonian to nearly Poissonian statistics under strong coupling asymmetry, corresponding to thermal rectification behavior. The method applies to noninteracting (bosonic or fermionic) nonequilibrium environments with arbitrary spectral densities, offering a powerful, nonperturbative framework for understanding heat transfer in open quantum systems.

arXiv:2504.10631•Physical Review Research•Feb 2, 2026•1 citation

Thermalization and Out-of-Equilibrium Dynamics in Open Quantum Many-Body Systems

Michael Buchhold

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 $\\mathcal{A}(\\omega)\\sim\\sqrt{\\omega}$. 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 $\\sigma^R=\\mbox{Im}\\Sigma^R$, 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 $\\sim\\tau^{-\\eta\\sub{D}}$ with $\\eta_D=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 $t\\rightarrow\\infty$ 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 $t\\rightarrow\\infty$ 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 $\\sim t^{\\alpha}$ 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 $T_{t}>T_{\\infty}$, decreasing towards its final value $T_{\\infty}$. 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 $T_{t}-T_{\\infty}\\sim t^{-\\eta\\sub{D}}$ 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 $k\\sub{F}$ 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 $\\partial_tn_q\\sim\\gamma |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 $x_t\\sim t^{4/5}$. For larger distances, the system is described by an non-equilibrium phonon distribution $n_q\\sim |q|$, which leads to a new, non-equilibrium behavior of large distance observables. For instance, the phonon decay rate scales universally as $\\gamma_q\\sim |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 $\\eta_T=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

•Jun 29, 2015•1 citation

A Quantum of Thermodynamics

Henrik Wilming

To understand in detail the relation between unitary quantum theory that describes our world at the microscopic scale and thermodynamics, which was long believed only to apply to macroscopic objects, is one of the most interesting and long-standing problems in physics. Recently, this problem has received renewed attention, in particular from the community of quantum information theory, but also from the field of statistical mechanics, inspired from stochastic thermodynamics. These results suggest that thermodynamics is also relevant for individual quantum systems, provided that they can be brought into contact with thermal baths. In this thesis, I use recently developed tools to provide new results both on fundamental questions, but also on questions which are of practical relevance for potential miniaturized thermal machines. In terms of fundamental questions, the results in this thesis contribute to understanding how the basic laws of thermodynamics and statistical mechanics can be understood directly from unitary quantum mechanics. In particular, I discuss and answer the following questions: i) How can we quantify the third law of thermodynamics using information theoretic methods? ii) How can we quantify the "thermodynamic value" of a state-transition in quantum systems? iii) How can we axiomatically characterize the non-equilibrium free energy and relative entropy? iv) How can we justify statistical ensembles from an operational perspective, without having to introduce either some probability measures or an information theoretic entropy measure? v) How can we understand the equilibration of closed quantum systems, how long does it take and how difficult is it to avoid? In the second part of the thesis I discuss in detail how experimental restrictions, which become important at the quantum scale, influence the ultimate thermodynamic bounds for thermal machines. In particular, the results in this thesis provide thermodynamic bounds on work-extraction and efficiencies of thermal machines in situations where 1.) an experimenter only has limited field strengths available, 2.) an experimenter cannot control the interactions between particles, but external fields arbitrarily well, and 3.) situations in which a small quantum system can only be strongly coupled to heat baths. These bounds are tight and I provide explicit examples illustrating the different behaviours. Finally, I come back to a classic problem in statistical physics: The emergence of spontaneous symmetry breaking. Here, I provide general and rigorous new results that show how symmetry-breaking stationary states emerge from fluctuations in order parameters in dissipative lattice models.

•Jan 1, 2018•0 citations

On the relaxation to equilibrium of a quantum oscillator interacting with a radiation field

Pierre-A. Vuillermot

In this article we investigate from the point of view of spectral theory the problem of relaxation to thermodynamical equilibrium of a quantum harmonic oscillator interacting with a radiation field. Our starting point is a system of infinitely many Pauli master equations governing the time evolution of the occupation probabilities of the available quantum states. The system we consider is derived from the evolution equation for the reduced density operator obtained after the initial interaction of the oscillator with the radiation field, the latter acting as a heat bath. We provide a complete spectral analysis of the infinitesimal generator of the equations, showing thereby that it generates an infinite-dimensional system of hyperbolic type. This implies that every global solution to the equations converge exponentially rapidly toward the corresponding Gibbs equilibrium state. We also provide a complete spectral analysis of a linear pencil naturally associated with the Pauli equations. All of our considerations revolve around the notion of compactness, more specifically around the existence of a new compact embedding result involving a space of sequences of Sobolev type into a weighted space of square-integrable summable sequences.

arXiv:2404.11329•Communications in Mathematical Sciences•Jan 1, 2025•1 citation
The present thesis is concerned with the study of a quantum physical system composed of a small particle system (such as a spin chain) and several quantized massless boson fields (as photon gasses or phonon fields) at positive temperature. The setup serves as a simplified model for matter in interaction with thermal "radiation" from different sources. Hereby, questions concerning the dynamical and thermodynamic properties of particle-boson configurations far from thermal equilibrium are in the center of interest. We study a specific situation where the particle system is brought in contact with the boson systems (occasionally referred to as heat reservoirs) where the reservoirs are prepared close to thermal equilibrium states, each at a different temperature. We analyze the interacting time evolution of such an initial configuration and we show thermal relaxation of the system into a stationary state, i.e., we prove the existence of a time invariant state which is the unique limit state of the considered initial configurations evolving in time. As long as the reservoirs have been prepared at different temperatures, this stationary state features thermodynamic characteristics as stationary energy fluxes and a positive entropy production rate which distinguishes it from being a thermal equilibrium at any temperature. Therefore, we refer to it as non-equilibrium stationary state or simply NESS. The physical setup is phrased mathematically in the language of C*-algebras. The thesis gives an extended review of the application of operator algebraic theories to quantum statistical mechanics and introduces in detail the mathematical objects to describe matter in interaction with radiation. The C*-theory is adapted to the concrete setup. The algebraic description of the system is lifted into a Hilbert space framework. The appropriate Hilbert space representation is given by a bosonic Fock space over a suitable L2-space. The first part of the present work is concluded by the derivation of a spectral theory which connects the dynamical and thermodynamic features with spectral properties of a suitable generator, say K, of the time evolution in this Hilbert space setting. That way, the question about thermal relaxation becomes a spectral problem. The operator K is of Pauli-Fierz type. The spectral analysis of the generator K follows. This task is the core part of the work and it employs various kinds of functional analytic techniques. The operator K results from a perturbation of an operator L0 which describes the non-interacting particle-boson system. All spectral considerations are done in a perturbative regime, i.e., we assume that the strength of the coupling is sufficiently small. The extraction of dynamical features of the system from properties of K requires, in particular, the knowledge about the spectrum of K in the nearest vicinity of eigenvalues of the unperturbed operator L0. Since convergent Neumann series expansions only qualify to study the perturbed spectrum in the neighborhood of the unperturbed one on a scale of order of the coupling strength we need to apply a more refined tool, the Feshbach map. This technique allows the analysis of the spectrum on a smaller scale by transferring the analysis to a spectral subspace. The need of spectral information on arbitrary scales requires an iteration of the Feshbach map. This procedure leads to an operator-theoretic renormalization group. The reader is introduced to the Feshbach technique and the renormalization procedure based on it is discussed in full detail. Further, it is explained how the spectral information is extracted from the renormalization group flow. The present dissertation is an extension of two kinds of a recent research contribution by Jakšić and Pillet to a similar physical setup. Firstly, we consider the more delicate situation of bosonic heat reservoirs instead of fermionic ones, and secondly, the system can be studied uniformly for small reservoir temperatures. The adaption of the Feshbach map-based renormalization procedure by Bach, Chen, Fröhlich, and Sigal to concrete spectral problems in quantum statistical mechanics is a further novelty of this work.
In this thesis interacting relativistic and non relativistic Quantum Field Theories (QFT) are studied, investigating their equilibrium configurations and their long time evolution in the framework of perturbative Algebraic Quantum Field Theory, studying also convergence in a specific model. In the first part, concerning the relativistic case, we study, in formal perturbation theory, interacting KMS states generalising results by Fredenhagen and Lindner to fermionic field theories. As an application their stability is investigated together with the characterisation of the Debye screening effect in plasmas via the explicit solution of the Semiclassical Maxwell equations. Connected to these results, we also initiate a systematic study of secular growths appearing in the truncated perturbative series and affecting the validity of perturbation theory in its long time predicting capability. We restrict to theories on Minkowski spacetime and provide, for interactions depending on time via a switching process, a sufficient condition preventing secular growths in the truncated perturbative series. In the second part, we focus on a model of non relativistic quantum field theory. We construct equilibrium states at positive temperature in the presence of a condensation phase for a gas of non relativistic Bose particles on an infinite space interacting through a localised two body interaction. The convergence is proven, both for the interacting KMS state and for the entropy of the interacting KMS state relative to the free KMS state. The strategy of our proof combines ideas of Algebraic QFT, Constructive QFT and methods of statistical mechanics in the operator algebraic setting. Limits where the localisation of the two-body interaction is removed are eventually discussed, in combination with limits of the other physical parameters of the theory such as temperature, strength of the interaction and chemical potential. Finally, in the appendix, we present results on the Araki-Uhlmann relative entropy in fermionic quantum field theories, as well as an entropy-area relation, obtained via semiclassical methods, for the cosmological horizon in the static patch of de Sitter spacetime sourced by a coherent excitation of a free quantum scalar field. This completes the collection of original results obtained during the PhD.
We develop a framework to analyze the dynamics of a finite-dimensional quantum system S in contact with a reservoir R . The full, interacting SR dynamics is unitary. The reservoir has a stationary state but otherwise dissipative dynamics. We identify a main part of the full dynamics, which approximates it for small values of the SR coupling constant, uniformly for all times t≥0 . The main part consists of explicit oscillating and decaying parts. We show that the reduced system evolution is Markovian for all times. The technical novelty is a detailed analysis of the link between the dynamics and the spectral properties of the generator of the SR dynamics, based on Mourre theory. We allow for SR interactions with little regularity, meaning that the decay of the reservoir correlation function only needs to be polynomial in time, improving on the previously required exponential decay.In this work we distill the structural and technical ingredients causing the characteristic features of oscillation and decay of the SR dynamics. In the companion paper [27] we apply the formalism to the concrete case of an N -level system linearly coupled to a spatially infinitely extended thermal bath of non-interacting Bosons.
Heat exchanged between an open quantum system and its environment exhibits fluctuations that carry crucial signatures of the underlying dynamics. Within the well-established two-point measurement scheme, we identify a “heat operator,” whose moments with respect to the vacuum state of a thermofield-doubled Hilbert space correspond to the stochastic moments of the heat exchanged with a bath. This recasts heat statistics as a unitary time evolution problem, which we solve by combining chain-mapped reservoirs with tensor-network propagation. In a multibath setup, all total and bath-resolved heat moments then follow from a single pure state evolution. We employ this approach to compute transient and steady-state heat fluctuations in Ohmic spin-boson models in and out of equilibrium, accessing the challenging low-temperature and long memory time regimes of the environment. In the nonequilibrium case, we show a crossover in the Fano factor from super-Poissonian to nearly Poissonian statistics under strong coupling asymmetry, corresponding to thermal rectification behavior. The method applies to noninteracting (bosonic or fermionic) nonequilibrium environments with arbitrary spectral densities, offering a powerful, nonperturbative framework for understanding heat transfer in open quantum systems.
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\\mathcal{A}(\\omega)\\sim\\sqrt{\\omega}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\\sigma^R=\\mbox{Im}\\Sigma^RsigmaR=mboxImSigmaR, determined by the self-energy at equilibrium. However, for long times tau\\tautau, 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\\sim\\tau^{-\\eta\\sub{D}}simtau−etasubD with etaD=0.58\\eta_D=0.58etaD​=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=0t=0t=0, such that for t>0t>0t>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 trightarrowinftyt\\rightarrow\\inftytrightarrowinfty 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 trightarrowinftyt\\rightarrow\\inftytrightarrowinfty is a thermal state with a finite temperature T>0T>0T>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\\sim t^{\\alpha}simtalpha in space with 0<alpha<10<\\alpha<10<alpha<1 depending on the quench. Until complete thermalization (i.e. for times t<inftyt<\\inftyt<infty), the thermal region contains more energy than the prethermal and prequench region, which is expressed in a larger temperature Tt>TinftyT_{t}>T_{\\infty}Tt​>Tinfty​, decreasing towards its final value TinftyT_{\\infty}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−etasubDT_{t}-T_{\\infty}\\sim t^{-\\eta\\sub{D}}Tt​−Tinfty​simt−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 ksubFk\\sub{F}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)U(1)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∣\\partial_tn_q\\sim\\gamma |q|partialt​nq​simgamma∣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/5x_t\\sim t^{4/5}xt​simt4/5. For larger distances, the system is described by an non-equilibrium phonon distribution nqsim∣q∣n_q\\sim |q|nq​sim∣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\\gamma_q\\sim |q|^{5/3}gammaq​sim∣q∣5/3, with a new non-equilibrium exponent eta=5/3\\eta=5/3eta=5/3, which differs from equilibrium. This new, universal behavior is guaranteed by the U(1)U(1)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\\eta=5/3eta=5/3 and for the spreading of thermalization etaT=4/5\\eta_T=4/5etaT​=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
To understand in detail the relation between unitary quantum theory that describes our world at the microscopic scale and thermodynamics, which was long believed only to apply to macroscopic objects, is one of the most interesting and long-standing problems in physics. Recently, this problem has received renewed attention, in particular from the community of quantum information theory, but also from the field of statistical mechanics, inspired from stochastic thermodynamics. These results suggest that thermodynamics is also relevant for individual quantum systems, provided that they can be brought into contact with thermal baths. In this thesis, I use recently developed tools to provide new results both on fundamental questions, but also on questions which are of practical relevance for potential miniaturized thermal machines. In terms of fundamental questions, the results in this thesis contribute to understanding how the basic laws of thermodynamics and statistical mechanics can be understood directly from unitary quantum mechanics. In particular, I discuss and answer the following questions: i) How can we quantify the third law of thermodynamics using information theoretic methods? ii) How can we quantify the "thermodynamic value" of a state-transition in quantum systems? iii) How can we axiomatically characterize the non-equilibrium free energy and relative entropy? iv) How can we justify statistical ensembles from an operational perspective, without having to introduce either some probability measures or an information theoretic entropy measure? v) How can we understand the equilibration of closed quantum systems, how long does it take and how difficult is it to avoid? In the second part of the thesis I discuss in detail how experimental restrictions, which become important at the quantum scale, influence the ultimate thermodynamic bounds for thermal machines. In particular, the results in this thesis provide thermodynamic bounds on work-extraction and efficiencies of thermal machines in situations where 1.) an experimenter only has limited field strengths available, 2.) an experimenter cannot control the interactions between particles, but external fields arbitrarily well, and 3.) situations in which a small quantum system can only be strongly coupled to heat baths. These bounds are tight and I provide explicit examples illustrating the different behaviours. Finally, I come back to a classic problem in statistical physics: The emergence of spontaneous symmetry breaking. Here, I provide general and rigorous new results that show how symmetry-breaking stationary states emerge from fluctuations in order parameters in dissipative lattice models.
In this article we investigate from the point of view of spectral theory the problem of relaxation to thermodynamical equilibrium of a quantum harmonic oscillator interacting with a radiation field. Our starting point is a system of infinitely many Pauli master equations governing the time evolution of the occupation probabilities of the available quantum states. The system we consider is derived from the evolution equation for the reduced density operator obtained after the initial interaction of the oscillator with the radiation field, the latter acting as a heat bath. We provide a complete spectral analysis of the infinitesimal generator of the equations, showing thereby that it generates an infinite-dimensional system of hyperbolic type. This implies that every global solution to the equations converge exponentially rapidly toward the corresponding Gibbs equilibrium state. We also provide a complete spectral analysis of a linear pencil naturally associated with the Pauli equations. All of our considerations revolve around the notion of compactness, more specifically around the existence of a new compact embedding result involving a space of sequences of Sobolev type into a weighted space of square-integrable summable sequences.