Asymmetric lineshapes of Efimov resonances in mass-imbalanced ultracold gases

The resonant profile of the rate coefficient for three-body recombination into a shallow dimer is investigated for mass-imbalanced systems. In the low-energy limit, three atoms collide with zero-range interactions, in a regime where the scattering lengths of the heavy-heavy and the heavy-light subsystems are positive and negative, respectively. For this physical system, the adiabatic hyperspherical representation is combined with a fully semi-classical method and we show that the shallow dimer recombination spectra display an asymmetric lineshape that originates from the coexistence of Efimov resonances with St\"uckelberg interference minima. These asymmetric lineshapes are quantified utilizing the Fano profile formula. In particular, a closed form expression is derived that describes the width of the corresponding Efimov resonances and the Fano lineshape asymmetry parameter $q$. The profile of Efimov resonances exhibits a $q-$reversal effect as the inter- and intra-species scattering lengths vary. In the case of a diverging asymmetry parameter, i.e. $|q|\to \infty$, we show that the Efimov resonances possess zero width and are fully decoupled from the three-body and atom-dimer continua, and the corresponding Efimov metastable states behave as bound levels.


I. INTRODUCTION
The Efimov effect is one of the most counter-intuitive phenomena in few-body physics where an infinity of three-body bound states is formed even when the scattering length of the two-body subsystems is negative [1][2][3][4]. This phenomenon was theoretically predicted by V. Efimov to occur for 3 equal mass particles that interact via zero-range potentials, with trimer binding energies that scale geometrically [5]. The existence of these exotic trimer states was experimentally confirmed by Kraemer et al. in an ultracold gas of Cs atoms [6]. This suggested new possibilities for theoretical and experimental investigations [1,2,[7][8][9][10][11][12][13] to address various physical aspects of the Efimov states, such as the discrete scale invariance of the trimer binding energies [14] or the sensitivity of the ground state energy on the short-range physics. In particular, the latter stems from the fact that within the zero-range model the trimer spectrum is unbound from below due to Thomas collapse [15] and an auxiliary parameter, i.e. three-body parameter, was introduced in order to specify the ground state energy rendering the entire spectrum system dependent [5]. However, experimental and theoretical advances demonstrated that for ultracold atoms the Efimov spectrum exhibits a certain class of universality, i.e. van der Waals universality [16][17][18][19][20][21][22][23][24][25][26]. Namely, it was shown that the lowest Efimov state appears at scattering lengths a (1) − ≈ −10ℓ vdW with ℓ vdW being the length scale of van der Waals interactions between two neutral atoms.
Mass-imbalanced ultracold gases are an ideal platform to explore more deeply the idiosyncrasies of Efimov * pgiannak@pks.mpg.de spectra. In particular, three-body collisions of ultracold atoms with unequal masses offer more favorable experimental conditions that enable observation of multiple successive trimer states and measurement of their geometrical energy scaling, i.e. the smoking gun of Efimov physics [27][28][29]. Apart from that, mass-imbalanced ensembles offer a large parameter space, such as the particles' mass-ratio, the sign and magnitude of the interand intra-species scattering lengths, which provide fertile ground to investigate pristine attributes of the Efimov states. Specifically, theoretical and experimental efforts have mapped out a large portion of the parameter space addressing the underlying physics of recombination processes in heavy-heavy-light (HHL) systems [27,28,[30][31][32][33]. The particular case which stands out corresponds to HHL systems that possess inter-and intraspecies scattering lengths of opposite sign, i.e. a HL < 0 and a HH > 0 respectively. For example, in the experimental works of Refs. [30,34], it was demonstrated in the regime of broad Fano-Feshbach resonances [35] that the lowest Efimov state is in good agreement with the predictions of the universal zero-range and van der Waals theory. However, subsequent experimental investigations show that deviations from the universal theory are more pronounced for narrow Fano-Feshbach resonances [36]. Furthermore, within the zero-range theory Ref. [37] illustrated that the diabaticity of the three-body collisions imposes additional limitations on the universal properties of Efimov spectrum, where mostly adiabatic collisions yield trimer states independent of the three-body parameter, as was pointed out in the case of Refs. [30,34].
Additionally, Ref. [37] showed that three-body recombination into a shallow heavy-heavy dimer possesses a unique property that only mass-imbalanced systems exhibit, namely the co-existence of Efimov resonances with Stückelberg suppression effects in the same range of scat-tering lengths. In this work, we further study this particular attribute of HHL systems and demonstrate that the corresponding Efimov resonances in the recombination rate coefficient plotted versus scattering length can display an asymmetric profile, which can be quantified by the Fano profile formula. In particular, our analysis employs the adiabatic hyperspherical framework for zerorange two-body interactions which is combined with a fully semi-classical theory [37]. Also, a simplified version of the semi-classical approach is shown where the lowest hyperspherical curves are approximated by universal potential tails at large hyperradii, as in Ref. [38]. This permits us to derive closed form relations for the S-matrix elements which are expressed in terms of the width of the Efimov resonance and Fano's lineshape asymmetry parameter q. As an example, the asymmetric profiles of the Efimov resonances in the recombination coefficient of 6 Li − 133 Cs − 133 Cs and 6 Li − 87 Rb − 87 Rb are analyzed, both of which showcase a q-reversal phenomenon as function of the inter-and intraspecies scattering length ratio. Furthermore, we observe that for a diverging q parameter the Efimov resonances behave as bound states that are embedded in the continuum [39]. This occurs since the decay width of the resonances vanishes as q → ∞ and the corresponding Efimovian quasi-bound states decouple from the three-body and atom-dimer continua.
The structure of this work is as follows: In Section II the Hamiltonian of the three-body system and the parameters of interest are given. Sections II A and II B provide a detailed review of the methods that are employed in our analysis. More specifically, Section II A discusses the adiabatic hyperspherical representation and the fully semi-classical treatment of the coupled hyperradial equations. In Section II B a simplified version of the semi-classical theory is given that permits us to express the S− matrix elements of recombination processes into shallow dimers in terms of the inter-and intra-species scattering lengths. Finally, Section III focuses on the asymmetric profile of Efimov resonances in the spectrum of the three-body recombination coefficient for HHL systems.

II. GENERAL CONSIDERATIONS AND METHODS
Consider a three-body system that consists of two heavy (H) alkali atoms and a light (L) one at lowenergies. The particles mutually interact through s-wave pairwise interactions that are modeled via Fermi-Huang's zero-range pseudopotential. Our greatest interest here is in the regime where the mass-imbalanced system can recombine into a shallow heavy-heavy dimer with a recoiling light atom. This scenario arises for inter-and intraspecies interactions of opposite sign meaning that the scattering length between a heavy-light or heavy-heavy pair of particles is a HL < 0 or a HH > 0, respectively. Furthermore, relaxation and recombination processes into deep dimer channels will be neglected, and we will focus on the physics that arises due to energies near the breakup threshold, i.e. the zero-energy limit. For this purpose, we focus on the two lowest potential curves of HHL systems, which suffice to describe three-body recombination processes into shallow dimers as was shown in Ref. [37] permitting the derivation of closed form expressions for the S−matrix.
A. The adiabatic hyperspherical representation and the semi-classical approach The total three-body Hamiltonian for the HHL system of interest is given by the following expression: where V ij represents the Fermi-Huang pseudopotential. a ij and µ ij refer to the scattering length and two-body reduced mass of the ij-pair of particles, respectively. ∇ 2 i denotes the Laplacian for r i , and m i indicates the mass of the i-th particle. Note that the scattering lengths a ij between the atoms are chosen to be larger than any other length scale of the system permitting us to focus on the universal characteristics of the three-body system under consideration.
Utilizing the Jacobi vector choice of Ref. [40], Eq. (1) can be separated into the Hamiltonians of center of mass and relative degrees of freedom. Since the s-wave interactions involve only the relative distance between a pair of particles, the center-of-mass Hamiltonian is fully decoupled, meaning that the relative one retains all the relevant information of the three-body system. Therefore, we focus only on the relative Hamiltonian which gives, after transforming it into hyperspherical coordinates (for details see [2]), the following expression: where indicates the three-body reduced mass, R is the hyperradius, and Ω is a collective coordinate denoting the five hyperangles [41,42]. H ad (R; Ω) represents the part of the Hamiltonian which contains the hyperangular centrifugal potential as well as the two-body interactions expressed in the hyperspherical coordinates.
whereΛ denotes the grand angular momentum operator. ν (R a HH ) with a HH > 0 and a HL < 0. The red (blue) line saturates at large hyperradii in the atom+dimer (three-body break-up) threshold. The quantities Φ U L , and Φ U L indicate the JWKB phase accumulation in the upper potential curve. For the lower potential the corresponding phase is denoted by Φ L L . The vertical dashed line represents the hyperradius where the non-adiabatic coupling P −matrix element P12 maximizes. The horizontal dotted line refers to the three-body collisional energyĒ in units of

HH
, and the three-body parameter, r 3b a HH depicted by the blue region.
In the spirit of the adiabatic hyperspherical representation, the properly symmetrized three-body wave function is provided by the following ansatz: where F ν (R) and φ ν (R; Ω) indicate the ν-th hyperradial and hyperangular part of the wave function, respectively. In particular, φ ν (R; Ω) components of Ψ(R, Ω) are obtained by diagonalizing Eq. (3) at a fixed hyperradius R.
where the eigenvalues U ν (R) are the so-called adiabatic hyperspherical potential curves. Substitution of Eqs. (4) and (5) into the Schrödinger equation of the Hamiltonian H rel and integration over all the hyperangles Ω yields a set of coupled ordinary second-order differential equations that solely depend on the hyperradius R.
where V νν ′ (R) indicate the non-adiabatic coupling matrix elements/operators that are given by the following expressions: where ⟨. . .⟩ Ω denotes that the integration over the hyperangles only.
Owing to the zero-range interactions, the nonadiabatic coupling matrix elements, P νν ′ (R) and Q νν ′ (R), as well as the hyperspherical potential curves U ν (R) can be calculated semi-analytically, [1,40,43,44]. However, the resulting hyperspherical potential curves U ν (R), especially the lowest one, possess attractive singularities at the origin, i.e. the Thomas collapse. Therefore, an auxiliary parameter is introduced in order to truncate the attractive singularity in the potential curves, which in its simplest form consists of a hard wall placed at a small hyperradius, R ≈ r 3b . The three-body parameter r 3b is arbitrary (from the point of view of zerorange theory) and is usually fixed via experimental observations. In addition, the zero-range approximation greatly simplifies the computational cost since only hyperradial equations in Eq. (6) require numerical solution using standardized R-matrix methods [45][46][47]. Fig. 1 depicts the two lowest hyperspherical potential curves U 1 3 ν (R a HH ) as obtained from zero-range approximation. The upper (blue) potential that vanishes at large hyperradii R in the break-up threshold and the lower (red) potential which in the limit of large R approaches the energy of the HH dimer. The light blue region denotes the hard wall boundary condition at r 3b a HH that removes the attractive singularity of the lower curve. The potential curves of Fig. 1 suffice in order to intuitively understand the recombination of three free particles into a universal pair of atoms with a recoiling one. Consider the three-body system at a collisional energyĒ (in units of

HH
) indicated by the dotted line in Fig. 1. In particular, we are interested in the lowenergy limit in order to validate the two-channel approximation and highlight the threshold behavior of threebody collisions in HHL settings. Viewing this three-body system heuristically as a time-dependent collision, starting from infinite long distances, the three-particles propagate inwards in the upper potential curve and tunnel with some probability under the repulsive barrier and then probe the corresponding classical allowed region at short hyperradii. In this region, the non-adiabatic P −matrix element P 12 between the upper and lower potential curve plays a key role in inducing transitions. More specifically, at distances R LZ a HH (vertical dashed line) the corresponding P -matrix maximizes indicating the strong coupling regime. This means that the particles transition with a certain probability from the upper to the lower curve and subsequently propagate outwards, fragmenting into a two-body molecule plus a spectator atom. This recombination process is quantified mainly by evaluating the S 12 2 element of the scattering S−matrix. As was shown in Ref. [37], the S 12 2 matrix element can be obtained analytically within the two-channel approximation by combining the Landau-Zener physics with the Jeffreys-Wentzel-Kramers-Brillouin (JWKB) approach.
The main constituents of this semi-classical approach are depicted in Fig. 1. More specifically, we assume that the P −matrix element of the potential curves in Fig. 1 possesses a Lorentzian lineshape in the vicinity or R ≈ R LZ , and we include the Langer correction in JWKB integrals [48]. Under these considerations, the S 12 2 matrix element for the hyperspherical potential curves in Fig. 1 reads: where e −2τ indicates the tunneling probability in a single collision with the repulsive barrier of the upper potential curve in Fig. 1. The JWKB phases in the upper curve are indicated by the terms Φ U L , and Φ U R . More specifically, Φ U L is the phase accumulation from the far left classical turning point up to R ≈ R LZ whereas Φ U R is the JWKB integral from R ≈ R LZ up to the inner classical turning point of the repulsive barrier. Similarly, in the lower potential curve Φ L L corresponds to the phase accumulation between the hard wall (blue shaded region) located at r 3b a HH and R ≈ R LZ . Furthermore, p corresponds to the Landau-Zener non-adiabatic probability to transition from the upper to the lower hyperspherical potential curve in a single pass through the avoided crossing region. The non-adiabatic probability p is evaluated from the P −matrix elements which as we mentioned above are approximated to have a Lorentzian lineshape versus the hyperradius and a maximum at R ≈ R LZ [49]. λ is the Stokes phase and it is a correction added to the components of the hyperradial wave function, i.e. F ν (R) with ν = 1, 2, as they propagate through the non-adiabatic transition region [50,51]. The Stokes correction phase depends on the non-adiabatic probability p and it obeys the relation where δ = − ln p 2 . Fig. 2 shows the Stokes phase versus p where in the diabatic (adiabatic) limit, i.e. p = 1 (0), the Stokes phase tends to λ = −π 4 (0). Eq. (8) captures the two main effects that occur in HHL systems. More specifically, the roots of the numerator of Eq. (8) indicate the Stückelberg suppression effects minimizing the probability of the HHL system to recombine into a shallow dimer. On the other hand, the roots of the denominator in Eq. (8) denote the Efimov resonance phenomenon that enhances the recombination into a shallow dimer. An additional insight obtained by Eq. (8) is that the Stückelberg suppression effects depend on the threebody parameter due to the phase Φ L L accumulation in the lower potential curve. In principle, also the Efimov resonances depend on r 3b , however, as Eq. (8) suggests in the limit of adiabatic collisions, i.e. p ≪ 1, only the phase accumulation in the upper potential curve survives which are independent of the three-body parameter meaning that such collisions possess a universal character.
The degree of diabaticity p is depicted in Fig. 3, in the zero-energy limit, as a function of a HH a HL for different mass ratios m H m L covering in this manner the regime from strong-to-weak mass-imbalanced atomic ensembles. Note that we consider values of the ratio a HH a HL that correspond to a HL and a HH both being larger than the van der Waals length scales of the HL and HH pairs of atoms, respectively, ensuring the validity of the zerorange theory. In particular, we observe in Fig. 3 that for large mass ratio, i.e. m H m L = 21, the corresponding three-body collision is more diabatic than in the case of weak mass-imbalance, i.e. m H m L = 6.3. This means that HHL systems with strong mass-imbalance, i.e. m H m L > 21, can easily transition from the threebody continuum to the shallow dimer+atom channel implying that the corresponding recombination process is strongly affected by the three-body parameter r 3b . This behavior of the non-adiabatic probability p on m H m L can be understood in terms of the ratio of the P -matrix elements and the energy difference of the hyperspherical potential curves, i.e. ∆, at R = R LZ . According to Ref. [49], the probability p is given by the relation p = e −π∆ [4vP12(R LZ )] where v refers to the semiclassical velocity of the particles at R = R LZ . Thus, for a HL → −∞, the ratio of the energy gap ∆ and P 12 (R LZ ) increases as m H m L decreases yielding in return a decreasing probability p and vice versa.
As an example, Fig. 4 illustrates the scaled S−matrix element S12 2 (ka HL ) 4 for 6 Li − 87 Rb − 87 Rb three-body system at low-energies E = ̵ h 2 k 2 2µ . Fig. 4(a) corresponds to the semi-classical model using Eq. (8) and Fig. 4(b) refers to the case where the hyperradial equations are solved numerically within the R-matrix approach. Both panels are in excellent agreement and the qualitative features, i.e. the enhancement and suppression of S12 2 (ka HL ) 4 , are similar to those shown in Ref. [37]. In particular, as discussed in Ref. [37] the enhancement of S12 2 (ka HL ) 4 is associated with an Efimov resonance. Namely, the upper potential curve in Fig. 1 can support a quasi-bound state behind the repulsive barrier at specific values of the ratios a HL a HH and r 3b a HH . Therefore, for colliding energies E that match the energy of the quasi-bound three-body state, the atoms can easily tunnel under the barrier, where they can probe the nonadiabatic transition region and eventually hop with some probability to the Rb 2 + Li channel. Therefore, the presence of quasi-bound state in the upper potential curve in Fig. 1 causes S12 2 (ka HL ) 4 to be more pronounced. On the other hand, the suppression of S12 2 (ka HL ) 4 is a manifestation of Stückelberg physics due to destructive interference of  the alternative pathways, which prevents the three particles from exiting to infinity along the Rb 2 + Li channel. However, in Ref. [37] the 6 Li − 133 Cs − 133 Cs system was investigated and the corresponding S12 2 (ka HL ) 4 possesses one main qualitative difference from 6 Li − 87 Rb − 87 Rb. Specifically, the 6 Li − 133 Cs − 133 Cs system exhibits narrower Efimov resonances [see Fig.2(b) and (c) in Ref. [37]] than those shown Fig. 4 for 6 Li − 87 Rb − 87 Rb. This difference mainly arises from the fact that the collisions in 6 Li − 133 Cs − 133 Cs are more diabatic than in 6 Li − 87 Rb − 87 Rb. As Fig. 3 suggests the non-adiabatic probability p for 6 Li − 133 Cs − 133 Cs is much closer to unit that for the case of 6 Li − 87 Rb − 87 Rb. The lower values of p for 6 Li − 87 Rb − 87 Rb indicate the weak coupling of the quasi-bound Efimov state to the atom-dimer continuum which in return is manifested as a broad resonance in S12 2 (ka HL ) 4 matrix element.

B. A simplified semi-classical model
In the following, we focus on the derivation of a simplified semi-classical model based on the prescription given in Ref. [38]. Our goal is to unveil the scaling behavior of S 12 matrix element with respect to the length scales and the degree of diabaticity p that govern HHL systems by incorporating only the necessary approximations. Therefore, for our purposes from this point and on we assume that the colliding energy of the three atoms tends to zero, i.e. E = ̵ h 2 k 2 2µ → 0. As in Fig. 1, Fig. 5 illustrates the two-lowest hyperspherical potential curves which are properly parameterized using only the limiting tails of the curves of Fig. 1 in a piecewise manner. Namely, for the upper potential curve (blue line) in Fig. 5 the universal tail is shown for hyperradii ranging from the non-adiabatic transition region, i.e. R LZ , up to R ∼ γ a HL . Also, for R > γ a HL we consider only the tail of the repulsive barrier of the potential curve shown in Fig. 1 which falls off as U 1 (R) ∼ ̵ h 2 2µR 2 (15 4) with the outer classical turning point being located at R ∼ 2 k. In addition, the effects of motion along the upper potential curve for R < R LZ is mapped to an arbitrary phase Φ. For the lower curve, at small hyperadii we employ the universal tail whereas for R > R LZ we assume that the potential curve is constant with energy equal to the heavy-heavy dimer. Note that the parameters s 0 and s * 0 correspond to the universal Efimov scaling coefficients for two and three resonant two-body interactions, respectively, and they are tabulated in Ref. [52] for several HHL systems.
Based on the piecewise potential curves of Fig. 5 and considering the low-energy limit, i.e. k → 0, the tunneling amplitude e −τ , and the semi-classical phases Φ U L , Φ U R and (a) (b) Φ L L are given by the following expressions: where the dimensionless parameters β and γ define the interval of hyperradius R such that the upper potential curve has the form U 1 (R) = − ̵ h 2 2µR 2 (s 2 0 + 1 4). In general, β and γ are considered as free parameters and they can be fixed by a fitting procedure to experimental or numerical data. Also, recall that the above JWKB integrals include the Langer corrections. After substitution of Eqs. (10a) to (10c) into Eq. (8), the S−matrix element S 12 reads where the terms (1 − (γka HL 2) 8 16) ≈ 1, and (1 − (γka HL 2) 4 4) ≈ 1 since we focus on the low-energy regime, i.e. k → 0. The phases ψ 1 and ψ 2 obey the expressions ψ 1 = Φ − s * 0 ln β − π 4 and ψ 2 = Φ + s 0 ln(γ β), respectively. Eq. (11) captures the main properties of the S−matrix element S 12 shown in Fig. 4. The numerator of Eq. (11) describes the positions of the Stückelberg interference minima which, as shown in Fig. 4, scale logarithmically with respect to the ratio r 3b a HH . Also, the spacing between successive minima is constant on a logarithmic scale, and related to the universal Efimov scaling coefficient s * 0 . On the other hand, the roots of the denominator of Eq. (11) trace out the maxima of S12 2 (ka HL ) 4 in Fig. 4, i.e. the Efimov resonances, where the position of the successive resonances is defined by the s 0 universal factor. We note Eq. (11) due to its simple structure can be used as a fitting formula for experimental measurements by treating the (ψ 1 , ψ 2 , γ) or (Φ, β, γ) as fitting parameters. Fig. 4 demonstrates that recombination resonant features are intertwined with Stückelberg interference minima.

III. ASYMMETRIC LINESHAPES IN THREE-BODY RECOMBINATION COEFFICIENT
This constitutes a unique feature of massimbalanced systems since for homonuclear three-body collisions the corresponding S-matrix element exhibits either Efimov resonances or Stückeleberg suppression effects for negative or positive scattering lengths, respectively. Therefore, this section focuses on the lineshape of the S 12 2 squared matrix element plotted as a function of the ratio r 3b a HH at fixed values of a HL a HH . In order to demonstrate the asymmetric lineshape of the Efimov resonances in HHL systems suffices to consider a range of r 3b a HH values in the neighborhood of a Stückelberg minimum assuming a total colliding energy E ≈ 0.
Under these considerations, utilizing the Fano profile formula Eq. (11) can be expressed in terms of the width of the resonance, Γ, and the Fano q−parameter which describes the asymmetry of the profile of the S 12 2 .
where x = 2(ln r 3b a HH −x r ) Γ with x r referring to the values of the ratio ln r 3b a HH that minimize the denominator of Eq. (11) at fixed a HL a HH . Note that Eq. (12) has same functional form as the conventional Fano formula, i.e. σ = σ 0 (ǫ + q) 2 (ǫ 2 + 1) [53], where the ratio ln r 3b a HH is the independent variable instead of the energy. For Eq. (12) the Fano lineshape asymmetry parameter q and the width Γ are given by the following expressions: Note that Γ is dimensionless here, in contrast to the usual Fano lineshape where Γ has units of energy (or frequency). The three-body recombination coefficient of HHL systems can be expressed in terms of the S 12 matrix  Table I. element yielding the relation where k = 2µE ̵ h 2 with E being the total colliding energy of the three-body system.
For a total colliding energy E ≈ 0, Fig. 6(a) and (b) depict the scaled recombination coefficient m H K3 ̵ ha 4 HL versus the ratio ln r 3b a HH in the vicinity of a Stückelberg minimum for two three-body systems, i.e. 6 Li − 133 Cs − 133 Cs and 6 Li − 87 Rb − 87 Rb, respectively. More specifically, the symbols in both panels correspond to the full semiclassical calculations whereas the solid lines are obtained by Eq. (11), i.e. the simplified semi-classical model, using γ, ψ 1 and ψ 2 as fitting parameters.
Note that Table I summarizes the values of these parameters for both HHL systems exhibiting universal characteristics since they are independent of scattering length ratio a HL a HH . Therefore, in order to extract the values of the γ, ψ 1 and ψ 2 parameters suffices to fit only the semiclassical calculations for a HL a HH = 47.4 and a HL a HH = 101.1 in panels (a) and (b), respectively. However, the phases ψ 1 and ψ 2 and the amplitude γ depend on the mass ratio of the HHL system since the corresponding hyperspherical potential curves are strongly influenced by variations of m H m L . Evidently, both panels showcase the asymmetric profile of the Efimov resonance as a distinctive feature of HHL systems where Eq. (11) is in excellent agreement with the corresponding semi-classical calculations. In particular, in Fig. 6  show the width of the Efimov resonance Γ and the asymmetry parameter q versus the scattering length ratio a HL a HH for the 6 Li − 133 Cs − 133 Cs ( 6 Li − 87 Rb − 87 Rb) system, respectively. Note that the total colliding energy is set to zero. Also, Γ and q are obtained via Eqs. (13a) and (13b), respectively, using the universal parameters shown in Table I. function of the ratio a HL a HH is known as the q-reversal effect, where the asymmetry parameter q changes sign at a HL a HH ∼ 70. For 6 Li− 87 Rb− 87 Rb shown in Fig. 6(b) a similar behavior is observed demonstrating that the occurrence of q-reversal is independent of the particles' mass ratio. The q-reversal phenomenon is a manifestation of quantum interference and in HHL systems it occurs when s * 0 x r + ψ 1 + λ = nπ 2 with n being an integer. Additionally, Fig. 7 demonstrates the validity of the Fano lineshape formula given in Eq. (12). More specifically, Fig. 7 illustrates a comparison of the scaled recombination coefficient between the fitting of Eq. (11) (red and black dots) and the Fano lineshape formula from Eq. (12) (red and black solid lines) at low collisional energies. In particular, the red (black) symbols and lines refer to the 6 Li − 133 Cs − 133 Cs ( 6 Li − 87 Rb − 87 Rb) for a scattering length ratio aCsLi aCsCs = 67.8 ( a RbLi a RbRb = 126.2). We observe that the Fano lineshape formula from Eq. (12) is in good agreement with the corresponding semi-classical calculations of Eq. (11).
The width Γ of the Efimov resonances and the lineshape asymmetry q are shown in Fig. 8 Table I. In panels (b) and (d), we observe the q−reversal effect where at a HL a HH = 70 and a HL a HH = 140 the lineshape asymmetry q diverges. This implies that for large q parameters the recombination coefficient approaches a symmetric lineshape that is centered at x r . Furthermore, we observe that at q → ∞ the corresponding widths of the Efimov resonances tend to zero, i.e. Γ → 0, as is illustrated in Fig. 8(a) for 6 Li − 133 Cs − 133 Cs and Fig. 8(c) for 6 Li − 87 Rb − 87 Rb. This means that in this range of parameters the Efimovian quasi-bound state stabilizes into a bound one which is fully decoupled from the three-body and the atom-dimer continua. This counterintuitive phenomenon is known as bound state in the continuum and such states have been observed in various fields of physics [39].

IV. SUMMARY
In summary, the properties of three-body recombination processes into shallow dimers for HHL systems are investigated. Focusing on the low-energy regime, we consider inter-and intraspecies interactions that possess negative and positive scattering lengths, respectively, thereby highlighting the threshold behavior of such HHL systems. For this three-body system, we have reviewed the theoretical methods used in Ref. [37] and in particular, the semi-classical approach providing additional details on the Stokes phase and the degree of diabaticity p. Furthermore, a simplified version of the semi-classical method is derived by approximating the hyperspherical curves with piecewise potential tails as in Ref. [38]. The simplified semi-classical model provides closed form expressions of the S−matrix elements which describe the process of three free-particles recombining into the shal-low dimer+atom channel. Namely, we show that Eq. (11) captures all the main attributes of the recombination spectra for HHL systems, such as the asymmetric lineshape in the three-body recombination coefficient, the logarithmic scaling of the Efimov resonances and Stückelberg interference minima. In particular, Fig. 6 demonstrates that Eq. (11) can be used as a fitting formula for the recombination spectra in HHL systems since the parameters ψ 1 , ψ 2 and γ are insensitive on the scattering length ratio a HL a HH . Focusing on the resonant profile of the recombination coefficient, Eq. (11) is parameterized in terms of the width of the resonance Γ and the lineshape asymmetry q. This parameterization enables us to identify two emergent phenomena that occur only in heteronuclear three-body collisions: (i) the q-reversal effect which describes the change in the asymmetry of the profile of the three-body recombination coefficient as a function of the scattering length ratio a HL a HH and (ii) the modification of an Efimov resonance into a bound state embedded in the three-body and atom-dimer continua for q → ∞.