In this Rapid Communication we study the transient dynamics of a Bose superfluid subsequent to an interaction quench. Essential for equilibration is a source of dissipation which we include following the approach of Caldeira and Leggett. Here we solve the equations of motion exactly by integrating out an environmental bath. We thereby derive precisely the time dependent density correlation functions with the appropriate analytic and asymptotic properties. The resulting structure factor exhibits the expected damping and thereby differs from that of strict Bogoliubov theory. These damped sound modes, which reflect the physics beyond mean-field approaches, are characterized and the structure factors are found to compare favorably with experiment.
In this Rapid Communication we study the transient dynamics of a Bose superfluid subsequent to an interaction quench. Essential for equilibration is a source of dissipation which we include following the approach of Caldeira and Leggett. Here we solve the equations of motion exactly by integrating out an environmental bath. We thereby derive precisely the time dependent density correlation functions with the appropriate analytic and asymptotic properties. The resulting structure factor exhibits the expected damping and thereby differs from that of strict Bogoliubov theory. These damped sound modes, which reflect the physics beyond mean-field approaches, are characterized and the structure factors are found to compare favorably with experiment.
We address the physics of equilibration in ultracold atomic gases following a quench of the interaction parameter. Our work is based on a bath model which generates damping of the bosonic excitations. We illustrate this dissipative behavior through the momentum distribution of the excitations nk, observing that larger k modes have shorter relaxation times $\ensuremath{\tau}(k)$; they will equilibrate faster, as has been claimed in recent experimental work. We identify three time regimes. At short times nk exhibits oscillations; these are damped out at intermediate times where the system appears to be in a false or slowly converging equilibrium. Finally, at longer times, full equilibration occurs. This false equilibrium is, importantly, associated with the k dependence in $\ensuremath{\tau}(k)$ and has implications for experiment.
The exact description of the time evolution of open correlated quantum systems remains one of the major challenges of condensed matter theory, especially for asymptotic long times where most numerical methods fail. Here, the post-quench dynamics of the N-component Bose-Anderson impurity model is studied in the $N\ensuremath{\rightarrow}\ensuremath{\infty}$ limit. The equilibrium phase diagram is similar to that of the Bose-Hubbard model in that it contains local versions of Mott and Bose phases. Using a numerically exact procedure, we are able to study the real-time evolution including asymptotic long-time regimes. The formation of long-lived transient phases is observed for quench paths crossing foreign phases. For quenches inside the local Bose condensed phase, a dynamical phase transition is reported that separates the evolution towards a new equilibrium state and a regime characterized at large times by a persistent phase rotation of the order parameter. We explain how such nondecaying modes can exist in the presence of a dissipative bath. We discuss the extension of our results to the experimental relevant finite-N case and their implication for the existence of nondecaying modes in generic quantum systems in the presence of a bath.
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.
Yin and Radzihovsky [X. Yin and L. Radzihovsky, Phys. Rev. A 88, 063611 (2013)] recently developed a self-consistent extension of a Bogoliubov theory, in which the condensate number density nc is treated as a mean field that changes with time, in order to analyze a JILA experiment by Makotyn et al. [P. Makotyn et al., Nat. Phys. 10, 116 (2014)] on a 85Rb Bose gas following a deep quench to a large scattering length. We apply this theory to construct a closed set of equations that highlight the role of ${\stackrel{\ifmmode \dot{}\else \.{}\fi{}}{n}}_{c}$, which is to induce an effective interaction between quasiparticles. We show analytically that such a system supports a steady state characterized by a constant condensate density and a steady but periodically changing momentum distribution, whose time average is described exactly by the generalized Gibbs ensemble. We discuss how the ${\stackrel{\ifmmode \dot{}\else \.{}\fi{}}{n}}_{c}$-induced effective interaction, which cannot be ignored on the grounds of the adiabatic approximation for modes near the gapless Goldstone mode, can significantly affect condensate populations and Tan's contact for a Bose gas that has undergone a deep quench.
We study the dynamics of a two-mode Bose-Einstein condensate in the vicinity of a mean-field dynamical instability. Convergence to mean-field theory (MFT), with increasing total number of particles N, is shown to be logarithmically slow. Using a density-matrix formalism rather than the conventional wave-function methods, we derive an improved set of equations of motion for the mean-field plus the fluctuations, which goes beyond MFT and provides accurate predictions for the leading quantum corrections and the quantum break time. We show that the leading quantum corrections appear as decoherence of the reduced single-particle quantum state; we also compare this phenomenon to the effects of thermal noise. Using the rapid dephasing near an instability, we propose a method for the direct measurement of scattering lengths.
We investigate the superfluid properties of a Bose–Einstein condensate (BEC) trapped in a one-dimensional periodic potential. We study, both analytically (in the tight binding limit) and numerically, the Bloch chemical potential, the Bloch energy and the Bogoliubov dispersion relation, and we introduce two different, density dependent, effective masses and group velocities. The Bogoliubov spectrum predicts the existence of sound waves, and the arising of energetic and dynamical instabilities at critical values of the BEC quasi-momentum which dramatically affect its coherence properties. We investigate the dependence of the dipole and Bloch oscillation frequencies in terms of an effective mass averaged over the density of the condensate. We illustrate our results with several animations obtained solving numerically the time-dependent Gross–Pitaevskii equation.