Shell-shaped quantum droplet in a three-component ultracold Bose gas

Abstract

We present a scheme to generate shell-shaped droplet in a three-component (1,2,3) ultracold Bose gas. Here binary mixtures (1,2) and (2,3) form quantum droplets due to inter-species attractions, and the two droplets are mutually immiscible due to strong 1-3 repulsion. Importantly, the shared component-2 serves as a glue to link the two droplets together as a globally self-bound object. In this system, the outer droplet naturally develops a shell structure, and its radius and width can be conveniently tuned through the size of core droplet. Moreover, to reach an equilibrium with the shell, the core droplet displays very different spin densities as compared to the vacuum case. These results have been demonstrated in a realistic ${}^{23}$Na-${}^{39}$K-${}^{41}$K mixture. Our scheme liberates the shell-shaped Bose gas from stringent conditions with microgravity or fine-tuned traps, and can be readily implemented in cold atoms laboratories on Earth. This paves the way for future exploration of quantum droplets in curved space with non-trivial real-space topologies.

Quantum systems in curved space exhibit many distinctive features due to their non-trivial real-space topologies. For instance, a periodic boundary allows a persistent superflow of toroidal Bose-Einstein condensates (BECs) ring_1 ; ring_2 , and a local curvature gives rise to new topological defectsNelson ; Ho ; Zhou , interesting few-body physicsZhang_Ho ; Shi_Zhai and even non-Hermitian phenomenaZhouQi . As an outstanding case with non-trivial real-space topology, the shell-shaped BEC has attracted a great amount of attention in the field of ultracold atomsreview , and various fascinating properties have been revealed in terms of thermodynamicsTononi ; Rhyno , collective modesLannert ; Padavic ; Sun , expansion dynamicsTononi2 and vortex formationPadavic2 ; Fetter . These studies have been largely stimulated by the experimental plan of creating shell-shaped BECs, or atomic bubbles, using the shell-shaped potentials under radio-frequency dressingproposal_rf_dress_1 ; proposal_rf_dress_2 . However, the Earth’s gravity prevents the formation of a closed shell in such setupColombe , and consequently a microgravity environment is required. Indeed, the first shell-shaped BEC was recently realized in NASA Cold Atom Laboratory aboard the International Space StationNASA_expt . Shortly after that, it was proposed alternatively that binary bosons in immiscible regimeHo_Shenoy ; Pu_Bigelow can also achieve the shell structureWolf . This idea has been successfully implemented in a recent experiment in the presence of Earth’s gravityDajun_expt , where a magic-wavelength optical trap was applied to avoid different gravitational sags between two species.

Figure 1: (a) Interaction model for realizing shell-shaped geometry in three-component(’1’, ’2’, ’3’) bosons. Here (1,2) and (2,3) both form quantum droplet due to inter-species attraction, while the two droplets are immiscible due to 1-3 repulsion. The outer droplet is then shell-shaped. The common component-2 links the whole system together as a self-bound object. (b) Typical half-sphere distribution and according density profile along the radius direction. Here we consider a realistic ${}^{23}$Na-${}^{39}$K-${}^{41}$K (’1’-’2’-’3’) mixture near $B\sim 150$G with $a_{23}=-70a_{0}$ (see text). The atom numbers are $(N_{1},N_{2},N_{3})/10^{5}=(5,10.7,1.5)$.

In this work, we introduce a new type of shell-shaped BEC without resorting to any trapping potential, and therefore it does not rely on microgravity environment or fine-tuned traps as in previous experimentsDajun_expt ; NASA_expt . Our scheme is motivated by the recent development of quantum droplet in ultracold atoms, which has been realized in both dipolar gasesferrier-barbut_observation_2016 ; schmitt_self-bound_2016 ; ferrier-barbut_liquid_2016 ; chomaz_quantum-fluctuation-driven_2016 ; tanzi_observation_2019 ; bottcher_transient_2019 ; Chomaz_long_lived_2019 and boson mixturescabrera_quantum_2018 ; cheiney_bright_2018 ; semeghini_self-bound_2018 ; derrico_observation_2019 ; burchianti_dual-species_2020 ; guo_lee-huang-yang_2021 . These droplets are stabilized by a mean-field attraction and the Lee-Huang-Yang repulsion from quantum fluctuationsPetrov , and thus can be self-bound in vacuum. Here we remark that the self-bound nature of quantum droplet offers an ideal opportunity for creating perfect shell geometry on Earth. This is because in the absence of any external trap, different species in the droplet will fall freely with the same speed and therefore no relative displacement will be produced due to the gravity.

Our scheme of creating shell-shaped quantum droplet is illustrated in Fig.1. Specifically, we consider a three-component (1,2,3) boson mixture with contact interactions. Here (1,2) and (2,3) both form quantum droplets due to inter-species attractions, and the two droplets are immiscible due to strong 1-3 repulsion. Importantly, the shared component-2 acts as a glue to link two droplets together as a globally self-bound object. In this way, the outer droplet is repelled by the inner/core one and naturally develops a shell structure in free space. By increasing the core size, we find the outer shell can efficiently expand to larger radius and become gradually thinner. Moreover, to reach an equilibrium with the outer shell, the core droplet exhibits equilibrium densities that are very different from those in vacuum (without shell). As a result, the spin densities of the core are highly tunable, rather than locked to a specific ratio as in vacuum. We have demonstrated these results using a realistic mixture of ${}^{23}$Na-${}^{39}$K-${}^{41}$K near $B\sim 150$G. The shell-shaped droplet presented here opens a new avenue for studying ultracold gases in curved space, where quantum fluctuations and non-trivial real-space topologies can both play important roles and lead to many intriguing phenomena.

We write down the Hamiltonian of three-component boson mixtures $H=\int d{\bf r}H({\bf r})$, with ($\hbar=1$)

$$H({\bf r})=\sum_{i=1}^{3}\phi_{i}^{\dagger}({\bf r})\left(-\frac{\nabla^{2}}{2m_{i}}\right)\phi_{i}({\bf r})+\sum_{ij}\frac{g_{ij}}{2}\phi_{i}^{\dagger}\phi_{j}^{\dagger}\phi_{j}\phi_{i}({\bf r}).$$ (1)

Here ${\bf r}$ is the coordinate; $m_{i}$ and $\phi_{i}$ are the mass and field operator of boson species $i$, respectively; $g_{ij}=2\pi a_{ij}/m_{ij}$ is the coupling constant between species $i$ and $j$, with scattering length $a_{ij}$ and reduced mass $m_{ij}=m_{i}m_{j}/(m_{i}+m_{j})$.

For a homogeneous dilute gas with densities $\{n_{i}\}$ $(i=1,2,3)$, the total energy density is composed by the mean-field part $\epsilon_{\rm mf}=\frac{1}{2}\sum_{ij}g_{ij}n_{i}n_{j}$ and a correction from quantum fluctuations:

$$\epsilon_{\rm qf}=\int\frac{d^{3}{\bf k}}{2(2\pi)^{3}}\left[\sum_{i}(E_{i{\bf k}}-\epsilon_{i{\bf k}}-g_{ii}n_{i})+\sum_{ij}\frac{2m_{ij}g_{ij}^{2}n_{i}n_{j}}{{\bf k}^{2}}\right]$$ (2)

with $\epsilon_{i{\bf k}}=k^{2}/2m_{i}$, and $E_{i{\bf k}}$ the $i$-th Bogoliubov modema_borromean_2021 . For a general system with inhomogeneous densities, we employ an extended Gross-Pitaevskii(GP) equations incorporating quantum fluctuations:

$$\begin{split}i\partial_{t}\phi_{i}=\left(-\frac{\nabla^{2}}{2m_{i}}+\sum_{j}g_{ij}n_{j}+\frac{\partial\epsilon_{\rm qf}}{\partial n_{i}}\right)\phi_{i}.\end{split}$$ (3)

The ground state can be obtained by the imaginary time evolution of above coupled equations.

Different from the miscible three-component droplet in Ref.ma_borromean_2021 , in this work we focus on the situation where components 1 and 3 are immiscible with strong repulsion, while (1,2) and (2,3) are still miscible and can form binary droplets, see Fig.1(a). To achieve this in practice, we consider the boson mixture ${}^{23}$Na-${}^{39}$K-${}^{41}$K all at hyperfine state $|F=1,m_{F}=-1\rangle$ (denoted as 1-2-3). Near $B\sim 150{\rm G}$, we have $(a_{11},a_{22},a_{33},a_{12},a_{13})=(52,30,63,-50,213)a_{0}$ ($a_{0}$ is the Bohr radius)lysebo_feshbach_2010 ; viel_feshbach_2016 ; Schulze_2018 , and $a_{23}$ is highly tunable via a Feshbach resonance at $B_{0}=149.8$G with width $\Delta B\sim 25$mGtanzi_feshbach_2018 . Under above parameters, one can easily check that the required mean-field instabilities, i.e., collapse for 1-2 and 2-3 and phase separation for 1-3, will occur as long as $a_{23}$ is sufficiently attractive.

Our numerical simulations based on (3) indeed produce two mutually immiscible droplets, (1,2) and (2,3), as the ground state for above system, see typical density distribution in Fig.1(b). Remarkably, the whole system is still self-bound due to the presence of component-2. In such an immiscible phase, in principle each droplet can either stay inside as a core or outside as a shell. However, for the ground state we find the one with a higher density of 2 always stays inside and the lower one stays outside, i.e., $n_{2}^{\rm core}>n_{2}^{\rm shell}$. This can be attributed to the minimized surface energy under this configuration. For the parameter of $a_{23}$ chosen in this work, the core and shell droplets are respectively (2,3) and (1,2).

We note that similar phenomena of immiscible droplets were also found previously in Helium mixturesPi1999 ; Barranco2006 , where a heavy ${}^{4}$He droplet was coated with a normal ${}^{3}$He liquid, and in dipole-dipole mixtures with anisotropic density profilesbisset_quantum_2021 ; smith_quantum_2021 . In all these studies, a long-range attraction between two droplets is required for their self-binding. This is drastically different from the binding mechanism in our case, which does not rely on any long-range force.

Now we explore the equilibrium expansion of shell (1,2) droplet by increasing the size of core (2,3) droplet. Let us start from the situation with a small core and a thick shell, see Fig.2(a1), where both droplets display flat-top densities. As increasing the size (or atom number) of the core, the shell is repelled to larger radius and becomes gradually thinner, see Fig.2(a2,a3); meanwhile, its flat-top profile gradually disappears and gives way to a Gaussian distribution with quite narrow width and low density. In Fig.2(b), we further extract the maximal density of each component during this process. One can see that the shell densities $\{n_{2}^{\rm shell},n_{1}^{\rm shell}\}$ decrease rapidly as the core size increases, and eventually a very thin and dilute shell can be created.

Figure 2: Expansion of shell (1,2) droplet as increasing the size of core (2,3) droplet. Here we consider the ${}^{23}$Na-${}^{39}$K-${}^{41}$K (’1’-’2’-’3’) mixture at $a_{23}=-200a_{0}$. (a1,a2,a3) Density profiles of ground states for different atom numbers $(N_{1},N_{2},N_{3})/10^{5}=(1,1.73,0.01)$ (a1), $(1,1.86,0.1)$ (a2) and $(1,2.42,0.5)$ (a3). (b) Maximal densities of the core ($n^{\rm core}_{2},n^{\rm core}_{3}$) and shell ($n^{\rm shell}_{1},n_{2}^{\rm shell}$) as functions of $N_{\rm 3}$. The dashed horizontal line with according color denotes the equilibrium density for the shell ($n_{i}^{\rm shell}=n_{i}^{(0)}$) or for the core ($n_{i}^{\rm core}$).

Because of the shared component-2 as a link, the core and shell droplets, although spatially separated, are strongly correlated with each other. This is clearly reflected in the balance condition as derived below, which leads to very different equilibrium densities for the core if the shell is present. To start with, let us first consider the shell. Since it can transfer atom or energy to vacuum, its equilibration is, similar to an isolated droplet in vacuumPetrov , determined by zero pressure and minimized energy:

$$P_{\rm shell}=0;\ \ \ \epsilon_{\rm shell}=\epsilon_{\rm min}.$$ (4)

In the thermodynamic limit, this gives the same equilibrium density $n_{i}^{\rm shell}=n_{i}^{(0)}$ as the vacuum droplet, as well as a locked density ratio $n_{i}^{\rm shell}/n_{j}^{\rm shell}=n_{i}^{(0)}/n_{j}^{(0)}=\sqrt{g_{jj}/g_{ii}}$. However, the condition (4) does not apply to the core, since it can only transfer atoms to the shell but not directly to vacuum. In this case, the core should have the same pressure and chemical potential as the shell. In view of the presence of component-2 in both droplets, we arrive at the following balance condition:

$$P_{\rm core}=0;\ \ \ \ \mu^{\rm core}_{2}=\mu^{\rm shell}_{2}.$$ (5)

Two remarks are in order for the second condition in (5). First, it obviously differs from the second condition in (4), and therefore a different equilibrium density can be resulted for the core as compared to the vacuum case. Secondly, it exactly expresses the correlation between core and shell droplets. Under this condition, the core densities ($\{n^{\rm core}_{i}\}$) depend crucially on the shell chemical potential ($\mu^{\rm shell}_{2}$). As a result, $\{n^{\rm core}_{i}\}$ can be highly tunable by the shell parameters but no longer locked to a specific ratio as in vacuum case.

We now analytically derive $\{n^{\rm core}_{i}\}$ based on (4,5). Assuming a thermodynamically large ($i,j$) droplet with uniform densities $\{n_{i},n_{j}\}$, its total energy density can be written as

$$\epsilon=\frac{1}{2}\left(g_{ii}n_{i}^{2}+g_{jj}n_{j}^{2}\right)+g_{ij}n_{i}n_{j}+\epsilon_{\rm qf}.$$ (6)

Here $\epsilon_{\rm qf}$ is the correction from quantum fluctuationsPetrov :

$$\epsilon_{\rm qf}=\frac{8m_{i}^{3/2}}{15\pi^{2}}(g_{ii}n_{i})^{5/2}f\left(\alpha,\beta,\gamma\right)$$ (7)

where $f>0$ is a dimensionless function of

$$\alpha\equiv\frac{m_{j}}{m_{i}},\ \ \beta\equiv\frac{g_{ij}^{2}}{g_{ii}g_{jj}},\ \ \gamma\equiv\frac{g_{jj}n_{j}}{g_{ii}n_{i}}.$$ (8)

The chemical potential $\mu_{j}=\partial\epsilon/\partial n_{j}$ is then

$$\mu_{j}=g_{ij}n_{i}+g_{jj}n_{j}+\frac{8\left(m_{i}g_{ii}n_{i}\right)^{3/2}}{15\pi^{2}}\frac{\partial f}{\partial\gamma}g_{jj},$$ (9)

and the pressure $P=\mu_{i}n_{i}+\mu_{j}n_{j}-\epsilon$ can be obtained as

$$P=\frac{1}{2}\left(g_{ii}n_{i}^{2}+g_{jj}n_{j}^{2}\right)+g_{ij}n_{i}n_{j}+\frac{3}{2}\epsilon_{\rm qf}.$$ (10)

Assuming $n_{j}/n_{i}=c_{ji}$, the $P=0$ condition leads to the equilibrium density

$$n_{i}=\frac{25\pi}{1024a_{ii}^{3}}\left(\frac{\frac{1}{2}\left(g_{jj}c_{ji}^{2}+g_{ii}\right)+g_{ij}c_{ji}}{g_{ii}f}\right)^{2}.$$ (11)

In this way, both $\mu_{j}$ (Eq.9) and $\epsilon$ (Eq.6) can be expressed as functions of a single unknown parameter $c_{ji}$, which is to be determined by another condition apart from $P=0$.

For the shell (1,2) droplet, the second condition in (5), i.e., the energy minimization, results in a locked density ratio $c^{\rm shell}_{21}\approx\sqrt{g_{11}/g_{22}}$ as in the vacuum casePetrov . The densities $n^{\rm shell}_{1}$ and $n^{\rm shell}_{2}$ can be further obtained via (11), which are the same as $n^{(0)}_{1,2}$ for the vacuum droplet. The chemical potential $\mu_{2}^{\rm shell}$ can be obtained via (9) as

$$\mu_{2}^{\rm shell}=(g_{12}+\sqrt{g_{11}g_{22}})n_{1}(1-\frac{2}{3}\frac{\gamma}{f}\frac{\partial f}{\partial\gamma}).$$ (12)

For the core (2,3) droplet, by enforcing the second condition $\mu_{2}^{\rm core}=\mu_{2}^{\rm shell}$ in (5), and recalling $\mu_{2}^{\rm core}$ is just parametrized by the density ratio $c^{\rm core}_{23}$, we can then solve $c^{\rm core}_{23}$, and finally the densities $n^{\rm core}_{2}$ and $n^{\rm core}_{3}$ can be obtained via (11). In Fig.2 (b), $\{n^{\rm shell}_{1},n^{\rm shell}_{2}\}$ and $\{n^{\rm core}_{2},n^{\rm core}_{3}\}$ are respectively shown by horizontal dashed and dotted lines. These analytical results fit well to the actual densities extracted from our numerics when the core and the shell are both in thermodynamical limit.

Figure 3: Tunable equilibrium densities of the core (2,3) droplet in ${}^{23}$Na-${}^{39}$K-${}^{41}$K (’1’-’2’-’3’) mixture. (a) Equilibrium densities $n_{2}^{\rm core}$ and $n_{3}^{\rm core}$ as functions of $a_{23}$, as compared to $n_{2}^{(0)}$ and $n_{3}^{(0)}$ in vacuum. (b) Density ratio $c^{\rm core}_{23}=n_{2}^{\rm core}/n_{3}^{\rm core}$ and its vacuum counterpart $c^{(0)}_{23}=n_{2}^{(0)}/n_{3}^{(0)}$ as functions of $a_{23}$. The stars ($\star$) show numerical results by simulating the GP equations (3) in imaginary time for large atom number $N_{i}=10^{5}\sim 10^{6}$. Inset shows $\mu_{2}^{(0)}$ for isolated (2,3) and (1,2) droplets in vacuum.

In Fig.3 (a,b), we further show how the core densities $\{n_{2}^{\rm core},n_{3}^{\rm core}\}$ and its ratio $c^{\rm core}_{23}$ vary with $a_{23}$, in comparison to the same quantities for a vacuum droplet. One can see that both $n_{2}^{\rm core}$ and $n_{3}^{\rm core}$ change dramatically with $a_{23}$, and they both deviate visibly from the vacuum values $n_{2}^{(0)}$ and $n_{3}^{(0)}$. Their ratio, $c^{\rm core}_{23}$, also shows a sensitive dependence on $a_{23}$, in distinct contrast to the locked $c^{(0)}_{23}\sim\sqrt{g_{33}/g_{22}}$ in vacuum. Moreover, we see that the deviations between $c^{\rm core}_{23}$ and $c^{(0)}_{23}$ become more pronounced if $a_{23}$ gets more attractive. This can be attributed to the larger mismatch of $\mu_{2}^{(0)}$ for isolated (2,3) and (1,2) droplets in vacuum, see the inset of Fig.3(b). As a result, to balance the two chemical potentials as required by (5), the (2,3) core has to adjust its densities to change $\mu_{2}^{\rm core}$, such that it can match $\mu_{2}^{(0)}$ in the (1,2) shell. This is how the correlation is built up between core and shell, and the modified core densities are a direct evidence of such correlation. In this way, the core-shell structure presents a rare situation where the equilibrium density of a quantum droplet can be efficiently tuned by its environment.

In summary, we have demonstrated a scheme to create shell-shaped quantum droplets in three-component ultracold boson mixtures. Different from existing experiments on shell-shaped BECNASA_expt ; Dajun_expt , this scheme does not rely on any trapping potential and therefore can be readily implemented in cold atoms laboratories on Earth. In addition, we have revealed the equilibrium expansion of the shell, the core-shell balance condition, as well as the high tunability of core densities. All these results can be practically detected in ${}^{23}$Na-${}^{39}$K-${}^{41}$K boson mixture near $B\sim 150$G. Other atomic candidates may also be possible in view of a growing number of Bose mixtures nowadays available in cold atoms laboratories with tunable interactionsK-Rb1 ; K-Rb ; Na-Rb ; K-Cs ; Schulze_2018 ; tanzi_feshbach_2018 .

Our work opens up a new avenue to study quantum droplets in curved space, where quantum fluctuations and non-trivial real-space topologies both play important roles and their interplay can lead to many intriguing phenomena. For instance, it would be interesting to investigate the collective excitation and vortex formation of shell droplet. Compared to the corresponding properties of binary droplets in flat geometriesPetrov ; Li2018:TwodimensionalVortexQuantum ; Zhang2019:SemidiscreteQuantumDroplets ; Tengstrand2019:RotatingBinaryBoseEinstein ; Caldara2022:VorticesQuantumDroplets ; Gu ; Oktel , the shell structure is expected to bring substantially new physics due to its compactness and local curvature. Moreover, dynamical expansion of core and shell droplets is another interesting subject to explore, which, for example, can be driven by suddenly quenching the coupling parameters of core or shell (such as $a_{23}$ in our model system). Indeed, the quantum droplet has been shown to exhibit unique dynamical properties that are absent in gaseous statecollision1 ; collision2 ; dynamic1 ; dynamic2 ; dynamic3 ; Ma2 . Furthermore, it is worthwhile to study quantum fluctuation effect when the shell become thin enough, in which limit the fluctuation energy in Eq.2 is no longer applicable since the system behaves as effectively 2D. Given that a 2D droplet can be supported at an arbitrarily low densityPetrov2 , we expect the actual shell can extend to very large radius with extremely thin width. Finally, it is noted that the present scheme of creating a shell droplet in 3D can be directly generalized to quasi-2D system in creating a toroidal (ring) droplet. We hope our current work could stimulate more studies in future on these fascinating curved-space droplets.