Fermions meet two bosons -- the heteronuclear Efimov effect revisited

In this article, we revisit the heteronuclear Efimov effect in a Bose-Fermi mixture with large mass difference in the Born-Oppenheimer picture. As a specific example, we consider the combination of bosonic $^{133}\mathrm{Cs}$ and fermionic $^6\mathrm{Li}$. In a system consisting of two heavy bosons and one light fermion, the fermion-mediated potential between the two heavy bosons gives rise to an infinite series of three-body bound states. The intraspecies scattering length determines the three-body parameter and the scaling factor between consecutive Efimov states. In a second scenario, we replace the single fermion by an entire Fermi Sea at zero temperature. The emerging interaction potential for the two bosons exhibits long-range oscillations leading to a weakening of the binding and a breakup of the infinite series of Efimov states. In this scenario, the binding energies follow a modified Efimov scaling law incorporating the Fermi momentum. The scaling factor between deeply bound states is governed by the intraspecies interaction, analogous to the Efimov states in vacuum.


Introduction
The late Mahir S. Hussein, to whom this special issue is dedicated, is best known for his ground-breaking contributions in nuclear physics, e.g., to the understanding of the fusion and break-up of loosely bound nuclei [1]. From a very early stage on, Mahir also recognized the importance of ultracold atomic gases, in particular Bose-Einstein condensates, for reaching a deeper understanding of strongly-correlated many-body quantum systems. In particular, based on his solid background in nuclear physics, he was among the first to realize the immense potential of magnetically tuned Feshbach resonances for tuning effective interactions over a wide range [2]. Therefore, it appeared appropriate to us to discuss one specific, rather spectacular example of the application of tunable interparticle interactions on the quantum level, namely the so-called Efimov effect. This effect, originally proposed in the context of nuclear binding [3], nicely demonstrates the fruitful crossfertilization of concepts between nuclear physics and ultracold atomic quantum gases, which was so dear to Mahir S. Hussein.
The unexpected emergence of an infinite series of bound three-particles states via resonant pairwise interaction was first predicted theoretically by Vitaly Efimov in 1970 [3]. Since then, the Efimov effect has become a prime example for studying universality in quantum few-body systems. In the universal regime, i.e. when the scattering length a between two particles exceeds the characteristic range of interparticle interaction r 0 , there exists a series of three-body bound states which follow a simple discrete scaling law. Remarkably, these three-body bound states can form, even though the twobody interactions are too weak to support a two-body bound state. Initially, Efimov's prediction raised serious doubts about the applicability of his concepts, but theorists trying to prove him wrong had to finally concede that he might actually be right. While Efimov proposed the effect to be observed in nuclear systems, such as in 3 He + or in the Hoyle state of 12 C, the existence of a series of three-bound states requires the two-body interactions to be close to resonance. It was only in 2006, 36 years after Efimov's theoretical prediction, that first experimental evidence was eventually found in a gas of ultracold 133 Cs atoms [4]. The precise tuning of twobody interactions from infinite repulsion to attraction via Feshbach resonances opened up unique opportunities to study the range of universality in the Efimov scenario. Following the investigations in homonuclear systems, the Efimov effect was later also observed in heteronuclear systems [5][6][7][8][9][10][11] where, e.g., a large mass ratio [8,9] leads to a denser Efimov spectrum which allowed for the observation of up to three consecutive Efimov resonances.
Besides its relevance in the field of few-body physics, the Efimov effect plays an important role in the understanding of many-body systems in the presence of threebody bound states. Similarly, one can consider the Efimov effect from a few-body perspective and study the influence of a surrounding many-body background to the trimer. The latter has been studied in different scenarios with either a Fermi Sea [12][13][14][15] or a BEC [16,17] serving as a background. The former has been studied in [18][19][20], following recent experimental advances in creating Bose [21][22][23][24][25] and Fermi polarons [26][27][28][29][30][31][32]. More details on different aspects of the Efimov effect are covered extensively in several reviews, which have appeared in the years after its first realization [33][34][35][36][37][38][39].
In this work, we will focus on the specific case of two heavy bosons and a light fermion, which allows to apply the Born-Oppenheimer (BO) approximation. We investigate two limiting cases of the Efimov scenario, first, in vacuum, and second, in the presence of a Fermi Sea. While the first case reproduces the wellknown features of the Efimov effect, the second case provides novel insights serving as a precursor to understand effective interactions of Fermi polarons [40], i.e., strongly correlated impurities in a Fermi sea. As a specific example, we consider the heteronuclear mixture of bosonic 133 Cs and fermionic 6 Li, for which experiments have been performed by us and others. Providing a simple and intuitive access to understand the Efimov effect in a mass-imbalanced system, the BO approximation does not only capture the existence of an infinite series of three-body bound states [33,[41][42][43], but can also account for intraspecies interactions via short-range van der Waals (vdW) potentials [44]. Our paper is structured as follows: after giving a basic introduction to for a < 0 and a > 0, respectively. The intraspecies scattering length a BB between the particles B is resonant. Short-range two-body interactions (red shaded area) limit the discrete scaling law between Efimov trimers. the Efimov scenario in Sec. 2, we describe how we solve the Schrdinger equation by employing the BO approximation in Sec. 3. In Sec. 4, we present our results on the Efimov energy spectrum for a system consisting of two Cs atoms and one Li atom, taking into account finite Cs-Cs s-wave interaction. We extend our results to a system of two Cs atoms immersed in a Fermi Sea (Sec. 5), before concluding in Sec. 6.

The Efimov Scenario
The Efimov scenario for the heteronuclear case can be visualized in an energy diagram ( Fig. 1) which is plotted against the inverse interspecies scattering length 1/a between the identical bosons B and the distinguishable particle X. For simplicity, we first consider the scenario in which the two bosons resonantly interact (a BB → ∞). Then the energy diagram can be divided into three different regions representing the threebody scattering states, Efimov trimers, and atom-dimer states. For positive energies (E > 0), the three atoms are unbound and possess a finite kinetic energy. For negative energies (E < 0), one needs to distinguish between the sign of the scattering length. On the positive scattering length side (a > 0), the system supports a weakly bound dimer state BX with an energy of 1/2µa 2 , where µ is the reduced mass. Above the atom-dimer threshold, the dimer state BX coexists with a free atom B. The Efimov trimers exist in the region below the three-body dissociation threshold at a < 0 and the atom-dimer threshold at a > 0. In this region, an infinite number of Efimov trimers with energies E n exist which cross the three-body dissociation threshold at values of a (n) − and a (n) * on the negative and positive scattering length side, respectively. These crossings can be described by the discrete scaling laws where the scaling factor λ = e π/s0 with the dimensionless parameter s 0 is dependent on the mass ratio, the number of resonant interactions and the quantum statistics of the particles [39].
In real systems, interactions have a finite range, given by non-zero interparticle distances in the underlying interaction potentials. Therefore, in the infinite progression of Efimov trimers a ground state has to be considered. This can be done by introducing a threebody parameter (3BP) which defines the position of the energy E 0 or the scattering length a (0) − of the lowest Efimov trimer. The 3BP is determined by shortrange two-body interactions and can be expressed in terms of vdW units (see Eq. 4). As indicated in Fig.  1, the universal Efimov scaling therefore only holds if |a| max(r BX vdW , r BB vdW ) due to finite range effects.

Born-Oppenheimer approximation
In the Born-Oppenheimer approximation, we can solve the three-body Schrdinger equation for two heavy bosons with mass M and a light atom with mass m in a twostep approach. We assume that the mass ratio is M/m 1, such that the light atom immediately follows the motion of the heavy ones. In the first step, the Schrdinger equation is solved for the light atom with a potential created by the two heavy scatterers at fixed distances. The resulting energy serves as an interaction potential V E (R), induced by the presence of the light particle, in the Schrdinger equation for the heavy bosons where R denotes the distance between the two heavy bosons. The boson-boson interaction is modeled by a van der Waals potential with a hard core of the form [45,46] V BB (R) = ∞, The C 6 coefficient allows us to naturally introduce the van der Waals radius r vdW and energy E vdW via [47] and The cut-off radius R 0 in Eq. 3 determines the 3BP and can be directly related to the boson-boson scattering length a BB via [45] where J ν (x) and N ν (x) are Bessel functions of the first and second kind, respectively. As Eq. 6 has more than one solution, the value of R 0 does not only determine a BB , but also the number of bound states supported by the vdW potential V BB (R). However, we note that the exact number of dimer states is irrelevant for our purposes as it has no significant effect on the long-range part of the Efimov wavefunctions [44]. The induced interaction potential V E (R) possesses a symmetric and an antisymmetric solution [43]. In the limit of a → ∞, only the former is relevant for the discussion of bound states. It reads where c ≈ 0.567 is the solution of c = e −c and is connected to the scaling factor s 0 via With the interaction potentials V E (R) and V BB (R) we can solve Eq. 2 in order to understand the Efimov scenario for a Cs-Cs-Li system and its scaling behavior taking into account boson-boson interactions.

Two bosons meet one fermion
In Fig. 2 we show the solution of Eq. 2 using the Li-Cs mass ratio of M Cs /m Li = 22.1 and a short-range cutoff R 0 for which the energy spectrum supports up to two Cs 2 dimer states. It is instructive to consider the two cases For case a), i.e. pure two-body vdW interaction, we identify two weakly bound dimer states (orange dashed lines). The energy of the least bound state approaches the binding energy of the Cs 2 dimer E b = 1/M a 2 CsCs for positive and increasing a CsCs . The energy of the most deeply bound state crosses the Cs-Cs resonance and shows a step-like behavior around a CsCs = r vdW . This behavior marks the crossover between a vdWdominated (a CsCs < r vdW ) dimer and a halo state (a CsCs > r vdW ). For case b) (blue lines), we find that the energy of the most deeply bound state closely follows case a), such that we can assign this state to the Cs 2 dimer state. The next bound state E n=0 does not follow the respective Cs 2 dimer state anymore, but persists also across the Cs-Cs resonance. This clearly shows the effect of the mediated interaction necessary to form Efimov trimer states. The following bound states with energies E n≥1 correspond to the infinite progression of Efimov states. They show a gradual step around a CsCs = r vdW which we assign to an overlap of the wavefunctions of the Efimov trimers with the Cs 2 halo dimer. This crossover between states originating from a short-range molecular vdW potential and a long-range ∝ −1/R 2 potential is one of the main results of the finite range BO approximation, providing an intuitive access to the three-body problem. From the calculated energy spectrum, we can additionally extract the scal- ing factor between adjacent energy levels. For λ 2 n := E n /E n−1 the scaling factor between the two most deeply bound states amounts to λ 2 n=1 = 42.9 at resonance. When crossing the resonance towards negative a CsCs , the scaling factor approaches the universal value of λ 2 n→∞ = (5.63) 2 = 31.7 for a pure ∝ −1/R 2 potential. The deviation close to resonance can again be explained by the existence of the weakly bound Cs 2 dimer state. We find that already the second deepest scaling factor λ 2 n=2 ≈ 31.9 is close to the universal value.

Two bosons meet the Fermi Sea
Let us now consider two heavy bosons immersed in a Fermi Sea. In the BO approximation the effective Limediated potential between the two bosons can be calculated by [48]: Here ∆E(R) denotes the energy reduction of the interacting system compared to the free system with the bound state wavevectors and scattering phase shifts described by κ ± and δ ± , respectively [48]. In the asymptotic limit, the energy reduction ∆E(R → ∞) is equivalent to the chemical potential of two free, heavy atoms in a Fermi Sea [49]. The length scale of the Fermi Sea is governed by the Fermi wavevector k F which is related to the atomic density n via k F = (6π 2 n) 1/3 . In Fig. 3 (left panel) we plot the total potential V (R) = V eff (R) + V BB (R) for which we solve the Schrdinger equation for unitarity a LiCs = a CsCs = ∞. The grey shaded area marks the hard wall for R < R 0 . If we set k F r vdW = 0 (blue line), we recover the potential from the Efimov scenario consisting of a short-range ∝ −1/R 6 vdW potential (grey dashed line) and a longrange Efimov potential ∝ −1/R 2 (grey dotted line). For increasing k F , the effective potential starts to grow a repulsive barrier around R ≈ k −1 F showing damped oscillations. In the inset of Fig. 3, this behavior can be seen more clearly for even larger values of k F . The form of the potential is reminiscent of the form of Friedel oscillations [50] which arise due to the sharp edge of the Fermi distribution. In fact, in the limit of a LiCs R, k −1 F , V eff (R) takes the same form of the Friedel oscillations or of the RKKY interaction [51] in the context of magnetic interactions.
We now use the effective potential V eff (R) to replace the Efimov potential V E (R) in Eq. 2 and solve for the eigenenergies of the system (right panel of Fig.  3). We note that the chosen values of k F r vdW = 0.01 and k F r vdW = 0.05 correspond to densities between 10 11 cm −3 and 10 13 cm −3 . For k F r vdW = 0.01, the Cs 2 dimer state and the deepest two Efimov states remain unchanged. For the following state, we see that the effective potential starts to deviate from the Efimov potential around R ≈ 70 r vdW leading to a smaller binding energy. As the potential gets more repulsive for increasing R, the formation of bound states is completely suppressed. Similarly, for k F r vdW = 0.05, the weakening of the binding can be understood from a deviation of the Efimov potential around R ≈ 10 r vdW and the suppression of the infinite series of Efimov states begins one state earlier. Analogous to the case of two Cs atoms and one Li atom, we want to investigate the role of the intraspecies interaction, as shown in Fig. 4. The energies in the case k F r vdW = 0 again represent the previous Efimov bound state energies. Introducing a Fermi Sea to the system (k F r vdW = 0.05) we see that, on the positive scattering length side, the n = 1 state shows first deviations from the Efimov scaling when a CsCs takes values larger than r vdW at the step going from the vdW-dominated to the long-range regime. Crossing the resonance towards negative scattering lengths, this deviation gets larger again at the step −a CsCs ≈ r vdW . Following the energy line further to the positve scattering length side, the state rapidly dissociates before it reaches the long-range regime. In the same way, with a smaller density of the Fermi Sea of k F r vdW = 0.01, the system supports one more bound state before the binding energy vanishes.
We calculate the density dependence of the binding energies (Fig. 5) and find that they remain nearly constant before they rapidly go to the continuum. The number of bound states is given by the Fermi wavevector k F and the shape of the lines suggest a similar scaling as we have seen in Fig. 1 following Eq. 1. Indeed, in the presence of the Fermi Sea, the two heavy bosons follow a new discrete scaling law including the additional length scale k F [12,14]: In the case of finite intraspecies interactions we find that this scaling is fulfilled for high-lying bound states n > 1 (inset of Fig. 5). Analogous to the scenario of two Cs atoms and one Li atom in Sec. 4, the first excited bound state (n = 1, orange line) shows only a small deviation from the scaling, whereas the scaling is broken for the ground state due to finite range effects.

Conclusions
In summary, we have calculated the binding energies of a Cs-Cs-Li system and of a system of two Cs atoms in a Li Fermi Sea and studied the influence of the intraspecies scattering length using the BO approximation. In the Cs-Cs-Li system, the intraspecies interaction leads to a step-like behavior in the energy spectrum and the existence of weakly bound Cs 2 dimers influence the scaling factor. Immersing the two Cs atoms in a Li Fermi Sea suppresses the formation of bound states for sufficiently high k F and breaks the discrete Efimov scaling law. Instead, a new scaling law can be formulated which takes into account the wavevector k F . This additional length scale of the Fermi Sea may also be used to define a new window of universality which is not only determined by short-range interactions, but also by the Fermi wavevector. In an experiment, the shifted position of the Efimov states in the Fermi Sea, which can also be interpreted as bipolaronic states, may be observed by means of three-body loss measurements, similar to the previous Efimov experiments [52] while now lower temperatures and higher densities are required. The influence of the intraspecies scattering length can be studied in the Li-Cs system which features two interspecies Feshbach resonances around 843 G and 889 G with negative and positive sign of the intraspecies scattering length, respectively. However, we note that our simple BO approximation does only provide qualitative results. For the Cs-Cs-Li system, more quantitative results beyond the BO approximation can be obtained by means of a spinless vdW theory [53] where the threebody problem is solved in the hyperspherical formalism with two-body interactions modeled by a single channel Lennard Jones potential. Also finite temperature effects as well as scattering of trimers by the Fermi Sea which may lead to the excitation of particle-hole pairs and a change of the effective interaction potential [13] have to be considered for a more realistic description of the system.