The central theme of this thesis is that of coherence in bosonic matter-wave fields. Analogous to optical coherence; the measure of the capacity of a photonic beam to produce fringes upon interference; matter-wave coherence is in many ways the defining characteristic of the BEC phase of identical massive bosons. Its study therefore plays a critical role in the quest to better understand the physics of ultra-cold quantum gases, which has been the goal of this research. This work approaches coherence from a number of different angles and is for this reason presented in two parts. The first of these delves into the physics of perfectly coherent matterwave fields which are adequately described by ‘mean field’ equations, while the second is devoted to origins and repercussions of the diminished coherence evident in actual atom laser experiments. The presentation of the original work of this project begins with a detailed study of the interaction between an atomic and a molecular condensate, coherently coupled at zero temperature via one of several possible mechanisms (a Feshbach resonance, for example) in three spatial dimensions. Despite the presence of repulsive scattering interaction between the atoms, the mean-field equations of motion are known to support bright solitary waves — time-independent solutions which remain localised even in the absence of a trapping potential. (These solutions are the matter-wave analogue of spatio-temporal solitons in quadratic nonlinear media, but can also represent stationary states of the Schr¨odinger-Newton equations describing bosons interacting via a gravitational field.) Ignoring all interactions save the repulsive inter-atomic scattering mentioned, a variational ansatz approach is employed to analytically approximate these stationary solutions using Gaussian density distributions of the atomic and molecular clouds. By evolving these approximate solutions numerically, a comprehensive map of the accessible parameter space is constructed, showing primarily stable localised propagation. Parameter space regions where the approximate Gaussian solutions display strong ‘breathing’ oscillations are identified, and a modified ansatz having exponentially decaying tails is found to dramatically reduce the amplitude of these oscillations. Finally, by seeding a numerical relaxation algorithm with our Gaussian solutions, we find that the only parameter space region in which the approximate solutions diffuse under propagation is also devoid of exact localised stationary solutions. The discussion then proceeds beyond the scope of mean field theory, and turns to a detailed treatment of the quantum mechanical observable corresponding to the centre of mass of systems of identical particles. Firstly, careful consideration is given to its operational definition in the presence of fluctuating particle numbers. Possible pitfalls due to naive approximations when the number fluctuations are large are also addressed. Attention is then given to the fundamental quantum noise associated with this variable, which provides a physical mechanism for diminished first order coherence in degenerate Bose gases, even at zero temperature. Characteristic lower bounds (standard quantum limits) to the resulting mean position uncertainty are identified both for systems of identical bosons and identical fermions with comparable density profiles. Interestingly, fermions are found to have an intrinsically lower mean-position uncertainty than bosons, a fact which is attributed to the action of Pauli repulsion between atoms confined to a fixed density distribution. The results of numerical simulations of the exact quantum evolution of mean position variance are then presented. These demonstrate firstly the quantum diffusion of the mean position of one dimensional matter-wave solitons and secondly the evolution of variance in the centre of mass of a Bose gas undergoing forced evaporative cooling. Finally, attention is paid to the fact that the Penrose-Onsager criterion, which essentially equates Bose condensation with perfect first-order matter-wave coherence, fails to correctly detect the presence of BEC when a priori knowledge of the exact condensate mode is imperfect or unavailable. Insight is gained by drawing analogies between this criterion and entropic measures of entanglement between particular system subspaces. This motivates employment of the concept of entanglement of formation as a means to providing a more accurate measure of condensation in the general case of mixed system states. Unfortunately, this new measure remains adversely sensitive to superpositions of condensate modes. For this reason, alternative criteria based on higher order particle number correlation functions are also proposed.
The central theme of this thesis is that of coherence in bosonic matter-wave fields. Analogous to optical coherence; the measure of the capacity of a photonic beam to produce fringes upon interference; matter-wave coherence is in many ways the defining characteristic of the BEC phase of identical massive bosons. Its study therefore plays a critical role in the quest to better understand the physics of ultra-cold quantum gases, which has been the goal of this research. This work approaches coherence from a number of different angles and is for this reason presented in two parts. The first of these delves into the physics of perfectly coherent matterwave fields which are adequately described by ‘mean field’ equations, while the second is devoted to origins and repercussions of the diminished coherence evident in actual atom laser experiments. The presentation of the original work of this project begins with a detailed study of the interaction between an atomic and a molecular condensate, coherently coupled at zero temperature via one of several possible mechanisms (a Feshbach resonance, for example) in three spatial dimensions. Despite the presence of repulsive scattering interaction between the atoms, the mean-field equations of motion are known to support bright solitary waves — time-independent solutions which remain localised even in the absence of a trapping potential. (These solutions are the matter-wave analogue of spatio-temporal solitons in quadratic nonlinear media, but can also represent stationary states of the Schr¨odinger-Newton equations describing bosons interacting via a gravitational field.) Ignoring all interactions save the repulsive inter-atomic scattering mentioned, a variational ansatz approach is employed to analytically approximate these stationary solutions using Gaussian density distributions of the atomic and molecular clouds. By evolving these approximate solutions numerically, a comprehensive map of the accessible parameter space is constructed, showing primarily stable localised propagation. Parameter space regions where the approximate Gaussian solutions display strong ‘breathing’ oscillations are identified, and a modified ansatz having exponentially decaying tails is found to dramatically reduce the amplitude of these oscillations. Finally, by seeding a numerical relaxation algorithm with our Gaussian solutions, we find that the only parameter space region in which the approximate solutions diffuse under propagation is also devoid of exact localised stationary solutions. The discussion then proceeds beyond the scope of mean field theory, and turns to a detailed treatment of the quantum mechanical observable corresponding to the centre of mass of systems of identical particles. Firstly, careful consideration is given to its operational definition in the presence of fluctuating particle numbers. Possible pitfalls due to naive approximations when the number fluctuations are large are also addressed. Attention is then given to the fundamental quantum noise associated with this variable, which provides a physical mechanism for diminished first order coherence in degenerate Bose gases, even at zero temperature. Characteristic lower bounds (standard quantum limits) to the resulting mean position uncertainty are identified both for systems of identical bosons and identical fermions with comparable density profiles. Interestingly, fermions are found to have an intrinsically lower mean-position uncertainty than bosons, a fact which is attributed to the action of Pauli repulsion between atoms confined to a fixed density distribution. The results of numerical simulations of the exact quantum evolution of mean position variance are then presented. These demonstrate firstly the quantum diffusion of the mean position of one dimensional matter-wave solitons and secondly the evolution of variance in the centre of mass of a Bose gas undergoing forced evaporative cooling. Finally, attention is paid to the fact that the Penrose-Onsager criterion, which essentially equates Bose condensation with perfect first-order matter-wave coherence, fails to correctly detect the presence of BEC when a priori knowledge of the exact condensate mode is imperfect or unavailable. Insight is gained by drawing analogies between this criterion and entropic measures of entanglement between particular system subspaces. This motivates employment of the concept of entanglement of formation as a means to providing a more accurate measure of condensation in the general case of mixed system states. Unfortunately, this new measure remains adversely sensitive to superpositions of condensate modes. For this reason, alternative criteria based on higher order particle number correlation functions are also proposed.
Bose-Einstein condensate (BEC) is a new state of matter, experimentally realized in a dilute alkali gas of weakly interacting bosons. The gas is cooled to very low temperatures, of the order of nano-Kelvin, at which point a large fraction of the atoms occupy the lowest quantum state. Under such conditions, the quantum effects become apparent on a macroscopic scale. Many fundamental ideas
of quantum mechanics have been realized and proven to be true in experiments,
by studying the Bose-Einstein condensates. Amongst them, the observation of
macroscopic matter wave interference, the creation of vortices, dark and bright
solitons etc, find special mention. The dynamics of this system in presence of an
optical lattice reveals a superfluid to Mott insulator quantum phase transition.
This has opened a new avenue for research on BEC and its applications. From
a theoretical standpoint, the dynamics of BEC can be described by means of an
effective mean-field theory. This formalism is much simpler than the many-body
Schr¨odinger equation and based on a classical nonlinear evolution equation, known
as Gross-Pitaevskii equation. The nonlinear nature of the GPE leads to several interesting phenomena like soliton propagation, vortices, four wave mixing etc.,
which are closely related to nonlinear fiber optics.
The present thesis deals with nonlinear excitations and different phases of BEC
and is organized as follows: after a brief introduction, in chapter two, we start
with BEC in a one-dimensional shallow optical lattice, experiencing both cubic
and quintic nonlinearity. A superfluid density wave, of twice the periodicity of
lattice potential, is identified. This indicates matter redistribution, where atoms
from neighbouring sites have been depleted leading to a density wave behaviour.
Interestingly, it ceases to exist when only one of these interactions is operative.
We predict the loss of superfluidity through a classical dynamical phase transition,
where modulational instability leads to the loss of phase coherence. In presence
of pure two-body interaction, the superfluid density shows a periodicity commensurate
with the lattice potential. Apart from the superfluid state, a density wave
insulating phase is also found to exist, possessing two frequency modulations commensurate
with the lattice potential. A dynamical phase transition connects the
superfluid phase with an insulating phase, which in turn, is connected with the
density wave insulating phase. The solutions are found to be marginally stable, as
per the Vakhitov-Kolokolov (VK) criterion.
Chapter three describes the dynamics of the atomic-molecular Bose-Einstein
condensate (AMBEC) by mean-field Gross-Pitaevskii equation (GPE), taking into
account the atom-molecular interaction, as well as atom-molecule inter-conversion.
The two-body atom-atom, atom-molecule and molecule-molecule interactions are
analogous to the cubic nonlinearity of the standard GPE, whereas, the quadratic
nonlinear term arising from the inter-conversion between atoms and molecules is
identical to the scenario in strongly coupled BEC in one dimension. The presence
of both weak and strong interactions adds further richness to this system. We
find the molecular BEC possesses a complex dark (grey) soliton, while the atomic condensate is a bright one. The energy and momentum of this system is computed
in order to see the effect of the both weak and strong interactions on Lieb mode
profile and compare it with the pure atomic condensate profile. Another class
of solutions exist in this system, where both atomic and molecular condensates
possess bright solitons.
In chapter four, we investigate the effect of a harmonic trap on the dynamics
of grey solitons in Bose-Einstein condensates. It is observed that, these nonlinear
excitations of the cigar-shaped Bose-Einstein condensate have strong coupling
with the trap at low momenta and hence can be effectively isolated from the
Bogoliubov sound modes, which respond weakly to harmonic confinement. This
strong coupling with the trap also makes the grey soliton amenable for control and
manipulation through trap modulation and temporal variation of the two-body interaction.
We have shown that in presence of an expulsive oscillator trap, the grey
soliton gets amplified, while having a compression in its width. Exact solutions are
also obtained in case of the two component Bose-Einstein condensates (TBECs),
where one component is a grey soliton and the other a bright one. The presence
of the trap and inter-species interaction significantly affects the grey soliton dynamics.
Similar to the single component, in the present case, we have shown that
both grey and bright soliton can be compressed simultaneously. A detailed study
of the energy and momentum of the soliton is reported, in order to find the effect
of the trap on the Lieb mode.
In chapter five, we study a BEC-defect complex in a quasi-one dimensional
scenario, where the defect atom moves in the condensate with constant velocity.
Explicit solutions of the coupled mean-field equations, describing defect-grey soliton
dynamics are obtained, demonstrating the coexistence of grey soliton and a
localized defect. Unlike the case of dark soliton, where the defect trapping center
has vanishing superfluid density, the moving grey soliton necessarily possesses a finite superfluid component at the defect location. The presence of the defect can
lower the grey soliton’s maximum velocity, as compared to the defect free case,
where the maximum velocity is the sound velocity. We have also carried out the
stability analysis of the obtained solutions using renormalized momentum formalism
and Vakhitov-Kolokolov criteria and find that the obtained solutions are stable
through both the methods. The grey soliton’s energy gets substantially modified
through its interaction with the defect, opening up the possibility of its control
through defect dynamics.
Chapter six describes a parametrically forced Bose-Einstein condensate in the
presence of a general moving optical lattice. The interaction between the atoms in
the condensate and the time dependent lattice potential leads to a novel propagating
superfluid matter wave, which can be controlled through chirp management.
This system, when placed in a trap, accelerates and undergoes rapid nonlinear
compression, controlled by the chirp. The density achieves its maximum, when
the matter wave changes direction. It is also shown that this system exhibits dynamical
superfluid insulator transition, where superfluidity breaks down and the
condensate transits to an insulating phase. The stability has been investigated
using Vakhitov-Kolokolov criterion and found that it depends on the chirp, as well
as on the nature of the atom-atom interaction. The exact expression for energy
is obtained and analyzed in detail to gain physical understanding of the chirp
management of the sinusoidal excitations and also the dynamical phase transition.
When a gas of bosonic atoms is cooled sufficiently near to absolute zero, it undergoes a startling transition from a collection of randomly-behaved individual atoms to a phase in which thousands, even millions, of atoms act in collective unison. In this new phase, known as a Bose-Einstein condensate (BEC), the atoms all occupy the exact same quantum mechanical ground state and can be collectively described by a single phase-coherent wavefunction of macroscopic size. But how exactly does this harmonious behaviour arise within a collection of incoherent atoms as their temperature decreases? Does the condensate wavefunction grow from a single central nucleus, or does it emerge simultaneously throughout the system? And how does the situation differ when the gas is cooled suddenly, away from thermodynamic equilibrium? Can non-equilibrium behaviour be inherited by the condensate that forms? In this thesis we aim to answer these and other questions by studying the formation dynamics of Bose-Einstein condensates and the development of phase coherence. The research described herein makes use of several theoretical tools, which we lay out on the workbench in Chapter 2: Semi-classical Hartree-Fock theory provides a mean-field description of a condensate and thermal cloud at equilibrium; Quantum kinetic theory is used to calculate the dynamics of condensate growth into an equilibrium state, from a non-equilibrium thermal cloud; The stochastic projected Gross-Pitaevskii equation (SPGPE) employs a classical-field description to encapsulate the non-equilibrium dynamics of the condensate and low-energy region of the atom cloud, coupled to an equilibrium thermal cloud. This thesis is composed of four principal studies, each constituting a separate chapter.In Chapter 3 we study the formation of a Bose-Einstein condensate in a cigar-shaped 3D harmonic trap, induced by the controlled addition of an attractive ``dimple'' potential along the weak axis in a recent experiment [Garrett et al., Phys. Rev. A 83, 013630 (2011)]. In this manner condensation is induced without cooling due to a localized increase in the phase space density. We perform a quantitative analysis of the thermodynamic transformation in both the sudden and adiabatic regimes for a range of dimple widths and depths. We find good agreement with equilibrium calculations based on self-consistent semi-classical Hartree-Fock theory describing the condensate and thermal cloud. We observe there is an optimal dimple depth that results in a maximum in the condensate fraction. We also study the non-equilibrium dynamics of condensate formation in the sudden turn-on regime, finding good agreement for the observed time dependence of the condensate fraction with calculations based on quantum kinetic theory. In Chapter 4 we study the onset of phase coherence during the formation of a BEC following a sudden quench to below the transition temperature. Using the SPGPE, we simulate a recent experiment performed at ANU in which a sudden RF quench is applied to a dilute gas of metastable Helium atoms confined in an elongated 3D harmonic magnetic trap. At regularly-spaced intervals following the quench, broadband RF pulses are applied to homogeneously outcouple small (<5%) pulses of atoms, which expand and fall to a microchannel plate detector. The second-order correlation function at the detector, postulated to provide a measure of phase coherence of the in-situ cloud from which they were outcoupled, is measured versus time and compared with the condensate growth curve. A delay was observed in the relaxation of the second-order correlation function, relative to the peak in the number of condensate atoms. In our simulations we incorporate the effects of periodic outcoupling both the condensate and thermal cloud, and monitor the time-dependence of first- and second-order correlation functions in both position and momentum space. We observe long-range phase coherence in the final equilibrium state, though mild persistent phase fluctuations are present. We also observe the curious phenomenon of ``superbunching'' - an enhancement of density fluctuations - at specific values of the momentum, which we relate to anomalous correlations of fluctuations. Despite obtaining a good quantitative fit to the experimentally-measured condensate growth curve in our SPGPE simulations, we do not observe any evidence of the delay in coherence observed in the experiment, which remains unexplained. In Chapter 5 we use the SPGPE to study the crossover from a true BEC with long-range phase coherence to a quasicondensate dominated by phase fluctuations, for the case of a Bose gas in an highly-elongated 3D harmonic trap. We use parameters corresponding to an experiment in which a sudden quench resulted in the formation of a non-equilibrium condensate, the dynamics of which we study in the chapter that follows. Sweeping through a range of temperatures and chemical potentials within the crossover region, we reveal the behaviour of phase fluctuations through analysis of correlation functions. We compare the ideal-gas approximation for the transition temperature with the more accurate value obtained via the Binder cumulant. By comparing the phase coherence length with the spatial extent across which density fluctuations are suppressed, we identify the BEC/quasicondensate crossover point. We compare directly with predictions of mean-field theory [Petrov et al., Phys. Rev. Lett. 87, 050404 (2001)] and find that phase fluctuations persist to lower temperatures in our SPGPE simulations. We calculate anomalous correlations of the fluctuation field and show that the resulting correction to the mean-field condensate self-interaction energy is significant within the quasicondensate regime. In Chapter 6 we study the formation of a BEC into a non-equilibrium state by simulating an experiment in which a dilute gas of Rubidium-87 atoms in a highly-elongated 3D harmonic magnetic trap is rapidly quenched to below the transition temperature. Following the quench, a condensate is observed to form with an excessively elongated shape, and subsequently undergo damped quadrupole shape oscillations [Shvarchuck et al., Phys. Rev. Lett. 89, 270404 (2002)]. We simulate this experiment by modifying the SPGPE to incorporate a hydrodynamic thermal cloud, parametrized by axially and temporally-dependent temperature and chemical potential. Comparison against simulations without this axial dependence suggests the hydrodynamic behaviour of the thermal cloud plays a crucial role in the formation of a non-equilibrium condensate. Analyses of the first-order correlation function and the condensate wavefunction reveal the spontaneous formation of quantum vortices, as well as strong phase fluctuations at equilibrium indicating quasicondensate behaviour.
Research in the field of ultra-cold quantum gases provides a testbed to examine Hamiltonians from other areas of physics such as the Bose-Hubbard Hamiltonian, as well as interesting phenomena such as superconductivity and superfluidity. Of emerging interest in the field is the analysis of non-equilibrium behaviour. This thesis focuses on the non-equilibrium dynamics of Bose gases in two different elongated systems: a continuously-pumped atom laser, and the formation of solitons as a result of a rapid quench through the Bose-Einstein condensation phase transition. We begin by introducing the theoretical tools required to model these systems, including the Stochastic Projected Gross-Pitaevskii formalism that describes the dynamics of a condensate in contact with a thermal cloud. We apply the Stochastic Projected Gross-Pitaevskii formalism to model a Raman-outcoupled atom laser continuously outcoupled from a replenished condensate reservoir at finite temperature. We investigate the first and second order coherence of the output as a function of the temperature of the condensate reservoir. For large momentum kicks given to outcoupled atoms, we show for the first time that the coherence of the atom laser beam is linked to that of the condensate reservoir. We observe a decrease in first order coherence of the condensate reservoir as its temperature is increased, and show that this decrease is caused by the dephasing of thermal fluctuations. We also investigate the effects of varying the outcoupling rate and the momentum kick on the first order coherence of the beam. We demonstrate that the first order coherence of the atom laser beam deviates from that of the condensate reservoir for small momentum kicks which we identify as due to interactions between the condensate and atom laser beam. Thermal density fluctuations in the condensate cause the atom laser beam to experience a stochastic force and the temporal coherence of the atom laser beam is reduced. We find that the condensate reservoir is significantly depleted for large outcoupling rates, resulting in a reduction in the steady state condensate fraction and in the first order coherence of the atom laser beam. In a second study, the momentum linewidth of a guided atom laser outcoupled using a radio-frequency transition is investigated at zero and finite temperature. We examine the momentum distribution of the beam by measuring the transmission of atoms over a potential barrier following the experimental scheme of the group of Professor Alain Aspect at Institut d'Optique. We observe a noisy momentum distribution for outcoupling rates in the weak-outcoupling limit at zero temperature. We show that the noise is caused by bulk excitations within the condensate reservoir as a result of a density depression caused by the non-adiabatic switch-on of outcoupling. We show that these excitations can be prevented by increasing the outcoupling rate to its final value linearly over time. In addition, we demonstrate that a thermally-broadened linewidth could potentially be measured using this experimental technique, and discuss challenges associated with the measurement scheme. Finally, we focus on the non-equilibrium dynamics of condensate growth in elongated geometries. Following a rapid quench through the Bose-Einstein condensation phase transition, we observe the spontaneous appearance of density dips which we identify as solitons. We develop an algorithm for counting these solitons in order to observe the temporal dependence of the soliton number as a function of the quench rate. To compare with the predictions of the Kibble-Zurek mechanism describing topological defect formation as a result of the dynamics of phase transitions in homogeneous systems, we calculated the scaling exponent of the soliton number with respect to the quench rate. Our results were qualitatively consistent with the Kibble-Zurek predictions, but quantitative comparisons were inconclusive. In trapped systems the phase transition proceeds inhomogeneously, with a phase front that propagates in space with a velocity. If this is lower than a critical value it has been predicted that this can inhibit defect formation. However, our results do not show conclusive evidence of this in the parameter regime that we study. We find that the number of solitons reduces as the condensate reaches equilibrium as a result of the damping included in the Stochastic Projected Gross-Pitaevskii formalism.
In this thesis we study the physics of quantum many-body systems confined to one-dimensional geometries. The work was motivated by the recent success of experimentalists in developing atom traps, capable of restricting the motion of the individual atoms to a single spatial dimension. Specifically, we look at aspects of the one-dimensional Bose gas including; excitation spectrum, correlation functions, and dynamical behaviour. In Chapter \ref{ch:excitation1D} we consider the Lieb-Liniger model of interacting bosons in one-dimension. We numerically solve the equations arising from the Bethe ansatz solution for the exact many-body wave function in a finite-size system of up to twenty particles for attractive interactions. We discuss novel features of the solutions, including deviations from the well-known string solutions due to finite size effects. We present excited state string solutions in the limit of strong interactions and discuss their physical interpretation, as well as the characteristics of the quantum phase transition that occurs as a function of interaction strength in the mean-field limit. Our results are compared to those obtained via exact diagonalization of the Hamiltonian in a truncated basis. In Chapter \ref{ch:g2} we analytically calculate the spatial nonlocal pair correlation function for an interacting uniform one dimensional Bose gas at finite temperature and propose an experimental method to measure nonlocal correlations. Our results span six different physical realms, including the weakly and strongly interacting regimes. We show explicitly that the characteristic correlation lengths are given by one of four length scales: the thermal de Broglie wavelength, the mean interparticle separation, the healing length, or the phase coherence length. In all regimes, we identify the profound role of interactions and find that under certain conditions the pair correlation may develop a global maximum at a finite interparticle separation due to the competition between repulsive interactions and thermal effects. In Chapter \ref{ch:casimirdrag} we study the drag force below the critical velocity for obstacles moving in a superfluid. The absence of drag is well established in the context of the mean-field Gross-Pitaevskii theory. We calculate the next order correction due to quantum and thermal fluctuations and find a non-zero force acting on a delta-function impurity moving through a quasi-one-dimensional Bose-Einstein condensate at all subcritical velocities and at all temperatures. The force occurs due to an imbalance in the Doppler shifts of reflected quantum fluctuations from either side of the impurity. Our calculation is based on a consistent extension of Bogoliubov theory to second order in the interaction strength, and finds new analytic solutions to the Bogoliubov-de Gennes equations for a gray soliton. In Chapter \ref{ch:solitons} we study the effect of quantum noise on the stability of a soliton. We find the soliton solutions exactly define the reflectionless potentials of the Bogoliubov-de Gennes equations. This results in complete stability of the solitons in a purely one dimensional system. We look at the modifications to the density profile of a black soliton due to quantum fluctuations.
Abstract Recent experimental progress has enabled cold atomic gases to be studied at nanokelvin temperatures, creating new states of matter where quantum degeneracy occurs—Bose–Einstein condensates and degenerate Fermi gases. Such quantum states are of macroscopic dimensions. As its title suggests, this book presents the phase space theory approach to treating the physics of degenerate quantum gases, an approach already widely used in quantum optics. However, degenerate quantum gases involve massive bosonic and fermionic atoms, not massless photons. The book begins with a review of Fock states for systems of identical atoms, where large numbers of atoms occupy the various single-particle states or modes. First, separate modes are considered, and here the quantum density operator is represented by a phase space distribution function of phase space variables which replace mode annihilation and creation operators, the dynamical equation for the density operator determines a Fokker–Planck equation for the distribution function, and measurable quantities such as quantum correlation functions are given as phase space integrals. Finally, the phase space variables are replaced by time-dependent stochastic variables satisfying Langevin stochastic equations obtained from the Fokker–Planck equation, with stochastic averages giving the measurable quantities. Second, a quantum field approach is then treated, the density operator being represented by a distribution functional of field functions which replace field annihilation and creation operators, the distribution functional satisfying a functional Fokker–Planck equation, etc. A novel feature of this book is that the phase space variables for fermions are Grassmann variables, not c-numbers. However, the book shows that Grassmann distribution functions and functionals still provide equations for obtaining both analytic and numerical solutions.
In this thesis we consider the application of phase-space methods to Bose-Einstein condensates; the work comprises of three main parts: Part I: A phase-space method for the Bose-Hubbard model; Part II: An analogue model of an expanding universe in Bose-Einstein condensates. and Part III: An analogue model of an acoustic Black Hole in Bose-Einstein condensates. In part I we present a phase-space method for the Bose-Hubbard model based on the Qfunction representation. In particular, we consider two model Hamiltonians in the mean-field approximation; the first is the standard “one site” model where quantum tunneling is approximated entirely using mean-field terms; the second “two site” model explicitly includes tunneling between two adjacent sites while treating tunneling with other neighbouring sites using the meanfield approximation. The ground state is determined by minimising the classical energy functional subject to quantum mechanical constraints, which take the form of uncertainty relations. For each model Hamiltonian we compare the ground state results from the Q-function method with the exact numerical solution. The results from the Q-function method, which are easy to compute, give a good qualitative description of the main features of the Bose-Hubbard model including the superfluid to Mott insulator. We find the quantum mechanical constraints dominate the problem and show there are some limitations of the method particularly in the weak lattice regime. Analogue models of gravity have been motivated by the possibility of investigating phenomena not readily accessible in their cosmological counterparts. In particular, the prediction of quasiparticle creation in ultra-cold Bose gases in specific configurations can be viewed as an analogue to either cosmological particle creation or the Hawking effect. In part II of this thesis we investigate the analogue of cosmological particle creation in an expanding universe by numerically simulating a Bose-Einstein condensate with a time-dependent scattering length. In particular, we simulate a 2D homogeneous condensate using the classical field method via the truncated Wigner approximation. We show that for several different expansion scenarios the calculated particle production is consistent with the underlying theory. For inflationary models we find the particle production for long wavelength modes coincides with the analytic theory within the acoustic approximation, whereas the particle production is suppressed for short wavelength (ie. free-particle like) modes. Moreover, particle production is enhanced for faster expansions, approaching the analytic result for the sudden expansion in the limit of a very fast expansion. For the case of a cyclic expansion, particle production peaks for a mode frequency that is approximately half of the driving frequency as expected for parametric resonance. In part III of this thesis we investigate an acoustic black hole in a Bose-Einstein condensate, formed by two de Laval nozzles in a ring configuration — a system we refer to as the quantum de Laval nozzle. Our model is formulated in one dimension with a sinusoidal potential. For nonzero superfluid flow, this system can exhibit stable transonic flow with both black and white hole sonic horizons. Stationary states are found by solving the time-independent Gross-Pitaevskii equation subject to a phase quantisation constraint. By solving the projected Bogoliubov-de Gennes equations for the system, we also find the discrete spectrum and quasiparticle modes. There are dynamical instabilities for certain values of winding number and potential depth, for which it is possible to construct pairs of normalisable modes. We further investigate the dynamics of the system using a classical field method based on the truncated Wigner approximation. For a low winding number and unstable configuration, we find exponential growth for the pair of unstable modes, whereas there is no growth in these modes for a stable configuration. This can be interpreted as non-degenerate parametric amplification, valid for short times. In contrast, for a large winding number, there is significant growth in modes for both stable and unstable configurations. This is indicative of higher order processes neglected in the quasiparticle picture, which is further reinforced by that fact that large winding number solutions require large nonlinearities. Finally, we consider the connection of our results with the usual semi-classical prediction of the Hawking effect. For an unstable configuration, the normalised unstable modes couple equal and opposite real frequencies, so that the growth in these modes represents the closest analogy with the Hawking effect for our quantum system.