Blog

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Scientific Reports volume 12, Article number: 18467 (2022) Cite this article
322 Accesses
2 Altmetric
Metrics details
A theoretical analysis of binary collisions of quantum droplets under feasible experimental conditions is reported. Droplets formed from degenerate dilute Bose gases made up from binary mixtures of ultracold atoms are considered. Reliable expressions for the surface tension of the droplets are introduced based on a study of low energy excitations of their ground state within the random phase approximation. Their relevance is evaluated considering an estimation of the expected excitation energy having in mind the Thouless variational theorem. The surface tension expressions allow calculating the Weber number of the droplets involved in the collisions. Several regimes on the outcomes of the binary frontal collisions that range from the coalescence of the quantum droplets to their disintegration into smaller droplets are identified. Atoms losses of the droplets derived from self-evaporation and three-body scattering are quantified for both homo- and hetero-nuclear mixtures. Their control is mandatory for the observation of some interesting effects arising from droplets collisions.
In most circumstances, degenerate states of dilute atomic gases are realized using samples trapped by different means. The question of whether they can also exist in isolation, without a trapping potential, was theoretically answered in the affirmative within two scenarios. The first one corresponds to a homogeneous atomic gas where the two-body scattering length is negative and much larger than the effective two-body interaction radius; then the contribution to the ground state energy due to the three-body correlations could become dominant and provide a self-trapping mechanism¹. The second scenario corresponds to a condensed Bose–Bose mixture where the interspecies attraction becomes stronger than the geometrical average of the intra-species repulsions². The stabilization mechanism is then provided by quantum fluctuations that are a direct manifestation of beyond mean-field effects³. The implementation of the latter mechanism was the basis for the experimental observation of self-bound droplets of ultra cold atoms in homo-nuclear^4,5 and hetero-nuclear⁶ mixtures of alkali atoms, and in dipolar condensates^7,8,9,10 that required driving the Bose-Einstein condensate (BEC) into the strongly interacting regime.
Collisions of quantum droplets is the only mechanism available in current experiments that allows the direct observation of surface excitations. Up to now, this mechanism has been used to study the crossover between compressible and incompressible regimes¹¹. It is expected that in the incompressible regime surface waves occur and surface tension should be a relevant parameter to understand the conditions required for different results of quantum droplet collisions.
Following the ideas implemented by the nuclear physics community to study the expectations of the nuclear liquid drop model^12,13,14, in this work we introduce expressions for the low energy surface excitations of a binary mixture of Bose–Bose ultra cold atoms. These expressions are taken as variational ansatz that include effects directly due to the thickness of the surface of the droplet and its symmetries.
Since we are interested in excitations from the ground state of the droplet, the variational ansatz depends on that function. We evaluate numerically the ground state for a finite number of atoms, ranging from (10^4)–(10^7) within the extended Gross–Pitaevskii equation (EGPE) that incorporates the Lee–Huang–Yang term (LHY)³ in three dimensions. The analysis is performed for binary mixtures of homo- and hetero-nuclear atomic Bose gases that give rise to quantum droplets. In the case of homo-nuclear BEC mixtures, we also obtain the ground state within an effective formalism that considers finite-range effects as developed in Refs^15,16. General features of the resulting ground states are then evaluated both numerically and using analytical expressions built from numerical ground state functions, since the corresponding densities support a fit by a Boltzmann profile. These expressions allow a direct calculation of the energy of the surface excitations derived from the different ansatz. The resulting energies, besides being compared to previous phenomenological expectations, permit both the selection of the most adequate ansatz to describe the excitations according to Thouless variational theorem, and the evaluation of the corresponding surface tension.
Atom losses due to self-evaporation—directly dependent on the excitation modes- and three-body recombination effects are analyzed within the EGPE scheme. The time evolution of such losses for experimental relevant conditions is evaluated (Supplementary Material). Finally, the outcomes of frontal binary collisions are explored. Similarly to previous experimental and theoretical studies^11,18, we find out several regimes for binary front collisions that range from the coalescence of the quantum droplets to their disintegration into smaller droplets. For classical liquid droplets their collisions are usually investigated in terms of the impact parameter and of the Weber number. The latter is a measure of the relative importance of the inertia of a fluid in terms of the kinetic energy compared to its surface tension¹⁹. In this work, collisions are studied in terms of a generalization of the Weber number (mathrm {We}) to the quantum regime. A direct correlation between (mathrm {We}) and the collision outcomes is found both for homo- and hetero-nuclear Bose mixtures and in the incompressible and compressible regimes.
We are interested in obtaining expressions for the surface tension of a quantum droplet. To that end we shall follow the formalism developed by Berstch in the study of capillary waves in a superfluid¹², and which was later extended to Fermi systems with spherical symmetry¹³ having in mind the liquid drop model of nuclei. This formalism is based on the random phase approximation (RPA) and the Thouless variational theorem²⁰. RPA provides a strategy to decouple the collective motion of a many-body system from their individual motion²¹. In our case, RPA is equivalent to the linearization of the time-dependent equations that define the dynamics of the system in the vicinity of a variational minimum. A contemporary view of the theoretical basis that supports the Thouless theorem and its relevance for many-body theory is revisited in Ref.²². To study quantum droplets arising from vacuum fluctuations in a binary mixture of Bose gases, it is useful to implement Thouless variational principle to such a mixture.
In an analysis of the collective modes of nuclei in RPA, Thouless²⁰ showed that it is always possible to find a solution of the independent particle model that makes these collective modes as stable, i.e., they have real eigenfrequencies. This occurs because the solution that gives the lowest expectation value of the Hamiltonian also stabilizes the collective modes in RPA. The stability condition allows the construction of the variational principle and assures that the eigenvectors form a complete set in most cases; this assertion requires the proper definition of a scalar product. Thouless variational theorem guarantees that actual excitation energies are a lower bound to the excitation energies that result from using, as an ansatz, known expressions of excitations with harmonic time dependence and satisfying adequate boundary conditions.
Consider a binary mixture of Bose gases composed by (N^{(a)}) and (N^{(b)}) atoms of each species. Let us assume that the state of the system is approximately described by the Hartree many-body wavefunction
As order parameters we consider
(psi ^{(alpha )}) evolves according to the equation
In Eq. (3a), (V^alpha _{ext}(vec{r})) is the external potential and (rho ^{(alpha )}=psi ^{(alpha )*}psi ^{(alpha )}) the density of (alpha )-atoms. The real functionals ({mathcal {U}}_{alpha alpha }) and ({mathcal {V}}_{alpha beta }) in the integral operator (hat{U}^alpha ) may incorporate, in an effective way, interactions beyond the standard mean field approximation. The functional ({mathcal {V}}_{alpha beta }) represents the interaction between different species, i.e. (alpha ne beta ). In the case of EPGE with LHY, the density dependent interactions are superpositions of contact terms²,
The stationary ground state solutions of Eq. (3) define the chemical potentials (mu _alpha ),
Let us consider collective Bogolubov excitations with a harmonic time dependence,
In this equation q is a label that determine the characteristics of the excitation derived, e. g., from its geometry. The excitation functions belong to the space expanded by a normalized and complete basis set ({u_q^a,v_q^a,u_q^b,u_q^b}), according to the scalar product
Excitations with (omega _q = 0) are assumed to be orthogonal to all other excitations with a different set of quantum numbers.
If (phi ^{(alpha )}_0) is a real function, the linear response approximation around the stationary function (psi ^{(alpha )}_0 (vec{r},t)) is equivalent to the RPA approximation. Taking into account the real character of functionals ({mathcal {U}}_{alpha alpha }) and ({mathcal {V}}_{alpha beta }), and factorizing the terms that are either proportional to (e^{ iomega _q t}) or (e^{-iomega _q t}), Eq. (3) in the RPA approximation is equivalent to the equations
(delta rho ^{(alpha )}_{q;eff}= phi ^{(alpha )}_0eta _q^{(alpha )}) could be directly interpreted as an the effective variation of the density for dynamical purposes since it multiplies the variational derivatives of ({mathcal {U}}_{alpha alpha }) and ({mathcal {V}}_{alpha beta }). In fact, using ({hat{Delta }}_0^alpha phi ^{(alpha )}_0=0), the continuity equation
is derived; it implies a conservation relation for both the real and imaginary parts of (int d^3x delta rho ^{(alpha )}_{q;eff}e^{-iomega _q t}).
Thouless variational theorem²⁰ is now applied. In the system under consideration, the application of this theorem is based on the stability inherited to (psi ^{(alpha )}_0 (vec{r},t)) from the stationary character of the ground state solutions of Eq. (3). The energy (hbar omega _q) is estimated after from Eq. (8) solving for (omega _q), and averaging the estimates obtained from each atoms species, (alpha = a,b). The variational condition is then written as
The exact expressions for (u_q) and (v_q) would satisfy Eqs. (8) and they are unknown. In a variational calculation, an ansatz for the variational functions (eta _q^{(alpha )}) and (zeta _q^{(alpha )}), Eq. (8c), is done incorporating in as much as possible the expected physical characteristics of the particular system under consideration. Our ansatz for the excitation functions depends functionally on (phi _0) and its derivatives and takes as variational parameters a relative weight between (eta _q^{(alpha )}) and (zeta _q^{(alpha )}) for each atomic species (alpha =a,b),
By requiring that (gamma _alpha ) yields a minimum for (omega _{eta zeta }) it results
and the corresponding optimal variational excitation energy is then
Since quantum droplets exhibit different regimes, e. g. compressible and incompressible, we decided to explore the relevance of two pairs of excitation functionals (eta ^{(alpha )}_q ) and ({tilde{zeta }}^{(alpha )}_q ), both of them with the simple structure,
The functions (vec{F}^alpha (vec{r})) and ( Lambda ^alpha (vec{r})) incorporate geometric properties of the quantum fluid. The parameter (xi ^{(alpha )}) corresponds to a characteristic scale length for the (alpha )-species. For these expressions for (eta ^{(alpha )}_q) and (zeta ^{(alpha )}_q=gamma _alpha {tilde{zeta }}^{(alpha )}_q),
Thus, each component of ( vec{F}_q^alpha ) would determine the variations of the density (vert phi _0^{(alpha )}vert ^2) along different spatial directions, and the “velocity” field (vec{v}^{(alpha )}(vec{r})) is determined by the gradient of (Lambda _q^alpha ), i.e. (vec{v}^{(alpha )}(vec{r})propto vec{nabla }Lambda _q^alpha ).
If a sphere delimits the boundary of the (alpha )-fluid there are two alternative configurations. In the case the fluid is contained within the sphere (droplet) or outside it (bubble). For a droplet (Lambda ^alpha _q) is taken as the solution of Laplace equation in spherical coordinates that is regular at the origin
with ({mathcal {A}}_{ell m}) a normalization constant to be evaluated according to Eq. (7). Even though in the present study we focus on quantum droplets, the bubble configuration could result useful for the description of the mixed-bubble regime in which bubbles of the mixed phase coexist with a pure phase of one of the components²³. In such a case, one would take (Lambda _{ell m}^alpha (vec{r}) ={mathcal {A}}_{ell m} Y_{ell m}(theta ,varphi ) (r/xi ^{(alpha )})^{-(ell +1)}).
Concerning the function (vec{F}^alpha _{ell m}(vec{r})), two reasonable ansatz for spherical boundaries are identified,
Ansatz 1.
It is equivalent to an initially infinitesimal movement of the droplet surface in the radial direction with an amplitude proportional to the angular solution of the wave equation, i.e., a spherical harmonic. This ansatz resembles that discussed by Berstch¹² for a quantum fluid whose surface is defined by the (z=0) plane and for which (eta = partial _zphi _0 e^{ikx}).
Ansatz 2.
This ansatz is equivalent to (delta rho _{q;eff}^{(alpha )} sim vec{v}^{(alpha )}cdot nabla vert phi _0^{(alpha )}vert ^2sim dvert phi _0^{(alpha )}vert ^2/dt), i. e., the surface of the droplet is displaced along the velocity field. Note that the second ansatz privileges an hydrodynamic interpretation, Eq. (19), while the first ansatz at (t=0) yields an infinitesimal deformation of the droplet perpendicular to its surface, Eq. (18).
A direct calculation shows that, in case the dynamics is governed by contact interactions as it occurs for the EGPE, the contribution of the interaction terms of these excitation functions to A and (B_alpha ), Eq. (13), is null. However, information about such interaction terms is still contained in the stationary solution (phi _0^{(alpha )}). For quantum droplets that can be generated without an external potential, the ground state order parameters depend only on the radial coordinate and are proportional to each other,
the last equality being a consequence of Eq. (2). Under these conditions, the excitation energies become
Ansatz 1.
Ansatz 2.
whenever an equal length scale
is chosen for both species. In these equations we notice the presence of the total mass of the droplet
In his seminal work²⁴, Lord Rayleigh considered a classical spherical droplet with constant density of particles (rho _0) with mass m and a radius (R_0). For the frequency of excitation (omega _{exc}), he showed that
and he identified (sigma ) as the surface tension. Lord Rayleigh obtained such an structure for (omega _{exc}) analyzing the lowest energy surface excitations of a classical droplet via a variational theorem. The term
is directly related to the total mass of the droplet (M_T).
We observe that the dependence on the multipole parameter (ell ) both for ansatz 1 and 2, shows that monopole excitations are not surface excitations, and that the (ell =1) contribution is null. This is expected since this term corresponds to a translation of the whole droplet. However for (ell ge 2), the estimation value for (epsilon =hbar omega _{exc}) of the spherical quantum droplet, and as a consequence of the surface tension, is different for each ansatz:
For ansatz 1, the (ell )-dependence of the coefficient resembles that of Rayleigh and that of the liquid drop model of nuclei in the limit of (Nrightarrow infty ) nucleons²⁵.
For ansatz 2, the (ell )-dependence of the coefficient is similar to that predicted for atomic nuclei by Bertsch¹³ and by Casas and Stringari¹⁴ using RPA and the density-density Green’s function formalism; in both cases, the structure of the second ansatz was taken for the excitation modes.
A variational criteria based on the excitation energy can be applied to select between the two excitation ansatz. Once adequate expressions for the order parameters (phi ^{(a)}_0) are given, the approximate excitation energies are evaluated using both ansatz. The variational criteria establishes that the best approximation to the exact description of those excitations is that which provides the lowest excitation energies.
Following Lord Rayleigh ideas, the corresponding surface tension estimation would be given by the expressions
Ansatz 1.
Ansatz 2.
with
The density of a spherical liquid droplet conformed by N atoms is characterized by a radius (R_0), a skin width dR and a saturation density. The generic shape of such a droplet can be approximately described by a Boltzmann function,
For Bose–Bose mixtures of experimental interest and number of atoms ranging from (10^4)–(10^7), we have performed numerical calculations for the ground state derived from an EGPE for homo- and hetero-nuclear mixtures, and, for the homo-nuclear case, from the effective formalism that considers finite-range effects^15,16. The resulting density profiles fits the Boltzmann density with a reliability above 0.998 in all considered cases. The resulting parameters are reported in the Supplementary Material. Using Eq. (28), the surface tension Eqs. (26–27) can be written in terms of the Fermi-Dirac integrals
So that, for (ell > 1)
expressions that correspond to ansatz 1 and 2 respectively.
Consider the collision of quantum droplets that initially were radially spaced at a distance (2d_0); they are assumed to be kicked off as described by a plane wave wavefunction factor,
The mean value
is taken as a measure of the initial kinetic energy of the droplets. Assuming a non significant overlap between the droplets at (t=0), it can be decomposed as the sum of the translational kinetic energy of each droplet as a whole ({mathcal {K}}_{trans}), and an internal kinetic energy of the atoms within the droplets ({mathcal {K}}_{int}). Thus, ({mathcal {K}}approx {mathcal {K}}_{trans} + {mathcal {K}}_{int}), with
The dimensionless quantity suitable for characterizing binary collisions of classical incompressible droplets is the Weber number¹⁹. It is a dimesionless quantity for each one of the droplets, defined as the ratio between the translational kinetic energy ({mathcal {K}}_{trans}) before the collision and the surface energy of excitation evaluated in terms of the surface tension of the droplet. This definition can be extended to quantum droplets as
where (R_0) is the initial radius of the droplet. For the collision of two identical quantum droplets like those described in this work, this expression becomes
Ansatz 1
Ansatz 2
Notice that the expressions of (mathrm {We}_ell ) for the ground state of the quantum droplets given by Eqs. (36–37) make evident the relevance of the relation between the de Broglie wavelength and the skin width of the droplet, (k_0 dR), as a measure of the momentum that could be transfered to the droplet during the collision. The second relevant parameter is the ratio (R_0/dR) that determines the deformation of the droplet due to a multipole excitation.
The theoretical structure developed in the previous Section will be applied to quantum droplets of homo-nuclear ((^{39})K) and hetero-nuclear ( (^{41})K and (^{87})Rb) mixtures as those predicted by^2,5. The physical parameters of the illustrative examples in this and the following sections are compatible with experiments reported in Refs^4,6. We consider a mixture of two ultracold Bose gases with atom masses (m_a) and (m_b), densities (n_a) and (n_b), scattering lengths (mathrm {a}_{alpha beta }), (alpha ,beta =a,b), and coupling constants (g_{alpha beta }= 2pi hbar ^2mathrm {a}_{alpha beta }(1/m_a +1/m_b)) with (g_{aa}> 0), (g_{bb} > 0), and (g_{ab}<0). The first correction to the mean-field contribution to the ground state energy of a homogeneous weakly repulsive Bose gas is given by the LHY term³. This term is negligible in many circumstances. However, there are situations where the LHY and the mean field contributions can be of the same order without leaving the weakly interacting regime, (nmathrm {a}^3ll 1). Whenever (g^2_{ab}sim g_{aa} g_{bb}) the EGPE which incorporates the LHY term is⁶
The inter- and intra-species coupling constants (g_{alpha beta }) can be independently controlled. The mean-field energy term, which is proportional to the square of the density (n^2), that is negative and the LHY contribution, which is proportional to (n^{5/2}) is positive. Then, the quantum LHY repulsion could neutralize the mean-field attraction and stabilize the system against collapse². A less intuitive, albeit more precise description of the self-trapping mechanism can be obtained via Quantum Monte Carlo simulations¹⁵, they predict similar general properties for the quantum droplets, though different specific scenarios of their occurrence. The general properties are: generic features of the collective modes⁹, the presence of striped states in confined geometries²⁷, and the existence of arrays of phase coherent droplets with transient super solid properties^10,28,29.
The self-trapping regime requires a minimum number of atoms (N_c^{(alpha )}) which can be evaluated by calculating the ground state numerically. There is a transient regime characterized by a set of additional critical numbers (N_{ql}^{(alpha )}). For (N_{c}^{(alpha )}<N^{(alpha )}<N_{ql}^{(alpha )}) the atomic density exhibits a surface with a thickness similar to the radius of the droplet. If the quantum droplet has no collective modes with energy lower than the particle emission threshold self-evaporation can occur^2,30. For higher N values, the ground state corresponds to that expected for a liquid either a classical liquid, a helium nano droplet, or an atomic nuclei within the liquid drop model. Then, the self-trapping quantum state is characterized by a core with a quasi uniform density and a thin free surface. According to EGPE Ref.²: (1) A spherical shape with radius (R_0approx ( 3 N/4pi n_alpha ^{(0)}xi ^3)^{1/3}xi ), is expected with
and (delta g = g_{ab} + sqrt{g_{aa}g_{bb}}); (2) a surface thickness of order (xi ); (3) a saturation density in the limit of (N^{(alpha )}rightarrow infty )
As a consequence, the well-known shape of a Boltzmann function Eq. (28) is a feasible mathematical model of the spherical density profile of a droplet for a given number of particles with (R_0) and dR the radius and suface thickness and the density at the origin, (rho _B(0;N) = {{{mathcal {B}}}}_1/(1 + mathrm {exp}(-R_0/dR))). Details for specific homo-nuclear and hetero-nuclear BEC mixtures of experimental interest for a finite number of atoms are given in Section SMA of the Supplementary Material. It is worth to mention that the reliability of the Boltzmann density model was found to be above 0.998 in all numerical simulations.
Within the Bogolubov scheme, the excitation modes of the quantum droplets are determined by the normalized and complete basis set ({u_q^{(a)},v_q^{(a)},u_q^{(b)},u_q^{(b)}}). An ideal collective excitation of the wave function (psi ^{(alpha )}_{exc} ) shows a harmonic time dependence, Eq. (6).
Figure 1 illustrates the resulting ideal evolution of the atomic densities for hetero-nuclear mixtures applied to the ground state wave function (phi ^{(alpha )}_0) of the EGPE being (u_q^{(a)}) and (v_q^{(a)}) the functions determined by either Eq. (18) (ansatz1:A1) or Eq. (19) (ansatz2:A2). The number of atoms in the illustrative case of Fig. 1 corresponds approximately to the minimum required for the ground state droplet to achieve the saturation density. In the Supplementary Material SMB other ideal evolutions for cases of particular relevance in the compressible regime are presented. The unit of time is taken as suggested in Ref.²,
and the unit of energy is (hbar /tau ). Graphs are obtained from isodensity surfaces chosen to approximately delimit the droplets. Apparently the droplets seem to divide in several sub-droplets during particular time intervals. However this interpretation is not precise since it is just based in the isodensity surfaces. The more detailed numerical analysis shows that there is always a dilute gas between the apparent fragments. In Fig. 1 we observe clear differences and similarities between the structures of the excitation modes for different ansatz.
It is worth mentioning that these ideal evolutions do not take into account atomic losses, which in the experimental realm involve at least two processes: self-evaporation and three-body dispersion. They will be discussed in a following Section and in the Supplementary Material. The case illustrated in Fig. 1 is not expected to show significant self-evaporation. Note, however, that in this study the relevance of ideal droplet evolution under excitation modes is conditional on its ability to illustrate the structure of atomic densities under surface excitations as a result of quantum droplet collisions.
Comparative illustrations of the ideal evolution of the density of a quantum droplet formed by (N^{(Rb)}) = 96818 and (N^{(K)}) = 84578 atoms considering the excited quadrupole state, (ell =2) and (m_ell =0), generated from the ground state droplet. Ansatz 1 corresponds to rows (a–d) and (i–l) and ansatz 2 to rows (e–h) and (m–p). In this case (omega _1 =0.21tau ^{-1}) and (omega _2 =0.22tau ^{-1}). The 3D graphs where generated by plotting isodensity surfaces at fifteen percent of the maximum density at each time.
The ground state wave functions can be used to calculate the energies associated to surface excitations within a description where just contact interactions are considered, i.e., according to ansatz 1 (A1), Eq. (21) and ansatz 2 (A2), Eq. (22). The numerical results are illustrated in Fig. 2 for the physical parameters given in Section SMA of the Supplementary Material.
For homo-nuclear mixtures, the phenomenological expression
based in RPA is also illustrated. The variational criterium identifies ansatz 2 as a better evaluation of the energy for small values of N. However, for a given (ell ) value, there is a critical value (N_ell ^c) for which ansatz 1 gives a lower excitation energy if (N >N_ell ^c). The phenomenological estimation (epsilon _ell ^{(P)}) cannot be considered closer to the exact one even though it may be lower than the result of ansatz 1 or 2, both because it is not supported by a variational theorem and because it is expected to be valid just for (Ngg 1)^2,31; note that the ansatz 1 results approximately coincide with (epsilon _ell ^{(P)}) in this limit.
For low N calculations where ansatz 2 gives lower energies, droplets have a radius (R_0) similar to its width dR, thus the quantum incompressible regime has not been reached. The dependence on (ell ) of the energy involving the factor (ell (ell -1)(ell +2)) that was predicted by Lord Rayleigh for classical liquids, is also predicted by ansatz 1 and the phenomenological result Eq. (41), but not by ansatz 2. As illustrated in Fig. 2, the similarities between the results for the excitation energies obtained for LHY calculations and the effective equation resulting from density Monte Carlo calculations are remarkable. Notice however that there are relevant subtleties. In particular, for given N and (ell ) values the excitation energies of the suggested modes that surpass the absolute value of the chemical potential (mu ) differ. So that, N values for which self-evaporation occurs are similar but not identical to those predicted in Ref.². The differences arise not just from the dynamical equation but also from the natural excitation modes description.
Summarizing, for homo-nuclear mixtures using the variational criteria, both for LHY and Monte Carlo inspired calculations, it results that for (N<N^ell _c) ((N>N^ell _c)) low energy excitations are better described by ansatz 2 (ansatz 1). We have (N^{ell =2}_csim N_{ql}) with (N_{ql}) the critical value between compressible and incompressible regimes.
Approximate excitation energies (epsilon _ell ) as a function of the number of atoms N for (a–c) a homo–nuclear mixture of (^{39})K atoms; (d) a hetero-nuclear mixture of (^{41})K and (^{87})Rb atoms. The negative value of the chemical potential (mu ) is also shown. (epsilon _ell ) is evaluated according to different models: (a) a phenomenological expression expected to be valid just for (Ngg N_c)² (P); (b, d) the variational scheme for the LHY-EGPE, Eqs. (21–22); (c) the variational scheme for the effective interaction obtained from Monte-Carlo (MC)¹⁶. In (b–d) the results from ansatz 1, A1, (dashed) and ansatz 2, A2, (continous) are shown. The natural unit (hbar /tau ) of energy is used with (tau =1.308)  ms (homo-nuclear) and (tau =1.191) ms (hetero-nuclear).
In the hetero-nuclear case, the behaviour of the two chemical potentials is different for the two species as expected from their differences in mass and coupling interactions. In Fig. 2d, the values of the negative of the chemical potentials (-mu _K) and (-mu _{Rb}) of the ground states obtained from LHY-EGPE are shown. Notice that (vert mu _Kvert ) always increases as (N^{(K)}=sqrt{g_{RbRb}/g_{KK}} N^{(Rb)}) increases, while for (^{87})Rb atoms (vert mu _{Rb}vert ) at (N^{(Rb)}=N_c^{(Rb)}) is higher than its asymptotic value(~ 0.1hbar /tau ) for (N^{(Rb)}rightarrow infty ). The crossing point of those chemical potentials corresponds to (N^{(Rb)}-N_c^{(Rb)} approx 8^4), that is (N^{(Rb)}approx 16326approx 1.3cdot N_c^{(Rb)} ) which is smaller than the maximum saturation density. That is, for (N_c^{(alpha )}le N^{(alpha )}le 1.3N_c^{(alpha )}) it is energetically favorable to release a single (^{41})K atom from the fluid, than a (^{87})Rb atom; the reverse condition applies for higher (N^{(alpha )}) values. That is, the differences in the value of the chemical potential indicates that the self-evaporation rates of each species differ and depend on the number of atoms in the droplet. In a similar way than the homo-nuclear case, ansatz 2 gives a better estimation of the energy for the compressible regime at small values of the number of atoms, (N^{(Rb)}< 70227) for quadrupole excitations and (N^{(Rb)}< 10^7) for octupole excitations. Ansatz 1 must be taken as a better variational estimation for high (N^{(a)}) values. According to the variational criteria, the corresponding Bogolubov functions should be selected to describe the excited state.
In the incompressible regime, the direct numerical analysis of the Bogolubov de Gennes equations Eq. (8) reveals that the energies estimated with the variational criteria reported in this work coincide with the exact ones at the two digit level.
The experimental observation of the dynamics of quantum droplets before and after a collision is highly constrained by atom losses. Self-evaporation and three-body scattering have been identified as relevant atom losses mechanisms.
Self-evaporation arises whenever the number of atoms in the droplet is such that its excitation energy (epsilon ) lies in the continuum, i.e., (epsilon >-mu ), see Fig. 2. Previously, self-evaporation was studied²⁶ using a monopole excitation induced by an isotropic harmonic potential, only radially dependent. Here we benchmark the quadrupole excitation functions introduced in “A pair of variational ansatz for surface excitations” section, and the time evolution within the EGPE scheme of an excited droplet considering angular deformations. The initial state of the atomic cloud was taken as given by Bogolubov expression, Eq. (6) at (t=0). The evaporation rates are evaluated in terms of the atoms that leave the droplet as a function of time; the latter are identified as those located outside the minimum radius of the sphere that includes all the atoms in the ideal evolution of the approximate excitation mode given by the Bogolubov expression. In the Section SMB of the Supplementary Material, the evolution of self-evaporation is illustrated for homo-nuclear and hetero-nuclear mixtures and some interesting processes are highlighted.
For the homo-nuclear case, the biggest atom loss is due to three-body scattering with typical time scale in the 10 ms range. Notice that additional mechanisms that threaten the stability of a droplet arise from changes on the quotient (N^{(a)}/N^{(b)}) as well as the losing the spatial overlap of the two atomic wave functions that is required for self-trapping. Our numerical simulations confirm that for hetero-nuclear mixtures of (^{87})Rb and (^{41})K, there is a time window in the tens of milliseconds scale where the droplet state is preserved, as seen in Section SMC of the Supplementary Material. Therefore, collisions of two droplets could be observed. In fact, experiments of collisions in the compressible regime have already been performed, Ref.¹¹.
The observation of the frontal collision of two quantum droplets were first reported in Ref.¹¹. The corresponding experiments involved homo-nuclear mixtures of (^{39})K atoms in two different hyperfine states. The two droplets were generated from degenerate gases located nearby the minima of a double well potential. Subsequently, the potential barrier was turned off and the two droplets migrate towards the center of a harmonic confining potential. After a controlled time interval (Delta t) the potential was turned off, and the droplets collided in a potential free environment; (Delta t) was used to tune the initial relative speed of the droplets³². The analysis of experimental data in Ref.¹¹ in terms of the relative velocity was inspired by a classical Weber number definition, (mathrm {We} = 2rho R_0 v^2/sigma ). It was found that if the speed were higher than a critical value (v_c) as a result of the collision the droplets coalesced; (v_c) depends in a non monotonic way on the number of atoms in each initial droplet. This is a consequence of the dependence of (rho ), (R_0) and (sigma ) on N, which we have studied in the previous sections. The experiments involved (Nsimeq 10^4) so that the droplets were in the compressible regime; due to large three-body atom losses the whole dynamics was required to involve small time intervals of around 15 ms.
Illustrative examples of frontal collisions of two droplets leading to its coalescence. Identical hetero-nuclear droplets of Rb and K atoms are considered in each column, where the density of Rb atoms integrated along the z-direction ({{tilde{rho }}}^{(Rb)}(x,y) = int dz rho ^{(Rb)}(x,y,z)) is depicted at different evolution times. The initial translational kinetic energy is varied to guarantee that for (N^{Rb} = 5 times 10^5) (first column), (N^{Rb} = 1 times 10^6) (second column) and (N^{Rb} = 5 times 10^6) (third column) the Weber number equals (mathrm {We}_2^{(A1)}= 6.2 ).
Illustrative examples of frontal collisions of two droplets leading to an initial coalescence followed by a disintegration of the resulting droplet into two droplets that are then expelled in opposite directions. Initially, identical hetero-nuclear droplets of Rb and K atoms are considered in each column, where the density of Rb atoms integrated along the z-direction ( {tilde{rho }}^{(Rb)} (x, y) = int dz rho ^{(Rb)} (x, y, z)) is depicted at different evolution times. The initial translational kinetic energy is varied to guarantee that for (N^{(Rb)} = 5 times 10^5) (first column), (N^{(Rb)} = 1 times 10^6) (second column) and (N^{(Rb)} = 5 times 10^6) (third column), the Weber number equals (hbox {We}^{(A1)}_2 =14.5).
Illustrative examples of frontal collisions of two droplets leading to an initial coalescence followed by a disintegration of the resulting droplet into two droplets that are then expelled in opposite directions and a third droplet is formed at the center.. Initially, identical hetero-nuclear droplets of Rb and K atoms are considered in each column, where the density of Rb atoms integrated along the z-direction ({{tilde{rho }}}^{(Rb)}(x,y) = int dz rho ^{(Rb)}(x,y,z)) is depicted at different evolution times. The initial translational kinetic energy is varied to guarantee that for (N^{Rb} = 5 times 10^5) (first column), (N^{Rb} = 1 times 10^6) (second column) and (N^{Rb} = 5 times 10^6) the Weber number equals (mathrm {We}_2^{(A1)}= 31.9).
At first, we have studied the effects of binary collisions of quantum droplets in an EGPE scheme described by Eqs. (38) without including three body atom losses. We performed simulations for homo-nuclear ((^{39})K) and hetero-nuclear ((^{41})K–(^{87})Rb) systems, with the parameters reported in Section SMA of the Supplementary Material. Different initial kinetic energies and number of particles were studied always in the incompressible regime where no self-evaporation is expected. The initial condition corresponds to identical ground state droplets with an overall movement described by Eq. (32). Later on, simulations including the experimental three-body scattering coefficients (K_{alpha beta gamma } ne 0) were performed.
In general, for the kinetic energies here considered, the simulations yield different possible final results of a frontal collision,
Coalescence: Formation of a single droplet from two colliding droplets. The resulting droplet is found to be created mainly in a quadrupole excitation, although higher order modes could also be present (see Fig. 3 ).
Disintegration into two droplets after coalescence: after the collision, the oscillations of the resultant coalesced droplet are so strong that it breaks into two similar droplets. These two droplets exit in opposite directions and leave the collision zone with a translational movement described mainly as a dipole like oscillation (See Fig. 4).
Disintegration into three droplets after coalescence: when the initial droplets are large and the collision is energetic enough, after the breaking and the expulsion of two droplets in opposite directions, a third droplet is formed at the center. The latter oscillates mainly in a quadrupole mode as in the first case (illustrated in Fig. 5). Disintegration into more droplets is expected for even larger translational kinetic energies.
Similar results for homo-nuclear droplets are described in Refs^11,18. In the simulations no quasi-elastic collisions were observed, though they have been predicted to occur in one dimensional simulations³⁴.
The Weber number given by Eq. (35) must be evaluated according to the expression of the surface tension adequate for the regime at which the initial droplets are prepared: in the incompressible regime ansatz 1 is used while in the compressible regime ansatz 2 is the correct one. The droplets dynamics seen in the collisions tend to be mainly quadrupolar as can be infered by comparing Fig. (1) to Figs. (3–5), so that the quadrupolar Weber number is proposed as a figure of merit. In the illustrated simulations, (N^{Rb}) is compatible with excitations described by ansatz 1, and the Weber number describes the qualitative outcomes. This assertion is also supported by the numerical calculation of the multipole moment of the quantum liquid after a frontal collision. Figure 6, illustrates the Weber number (We_{ell =2}) as a function of the number of atoms and the wave vector (k_0). When three-body losses are included, the Weber number is still a good parameter to delimit the immediate breakdown of a droplet after the collision. However the transient regime is linked to a different band of (mathrm {We}_2) values. A rich variety of transient configurations that culminate in one of the three above results are found. Their observation depends highly on the inclusion of three-body atom losses. In the ideal case (K_{alpha beta gamma } = 0), frequently after the break into two droplets like that described in (2), the dipole like oscillation of each individual droplet is faster than their relative translation and afterwards coalescence occurs-described by (1), the final single droplet remains oscillating mainly in a quadrupole way. The transient regime usually takes place for quadrupole Weber numbers between two given values (mathrm {We}_2^{min}) and (mathrm {We}_2^{max}) whose values are exemplified in Fig. 6 for collisions between the droplets. For (mathrm {We}_2 < mathrm {We}_2^{min}) the scenario described by (1) takes place without a transient stage. For (mathrm {We}_2 > mathrm {We}_2^{max}) the scenario described by either (2) or (3) takes place without a transient stage. Figure 6 is supported by extensive numerical simulations carried out for the values of the parameters (k_0xi ) and N reported in it.
When three-body losses are included, the observation of the results of the collision is strongly dependent on the initial conditions. For a given (vec{k}_0), the parameter (vec{d}_0) should be chosen so that atom losses by three-body scattering are small-less than (2%) before the collision takes place- and, simultaneously, the droplets overlap at (t = 0) is avoided. The mean life time of the hetero-nuclear droplets is in the range of tens of milliseconds as expected from calculations in “Atom losses” section, and confirmed by experimental evidence⁶. So that, the long time evolutions like those described in the ideal case are difficult to observe for low (k_0) values. Nevertheless, in spite of the short times involved, the regimes at which the scenarios (1)–(3) occur were observed in the simulations. An interesting effect due to the atomic gas that surrounds each one of the droplets is also predicted: the surrounding evaporated gas created nearby a given droplet interact with the atoms from the other droplet through contact terms. This results in a reconfiguration of a merged quantum droplet. The effect is similar to that exemplified by Fig. 1c,d. Taking into account three-body losses, the Weber number is still a good parameter to delimit the immediate breakdown of a coalesced droplet after the collision.
Quadrupole Weber number (mathrm {We}_2) as a function of the number of atoms of a given species is illustrated for (a) homunuclear and (b) hetero-nuclear droplets (c) hetero-nucleat droplets including atom losses. In the first case N corresponds to half the total number of atoms in the droplet and in case (b), (c) to (N^{(Rb)}) in an Rb and K admixture. The outcome of a frontal collision will give rise to the coalescence of the droplets for (mathrm {We}_2) below the lower line (C region). The disintegration of the originally droplet formed during the collision is observed for (mathrm {We}_2) above the higher line (D region). For values of (mathrm {We}_2) between these lines involve a long transient stage where the coalesced droplet oscillates strongly before its disintegration (OCD region). The dots in (b) correspond to the values of (k_0xi ) and (N^{(Rb)}) illustrated in previous figures.
In this work we have studied theoretically the static and dynamical properties of quantum droplets constituted by binary mixtures of homo-nuclear and hetero-nuclear ultracold atoms. One of the more interesting features of quantum droplets corresponds to self-trapping. Our calculations show that the effective equations as those built as extensions to the Gross–Pitaevskii formalism provide robust interpretations of the droplets behavior both in the compressible and the incompressible regime. Our calculations confirm the slight differences between the expectations derived from Monte Carlo studies and the LHY-EGPE. They manifest not only on the minimum number of atoms required for self-trapping, but also on the radius (R_0) and thickness of the surface dR of the droplets; their dependence on the scattering lengths is discussed at depth in Ref.^16,17. We have also shown that the analytical model Eq. (28) is a reasonable representation of the density of atoms (rho _0^{(alpha )}).
We found that in the incompressible and compressible regimes the elementary Bogolubov excitations, Eqs. (15–17), are better described in a variational basis by Eqs. (18) and (19) respectively. We provide an analytic expresision for these excitations for any given multipole moment (ell ). The spherical symmetry of the droplets is directly incorporated in the selection of the basis. The corresponding expressions for the surface tension reflect such a symmetry, and exhibit a different behavior on the multipole order. For ansatz 1, the (ell )-dependent coefficient that accompanies the expression of the energy of the excitations coincides with that expected for classical liquids and that of the liquid drop model of nuclei in the limit of (Nrightarrow infty ) nucleons. Explicit calculations of the excitation energies show that this ansatz is the more adequate one to describe the incompressible regime of the quantum droplets. The excitation modes and the surface tension derived from ansatz 2 are more adequate for the quantum fluid in the compressible regime. Differences in the predicted shape of the droplets between the effective range formulation and the EGPE manifest also on slight differences in the excitation energies of the Bogolubov modes.
Self-evaporation is a particularly interesting phenomenon predicted for quantum droplets that is not easily observed due to its competition with atom losses arising from three-body scattering. We have identified two other relevant mechanisms that accelerate the disintegration of the quantum droplet: loss of the adequate proportion of each atomic species and differences in spatial configuration of each species.
In spite of the difficulties to experimentally study the dynamics of quantum droplets, we have theoretically shown that if the droplets are excited via frontal binary collisions it is still possible to observe a non trivial dynamics. The Weber number resulting from the surface tension expressions introduced in our study exhibit a simple dependence on the radii of the droplets (R_0) and the relative droplets momentum (hbar k_0) which are measured in terms of the surface width dR, Eqs. (36, 37). These make explicit the surface character of the emerging excitations. We have also shown that the Weber number evaluated with the adequate ansatz define different regimes of excitation that include conditions for (1) the formation of pairs of quantum droplets excited as oscillating dipoles that detach from each other leaving the impact region, (2) the formation of a single droplet (mainly in a quadrupole excitation mode) that merges from two elementary droplets; (3) if three-body losses are included, space regions are created where the lost atoms interact via de contact terms and lead to a reconfiguration of otherwise disjoint droplets (this effect is similar to that exemplified by Fig. SM 3 c-d of the Supplementary Material but it is enhanced by the evaporated atoms). Contrary to analysis that evaluate the surface tension of the planar interface described by local energy functionals and later on incorporate curvature through Tolman like terms³³, no extra parameters are required within our formalism to identify the different collision regimes.
The Bogolubov functions introduced in this work and selected via variational criteria, give a reliable representation of the low energy excitations of quantum droplets constituted by binary mixtures with a finite number of homo- and hetero-nuclear atoms. To our knowledge, the dynamical definition of the surface tension of a quantum droplet that we introduce via the Rayleigh interpretation (25) has not been reported yet in the literature. Our ansatze give the energy estimation of a surface excitation up to a multipole order (ell ).
As the results of this work do not depend on the strength of the interactions, albeit the ground state of the droplets, these could be extrapolated to many different droplets with similar geometries. The basis of our work could be applied to different systems such as dipolar droplets, droplets with impurities, or the novel mixed bubble phase²³.
We have focused our study to low energy surface excitations that can be induced by collisions between ground state droplets. Other interesting excitation modes have been predicted, in particular, vortical quantum droplets that could provide a three-dimensional realization of 3D topological solitons in nonlinear media³⁵. In this case, the bulk excitation yield self-trapped vortex tori. It would result very interesting to complement the careful numerical study performed in³⁵ by a variational study of these high energy excitations.
Our variational expressions are useful to describe the qualitative and quantitative expectations of frontal binary collisions of the droplets as a function of the most relevant parameters, such as number of atoms, interaction strengths, translation kinetic energy. This panorama is valid even in the presence of atom losses where the Weber number- here defined and calculated—is a good figure of merit to discriminate the expected outcomes of such collisions. Two interesting effects concerning atom losses are predicted, the first corresponds to the saturation of the exponent that describes the time dependence of atom losses via three body scattering (details in SMC). The second effect arises in atom collisions when the surrounding evaporated gas created nearby a given droplet interact with the atoms from the other droplet through contact terms, resulting in a reconfiguration of a merged quantum droplet.
We hope that these theoretical analyses encourages the research of quantum droplets in both the compressible and the incompressible regimes, the latter fairly unexplored due to the experimental difficulties related to its generation and characterization. These theoretical and experimental studies would pave the way to further studies of novel excited states induced in ultracold atoms, and their analog to other quantum matter systems.
The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.
Bulgac, A. Dilute quantum droplets. Phys. Rev. Lett. 89, 050402 (2002).
Article  ADS  PubMed  Google Scholar 
Petrov, D. S. Quantum mechanical stabilization of a collapsing Bose–Bose mixture. Phys. Rev. Lett. 115, 155302 (2015).
Article  ADS  CAS  PubMed  Google Scholar 
Lee, T. D., Huang, K. & Yang, C. N. Eigenvalues and eigenfunctions of a Bose system of hard spheres and its low-temperature properties. Phys. Rev. 106, 1135 (1957).
Article  ADS  MathSciNet  CAS  MATH  Google Scholar 
Cabrera, C. et al. Quantum liquid droplets in a mixture of Bose–Einstein condensates. Science 359, 301 (2018).
Article  ADS  MathSciNet  CAS  PubMed  Google Scholar 
Semeghini, G. et al. Self-bound quantum droplets of atomic mixtures in free space. Phys. Rev. Lett. 120, 235301 (2018).
Article  ADS  CAS  PubMed  Google Scholar 
D’Errico, C. et al. Observation of quantum droplets in a hetero-nuclear bosonic mixture. Phys. Rev. Res. 1, 033155 (2019).
Article  Google Scholar 
Ferrier-Barbut, I., Kadau, H., Schmitt, M., Wenzel, M. & Pfau, T. Observation of quantum droplets in a strongly dipolar Bose gas Phys. Rev. Lett. 116, 215301 (2016).
Article  ADS  Google Scholar 
Schmitt, M., Wenzel, M., Böttcher, F., Ferrier-Barbut, I. & Pfau, T. Self-bound droplets of a dilute magnetic quantum liquid. Nature (London) 539, 259 (2016).
Article  ADS  CAS  Google Scholar 
Chomaz, L. et al. Quantum-fluctuation-driven crossover from a dilute Bose–Einstein condensate to a macrodroplet in a dipolar quantum fluid. Phys. Rev. X 6, 041039 (2016).
Google Scholar 
Tanzi, L. et al. Observation of a dipolar quantum gas with metastable supersolid properties. Phys. Rev. Lett. 122, 130405 (2019).
Article  ADS  CAS  PubMed  Google Scholar 
Ferioli, G. et al. Collisions of self-bound quantum droplets. Phys. Rev. Lett. 122, 090401 (2019).
Article  ADS  CAS  PubMed  Google Scholar 
Bertsch, G. F. Capillary waves in a quantum liquid. Phys. Rev. A 9, 819 (1974).
Article  ADS  Google Scholar 
Bertsch, G. F. Dynamics of heavy ion collisions. In Nuclear Physics of Heavy Ions (eds. Balian, R., Rho, M., & Ripka, G.) (North Holland Publishing Company, Amsterdam, 1978).
Casas, M. & Stringari, S. Elementary excitations of (^4)He clusters. J. Low Temp. Phys. 79, 135 (1990).
Article  ADS  CAS  Google Scholar 
Böttcher, F. et al. Dilute dipolar quantum droplets beyond the extended Gross–Pitaevskii equation. Phys. Rev. Res. 1, 033088 (2019).
Article  Google Scholar 
Cikojević, V., Dzelalija, K., Stipanović, P., Vranjec Markić, L. & Boronat, J. Ultradilute quantum liquid drops. Phys. Rev. B 97, 140502(R) (2018).
Article  ADS  Google Scholar 
Cikojević, V., Poli, E., Ancilotto, F., Vranjes̃-Markić, L. & Boronat, J. Dilute quantum liquid in a K-Rb Bose mixture. Phys. Rev. A 104, 033319 (2021).
Article  ADS  Google Scholar 
Cikojevic, V. et al. Dynamics of equilibration and collisions in ultradilute quantum droplets. Phys. Rev. Res. 3, 043139 (2021).
Article  CAS  Google Scholar 
Frohn, A. & Roth, N. Dynamics of Droplets (Springer, 2000).
Book  MATH  Google Scholar 
Thouless, D. J. Vibrational states of nuclei in the random phase approximation. Nuclear Phys. 22, 78 (1961).
Article  ADS  MathSciNet  CAS  MATH  Google Scholar 
Bohm, D. & Pines, D. A collective description of electron interactions: III. Coulomb interactions in a degenerate electron gas. Phys. Rev. 92, 609 (1953).
Article  ADS  MathSciNet  CAS  MATH  Google Scholar 
Wouters, S., Nakatani, N., Van Neck, D. & Kin-Lic Chan, G. The Thouless theorem for matrix product states and subsequent post-density matrix renormalization group methods. Phys. Rev. B 88, 075122 (2013).
Article  ADS  Google Scholar 
Naidon, P. & Petrov, D. S. Mixed bubbles in Bose–Bose mixtures. Phys. Rev. Lett. 126, 115301 (2021).
Article  ADS  MathSciNet  CAS  PubMed  Google Scholar 
Rayleigh, L. On the capillary phenomena of jets. Proc. R. Soc. Lond. 29, 71 (1879).
Article  Google Scholar 
Bohr, A. & Mottelson, B. R. Nuclear Structure Vol. II (W. A. Benjamin, New York, 1975).
MATH  Google Scholar 
Fort, C. & Modugno, M. Self-evaporation dynamics of quantum droplets in (^{41})K – (^{87})Rb mixture. Appl. Sci. 11, 866 (2021).
Article  CAS  Google Scholar 
Wenzel, M., Böttcher, F., Langen, T., Ferrier-Barbut, I. & Pfau, T. Striped states in a many-body system of tilted dipoles. Phys. Rev. A 96, 053630 (2017).
Article  ADS  Google Scholar 
Böttcher, F. et al. Transient supersolid properties in an array of dipolar quantum droplets. Phys. Rev. X 9, 011051 (2019).
Google Scholar 
Chomaz, L. et al. Long-lived and transient supersolid behaviors in dipolar quantum Gases. Phys. Rev. X 9, 021012 (2019).
CAS  Google Scholar 
Ferioli, G. et al. Dynamical formation of quantum droplets in a (^{39})K mixture. Phys. Rev. Res. 2, 013269 (2020).
Article  CAS  Google Scholar 
Barranco, M. et al. Helium nanodroplets: An overview. J. Low Temp. Phys. 142, 1 (2006).
Article  ADS  CAS  Google Scholar 
Dupont, N. et al. Quantum state control of a Bose–Einstein condensate in an optical lattice phys. Rev. X Quantum 2, 040303 (2021).
ADS  Google Scholar 
Tolman, R. C. The effect of droplet size on surface tension. J. Chem. Phys. 17, 333 (1949).
Article  ADS  CAS  Google Scholar 
Astrakharchik, G. E. & Malomed, B. A. Dynamics of one-dimensional quantum droplets. Phys. Rev. A 98, 013631 (2018).
Article  ADS  CAS  Google Scholar 
Kartashov, Y. V., Malomed, B. A., Tarruell, L. & Torner, L. Three-dimensional droplets of swirling superfluids. Phys. Rev. A 98, 013612 (2018).
Article  ADS  CAS  Google Scholar 
Download references
This work was partially supported by the grants DGAPA-PAPIIT, UNAM: IN-103020, IN-109619, UNAM-AG810720, and CONACYT Ciencia Básica: A1-S-30934 and Laboratorios Nacionales 315838. E. Alba-Arroyo thanks CONACyT for his posgraduate studies fellowship 464666 and DGAPA-PAPIIT IN-103020 graduation support.
Departamento de Física Cuántica y Fotónica, Instituto de Física, Universidad Nacional Autónoma de México, 04510, Mexico City, Mexico
J. E. Alba-Arroyo, S. F. Caballero-Benitez & R. Jáuregui
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
J.E.A.-A. performed the numerical calculations, contributed to the theoretical description and to the writing of the manuscript. S.F.C.-B. contributed to the development of numerical algorithms, to the discussion on the theoretical description and to the writing of the manuscript. R.J. contributed to the theoretical description and to the writing of the manuscript.
Correspondence to R. Jáuregui.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Alba-Arroyo, J.E., Caballero-Benitez, S.F. & Jáuregui, R. Weber number and the outcome of binary collisions between quantum droplets. Sci Rep 12, 18467 (2022). https://doi.org/10.1038/s41598-022-22904-8
Download citation
Received: 06 August 2022
Accepted: 20 October 2022
Published: 02 November 2022
DOI: https://doi.org/10.1038/s41598-022-22904-8
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Advertisement
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
© 2022 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

source

In Arizona Senate Race, Libertarian Drops Out and Backs Blake Masters  The New York Times
source

If you don’t remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Please enter a term before submitting your search.
Access provided by

DOI: https://doi.org/10.1016/j.lanepe.2022.100539

Permitted

For non-commercial purposes:

• Read, print & download
• Redistribute or republish the final article
• Text & data mine
• Translate the article (private use only, not for distribution)
• Reuse portions or extracts from the article in other works

Not Permitted

• Sell or re-use for commercial purposes
• Distribute translations or adaptations of the article

Elsevier’s open access license policy
COVID-19 Resource Centre
Access the latest 2019 novel coronavirus disease (COVID-19) content from across The Lancet journals as it is published.
We use cookies to help provide and enhance our service and tailor content and ads. To update your cookie settings, please visit the Cookie Preference Center for this site.
Copyright © 2022 Elsevier Inc. except certain content provided by third parties. The content on this site is intended for healthcare professionals.

Privacy Policy   Terms and Conditions   Accessibility

source

Play Now
Football Pick’em
Play Now
College Football Pick’em
Week 9 Recap
Talking the business of sports
Tennessee @ Kansas City
Current Records: Tennessee 5-2; Kansas City 5-2
The Tennessee Titans are 4-1 against the Kansas City Chiefs since December of 2016, and they’ll have a chance to extend that success on Sunday. The Titans are on the road again on Sunday and play against Kansas City at 8:20 p.m. ET Nov. 6 at GEHA Field at Arrowhead Stadium. These two teams come into the matchup bolstered by wins in their previous games.
Tennessee was able to grind out a solid victory over the Houston Texans this past Sunday, winning 17-10. RB Derrick Henry had a stellar game for Tennessee as he rushed for two TDs and 219 yards on 32 carries.
Meanwhile, Kansas City’s game against the San Francisco 49ers two weeks ago was close at halftime, but Kansas City turned on the heat in the second half with 30 points. The Chiefs put the hurt on San Francisco with a sharp 44-23 win. Kansas City’s win was all the more impressive since San Francisco was averaging only 14.83 points allowed on the season. Kansas City relied on the efforts of WR Mecole Hardman, who punched in two rushing touchdowns in addition to snatching one receiving TD, and QB Patrick Mahomes, who passed for three TDs and 423 yards on 34 attempts.
Tennessee is the clear underdog, so they’re hoping it’s the size of the fight that determines this one. They are currently five-for-five against the spread in their most recent games, however, a trend bettors might want to take into account.
The wins brought both teams up to an identical 5-2. A couple offensive numbers to keep in mind before kickoff: The Titans are worst in the league in yards per game, with only 285.6 on average. The Chiefs’ offense has more to brag about, as they they enter the contest with 27 overall offensive touchdowns, which is the best in the NFL. So the cards are definitely stacked in their favor.
The Chiefs are a big 12.5-point favorite against the Titans, according to the latest NFL odds.
Over/Under: -110
See NFL picks for every single game, including this one, from SportsLine’s advanced computer model. Get picks now.
Tennessee have won four out of their last five games against Kansas City.
© 2004-2022 CBS Interactive. All Rights Reserved.
CBS Sports is a registered trademark of CBS Broadcasting Inc. Commissioner.com is a registered trademark of CBS Interactive Inc.
Images by Getty Images and US Presswire

source

Eunisses Hernandez, an incoming council member, is grappling with what it means to gain power in an unequal system
In 2014, Eunisses Hernandez, a young Latina activist, watched as a Los Angeles politician was sworn in as the first Latino leader of California’s state senate in more than a century.
“I got chills the whole time,” Hernandez recalled on a recent Monday afternoon at a restaurant in Highland Park. The background of Kevin de León, the child of working-class immigrants who had now won serious political power in the state, reminded her of her own family, and of herself. She was confident that De León, a longtime activist for immigrants’ rights, was “going to help us pass better laws”.
Today, Hernandez, is about to take her own seat on Los Angeles’ city council, and is leading the calls for De León’s resignation, after leaked audio captured racist conversation between De León and other powerful Latino Democrats in the city.
The scandal has prompted a reckoning over anti-Black racism and colorism within Latino communities, as well as an ongoing political crisis in America’s second-largest city, where even Joe Biden has called for all of the local Democratic officials involved in the racist conversation to resign. Nury Martinez, the former president of the Los Angeles city council, has stepped down, but De León and Gil Cedillo, the incumbent Democrat Hernandez is set to replace after she defeated him in the June primary, have not.
For Los Angeles’ next generation of Latino politicians, including Hernandez, the crisis is forcing an even deeper reckoning with what it means to gain political power in a deeply unequal system.
“It’s a different place to be in when the people who have these values and views, anti-Blackness and racism and homophobia, are the ones in leadership that look like us,” Hernandez said.
De León described himself as a progressive and was a prominent voice in California’s environmental justice movement. Cedillo, the incumbent Hernandez defeated this June, had been endorsed by Bernie Sanders, and spent years fighting to allow the state to issue driver’s licenses to undocumented workers and provide financial aid to undocumented students.
The racism advanced by politicians like Donald Trump is “in your face”, Hernandez argued. “This type of racism and anti-Blackness and homophobia is covered by faces that have been adored and loved by community, and that people have looked up to.”
The leaked audio was recorded during a meeting between the three influential Latino city council members and Ron Herrera, a prominent Latino labor leader, last year. The substance of the conversation was how to redraw the political districts in the city to consolidate their own personal power, as well as what they perceived as Latino political power in the city.
The tape captured the lawmakers making disparaging remarks about Black people, members of the gay community and mocked the appearance of Indigenous people from Oaxaca, a community that has hundreds of thousands of people in Los Angeles. “I’m glad they’re wearing shoes,” Cedillo was recorded as saying. Cedillo and De León did not immediately respond to requests for comment.
Emerging just weeks before she is set to take up her seat in the city council, the leaked recording prompted “a rollercoaster of emotions” for Hernandez, from anger to sadness to disappointment to anger again. “I felt that they took us back several decades,” she said of the lawmakers involved.
Before the audio leaked, she said, she had just heard Herrera, the labor leader, speak publicly about how to build Black and Latino solidarity in Los Angeles.
Hernandez, the 32-year-old daughter of Mexican immigrants, grew up in Highland Park, a predominantly Latino neighborhood in east Los Angeles that in recent years has seen intense gentrification, soaring housing prices and the displacement of longtime Latino residents.
Hernandez challenged Cedillo, the district’s longtime Democratic incumbent this spring, with a campaign focused on the issue of housing affordability. She accused Cedillo of supporting the interests of real estate developers over his constituents and pledged to put the needs of renters first.
She was targeted with attack ads over her support of police abolition, but won decisively in the primary anyway, making her Los Angeles’ first abolitionist city council member and a key example of a potential leftward swing in Los Angeles politics.
For someone who spent the past decade as a progressive organizer trying to “build bridges” between Black and Latino residents in order to change drug laws, close local jails and secure public funding for alternatives to incarceration, Hernandez said, the bigoted remarks from people seen as Latino champions is particularly frustrating.
Today, Latino residents make up nearly 50% of Los Angeles’ total population: non-Hispanic Black residents make up 9%. Interracial tensions have played out for decades as a result the lack of investment in both Black and Latino communities in Los Angeles, Hernandez said.
There has long been anti-Blackness in Latino communities, Hernandez said, and it’s something she’s encountered “in a lot of spaces that I come from”.
“We’ve always been put against each other. The bureaucracy runs on scarcity mode, ‘If we give it to them, we can’t give it to you,’” Hernandez said.
Much of the work for organizers on the ground “is about breaking that narrative”, she said, and showing how Black and Latino residents can combine power to achieve a shared goal, as Hernandez and Black state lawmaker Isaac Bryan did on their campaign for Measure J, in which voters approved spending 10% of the county’s unrestricted funds on social services, mental health programs, youth intervention programs, and other alternatives to incarceration.
The leaked audio had other lessons for young politicians, including open disparagement of the city council’s two most progressive members, and a suggestion that the district represented by Nithya Raman, a progressive on the council who had just unseated a Democratic incumbent in 2020, was “the one to put in the blender and chop up left and right” when district boundaries were redrawn. (California’s attorney general announced an investigation in Los Angeles’ redistricting in response to the audio leak.)
Raman’s district was dramatically changed by the redistricting process, forcing her to get to know an entirely new swath of constituents, an outcome Hernandez said she saw as “intentional” and designed, “to let other people know: ‘Better watch yourself, because this is what’s going to happen to you.’”
Hernandez said she believes one outcome to the scandal might be election victories for more leftwing city council candidates, and a chance to grow the council’s progressive bloc from two members to five. She hopes progressive candidates can have more of an impact in the city’s approach to public safety and housing.
Hernandez said she believed the revelations from the audio recording had political relevance far outside Los Angeles, particularly in other nominally liberal cities struggling with homelessness and massive spending on police and law enforcement.
Her home town is “one of the most beautifully diverse cities of any place in the world”, and one of the wealthiest, and it’s in a state with a Democratic supermajority. Yet Los Angeles is also suffering “one of the worst housing crises”, Hernandez said, and there are an estimated 41,000 people living on the streets. It was time for voters to be more cautious about who they trusted. “Look where we are,” she said.

source

Play Now
Football Pick’em
Play Now
College Football Pick’em
Week 9 Recap
Talking the business of sports
The Houston Astros are one win away from the second World Series championship in franchise history. They held on to beat the Philadelphia Phillies in Game 5 on Thursday, and now hold a 3-2 series lead. Historically, teams with a 3-2 lead in a best-of-seven have gone on to win the series 70 percent of the time. The Astros are in driver’s seat.
Now, that said, the Astros were in this exact same position three years ago. They took a 3-2 series lead home to Minute Maid Park in the 2019 World Series, but the Washington Nationals won Games 6 and 7 to steal away the title. The Nationals showed the world what the Phillies are trying to do can be done. It won’t be easy, but it is doable.
You can view the full postseason schedule by clicking here, and here’s the 2022 playoff bracket. Now, let’s get to the important stuff for Game 6.
Date: Saturday, Nov. 5 | Time: 8:03 p.m. ET
Location: Minute Maid Park (Houston, Texas)
TV channel: FOX | Live stream: fuboTV (try for free)  
Odds: HOU -145; PHI +122 | 7 o/u (via Caesars Sportsbook)
Starting pitchers: LHP Framber Valdez (17.6, 2.82 ERA) vs. RHP Zack Wheeler (12-7, 2.82 ERA)
Game 6 is rematch of the Game 2 pitching matchup and the Astros opened that game with three consecutive doubles (on only four pitches!) against Wheeler. They led 3-0 six batters into the game and, needless to say, the Phillies need Wheeler to be a lot better than that to have a chance at a comeback. Valdez was dominant in Game 2, holding the Phillies to one run in 6 1/3 innings, and only five of the 24 batters he faced hit the ball out of the infield in the air. Kyle Schwarber and Bryce Harper has been great all series. Nick Castellanos, Rhys Hoskins, and J.T. Realmuto have not, and the Phillies will need something from their big righty bats against Valdez to keep their season alive Saturday.
Wheeler allowed five runs (four earned) in five innings in Game 2 and I have a hard time believing a pitcher this good will be that bad in back-to-back starts. I’m going to say the Phillies escape with a narrow win in Game 6 and force a winner-take-all Game 7 on Sunday. I admit this is a hope more than a prediction. I am not ready for the baseball season to end and Game 7s rule, so I hope we get one this series.
Pick: Phillies 3, Astros 2
© 2004-2022 CBS Interactive. All Rights Reserved.
CBS Sports is a registered trademark of CBS Broadcasting Inc. Commissioner.com is a registered trademark of CBS Interactive Inc.
Images by Getty Images and US Presswire

source