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

Thermalization Using Quantum Field Dynamics

M. Sallé, Jan Smit A, Jeroen C. Vink

Published Oct 30, 2000• 0 citations
OpenAlex

Abstract

We describe a Hartree ensemble method to approximately solve the Heisenberg equations for the ϕ 4 model in 1 + 1 dimensions. We compute the energies and number densities of the quantum particles described by the ϕ field and find that the particles initially thermalize with a Bose-Einstein distribution for the particle density. Gradually, however, the distribution changes towards classical equipartition. Using suitable initial conditions quantum thermalization is achieved much faster than the onset of this undesirable equipartition. We also show how the numerical efficiency of our method can be significantly improved. 1. Hartree ensemble approximation Non-perturbative equilibrium properties of quantum field theories are studied with great success using numerical lattice techniques. Similarly one would like to investigate quantum field dynamics non-perturbatively. This has obvious applications to e.g. heavy ion collisions and early universe physics. Furthermore, it could be used to study equilibrium properties in cases where Monte Carlo methods do not work, such as in QCD at finite density or chiral gauge theories. Using Monte Carlo techniques to compute the Minkowski time path integral is, of course, not possible. Neither can one directly solve the Heisenberg equations for the quantum field operators. Hence one has to resort to approximations, such as classical dynamics [2] (for recent work see [3]), large n or Hartree (see e.g. [4]). Here we discuss a different type of Hartree or gaussian approximation than used previously, in which we use an ensemble of gaussian wavefunctions to compute Green functions. Ref. [1] also contains a discussion of this method and a detailed presentation will appear elsewhere. We use the lattice ϕ4 model in 1+1 dimensions as a test model, which has the following Heisenberg operator equations, ˙ˆϕ = ˆπ, ˙ ˆπ = ( ∆ − µ

Related papers

1 Finite Temperature Simulations from Quantum Field Dynamics?

M. Sallé, Jan Smit A, Jeroen C. Vink

We describe a Hartree ensemble method to approximately solve the Heisenberg equations for the ϕ 4 model in 1 + 1 dimensions. We compute the energies and number densities of the quantum particles described by the ϕ field and find that the particles initially thermalize with a Bose-Einstein distribution for the particle density. Gradually, however, the distribution changes towards classical equipartition. Using suitable initial conditions quantum thermalization is achieved much faster than the onset of this undesirable equipartition. We also show how the numerical efficiency of our method can be significantly improved. 1. Hartree ensemble approximation Non-perturbative equilibrium properties of quantum field theories are studied with great success using numerical lattice techniques. Similarly one would like to investigate quantum field dynamics non-perturbatively. This has obvious applications to e.g. heavy ion collisions and early universe physics. Furthermore, it could be used to study equilibrium properties in cases where Monte Carlo methods do not work, such as in QCD at finite density or chiral gauge theories. Using Monte Carlo techniques to compute the Minkowski time path integral is, of course, not possible. Neither can one directly solve the Heisenberg equations for the quantum field operators. Hence one has to resort to approximations, such as classical dynamics [2] (for recent work see [3]), large n or Hartree (see e.g. [4]). Here we discuss a different type of Hartree or gaussian approximation than used previously, in which we use an ensemble of gaussian wavefunctions to compute Green functions. Ref. [1] also contains a discussion of this method and a detailed presentation will appear elsewhere. We use the lattice ϕ4 model in 1+1 dimensions as a test model, which has the following Heisenberg operator equations, ˙ˆϕ = ˆπ, ˙ ˆπ = ( ∆ − µ

•Jan 1, 2000•1 citation

Thermalization in a Hartree ensemble approximation to quantum field dynamics

M. Sallé, Jan Smit A, Jeroen C. Vink

For homogeneous initial conditions, Hartree (Gaussian) dynamical approximations are known to have problems with thermalization because of insufficient scattering. We attempt to improve on this by writing an arbitrary density matrix as a superposition of Gaussian pure states and applying the Hartree approximation to each member of such an ensemble. Particles can then scatter via their back reaction on the typically inhomogeneous mean fields. Starting from initial states that are far from equilibrium we numerically compute the time evolution of particle distribution functions and observe that they indeed display approximate thermalization on intermediate time scales by approaching a Bose-Einstein form. However, for very large times the distributions drift towards classical-like equipartition.

arXiv:hep-ph/0012346•Physical review. D. Particles, fields, gravitation, and cosmology/Physical review. D. Particles and fields•Jun 25, 2001•47 citations

SCALAR FIELD DYNAMICS: CLASSICAL, QUANTUM AND IN BETWEEN ∗

M. Sallé, Jan Smit A, Jeroen C. Vink +1 more

Using a Hartree ensemble approximation, we investigate the dynamics of the \\phi^4 model in 1+1 dimensions. We find that the fields initially thermalize with a Bose-Einstein distribution for the fields. Gradually, however, the distribution changes towards classical equipartition. Using suitable initial conditions quantum thermalization is achieved much faster than the onset of this undesirable equipartition. We also show how the numerical efficiency of our method can be significantly improved.

•Jan 1, 2000•1 citation

SCALAR FIELD DYNAMICS: CLASSICAL, QUANTUM AND IN BETWEEN

M. Sallé, Jan Smit A, Jeroen C. Vink

Using a Hartree ensemble approximation, we investigate the dynamics of the \f^4 model in 1+1 dimensions. We find that the fields initially thermalize with a Bose-Einstein distribution for the fields. Gradually, however, the distribution changes towards classical equipartition. Using suitable initial conditions quantum thermalization is achieved much faster than the onset of this undesirable equipartition. We also show how the numerical efficiency of our method can be significantly improved.

arXiv:hep-ph/0009120•Strong and Electroweak Matter (SEWM)•Jul 1, 2001•0 citations

An exact stochastic field method for the interacting Bose gas at thermal equilibrium

Iacopo Carusotto, Yvan Castin

We present a new exact method to numerically compute the thermodynamical properties of an interacting Bose gas in the canonical ensemble. As in our previous paper (Carusotto I, Castin Y and Dalibard J 2001 Phys. Rev. A 63 023606), we write the density operator ρ as an average of Hartree dyadics | N :ϕ 1 ⟩⟨ N :ϕ 2 | and we find stochastic evolution equations for the wavefunctions ϕ 1,2 such that the exact imaginary-time evolution of ρ is recovered after averaging over noise. In this way, the thermal equilibrium density operator can be obtained for any temperature T . The method is then applied to study the thermodynamical properties of a homogeneous one-dimensional N -boson system: although Bose-Einstein condensation cannot occur in the thermodynamical limit, a macroscopic occupation of the lowest mode of a finite system is observed at sufficiently low temperatures. If k B T >>µ, the main effect of interactions is to suppress density fluctuations and to reduce their correlation length. Different effects such as a spatial antibunching of the atoms are predicted for the opposite k B T ⩽µ regime. Our exact stochastic calculations have been compared with existing approximate theories.

arXiv:cond-mat/0108042•Journal of Physics B Atomic Molecular and Optical Physics•Nov 23, 2001•37 citations

Exact and truncated dynamics in nonequilibrium field theory

Gert Aarts A, Gian Franco Bonini, CHRISTOF WETTERICH

Nonperturbative dynamics of quantum fields out of equilibrium is often described by the time evolution of a hierarchy of correlation functions, using approximation methods such as Hartree, large ${N}_{f},$ and $n\mathrm{PI}$-effective action techniques. These truncation schemes can be implemented equally well in a classical statistical system, where results can be tested by comparison with the complete nonlinear evolution obtained by numerical methods. For a $(1+1)$-dimensional scalar field we find that the early-time behavior is reproduced qualitatively by the Hartree dynamics. The inclusion of direct scattering improves this to the quantitative level. We show that the emergence of nonthermal temperature profiles at intermediate times can be understood in terms of the fixed points of the evolution equations in the Hartree approximation. The form of the profile depends explicitly on the initial ensemble. While the truncated evolution equations do not seem to be able to get away from the fixed point, the full nonlinear evolution shows thermalization with a (surprisingly) slow relaxation.

arXiv:hep-ph/0007357•Physical review. D. Particles, fields, gravitation, and cosmology/Physical review. D. Particles and fields•Dec 21, 2000•77 citations
We describe a Hartree ensemble method to approximately solve the Heisenberg equations for the ϕ 4 model in 1 + 1 dimensions. We compute the energies and number densities of the quantum particles described by the ϕ field and find that the particles initially thermalize with a Bose-Einstein distribution for the particle density. Gradually, however, the distribution changes towards classical equipartition. Using suitable initial conditions quantum thermalization is achieved much faster than the onset of this undesirable equipartition. We also show how the numerical efficiency of our method can be significantly improved. 1. Hartree ensemble approximation Non-perturbative equilibrium properties of quantum field theories are studied with great success using numerical lattice techniques. Similarly one would like to investigate quantum field dynamics non-perturbatively. This has obvious applications to e.g. heavy ion collisions and early universe physics. Furthermore, it could be used to study equilibrium properties in cases where Monte Carlo methods do not work, such as in QCD at finite density or chiral gauge theories. Using Monte Carlo techniques to compute the Minkowski time path integral is, of course, not possible. Neither can one directly solve the Heisenberg equations for the quantum field operators. Hence one has to resort to approximations, such as classical dynamics [2] (for recent work see [3]), large n or Hartree (see e.g. [4]). Here we discuss a different type of Hartree or gaussian approximation than used previously, in which we use an ensemble of gaussian wavefunctions to compute Green functions. Ref. [1] also contains a discussion of this method and a detailed presentation will appear elsewhere. We use the lattice ϕ4 model in 1+1 dimensions as a test model, which has the following Heisenberg operator equations, ˙ˆϕ = ˆπ, ˙ ˆπ = ( ∆ − µ
We describe a Hartree ensemble method to approximately solve the Heisenberg equations for the ϕ 4 model in 1 + 1 dimensions. We compute the energies and number densities of the quantum particles described by the ϕ field and find that the particles initially thermalize with a Bose-Einstein distribution for the particle density. Gradually, however, the distribution changes towards classical equipartition. Using suitable initial conditions quantum thermalization is achieved much faster than the onset of this undesirable equipartition. We also show how the numerical efficiency of our method can be significantly improved. 1. Hartree ensemble approximation Non-perturbative equilibrium properties of quantum field theories are studied with great success using numerical lattice techniques. Similarly one would like to investigate quantum field dynamics non-perturbatively. This has obvious applications to e.g. heavy ion collisions and early universe physics. Furthermore, it could be used to study equilibrium properties in cases where Monte Carlo methods do not work, such as in QCD at finite density or chiral gauge theories. Using Monte Carlo techniques to compute the Minkowski time path integral is, of course, not possible. Neither can one directly solve the Heisenberg equations for the quantum field operators. Hence one has to resort to approximations, such as classical dynamics [2] (for recent work see [3]), large n or Hartree (see e.g. [4]). Here we discuss a different type of Hartree or gaussian approximation than used previously, in which we use an ensemble of gaussian wavefunctions to compute Green functions. Ref. [1] also contains a discussion of this method and a detailed presentation will appear elsewhere. We use the lattice ϕ4 model in 1+1 dimensions as a test model, which has the following Heisenberg operator equations, ˙ˆϕ = ˆπ, ˙ ˆπ = ( ∆ − µ
For homogeneous initial conditions, Hartree (Gaussian) dynamical approximations are known to have problems with thermalization because of insufficient scattering. We attempt to improve on this by writing an arbitrary density matrix as a superposition of Gaussian pure states and applying the Hartree approximation to each member of such an ensemble. Particles can then scatter via their back reaction on the typically inhomogeneous mean fields. Starting from initial states that are far from equilibrium we numerically compute the time evolution of particle distribution functions and observe that they indeed display approximate thermalization on intermediate time scales by approaching a Bose-Einstein form. However, for very large times the distributions drift towards classical-like equipartition.
We present a new exact method to numerically compute the thermodynamical properties of an interacting Bose gas in the canonical ensemble. As in our previous paper (Carusotto I, Castin Y and Dalibard J 2001 Phys. Rev. A 63 023606), we write the density operator ρ as an average of Hartree dyadics | N :ϕ 1 ⟩⟨ N :ϕ 2 | and we find stochastic evolution equations for the wavefunctions ϕ 1,2 such that the exact imaginary-time evolution of ρ is recovered after averaging over noise. In this way, the thermal equilibrium density operator can be obtained for any temperature T . The method is then applied to study the thermodynamical properties of a homogeneous one-dimensional N -boson system: although Bose-Einstein condensation cannot occur in the thermodynamical limit, a macroscopic occupation of the lowest mode of a finite system is observed at sufficiently low temperatures. If k B T >>µ, the main effect of interactions is to suppress density fluctuations and to reduce their correlation length. Different effects such as a spatial antibunching of the atoms are predicted for the opposite k B T ⩽µ regime. Our exact stochastic calculations have been compared with existing approximate theories.
Nonperturbative dynamics of quantum fields out of equilibrium is often described by the time evolution of a hierarchy of correlation functions, using approximation methods such as Hartree, large Nf,{N}_{f},Nf​, and nPIn\mathrm{PI}nPI-effective action techniques. These truncation schemes can be implemented equally well in a classical statistical system, where results can be tested by comparison with the complete nonlinear evolution obtained by numerical methods. For a (1+1)(1+1)(1+1)-dimensional scalar field we find that the early-time behavior is reproduced qualitatively by the Hartree dynamics. The inclusion of direct scattering improves this to the quantitative level. We show that the emergence of nonthermal temperature profiles at intermediate times can be understood in terms of the fixed points of the evolution equations in the Hartree approximation. The form of the profile depends explicitly on the initial ensemble. While the truncated evolution equations do not seem to be able to get away from the fixed point, the full nonlinear evolution shows thermalization with a (surprisingly) slow relaxation.