A quantum vortex is an excitation of a superfluid that carries angular momentum, in which the superfluid is expelled from and circulates around a region known as the vortex core. A quantum vortex dipole consists of two bound quantum vortices with opposite circulation, which together carry linear momentum. The topic of this thesis is two-fold, with a common theme of quantum vortex dipoles.\nIn the first part of the thesis, the motion of a quantum vortex dipole incident upon a step-change in the background superfluid density of an otherwise uniform two-dimensional Bose-Einstein condensate is investigated, both analytically and numerically.\nDue to the conservation of fluid momentum and energy, the incident and refracted angles of the dipole satisfy a relation analogous to Snell’s law, when crossing the interface between regions of different density. The predictions of the analogue Snell’s law relation are confirmed for a wide range of incident angles by systematic numerical simulations of the Gross-Pitaevskii equation (GPE). Near the critical angle for total internal reflection, we identify a regime of anomalous Snell’s law behaviour where the finite size of the dipole causes transient capture by the interface. Remarkably, despite the extra complexity of the surface interaction, the incoming and outgoing dipole paths obey Snell’s law.\nIn the second part of the thesis, the point-vortex model is used to find the stationary states of an arbitrary neutral number of vortices of alternating sign arranged on a ring in the simply bounded domain. These stationary states are collectively referred to as necklace states, due to the symmetry of their appearance. Curiously, the necklace states are found to have a simple relation to the metallic means, which are generalisations of the golden ratio. The necklace state is numerically simulated with different initial perturbations in both the point-vortex model and the GPE with a comparison of the results. Perturbations from the necklace solution evolve differently in the point-vortex model and GPE due to the healing length scale in the GPE enforcing a departure from the point-vortex model predictions. We introduce a simple numerical scheme to compensate for the healing length and confirm the necklace as a steady of the GPE. We find that the agreement between the point-vortex model and GPE for evolution of significant perturbations persists for longer as the healing length reduces compared to system size.
A quantum vortex is an excitation of a superfluid that carries angular momentum, in which the superfluid is expelled from and circulates around a region known as the vortex core. A quantum vortex dipole consists of two bound quantum vortices with opposite circulation, which together carry linear momentum. The topic of this thesis is two-fold, with a common theme of quantum vortex dipoles.\nIn the first part of the thesis, the motion of a quantum vortex dipole incident upon a step-change in the background superfluid density of an otherwise uniform two-dimensional Bose-Einstein condensate is investigated, both analytically and numerically.\nDue to the conservation of fluid momentum and energy, the incident and refracted angles of the dipole satisfy a relation analogous to Snell’s law, when crossing the interface between regions of different density. The predictions of the analogue Snell’s law relation are confirmed for a wide range of incident angles by systematic numerical simulations of the Gross-Pitaevskii equation (GPE). Near the critical angle for total internal reflection, we identify a regime of anomalous Snell’s law behaviour where the finite size of the dipole causes transient capture by the interface. Remarkably, despite the extra complexity of the surface interaction, the incoming and outgoing dipole paths obey Snell’s law.\nIn the second part of the thesis, the point-vortex model is used to find the stationary states of an arbitrary neutral number of vortices of alternating sign arranged on a ring in the simply bounded domain. These stationary states are collectively referred to as necklace states, due to the symmetry of their appearance. Curiously, the necklace states are found to have a simple relation to the metallic means, which are generalisations of the golden ratio. The necklace state is numerically simulated with different initial perturbations in both the point-vortex model and the GPE with a comparison of the results. Perturbations from the necklace solution evolve differently in the point-vortex model and GPE due to the healing length scale in the GPE enforcing a departure from the point-vortex model predictions. We introduce a simple numerical scheme to compensate for the healing length and confirm the necklace as a steady of the GPE. We find that the agreement between the point-vortex model and GPE for evolution of significant perturbations persists for longer as the healing length reduces compared to system size.
A quantum vortex dipole, comprised of a closely bound pair of vortices of equal strength with opposite circulation, is a spatially localized travelling excitation of a planar superfluid that carries linear momentum, suggesting a possible analogy with ray optics. We investigate numerically and analytically the motion of a quantum vortex dipole incident upon a step-change in the background superfluid density of an otherwise uniform two-dimensional Bose-Einstein condensate. Due to the conservation of fluid momentum and energy, the incident and refracted angles of the dipole satisfy a relation analogous to Snell’s law, when crossing the interface between regions of different density. The predictions of the analogue Snell’s law relation are confirmed for a wide range of incident angles by systematic numerical simulations of the Gross-Piteavskii equation. Near the critical angle for total internal reflection, we identify a regime of anomalous Snell’s law behaviour where the finite size of the dipole causes transient capture by the interface. Remarkably, despite the extra complexity of the interface interaction, the incoming and outgoing dipole paths obey Snell’s law.
In this thesis, we analytically and numerically study two different scenarios of vortex dynamics in a two-dimensional superfluid. The first system we explore is an expanding cluster of like-sign vortices, which recently suggested it could explain observations in a superfluid optomechanics experiment. We employ a recent theoretical model that approximates large collections of vortices as a fluid and present an analytic solution for its dynamics. Analysis suggests that as a cluster expands due to thermal friction, it forms a Rankine vortex (i.e., a top-hat distribution) that expands self-similarly in time regardless of initial vorticity distribution. Numerical simulations of point vortices confirm this universal expanding regime and further demonstrates vortex lattice formation that goes beyond the vortex fluid theory. We analyse correlations that emerge in the resulting lattices and observe decaying disorder as the cluster expands. Furthermore, we analyse experimental data from a quasi-two-dimensional 87Rb Bose-Einstein condensate, which supports the vortex fluid theory’s predictions.The second scenario considers the dynamics of a scattering event where a single quantum vortex interacts with a uniform, circular pinning potential in a steady two-dimensional superfluid flow. We first develop an approximate hydrodynamical model to understand the interaction between the vortex and the pinning potential in the conservative regime and find it is in good agreement with the mean-field Gross-Pitaevskii equation in the low-velocity limit. We further solve the Gross-Pitaevskii equation for higher background velocities to study how the vortex moves around pinning sites, interacts with it, and ultimately becomes pinned. We uncover four different regimes of vortex pinning, each marked with a characteristic emission of sound energy. Finally, we study how vortices become unpinned from a pinning site, revealing large vortex dipoles form to stabilise the system.Our work in this thesis takes significant steps towards advancing knowledge of both the collective dynamics of vortices and their individual interaction with defects upon surface substrates. Our numerical and analytical results opens up new avenues of exploration be it further understanding dissipation mechanisms in two-dimensional superfluids, or elucidating the nature of vortex pinning in systems such as neutron stars.
This thesis investigates the behaviour of a quantum impurity coupled to a Bose–Einstein condensate. The condensate dynamics are described by the Gross–Pitaevskii equation, while the impurity is modelled by a coupled Schrödinger equation, forming a Gross–Clark-system. The work focuses on impurity localisation, vortex structure, and the dynamical consequences of impurity-condensate coupling in effectively two-dimensional geometries. The thesis begins by establishing the Gross–Pitaevskii model and its dimensionless formulation, together with the structure and stability of quantised vortices. Both numerical and approximate analytical descriptions of vortex profiles are developed and used throughout the thesis. The localisation of impurities is then examined in a range of static settings, including idealised trapping potentials and vortex cores, revealing how condensate density depletion can generate effective confinement for the impurity. Self-consistent impurity-vortex bound states are obtained by solving coupled stationary equations. The spectrum of localised impurity states is analysed, and the dependence of localisation on interaction parameters and angular momentum is characterised. These static results provide the foundation for subsequent dynamical studies. Real-time simulations are used to investigate the stability of impurity-bound configurations. For singly quantised vortices, impurity ground states are shown to be dynamically robust, while excited impurity states relax under time evolution despite remaining localised. A quantitative overlap measure is introduced to characterise the dynamical stability of impurity states. The thesis then examines multiply quantised vortices, demonstrating that impurity coupling can modify vortex dynamics and, in certain cases, suppress intrinsic instabilities. Overall, this work provides a systematic analysis of impurity localisation and impurity–vortex interactions in Bose-Einstein condensates, combining theoretical modelling, numerical methods, and dynamical simulations to clarify the role of impurities in coupled quantum fluids.
In this thesis we solve the Gross-Pitaevskii equation numerically in order to model the response of trapped Bose-Einstein condensed gases to perturbations by electromagnetic fields. First, we simulate output coupling of pulses from the condensate and compare our results to experiments. The excitation and separation of eigenmodes on flow through a constriction is also studied. We then move on to the main theme of this thesis: the important subject of quantised vortices in Bose condensates, and the relation between Bose-Einstein condensation and superfluidity. We propose methods of producing vortex pairs and rings by controlled motion of objects. Full three-dimensional simulations under realistic experimental conditions are performed in order to test the validity of these ideas. We link vortex formation to drag forces on the object, which in turn is connected with energy transfer to the condensate. We therefore argue that vortex formation by moving objects is intimately related to the onset of dissipation in superfluids. We discuss this idea in the context of a recent experiment, using simulations to provide evidence of vortex formation in the experimental scenario. Superfluidity is also manifest in the property of persistent currents, which is linked to vortex stability and dynamics. We simulate vortex line and ring motion, and find in both cases precessional motion and thermodynamic instability to dissipation. Strictly speaking, the Gross-Pitaevskii equation is valid only for temperatures far below the BEC transition. We end the thesis by describing a simple finite-temperature model to describe mean-field coupling between condensed and non-condensed components of the gas. We show that our hybrid Monte-Carlo/FFT technique can describe damping of the lowest energy excitations of the system. Extensions to this model and future research directions are discussed in the conclusion. (author)
The dynamics and interaction of quantized vortices in Bose-Einstein condensates (BECs) are investigated by using the two-dimensional Gross-Pitaevskii equation (GPE) with/without an angular momentum rotation term. If all vortices have the same winding number, they would rotate around the trap center but never collide. In contrast, if the winding numbers are different, their interaction highly depends on the initial distance between vortex centers. The analytical results are presented to describe the dynamics of the vortex centers when β=0. While if β≠0, there is no analytical result but some conclusive numerical findings are provided for the further understanding of vortex interaction in BECs. Finally, the dynamic laws describing the relation of vortex interaction in nonrotating and rotating BECs are presented.
Recently, Freilich et al. [Science 329, 1182 (2010)] experimentally discovered stationary states of vortex dipoles, pairs of vortices of opposite circulation, in dilute Bose-Einstein condensates. To explain their observations, we perform simulations based on the Gross-Pitaevskii equation and obtain excellent quantitative agreement on the size of the stationary dipole. We also investigate how their imaging method, in which atoms are repeatedly extracted from a single condensate, affects the vortex dynamics. We find that it mainly induces isotropic size oscillations of the condensate without otherwise disturbing the vortex trajectories. Thus, the imaging technique appears to be a promising tool for studying real-time superfluid dynamics.