Dilute Bose-Einstein Condensate Research Articles

Superfluid Fraction in an Interacting Spatially Modulated Bose-Einstein Condensate.

At zero temperature, a Galilean-invariant Bose fluid is expected to be fully superfluid. Here we investigate theoretically and experimentally the quenching of the superfluid density of a dilute Bose-Einstein condensate due to the breaking of translational (and thus Galilean) invariance by an external 1D periodic potential. Both Leggett's bound fixed by the knowledge of the total density and the anisotropy of the sound velocity provide a consistent determination of the superfluid fraction. The use of a large-period lattice emphasizes the important role of two-body interactions on superfluidity.

Complex Langevin approach to interacting Bose gases

Quantitative numerical analyses of interacting dilute Bose-Einstein condensates are most often based on semiclassical approximations. Since the complex-valued field-theoretic action of the Bose gas does not offer itself to standard Monte Carlo techniques, simulations beyond their scope almost exclusively rely on quantum-mechanical techniques, which are inherently limited to small particle numbers. Here we explore an alternative approach based on a Langevin-type sampling in an extended state space, known as complex Langevin (CL) algorithm. While the use of the CL technique has a long-standing history in high-energy physics, in particular in the simulation of QCD at finite baryon density, applications to ultracold atoms are still in their infancy. Here we examine the applicability of the CL approach for a one- and two-component, three-dimensional non-relativistic Bose gas in thermal equilibrium, below and above the Bose-Einstein phase transition. By comparison with analytic descriptions at the Gaussian level, including Bogoliubov and Hartree-Fock theory, we find that the method allows computing reliably and efficiently observables in the regime of experimentally accessible parameters. Close to the transition, quantum corrections lead to a shift of the critical temperature which we reproduce within the numerical range known in the literature. With this work, we aim to provide a first test of CL as a potential out-of-the-box tool for the simulation of experimentally realistic situations, including trapping geometries and multicomponent/-species models.

Topologically protected vortex knots and links

In 1869, Lord Kelvin found that the way vortices are knotted and linked in an ideal fluid remains unchanged in evolution, and consequently hypothesized atoms to be knotted vortices in a ubiquitous ether, different knotting types corresponding to different types of atoms. Even though Kelvin’s atomic theory turned out incorrect, it inspired several important developments, such as the mathematical theory of knots and the investigation of knotted structures that naturally arise in physics. However, in previous studies, knotted and linked structures have been found to untie via local cut-and-paste events referred to as reconnections. Here, in contrast, we construct knots and links of non-Abelian vortices that are topologically protected in the sense that they cannot be dissolved employing local reconnections and strand crossings. Importantly, the topologically protected links are supported by a variety of physical systems such as dilute Bose-Einstein condensates and liquid crystals. We also propose a classification scheme for topological vortex links, in which two structures are considered equivalent if they differ from each other by a sequence of topologically allowed reconnections and strand crossings, in addition to the typical continuous transformations. Interestingly, this scheme produces a remarkably simple classification.

Resolving the puzzle of sound propagation in a dilute Bose–Einstein condensate

A unified model of a dilute Bose–Einstein condensate is proposed, combining the logarithmic and Gross–Pitaevskii (GP) nonlinear terms in a wave equation, where the GP term describes two-body interactions, as suggested by the standard perturbation theory; while the logarithmic term is essentially nonperturbative, and takes into account quantum vacuum effects. The model is shown to have excellent agreement with sound propagation data in the condensate of cold sodium atoms known since the now classic works by Andrews and collaborators. The data also allowed us to place constraints on two of the unified model’s parameters, which describe the strengths of the logarithmic and GP terms. Additionally, we suggest an experiment constraining the value of the third parameter (the characteristic density scale of the logarithmic part of the model), using the conjectured attraction–repulsion transition of many-body interaction inside the condensate.

Dropleton-soliton crossover mediated via trap modulation

We report a droplet to a soliton crossover by tuning the external confinement potential in a dilute Bose-Einstein condensate by numerically solving the modified Gross-Pitaevskii equation. The testimony of such a crossover is presented via studying the fractional density of the condensate which smoothly migrates from being a flat-head curve at weak confinement to a bright soliton at strong confinement. Such a transition occurs across a region of the potential whose strength varies over an order of magnitude and thus should be fit to be termed as a crossover. We supplement our studies via exploring the size of the bound pairs and the ramifications of the particle density therein. Eventually, all of these aid us in arriving at a phase diagram in a space defined by the trap strength and the particle number that shows the formation of two phases consisting of droplets and solitons, along with a regime of coexistence of these two.

Improved optical standing-wave beam splitters for dilute Bose–Einstein condensates

Bose–Einstein condensate (BEC)-based atom interferometry exploits low temperatures and long coherence lengths to facilitate high-precision measurements. Progress in atom interferometry promises improvements in navigational devices like gyroscopes and accelerometers, as well as applications in fundamental physics such as accurate determination of physical constants. Previous work demonstrates that beam splitters and mirrors for coherent manipulation of dilute BEC momentum in atom interferometers can be implemented with sequences of non-resonant standing-wave light pulses. While previous work focuses on the optimization of the optical pulses’ amplitude and duration to produce high-order momentum states with high fidelity, we explore how varying the shape of the optical pulses affects optimal beam-splitter performance, as well as the effect of pulse shape on the sensitivity of optimized parameters in achieving high fidelity in high-momentum states. In simulations of two-pulse beam splitters utilizing optimized square, triangle, and sinc-squared pulse shapes applied to dilute BECs, we, in some cases, reduce parameter sensitivity by an order of magnitude while maintaining fidelity.

Bogoliubov–de Gennes theory of the snake instability of gray solitons in higher dimensions

Gray solitons are a one-parameter family of solutions to the one-dimensional non-linear Schr\"odinger equation (NLSE) with positive cubic nonlinearity, as found in repulsively interacting dilute Bose-Einstein condensates or electromagnetic waves in the visible spectrum in waveguides described by Gross-Pitaevskii mean field theory. In two dimensions these solutions to the NLSE appear as a line or plane of depressed condensate density or light intensity, but numerical solutions show that this line is dynamically unstable to `snaking': the initially straight line of density or intensity minimum undulates with exponentially growing amplitude. To assist future studies of quantum mechanical instability beyond mean field theory, we here pursue an approximate analytical description of the snake instability within Bogoliubov-de Gennes perturbation theory. Within this linear approximation the two-dimensional result applies trivially to three dimensions as well, describing buckling modes of the low-density plane. We extend the analytical results of Kuznetsov and Turitsyn [Sov. Phys. JETP \textbf{67}, 1583 (1988)] to shorter wavelengths of the `snake' modulation and show to what extent the snake mode can be described accurately as a parametric instability, in which the position and grayness parameter of the initial soliton simply become dependent on the transverse dimension(s). We find that the parametric picture remains accurate up to second order in the snaking wave number, if the snaking soliton is also dressed by an outward-propagating sound wave, but that beyond second order in the snaking wave number the parametric description breaks down.

Mobile impurity in a Bose-Einstein condensate and the orthogonality catastrophe

We analyze the properties of an impurity in a dilute Bose-Einstein condensate (BEC). First the quasiparticle residue of a static impurity in an ideal BEC is shown to vanish with increasing particle number as a stretched exponential, leading to a bosonic orthogonality catastrophe. Then we introduce a variational ansatz, which recovers this exact result and describes the macroscopic dressing of the impurity including its back-action onto the BEC as well as boson-boson repulsion beyond the Bogoliubov approximation. This ansatz predicts that the orthogonality catastrophe also occurs for mobile impurities, whenever the BEC becomes ideal. Finally, we show that our ansatz agrees well with experimental results.

Disordered Bose–Einstein condensate in hard walls trap

We discuss the effects of quenched disorder in a dilute Bose–Einstein condensate confined in a hard walls trap. Starting from the disordered Gross–Pitaevskii functional, we obtain a representation for the quenched free energy as a series of integer moments of the partition function. Positive and negative disorder-dependent effective coupling constants appear in the integer moments. Going beyond the mean-field approximation, we compute the static two-point correlation functions at first-order in the positive effective coupling constants. We obtain the combined contributions of effects due to boundary conditions and disorder in this weakly disordered condensate. The ground state renormalized density profile of the condensate is presented. We also discuss the appearance of metastable and true ground states for strong disorder, when the effective coupling constants become negative.

Three-dimensional splitting dynamics of giant vortices in Bose-Einstein condensates

We study the splitting dynamics of giant vortices in dilute Bose-Einstein condensates by numerically integrating the three-dimensional Gross-Pitaevskii equation in time. By taking advantage of tetrahedral tiling in the spatial discretization, we decrease the error and increase the reliability of the numerical method. An extensive survey of vortex splitting symmetries is presented for different aspect ratios of the harmonic trapping potential. The symmetries of the splitting patterns observed in the simulated dynamics are found to be in good agreement with predictions obtained by solving the dominant dynamical instabilities from the corresponding Bogoliubov equations. Furthermore, we observe intertwining of the split vortices in prolate condensates and a split-and-revival phenomenon in a spherical condensate.

Quasi-steady radiation of sound from turbulent sonic ergoregions

Sonic Hawking radiation has recently been observed in dilute Bose–Einstein condensates (BECs), but it remains an open question whether this landmark achievement of atomic physics can lead to new insights into the effects on Hawking radiation of nonlinear back-reaction and new short-distance physics, as was originally hoped by Unruh when he introduced the sonic analogy. Furthermore, studies of sonic analog black holes have until now concentrated on (1+1)-dimensional scenarios, but Unruh’s sonic analogy for curved spacetime is only valid in more than one spatial dimension. We therefore model the evolution of a (2+1)-dimensional sonic black hole in a dilute BEC, over a long enough time to let the initial Corley–Jacobson instabilities saturate in vortex production and give way to a long-lived quasi-stationary state. In this quasi-equilibrium state we find the initial laminar ergoregion replaced by a turbulent zone that steadily radiates sound, but with a non-thermal power spectrum.

Gravitationally bound Bose condensates with rotation

We develop a self-consistent, Gravitoelectromagnetic (GEM) formulation of a slowly rotating, self-gravitating and dilute Bose-Einstein condensate (BEC), intended for astrophysical applications in the context of dark matter halos. GEM self-consistently incorporates the effects of frame dragging to lowest order in $v/c$ via the Gravitomagnetic field. BEC dark matter has attracted attention as an alternative to Cold dark matter (CDM) and Warm dark matter (WDM) for some time now. The BEC is described by the Gross-Pitaevskii-Poisson (GPP) equation with an arbitrary potential allowing for either attractive or repulsive interactions. Owing to the difficulty in obtaining exact solutions to the GEM equations of motion without drastic approximations, we employ the variational method to examine the conditions under which rotating condensates, stable against gravitational collapse, may form in models with attractive and repulsive quartic interactions. We also describe the approximate dynamics of an imploding and rotating condensate by employing a collective coordinate description in terms of the condensate radius.

Vortices in Bose-Einstein condensates with PT -symmetric gain and loss

We investigate vortex excitations in dilute Bose-Einstein condensates in the presence of complex $\mathcal{PT}$-symmetric potentials. These complex potentials are used to describe a balanced gain and loss of particles and allow an easier calculation of stationary states in open systems than in a full dynamical calculation including the whole environment. We examine the conditions under which stationary vortex states can exist and consider transitions from vortex to nonvortex states. In addition, we study the influences of $\mathcal{PT}$ symmetry on the dynamics of nonstationary vortex states placed at off-center positions.

Quantum-Fluctuation-Driven Crossover from a Dilute Bose-Einstein Condensate to a Macrodroplet in a Dipolar Quantum Fluid

In a joint experimental and theoretical effort, we report on the formation of a macro-droplet state in an ultracold bosonic gas of erbium atoms with strong dipolar interactions. By precise tuning of the s-wave scattering length below the so-called dipolar length, we observe a smooth crossover of the ground state from a dilute Bose-Einstein condensate (BEC) to a dense macro-droplet state of more than $10^4$ atoms. Based on the study of collective excitations and loss features, we quantitative prove that quantum fluctuations stabilize the ultracold gas far beyond the instability threshold imposed by mean-field interactions. Finally, we perform expansion measurements, showing the evolution of the normal BEC towards a three-dimensional self-bound state and show that the interplay between quantum stabilization and three-body losses gives rise to a minimal expansion velocity at a finite scattering length.

Hybrid OpenMP/MPI programs for solving the time-dependent Gross–Pitaevskii equation in a fully anisotropic trap

We present hybrid OpenMP/MPI (Open Multi-Processing/Message Passing Interface) parallelized versions of earlier published C programs (D. Vudragovic et al., Comput. Phys. Commun. 183, 2021 (2012)) for calculating both stationary and non-stationary solutions of the time-dependent Gross-Pitaevskii (GP) equation in three spatial dimensions. The GP equation describes the properties of dilute Bose-Einstein condensates at ultra-cold temperatures. Hybrid versions of programs use the same algorithms as the C ones, involving real- and imaginary-time propagation based on a split-step Crank-Nicolson method, but consider only a fully-anisotropic three-dimensional GP equation, where algorithmic complexity for large grid sizes necessitates parallelization in order to reduce execution time and/or memory requirements per node. Since distributed memory approach is required to address the latter, we combine MPI programing paradigm with existing OpenMP codes, thus creating fully flexible parallelism within a combined distributed/shared memory model, suitable for different modern computer architectures. The two presented C/OpenMP/MPI programs for real- and imaginary-time propagation are optimized and accompanied by a customizable makefile. We present typical scalability results for the provided OpenMP/MPI codes and demonstrate almost linear speedup until inter-process communication time starts to dominate over calculation time per iteration. Such a scalability study is necessary for large grid sizes in order to determine optimal number of MPI nodes and OpenMP threads per node.

Propagation of a quantum fluid of light in a cavityless nonlinear optical medium: General theory and response to quantum quenches

Making use of a generalized quantum theory of paraxial light propagation where the radiation-axis and the temporal coordinates play exchanged roles, we discuss the potential of bulk nonlinear optical media in cavityless configurations for quantum statistical mechanics studies of the conservative many-body dynamics of a gas of interacting photons. To illustrate the general features of this point of view, we investigate the response of the fluid of light to the quantum quenches in the photon-photon interaction constant experienced at the front and the back faces of a finite slab of weakly nonlinear material. Extending the standard Bogoliubov theory of dilute Bose-Einstein condensates, peculiar features are predicted for the statistical properties of the light emerging from the nonlinear medium.

The viscosity of dilute Bose–Einstein condensates

The viscosity of a dilute Bose–Einstein condensate is obtained for a range of temperatures and densities. The kinetic equation used to derive the viscosity is based on Bogoliubov mean field theory and is a Boltzmann-like equation for the Bogoliubov excitations in the condensate. The viscosity can be assigned a numerical value for any gas for which the mass and scattering length are known. The viscosity decreases slowly as the temperature is decreased, but at low temperature (around 0.01 Tc), begins to show a much more rapid decrease with decreasing temperature. This interesting behavior of the viscosity is related to the behavior of the eigenvalues of the bogolon collision operator at these low temperatures.

Symmetry breaking of localized discrete matter waves induced by spin–orbit coupling

We study localized nonlinear excitations of a dilute Bose–Einstein condensate (BEC) with spin–orbit coupling in a deep optical lattice (OL). For this we introduce a tight-binding model that includes the spin–orbit coupling (SOC) at the discrete level in the form of a generalized discrete nonlinear Schrödinger equation. Existence and stability of discrete solitons of different symmetry types is demonstrated. Quite interestingly, we find three distinctive regions in which discrete solitons undergo spontaneously symmetry breaking, passing from on-site to inter-site and to asymmetric, simply by varying the interatomic interactions. Existence ranges of discrete solitons with inter-site symmetry depend on SOC and shrink to zero as the SOC parameter is increased. Asymmetric discrete solitons appear as novel excitations specific of the SOC. Possible experimental implementation of these results is briefly discussed.

Optomechanical signature of a frictionless flow of superfluid light

We propose an experimental setup that should make it possible to reveal the frictionless flow of a superfluid of light from the suppression of the drag force that it exerts onto a material obstacle. In the paraxial-propagation geometry considered here, the photon-fluid dynamics is described by a wave equation analogous to the Gross-Pitaevskii equation of dilute Bose-Einstein condensates and the obstacle consists in a solid dielectric slab immersed into a nonlinear optical liquid. By means of an ab initio calculation of the electromagnetic force experienced by the obstacle, we anticipate that superfluidity is detectable in state-of-the-art experiments from the disappearance of the optomechanical deformation of the obstacle.

Association of atoms into universal dimers using an oscillating magnetic field.

In a system of ultracold atoms near a Feshbach resonance, pairs of atoms can be associated into universal dimers by an oscillating magnetic field with a frequency near that determined by the dimer binding energy. We present a simple expression for the transition rate that takes into account many-body effects through a transition matrix element of the contact. In a thermal gas, the width of the peak in the transition rate as a function of the frequency is determined by the temperature. In a dilute Bose-Einstein condensate of atoms, the width is determined by the inelastic scattering rates of a dimer with zero-energy atoms. Near an atom-dimer resonance, there is a dramatic increase in the width from inelastic atom-dimer scattering and from atom-atom-dimer recombination. The recombination contribution provides a signature for universal tetramers that are Efimov states consisting of two atoms and a dimer.