Bright Soliton to Quantum Droplet Transition in a Mixture of Bose-Einstein Condensates

Measurement of the Excitation Spectrum of a Dipolar Gas in the Macrodroplet Regime

The excitation spectrum of a cigar-shaped strongly dipolar quantum gas at the crossover from a Bose-Einstein condensate to a trapped macrodroplet is predicted to exhibit peculiar features-a strong upward shift of low momentum excitation energies together with a strong multiband response for high momenta. By performing Bragg spectroscopy over a wide range of momenta, we observe both key elements and also confirm the predicted stiffening of excitation modes when approaching the macrodroplet regime. Our measurements are in good agreement with numerical calculations taking into account finite size effects.

Discrete and Semi-Discrete Multidimensional Solitons and Vortices: Established Results and Novel Findings

This article presents a concise survey of basic discrete and semi-discrete nonlinear models, which produce two- and three-dimensional (2D and 3D) solitons, and a summary of the main theoretical and experimental results obtained for such solitons. The models are based on the discrete nonlinear Schrödinger (DNLS) equations and their generalizations, such as a system of discrete Gross-Pitaevskii (GP) equations with the Lee-Huang-Yang corrections, the 2D Salerno model (SM), DNLS equations with long-range dipole-dipole and quadrupole-quadrupole interactions, a system of coupled discrete equations for the second-harmonic generation with the quadratic (χ(2)) nonlinearity, a 2D DNLS equation with a superlattice modulation opening mini-gaps, a discretized NLS equation with rotation, a DNLS coupler and its PT-symmetric version, a system of DNLS equations for the spin-orbit-coupled (SOC) binary Bose-Einstein condensate, and others. The article presents a review of the basic species of multidimensional discrete modes, including fundamental (zero-vorticity) and vortex solitons, their bound states, gap solitons populating mini-gaps, symmetric and asymmetric solitons in the conservative and PT-symmetric couplers, cuspons in the 2D SM, discrete SOC solitons of the semi-vortex and mixed-mode types, 3D discrete skyrmions, and some others.

Quantum Criticality of Liquid-Gas Transition in a Binary Bose Mixture

Quantum liquid, in the form of a self-bound droplet, is stabilized by a subtle balance between the mean-field contribution and quantum fluctuations. While a liquid-gas transition is expected when such a balance is broken, it remains elusive whether liquid-gas critical points exist in the quantum regime. Here, we study the quantum criticality in a binary Bose mixture undergoing the liquid-gas transition. We show that, beyond a narrow stability window of the self-bound liquid, a liquid-gas coexistence persists, which eventually transits into a homogeneous mixture. Importantly, we identify two distinct critical points where the liquid-gas coexistence terminates. These critical points are characterized by rich critical behaviors in their vicinity, including divergent susceptibility, unique phonon-mode softening, and enhanced density correlations. The liquid-gas transition and the critical points can be readily explored in ultracold atoms confined to a box potential. Our work highlights the thermodynamic approach as a powerful tool in revealing the quantum liquid-gas criticality, and paves the way for further studies of critical phenomena in quantum liquids.

Formations and dynamics of two-dimensional spinning asymmetric quantum droplets controlled by a PT-symmetric potential

In this paper, vortex solitons are produced for a variety of 2D spinning quantum droplets (QDs) in a PT-symmetric potential, modeled by the amended Gross-Pitaevskii equation with Lee-Huang-Yang corrections. In particular, exact QD states are obtained under certain parameter constraints, providing a guide to finding the respective generic family. In a parameter region of the unbroken PT symmetry, different families of QDs originating from the linear modes are obtained in the form of multipolar and vortex droplets at low and high values of the norm, respectively, and their stability is investigated. In the spinning regime, QDs become asymmetric above a critical rotation frequency, most of them being stable. The effect of the PT-symmetric potential on the spinning and nonspinning QDs is explored by varying the strength of the gain-loss distribution. Generally, spinning QDs trapped in the PT-symmetric potential exhibit asymmetry due to the energy flow affected by the interplay of the gain-loss distribution and rotation. Finally, interactions between spinning or nonspinning QDs are explored, exhibiting elastic collisions under certain conditions.

Ultrawide Dark Solitons and Droplet-Soliton Coexistence in a Dipolar Bose Gas with Strong Contact Interactions

We look into dark solitons in a quasi-1D dipolar Bose gas and in a quantum droplet. We derive the analytical solitonic solution of a Gross-Pitaevskii-like equation accounting for beyond mean-field effects. The results show there is a certain critical value of the dipolar interactions, for which the width of a motionless soliton diverges. Moreover, there is a peculiar solution of the motionless soliton with a nonzero density minimum. We also present the energy spectrum of these solitons with an additional excitation subbranch appearing. Finally, we perform a series of numerical experiments revealing the coexistence of a dark soliton inside a quantum droplet.

Superexchange Liquefaction of Strongly Correlated Lattice Dipolar Bosons

We propose a mechanism for liquid formation in strongly correlated lattice systems. The mechanism is based on an interplay between long-range attraction and superexchange processes. As an example, we study dipolar bosons in one-dimensional optical lattices. We present a perturbative theory and validate it in comparison with full density-matrix renormalization group simulations for the energetic and structural properties of different phases of the system, i.e., self-bound Mott insulator, liquid, and gas. We analyze the nonequilibrium properties and calculate the dynamic structure factor. Its structure differs in compressible and insulating phases. In particular, the low-energy excitations in compressible phases are linear phonons. We extract the speed of sound and analyze its dependence on dipolar interaction and density. We show that it exhibits a nontrivial behavior owing to the breaking of Galilean invariance. We argue that an experimental detection of this previously unknown quantum liquid could provide a fingerprint of the superexchange process and open intriguing possibilities for investigating non-Galilean invariant liquids.

Reflectionless potentials and resonant scattering of flat-top and thin-top solitons

We identify a class of potentials for which the scattering of flat-top solitons and thin-top solitons of the nonlinear Schrödinger equation with dual nonlinearity can be reflectionless. The scattering is characterized by sharp resonances between regimes of full transmission and full quantum reflection. Perturbative expansion in terms of the magnitude of radiation losses leads to the general form of reflectionless potentials. Simulating the scattering of flat-top solitons and thin-top solitons confirms the reflectionless feature of these potentials.

Dipolar physics: a review of experiments with magnetic quantum gases

Since the achievement of quantum degeneracy in gases of chromium atoms in 2004, the experimental investigation of ultracold gases made of highly magnetic atoms has blossomed. The field has yielded the observation of many unprecedented phenomena, in particular those in which long-range and anisotropic dipole-dipole interactions (DDIs) play a crucial role. In this review, we aim to present the aspects of the magnetic quantum-gas platform that make it unique for exploring ultracold and quantum physics as well as to give a thorough overview of experimental achievements. Highly magnetic atoms distinguish themselves by the fact that their electronic ground-state configuration possesses a large electronic total angular momentum. This results in a large magnetic moment and a rich electronic transition spectrum. Such transitions are useful for cooling, trapping, and manipulating these atoms. The complex atomic structure and large dipolar moments of these atoms also lead to a dense spectrum of resonances in their two-body scattering behaviour. These resonances can be used to control the interatomic interactions and, in particular, the relative importance of contact over dipolar interactions. These features provide exquisite control knobs for exploring the few- and many-body physics of dipolar quantum gases. The study of dipolar effects in magnetic quantum gases has covered various few-body phenomena that are based on elastic and inelastic anisotropic scattering. Various many-body effects have also been demonstrated. These affect both the shape, stability, dynamics, and excitations of fully polarised repulsive Bose or Fermi gases. Beyond the mean-field instability, strong dipolar interactions competing with slightly weaker contact interactions between magnetic bosons yield new quantum-stabilised states, among which are self-bound droplets, droplet assemblies, and supersolids. Dipolar interactions also deeply affect the physics of atomic gases with an internal degree of freedom as these interactions intrinsically couple spin and atomic motion. Finally, long-range dipolar interactions can stabilise strongly correlated excited states of 1D gases and also impact the physics of lattice-confined systems, both at the spin-polarised level (Hubbard models with off-site interactions) and at the spinful level (XYZ models). In the present manuscript, we aim to provide an extensive overview of the various related experimental achievements up to the present.

Droplet to soliton crossover at negative temperature in presence of bi-periodic optical lattices

It is shown that the phenomenon of negative temperature essentially occurs in Bose-Einstein condensate due to the realization of the upper bound energy state utilizing a combination of expulsive harmonic oscillator and optical lattice potentials. We study the existence of quantum droplets at negative temperature and droplet-to-soliton crossover in the binary Bose-Einstein condensate mixture in the presence of bi-periodic optical lattices and expulsive-BOL confinements. Based on the beyond mean field approximation, we employ the extended Gross-Pitäevskii equation and calculate the exact analytical form of wavefunction solutions for BOL, expulsive-BOL confinements. An interesting transition of quantum droplets from positive to negative temperatures and the droplet-to-soliton crossover by modulating the disorder in BOL potential are illustrated. The affirmation of such crossover is performed by exploring the profile of atomic condensate density which smoothly transits from being a flat top density in optical lattice confinement to a bright soliton for BOL trap. Further, we confirm the crossover by exploring the energy per particle and the variation in the root mean square size of the condensate with respect to the potential depth of the BOL trap. Eventually, all of this aid us to construct a phase diagram in a space between the amplitude of BOL potential depth and particle number which reveals the formation of droplet and soliton phases. In expulsive-BOL confinement, it is seen that the impact of the expulsive trap is insignificant on atomic condensate density in the droplet phase and it becomes prominent in the soliton region. Further, the variation of total energy reveals that the amplitude of the expulsive oscillator strengthens the droplet phase and leads to an increase in the negative temperature of the considered system.

Tuning Long-Range Fermion-Mediated Interactions in Cold-Atom Quantum Simulators

Engineering long-range interactions in cold-atom quantum simulators can lead to exotic quantum many-body behavior. Fermionic atoms in ultracold atomic mixtures can act as mediators, giving rise to long-range Ruderman-Kittel-Kasuya-Yosida-type interactions characterized by the dimensionality and density of the fermionic gas. Here, we propose several tuning knobs, accessible in current experimental platforms, that allow one to further control the range and shape of the mediated interactions, extending the existing quantum simulation toolbox. In particular, we include an additional optical lattice for the fermionic mediator, as well as anisotropic traps to change its dimensionality in a continuous manner. This allows us to interpolate between power-law and exponential decays, introducing an effective cutoff for the interaction range, as well as to tune the relative interaction strengths at different distances. Finally, we show how our approach allows one to investigate frustrated regimes that were not previously accessible, where symmetry-protected topological phases as well as chiral spin liquids emerge.

Stability limits for modes held in alternating trapping-expulsive potentials

We elaborate a scheme of trapping-expulsion management (TEM), in the form of the quadratic potential periodically switching between confinement and expulsion, as a means of stabilization of two-dimensional dynamical states against the backdrop of the critical collapse driven by the cubic self-attraction with strength g. The TEM scheme may be implemented, as spatially or temporally periodic modulations, in optics or BEC, respectively. The consideration is carried out by dint of numerical simulations and variational approximation (VA). In terms of the VA, the dynamics amounts to a nonlinear Ermakov equation, which, in turn, is tantamount to a linear Mathieu equation. Stability boundaries are found as functions of g and parameters of the periodic modulation of the trapping potential. Below the usual collapse threshold, which is known, in the numerical form, as g<g_{c}^{(num)}≈5.85 (in the standard notation), the stability is limited by the onset of the parametric resonance. This stability limit, including the setup with the self-repulsive sign of the cubic term (g<0), is accurately predicted by the VA. At g>g_{c}^{(num)}, the collapse threshold is found with the help of full numerical simulations. The relative increase of g_{c} above g_{c}^{(num)} is ≈1.5%. It is a meaningful result, even if its size is small, because the collapse threshold is a universal constant which is difficult to change.

Ultradilute Quantum Droplets in the Presence of Higher-Order Quantum Fluctuations

We investigate the effects of higher-order quantum fluctuations on the bulk properties of self-bound droplets in three-, two- and one-dimensional binary Bose mixtures using the Hartree–Fock–Bogoliubov theory. We calculate higher-order corrections to the equation of state of the droplet at both zero and finite temperatures. We show that our results for the ground-state energy are in a good agreement with recent quantum Monte Carlo simulations in any dimension. Our study extends to the finite temperature case where it is found that thermal fluctuations may destabilize the droplet state and eventually destroy it. In two dimensions, we reveal that the droplet occurs at temperatures well below the Berezinskii–Kosterlitz–Thouless transition temperature.

Dynamics of quantum droplets in an external harmonic confinement

Recent theoretical and experimental results show that one-dimensional (1D) weakly interacting atomic Bose-Bose mixtures with repulsive interspecies mean field (MF) interaction are stabilized by attractive quadratic beyond-mean-field (BMF) effects into self-bound quantum droplet (QD) in free space. Here, we construct an exact analytical model to investigate the structure and dynamics of QDs in presence of external harmonic confinement by solving the 1D extended Gross–Pitäevskii equation (eGPE) with temporal variation of MF and BMF interactions. The model provides the analytical form of wavefunction, phase, MF and BMF nonlinearities. The generation of QDs and interesting droplet to soliton transition in presence of regular/expulsive parabolic traps by taking the comparable MF and BMF interactions are illustrated. We derive the phase diagram of the droplet-soliton phase transition between amplitude of MF, BMF interactions and harmonic oscillator frequency. The strength and form of oscillator frequency are identified as key parameter for tuning the compression, fragmentation and transport of droplets. Finally, the stability of the obtained solutions are confirmed from Vakhitov–Kolokolov (VK) criterion and are found stable.

Tunable Three-Body Interactions in Driven Two-Component Bose-Einstein Condensates

We propose and demonstrate the appearance of an effective attractive three-body interaction in coherently driven two-component Bose-Einstein condensates. It originates from the spinor degree of freedom that is affected by a two-body mean-field shift of the driven transition frequency. Importantly, its strength can be controlled with the Rabi-coupling strength and it does not come with additional losses. In the experiment, the three-body interactions are adjusted to play a predominant role in the equation of state of a cigar-shaped trapped condensate. This is confirmed through two striking observations: a downshift of the radial breathing mode frequency and the radial collapses for positive values of the dressed-state scattering length.

Interaction Control and Bright Solitons in Coherently Coupled Bose-Einstein Condensates

We demonstrate fast control of the interatomic interactions in a Bose-Einstein condensate by coherently coupling two atomic states with intra- and interstate scattering lengths of opposite signs. We measure the elastic and inelastic scattering properties of the system and find good agreement with a theoretical model describing the interactions between dressed states. In the attractive regime, we observe the formation of bright solitons formed by dressed-state atoms. Finally, we study the response of the system to an interaction quench from repulsive to attractive values, and observe how the resulting modulational instability develops into a bright soliton train.

Beyond-Mean-Field Effects in Rabi-Coupled Two-Component Bose-Einstein Condensate

We theoretically calculate and experimentally measure the beyond-mean-field (BMF) equation of state in a coherently coupled two-component Bose-Einstein condensate (BEC) in the regime where averaging of the interspecies and intraspecies coupling constants over the hyperfine composition of the single-particle dressed state predicts the exact cancellation of the two-body interaction. We show that with increasing the Rabi-coupling frequency Ω, the BMF energy density crosses over from the nonanalytic Lee-Huang-Yang scaling ∝n^{5/2} to an expansion in integer powers of density, where, in addition to a two-body BMF term ∝n^{2}sqrt[Ω], there emerges a repulsive three-body contribution ∝n^{3}/sqrt[Ω]. We experimentally evidence these two contributions, thanks to their different scaling with Ω, in the expansion of a Rabi-coupled two-component ^{39}K condensate in a waveguide. By studying the expansion with and without Rabi coupling, we reveal an important feature relevant for observing BMF effects and associated phenomena in mixtures with spin-asymmetric losses: Rabi coupling helps preserve the spin composition and thus prevents the system from drifting away from the point of the vanishing mean field.

Many-body and temperature effects in two-dimensional quantum droplets in Bose-Bose mixtures

We study the equilibrium properties of self-bound droplets in two-dimensional Bose mixtures employing the time-dependent Hartree–Fock–Bogoliubov theory. This theory allows one to understand both the many-body and temperature effects beyond the Lee–Huang–Yang description. We calculate higher-order corrections to the excitations, the sound velocity, and the energy of the droplet. Our results for the ground-state energy are compared with the diffusion Monte Carlo data and good agreement is found. The behavior of the depletion and anomalous density of the droplet is also discussed. At finite temperature, we show that the droplet emerges at temperatures well below the Berezinskii–Kosterlitz–Thouless transition temperature. The critical temperature strongly depends on the interspecies interactions. Our study is extended to the finite size droplet by numerically solving the generalized finite-temperature Gross-Pitaevskii equation which is obtained self-consistently from our formalism in the framework of the local density approximation.

Borromean Droplet in Three-Component Ultracold Bose Gases

We investigate droplet formation in three-component ultracold bosons. In particular, we identify the formation of a Borromean droplet, where only the ternary bosons can form a self-bound droplet while any binary subsystems cannot, as the first example of Borromean binding due to a collective many-body effect. Its formation is facilitated by an additional attractive force induced by the density fluctuation of a third component, which enlarges the mean-field collapse region in comparison to the binary case and renders the formation of a Borromean droplet after incorporating the repulsive force from quantum fluctuations. Outside the Borromean regime, we demonstrate an interesting phenomenon of droplet phase separation due to the competition between ternary and binary droplets. We further show that the transition between different droplets and gas phase can be conveniently tuned by boson numbers and interaction strengths. The study reveals the rich physics of a quantum droplet in three-component boson mixtures and sheds light on the more intriguing many-body bound state formed in multicomponent systems.

Observation of Scale Invariance in Two-Dimensional Matter-Wave Townes Solitons

We report near-deterministic generation of two-dimensional (2D) matter-wave Townes solitons and a precision test on scale invariance in attractive 2D Bose gases. We induce a shape-controlled modulational instability in an elongated 2D matter wave to create an array of isolated solitary waves of various sizes and peak densities. We confirm scale invariance by observing the collapse of solitary-wave density profiles onto a single curve in a dimensionless coordinate rescaled according to their peak densities and observe that the scale-invariant profiles measured at different coupling constants g can further collapse onto the universal profile of Townes solitons. The reported scaling behavior is tested with a nearly 60-fold difference in soliton interaction energies and allows us to discuss the impact of a non-negligible magnetic dipole-dipole interaction (MDDI) on 2D scale invariance. We confirm that the effect of MDDI in our alkali cesium quasi-2D samples effectively conforms to the same scaling law governed by a contact interaction to well within our experiment uncertainty.

Realization of a Townes Soliton in a Two-Component Planar Bose Gas

Most experimental observations of solitons are limited to one-dimensional (1D) situations, where they are naturally stable. For instance, in 1D cold Bose gases, they exist for any attractive interaction strength g and particle number N. By contrast, in two dimensions, solitons appear only for discrete values of gN, the so-called Townes soliton being the most celebrated example. Here, we use a two-component Bose gas to prepare deterministically such a soliton: Starting from a uniform bath of atoms in a given internal state, we imprint the soliton wave function using an optical transfer to another state. We explore various interaction strengths, atom numbers, and sizes and confirm the existence of a solitonic behavior for a specific value of gN and arbitrary sizes, a hallmark of scale invariance.

Rotating Multidimensional Quantum Droplets

We predict a new type of two- and three-dimensional stable quantum droplets persistently rotating in broad external two-dimensional and weakly anharmonic potential. Their evolution is described by the system of the Gross-Pitaevskii equations with Lee-Huang-Yang quantum corrections. Such droplets resemble whispering-gallery modes localized in the polar direction due to nonlinear interactions and, depending on their chemical potential and rotation frequency, they appear in rich variety of shapes, ranging from nearly flat-top or strongly localized rotating wave packets, to crescentlike objects extending nearly over the entire range of polar angles. Above critical rotation frequency quantum droplets transform into vortex droplets (in two dimensions) or vortex tori (in three dimensions), whose topological charge gradually increase with the increase of the modulus of chemical potential, and therefore they belong to the family of nonlinear modes connecting fundamental and vortex quantum droplets. Rotating quantum droplets are exceptionally robust objects, stable practically in the entire range of their existence.

Observation of a Lee-Huang-Yang Fluid

We observe monopole oscillations in a mixture of Bose-Einstein condensates, where the usually dominant mean-field interactions are canceled. In this case, the system is governed by the next-order Lee-Huang-Yang (LHY) correction to the ground state energy, which describes the effect of quantum fluctuations. Experimentally such a LHY fluid is realized by controlling the atom numbers and interaction strengths in a ^{39}K spin mixture confined in a spherical trap potential. We measure the monopole oscillation frequency as a function of the LHY interaction strength as proposed recently by Jrgensen et al. [Phys. Rev. Lett. 121, 173403 (2018)PRLTAO0031-900710.1103/PhysRevLett.121.173403] and find excellent agreement with simulations of the complete experiment including the excitation procedure and inelastic losses. This confirms that the system and its collective behavior are initially dominated by LHY interactions. Moreover, the monopole oscillation frequency is found to be stable against variations of the involved scattering lengths in a broad region around the ideal values, confirming the stabilizing effect of the LHY interaction. These results pave the way for using the nonlinearity provided by the LHY term in quantum simulation experiments and for investigations beyond the LHY regime.

Mixed Bubbles in Bose-Bose Mixtures

Repulsive Bose-Bose mixtures are known to either mix or phase separate into pure components. Here we predict a mixed-bubble regime in which bubbles of the mixed phase coexist with a pure phase of one of the components. This is a beyond-mean-field effect that occurs for unequal masses or unequal intraspecies coupling constants and is due to a competition between the mean-field term, quadratic in densities, and a nonquadratic beyond-mean-field correction. We find parameters of the mixed-bubble regime in all dimensions and discuss implications for current experiments.

Self-Evaporation Dynamics of Quantum Droplets in a 41K-87Rb Mixture

We theoretically investigate the self-evaporation dynamics of quantum droplets in a 41K-87Rb mixture, in free-space. The dynamical formation of the droplet and the effects related to the presence of three-body losses are analyzed by means of numerical simulations. We identify a regime of parameters allowing for the observation of the droplet self-evaporation in a feasible experimental setup.

Universal Dimerized Quantum Droplets in a One-Dimensional Lattice

The ground-state properties of two-component bosonic mixtures in a one-dimensional optical lattice are studied both from few- and many-body perspectives. We rely directly on a microscopic Hamiltonian with attractive intercomponent and repulsive intracomponent interactions to demonstrate the formation of a quantum liquid. We reveal that its formation and stability can be interpreted in terms of finite-range interactions between dimers. We derive an effective model of composite bosons (dimers) which correctly captures both the few- and many-body properties and validate it against exact results obtained by the density matrix renormalization group method for the full Hamiltonian. The threshold for the formation of the liquid coincides with the appearance of a bound state in the dimer-dimer problem and possesses a universality in terms of the two-body parameters of the dimer-dimer interaction, namely, scattering length and effective range. For sufficiently strong effective dimer-dimer repulsion we observe fermionization of the dimers which form an effective Tonks-Girardeau state and identify conditions for the formation of a solitonic solution. Our predictions are relevant to experiments with dipolar atoms and two-component mixtures.

New states of matter with fine-tuned interactions: quantum droplets and dipolar supersolids

Quantum fluctuations can stabilize Bose-Einstein condensates (BEC) against the mean-field collapse. Stabilization of the condensate has been observed in quantum degenerate Bose-Bose mixtures and dipolar BECs. The fine-tuning of the interatomic interactions can lead to the emergence of two new states of matter: liquid-like self-bound quantum droplets and supersolid crystals formed from these droplets. We review the properties of these exotic states of matter and summarize the experimental progress made using dipolar quantum gases and Bose-Bose mixtures. We conclude with an outline of important open questions that could be addressed in the future.

Consistent Theory of Self-Bound Quantum Droplets with Bosonic Pairing

We revisit the Bogoliubov theory of quantum droplets proposed by Petrov [Phys. Rev. Lett. 115, 155302 (2015)PRLTAO0031-900710.1103/PhysRevLett.115.155302] for an ultracold Bose-Bose mixture, where the mean-field collapse is stabilized by the Lee-Huang-Yang quantum fluctuations. We show that a loophole in Petrov's theory, i.e., the ignorance of the softening complex Bogoliubov spectrum, can be naturally removed by the introduction of bosonic pairing. The pairing leads to weaker mean-field attractions, and also a stronger Lee-Huang-Yang term in the case of unequal intraspecies interactions. As a result, the equilibrium density for the formation of self-bound droplets significantly decreases in the deep droplet regime, in agreement with a recent observation from diffusion Monte Carlo simulations. Our construction of a consistent Bogoliubov theory paves the way to understand the puzzling low critical number of small quantum droplets observed in the experiment [C. Cabrera et al., Science 359, 301 (2018)SCIEAS0036-807510.1126/science.aao5686].

Normal-Superfluid Phase Separation in Spin-Half Bosons at Finite Temperature

For pseudospin-half bosons with interspin attraction and intraspin repulsion, the normal phase and Bose condensed phase can coexist at finite temperature. The homogeneous system is unstable against the spinodal decomposition within a medium density interval, and, consequently, a normal-superfluid phase separation takes place. The isothermal equation of state shows a characteristic plateau in the P-V (pressure-volume) diagram, which is reminiscent of a classical gas-liquid transition, although, unlike the latter, the coexistence lines never terminate at a critical point as temperature increases. In a harmonic trap, the phase separation can be revealed by the density profile of the atomic cloud, which exhibits a sudden jump across the phase boundary.

Quantum Fluctuations of the Center of Mass and Relative Parameters of Nonlinear Schrödinger Breathers

We study quantum fluctuations of macroscopic parameters of a nonlinear Schrödinger breather-a nonlinear superposition of two solitons, which can be created by the application of a fourfold quench of the scattering length to the fundamental soliton in a self-attractive quasi-one-dimensional Bose gas. The fluctuations are analyzed in the framework of the Bogoliubov approach in the limit of a large number of atoms N, using two models of the vacuum state: white noise and correlated noise. The latter model, closer to the ab initio setting by construction, leads to a reasonable agreement, within 20% accuracy, with fluctuations of the relative velocity of constituent solitons obtained from the exact Bethe-ansatz results [Phys. Rev. Lett. 119, 220401 (2017)PRLTAO0031-900710.1103/PhysRevLett.119.220401] in the opposite low-N limit (for N≤23). We thus confirm, for macroscopic N, the breather dissociation time to be within the limits of current cold-atom experiments. Fluctuations of soliton masses, phases, and positions are also evaluated and may have experimental implications.

Static and Dynamic Properties of a Few Spin 1/2 Interacting Fermions Trapped in a Harmonic Potential

We provide a detailed study of the properties of a few interacting spin 1 / 2 fermions trapped in a one-dimensional harmonic oscillator potential. The interaction is assumed to be well represented by a contact delta potential. Numerical results obtained by means of direct diagonalization techniques are combined with analytical expressions for both the non-interacting and strongly interacting regime. The N = 2 case is used to benchmark our numerical techniques with the known exact solution of the problem. After a detailed description of the numerical methods, in a tutorial-like manner, we present the static properties of the system for N = 2 , 3 , 4 and 5 particles, e.g., low-energy spectrum, one-body density matrix, ground-state densities. Then, we consider dynamical properties of the system exploring first the excitation of the breathing mode, using the dynamical structure function and corresponding sum-rules, and then a sudden quench of the interaction strength.

van der Waals Universality near a Quantum Tricritical Point

We study the three-body scattering hypervolume D of atoms whose scattering length a is on the order of or smaller than the typical range r_{vdW} of the van der Waals attraction. We find that the real part of D behaves universally in this weakly interacting regime (|a|/r_{vdW}≲1) in the absence of trimer resonances. This universality originates from hard-spherelike collisions that dominate elastic three-body scattering. We use this result to make quantitative predictions for the thermodynamics and elementary excitations of an atomic Bose-Einstein condensate in the vicinity of a quantum tricritical point, including quantum droplets stabilized by effective three-body interactions.

A Dual-Species Bose-Einstein Condensate with Attractive Interspecies Interactions

We report on the production of a 41 K- 87 Rb dual-species Bose–Einstein condensate with tunable interspecies interaction and we study the mixture in the attractive regime; i.e., for negative values of the interspecies scattering length a 12 . The binary condensate is prepared in the ground state and confined in a pure optical trap. We exploit Feshbach resonances for tuning the value of a 12 . After compensating the gravitational sag between the two species with a magnetic field gradient, we drive the mixture into the attractive regime. We let the system evolve both in free space and in an optical waveguide. In both geometries, for strong attractive interactions, we observe the formation of self-bound states, recognizable as quantum droplets. Our findings prove that robust, long-lived droplet states can be realized in attractive two-species mixtures, despite the two atomic components possibly experiencing different potentials.

Singular Mean-Field States: A Brief Review of Recent Results

This article provides a focused review of recent findings which demonstrate, in some cases quite counter-intuitively, the existence of bound states with a singularity of the density pattern at the center; the states are physically meaningful because their total norm converges. One model of this type is based on the 2D Gross–Pitaevskii equation (GPE), which combines the attractive potential ∼ r − 2 and the quartic self-repulsive nonlinearity, induced by the Lee–Huang–Yang effect (quantum fluctuations around the mean-field state). The GPE demonstrates suppression of the 2D quantum collapse, driven by the attractive potential, and emergence of a stable ground state (GS), whose density features an integrable singularity ∼ r − 4 / 3 at r → 0 . Modes with embedded angular momentum exist too, but they are unstable. A counter-intuitive peculiarity of the model is that the GS exists even if the sign of the potential is reversed from attraction to repulsion, provided that its strength is small enough. This peculiarity finds a relevant explanation. The other model outlined in the review includes 1D, 2D, and 3D GPEs, with the septimal (seventh-order), quintic, and cubic self-repulsive terms, respectively. These equations give rise to stable singular solitons, which represent the GS for each dimension D, with the density singularity ∼ r − 2 / ( 4 − D ) . Such states may be considered the results of screening a “bare” delta-functional attractive potential by the respective nonlinearities.

Strongly Correlated Quantum Droplets in Quasi-1D Dipolar Bose Gas

We exploit a few- to many-body approach to study strongly interacting dipolar bosons in the quasi-one-dimensional system. The dipoles attract each other while the short range interactions are repulsive. Solving numerically the multiatom Schrödinger equation, we discover that such systems can exhibit not only the well-known bright soliton solutions but also novel quantum droplets for a strongly coupled case. For larger systems, basing on microscopic properties of the found few-body solution, we propose a new equation for a density amplitude of atoms. It accounts for fermionization for strongly repelling bosons by incorporating the Lieb-Liniger energy in a local density approximation and approaches the standard Gross-Pitaevskii equation (GPE) in the weakly interacting limit. Not only does such a framework provide an alternative mechanism of the droplet stability, but it also introduces means to further analyze this previously unexplored quantum phase. In the limiting strong repulsion case, yet another simple multiatom model is proposed. We stress that the celebrated Lee-Huang-Yang term in the GPE is not applicable in this case.

Singular solitons

We demonstrate that the commonly known concept which treats solitons as nonsingular solutions produced by the interplay of nonlinear self-attraction and linear dispersion may be extended to include modes with a relatively weak singularity at the central point, which keeps their integral norm convergent. Such states are generated by self-repulsion, which should be strong enough, represented by septimal, quintic, and usual cubic terms in the framework of the one-, two-, and three-dimensional (1D, 2D, and 3D) nonlinear Schrödinger equations (NLSEs), respectively. Although such solutions seem counterintuitive, we demonstrate that they admit a straightforward interpretation as a result of screening of an additionally introduced attractive δ-functional potential by the defocusing nonlinearity. The strength ("bare charge") of the attractive potential is infinite in 1D, finite in 2D, and vanishingly small in 3D. Analytical asymptotics of the singular solitons at small and large distances are found, entire shapes of the solitons being produced in a numerical form. Complete stability of the singular modes is accurately predicted by the anti-Vakhitov-Kolokolov criterion (under the assumption that it applies to the model), as verified by means of numerical methods. In 2D, the NLSE with a quintic self-focusing term admits singular-soliton solutions with intrinsic vorticity too, but they are fully unstable. We also mention that dissipative singular solitons can be produced by the model with a complex coefficient in front of the nonlinear term.

Rotating a Supersolid Dipolar Gas

Distinctive features of supersolids show up in their rotational properties. We calculate the moment of inertia of a harmonically trapped dipolar Bose-Einstein condensed gas as a function of the tunable scattering length parameter, providing the transition from the (fully) superfluid to the supersolid phase and eventually to an incoherent crystal of self-bound droplets. The transition from the superfluid to the supersolid phase is characterized by a jump in the moment of inertia, revealing its first order nature. In the case of elongated trapping in the plane of rotation, we show that the moment of inertia determines the value of the frequency of the scissors mode, which is significantly affected by the reduction of superfluidity in the supersolid phase. The case of an in-plane isotropic trapping is instead well suited to study the formation of quantized vortices, which are shown to be characterized, in the supersolid phase, by a sizeable deformed core, caused by the presence of the surrounding density peaks.

Modulational Instability, Inter-Component Asymmetry, and Formation of Quantum Droplets in One-Dimensional Binary Bose Gases

Quantum droplets are ultradilute liquid states that emerge from the competitive interplay of two Hamiltonian terms, the mean-field energy and beyond-mean-field correction, in a weakly interacting binary Bose gas. We relate the formation of droplets in symmetric and asymmetric two-component one-dimensional boson systems to the modulational instability of a spatially uniform state driven by the beyond-mean-field term. Asymmetry between the components may be caused by their unequal populations or unequal intra-component interaction strengths. Stability of both symmetric and asymmetric droplets is investigated. Robustness of the symmetric solutions against symmetry-breaking perturbations is confirmed.

One-dimensional mixtures of several ultracold atoms: a review

Recent theoretical and experimental progress on studying one-dimensional systems of bosonic, fermionic, and Bose-Fermi mixtures of a few ultracold atoms confined in traps is reviewed in the broad context of mesoscopic quantum physics. We pay special attention to limiting cases of very strong or very weak interactions and transitions between them. For bosonic mixtures, we describe the developments in systems of three and four atoms as well as different extensions to larger numbers of particles. We also briefly review progress in the case of spinor Bose gases of a few atoms. For fermionic mixtures, we discuss a special role of spin and present a detailed discussion of the two- and three-atom cases. We discuss the advantages and disadvantages of different computation methods applied to systems with intermediate interactions. In the case of very strong repulsion, close to the infinite limit, we discuss approaches based on effective spin chain descriptions. We also report on recent studies on higher-spin mixtures and inter-component attractive forces. For both statistics, we pay particular attention to impurity problems and mass imbalance cases. Finally, we describe the recent advances on trapped Bose-Fermi mixtures, which allow for a theoretical combination of previous concepts, well illustrating the importance of quantum statistics and inter-particle interactions. Lastly, we report on fundamental questions related to the subject which we believe will inspire further theoretical developments and experimental verification.

Semidiscrete Quantum Droplets and Vortices

We consider a binary bosonic condensate with weak mean-field (MF) residual repulsion, loaded in an array of nearly one-dimensional traps coupled by transverse hopping. With the MF force balanced by the effectively one-dimensional attraction, induced in each trap by the Lee-Hung-Yang correction (produced by quantum fluctuations around the MF state), stable on-site- and intersite-centered semidiscrete quantum droplets (QDs) emerge in the array, as fundamental ones and self-trapped vortices, with winding numbers, at least, up to five, in both tightly bound and quasicontinuum forms. The application of a relatively strong trapping potential leads to squeezing transitions, which increase the number of sites in fundamental QDs and eventually replace vortex modes by fundamental or dipole ones. The results provide the first realization of stable semidiscrete vortex QDs, including ones with multiple vorticity.

Nondegenerate bound-state solitons in multicomponent Bose-Einstein condensates

We investigate nondegenerate bound-state solitons systematically in multicomponent Bose-Einstein condensates, through developing the Darboux transformation method to derive exact soliton solutions analytically. In particular, we show that bright solitons with nodes correspond to the excited bound states in effective quantum wells, in sharp contrast to the bright solitons and dark solitons reported before (which usually correspond to ground state and free state, respectively). We further demonstrate that bound-state solitons with nodes are induced by incoherent superposition of solitons in different components. Moreover, we reveal that the interactions between these bound-state solitons are usually inelastic, caused by the incoherent interactions between solitons in different components and the coherent interactions between solitons in the same component. Additionally, the detailed spectral stability analysis demonstrates the stability of nondegenerate bound-state solitons. The bound-state solitons can be used to study many different physical problems, such as beating dynamics, spin-orbit coupling effects, quantum fluctuations, and even quantum entanglement states.

Metastability of Quantum Droplet Clusters

We show that metastable ring-shaped clusters can be constructed from two-dimensional quantum droplets in systems described by the Gross-Pitaevskii equations augmented with Lee-Huang-Yang quantum corrections. The clusters exhibit dynamical behavior ranging from contraction to rotation with simultaneous periodic pulsations, or expansion, depending on the initial radius of the necklace pattern and phase shift between adjacent quantum droplets. We show that, using an energy-minimization analysis, one can predict equilibrium values of the cluster radius that correspond to rotation without radial pulsations. In such a regime, the clusters evolve as metastable states, withstanding abrupt variations in the underlying scattering lengths and keeping their azimuthal symmetry in the course of evolution, even in the presence of considerable perturbations.

Directional soliton and breather beams

Solitons and breathers are nonlinear modes that exist in a wide range of physical systems. They are fundamental solutions of a number of nonlinear wave evolution equations, including the unidirectional nonlinear Schrödinger equation (NLSE). We report the observation of slanted solitons and breathers propagating at an angle with respect to the direction of propagation of the wave field. As the coherence is diagonal, the scale in the crest direction becomes finite; consequently, beam dynamics form. Spatiotemporal measurements of the water surface elevation are obtained by stereo-reconstructing the positions of the floating markers placed on a regular lattice and recorded with two synchronized high-speed cameras. Experimental results, based on the predictions obtained from the (2D + 1) hyperbolic NLSE equation, are in excellent agreement with the theory. Our study proves the existence of such unique and coherent wave packets and has serious implications for practical applications in optical sciences and physical oceanography. Moreover, unstable wave fields in this geometry may explain the formation of directional large-amplitude rogue waves with a finite crest length within a wide range of nonlinear dispersive media, such as Bose-Einstein condensates, solids, plasma, hydrodynamics, and optics.

Liquid State of One-Dimensional Bose Mixtures: A Quantum Monte Carlo Study

By using exact quantum Monte Carlo methods we calculate the ground-state properties of the liquid phase in one-dimensional Bose mixtures with contact interactions. We find that the liquid state can be formed if the ratio of coupling strengths between interspecies attractive and intraspecies repulsive interactions exceeds a critical value. As a function of this ratio we determine the density where the energy per particle has a minimum and the one where the compressibility diverges, thereby identifying the equilibrium density and the spinodal point in the phase diagram of the homogeneous liquid. Furthermore, in the stable liquid state, we calculate the chemical potential, the speed of sound, as well as structural and coherence properties, such as the pair correlation function, the static structure factor, and the one-body density matrix, thus providing a detailed description of the bulk region in self-bound droplets.

Collisions of Self-Bound Quantum Droplets

We report on the study of binary collisions between quantum droplets formed by an attractive mixture of ultracold atoms. We distinguish two main outcomes of the collision, i.e., merging and separation, depending on the velocity of the colliding pair. The critical velocity v_{c} that discriminates between the two cases displays a different dependence on the atom number N for small and large droplets. By comparing our experimental results with numerical simulations, we show that the nonmonotonic behavior of v_{c}(N) is due to the crossover from a compressible to an incompressible regime, where the collisional dynamics is governed by different energy scales, i.e., the droplet binding energy and the surface tension. These results also provide the first evidence of the liquidlike nature of quantum droplets in the large N limit, where their behavior closely resembles that of classical liquid droplets.

Dilute Fluid Governed by Quantum Fluctuations

Understanding the effects of interactions in complex quantum systems beyond the mean-field paradigm constitutes a fundamental problem in physics. We show how the atom numbers and interactions in a Bose-Bose mixture can be tuned to cancel mean-field interactions completely. The resulting system is entirely governed by quantum fluctuations-specifically the Lee-Huang-Yang correlations. We derive an effective one-component Gross-Pitaevskii equation for this system, which is shown to be very accurate by comparison with a full two-component description. This allows us to show how the Lee-Huang-Yang correlation energy can be accurately measured using two powerful probes of atomic gases: collective excitations and radio-frequency spectroscopy. Importantly, the behavior of the system is robust against deviations from the atom number and interaction criteria for canceling the mean-field interactions. This shows that it is feasible to realize a setting where quantum fluctuations are not masked by mean-field forces, allowing investigations of the Lee-Huang-Yang correction at unprecedented precision.

Observation of Photon Droplets and Their Dynamics

We present experimental evidence of photon droplets in an attractive (focusing) nonlocal nonlinear medium. Photon droplets are self-bound, finite-sized states of light that are robust to size and shape perturbations due to a balance of competing attractive and repulsive forces. It has recently been shown theoretically, via a multipole expansion of the nonlocal nonlinearity, that the self-bound state arises due to competition between the s-wave and d-wave nonlinear terms, together with diffraction. The theoretical photon droplet framework encompasses both a solitonlike stationary ground state and the nonsolitonlike dynamics that ensue when the system is displaced from equilibrium, i.e., driven into an excited state. We present numerics and experiments supporting the existence of these photon droplet states and measurements of the dynamical evolution of the photon droplet orbital angular momentum.

Quantum Spin Dynamics of Individual Neutral Impurities Coupled to a Bose-Einstein Condensate

We report on spin dynamics of individual, localized neutral impurities immersed in a Bose-Einstein condensate. Single cesium atoms are transported into a cloud of rubidium atoms and thermalize with the bath, and the ensuing spin exchange between localized impurities with quasispin F_{i}=3 and bath atoms with F_{b}=1 is resolved. Comparing our data to numerical simulations of spin dynamics, we find that, for gas densities in the Bose-Einstein condensate regime, the dynamics is dominated by the condensed fraction of the cloud. We spatially resolve the density overlap of impurities and gas by the spin population of impurities. Finally, we trace the coherence of impurities prepared in a coherent superposition of internal states when coupled to a gas of different densities. For our choice of states, we show that, despite high bath densities and, thus, fast thermalization rates, the impurity coherence is not affected by the bath, realizing a regime of sympathetic cooling while maintaining internal state coherence. Our work paves the way toward the nondestructive probing of quantum many-body systems via localized impurities.

Self-Bound Quantum Droplets of Atomic Mixtures in Free Space

Self-bound quantum droplets are a newly discovered phase in the context of ultracold atoms. In this Letter, we report their experimental realization following the original proposal by Petrov [Phys. Rev. Lett. 115, 155302 (2015)PRLTAO0031-900710.1103/PhysRevLett.115.155302], using an attractive bosonic mixture. In this system, spherical droplets form due to the balance of competing attractive and repulsive forces, provided by the mean-field energy close to the collapse threshold and the first-order correction due to quantum fluctuations. Thanks to an optical levitating potential with negligible residual confinement, we observe self-bound droplets in free space, and we characterize the conditions for their formation as well as their size and composition. This work sets the stage for future studies on quantum droplets, from the measurement of their peculiar excitation spectrum to the exploration of their superfluid nature.

Scissors Mode of Dipolar Quantum Droplets of Dysprosium Atoms

We report on the observation of the scissors mode of a single dipolar quantum droplet. The existence of this mode is due to the breaking of the rotational symmetry by the dipole-dipole interaction, which is fixed along an external homogeneous magnetic field. By modulating the orientation of this magnetic field, we introduce a new spectroscopic technique for studying dipolar quantum droplets. This provides a precise probe for interactions in the system, allowing us to extract a background scattering length for ^{164}Dy of 69(4)a_{0}. Our results establish an analogy between quantum droplets and atomic nuclei, where the existence of the scissors mode is also only due to internal interactions. They further open the possibility to explore physics beyond the available theoretical models for strongly dipolar quantum gases.