After the first experimental realization of a Bose-Einstein condensate (BEC) in 1995, BECs have become a subject of intense experimental and theoretical study. In this dissertation, I present our results on the classical and quantum dynamics of BECs at zero temperature under different scenarios. First, I consider the analog of slow light in the collision of two BECs near a Feshbach resonance. The scattering length then becomes a function of the collision energy. I derive a generalization of the Gross-Pitaevskii equation for incorporating this energy dependence. In certain parameter regimes, the group velocity of a BEC traveling through another BEC decreases. I also study the feasibility of an experimental realization of this phenomena. Second, I analyze an experiment in which a BEC in a ring-shaped trap is stirred by a rotating barrier. The phase drop across and current flow through the barrier is measured from spiral-shaped density profiles created by interfering the BEC in the ring-shaped trap and a concentric reference BEC after release from all trapping potentials. I show that a free-particle expansion is sufficient to explain the origin of the spiral pattern and relate the phase drop to the geometry of a spiral. I also bound the expansion times for which the phase drop can be accurately determined and study the effect of inter-atomic interactions on the expansion time scales. Third, I study the dynamics of few-mode BECs when they become dynamically unstable after preparing an initial state at a saddle point of the Hamiltonian. I study the dynamics within the truncated Wigner approximation (TWA) and find that, due to phase-space mixing, the expectation value of an observable relaxes to a steady-state value. Using the action-angle formalism, we derive analytical expressions for the steady-state value and the time evolution towards this value. I apply these general results to two systems: a condensate in a double-well potential and a spin-1 (spinor) condensate. Finally, I study quantum corrections beyond the TWA in the semiclassical limit. I derive general expressions for the dynamics of an observable by using the van Vleck-Gutzwiller propagator and find that the interference of classical paths leads to non-perturbative corrections. As a case study, I consider a single-mode nonlinear oscillator; this system displays collapse and revival of observables. I find that the interference of classical paths, which is absent in the TWA, leads to revivals.
After the first experimental realization of a Bose-Einstein condensate (BEC) in 1995, BECs have become a subject of intense experimental and theoretical study. In this dissertation, I present our results on the classical and quantum dynamics of BECs at zero temperature under different scenarios. First, I consider the analog of slow light in the collision of two BECs near a Feshbach resonance. The scattering length then becomes a function of the collision energy. I derive a generalization of the Gross-Pitaevskii equation for incorporating this energy dependence. In certain parameter regimes, the group velocity of a BEC traveling through another BEC decreases. I also study the feasibility of an experimental realization of this phenomena. Second, I analyze an experiment in which a BEC in a ring-shaped trap is stirred by a rotating barrier. The phase drop across and current flow through the barrier is measured from spiral-shaped density profiles created by interfering the BEC in the ring-shaped trap and a concentric reference BEC after release from all trapping potentials. I show that a free-particle expansion is sufficient to explain the origin of the spiral pattern and relate the phase drop to the geometry of a spiral. I also bound the expansion times for which the phase drop can be accurately determined and study the effect of inter-atomic interactions on the expansion time scales. Third, I study the dynamics of few-mode BECs when they become dynamically unstable after preparing an initial state at a saddle point of the Hamiltonian. I study the dynamics within the truncated Wigner approximation (TWA) and find that, due to phase-space mixing, the expectation value of an observable relaxes to a steady-state value. Using the action-angle formalism, we derive analytical expressions for the steady-state value and the time evolution towards this value. I apply these general results to two systems: a condensate in a double-well potential and a spin-1 (spinor) condensate. Finally, I study quantum corrections beyond the TWA in the semiclassical limit. I derive general expressions for the dynamics of an observable by using the van Vleck-Gutzwiller propagator and find that the interference of classical paths leads to non-perturbative corrections. As a case study, I consider a single-mode nonlinear oscillator; this system displays collapse and revival of observables. I find that the interference of classical paths, which is absent in the TWA, leads to revivals.
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.
In this Thesis we study steady state solutions and dynamical evolutions of two– component atomic Bose–Einstein Condensates. We initially investigate the equilibrium properties of condensate mixtures in harmonic trapping potentials at zero temperature. Subsequently we simulate the coupled growth of these condensates by inclusion of damping terms. Finally, we investigate the evolution of coupled Bose gases via the so-called classical–field method. A recent experiment [D. J. McCarron et al., Phys. Rev. A, 84, 011603(R) (2011)] achieved Bose–Einstein Condensation of a two–species 87Rb–133Cs phase segregated mixture in harmonic trapping potentials. Depending on relative atom numbers of the two species, three distinct regimes of density distributions were observed. For these experimental parameters, we investigate the corresponding time–independent ground state solutions through numerical simulations of the coupled Gross–Pitaevskii equations. By including experimentally relevant shifts between the traps, we observe a range of structures including ‘ball and shell’ formations and axially/radially separated states. These are found to be very sensitive to the trap shifts. For all three experimental regimes, our numerical simulations reveal good qualitative agreement. The observed experimental profiles cannot be guaranteed to be fully equilibrated. This, coupled with the rapid sympathetic cooling of the experimental system, leads to a situation where growth may play a determining factor in the density structures formed. To investigate this further, we introduce phenomenological damping to describe the associated condensate growth/decay, revealing a range of transient structures. However, such a model always predicts the predominance of one condensate species over longer evolution times. Work undertaken by collaborators with the more elaborate Stochastic Projected Gross–Pitaevskii equations, which can describe condensate formation by coupling to a heat bath, predicts the spontaneous formation of dark–bright solitons. Motivated by this, we show how the presence of solitons can affect the condensate distribution, thus highlighting the overall dynamical role in the emerging patterns. Finally, we use classical field methods to analyse the evolution of non trapped Bose gases from strongly nonequilibrium initial distributions. The contrast between miscible (overlapping) and immiscible (phase segregated) components gives rise to important distinctions for condensate fractions and the formation of domains and vortices. In addition, splitting the particles of a single component thermalised state into two components is investigated. We then study the effects of suddenly quenching the strength of the interspecies interactions. Under suitable conditions, this quench generates isotropic vortex tangles. While this tangle subsequently decays over time, we propose how a repeat sequence of quenches at regular intervals could be employed to drive the tangle, thereby potentially providing a novel route to the generation of quantum turbulence.
In the years following the first Bose-Einstein condensation of dilute atomic gases in 1995, several experiments demonstrated many fundamental properties of two-component BEC. Subsequent to these pioneering experiments, almost a decade passed with little experimental activity, relative to investigations of single-component condensates. In somewhat of a renaissance, the initial observations of two-component BEC have been revisited in recent years, consolidating and clarifying results of the early works, and illuminating new and surprising behaviour of these rich quantum fluids. This thesis presents new experimental results on the relative phase evolution of two-component Bose-Einstein condensates (BECs). We approach our study of this system from the dual perspective of binary superfluidity and trapped atom interferometry. We observe nonequilibrium dynamics of the condensate that are deeply connected to the physics of interpenetrating superfluids. These dynamics have consequences for interferometers whose two paths are different internal spin states of the condensate atoms. Using trapped condensates for precision measurement draws upon the benefits of confined atom interferometry - long interrogation times, precise position control to measure near-field potentials - in addition to the unique coherence properties of BECs. Precision measurementwith trapped particles, however, is obstructed by the inherent inhomogeneities associated with the confinement itself. Spatial variation of the external potentials and atomic density can give rise to dephasing during the free evolution and imperfect splitting and recombination of the interferometric paths. In this thesis, we consider the dephasing of the two-componentmean field order parameter, which tends to diminish the interference signal attained with the condensate. The promise of fuhs for atom interferometry lies within using nonclassical many body states to reduce the fundamental noise limit imposed upon measuring the projection of atomic pseudospin. Before the many body quantum properties of two-component condensates can be used to their full potential, the ideal classical limit must be achieved and investigated. As in optical interferometry, the classical ideality of an interference measurement is contingent upon control of the spatialmodes - thematter wave fields, or single particle wave functions. Whilst much attention is currently focused on preparing nonclassical many body states for entanglement and squeezing of the condensate pseudospin, this thesis focuses on the ancillary yet important issue of spatial inhomogeneity of the condensate mean field. Reversal of the spatially inhomogeneous mean field evolution is demonstrated, employing the techniques of nuclear magnetic resonance. A new method is demonstrated to probe pseudospin -1/2 condensates involving microwave adiabatic passage. The technique permits simultaneous, spatially resolved absorption imaging of the spatial mode of each component. Controlled manipulation of the components into Zeeman states with different magnetic moments allows spatial separation using a Stern-Gerlach separation technique, providing a spatially resolved projection of the local condensate pseudospin from a single absorption image. We propose an interferometric technique which employs this imaging method to directly image the relative phase. The miscibility of the two components depends on the relative magnitude of the interactions between various collision partners. The boundary of this miscibility is found to be less pronounced than expected from the commonly applied results of homogeneous super-fluids.The effects of kinetic energy are shown to be significant at the threshold of miscibility, leading to a breakdown in the Thomas-Fermi approximation. Magnetic dipole transitions amongst alkali ground states are often used to prepare, manipulate, and interrogate pseudospin -1/2 condensates. We present a thorough investigation of the two-photonmagnetic dipole transitions between the |F = 1,mF = -1〉 and |F = 2,mF = +1〉 states, and the |F = 1,mF = +1〉 and |F = 2,mF = -1〉 states of 87Rb. Taking into account all 8 Zeeman levels of the electronic ground state reveals markedly different radiative shifts and effects beyond the rotating wave approximation, compared to the three-level formalism commonly used for these transitions.This has implications for Ramsey interferometry using the two-photon magnetic dipole transitions. The unchartered experimental territory offered by two-component condensates is vast, and future directions, supplementary to the work in this thesis, are elucidated.
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.
We study the collision dynamics of two Bose-Einstein condensates with their dynamical wave functions modeled by a set of coupled, time-dependent Gross-Pitaevskii equations. Beginning with an effective one-dimensional system, we identify regimes characterized by the relationship between inter- and intra-atomic interactions and the initial configuration of the system, akin to the equilibrium phase diagram of two interacting Bose condensates. We consider a dynamical setup in which two wave packets are initially at rest with a small separation about the center of an asymmetric harmonic trap. Upon release, we observe a rapid approach to dynamical equilibrium in the limits of very large and very small inter-particle repulsion, characterized by periodic transmission or reflection of the condensates as distinguishable units, whereas the intermediate, critical regime is characterized by extended transient dynamics, density fracturing, and dynamical mixing. We therefore have reason to believe that non-trivial behavior may exist during the collisions of Bose-Einstein condensates as a result of their non-linear interactions, and these effects may be observable in a suitable laboratory environment.
This thesis examines two aspects of the non-equilibrium dynamics of Bose-Einstein condensates (BECs). Our particular interests are the relaxation of near-integrable systems, and the inter-particle correlations present in BECs, both in equilibrium and following an abrupt change in the Hamiltonian. In the first part of this thesis we investigate the stability of solitary wave solutions to a system of two coupled one-dimensional non-linear Schrodinger equations (CNLSE). This system has an integrable point known as the Manakov[1] model which has analytic dark-bright soliton solutions[2]. We break the integrability of this model by varying the inter-species interaction strength, and employ variational and numerical methods in order to find dark-bright solitary wave solutions, away from the integrable point. The time evolution of these states is calculated via numerical integration of the CNLSE, which allows us to assess the stability of these solutions in isolation and to quantify their robustness against collisions with a dark soliton. We find that there is a broad region of the parameter space in which solutions for black-bright (stationary) solitary waves can be found. In this domain, the integrability breaking is revealed only during collisions. Prior to a collision, the numerical solutions are indistinguishable from true solitons, however, the collisional stability of these solitary waves is significantly affected by the extent to which integrability is broken. We find that there is a smooth transition between the non-dispersive particle-like nature of soliton interactions for integrable systems, and the destructive collisions which take place when the system is perturbed far from the integrable point. In the case of moving (grey-bright) solitary waves, we find a more restricted region of the parameter space which admits solitary wave solutions. In this domain we are able to find approximately stable solitary wave solutions using a variational ansatz. The collisional stability of these grey-bright solitary waves is also found to depend upon the extent to which integrability is broken. We conclude that the observation of stable, long-lived dark-bright solitary waves is not sufficient to indicate near-integrability of the two-component system. However, near-integrability may inferred if dark-bright solitary waves are observed to survive multiple collisions. These results are not only of theoretical interest, but may also inform future experiments which explore the behaviour of solitary waves in two-component Bose-Einstein condensates. Following this we investigate the properties of magnetic solitary wave solutions to the CNLSE. We find that the analytic ‘magnetic soliton’ solutions derived in reference [3] are dynamically unstable and relax to non-stationary magnetic solitary wave states over a short timescale. We investigate the collisional dynamics of these magnetic solitary waves and find that these excitations are remarkably robust. In the second part of this thesis we calculate the fluctuations in the number difference between two sub-regions of a Bose-Einstein condensate. This study is motivated by recent experiments performed by the Truscott group at the Australian National University (ANU), which indicate sub- Poissonian fluctuations even for condensates with large thermal depletion. This result cannot easily be reconciled with existing literature, which predicts strongly super-Poissonian fluctuations below the critical temperature (Tc). We develop one-dimensional, and three-dimensional, models to describe these experiments and employ the Bogoliubov formalism to explore the quantum fluctuations of the condensate, which have not previously been studied in detail. We calculate the relative number fluctuations between two halves of the condensate, and also between multiple pairs of spatial bins. For harmonically confined condensates at zero temperature, quantum correlations result in a suppression of the relative number fluctuations between the two halves of the condensate below the shot-noise level. However, the presence of even a small thermal component results in strongly super-Poissonian fluctuations below Tc; this is in agreement with the existing theoretical and experimental results. In addition, we find that the fluctuations of the harmonically confined condensate are non-uniform, with the relative number fluctuations being peaked towards the edges of the condensate. The recent studies of the Truscott group measure the fluctuations after free expansion. We therefore calculate the propagation of correlations following a sudden removal of the harmonic confinement. We find that free expansion results in a strong decrease in the relative number fluctuations. Condensates which exhibit super-Poissonian fluctuations when harmonically confined may therefore display number- squeezing after expansion. A second set of experiments by the Truscott group measured the density fluctuations present in small ‘atom laser’ pulses of atoms out-coupled from their BEC. We also model this procedure and find that immediately after the out-coupling process, the number fluctuations of the out-coupled atomic pulse are approximately Poissonian, irrespective of the correlations present in the atomic reservoir. However, the competing effects of interactions between the atomic pulse and the condensate reservoir, and free expansion of the pulse after out-coupling, can result in the measured fluctuations being either number squeezed or super-Poissonian. Our results will guide future experiments planned by the Truscott group.