Over the past decades ultracold atomic gases have emerged as a highly controllable platform for probing quantum many-body phenomena, including Bose-Einstein condensation and superfluidity. These systems are particularly suited to exploring the physics of open quantum systems where interactions with the environment, traditionally considered a source of decoherence, can be engineered to reveal new types of non-equilibrium behaviour. This thesis focuses on the dissipative Josephson array, a setup of ultracold atoms subjected to a controlled loss mechanism. The primary goal is to better understand a series of intriguing experimental observations by analysing a range of theoretical models. As such, the focus is put on the dissipative Bose-Hubbard model, which isolates the lattice dynamics and provides insight into the system’s steady-state properties. Our analysis reveals qualitative similarities to the rather simple picture of a Kerr resonator, but also points to significant differences. To refine our understanding, we develop an effective model for the internal reservoirs based on incoherent driving, which accurately captures the observed non-equilibrium steady states. Additionally, signs of a previously unexplored region in phase space are observed. The work culminates in the identification of this phase: the standing dark soliton. Using a detailed 3D model, we reconstruct the system’s full phase diagram and perform a direct comparison with experimental data, demonstrating the relevance of the dark soliton to the observed phenomena.
Over the past decades ultracold atomic gases have emerged as a highly controllable platform for probing quantum many-body phenomena, including Bose-Einstein condensation and superfluidity. These systems are particularly suited to exploring the physics of open quantum systems where interactions with the environment, traditionally considered a source of decoherence, can be engineered to reveal new types of non-equilibrium behaviour. This thesis focuses on the dissipative Josephson array, a setup of ultracold atoms subjected to a controlled loss mechanism. The primary goal is to better understand a series of intriguing experimental observations by analysing a range of theoretical models. As such, the focus is put on the dissipative Bose-Hubbard model, which isolates the lattice dynamics and provides insight into the system’s steady-state properties. Our analysis reveals qualitative similarities to the rather simple picture of a Kerr resonator, but also points to significant differences. To refine our understanding, we develop an effective model for the internal reservoirs based on incoherent driving, which accurately captures the observed non-equilibrium steady states. Additionally, signs of a previously unexplored region in phase space are observed. The work culminates in the identification of this phase: the standing dark soliton. Using a detailed 3D model, we reconstruct the system’s full phase diagram and perform a direct comparison with experimental data, demonstrating the relevance of the dark soliton to the observed phenomena.
We show that a Bose-Hubbard model with local loss exhibits a first-order dark-state phase transition between a dark soliton and a uniform superfluid. Semiclassical simulations and stability analysis reveal that the soliton acts as a many-body dark state, whose suppressed fluctuations modify the transition relative to the comparable phenomenon of optical bistability in Kerr resonators. This mechanism quantitatively captures the bistable phase boundary observed in the experiment of Labouvie [], resolving long-standing discrepancies in previous theoretical modeling. Our results identify a general route by which spatial coherence and engineered loss can generate nonequilibrium phase transitions in ultracold gases.
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.
This thesis presents several experiments investigating the regimes of behaviour of a proto-typical open quantum many-body system far away from equilibrium. The experimental platform is based on ultracold atoms laser-excited to Rydberg states, which we engineer to emulate a strongly interacting, driven-dissipative quantum spin system. The high degree of control over the relevant microscopic processes and their parameters, as well as the ability to widely tune interaction and driving strengths provides the means to address fundamental questions on how different regimes of dynamical behaviour emerge in complex open quantum systems and how they can be characterised. \n \nIn the first part, we discuss how control over the relevant properties of Rydberg spin systems, like temperature, density, long evolution times, strong interactions, tunable driving, and dissipation can be achieved. We introduce the new experimental apparatus developed during this thesis, which combines both single-photon and two-photon Rydberg excitation schemes with high atom densities and long excitation times using an optical dipole trap. Additionally, we investigate theoretically how long coherence times exceeding motional timescales combined with strong many-body interactions can be achieved in such a system. We identify an optimal parameter regime for two-photon excitation where such conditions can be generated with coherence times and interaction strengths comparable to what is achievable for single-photon excitation. \n \nIn the second part we explore the macroscopic non-equilibrium behaviour of our system and devise suitable observables for characterising different regimes of behaviour. We discover that the slow atom loss inherent to our system provides a convenient observable for the many-body state of the system. Focussing on evolution times where the effect of the atom loss on the evolution is small, we discover that the rate of atom loss exhibits powerlaw scaling with the driving strength over several orders of magnitude. The measured scaling exponents reveal the non-equilibrium phase structure of the many-body system and allow us to distinguish dissipation-dominated, paramagnetic and critical regimes, as well as an instability. In the observed critical regime, collectively enhanced driving dominates over dissipation, leading to scaling associated to the critical point of the non-dissipative equilibrium Rydberg spin system. Based on the known microscopic processes of our system, we perform classical many-body rate equation simulations, which agree well with the observed phase structure as well as the position and associated scaling exponents of the individual regimes. These findings open up new means to study and classify out of equilibrium systems based on slow particle loss and powerlaw scaling. \n \nIn the last part of this thesis we reveal a second type of criticality in our system, which is genuinely non-equilibrium in character and arises from an absorbing state phase transition. The critical state of this phase transition dominates the system evolution at late times where particle loss is no longer negligible. We identify self-organisation as the mechanism driving the system to this critical state, which we observe through powerlaw scaling of the non-equilibrium steady-state. These experiments establish Rydberg atoms as a well controlled platform for implementing and exploring models of absorbing state phase transitions and self-organised criticality with unprecedented access to the underlying microscopic properties of the system.
Bose-Einstein condensates (BECs), a quantum state of matter formed when bosonic atoms are cooled close to absolute zero, have become the premier platform for investigating fundamental physics with atomic vapours. Experiments on Bose-Einstein condensates now achieve exquisite control over many aspects of the system, including interactions, trapping potential, and dynamics. This has precipitated a new wave of research into many-body quantum phenomena and, in particular, solitons. These structures are fundamental excitations of an interacting non-linear medium, of interest to a multitude of scientific disciplines from non-linear optics to financial markets. The highly controllable environment of BECs form an attractive playground for the study of solitons allowing the non-linearity to be dynamically tuned, facilitating deeper investigations into these structures. Consistently generating and analysing solitons in BEC experiments continues to be problematic. In particular, the non-linear dynamics of BECs, though required for the generation of solitons, produce particularly challenging control and optimization problems. These control problems must be solved before further investigations into the fundamental physics of soliton dynamics can be answered. This thesis makes three important advances in the control and measurement of BECs that will lead to better generation and observation of solitons. (1) a theoretical model for a control scheme capable of highly precise wavefunction engineering, (2) the experimental implementation of a machine learning algorithm for online optimisation, and (3) a continuous non-destructive imaging system capable of directly observing soliton dynamics in real-time. Together, these advances provide a suite of tools for manipulating and exploiting solitons in Bose-Einstein condensates. A novel technique was developed theoretically, offering control of the macroscopic wavefunction of a Bose-Einstein condensate with unprecedented spatial resolution and speed. The ability to control the atomic wavefunction at the fundamental length scale is key to the advancement of many quantum technologies such as quantum simulators. The magnetic resonance control scheme is demonstrated through simulation of a 87Rb condensate with the exemplar model generating a single dark soliton with corresponding π phase kink. The soliton represents a structure at the fundamental length scale of the system, and demonstrates the potential of the scheme for precision state engineering. The scheme is extended to generate higher-order soliton modes which are yet to be experimentally realised. A machine learning algorithm based on Gaussian processes was developed and implemented on the evaporative cooling stage of the production of a 87Rb Bose-Einstein condensate, successfully demonstrating fast optimisation to condensation. The Gaussian process develops a statistical model based on the data that enables the characterisation of the relationship between the experimental controls and resultant quality of the BEC. This relationship is often obfuscated through technical details of the apparatus, frustrating the use of theoretical models to design optimal evaporation ramps. These models often only consider ergodic dynamics with two-body s-wave interactions and no other loss rates with better ramps likely exploiting more complex interactions. The internal model generated from the Gaussian process utilised uncertainty in measured data, making the optimisation more robust to experimental noise than alternate methods. The algorithm is shown to produce high quality Bose-Einstein condensates in 10 times fewer experimental iterations than previously used online optimisation techniques. By exploiting information on the sensitivity of each control, the model can be used to aid experimental design. The convergence of the optimisation is further improved by eliminating a superfluous parameter identified by the model. The general usefulness of machine learning compared with bespoke optimisation algorithms has seen machine learning approach ubiquity. Finally, an experimentally straightforward technique for continuous non-destructive imaging of matter-wave solitons was developed and implemented, facilitating measurements of stochastic phenomena. The technique is readily practicable on any ultracold atom experiment with an existing absorption imaging system, simply requiring the probe laser be far-detuned from resonance. With a signal-to-noise of ∼ 33 at 1.25 GHz detuning, the technique is capable of producing 100 images with no observable heating or atom loss. Coupled with a fast optical phase locked loop, the technique can be used in conjunction with absorption imaging to generate a series of non-destructive images followed by a final high signal-to-noise absorption image solely through moving the laser on resonance for the final image. The high performance and utility of this imaging setup make it a powerful tool for ultra-cold atom experiments.
Non-equilibrium (NEQ) many-body dynamics are ubiquitous, ranging from glassy systems to the emergence of structure in the early universe. Ultracold atomic gases provide an ideal, highly controllable platform for studying these phenomena, and this thesis presents a comprehensive investigation into the NEQ dynamics of coupled bilayer two-dimensional (2D) Bose gases. In this thesis, we numerically investigate Josephson effects in bilayer 2D superfluids using classical-field simulations. We show that the system exhibits damped Josephson oscillations and macroscopic quantum self-trapping, and we discuss the associated excitations that emerge in the system. We then experimentally study the phase-ordering dynamics following a coupling quench that drives the system from a disordered to an ordered phase, realized using a well-controlled multiple-radio-frequency (MRF) dressed double-well potential. The resulting NEQ dynamics are probed using matter-wave interferometry, which provides direct access to correlation functions and vortex excitations of the relative-phase mode. Our results demonstrate that the system undergoes universal dynamics in the form of diffusion-type coarsening, in excellent agreement with classical-field numerical simulations. This work broadens the understanding of NEQ and universal dynamics, and the methods presented here serve as a novel and highly tunable approach for exploring a wide range of NEQ dynamics.
The Bose-Hubbard model provides a foundational paradigm for exploring quantum phase transitions in many-body systems. This thesis presents a systematic investigation into this model and its extensions, designed to provide a clear path for understanding both the rich underlying physics and the sophisticated numerical methods required for its study. Our approach is sequential: we begin with the standard one-dimensional Bose-Hubbard model, introduce complexity via the extended model with long-range interactions, and conclude by analyzing the impact of dissipation in a cavity-coupled open quantum system.To navigate these regimes, we implement and apply a suite of powerful computational techniques, including Density Matrix Renormalization Group (DMRG), finite-temperature Quantum Monte Carlo (QMC), and Mean-Field Theory (MFT). A significant component of this work is the detailed exposition of the implementation and logic of these numerical methods. Through their comparative application, we precisely characterize the equilibrium phase diagrams and benchmark the capabilities of each method. Crucially, our analysis of the open-system dynamics reveals that dissipation replaces the equilibrium phase diagram with one governed by non-equilibrium steady states (NESS), where the system evolves to either a stable solution or a persistent limit cycle.As such, this thesis is intended to serve a dual purpose: it is both a rigorous physical study of quantum phase transitions in a condensed matter system and a selfcontained pedagogical guide for researchers beginning to explore the phenomenology of quantum phase transitions through numerical schemes.