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:cond-mat/0410653 · Journal of Physics B Atomic Molecular and Optical Physics

Properties of the stochastic Gross–Pitaevskii equation: finite temperature Ehrenfest relations and the optimal plane wave representation

Ashton S. Bradley, P. B. Blakie, Crispin W. Gardiner

Published Nov 14, 2005• 43 citations
arXivHTMLPDF

Abstract

We present Ehrenfest relations for the high temperature stochastic Gross–Pitaevskii equation description of a trapped Bose gas, including the effect of growth noise and the energy cutoff. A condition for neglecting the cutoff terms in the Ehrenfest relations is found which is more stringent than the usual validity condition of the truncated Wigner or classical field method—that all modes are highly occupied . The condition requires a small overlap of the nonlinear interaction term with the lowest energy single particle state of the noncondensate band, and gives a means to constrain dynamical artefacts arising from the energy cutoff in numerical simulations. We apply the formalism to two simple test problems: (i) simulation of the Kohn mode oscillation for a trapped Bose gas at zero temperature, and (ii) computing the equilibrium properties of a finite temperature Bose gas within the classical field method. The examples indicate ways to control the effects of the cutoff, and that there is an optimal choice of plane wave basis for a given cutoff energy. This basis gives the best reproduction of the single particle spectrum, the condensate fraction and the position and momentum densities.

DOI
10.1088/0953-4075/38/23/008

Related papers

The stochastic Gross-Pitaevskii equation: II

C. W. Gardiner, Matthew J. Davis

We provide a derivation of a more accurate version of the stochastic Gross– Pitaevskii equation, as introduced by Gardiner et al (2002 J. Phys. B: At. Mol. Opt. Phys. 35 1555). This derivation does not rely on the concept of local energy and momentum conservation and is based on a quasiclassical Wigner function representation of a ‘high temperature ’ master equation for a Bose gas, which includes only modes below an energy cut-off ER that are sufficiently highly occupied (the condensate band). The modes above this cut-off (the non-condensate band) are treated as being essentially thermalized. The interaction between these two bands, known as growth and scattering processes, provides noise and damping terms in the equation of motion for the condensate band, which we call the stochastic Gross–Pitaevskii equation. This approach is distinguished by the control of the approximations made in its derivation and by the feasibility of its numerical implementation. 1.

•Jan 1, 2003•104 citations

The stochastic Gross–Pitaevskii equation: II

C W Gardiner, Matthew J. Davis

We provide a derivation of a more accurate version of the stochastic Gross–Pitaevskii equation, as introduced by Gardiner et al ( 2002 J. Phys. B: At. Mol. Opt. Phys. 35 1555 ). This derivation does not rely on the concept of local energy and momentum conservation and is based on a quasiclassical Wigner function representation of a 'high temperature' master equation for a Bose gas, which includes only modes below an energy cut-off E R that are sufficiently highly occupied (the condensate band). The modes above this cut-off (the non-condensate band) are treated as being essentially thermalized. The interaction between these two bands, known as growth and scattering processes, provides noise and damping terms in the equation of motion for the condensate band, which we call the stochastic Gross–Pitaevskii equation. This approach is distinguished by the control of the approximations made in its derivation and by the feasibility of its numerical implementation.

arXiv:cond-mat/0308044•Journal of Physics B Atomic Molecular and Optical Physics•Nov 10, 2003•154 citations

Anomalous Dissipation in Bose Gases: Analytic Theory of Energy Damping

R. G. McDonald

The stochastic projected Gross-Pitaevskii equation (SPGPE) is a finite-temperature theory of Bose-Einstein condensate dynamics, utilizing a classical-field approach to describing the dynamics of a low-energy coherent subspace in contact with an incoherent thermal reservoir. Interactions with the reservoir are encapsulated in two distinct processes, known as number damping and energy damping. Historically the energy-damping process has received little attention, at least in part due to difficulty in finding a sufficiently efficient and accurate numerical algorithm for solving the SPGPE. Rooney et al. [Phys. Rev. E 89, 013302 (2014)] developed such an algorithm, providing new opportunities to investigate the energy-damping process. In this work we provide an analytic treatment of the energy-damping process for a range of systems using functional Ito calculus.\n\nWe perform a linear fluctuation analysis on the SPGPE for a homogeneous system. We find the dispersion relation to be the usual Bogoliubov dispersion with an additional momentum dependent damping rate. The damping rate is a combination of number- and energy-damping terms with distinct momentum dependence, revealing that each damping process is dominant for excitations of different length scales. We also find the spectra of density and phase fluctuations in momentum space, again with distinct momentum dependence associated with the two damping processes. These results demonstrate that there always exists a regime of length scales where one can expect energy-damping to be the dominant process.\n\nWe derive stochastic Ehrenfest relations (SER) from the SPGPE using projected functional change of variables. The SER take the form of stochastic differential equations for matter wave moments of the system and include drift and diffusion terms corresponding to the number- and energy-damping processes. We find that for a well-chosen cutoff the SER can often admit approximate analytic solutions. Analytic solutions of the SER for a quasi-1D Thomas-Fermi system show good agreement with simulations of the 1D SPGPE.\n\nNext, we characterise the centre of mass equilibrium properties of a 3D harmonically trapped finite temperature system in the Thomas-Fermi regime using the SER for position, momentum, and particle number. We find that the energy-damping process is generally dominant over the number-damping process when considering position and momentum fluctuations. Considering the stochastic evolution of particle number in this system, we find that the steady-state fluctuation properties depend only on the number-damping process. This suggests that the effects of each reservoir interaction can be measured individually in the same system, allowing for experimental detection of the energy-damping reservoir interaction.\n\nWe finish by considering the motion of a singly-charged quantum vortex confined to a disc trap in an energy-damped system. We find a stochastic differential equation for the dissipative motion of the vortex. It is revealed that the presence of noise can cause a stabilizing effect on the vortex, extending its lifetime. We also consider the presence of a rotating thermal cloud, a scenario that leads to the vortex becoming the energetically favourable state. This allows for analysis of steady-state stochastic properties of the vortex motion.

•Jan 1, 2019•3 citations

Dynamics of hot Bose-Einstein condensates: stochastic Ehrenfest relations for number and energy damping

Rob G. McDonald, Peter S. Barnett, Fradom Atayee +1 more

Describing partially-condensed Bose gases poses a long-standing theoretical challenge. We present exact stochastic Ehrenfest relations for the stochastic projected Gross-Pitaevskii equation, including both number and energy damping mechanisms, and all projector terms that arise from the energy cutoff separating system from reservoir. We test the theory by applying it to the center of mass fluctuations of a harmonically trapped prolate system, finding close agreement between c-field simulations and analytical results. The formalism lays the foundation to analytically explore experimentally accessible hot Bose-Einstein condensates.

arXiv:1908.05809•SciPost Physics•Feb 18, 2020•6 citations

Equilibrium, relaxation, and fluctuations in homogeneous Bose-Einstein condensates: Linearized classical field analysis

Nils A. Krause, Ashton S. Bradley

The theory of open quantum systems is central to describing the dynamics and equilibration of dilute-gas Bose-Einstein condensates (BECs). We present an analysis of the linearized stochastic projected Gross-Pitaevskii equation (SPGPE) describing finite-temperature BECs. Our treatment provides an optimal choice for the cutoff that divides the Bose gas into the low-energy coherent region forming a classical wave and the high-energy thermal cloud treated as a reservoir. Moreover, it highlights the relevance of energy damping, the number-conserving scattering between thermal and coherent atoms. We analyze the equilibrium properties and near-equilibrium relaxation of a homogeneous BEC in one, two, and three dimensions at high phase-space density and calculate the autocorrelation function and power spectrum of the density and phase fluctuations. Simulations of the full nonlinear SPGPE are in close agreement and extend our arguments beyond the linear regime. Our work suggests the need for a reexamination of decay processes in BECs studied under the neglect of energy damping.

arXiv:2505.03216•Physical Review A•Aug 22, 2025•2 citations

The truncated Wigner method for Bose-condensed gases: limits of validity and applications

Alice Sinatra, Carlos Lobo, Yvan Castin

We study the truncated Wigner method applied to a weakly interacting spinless Bose-condensed gas which is perturbed away from thermal equilibrium by a time-dependent external potential. The principle of the method is to generate an ensemble of classical fields ψ( r ) which samples the Wigner quasi-distribution function of the initial thermal equilibrium density operator of the gas, and then to evolve each classical field with the Gross–Pitaevskii equation. In the first part of the paper we improve the sampling technique over our previous work (Sinatra et al 2000 J. Mod. Opt. 47 2629–44) and we test its accuracy against the exactly solvable model of the ideal Bose gas. In the second part of the paper we investigate the conditions of validity of the truncated Wigner method. For short evolution times it is known that the time-dependent Bogoliubov approximation is valid for almost pure condensates. The requirement that the truncated Wigner method reproduces the Bogoliubov prediction leads to the constraint that the number of field modes in the Wigner simulation must be smaller than the number of particles in the gas. For longer evolution times the nonlinear dynamics of the noncondensed modes of the field plays an important role. To demonstrate this we analyse the case of a three-dimensional spatially homogeneous Bose-condensed gas and we test the ability of the truncated Wigner method to correctly reproduce the Beliaev–Landau damping of an excitation of the condensate. We have identified the mechanism which limits the validity of the truncated Wigner method: the initial ensemble of classical fields, driven by the time-dependent Gross–Pitaevskii equation, thermalizes to a classical field distribution at a temperature T class which is larger than the initial temperature T of the quantum gas. When T class significantly exceeds T a spurious damping is observed in the Wigner simulation. This leads to the second validity condition for the truncated Wigner method, T class − T ≪ T , which requires that the maximum energy max of the Bogoliubov modes in the simulation does not exceed a few k B T .

arXiv:cond-mat/0201217•Journal of Physics B Atomic Molecular and Optical Physics•Aug 21, 2002•352 citations
We present Ehrenfest relations for the high temperature stochastic Gross–Pitaevskii equation description of a trapped Bose gas, including the effect of growth noise and the energy cutoff. A condition for neglecting the cutoff terms in the Ehrenfest relations is found which is more stringent than the usual validity condition of the truncated Wigner or classical field method—that all modes are highly occupied . The condition requires a small overlap of the nonlinear interaction term with the lowest energy single particle state of the noncondensate band, and gives a means to constrain dynamical artefacts arising from the energy cutoff in numerical simulations. We apply the formalism to two simple test problems: (i) simulation of the Kohn mode oscillation for a trapped Bose gas at zero temperature, and (ii) computing the equilibrium properties of a finite temperature Bose gas within the classical field method. The examples indicate ways to control the effects of the cutoff, and that there is an optimal choice of plane wave basis for a given cutoff energy. This basis gives the best reproduction of the single particle spectrum, the condensate fraction and the position and momentum densities.
We provide a derivation of a more accurate version of the stochastic Gross– Pitaevskii equation, as introduced by Gardiner et al (2002 J. Phys. B: At. Mol. Opt. Phys. 35 1555). This derivation does not rely on the concept of local energy and momentum conservation and is based on a quasiclassical Wigner function representation of a ‘high temperature ’ master equation for a Bose gas, which includes only modes below an energy cut-off ER that are sufficiently highly occupied (the condensate band). The modes above this cut-off (the non-condensate band) are treated as being essentially thermalized. The interaction between these two bands, known as growth and scattering processes, provides noise and damping terms in the equation of motion for the condensate band, which we call the stochastic Gross–Pitaevskii equation. This approach is distinguished by the control of the approximations made in its derivation and by the feasibility of its numerical implementation. 1.
We provide a derivation of a more accurate version of the stochastic Gross–Pitaevskii equation, as introduced by Gardiner et al ( 2002 J. Phys. B: At. Mol. Opt. Phys. 35 1555 ). This derivation does not rely on the concept of local energy and momentum conservation and is based on a quasiclassical Wigner function representation of a 'high temperature' master equation for a Bose gas, which includes only modes below an energy cut-off E R that are sufficiently highly occupied (the condensate band). The modes above this cut-off (the non-condensate band) are treated as being essentially thermalized. The interaction between these two bands, known as growth and scattering processes, provides noise and damping terms in the equation of motion for the condensate band, which we call the stochastic Gross–Pitaevskii equation. This approach is distinguished by the control of the approximations made in its derivation and by the feasibility of its numerical implementation.
The stochastic projected Gross-Pitaevskii equation (SPGPE) is a finite-temperature theory of Bose-Einstein condensate dynamics, utilizing a classical-field approach to describing the dynamics of a low-energy coherent subspace in contact with an incoherent thermal reservoir. Interactions with the reservoir are encapsulated in two distinct processes, known as number damping and energy damping. Historically the energy-damping process has received little attention, at least in part due to difficulty in finding a sufficiently efficient and accurate numerical algorithm for solving the SPGPE. Rooney et al. [Phys. Rev. E 89, 013302 (2014)] developed such an algorithm, providing new opportunities to investigate the energy-damping process. In this work we provide an analytic treatment of the energy-damping process for a range of systems using functional Ito calculus.\n\nWe perform a linear fluctuation analysis on the SPGPE for a homogeneous system. We find the dispersion relation to be the usual Bogoliubov dispersion with an additional momentum dependent damping rate. The damping rate is a combination of number- and energy-damping terms with distinct momentum dependence, revealing that each damping process is dominant for excitations of different length scales. We also find the spectra of density and phase fluctuations in momentum space, again with distinct momentum dependence associated with the two damping processes. These results demonstrate that there always exists a regime of length scales where one can expect energy-damping to be the dominant process.\n\nWe derive stochastic Ehrenfest relations (SER) from the SPGPE using projected functional change of variables. The SER take the form of stochastic differential equations for matter wave moments of the system and include drift and diffusion terms corresponding to the number- and energy-damping processes. We find that for a well-chosen cutoff the SER can often admit approximate analytic solutions. Analytic solutions of the SER for a quasi-1D Thomas-Fermi system show good agreement with simulations of the 1D SPGPE.\n\nNext, we characterise the centre of mass equilibrium properties of a 3D harmonically trapped finite temperature system in the Thomas-Fermi regime using the SER for position, momentum, and particle number. We find that the energy-damping process is generally dominant over the number-damping process when considering position and momentum fluctuations. Considering the stochastic evolution of particle number in this system, we find that the steady-state fluctuation properties depend only on the number-damping process. This suggests that the effects of each reservoir interaction can be measured individually in the same system, allowing for experimental detection of the energy-damping reservoir interaction.\n\nWe finish by considering the motion of a singly-charged quantum vortex confined to a disc trap in an energy-damped system. We find a stochastic differential equation for the dissipative motion of the vortex. It is revealed that the presence of noise can cause a stabilizing effect on the vortex, extending its lifetime. We also consider the presence of a rotating thermal cloud, a scenario that leads to the vortex becoming the energetically favourable state. This allows for analysis of steady-state stochastic properties of the vortex motion.
Describing partially-condensed Bose gases poses a long-standing theoretical challenge. We present exact stochastic Ehrenfest relations for the stochastic projected Gross-Pitaevskii equation, including both number and energy damping mechanisms, and all projector terms that arise from the energy cutoff separating system from reservoir. We test the theory by applying it to the center of mass fluctuations of a harmonically trapped prolate system, finding close agreement between c-field simulations and analytical results. The formalism lays the foundation to analytically explore experimentally accessible hot Bose-Einstein condensates.
The theory of open quantum systems is central to describing the dynamics and equilibration of dilute-gas Bose-Einstein condensates (BECs). We present an analysis of the linearized stochastic projected Gross-Pitaevskii equation (SPGPE) describing finite-temperature BECs. Our treatment provides an optimal choice for the cutoff that divides the Bose gas into the low-energy coherent region forming a classical wave and the high-energy thermal cloud treated as a reservoir. Moreover, it highlights the relevance of energy damping, the number-conserving scattering between thermal and coherent atoms. We analyze the equilibrium properties and near-equilibrium relaxation of a homogeneous BEC in one, two, and three dimensions at high phase-space density and calculate the autocorrelation function and power spectrum of the density and phase fluctuations. Simulations of the full nonlinear SPGPE are in close agreement and extend our arguments beyond the linear regime. Our work suggests the need for a reexamination of decay processes in BECs studied under the neglect of energy damping.
We study the truncated Wigner method applied to a weakly interacting spinless Bose-condensed gas which is perturbed away from thermal equilibrium by a time-dependent external potential. The principle of the method is to generate an ensemble of classical fields ψ( r ) which samples the Wigner quasi-distribution function of the initial thermal equilibrium density operator of the gas, and then to evolve each classical field with the Gross–Pitaevskii equation. In the first part of the paper we improve the sampling technique over our previous work (Sinatra et al 2000 J. Mod. Opt. 47 2629–44) and we test its accuracy against the exactly solvable model of the ideal Bose gas. In the second part of the paper we investigate the conditions of validity of the truncated Wigner method. For short evolution times it is known that the time-dependent Bogoliubov approximation is valid for almost pure condensates. The requirement that the truncated Wigner method reproduces the Bogoliubov prediction leads to the constraint that the number of field modes in the Wigner simulation must be smaller than the number of particles in the gas. For longer evolution times the nonlinear dynamics of the noncondensed modes of the field plays an important role. To demonstrate this we analyse the case of a three-dimensional spatially homogeneous Bose-condensed gas and we test the ability of the truncated Wigner method to correctly reproduce the Beliaev–Landau damping of an excitation of the condensate. We have identified the mechanism which limits the validity of the truncated Wigner method: the initial ensemble of classical fields, driven by the time-dependent Gross–Pitaevskii equation, thermalizes to a classical field distribution at a temperature T class which is larger than the initial temperature T of the quantum gas. When T class significantly exceeds T a spurious damping is observed in the Wigner simulation. This leads to the second validity condition for the truncated Wigner method, T class − T ≪ T , which requires that the maximum energy max of the Bogoliubov modes in the simulation does not exceed a few k B T .