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

Bright solitary waves and non-equilibrium dynamics in atomic Bose-Einstein condensates

T. P. Billam

Published Jan 1, 2012• 0 citations
OpenAlex

Abstract

In this thesis we investigate the static properties and non-equilibrium dynamics of bright solitary waves in atomic Bose-Einstein condensates in the zero-temperature limit, and we investigate the non-equilibrium dynamics of a driven atomic Bose-Einstein condensate at finite temperature. Bright solitary waves in atomic Bose-Einstein condensates are non-dispersive and soliton-like matter-waves which could be used in future atom-interferometry exper-iments. Using the mean-field, Gross-Pitaevskii description, we propose an exper-imental scheme to generate pairs of bright solitary waves with controlled velocity and relative phase; this scheme could form an important part of a future atom in-terferometer, and we demonstrate that it can also be used to test the validity of the mean-field model of bright solitary waves. We also develop a method to quanti-tatively assess how soliton-like static, three-dimensional bright solitary waves are; this assessment is particularly relevant for the design of future experiments. In reality, the non-zero temperatures and highly non-equilibrium dynamics occur-ring in a bright solitary wave interferometer are likely to necessitate a theoretical description which explicitly accounts for the non-condensate fraction. We show that a second-order, number-conserving description offers a minimal self-consistent treatment of the relevant condensate – non-condensate interactions at low temper-atures and for moderate non-condensate fractions. We develop a method to obtain a fully-dynamical numerical solution to the integro-differential equations of motion of this description, and solve these equations for a driven, quasi-one-dimensional test system. We show that rapid non-condensate growth predicted by lower-order descriptions, and associated with linear dynamical instabilities, can be damped by the self-consistent treatment of interactions included in the second-order descrip-tion. 1 Bright solitary waves and non-equilibrium dynamics in atomic Bose-Einstein condensates

Related papers

Non-equilibrium dynamics of Bose Einstein condensates

James Edward Kennett Mills

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.

•Oct 6, 2017•0 citations

Bright solitary matter waves: formation, stability and interactions

T. P. Billam, A. L. Marchant, Simon L. Cornish +2 more

In recent years, bright soliton-like structures composed of gaseous Bose-Einstein condensates have been generated at ultracold temperature. The experimental capacity to precisely engineer the nonlinearity and potential landscape experienced by these solitary waves offers an attractive platform for fundamental study of solitonic structures. The presence of three spatial dimensions and trapping implies that these are strictly distinct objects to the true soliton solutions. Working within the zero-temperature mean-field description, we explore the solutions and stability of bright solitary waves, as well as their interactions. Emphasis is placed on elucidating their similarities and differences to the true bright soliton. The rich behaviour introduced in the bright solitary waves includes the collapse instability and symmetry-breaking collisions. We review the experimental formation and observation of bright solitary matter waves to date, and compare to theoretical predictions. Finally we discuss the current state-of-the-art of this area, including beyond-mean-field descriptions, exotic bright solitary waves, and proposals to exploit bright solitary waves in interferometry and as surface probes.

arXiv:1209.0560•Sep 4, 2012•21 citations

Dynamics of bright solitary matter-waves

Timothy P. Wiles

This thesis describes the formation of a 85Rb Bose-Einstein condensate and the subsequent creation of a bright solitary matter-wave in a quasi-one-dimensional optical waveguide, with experiments investigating the dynamics of a solitary wave in comparison with a repulsive Bose-Einstein condensate. In the final chapters of this thesis, progress towards the interaction of a solitary wave with a narrow barrier and an attractive atom-surface potential is presented. Beyond the above, a review of recent soliton and solitary wave theory is presented from the perspective of an experimentalist, culminating in the numerical analysis of a variety of key areas, namely, quantum reflection, solitary wave size and solitary wave profile. The modelling and experimental results relating to the merging of ultracold gases of 85Rb and 87Rb, in order to create isotopic mixtures, is also described within this thesis. Such a scheme could be used as an initial step in the process of forming molecules or undertaking sympathetic cooling. Finally, the creation of a complex LabVIEW based experimental control system is also described within.

•Jan 1, 2013•0 citations

Two-component Bose-Einstein condensates :equilibria and dynamics at zero temperature and beyond

RW Pattinson

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.

•Jan 1, 2014•4 citations

Controlled formation and reflection of a bright solitary matter-wave

A. L. Marchant, T. P. Billam, Timothy P. Wiles +3 more

Bright solitons are non-dispersive wave solutions, arising in a diverse range of nonlinear, one-dimensional systems, including atomic Bose-Einstein condensates with attractive interactions. In reality, cold-atom experiments can only approach the idealized one-dimensional limit necessary for the realization of true solitons. Nevertheless, it remains possible to create bright solitary waves, the three-dimensional analogue of solitons, which maintain many of the key properties of their one-dimensional counterparts. Such solitary waves offer many potential applications and provide a rich testing ground for theoretical treatments of many-body quantum systems. Here we report the controlled formation of a bright solitary matter-wave from a Bose-Einstein condensate of (85)Rb, which is observed to propagate over a distance of ∼1.1 mm in 150 ms with no observable dispersion. We demonstrate the reflection of a solitary wave from a repulsive Gaussian barrier and contrast this to the case of a repulsive condensate, in both cases finding excellent agreement with theoretical simulations using the three-dimensional Gross-Pitaevskii equation.

arXiv:1301.5759•Nature Communications•May 14, 2013•183 citations

Theoretical studies of bright solitons in trapped atomic Bose-Einstein condensates

Andrew David Martin

Bright solitary-waves may be created in dilute Bose-Einstein condensates of at tractively interacting atoms in one dimensional regimes. In integrable systems, such solitary waves are particle-like objects called solitons. We investigate the consequences of non-integrability on the solitary waves in trapped Bose-Einstein condensates caused by an axial harmonic trap, and non-integrability caused by three dimensional effects. To analyse the soliton-like nature of the solitary-waves in an axial harmonic trap, a particle analogy for the solitary-waves is formulated. Exact Soliton solutions exist in the absence of an external trapping potential, which behave in a particle-like manner, and we find the particle analogy we employ to be a good model also when a harmonic trapping potential is present up to a gradual shift in the trajectories when the harmonic trap period is short compared with the, collision time of the solitons. We find that the collision time of the solions is dependent on the relative phase of the solitons as they collide. In the case of two solitons, the particle model is integrable, and the dynamics are completely regular. In the case of a system of two solitary waves of equal norm, the solitons are shown to retain their phase difference for repeated collisions. The extension to three particles supports both regular and chaotic regimes. The trajectory shift observed for two solitons carrier over to the case of three solitons. This shift aside, the agreement between the particle model and the wave dynamics remains good, even in chaotic regimes. We predict that these chaotic regimes will be an indicator of rapid depletion of the condensate due to quantum transitions of the condensate particles into non-condensate modes. To analyse the residual effects of the three dimensional nature of the solitary waves, we use a nonlinear Schrödinger equation with an additional quintic term. We perform variational calculations, and confirm the collapse of a soliton when the number of particles contained therein is increased past a critical number. We investigate the effects of varying the axial trap frequency and scattering length on the critical number. We propose a method to model particle exchange between solitons by extending the variational treatment to two solitons.

•Jan 1, 2008•1 citation
In this thesis we investigate the static properties and non-equilibrium dynamics of bright solitary waves in atomic Bose-Einstein condensates in the zero-temperature limit, and we investigate the non-equilibrium dynamics of a driven atomic Bose-Einstein condensate at finite temperature. Bright solitary waves in atomic Bose-Einstein condensates are non-dispersive and soliton-like matter-waves which could be used in future atom-interferometry exper-iments. Using the mean-field, Gross-Pitaevskii description, we propose an exper-imental scheme to generate pairs of bright solitary waves with controlled velocity and relative phase; this scheme could form an important part of a future atom in-terferometer, and we demonstrate that it can also be used to test the validity of the mean-field model of bright solitary waves. We also develop a method to quanti-tatively assess how soliton-like static, three-dimensional bright solitary waves are; this assessment is particularly relevant for the design of future experiments. In reality, the non-zero temperatures and highly non-equilibrium dynamics occur-ring in a bright solitary wave interferometer are likely to necessitate a theoretical description which explicitly accounts for the non-condensate fraction. We show that a second-order, number-conserving description offers a minimal self-consistent treatment of the relevant condensate – non-condensate interactions at low temper-atures and for moderate non-condensate fractions. We develop a method to obtain a fully-dynamical numerical solution to the integro-differential equations of motion of this description, and solve these equations for a driven, quasi-one-dimensional test system. We show that rapid non-condensate growth predicted by lower-order descriptions, and associated with linear dynamical instabilities, can be damped by the self-consistent treatment of interactions included in the second-order descrip-tion. 1 Bright solitary waves and non-equilibrium dynamics in atomic Bose-Einstein condensates
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.
In recent years, bright soliton-like structures composed of gaseous Bose-Einstein condensates have been generated at ultracold temperature. The experimental capacity to precisely engineer the nonlinearity and potential landscape experienced by these solitary waves offers an attractive platform for fundamental study of solitonic structures. The presence of three spatial dimensions and trapping implies that these are strictly distinct objects to the true soliton solutions. Working within the zero-temperature mean-field description, we explore the solutions and stability of bright solitary waves, as well as their interactions. Emphasis is placed on elucidating their similarities and differences to the true bright soliton. The rich behaviour introduced in the bright solitary waves includes the collapse instability and symmetry-breaking collisions. We review the experimental formation and observation of bright solitary matter waves to date, and compare to theoretical predictions. Finally we discuss the current state-of-the-art of this area, including beyond-mean-field descriptions, exotic bright solitary waves, and proposals to exploit bright solitary waves in interferometry and as surface probes.
This thesis describes the formation of a 85Rb Bose-Einstein condensate and the subsequent creation of a bright solitary matter-wave in a quasi-one-dimensional optical waveguide, with experiments investigating the dynamics of a solitary wave in comparison with a repulsive Bose-Einstein condensate. In the final chapters of this thesis, progress towards the interaction of a solitary wave with a narrow barrier and an attractive atom-surface potential is presented. Beyond the above, a review of recent soliton and solitary wave theory is presented from the perspective of an experimentalist, culminating in the numerical analysis of a variety of key areas, namely, quantum reflection, solitary wave size and solitary wave profile. The modelling and experimental results relating to the merging of ultracold gases of 85Rb and 87Rb, in order to create isotopic mixtures, is also described within this thesis. Such a scheme could be used as an initial step in the process of forming molecules or undertaking sympathetic cooling. Finally, the creation of a complex LabVIEW based experimental control system is also described within.
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.
Bright solitons are non-dispersive wave solutions, arising in a diverse range of nonlinear, one-dimensional systems, including atomic Bose-Einstein condensates with attractive interactions. In reality, cold-atom experiments can only approach the idealized one-dimensional limit necessary for the realization of true solitons. Nevertheless, it remains possible to create bright solitary waves, the three-dimensional analogue of solitons, which maintain many of the key properties of their one-dimensional counterparts. Such solitary waves offer many potential applications and provide a rich testing ground for theoretical treatments of many-body quantum systems. Here we report the controlled formation of a bright solitary matter-wave from a Bose-Einstein condensate of (85)Rb, which is observed to propagate over a distance of ∼1.1 mm in 150 ms with no observable dispersion. We demonstrate the reflection of a solitary wave from a repulsive Gaussian barrier and contrast this to the case of a repulsive condensate, in both cases finding excellent agreement with theoretical simulations using the three-dimensional Gross-Pitaevskii equation.
Bright solitary-waves may be created in dilute Bose-Einstein condensates of at tractively interacting atoms in one dimensional regimes. In integrable systems, such solitary waves are particle-like objects called solitons. We investigate the consequences of non-integrability on the solitary waves in trapped Bose-Einstein condensates caused by an axial harmonic trap, and non-integrability caused by three dimensional effects. To analyse the soliton-like nature of the solitary-waves in an axial harmonic trap, a particle analogy for the solitary-waves is formulated. Exact Soliton solutions exist in the absence of an external trapping potential, which behave in a particle-like manner, and we find the particle analogy we employ to be a good model also when a harmonic trapping potential is present up to a gradual shift in the trajectories when the harmonic trap period is short compared with the, collision time of the solitons. We find that the collision time of the solions is dependent on the relative phase of the solitons as they collide. In the case of two solitons, the particle model is integrable, and the dynamics are completely regular. In the case of a system of two solitary waves of equal norm, the solitons are shown to retain their phase difference for repeated collisions. The extension to three particles supports both regular and chaotic regimes. The trajectory shift observed for two solitons carrier over to the case of three solitons. This shift aside, the agreement between the particle model and the wave dynamics remains good, even in chaotic regimes. We predict that these chaotic regimes will be an indicator of rapid depletion of the condensate due to quantum transitions of the condensate particles into non-condensate modes. To analyse the residual effects of the three dimensional nature of the solitary waves, we use a nonlinear Schrödinger equation with an additional quintic term. We perform variational calculations, and confirm the collapse of a soliton when the number of particles contained therein is increased past a critical number. We investigate the effects of varying the axial trap frequency and scattering length on the critical number. We propose a method to model particle exchange between solitons by extending the variational treatment to two solitons.