Next Article in Journal
Special Issue “Interaction of Ionizing Photons with Atomic and Molecular Ions”
Next Article in Special Issue
One-Dimensional Disordered Bosonic Systems
Previous Article in Journal
Oleg Zatsarinny (1953–2021): Memories by His Colleagues
Previous Article in Special Issue
Emergence of Anyons on the Two-Sphere in Molecular Impurities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymmetric Lineshapes of Efimov Resonances in Mass-Imbalanced Ultracold Gases

by
Panagiotis Giannakeas
1,* and
Chris H. Greene
2,3
1
Max-Planck-Institut für Physik Komplexer Systeme, Nöthnitzer Str. 38, D-01187 Dresden, Germany
2
Department of Physics and Astronomy, Purdue University, West Lafayette, IN 47907, USA
3
Purdue Quantum Science and Engineering Institute, Purdue University, West Lafayette, IN 47907, USA
*
Author to whom correspondence should be addressed.
Atoms 2021, 9(4), 110; https://doi.org/10.3390/atoms9040110
Submission received: 28 September 2021 / Revised: 30 November 2021 / Accepted: 3 December 2021 / Published: 7 December 2021
(This article belongs to the Special Issue Physics of Impurities in Quantum Gases)

Abstract

:
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ückelberg 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 | , 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.

1. 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 three 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,3,7,8,9,10,11,12] to address various physical aspects of the Efimov states, such as the discrete-scale invariance of the trimer binding energies [13] 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 [14] 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 [15,16,17,18,19,20,21,22,23,24,25]. 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 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 [26,27,28]. Apart from this, mass-imbalanced ensembles offer a large parameter space, such as the particles’ mass ratio, the sign and magnitude of the inter- and intra-species scattering lengths, which provide fertile ground to investigate the 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 [26,27,29,30,31,32]. The particular case that stands out corresponds to HHL systems that possess inter- and intra-species scattering lengths of opposite sign, i.e., a H L < 0 and a H H > 0 , respectively. For example, in the experimental works of Refs. [29,33], it was demonstrated in the regime of broad Fano–Feshbach resonances [34] 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 [35]. Furthermore, within the zero-range theory, Ref. [36] illustrated that the diabaticity of the three-body collisions imposes additional limitations on the universal properties of the Efimov spectrum, where mostly adiabatic collisions yield trimer states independent of the three-body parameter, as was pointed out in the case of Refs. [29,33].
Additionally, Ref. [36] 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 scattering 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 zero-range two-body interactions, which is combined with a fully semi-classical theory [36]. Moreover, 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. [37]. 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 a function of the inter- and intra-species 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 [38]. 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 2, the Hamiltonian of the three-body system and the parameters of interest are given. Section 2.1 and Section 2.2 provide a detailed review of the methods that are employed in our analysis. More specifically, Section 2.1 discusses the adiabatic hyperspherical representation and the fully semi-classical treatment of the coupled hyperradial equations. In Section 2.2, 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 3 focuses on the asymmetric profile of Efimov resonances in the spectrum of the three-body recombination coefficient for HHL systems.

2. General Considerations and Methods

Consider a three-body system that consists of two heavy (H) alkali atoms and a light (L) one at low energies. 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 intra-species interactions of opposite sign, meaning that the scattering length between a heavy–light or heavy–heavy pair of particles is a H L < 0 or a H H > 0 , respectively. Furthermore, relaxation and recombination processes into deep dimer channels are neglected, and we focus on the physics that arises due to energies near the break-up 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. [36], permitting the derivation of closed-form expressions for the S-matrix.

2.1. 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:
H tot = i = 1 3 ħ 2 2 m i i 2 + i > j V i j ( r i j ) , with V i j ( r i j ) = 4 π ħ 2 a i j 2 μ i j δ ( r i j ) r i j [ r i j × ] ,
where V i j represents the Fermi–Huang pseudopotential. a i j and μ i j refer to the scattering length and two-body reduced mass of the i j -pair of particles, respectively. i 2 denotes the Laplacian for r i , and m i indicates the mass of the i-th particle. Note that the scattering lengths a i j 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. [39], Equation (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:
H rel = ħ 2 2 μ R 5 / 2 2 R 2 R 5 / 2 + H ad ( R ; Ω ) ,
where μ = m 1 m 2 m 3 / ( m 1 + m 2 + m 3 ) m H / 1 + 2 m H / m L indicates the three-body reduced mass, R is the hyperradius, and Ω is a collective coordinate denoting the five hyperangles [40,41]. H ad ( R ; Ω ) represents the part of the Hamiltonian that contains the hyperangular centrifugal potential as well as the two-body interactions expressed in the hyperspherical coordinates.
H ad ( R ; Ω ) = ħ 2 2 μ Λ ^ 2 + 15 ħ 2 8 μ R 2 + i > j V i j ( R ; Ω ) ,
where Λ ^ denotes the grand angular momentum operator.
In the spirit of the adiabatic hyperspherical representation, the properly symmetrized three-body wave function is provided by the following ansatz:
Ψ ( R , Ω ) = 1 R 5 / 2 ν ϕ ν ( R ; Ω ) F ν ( R ) ,
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 Equation (3) at a fixed hyperradius R.
H ad ( R ; Ω ) ϕ ν ( R ; Ω ) = U ν ( R ) ϕ ν ( R ; Ω ) ,
where the eigenvalues U ν ( R ) are the so-called adiabatic hyperspherical potential curves.
Substitution of Equations (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.
d 2 d R 2 + 2 μ ħ 2 ( U ν ( R ) E ) F ν ( R ) = ν V ν ν ( R ) F ν ( R ) ,
where V ν ν ( R ) indicate the non-adiabatic coupling matrix elements/operators that are given by the following expressions:
V ν ν ( R ) = 2 P ν ν ( R ) d d R + Q ν ν ( R ) with
P ν ν ( R ) = ϕ ν ( R ; Ω ) | R ϕ ν ( R ; Ω ) Ω
Q ν ν ( R ) = ϕ ν ( R ; Ω ) | 2 R 2 ϕ ν ( R ; Ω ) Ω ,
where Ω denotes the integration over the hyperangles only.
Owing to the zero-range interactions, the non-adiabatic coupling matrix elements, P ν ν ( R ) and Q ν ν ( R ) , as well as the hyperspherical potential curves, U ν ( R ) , can be calculated semi-analytically [1,39,42,43]. 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 3 b . The three-body parameter r 3 b is arbitrary (from the point of view of zero-range theory) and is usually fixed via experimental observations. In addition, the zero-range approximation greatly simplifies the computational cost since only hyperradial equations in Equation (6) require a numerical solution using standardized R-matrix methods [44,45,46].
Figure 1 depicts the two lowest hyperspherical potential curves U ν 1 / 3 ( R / a H H ) 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 3 b / a H H that removes the attractive singularity of the lower curve. The potential curves of Figure 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 E ¯ (in units of ħ 2 m H a H H 2 ) indicated by the dotted line in Figure 1. In particular, we are interested in the low-energy limit in order to validate the two-channel approximation and highlight the threshold behavior of three-body 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 L Z / a H H (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. [36], 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 Figure 1. More specifically, we assume that the P-matrix element of the potential curves in Figure 1 possesses a Lorentzian lineshape in the vicinity or R R L Z , and we include the Langer correction in JWKB integrals [47]. Under these considerations, the | S 12 | 2 matrix element for the hyperspherical potential curves in Figure 1 reads:
| S 12 | 2 = e 2 τ p ( 1 p ) cos 2 ( Φ L U Φ L L π 4 + λ ) { ( 1 e 4 τ 16 ) [ p cos 2 ( Φ L L + Φ R U π 4 ) + ( 1 p ) cos 2 ( Φ L U + Φ R U + λ ) ] ( 1 e 2 τ 4 ) 2 p ( 1 p ) cos 2 ( Φ L U Φ L L π 4 + λ ) + e 4 τ 16 } 1 ,
where e 2 τ indicates the tunneling probability in a single collision with the repulsive barrier of the upper potential curve in Figure 1.
The JWKB phases in the upper curve are indicated by the terms Φ L U and Φ R U . More specifically, Φ L U is the phase accumulation from the far-left classical turning point up to R R L Z , whereas Φ R U is the JWKB integral from R R L Z 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 3 b a H H and R R L Z . 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 L Z [48]. λ 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 [49,50]. The Stokes correction phase depends on the non-adiabatic probability p and it obeys the relation
λ = arg Γ i δ π δ π ln δ π + δ π + π 4 ,
where δ = ln p 2 . Figure 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 ) . Equation (10) captures the two main effects that occur in HHL systems. More specifically, the roots of the numerator of Equation (10) 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 Equation (10) denote the Efimov resonance phenomenon that enhances the recombination into a shallow dimer. An additional insight obtained by Equation (10) is that the Stückelberg suppression effects depend on the three-body parameter due to the phase Φ L L accumulation in the lower potential curve. In principle, also the Efimov resonances depend on r 3 b ; however, as Equation (10) suggests, in the limit of adiabatic collisions, i.e., p 1 , only the phase accumulation in the upper potential curve survives, which is independent of the three-body parameter, meaning that such collisions possess a universal character.
The degree of diabaticity p is depicted in Figure 3, in the zero-energy limit, as a function of a H H / | a H L | 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 H H / | a H L | that correspond to | a H L | and a H H , both being larger than the van der Waals length scales of the HL and HH pairs of atoms, respectively, ensuring the validity of the zero-range theory. In particular, we observe in Figure 3 that, for a 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 three-body continuum to the shallow dimer–atom channel, implying that the corresponding recombination process is strongly affected by the three-body parameter r 3 b . 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 L Z . According to Ref. [48], the probability p is given by the relation p = e π Δ / [ 4 v P 12 ( R L Z ) ] , where v refers to the semi-classical velocity of the particles at R = R L Z . Thus, for a H L , the ratio of the energy gap Δ and P 12 ( R L Z ) increases as m H / m L decreases, yielding, in return, a decreasing probability p and vice versa.
As an example, Figure 4 illustrates the scaled S-matrix element | S 12 | 2 ( k a H L ) 4 for the 6 Li - 87 Rb - 87 Rb three-body system at low energies E = ħ 2 k 2 2 μ . Figure 4a corresponds to the semi-classical model using Equation (10), and Figure 4b 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 | S 12 | 2 ( k a H L ) 4 , are similar to those shown in Ref. [36]. In particular, as discussed in Ref. [36], the enhancement of | S 12 | 2 ( k a H L ) 4 is associated with an Efimov resonance. Namely, the upper potential curve in Figure 1 can support a quasi-bound state behind the repulsive barrier at specific values of the ratios | a H L | a H H and r 3 b a H H . 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 non-adiabatic transition region and eventually hop with some probability to the Rb 2 + Li channel. Therefore, the presence of a quasi-bound state in the upper potential curve in Figure 1 causes | S 12 | 2 ( k a H L ) 4 to be more pronounced. On the other hand, the suppression of | S 12 | 2 ( k a H L ) 4 is a manifestation of Stückelberg physics due to the destructive interference of the alternative pathways, which prevents the three particles from exiting to infinity along the Rb 2 + Li channel. However, in Ref. [36], the 6 Li - 133 Cs - 133 Cs system was investigated and the corresponding | S 12 | 2 ( k a H L ) 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 Figure 2b,c in Ref. [36]) than those shown Figure 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 Figure 3 suggests, the non-adiabatic probability p for 6 Li - 133 Cs - 133 Cs is much closer to the unit than 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 the | S 12 | 2 ( k a H L ) 4 matrix element.

2.2. 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. [37]. Our goal is to unveil the scaling behavior of the 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 on, we assume that the colliding energy of the three atoms tends to zero, i.e., E = ħ 2 k 2 / 2 μ 0 .
As in Figure 1, Figure 5 illustrates the two lowest hyperspherical potential curves, which are properly parameterized using only the limiting tails of the curves of Figure 1 in a piecewise manner. Namely, for the upper potential curve (blue line) in Figure 5, the universal tail U 1 ( R ) = ħ 2 2 μ R 2 ( s 0 2 + 1 / 4 ) is shown for hyperradii ranging from the non-adiabatic transition region, i.e., R L Z , up to R γ | a H L | . Moreover, for R > γ | a H L | , we consider only the tail of the repulsive barrier of the potential curve shown in Figure 1, which falls off as U 1 ( R ) ħ 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 L Z are mapped to an arbitrary phase Φ . For the lower curve, at small hyperadii, we employ the universal tail U 2 ( R ) = ħ 2 2 μ R 2 [ ( s 0 * ) 2 + 1 / 4 ] , whereas, for R > R L Z , 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. [51] for several HHL systems.
Based on the piecewise potential curves of Figure 5 and considering the low-energy limit, i.e., k 0 , the tunneling amplitude e τ and the semi-classical phases Φ L U , Φ R U and Φ L L are given by the following expressions:
e τ γ | a H L | 2 / k 4 / R 2 ( γ k a H L / 2 ) 2 , Φ L U = Φ ,
Φ R U β | a H L | γ | a H L | s 0 2 / R 2 s 0 ln γ | a H L | β a H H
and Φ L L r 3 b β | a H L | ( s 0 * ) 2 / R 2 s 0 * ln β a H H r 3 b .
where the dimensionless parameters β and γ define the interval of hyperradius R such that the upper potential curve has the form U 1 ( R ) = ħ 2 2 μ R 2 ( s 0 2 + 1 / 4 ) . In general, β and γ are considered free parameters and they can be fixed by a fitting procedure to experimental or numerical data. Moreover, recall that the above JWKB integrals include the Langer corrections. After substitution of Equation (12)–(14) into Equation (10), the S-matrix element S 12 reads
| S 12 | 2 ( k a H L ) 4 = γ 4 16 p cos 2 s 0 * ln r 3 b a H H + ψ 1 + λ { p 1 p × sin 2 [ s 0 ln | a H L | a H H + ψ 2 ( s 0 * ln r 3 b a H H + ψ 1 ) ] + cos 2 ( s 0 ln | a H L | a H H + ψ 2 + λ ) p cos 2 ( s 0 * ln r 3 b a H H + ψ 1 + λ ) } 1 ,
where the terms ( 1 ( γ k a H L / 2 ) 8 / 16 ) 1 and ( 1 ( γ k a H L / 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.
Equation (15) captures the main properties of the S-matrix element S 12 shown in Figure 4. The numerator of Equation (15) describes the positions of the Stückelberg interference minima, which, as shown in Figure 4, scale logarithmically with respect to the ratio r 3 b / a H H . Moreover, 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 Equation (15) trace out the maxima of | S 12 | 2 ( k a H L ) 4 in Figure 4, i.e., the Efimov resonances, where the position of the successive resonances is defined by the s 0 universal factor. We note that Equation (15), due to its simple structure, can be used as a fitting formula for experimental measurements by treating the ( ψ 1 , ψ 2 , γ ) or ( Φ , β , γ ) as fitting parameters.

3. Asymmetric Lineshapes in Three-Body Recombination Coefficients

Figure 4 demonstrates that recombination resonant features are intertwined with Stückelberg interference minima. This constitutes a unique feature of mass-imbalanced 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 3 b a H H at fixed values of | a H L | a H H . In order to demonstrate the asymmetric lineshape of the Efimov resonances in HHL systems, it suffices to consider a range of r 3 b a H H values in the neighborhood of a Stückelberg minimum, assuming a total colliding energy E 0 .
Under these considerations, utilizing the Fano profile formula, Equation (15) 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 .
| S 12 | 2 ( k a H L ) 4 = A ( x + q ) 2 x 2 + 1 , with A = γ 4 ( 1 p ) 16 sin 2 ( s 0 * x r + ψ 1 + λ ) { cos [ 2 ( s 0 ln | a H L | a H H s 0 * x r + ψ 2 ψ 1 ) ] + ( 1 p ) cos 2 ( s 0 * x r + ψ 1 + λ ) } 1 ,
where x = 2 ( ln r 3 b a H H x r ) / Γ , with x r referring to the values of the ratio ln r 3 b a H H that minimize the denominator of Equation (15) at fixed | a H L | a H H . Note that Equation (16) has the same functional form as the conventional Fano formula, i.e., σ = σ 0 ( ϵ + q ) 2 / ( ϵ 2 + 1 ) [52], where the ratio ln r 3 b a H H is the independent variable instead of the energy. For Equation (16), the Fano lineshape asymmetry parameter q and the width Γ are given by the following expressions:
q = 2 s 0 * Γ cot ( s 0 * x r + ψ 1 + λ ) and
Γ 2 2 = 1 p ( s 0 * ) 2 p [ cos 2 ( s 0 ln | a H L | a H H + ψ 2 + λ ) p cos 2 ( s 0 * x r + ψ 1 + λ ) + p 1 p sin 2 ( s 0 ln | a H L | a H H s 0 * x r + ψ 2 ψ 1 ) ] × { cos [ 2 ( s 0 ln | a H L | a H H s 0 * x r + ψ 2 ψ 1 ) ] + ( 1 p ) cos [ 2 ( s 0 * x r + ψ 1 + λ ) ] } 1 .
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 element, yielding the relation
K 3 = 64 ħ π 2 μ k 4 | S 12 | 2 ,
where k = 2 μ E / ħ 2 , with E being the total colliding energy of the three-body system.
For a total colliding energy E 0 , Figure 6a,b depict the scaled recombination coefficient m H K 3 ħ a H L 4 versus the ratio ln r 3 b a H H 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 semi-classical calculations, whereas the solid lines are obtained by Equation (15), i.e., the simplified semi-classical model, using γ , ψ 1 and ψ 2 as fitting parameters.
Note that Table 1 summarizes the values of these parameters for both HHL systems exhibiting universal characteristics, since they are independent of scattering length ratio | a H L | a H H . Therefore, in order to extract the values of the γ , ψ 1 and ψ 2 parameters, it suffices to fit only the semi-classical calculations for | a H L | / a H H = 47.4 and | a H L | / a H H = 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 in m H / m L . Evidently, both panels showcase the asymmetric profile of the Efimov resonance as a distinctive feature of HHL systems, where Equation (15) is in excellent agreement with the corresponding semi-classical calculations. In particular, in Figure 6a, we observe that, for scattering length ratios in the range 55 < | a H L | a H H < 69 , the Efimov resonance occurs to the left of the Stückelberg minimum and its width decreases towards | a H L | a H H 70 . For | a H L | a H H > 70 , the Efimov resonance emerges to the right of the Stückelberg minimum with an increasing width. This behavior of the resonant structure as a function of the ratio | a H L | a H H is known as the q-reversal effect, where the asymmetry parameter q changes sign at | a H L | a H H 70 . For 6 Li - 87 Rb - 87 Rb shown in Figure 6b, 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, Figure 7 demonstrates the validity of the Fano lineshape formula given in Equation (16). More specifically, Figure 7 illustrates a comparison of the scaled recombination coefficient between the fitting of Equation (15) (red and black dots) and the Fano lineshape formula from Equation (16) (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 | a CsLi | a CsCs = 67.8 ( | a RbLi | a RbRb = 126.2 ). We observe that the Fano lineshape formula from Equation (16) is in good agreement with the corresponding semi-classical calculations of Equation (15).
The width Γ of the Efimov resonances and the lineshape asymmetry q are shown in Figure 8 for 6 Li - 133 Cs - 133 Cs (see panels (a) and (b)) and 6 Li - 87 Rb - 87 Rb (see panels (c) and (d)). Γ and q are obtained via Equations (17) and (18) using the universal parameters of Table 1. In panels (b) and (d), we observe the q-reversal effect, where, at | a H L | a H H = 70 and | a H L | a H H = 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 Figure 8a for 6 Li - 133 Cs - 133 Cs and Figure 8c 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 that is fully decoupled from the three-body and the atom–dimer continua. This counter-intuitive phenomenon is known as a bound state in the continuum and such states have been observed in various fields of physics [38].

4. 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 intra-species 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. [36] 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. [37]. The simplified semi-classical model provides closed-form expressions of the S-matrix elements that describe the process of three free particles recombining into the shallow dimer–atom channel. Namely, we show that Equation (15) 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, Figure 6 demonstrates that Equation (15) can be used as a fitting formula for the recombination spectra in HHL systems since the parameters ψ 1 , ψ 2 and γ are insensitive to the scattering length ratio | a H L | a H H . Focusing on the resonant profile of the recombination coefficient, Equation (15) 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 H L | a H H , and (ii) the modification of an Efimov resonance into a bound state embedded in the three-body and atom–dimer continua for | q | .

Author Contributions

All authors contributed equally to the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The work of CHG was supported in part by the U.S. National Science Foundation, Grant No. PHY-1912350. The numerical calculations were performed using NSF XSEDE Resource Allocation No. TG-PHY150003.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

The authors would like to thank M.T. Eiles for fruitful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nielsen, E.; Fedorov, D.V.; Jensen, A.S.; Garrido, E. The three-body problem with short-range interactions. Phys. Rep. 2001, 347, 373–459. [Google Scholar] [CrossRef]
  2. Greene, C.H.; Giannakeas, P.; Pérez-Ríos, J. Universal few-body physics and cluster formation. Rev. Mod. Phys. 2017, 89, 035006. [Google Scholar] [CrossRef] [Green Version]
  3. Naidon, P.; Endo, S. Efimov Physics: A Review. Rep. Prog. Phys. 2017, 80, 056001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. D’Incao, J.P. Few-Body Physics in Resonantly Interacting Ultracold Quantum Gases. J. Phys. B At. Mol. Opt. Phys. 2018, 51, 043001. [Google Scholar] [CrossRef] [Green Version]
  5. Efimov, V. Hard-core interaction and the three-nucleon problem. Sov. J. Nucl. Phys. 1970, 10, 62. [Google Scholar]
  6. Kraemer, T.; Mark, M.; Waldburger, P.; Danzl, J.G.; Chin, C.; Engeser, B.; Lange, A.D.; Pilch, K.; Jaakkola, A.; Nägerl, H.C.; et al. Evidence for Efimov quantum states in an ultracold gas of caesium atoms. Nature 2006, 440, 315–318. [Google Scholar] [CrossRef] [Green Version]
  7. Riisager, K. Nuclear halo states. Rev. Mod. Phys. 1994, 66, 1105–1116. [Google Scholar] [CrossRef]
  8. Braaten, E.; Hammer, H.W. Universality in few-body systems with large scattering length. Phys. Rep. 2006, 428, 259–390. [Google Scholar] [CrossRef] [Green Version]
  9. Rittenhouse, S.T.; von Stecher, J.; D’Incao, J.; Mehta, N.; Greene, C.H. The hyperspherical four-fermion problem. J. Phys. B At. Mol. Opt. Phys. 2011, 44, 172001. [Google Scholar] [CrossRef] [Green Version]
  10. Blume, D. Few-body physics with ultracold atomic and molecular systems in traps. Rep. Prog. Phys. 2012, 75, 046401. [Google Scholar] [CrossRef] [Green Version]
  11. Wang, Y.; D’Incao, J.P.; Esry, B.D. Chapter 1—Ultracold Few-Body Systems. In Advances in Atomic, Molecular, and Optical Physics; Ennio Arimondo, P.R.B., Lin, C.C., Eds.; Academic Press: Berlin/Heidelberg, Germany, 2013; Volume 62, 115p. [Google Scholar]
  12. Wang, Y.; Julienne, P.; Greene, C.H. Few-body physics of ultracold atoms and molecules with long-range interactions. In Annual Review of Cold Atoms and Molecules; Madison, K.W., Bongs, K., Carr, L.D., Rey, A.M., Zhai, H., Eds.; World Scientific Press: Singapore, 2015; Volume 3, pp. 77–134. [Google Scholar]
  13. Huang, B.; Sidorenkov, L.A.; Grimm, R.; Hutson, J.M. Observation of the second triatomic resonance in Efimov’s scenario. Phys. Rev. Lett. 2014, 112, 190401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Thomas, L.H. The interaction between a neutron and a proton and the structure of H3$. Phys. Rev. 1935, 47, 903–909. [Google Scholar] [CrossRef]
  15. Wang, Y.; Julienne, P.S. Universal van der Waals physics for three cold atoms near Feshbach resonances. Nat. Phys. 2014, 10, 768–773. [Google Scholar] [CrossRef] [Green Version]
  16. Roy, S.; Landini, M.; Trenkwalder, A.; Semeghini, G.; Spagnolli, G.; Simoni, A.; Fattori, M.; Inguscio, M.; Modugno, G. Test of the Universality of the Three-Body Efimov Parameter at Narrow Feshbach Resonances. Phys. Rev. Lett. 2013, 111, 053202. [Google Scholar] [CrossRef] [PubMed]
  17. Wang, J.; D’Incao, J.; Esry, B.; Greene, C.H. Origin of the three-body parameter universality in Efimov physics. Phys. Rev. Lett. 2012, 108, 263001. [Google Scholar] [CrossRef] [Green Version]
  18. Gross, N.; Shotan, Z.; Kokkelmans, S.; Khaykovich, L. Observation of universality in ultracold 7Li three-body recombination. Phys. Rev. Lett. 2009, 103, 163202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Naidon, P.; Endo, S.; Ueda, M. Microscopic Origin and Universality Classes of the Efimov Three-Body Parameter. Phys. Rev. Lett. 2014, 112, 105301. [Google Scholar] [CrossRef] [Green Version]
  20. Naidon, P.; Endo, S.; Ueda, M. Physical origin of the universal three-body parameter in atomic Efimov physics. Phys. Rev. A 2014, 90, 022106. [Google Scholar] [CrossRef] [Green Version]
  21. Ferlaino, F.; Grimm, R. Forty years of Efimov physics: How a bizarre prediction turned into a hot topic. Physics 2010, 3, 9. [Google Scholar] [CrossRef] [Green Version]
  22. Giannakeas, P.; Greene, C.H. Van Der Waals Universality in Homonuclear Atom-Dimer Elastic Collisions. Few-Body Syst. 2017, 58, 20. [Google Scholar] [CrossRef] [Green Version]
  23. Mestrom, P.M.A.; Wang, J.; Greene, C.H.; D’Incao, J.P. Efimov–van Der Waals Universality for Ultracold Atoms with Positive Scattering Lengths. Phys. Rev. A 2017, 95, 032707. [Google Scholar] [CrossRef] [Green Version]
  24. Mestrom, P.M.A.; Colussi, V.E.; Secker, T.; Groeneveld, G.P.; Kokkelmans, S.J.J.M.F. Van Der Waals Universality near a Quantum Tricritical Point. Phys. Rev. Lett. 2020, 124, 143401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Kunitski, M.; Zeller, S.; Voigtsberger, J.; Kalinin, A.; Schmidt, L.P.H.; Schoeffler, M.; Czasch, A.; Schoellkopf, W.; Grisenti, R.E.; Jahnke, T.; et al. Observation of the Efimov state of the helium trimer. Science 2015, 348, 551–555. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Pires, R.; Ulmanis, J.; Häfner, S.; Repp, M.; Arias, A.; Kuhnle, E.D.; Weidemüller, M. Observation of Efimov Resonances in a Mixture with Extreme Mass Imbalance. Phys. Rev. Lett. 2014, 112, 250404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Ulmanis, J.; Häfner, S.; Pires, R.; Kuhnle, E.D.; Weidemüller, M.; Tiemann, E. Universality of weakly bound dimers and Efimov trimers close to Li–Cs Feshbach resonances. New J. Phys. 2015, 17, 055009. [Google Scholar] [CrossRef] [Green Version]
  28. Tung, S.K.; Jiménez-García, K.; Johansen, J.; Parker, C.V.; Chin, C. Geometric Scaling of Efimov States in a 6Li-133$Cs. Phys. Rev. Lett. 2014, 113, 240402. [Google Scholar] [CrossRef] [Green Version]
  29. Häfner, S.; Ulmanis, J.; Kuhnle, E.D.; Wang, Y.; Greene, C.H.; Weidemüller, M. Role of the intraspecies scattering length in the Efimov scenario with large mass difference. Phys. Rev. A 2017, 95, 062708. [Google Scholar] [CrossRef] [Green Version]
  30. Petrov, D.S.; Werner, F. Three-body recombination in heteronuclear mixtures at finite temperature. Phys. Rev. A 2015, 92, 022704. [Google Scholar] [CrossRef] [Green Version]
  31. Zhao, C.Y.; Han, H.L.; Wu, M.S.; Shi, T.Y. Universal three-body parameter of heavy-heavy-light systems with a negative intraspecies scattering length. Phys. Rev. A 2019, 100, 052702. [Google Scholar] [CrossRef] [Green Version]
  32. Sun, M.; Liu, C.; Shi, Z.Y. Efimov physics in the complex plane. arXiv 2021, arXiv:2109.11206. [Google Scholar]
  33. Ulmanis, J.; Häfner, S.; Pires, R.; Kuhnle, E.D.; Wang, Y.; Greene, C.H.; Weidemüller, M. Heteronuclear Efimov Scenario with Positive Intraspecies Scattering Length. Phys. Rev. Lett. 2016, 117, 153201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Chin, C.; Grimm, R.; Julienne, P.; Tiesinga, E. Feshbach resonances in ultracold gases. Rev. Mod. Phys. 2010, 82, 1225. [Google Scholar] [CrossRef]
  35. Johansen, J.; DeSalvo, B.J.; Patel, K.; Chin, C. Testing universality of Efimov physics across broad and narrow Feshbach resonances. Nat. Phys. 2017, 13, 731–735. [Google Scholar] [CrossRef] [Green Version]
  36. Giannakeas, P.; Greene, C.H. Ultracold Heteronuclear Three-Body Systems: How Diabaticity Limits the Universality of Recombination into Shallow Dimers. Phys. Rev. Lett. 2018, 120, 023401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. D’Incao, J.P.; Esry, B.D. Scattering Length Scaling Laws for Ultracold Three-Body Collisions. Phys. Rev. Lett. 2005, 94, 213201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Hsu, C.W.; Zhen, B.; Stone, A.D.; Joannopoulos, J.D.; Soljačić, M. Bound States in the Continuum. Nat. Rev. Mater. 2016, 1, 16048. [Google Scholar] [CrossRef] [Green Version]
  39. Rittenhouse, S.T.; Mehta, N.P.; Greene, C.H. Green’s functions and the adiabatic hyperspherical method. Phys. Rev. A 2010, 82, 022706. [Google Scholar] [CrossRef] [Green Version]
  40. Avery, J. Hyperspherical Harmonics: Applications in Quantum Theory; Kluwer Academic Publishers: Norwell, MA, USA, 1989. [Google Scholar]
  41. Smirnov, Y.F.; Shitikova, K.V. Method of K harmonics and the shell model. Sov. J. Part. Nucl. 1977, 8, 44. [Google Scholar]
  42. Kartavtsev, O.I.; Malykh, A.V. Low-energy three-body dynamics in binary quantum gases. J. Phys. B 2007, 40, 1429–1441. [Google Scholar] [CrossRef] [Green Version]
  43. Kartavtsev, O.I.; Malykh, A.V. Universal low-energy properties of three two-dimensional bosons. Phys. Rev. A 2006, 74, 042506. [Google Scholar] [CrossRef] [Green Version]
  44. Aymar, M.; Greene, C.H.; Luc-Koenig, E. Multichannel Rydberg spectroscopy of complex atoms. Rev. Mod. Phys. 1996, 68, 1015. [Google Scholar] [CrossRef]
  45. Mehta, N.; Esry, B.; Greene, C.H. Three-body recombination in one dimension. Phys. Rev. A 2007, 76, 022711. [Google Scholar] [CrossRef] [Green Version]
  46. Burke, J.P., Jr. Theoretical Investigation of Cold Alkali Atom Collisions. Ph.D. Thesis, University of Colorado, Boulder, CO, USA, 1999. [Google Scholar]
  47. Nielsen, E.; Macek, J. Low-energy recombination of identical bosons by three-body collisions. Phys. Rev. Lett. 1999, 83, 1566. [Google Scholar] [CrossRef]
  48. Clark, C.W. Calculation of Non-adiabatic Transition Probabilities. Phys. Lett. A 1979, 70, 295. [Google Scholar] [CrossRef]
  49. Child, M. Semiclassical theory of tunneling and curve-crossing problems: A diagrammatic approach. J. Mol. Spectrosc. 1974, 53, 280–301. [Google Scholar] [CrossRef]
  50. Zhu, C.; Nakamura, H. Theory of nonadiabatic transition for general two-state curve crossing problems. I. Nonadiabatic tunneling case. J. Chem. Phys. 1994, 101, 10630–10647. [Google Scholar] [CrossRef]
  51. Wang, Y.; Wang, J.; D’Incao, J.P.; Greene, C.H. Universal Three-Body Parameter in Heteronuclear Atomic Systems. Phys. Rev. Lett. 2012, 109, 243201. [Google Scholar] [CrossRef] [Green Version]
  52. Fano, U. Effects of Configuration Interaction on Intensities and Phase Shifts. Phys. Rev. 1961, 124, 1866–1878. [Google Scholar] [CrossRef]
Figure 1. An illustration of the lowest hyperspherical potential curves U ν 1 / 3 ( R / a H H ) with a H H > 0 and a H L < 0 . The red (blue) line saturates at large hyperradii in the atom–dimer (three-body break-up) threshold. The quantities Φ L U and Φ L U 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 P 12 maximizes. The horizontal dotted line refers to the three-body collisional energy E ¯ in units of ħ 2 m H a H H 2 , and the three-body parameter, r 3 b a H H , depicted by the blue region.
Figure 1. An illustration of the lowest hyperspherical potential curves U ν 1 / 3 ( R / a H H ) with a H H > 0 and a H L < 0 . The red (blue) line saturates at large hyperradii in the atom–dimer (three-body break-up) threshold. The quantities Φ L U and Φ L U 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 P 12 maximizes. The horizontal dotted line refers to the three-body collisional energy E ¯ in units of ħ 2 m H a H H 2 , and the three-body parameter, r 3 b a H H , depicted by the blue region.
Atoms 09 00110 g001
Figure 2. The Stokes correction phase as a function of the non-adiabatic probability p.
Figure 2. The Stokes correction phase as a function of the non-adiabatic probability p.
Atoms 09 00110 g002
Figure 3. The degree of diabaticity p as a function of the scattering length ratio a H H / | a H L | for different mass ratios m H / m L , covering the regime of strong-to-weak mass-imbalanced three-body systems.
Figure 3. The degree of diabaticity p as a function of the scattering length ratio a H H / | a H L | for different mass ratios m H / m L , covering the regime of strong-to-weak mass-imbalanced three-body systems.
Atoms 09 00110 g003
Figure 4. The scaled | S 12 | 2 ( k a H L ) 4 matrix element versus the ratios | a H L | a H H and r 3 b a H H for the 6 Li - 87 Rb - 87 Rb system at low energy E = ħ 2 k 2 2 μ . (a) Semi-classical approach and (b) R-matrix numerical calculations.
Figure 4. The scaled | S 12 | 2 ( k a H L ) 4 matrix element versus the ratios | a H L | a H H and r 3 b a H H for the 6 Li - 87 Rb - 87 Rb system at low energy E = ħ 2 k 2 2 μ . (a) Semi-classical approach and (b) R-matrix numerical calculations.
Atoms 09 00110 g004
Figure 5. An illustration of the approximate hyperspherical potential curves shown in Figure 1, where s 0 and s 0 * are the universal Efimov scaling coefficients. These piecewise curves are used in Equations (12), (13) and (15).
Figure 5. An illustration of the approximate hyperspherical potential curves shown in Figure 1, where s 0 and s 0 * are the universal Efimov scaling coefficients. These piecewise curves are used in Equations (12), (13) and (15).
Atoms 09 00110 g005
Figure 6. In the limit of E 0 , the scaled recombination coefficient m H K 3 ħ a H L 4 is shown as a function of ln r 3 b a H H for (a) 6 Li - 133 Cs - 133 Cs and (b) 6 Li - 87 Rb - 87 Rb . The symbols refer to the corresponding calculations in the semi-classical approach. The solid lines indicate the fitting of Equation (15) using the universal parameters shown in Table 1.
Figure 6. In the limit of E 0 , the scaled recombination coefficient m H K 3 ħ a H L 4 is shown as a function of ln r 3 b a H H for (a) 6 Li - 133 Cs - 133 Cs and (b) 6 Li - 87 Rb - 87 Rb . The symbols refer to the corresponding calculations in the semi-classical approach. The solid lines indicate the fitting of Equation (15) using the universal parameters shown in Table 1.
Atoms 09 00110 g006
Figure 7. A comparison of the scaled recombination coefficient obtained via the fitting of Equation (15) (points) and the Fano lineshape formula (solid lines) from Equation (16) for two HHL systems. The red points and lines correspond to 6 Li - 133 Cs - 133 Cs for a scattering length ratio | a CsLi | a CsCs = 67.8 . The black points and lines denote the 6 Li - 87 Rb - 87 Rb system at | a RbLi | a RbRb = 126.2 . Note that the total colliding energy is set to zero.
Figure 7. A comparison of the scaled recombination coefficient obtained via the fitting of Equation (15) (points) and the Fano lineshape formula (solid lines) from Equation (16) for two HHL systems. The red points and lines correspond to 6 Li - 133 Cs - 133 Cs for a scattering length ratio | a CsLi | a CsCs = 67.8 . The black points and lines denote the 6 Li - 87 Rb - 87 Rb system at | a RbLi | a RbRb = 126.2 . Note that the total colliding energy is set to zero.
Atoms 09 00110 g007
Figure 8. Panels (ad) show the width of the Efimov resonance Γ and the asymmetry parameter q versus the scattering length ratio | a H L | a H H 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. Moreover, Γ and q are obtained via Equations (17) and (18), respectively, using the universal parameters shown in Table 1.
Figure 8. Panels (ad) show the width of the Efimov resonance Γ and the asymmetry parameter q versus the scattering length ratio | a H L | a H H 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. Moreover, Γ and q are obtained via Equations (17) and (18), respectively, using the universal parameters shown in Table 1.
Atoms 09 00110 g008
Table 1. A summary of the universal parameters used in Equations (16)–(18) for the systems of 6 Li - 133 Cs - 133 Cs and 6 Li - 87 Rb - 87 Rb . Note that the values of s 0 and s 0 * are calculated in Ref. [51].
Table 1. A summary of the universal parameters used in Equations (16)–(18) for the systems of 6 Li - 133 Cs - 133 Cs and 6 Li - 87 Rb - 87 Rb . Note that the values of s 0 and s 0 * are calculated in Ref. [51].
HHL System s 0 s 0 * γ ψ 1 ψ 2
6 Li - 133 Cs - 133 Cs 1.9832.0034.420.460.13
6 Li - 87 Rb - 87 Rb 1.6331.6823.130.80.4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Giannakeas, P.; Greene, C.H. Asymmetric Lineshapes of Efimov Resonances in Mass-Imbalanced Ultracold Gases. Atoms 2021, 9, 110. https://doi.org/10.3390/atoms9040110

AMA Style

Giannakeas P, Greene CH. Asymmetric Lineshapes of Efimov Resonances in Mass-Imbalanced Ultracold Gases. Atoms. 2021; 9(4):110. https://doi.org/10.3390/atoms9040110

Chicago/Turabian Style

Giannakeas, Panagiotis, and Chris H. Greene. 2021. "Asymmetric Lineshapes of Efimov Resonances in Mass-Imbalanced Ultracold Gases" Atoms 9, no. 4: 110. https://doi.org/10.3390/atoms9040110

APA Style

Giannakeas, P., & Greene, C. H. (2021). Asymmetric Lineshapes of Efimov Resonances in Mass-Imbalanced Ultracold Gases. Atoms, 9(4), 110. https://doi.org/10.3390/atoms9040110

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop