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

A perturbative model for the feedback cooling of finite temperature Bose gases

Matthew J. Blacker

Published Jan 1, 2020• 0 citations
Publisher

Abstract

Ultra-cold atomic gases are an ideal platform for precision measurement devices and analogue quantum simulators, which will prove instrumental in unlocking the secrets of quantum gravity and high-temperature superconductivity. However, current experimental techniques cannot cool atomic gases to simulatenously both the low entropies and high particle numbers necessary for these applications. A promising alternative is feedback cooling: using continuous-measurement feedback control to damp out energy fluctuations and cool a gas. So far, feedback cooling has been primarily studied at zero temperature, with the only finite temperature simulation achieved via a computationally expensive numerical method restricted to bosons. This thesis develops a perturbative model for the feedback cooling of a finite temperature condensed Bose gas using Bogoliubov theory, with the aim of deriving dynamic equations for the system that are both analytically tractable, and allow for fast numerical prototyping of new feedback control schemes. Using the measurement-feedback model of Szigeti \textit{et al.} \cite{szigeti_continuous_2009,szigeti_feedback_2010}, we derive a low temperature perturbative model for feedback cooling of a Bose gas in an arbitrary trapping and control potential. Using this general model, we then derive a model for the dynamics of a Bose gas in a hard box trap being cooled with an energy damping control. We complete preliminary simulations of this model in the no-backaction conditional measurement limit, damping out density fluctuations in a gas of $90\%$ condensate fraction and cooling it to $93.5 \pm 1\%$. We find that, in this limit, the dynamics of the gas are largely independent of number but significantly depend on the inter-particle interaction. We also find an optimal energy damping control strength in this limit. However, our model is not very efficient for simulation, particularly for a large number of particles and measurement strength. As an alternative, we propose an approximation scheme in which steady-state analytic solutions could be obtained from the model. Finally, we propose two methods to develop a Bogoliubov model for the feedback cooling of fermions, which would be the first finite temperature model for the Fermi gas case.

DOI
10.25911/bfj1-rw62

Related papers

Feedback control of atomic Fermi gases

Goh, Matthew

Ultracold atomic Fermi gases are the leading platform for analogue quantum simulation, and provide a promising avenue to study the origin of high-temperature superconductivity in cuprates. However, current experimental approaches to cooling Fermi gases use evaporative cooling, which is limited by poor thermalisation properties of fermions and is non-number-conserving. This prevents the creation of useful analogue simulators of many collective phenomena. This thesis is the first theoretical investigation into the use of continuous-measurement feedback control as an alternative means of cooling an atomic Fermi gas. Since tractable simulation of Fermi gas dynamics requires simplifications to the full quantum field theory, we derive and simulate a fermionic equivalent to the Gross-Pitaevskii equation, generalising a model of feedback-controlled BECs by Haine et al. to multimode ultracold atomic Fermi gases. We demonstrate that in the absence of measurement effects, a suitable control can drive an interacting Fermi gas arbitrarily close to its ground state. However, although control schemes based upon damping spatial density fluctuations work well for single-spatial-mode BECs, we show that they perform poorly for Fermi gases with a large number of atoms due to counter-oscillation of multiple spatial modes, which must exist due to Pauli exclusion. We generalise a feedback-measurement model of BECs by Szigeti et al. to a multimode atomic Fermi gas, and perform stochastic simulations of measured, feedback-controlled fermions in the single-atom and many-atom mean-field limits. The effects of measurement backaction are an important consideration, since in a realistic experiment knowledge of the system state used for feedback must be obtained from measurement, leading to competition between measurement-induced heating and feedback cooling. We show that weaker and less precise measurements cool the system to a lower equilibrium excitation energy, but are unable to place practical lower bounds on measurement strength due to the lack of a system-filter separation. When measurement-induced heating is accounted for, we find that the equilibrium energy per particle scales superlinearly, suggesting that existing control schemes which work well for bosons would not be effective for fermions. In light of this, we propose several avenues of future investigation to overcome this limitation, leaving open the possibility of feedback control of atomic Fermi gases as a pathway to analogue quantum simulation.

•Jan 1, 2019•0 citations

Feedback cooling Bose gases to quantum degeneracy

Matthew L. Goh, Zain Mehdi, Taylor, Richard L. +5 more

Degenerate quantum gases are instrumental in advancing many-body quantum physics and underpin emerging precision sensing technologies. All state-of-the-art experiments use evaporative cooling to achieve the ultracold temperatures needed for quantum degeneracy, yet evaporative cooling is extremely lossy: more than 99.9% of the gas is discarded. Such final particle number limitations constrain imaging resolution, gas lifetime, and applications leveraging macroscopic quantum coherence. Here we show that atomic Bose gases can be cooled to quantum degeneracy using real-time feedback, an entirely new method that does not suffer the same limitations as evaporative cooling. Through novel quantum-field simulations and scaling arguments, we demonstrate that an initial low-condensate-fraction thermal Bose gas can be cooled to a high-purity Bose-Einstein condensate (BEC) by feedback control, with substantially lower atomic loss than evaporative cooling. Advantages of feedback cooling are found to be robust to imperfect detection, finite resolution of the control and measurement, time delay in the control loop, and spontaneous emission. Using feedback cooling to create degenerate sources with high coherence and low entropy enables new capabilities in precision measurement, atomtronics, and few- and many-body quantum physics.

arXiv:2206.05069•Jun 10, 2022•0 citations

Multi-mode cooling of a Bose-Einstein condensate with linear quantum feedback

Zain Mehdi, Matthew L. Goh, Matthew J. Blacker +2 more

Abstract We theoretically investigate measurement-based feedback control over the motional degrees of freedom of an oblate quasi-2D atomic Bose-Einstein condensate (BEC) subject to continuous density monitoring. We develop a linear-quadratic-Gaussian (LQG) model that describes the multi-mode dynamics of the condensate's collective excitations under continuous measurement and control. Crucially, the multi-mode cold-damping feedback control we consider uses a realistic state-estimation scheme that does not rely upon a particular model of the atomic dynamics. We present analytical results showing that collective excitations can be cooled to below single-phonon average occupation (ground-state cooling) across a broad parameter regime, and identify the conditions under which the lowest steady-state phonon occupation is asymptotically achieved. Further, we develop multi-objective optimization methods that explore the trade-off between cooling speed and the final energy of the cloud, and provide numerical simulations demonstrating the ground-state cooling of the lowest ten motional modes above the condensate ground state. Our investigation provides concrete guidance on the feedback control design and parameters needed to experimentally realize a feedback-cooled BEC.

arXiv:2506.02377•Quantum Science and Technology•Sep 23, 2025•0 citations

Out Of Equilibrium Phenomena In Ultra-Cold Gases

Stefan S. Natu

The study of out-of-equilibrium dynamics in ultra-cold gases is a new and exciting field, driven largely by the recent experimental advances in controlling and imaging cold clouds. The experimental and theoretical work thus far has been somewhat exploratory and largely numerical in nature, as the very paradigms for thinking about these systems are not well established. In this thesis I consider several different scenarios of ultra-cold bosonic and fermionic gases driven out of equilibrium and study their properties. In Chapter 1, I provide an overview of the phenomenology of ultra-cold gases, highlighting the timescales governing these systems and how the experimentalist can tune them. I discuss how cold gases can be cooled and trapped and discuss the basic physics behind optical lattices. I also discuss experimental probes of these gases, in particular the new high resolution imaging techniques developed recently at Chicago, Munich and Harvard. In Chapter 2, I discuss an early experiment (circa. 2008) which observed long lived spin dynamics in a thermal spin-1/2 Fermi gas. This experiment is an nice illustration of interesting physics resulting from the separation of timescales between spin and collisional dynamics. In my opinion, it is an excellent example of why cold gases are naturally suited to studying non-equilibrium dynamics. I simulate the experiment numerically using a collisionless Boltzmann equation and explain the observed spin dynamics both qualitatively and quantitively. In Chapter 3, I continue the discussion of spin waves in thermal gases by extending previous works on spin-1/2 gases to spin-1 Bose gases. In contrast to Chapter 2, the bulk of the work in this Chapter is analytic in nature. In particular, I find a spin wave instability in the thermal spin-1 Bose gas, which is the high temperature analog of the polar to ferromagnetic transition in a spin-1 Bose Einstein condensate. In Chapter 3, I turn my attention to bosonic systems and briefly review the the Bogoliubov mean-field theory. I calculate the momentum distribution and density-density correlation function of an interacting Bose gas within the Bogoliubov framework. Then I consider bosons in an optical lattice, and introduce the Bose Hubbard model. I calculate the mean-field phase diagram of the Bose Hubbard model and then consider fluctuations about the mean field, and derive the excitation spectrum of the lattice gas in the superfluid and insulating regimes. In Chapter 4, I ask what we learn by studying the dynamics of correlation functions following a sudden change in the interactions in a superfluid. Using the Bogoliubov theory developed in Chapter 3, I will show how the underlying excitation spectrum influences the long and short time behavior of the correlation functions. By considering a lattice dispersion, I study the analogous problem in a weak optical lattice and discuss how the lattice dispersion leads to additional features in the correlation functions. I will also discuss the timescale governing the revival of the condensate fraction in a quantum depleted gas. In Chapter 5, I derive equations of motion governing the dynamics of one and two body correlation functions in the single-band Bose Hubbard model, applicable to bosons in deep lattices. I then consider a simple quench from a Mott insulating initial state to a weakly interacting final state and produce analytic expressions describing the dynamics of correlations following such a quench. I discuss the timescale for the development of long range order following such a quench. I study the problem of chapter 4 using an equations of motion approach. This approach complements the Bogoliubov approach of Chapter 4. First, I derive exact expressions for a quench to a non-interacting state. I then consider how interactions redistribute quasi-momentum to first order in perturbation theory in different dimensions. In Chapter 6, I calculate the relevant timescales for local and global dynamics in trapped lattice Bose gases, a work done in collaboration with Dr. Kaden R. A Hazzard. Using a time-dependent Gutzwiller mean-field theory, I show that the timescale for local equilibration in these systems is fast in experimental terms. I then show that due to the spatial inhomogeneities inherent to cold gases, achieving global equilibrium can be quite complicated, sometimes taking longer than the lifetime of the experiment, an issue of practical importance to current day experiments. I continue this discussion in Chapter 7 which is a collaborative work with experimentalists David McKay and Prof. Brian DeMarco from the University of Toronto and the University of Illinois, Urbana Champaign. Using experimental and numerical methods, we show that the rapid timescales for local dynamics in interacting systems invalidates a frequently used cold atom technique for mapping out the momentum distribution of atoms in an optical lattice.

•Jan 28, 2013•0 citations

Feedback-cooled Bose-Einstein condensation: Near and far from equilibrium

Evan P. Yamaguchi, Hilary M. Hurst, Ian B. Spielman

Continuously measured interacting quantum systems almost invariably heat, causing loss of quantum coherence. Here, we study Bose-Einstein condensates (BECs) subject to repeated weak measurement of the atomic density and describe several protocols for generating a feedback signal designed to remove excitations created by measurement back-action. We use a stochastic Gross-Pitaevskii equation to model the system dynamics and find that a feedback protocol utilizing momentum-dependent gain and filtering can effectively cool both one-dimensional (1D) and two-dimensional systems. The performance of these protocols is quantified in terms of the steady-state energy, entropy, and condensed fraction. These are the first feedback cooling protocols demonstrated in two dimensions, and in one dimension our optimal protocol reduces the equilibrium energy by more than a factor of 100 as compared with a previous cooling protocol developed using the same methodology. We also use this protocol to quench-cool 1D BECs from noncondensed highly excited states and find that they rapidly condense into a far-from-equilibrium state with energy orders of magnitude higher than the equilibrium ground-state energy for that condensate fraction. We explain this in terms of the near integrability of our 1D system, whereby efficiently cooled low-momentum modes are effectively decoupled from the energetic ``reservoir'' of the higher-momentum modes. We observe that the quench-cooled condensed states can have nonzero integer winding numbers described by quantized supercurrents.

arXiv:2206.04156•Physical Review A•Jun 13, 2023•8 citations

Bose gases in and out of equilibrium within the stochastic gross-pitaevskii equation

Stuart P. Cockburn

In this thesis, we apply the stochastic Gross-Pitaevskii equation, using the formulation of Stoof, to the study of weakly interacting, finite temperature atomic Bose gases, both in and out of equilibrium. Explicitly maintaining the effects of phase and density fluctuations, this state-of-the-art method is well suited to the study of low dimensional Bose systems. As such, we consider the non-equilibrium dynamics of dark matter wave solitons, oscillating within harmonically trapped, one-dimensional phase-fluctuating condensates, through numerical simulations of the stochastic Gross-Pitaevskii equation. In treating the damping due to higher lying thermal modes of the system, we include the HartreeFock contribution of the equilibrium atomic density, which is a novel feature of our work. A stochastic approach is found to be essential in capturing the shot-to-shot variation in soliton trajectories observed experimentally. Neglecting stochastic effects, we find the well-known dissipative Gross-Pitaevskii equation to recover the average dynamics of the stochastic model, on retaining the microscopically justified form for the damping of this approach. Such a coupling is found to be necessary in modelling the dissipative nature of finite temperature soliton dynamics. Furthermore, we undertake a comparison between the finite temperature equilibrium state of the stochastic Gross-Pitaevskii equation, and that predicted within a number conserving Bogoliubov approach. Through a series of numerical experiments, we compare the equilibrium properties of a harmonically trapped Bose gas, considering our results in light of those of independent theoretical treatments, and additionally assessing the sensitivity of the numerical implementation of the former method to the numerical grid spacing. In the regime in which thermal effects dominate, the stochastic Gross-Pitaevskii approach proves to be a more consistent means of generating such equilibrium states.

•Jan 1, 2010•2 citations
Ultra-cold atomic gases are an ideal platform for precision measurement devices and analogue quantum simulators, which will prove instrumental in unlocking the secrets of quantum gravity and high-temperature superconductivity. However, current experimental techniques cannot cool atomic gases to simulatenously both the low entropies and high particle numbers necessary for these applications. A promising alternative is feedback cooling: using continuous-measurement feedback control to damp out energy fluctuations and cool a gas. So far, feedback cooling has been primarily studied at zero temperature, with the only finite temperature simulation achieved via a computationally expensive numerical method restricted to bosons. This thesis develops a perturbative model for the feedback cooling of a finite temperature condensed Bose gas using Bogoliubov theory, with the aim of deriving dynamic equations for the system that are both analytically tractable, and allow for fast numerical prototyping of new feedback control schemes. Using the measurement-feedback model of Szigeti \textit{et al.} \cite{szigeti_continuous_2009,szigeti_feedback_2010}, we derive a low temperature perturbative model for feedback cooling of a Bose gas in an arbitrary trapping and control potential. Using this general model, we then derive a model for the dynamics of a Bose gas in a hard box trap being cooled with an energy damping control. We complete preliminary simulations of this model in the no-backaction conditional measurement limit, damping out density fluctuations in a gas of 90%90\%90% condensate fraction and cooling it to 93.5±1%93.5 \pm 1\%93.5±1%. We find that, in this limit, the dynamics of the gas are largely independent of number but significantly depend on the inter-particle interaction. We also find an optimal energy damping control strength in this limit. However, our model is not very efficient for simulation, particularly for a large number of particles and measurement strength. As an alternative, we propose an approximation scheme in which steady-state analytic solutions could be obtained from the model. Finally, we propose two methods to develop a Bogoliubov model for the feedback cooling of fermions, which would be the first finite temperature model for the Fermi gas case.
Ultracold atomic Fermi gases are the leading platform for analogue quantum simulation, and provide a promising avenue to study the origin of high-temperature superconductivity in cuprates. However, current experimental approaches to cooling Fermi gases use evaporative cooling, which is limited by poor thermalisation properties of fermions and is non-number-conserving. This prevents the creation of useful analogue simulators of many collective phenomena. This thesis is the first theoretical investigation into the use of continuous-measurement feedback control as an alternative means of cooling an atomic Fermi gas. Since tractable simulation of Fermi gas dynamics requires simplifications to the full quantum field theory, we derive and simulate a fermionic equivalent to the Gross-Pitaevskii equation, generalising a model of feedback-controlled BECs by Haine et al. to multimode ultracold atomic Fermi gases. We demonstrate that in the absence of measurement effects, a suitable control can drive an interacting Fermi gas arbitrarily close to its ground state. However, although control schemes based upon damping spatial density fluctuations work well for single-spatial-mode BECs, we show that they perform poorly for Fermi gases with a large number of atoms due to counter-oscillation of multiple spatial modes, which must exist due to Pauli exclusion. We generalise a feedback-measurement model of BECs by Szigeti et al. to a multimode atomic Fermi gas, and perform stochastic simulations of measured, feedback-controlled fermions in the single-atom and many-atom mean-field limits. The effects of measurement backaction are an important consideration, since in a realistic experiment knowledge of the system state used for feedback must be obtained from measurement, leading to competition between measurement-induced heating and feedback cooling. We show that weaker and less precise measurements cool the system to a lower equilibrium excitation energy, but are unable to place practical lower bounds on measurement strength due to the lack of a system-filter separation. When measurement-induced heating is accounted for, we find that the equilibrium energy per particle scales superlinearly, suggesting that existing control schemes which work well for bosons would not be effective for fermions. In light of this, we propose several avenues of future investigation to overcome this limitation, leaving open the possibility of feedback control of atomic Fermi gases as a pathway to analogue quantum simulation.
Degenerate quantum gases are instrumental in advancing many-body quantum physics and underpin emerging precision sensing technologies. All state-of-the-art experiments use evaporative cooling to achieve the ultracold temperatures needed for quantum degeneracy, yet evaporative cooling is extremely lossy: more than 99.9% of the gas is discarded. Such final particle number limitations constrain imaging resolution, gas lifetime, and applications leveraging macroscopic quantum coherence. Here we show that atomic Bose gases can be cooled to quantum degeneracy using real-time feedback, an entirely new method that does not suffer the same limitations as evaporative cooling. Through novel quantum-field simulations and scaling arguments, we demonstrate that an initial low-condensate-fraction thermal Bose gas can be cooled to a high-purity Bose-Einstein condensate (BEC) by feedback control, with substantially lower atomic loss than evaporative cooling. Advantages of feedback cooling are found to be robust to imperfect detection, finite resolution of the control and measurement, time delay in the control loop, and spontaneous emission. Using feedback cooling to create degenerate sources with high coherence and low entropy enables new capabilities in precision measurement, atomtronics, and few- and many-body quantum physics.
Abstract We theoretically investigate measurement-based feedback control over the motional degrees of freedom of an oblate quasi-2D atomic Bose-Einstein condensate (BEC) subject to continuous density monitoring. We develop a linear-quadratic-Gaussian (LQG) model that describes the multi-mode dynamics of the condensate's collective excitations under continuous measurement and control. Crucially, the multi-mode cold-damping feedback control we consider uses a realistic state-estimation scheme that does not rely upon a particular model of the atomic dynamics. We present analytical results showing that collective excitations can be cooled to below single-phonon average occupation (ground-state cooling) across a broad parameter regime, and identify the conditions under which the lowest steady-state phonon occupation is asymptotically achieved. Further, we develop multi-objective optimization methods that explore the trade-off between cooling speed and the final energy of the cloud, and provide numerical simulations demonstrating the ground-state cooling of the lowest ten motional modes above the condensate ground state. Our investigation provides concrete guidance on the feedback control design and parameters needed to experimentally realize a feedback-cooled BEC.
The study of out-of-equilibrium dynamics in ultra-cold gases is a new and exciting field, driven largely by the recent experimental advances in controlling and imaging cold clouds. The experimental and theoretical work thus far has been somewhat exploratory and largely numerical in nature, as the very paradigms for thinking about these systems are not well established. In this thesis I consider several different scenarios of ultra-cold bosonic and fermionic gases driven out of equilibrium and study their properties. In Chapter 1, I provide an overview of the phenomenology of ultra-cold gases, highlighting the timescales governing these systems and how the experimentalist can tune them. I discuss how cold gases can be cooled and trapped and discuss the basic physics behind optical lattices. I also discuss experimental probes of these gases, in particular the new high resolution imaging techniques developed recently at Chicago, Munich and Harvard. In Chapter 2, I discuss an early experiment (circa. 2008) which observed long lived spin dynamics in a thermal spin-1/2 Fermi gas. This experiment is an nice illustration of interesting physics resulting from the separation of timescales between spin and collisional dynamics. In my opinion, it is an excellent example of why cold gases are naturally suited to studying non-equilibrium dynamics. I simulate the experiment numerically using a collisionless Boltzmann equation and explain the observed spin dynamics both qualitatively and quantitively. In Chapter 3, I continue the discussion of spin waves in thermal gases by extending previous works on spin-1/2 gases to spin-1 Bose gases. In contrast to Chapter 2, the bulk of the work in this Chapter is analytic in nature. In particular, I find a spin wave instability in the thermal spin-1 Bose gas, which is the high temperature analog of the polar to ferromagnetic transition in a spin-1 Bose Einstein condensate. In Chapter 3, I turn my attention to bosonic systems and briefly review the the Bogoliubov mean-field theory. I calculate the momentum distribution and density-density correlation function of an interacting Bose gas within the Bogoliubov framework. Then I consider bosons in an optical lattice, and introduce the Bose Hubbard model. I calculate the mean-field phase diagram of the Bose Hubbard model and then consider fluctuations about the mean field, and derive the excitation spectrum of the lattice gas in the superfluid and insulating regimes. In Chapter 4, I ask what we learn by studying the dynamics of correlation functions following a sudden change in the interactions in a superfluid. Using the Bogoliubov theory developed in Chapter 3, I will show how the underlying excitation spectrum influences the long and short time behavior of the correlation functions. By considering a lattice dispersion, I study the analogous problem in a weak optical lattice and discuss how the lattice dispersion leads to additional features in the correlation functions. I will also discuss the timescale governing the revival of the condensate fraction in a quantum depleted gas. In Chapter 5, I derive equations of motion governing the dynamics of one and two body correlation functions in the single-band Bose Hubbard model, applicable to bosons in deep lattices. I then consider a simple quench from a Mott insulating initial state to a weakly interacting final state and produce analytic expressions describing the dynamics of correlations following such a quench. I discuss the timescale for the development of long range order following such a quench. I study the problem of chapter 4 using an equations of motion approach. This approach complements the Bogoliubov approach of Chapter 4. First, I derive exact expressions for a quench to a non-interacting state. I then consider how interactions redistribute quasi-momentum to first order in perturbation theory in different dimensions. In Chapter 6, I calculate the relevant timescales for local and global dynamics in trapped lattice Bose gases, a work done in collaboration with Dr. Kaden R. A Hazzard. Using a time-dependent Gutzwiller mean-field theory, I show that the timescale for local equilibration in these systems is fast in experimental terms. I then show that due to the spatial inhomogeneities inherent to cold gases, achieving global equilibrium can be quite complicated, sometimes taking longer than the lifetime of the experiment, an issue of practical importance to current day experiments. I continue this discussion in Chapter 7 which is a collaborative work with experimentalists David McKay and Prof. Brian DeMarco from the University of Toronto and the University of Illinois, Urbana Champaign. Using experimental and numerical methods, we show that the rapid timescales for local dynamics in interacting systems invalidates a frequently used cold atom technique for mapping out the momentum distribution of atoms in an optical lattice.
Continuously measured interacting quantum systems almost invariably heat, causing loss of quantum coherence. Here, we study Bose-Einstein condensates (BECs) subject to repeated weak measurement of the atomic density and describe several protocols for generating a feedback signal designed to remove excitations created by measurement back-action. We use a stochastic Gross-Pitaevskii equation to model the system dynamics and find that a feedback protocol utilizing momentum-dependent gain and filtering can effectively cool both one-dimensional (1D) and two-dimensional systems. The performance of these protocols is quantified in terms of the steady-state energy, entropy, and condensed fraction. These are the first feedback cooling protocols demonstrated in two dimensions, and in one dimension our optimal protocol reduces the equilibrium energy by more than a factor of 100 as compared with a previous cooling protocol developed using the same methodology. We also use this protocol to quench-cool 1D BECs from noncondensed highly excited states and find that they rapidly condense into a far-from-equilibrium state with energy orders of magnitude higher than the equilibrium ground-state energy for that condensate fraction. We explain this in terms of the near integrability of our 1D system, whereby efficiently cooled low-momentum modes are effectively decoupled from the energetic ``reservoir'' of the higher-momentum modes. We observe that the quench-cooled condensed states can have nonzero integer winding numbers described by quantized supercurrents.
In this thesis, we apply the stochastic Gross-Pitaevskii equation, using the formulation of Stoof, to the study of weakly interacting, finite temperature atomic Bose gases, both in and out of equilibrium. Explicitly maintaining the effects of phase and density fluctuations, this state-of-the-art method is well suited to the study of low dimensional Bose systems. As such, we consider the non-equilibrium dynamics of dark matter wave solitons, oscillating within harmonically trapped, one-dimensional phase-fluctuating condensates, through numerical simulations of the stochastic Gross-Pitaevskii equation. In treating the damping due to higher lying thermal modes of the system, we include the HartreeFock contribution of the equilibrium atomic density, which is a novel feature of our work. A stochastic approach is found to be essential in capturing the shot-to-shot variation in soliton trajectories observed experimentally. Neglecting stochastic effects, we find the well-known dissipative Gross-Pitaevskii equation to recover the average dynamics of the stochastic model, on retaining the microscopically justified form for the damping of this approach. Such a coupling is found to be necessary in modelling the dissipative nature of finite temperature soliton dynamics. Furthermore, we undertake a comparison between the finite temperature equilibrium state of the stochastic Gross-Pitaevskii equation, and that predicted within a number conserving Bogoliubov approach. Through a series of numerical experiments, we compare the equilibrium properties of a harmonically trapped Bose gas, considering our results in light of those of independent theoretical treatments, and additionally assessing the sensitivity of the numerical implementation of the former method to the numerical grid spacing. In the regime in which thermal effects dominate, the stochastic Gross-Pitaevskii approach proves to be a more consistent means of generating such equilibrium states.