We consider a time-dependent extension of a perturbative mean-field approach to the homogeneous dirty boson problem by considering how switching on and off a weak disorder potential affects the stationary state of an initially {equilibrated} Bose-Einstein condensate by the emergence of a disorder-induced condensate deformation. We find that in the switch on scenario the stationary condensate deformation turns out to be a sum of an equilibrium part{, that actually corresponds to adiabatic switching on the disorder,} and a dynamically-induced part, where the latter depends on the particular driving protocol. If the disorder is switched off afterwards, the resulting condensate deformation acquires an additional dynamically-induced part in the long-time limit, while the equilibrium part vanishes. {We also present an appropriate generalization to inhomogeneous trapped condensates.} Our results demonstrate that the condensate deformation represents an indicator of the generically non-equilibrium nature of steady states of a Bose gas in a temporally controlled weak disorder.
We consider a time-dependent extension of a perturbative mean-field approach to the homogeneous dirty boson problem by considering how switching on and off a weak disorder potential affects the stationary state of an initially {equilibrated} Bose-Einstein condensate by the emergence of a disorder-induced condensate deformation. We find that in the switch on scenario the stationary condensate deformation turns out to be a sum of an equilibrium part{, that actually corresponds to adiabatic switching on the disorder,} and a dynamically-induced part, where the latter depends on the particular driving protocol. If the disorder is switched off afterwards, the resulting condensate deformation acquires an additional dynamically-induced part in the long-time limit, while the equilibrium part vanishes. {We also present an appropriate generalization to inhomogeneous trapped condensates.} Our results demonstrate that the condensate deformation represents an indicator of the generically non-equilibrium nature of steady states of a Bose gas in a temporally controlled weak disorder.
We theoretically study how the superfluid and condensate deformation of a weakly interacting ultracold Bose gas evolve during the ramp-up of an external weak disorder potential. Both resulting deformations turn out to consist of two distinct contributions, namely a reversible equilibrium one, already predicted by Huang and Meng in 1992, and a nonequilibrium dynamical one, whose magnitude depends on the details of the ramping protocol. For the specific case of the exponential ramp-up protocol, we are able to derive analytical time-dependent expressions for the above quantities. After a sufficiently long time, a steady state emerges that is generically out of equilibrium. We take the first step in investigating its properties by studying its relaxation dynamics. In addition, we analyze the two-time correlation function and elucidate its relation to the equilibrium and the dynamical part of the condensate deformation.
Abstract In an equilibrium three-dimensional (3D) disordered condensate, it is well established that disorder can generate an amount of normal fluid ρ n equaling to of ρ ex , where ρ ex is a sum of interaction-induced quantum depletion and disorder-induced condensate deformation. The concept that the superfluid is more volatile to the existence of disorder than the condensate is crucial to the understanding of the Bose glass phase. In this work, we show that, by bringing a weakly disordered 3D condensate to non-equilibrium regime via a quantum quench in the interaction, disorder can destroy superfluid significantly more, leading to a steady state of Hamiltonian H f in which the ρ n far exceeds of the ρ ex . This suggests the possibility of engineering Bose glass in the dynamic regime. Here, we refer to the dynamical Bose glass as the case where in the steady state of quenched condensate , the superfluid density goes to zero while the condensate density remains finite. As both the ρ n and ρ ex are measurable quantities, our results allow an experimental demonstration of the dramatized interplay between the disorder and interaction in the non-equilibrium scenario.
We study the non-equilibrium evolution of binary Bose-Einstein condensates in the presence of a weak random potential with Gaussian correlation function using the time-dependent perturbation theory. We apply this theory to construct a closed set of equations that highlight the role of the spectacular interplay between the disorder and the interspecies interactions in the time evolution of the density induced by disorder in each component. It is found that this latter increases with time favoring localization of both species. The time scale at which the theory remains valid depends on the respective system parameters. We show analytically and numerically that such a system supports a steady state that periodically changing during its time propagation. The obtained dynamical corrections indicate that disorder may transform the system into a stationary out-of-equilibrium states. Understanding this time evolution is pivotal for the realization of Floquet condensates.
The time-dependent non-equilibrium dynamics of a Bose–Einstein condensate (BEC) typically generates incoherent excitations out of the condensate due to the finite frequencies present in the time evolution. We present a detailed derivation of a general non-equilibrium Greenʼs function technique that describes the coupled time evolution of an interacting BEC and its single-particle excitations in a trap, based on an expansion in terms of the exact eigenstates of the trap potential. We analyze the dynamics of a Bose system in a small double-well potential with initially all particles in the condensate. When the trap frequency is larger than the Josephson frequency, , the dynamics changes at a characteristic time, , abruptly from the slow Josephson oscillations of the BEC to fast Rabi oscillations driven by quasiparticle excitations in the trap. For times , the Josephson oscillations are undamped, in agreement with the experiments. We analyze the physical origin of the finite scale as well as its dependence on the trap parameter Δ .
Weakly interacting Bose gases in a disorder environment have long been a challenging topic in the field of solid-state physics due to the intriguing interplay between superfluidity and localization. In this thesis we investigate theoretically the equilibrium properties of one- and three-dimensional harmonically trapped dirty Bose-Einstein condensates (BECs). In particular we focus on the decomposition of the total particle density into three components, namely, the condensate density, the thermal density, and the density of fragmented local Bose-Einstein condensates within the respective minima of the random potential landscape. In order to determine the different density components we develop a Hartree-Fock mean-field theory for a weakly interacting trapped BEC in a random environment, then specialize it to contact interaction and delta-correlated disorder. This non-perturbative theory is worked out on the basis of the replica method, which represents a well- established tool to deal with disorder problems. With this we calculate the corresponding free energy at finite temperature and derive from it the underlying self-consistency equations for the respective density components. Then, as a first step, we investigate a quasi one-dimensional Bose-Einstein condensed gas in a harmonic trapping potential with an additional delta- correlated disorder potential at zero temperature. To this end we solve the self-consistency equations within the Thomas-Fermi approximation and find the emergence of a Bose-glass region, where the condensate vanishes. We corroborate this analysis by an elaborate numerical treatment, where the corresponding one-dimensional time-independent Gross-Pitaevskii equation is numerically solved and disorder-ensemble averages of the condensate wave function are performed. The variance of the condensate wave function quantifies the number of bosons which condense in the local minima of the random potential. For weak disorder these mini-condensates turn out to occur preferentially at the border of the condensate, while for intermediate disorder strength this happens in the trap center. Additionally we use a variational ansatz in order to describe analytically the numerically observed redistribution of the fragmented mini-condensates with increasing disorder strength. In close analogy to the one-dimensional case we then treat the three-dimensional dirty Bose gas in an isotropic harmonic trap and a delta- correlated disorder potential at zero temperature using the corresponding self-consistency equations obtained via the Hartree-Fock mean-field theory within the Thomas-Fermi approximation. Additionally we use a variational ansatz, whose results turn out to coincide qualitatively with those obtained from the Thomas-Fermi approximation. Inparticular, a first order quantum phase transition from the superfluid phase, where the condensate density contributes to the total density, to the Bose-glass phase, where all particles are in the mini-condensates, is detected at a critical disorder strength, which agrees with findings in the literature. Furthermore, in a generaltriaxial harmonic trap, we investigate the geometric effect of different trap aspect ratios on the respective properties of the dirty BEC system. Finally, we consider the three- dimensional dirty BEC at finite temperature. This allows us to study the impact of both temperature and disorder fluctuations on the respective components of the density as well as their Thomas-Fermi radii. In particular, we find that the superfluid region, the Bose-glass region, and the thermal region coexist. Furthermore, depending on the system parameters, three phase transitions are detected, namely, one from the superfluid to the Bose-glass phase, one from the Boseglass to the thermal phase, and finally one directly from the superfluid to the thermal phase. All these results could be particularly useful for a quantitative analysis of world-wide on-going experiments with dirty bosons in quasi one-dimensional and three-dimensional harmonic traps.
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.