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

arXiv:2404.11329 · Communications in Mathematical Sciences

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

Pierre-A. Vuillermot

Published Jan 1, 2025• 1 citation
arXivHTMLPDF

Abstract

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.

DOI
10.4310/cms.250607101019

Related papers

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

Matthias Mück

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.

•Jan 1, 2005•7 citations

Coherence and decoherence processes of a harmonic oscillator coupled with finite temperature field: exact eigenbasis solution of Kossakowski-Linblad's equation

Buang Ann Tay

The eigenvalue problem of Kossakowski-Linblad’s kinetic equation associated with the reduced density matrix of a harmonic oscillator interacting with a thermal bath in equilibrium is solved. The solution gives rise to a complete orthogonal eigenbasis endowed with Hilbert space structure that has a weighted norm. We find that the eigenfunctions at finite temperature can be obtained from the eigenfunction at zero temperature through a hyperbolic rotation on the position variables. This transformation enables the extension of the simple harmonic oscillator density matrix to that of a finite temperature. We further investigate the decay of these extended states under our dissipative kinetic equation. Furthermore, the Hilbert space structure enables the proof of a H-theorem in this system. We apply the eigenbasis expansion of an initial state to analyze decohorence as well as coherence processes. We find that coherence process occurs at a longer time scale compared to decoherence process. The time scales of both processes are estimated with the eigenbasis expansion. In the same way we analyze the evolution of the coherent state. We show that in addition to the ordinary decay time, we found another time scale which is defined by the time when the motion of the peak of the coherent state become comparative to the width of the coherent state. In contrast to the ordinary decay time this new relaxation time depends on the initial value of the momentum of the oscillator. We also find that our eigenbasis is applicable to a class of non-linear interactions, with a slight extension of the form of transport coefficients due to the non-linear interactions.

•Jan 1, 2004•0 citations

Relaxation spectrum of quantum harmonic motion in the presence of nonlinear dissipative-diffusive couplings

Marcelo A. Despósito, E. Susana Hernández

We analyze the behavior of a damped quantum harmonic oscillator immersed in a heat bath with a nonlinear coupling. For this purpose, we construct a generalized master equation (GME) for the reduced density of the oscillator. We specialize the results for a quadratic coupling in the presence of a bosonic reservoir and an Ohmic dissipation model, analyzing the evolution equation for the mean value of the oscillator coordinate. In the asymptotic case we obtain a gain-loss master equation whose transition rates are polynomials in the number of quanta. No explicit solutions other than the equilibrium one can be written, since the standard methods do not apply. The equilibrium solution is investigated and the eigenvalues of the corresponding spectral problem are numerically computed to study the effect of the nonlinearity on the evolution of the system. By means of the Wigner transformation of the GME, we extract a third-order partial derivative equation that represents the semiclassical evolution of the damped oscillator. In the classical limit this is a nonlinear Fokker-Planck equation.

Physical Review A•Dec 1, 1992•5 citations

On the relaxation to quantum-statistical equilibrium of the Wigner-Weisskopf atom in a one-dimensional radiation field. V. Exact solution for infinite systems

Russell Davidson, John J. Kozak

In pursuit of the study begun in the preceding paper of the spontaneous emission of a Wigner-Weisskopf atom in a one-dimensional radiation field, the exact solution of this problem in the thermodynamic limit of a system of infinite extent is examined. The analysis reveals the existence in general of three contributions to the probability of the atom's being excited at any moment: a constant contribution and two sorts of time-dependent ones, exponential and nonexponential. For a specific choice of coupling function, this probability is written down in closed form. The link between the constant contribution and the appearance of ghost states and ergodicity is investigated, and criteria given to determine when this contribution does not arise. A numerical study is made of the expressions obtained, and the nonexponential part of the solution compared with approximations for it which have been previously obtained by various means. It is concluded that the predictions of the weak-coupling Prigogine-Résibois master equation are inadequate except when the model is ergodic in a certain sense, and some suggestions are made concerning the relevance of these studies to nonequilibrium statistical mechanics.

Journal of Mathematical Physics•Apr 1, 1973•15 citations

Relaxation to Quantum-Statistical Equilibrium of the Wigner-Weisskopf Atom in a One-Dimensional Radiation Field. III. The Quantum-Mechanical Solution

Russell Davidson, John J. Kozak

A study is undertaken to cast light on difficulties, which arose in the first two papers of this series, pertaining to the occurrence of negative probabilities in the weak-coupling solution of the generalized Prigogine-Résibois master equation for the model of the Wigner-Weisskopf atom in a one-dimensional radiation field. The Schrödinger equation is solved exactly for the model with the initial condition for spontaneous emission, and then the weak-coupling approximations to the solution, both for an infinite and for a finite system, are derived as inverse Laplace transform integrals. An extensive analysis, theoretical and numerical, of these is undertaken, and comparison is made with the corresponding results based on the master equation. In particular, quantitative estimates of the Poincaré recurrence times for finite systems are made. It is found that considerable differences exist between the statistical-mechanical and quantum-mechanical results, but that both manifest nonanalyticity in the coupling parameter as it tends to zero. Suggestions are given for further work toward the resolution of these discrepancies and a better understanding of the weak-coupling limit.

Journal of Mathematical Physics•Jun 1, 1971•39 citations

Quantum driven dissipative parametric oscillator in a blackbody radiation field

Leonardo A. Pachón, Paul W. Brumer

We consider the general open system problem of a charged quantum oscillator confined in a harmonic trap, whose frequency can be arbitrarily modulated in time, that interacts with both an incoherent quantized (blackbody) radiation field and with an arbitrary coherent laser field. We assume that the oscillator is initially in thermodynamic equilibrium with its environment, a non-factorized initial density matrix of the system and the environment, and that at t = 0 the modulation of the frequency, the coupling to the incoherent and the coherent radiation are switched on. The subsequent dynamics, induced by the presence of the blackbody radiation, the laser field, and the frequency modulation, is studied in the framework of the influence functional approach. This approach allows incorporating, in analytic closed formulae, the non-Markovian character of the oscillator-environment interaction at any temperature as well the non-Markovian character of the blackbody radiation and its zero-point fluctuations. Expressions for the time evolution of the covariance matrix elements of the quantum fluctuations and the reduced density-operator are obtained.

arXiv:1207.3104•Journal of Mathematical Physics•Jan 1, 2014•18 citations
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.
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.
The eigenvalue problem of Kossakowski-Linblad’s kinetic equation associated with the reduced density matrix of a harmonic oscillator interacting with a thermal bath in equilibrium is solved. The solution gives rise to a complete orthogonal eigenbasis endowed with Hilbert space structure that has a weighted norm. We find that the eigenfunctions at finite temperature can be obtained from the eigenfunction at zero temperature through a hyperbolic rotation on the position variables. This transformation enables the extension of the simple harmonic oscillator density matrix to that of a finite temperature. We further investigate the decay of these extended states under our dissipative kinetic equation. Furthermore, the Hilbert space structure enables the proof of a H-theorem in this system. We apply the eigenbasis expansion of an initial state to analyze decohorence as well as coherence processes. We find that coherence process occurs at a longer time scale compared to decoherence process. The time scales of both processes are estimated with the eigenbasis expansion. In the same way we analyze the evolution of the coherent state. We show that in addition to the ordinary decay time, we found another time scale which is defined by the time when the motion of the peak of the coherent state become comparative to the width of the coherent state. In contrast to the ordinary decay time this new relaxation time depends on the initial value of the momentum of the oscillator. We also find that our eigenbasis is applicable to a class of non-linear interactions, with a slight extension of the form of transport coefficients due to the non-linear interactions.
We analyze the behavior of a damped quantum harmonic oscillator immersed in a heat bath with a nonlinear coupling. For this purpose, we construct a generalized master equation (GME) for the reduced density of the oscillator. We specialize the results for a quadratic coupling in the presence of a bosonic reservoir and an Ohmic dissipation model, analyzing the evolution equation for the mean value of the oscillator coordinate. In the asymptotic case we obtain a gain-loss master equation whose transition rates are polynomials in the number of quanta. No explicit solutions other than the equilibrium one can be written, since the standard methods do not apply. The equilibrium solution is investigated and the eigenvalues of the corresponding spectral problem are numerically computed to study the effect of the nonlinearity on the evolution of the system. By means of the Wigner transformation of the GME, we extract a third-order partial derivative equation that represents the semiclassical evolution of the damped oscillator. In the classical limit this is a nonlinear Fokker-Planck equation.
In pursuit of the study begun in the preceding paper of the spontaneous emission of a Wigner-Weisskopf atom in a one-dimensional radiation field, the exact solution of this problem in the thermodynamic limit of a system of infinite extent is examined. The analysis reveals the existence in general of three contributions to the probability of the atom's being excited at any moment: a constant contribution and two sorts of time-dependent ones, exponential and nonexponential. For a specific choice of coupling function, this probability is written down in closed form. The link between the constant contribution and the appearance of ghost states and ergodicity is investigated, and criteria given to determine when this contribution does not arise. A numerical study is made of the expressions obtained, and the nonexponential part of the solution compared with approximations for it which have been previously obtained by various means. It is concluded that the predictions of the weak-coupling Prigogine-Résibois master equation are inadequate except when the model is ergodic in a certain sense, and some suggestions are made concerning the relevance of these studies to nonequilibrium statistical mechanics.
A study is undertaken to cast light on difficulties, which arose in the first two papers of this series, pertaining to the occurrence of negative probabilities in the weak-coupling solution of the generalized Prigogine-Résibois master equation for the model of the Wigner-Weisskopf atom in a one-dimensional radiation field. The Schrödinger equation is solved exactly for the model with the initial condition for spontaneous emission, and then the weak-coupling approximations to the solution, both for an infinite and for a finite system, are derived as inverse Laplace transform integrals. An extensive analysis, theoretical and numerical, of these is undertaken, and comparison is made with the corresponding results based on the master equation. In particular, quantitative estimates of the Poincaré recurrence times for finite systems are made. It is found that considerable differences exist between the statistical-mechanical and quantum-mechanical results, but that both manifest nonanalyticity in the coupling parameter as it tends to zero. Suggestions are given for further work toward the resolution of these discrepancies and a better understanding of the weak-coupling limit.
We consider the general open system problem of a charged quantum oscillator confined in a harmonic trap, whose frequency can be arbitrarily modulated in time, that interacts with both an incoherent quantized (blackbody) radiation field and with an arbitrary coherent laser field. We assume that the oscillator is initially in thermodynamic equilibrium with its environment, a non-factorized initial density matrix of the system and the environment, and that at t = 0 the modulation of the frequency, the coupling to the incoherent and the coherent radiation are switched on. The subsequent dynamics, induced by the presence of the blackbody radiation, the laser field, and the frequency modulation, is studied in the framework of the influence functional approach. This approach allows incorporating, in analytic closed formulae, the non-Markovian character of the oscillator-environment interaction at any temperature as well the non-Markovian character of the blackbody radiation and its zero-point fluctuations. Expressions for the time evolution of the covariance matrix elements of the quantum fluctuations and the reduced density-operator are obtained.