Binary Condensates Research Articles

Formation of compact objects at finite temperatures in a dark-matter-candidate self-gravitating bosonic system

We study self-gravitating bosonic systems, candidates for dark-matter halos, by carrying out a suite of direct numerical simulations (DNSs) designed to investigate the formation of finitetemperature, compact objects in the three-dimensional (3D) Fourier-truncated Gross-Pitaevskii-Poisson equation (GPPE). This truncation allows us to explore the collapse and fluctuations of compact objects, which form at both zero temperature and finite temperature. We show that the statistically steady state of the GPPE, in the large-time limit and for the system sizes we study, can also be obtained efficiently by tuning the temperature in an auxiliary stochastic Ginzburg-Landau-Poisson equation (SGLPE). We show that, over a wide range of model parameters, this system undergoes a thermally driven first-order transition from a collapsed, compact, Bose-Einstein condensate (BEC) to a tenuous Bose gas without condensation. By a suitable choice of initial conditions in the GPPE, we also obtain a binary condensate that comprises a pair of collapsed objects rotating around their center of mass.

Oscillatory stability of quantum droplets in -symmetric optical lattice

We study stability and collisions of quantum droplets (QDs) forming in a binary bosonic condensate trapped in parity-time ()-symmetric optical lattices. It is found that the stability of QDs in the -symmetric system depends strongly on the values of the imaginary part W 0 of the -symmetric optical lattices, self-repulsion strength g, and the condensate norm N. As expected, the -symmetric QDs are entirely unstable in the broken -symmetric phase. However, the -symmetric QDs exhibit oscillatory stability with the increase of N and g in the unbroken -symmetric phase. Finally, collisions between -symmetric QDs are considered. The collisions of droplets with unequal norms are completely different from that in free space. Besides, a stable -symmetric QDs collides with an unstable ones tend to merge into breathers after the collision.

Dark-soliton-like magnetic domain walls in a two-dimensional ferromagnetic superfluid

We report a stable magnetic domain wall in a uniform ferromagnetic spin-1 condensate, characterized by the magnetization having a dark soliton profile with nonvanishing superfluid density. We find exact stationary solutions for a particular ratio of interaction parameters with and without magnetic fields, and develop an accurate analytic solution applicable to the whole ferromagnetic phase. In the absence of magnetic fields, this domain wall relates various distinct solitary excitations in binary condensates through $\textrm{SO}(3)$ spin rotations, which otherwise are unconnected. Remarkably, studying the dynamics of a quasi-two-dimensional (quasi-2D) system we show that standing wave excitations of the domain wall oscillate without decay, being stable against the snake instability. The domain wall is dynamically unstable to modes that cause the magnetization to grow perpendicularly while leaving the domain wall unchanged. Real time dynamics in the presence of white noise reveals that this "spin twist" instability does not destroy the topological structure of the magnetic domain wall.

Ultracold Bose mixtures with spin-dependent fermion-mediated interactions

We develop a functional integral formulation for binary Bose-Einstein condensates coupled to polarized fermions. We find that spin-dependent fermion-mediated interactions have dramatic effects on the properties of the binary condensates. The quasiparticle spectrum features two branches. The upper branch, which is of density nature, gets modified by the induced interactions, while the lower branch, which is of spin nature, is left intact. The ground-state phase diagram consists of stable region of miscible phases and unstable region toward phase separation. In the stable region, it is further classified by the damping of excitations of the upper branch. We show that it is possible to find region of well-defined, long-lived quasiparticle excitations by tuning relevant parameters, such as boson-fermion mass ratio, boson-fermion number density ratio, and interspecies interactions between bosons as well. We explore the effects of quantum fluctuation due to the effective potential on the binary condensates. It turns out that both the density structure factor and spin density structure factor fulfill the Feynman relation, except that the latter is immune to the fermion-mediated interactions.

Breathing modes of repulsive polarons in Bose–Bose mixtures

We consider impurity atoms embedded in a two-component Bose–Einstein condensate in a quasi-one dimensional regime. We study the effects of repulsive coupling between the impurities and Bose species on the equilibrium of the system for both miscible and immiscible mixtures by numerically solving the underlying coupled Gross–Pitaevskii equations. Our results reveal that the presence of impurities may lead to a miscible–immiscible phase transition due to the interaction of the impurities and the two condensates. Within the realm of the Bogoliubov–de Gennes equations we calculate the quantum fluctuations due to the different types of interactions. The breathing modes and the time evolution of harmonically trapped impurities in both homogeneous and inhomogeneous binary condensates are deeply discussed in the miscible case using variational and numerical means. We show in particular that the self-trapping, the miscibility and the inhomogeneity of the trapped Bose mixture may strongly modify the low-lying excitations and the dynamical properties of impurities. The presence of phonons in the homogeneous Bose mixture gives rise to the damping of breathing oscillations of impurities width.

Stability and collisions of quantum droplets in [formula omitted]-symmetric dual-core couplers

We study the effect of the interplay between parity-time (PT) symmetry and optical lattice (OL) potential on dynamics of quantum droplets (QDs) forming in a binary bosonic condensate trapped in a dual-core system. It is found that the stability of symmetric QDs in such non-Hermitian system depends critically on the competition of gain and loss γ, inter-core coupling κ, and OL potential. In the absence of OL potential, the PT-symmetric QDs are unstable against symmetry-breaking perturbations with the increase of the total condensate norm N, and they retrieve the stability at larger N, in the weakly-coupled regime. As expected, the stable region of the PT-symmetric QDs shrinks when γ increases, i.e., the PT symmetry is prone to break the stability of QDs. There is a critical value of κ beyond which the PT-symmetric QDs are entirely stable in the unbroken PT-symmetric phase. In the presence of OL potential, the PT-symmetric on-site QDs are still stable for relatively small and large values of N. Nevertheless, it is demonstrated that the OL potential can assist stabilization of PT-symmetric on-site QDs for some moderate values of N. On the other hand, it is worth noting that the relatively small PT-symmetric off-site QDs are unstable, and only the relatively large ones are stable. Furthermore, collisions between stable PT-symmetric QDs are considered too. It is revealed that the slowly moving PT-symmetric QDs tend to merge into breathers, while the fast-moving ones display quasi-elastic collision and suffer fragmentation for small and large values of N, respectively.

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.

Dipolar condensed atomic mixtures and miscibility under rotation

By considering symmetric- and asymmetric-dipolar coupled mixtures (with dysprosium and erbium isotopes), we report a study on relevant anisotropic effects, related to spatial separation and miscibility, due to dipole-dipole interactions (DDIs) in rotating binary dipolar Bose-Einstein condensates. The binary mixtures are kept in strong pancake-like traps, with repulsive two-body interactions modeled by an effective two-dimensional (2D) coupled Gross-Pitaevskii equation. The DDI are tuned from repulsive to attractive by varying the dipole polarization angle. A clear spatial separation is verified in the densities for attractive DDIs, being angular for symmetric mixtures and radial for asymmetric ones. Also relevant is the mass-imbalance sensibility observed by the vortex-patterns in symmetric- and asymmetric-dipolar mixtures. In an extension of this study, here we show how the rotational properties and spatial separation of these dipolar mixture are affected by a quartic term added to the harmonic trap of one of the components.

Gap and embedded solitons in microwave-coupled binary condensates

It was recently found that, under the action of the spin-orbit coupling (SOC) and Zeeman splitting (ZS), binary BEC with intrinsic cubic nonlinearity supports families of gap solitons, provided that the kinetic energy is negligible in comparison with the SOC and ZS terms. We demonstrate that, also under the action of SOC and ZS, a similar setting may be introduced for BEC with two components representing different atomic states, resonantly coupled by microwave radiation, while the Poisson equation accounts for the feedback of the two-component atomic wave function onto the radiation. The microwave-mediated interaction induces an effective nonlinear trapping potential, which strongly affects the purport of the linear spectrum in this system. As a result, families of both gap and embedded solitons (those overlapping with the continuous spectrum) are found, being chiefly stable. The shape of the solitons features exact or broken skew symmetry. In addition to fundamental solitons (whose shape may or may not include a node), a family of dipole solitons is constructed too, which are even more stable than their fundamental counterparts. A nontrivial stability area is identified for moving solitons in the present system, which lacks Galilean invariance. Colliding solitons merge into a single one.

Rotating Binary Bose-Einstein Condensates and Vortex Clusters in Quantum Droplets.

Quantum droplets may form out of a gaseous Bose-Einstein condensate, stabilized by quantum fluctuations beyond mean field. Determining the ground state of a rotating binary condensate, we show that multiple singly quantized vortices may form in these droplets at moderate angular momenta in two dimensions. Droplets carrying these precursors of an Abrikosov lattice remain self-bound for certain timescales after switching off an initial harmonic confinement. Intriguingly, we find evidence of a metastable persistent current in these new types of binary condensates. We discuss how this finding can be used to experimentally generate vortex-carrying quantum droplets.

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.

Immiscible and miscible states in binary condensates in the ring geometry

We report detailed investigation of the existence and stability of mixed and demixed modes in binary atomic Bose–Einstein condensates with repulsive interactions in a ring-trap geometry. The stability of such states is examined through eigenvalue spectra for small perturbations, produced by the Bogoliubov–de Gennes equations, and directly verified by simulations based on the coupled Gross–Pitaevskii equations, varying inter- and intra-species scattering lengths so as to probe the entire range of miscibility–immiscibility transitions. In the limit of the one-dimensional (1D) ring, i.e. a very narrow one, stability of mixed states is studied analytically, including hidden-vorticity (HV) modes, i.e. those with opposite vorticities of the two components and zero total angular momentum. The consideration of demixed 1D states reveals, in addition to stable composite single-peak structures, double- and triple-peak ones, above a certain particle-number threshold. In the 2D annular geometry, stable demixed states exist both in radial and azimuthal configurations. We find that stable radially-demixed states can carry arbitrary vorticity and, counter-intuitively, the increase of the vorticity enhances stability of such states, while unstable ones evolve into randomly oscillating angular demixed modes. The consideration of HV states in the 2D geometry expands the stability range of radially-demixed states.

Spatial separation of rotating binary Bose-Einstein condensates by tuning the dipolar interactions

We are pointing out relevant anisotropic effects, related to spatial separation, miscibility and mass-symmetry, due to dipole-dipole interactions in rotating binary dipolar Bose-Einstein condensates, by considering symmetric ($^{164}$Dy-$^{162}$Dy) and asymmetric ($^{168}$Er-$^{164}$Dy, $^{164}$Dy-$^{87}$Rb) dipolar mixtures. The binary mixtures are kept in strong pancake-shaped trap, modeled by an effective two-dimensional coupled Gross-Pitaevskii equation. The anisotropy of the dipolar interactions, on miscibility and vortex-lattice structures, is studied by tuning the polarization angle of the dipoles $\varphi$, which can enhance the attractive part of the dipole-dipole interaction (DDI) for both inter- and intra-species. Within this procedure of changing to attractive the DDI, a clear spatial separation is verified in the densities at some critical polarization angle. The spatial separations, being angular for symmetric mixtures and radial for asymmetric ones, are verified for repulsive contact interactions when the inter- to intra-species ratio $\delta$ is larger than one, implying the system is less miscible. The corresponding result for the critical polarization angle as a function of $\delta$ is shown in the particular dipolar symmetric case. A striking outcome of the present study is the observed sensibility of the vortex-pattern binary distributions due to the mass-asymmetry. This is exemplified by the symmetric dipolar mixture, where the two isotopes are of the same species.

Three-dimensional solitons in Bose-Einstein condensates with spin-orbit coupling and Bessel optical lattices

Applying a variational approach and numerical analysis to the system of Gross-Pitaevskii equations, we find three-dimensional (3D) stable solitons in binary atomic Bose-Einstein condensates with spin-orbit coupling (SOC) and out-of-phase linear and nonlinear Bessel optical lattices. We discuss the stability of 3D solitons by utilizing their norm and energy. The introduction of Bessel potentials makes the evolution and collisions of solitons more stable and improves their resistance to collapse. Depending on the strength of the intra- and intercomponent spatial modulation of the nonlinearity and SOC, we find stable solitons of the semivortex and mixed-mode structures. Furthermore, we show that the solitons are stable against small perturbations in propagation and collisions.

Suppression of Quantum-Mechanical Collapse in Bosonic Gases with Intrinsic Repulsion: A Brief Review

It is known that attractive potential ~− 1 / r 2 gives rise to the critical quantum collapse in the framework of the three-dimensional (3D) linear Schrödinger equation. This article summarizes theoretical analysis, chiefly published in several original papers, which demonstrates suppression of the collapse caused by this potential, and the creation of the otherwise missing ground state in a 3D gas of bosonic dipoles pulled by the same potential to the central charge, with repulsive contact interactions between them, represented by the cubic term in the respective Gross–Pitaevskii equation (GPE). In two dimensions (2D), quintic self-repulsion is necessary for the suppression of the collapse; alternatively, this may be provided by the effective quartic repulsion produced by the Lee–Huang–Yang correction to the GPE. 3D states carrying angular momentum are constructed in the model with the symmetry reduced from spherical to cylindrical by an external polarizing field. Interplay of the collapse suppression and miscibility–immiscibility transition is considered in a binary condensate. The consideration of the 3D setting in the form of the many-body quantum system, with the help of the Monte Carlo method, demonstrates that, although the quantum collapse cannot be fully suppressed, the self-trapped states predicted by the GPE exist in the many-body setting as metastable modes protected against the collapse by a tall potential barrier.

Domain-area distribution anomaly in segregating multicomponent superfluids

The domain-area distribution in the phase transition dynamics of ${Z}_{2}$ symmetry breaking is studied theoretically and numerically for segregating binary Bose-Einstein condensates in quasi-two-dimensional systems. Due to the dynamic-scaling law of the phase ordering kinetics, the domain-area distribution is described by a universal function of the domain area, rescaled by the mean distance between domain walls. The scaling theory for general coarsening dynamics in two dimensions hypothesizes that the distribution during the coarsening dynamics has a hierarchy with the two scaling regimes, the microscopic and macroscopic regimes with distinct power-law exponents. The power law in the macroscopic regime, where the domain size is larger than the mean distance, is universally represented with the Fisher's exponent of the percolation theory in two dimensions. On the other hand, the power-law exponent in the microscopic regime is sensitive to the microscopic dynamics of the system. This conjecture is confirmed by large-scale numerical simulations of the coupled Gross-Pitaevskii equation for binary condensates. In the numerical experiments of the superfluid system, the exponent in the microscopic regime anomalously reaches to its theoretical upper limit of the general scaling theory. The anomaly comes from the quantum-fluid effect in the presence of circular vortex sheets, described by the hydrodynamic approximation neglecting the fluid compressibility. It is also found that the distribution of superfluid circulation along vortex sheets obeys a dynamic-scaling law with different power-law exponents in the two regimes. An analogy to quantum turbulence on the hierarchy of vorticity distribution and the applicability to chiral superfluid $^{3}\mathrm{He}$ in a slab are also discussed.

Multiply quantized and fractional skyrmions in a binary dipolar Bose-Einstein condensate under rotation

We consider a binary dipolar Bose-Einstein condensate with repulsive contact and dipolar interactions under rotation. Our results show that the interplay among short-range interaction, long-range interaction, and rotation can give rise to a rich variety of topological configurations, including giant skyrmions with multiply topological charges and skyrmion-vortex lattices. In particular, we find that for fixed rotation frequencies, tuning the short- and long-range interactions can derive novel ground-state phases, such as a meron pair composed of two fractional skyrmions and a skyrmion with topological charge $Q=2$ centered in giant skyrmions.

Characteristic temperature for the immiscible-miscible transition of binary condensates in optical lattices

We study a two-species Bose-Einstein condensates confined in quasi-two-dimensional (quasi-2D) optical lattices at finite temperatures, employing the Hartree-Fock-Bogoliubov theory with the Popov approximation. We examine the role of thermal fluctuations on the ground-state density distributions, and the quasiparticle mode evolution. At zero temperature, the geometry of the ground-state in the immiscible domain is side-by-side. Our results show that the thermal fluctuations enhance the miscibility of the condensates, and at a characteristic temperature the system becomes miscible with rotationally symmetric overlapping density profiles. This immiscible-miscible transition is accompanied by a discontinuity in the excitation spectrum, and the low-lying quasiparticle modes such as slosh mode becomes degenerate at the characteristic temperature.

Two-dimensional symbiotic solitons and vortices in binary condensates with attractive cross-species interaction

We consider a two-dimensional (2D) two-component spinor system with cubic attraction between the components and intra-species self-repulsion, which may be realized in atomic Bose-Einstein condensates, as well as in a quasi-equilibrium condensate of microcavity polaritons. Including a 2D spatially periodic potential, which is necessary for the stabilization of the system against the critical collapse, we use detailed numerical calculations and an analytical variational approximation (VA) to predict the existence and stability of several types of 2D symbiotic solitons in the spinor system. Stability ranges are found for symmetric and asymmetric symbiotic fundamental solitons and vortices, including hidden-vorticity (HV) modes, with opposite vorticities in the two components. The VA produces exceptionally accurate predictions for the fundamental solitons and vortices. The fundamental solitons, both symmetric and asymmetric ones, are completely stable, in either case when they exist as gap solitons or regular ones. The symmetric and asymmetric vortices are stable if the inter-component attraction is stronger than the intra-species repulsion, while the HV modes have their stability region in the opposite case.

Faraday and resonant waves in binary collisionally-inhomogeneous Bose–Einstein condensates

We study Faraday and resonant waves in two-component quasi-one-dimensional (cigar-shaped) collisionally inhomogeneous Bose–Einstein condensates subject to periodic modulation of the radial confinement. We show by means of extensive numerical simulations that, as the system exhibits stronger spatially-localised binary collisions (whose scattering length is taken for convenience to be of Gaussian form), the system becomes effectively a linear one. In other words, as the scattering length approaches a delta-function, we observe that the two nonlinear configurations typical for binary cigar-shaped condensates, namely the segregated and the symbiotic one, turn into two overlapping Gaussian wave functions typical for linear systems, and that the instability onset times of the Faraday and resonant waves become longer. Moreover, our numerical simulations show that the spatial period of the excited waves (either resonant or Faraday ones) decreases as the inhomogeneity becomes stronger. Our results also demonstrate that the topology of the ground state impacts the dynamics of the ensuing density waves, and that the instability onset times of Faraday and resonant waves, for a given level of inhomogeneity in the two-body interactions, depend on whether the initial configuration is segregated or symbiotic.