Next Article in Journal
Spectral Line Shapes in Plasmas II
Next Article in Special Issue
Universal Scattering of Ultracold Atoms and Molecules in Optical Potentials
Previous Article in Journal
Quasars: From the Physics of Line Formation to Cosmology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Two-Body and Three-Body Contacts for Three Bosons in the Unitary Regime: Analytic Expressions and Limiting Forms

by
Victor E. Colussi
Department of Applied Physics, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
Atoms 2019, 7(1), 19; https://doi.org/10.3390/atoms7010019
Submission received: 19 December 2018 / Revised: 30 January 2019 / Accepted: 2 February 2019 / Published: 5 February 2019
(This article belongs to the Special Issue Few-body Physics in Ultracold Quantum Gases)

Abstract

:
The two- and three-body contacts are central to a set of univeral relations between microscopic few-body physics within an ultracold Bose gas and its thermodynamical properties. They may also be defined in trapped few-particle systems, which is the subject of this work. In this work, we focus on the unitary three-body problem in a trap, where interactions are as strong as allowed by quantum mechanics. We derive analytic results for the two- and three-body contacts in this regime and compare them with existing limiting expressions and previous numerical studies.

1. Introduction

As the s-wave scattering length approaches unitarity | a | , properties of both macroscopic and microscopic strongly-interacting ultracold quantum systems can be described effectively by a reduced set of remaining finite quantities. For three trapped bosonic alkali atoms at unitarity, physics in the unitary regime is parametrized solely by the range of the potential, captured by the van der Waals length r vdW [1], and the trap length a ho = ( / m ω ) 1 / 2 , with frequency ω and single-particle mass m. Consequently, these scales parametrize the two and three-body contacts, respectively, that are central to many relevant observables in the system including the tail of the single-particle density, short-distance behavior of correlation functions, high-frequency tail of the rf transition rate, virial theorem, and total energy of the system [2,3]. Intuitively, the two and three-body contacts measure the mean number of clustered pairs and triplets of bosons, respectively [4]. The centrality of the extensive two-body contact C 2 to the thermodynamic properties of the two-component unitary Fermi gas was worked out by Shina Tan [5,6,7], and this quantity frames part of our current understanding of experiments in this regime [8]. For strongly-interacting Bose gases, however, the presence of the Efimov effect [9,10] and associated infinite sequence of bound trimers require the introduction of an additional parameter, the extensive three-body contact C 3 [2,3]. In the strongly-interacting regime, where interaction length scales are diverging, the contacts may also be used to interpolate in between perturbative results in neighboring regimes [11].
The centrality of the contacts means that they may be measured through various means. By measuring the high-momentum tail of the rf transition rate in a weakly-interacting Bose gas, Wild et al. [12] measured a nonzero C 2 and found that C 3 was consistent with zero. In the non-equilibrium regime of a nondegenerate Bose gas quenched to unitarity, both C 2 and C 3 were measured interferrometrically in Ref. [13], finding a gradual approach of C 3 to equilibrium. It has also been suggested that the high-momentum tail of the degenerate Bose gas quenched to unitarity measured in Ref. [14] is a nonzero measurement of C 3 [15], although there is currently no consensus and a recent experimental study by Eigen et al. [16] found no evidence of the contacts in the single-particle momentum distribution at early-times after the quench. Theoretically, a prediction for the gradual growth of C 3 in a uniform degenerate Bose gas quenched to unitary was made in Ref. [17], where a coherent beating phenomenon at the frequency of Efimov states was observed.
The contacts can also be measured in microscopic few-body systems, for instance, by measuring interaction energies or spectroscopically. Theoretically, the contacts corresponding to the lowest few eigenstates were calculated numerically for three and four trapped bosons in Ref. [18] at fixed trap length. In the present work, we provide analytic results for C 2 and C 3 for the trapped unitary three-body problem including analytic expressions and limiting forms. This includes a review of some of previous results derived already in the supplementary materials of Ref. [17] for C 3 , included here in more detail, in addition to new derivations for C 2 .
In this work, we begin in Section 2 by reviewing analytic solutions of the unitary three-boson problem in a trap in the zero-range model first derived in Ref. [19]. In Section 3, we use a relation between the short-distance behavior of the three-body correlation function and C 3 from Ref. [17] to derive analytic expressions for the three-body contacts for each trapped three-body eigenstate including limiting expressions in the large-energy limit. In Section 4, we perform an analogous derivation except for the two-body contacts for each trapped three-body eigenstate. Finally, in Section 5, we discuss our findings and compare against the available existing results in [18], which evaluated C 2 and C 3 including more realistic finite-range effects at fixed trap length. We find generally excellent agreement with the available results of that work.

2. Three Trapped Bosons at Unitarity

In this section, we outline analytic solutions of the unitary three-body problem in a trap given in Ref. [19]. Consider three non-interacting bosons located at r 1 , r 2 , r 3 in a harmonic trap whose wave function Ψ ( r 1 , r 2 , r 3 ) solves the non-interacting Schrödinger equation
l = 1 3 2 2 m r l 2 + 1 2 m ω 2 r i 2 Ψ = E Ψ .
The three-body wave function can be separated | Ψ = | Ψ CM | ψ into center of mass | Ψ CM and relative | ψ components. The center of mass coordinate is defined as C ( r 1 + r 2 + r 3 ) / 3 , and the relative motion is parametrized by the hyperradius R 2 ( r 2 + ρ 2 ) / 2 and hyperangles Ω { ρ ^ , r ^ , α } , where r r 2 r 1 and ρ = ( 2 r 3 r 1 r 2 ) / 3 are the Jacobi vectors with spherical angles ρ ^ and r ^ , respectively. The hyperangle α is defined as tan α = r / ρ , and is related to each individual Jacobi vector as r = 2 R sin α , ρ = 2 R cos α where α [ 0 , π / 2 ] .
In ultracold quantum gases of alkali atoms, it is generally the case that the typical momentum of the problem are such that k r vdW 1 , where r vdW is the van der Waals length which furnishes a natural short-range for pairwise interactions [1]. In this limit, the zero-range approximation can be made, and pairwise interactions can be accounted for via the Bethe-Peierls boundary conditions
ψ ( R , Ω ) ) = | r | 0 1 a 1 r A ( ρ ) ,
where the unknown function A is determined by solving the three-body Schrödinger equation (Equation (1)) with the above boundary condition.
In the unitary regime, | a | in Equation (2) and the relative three-body wave function becomes singular as 1 / r in the r 0 limit. Following Ref. [19] and the original approach of Efimov [9], the relative three-body wave function in this regime can be decomposed as ψ ( R , Ω ) = R 2 N F ( R ) ϕ ( Ω ) in terms of hyperradial F ( R ) and hyperrangular ϕ ( Ω ) ) = ( 1 + Q ^ ) φ ( α ) / sin 2 α 4 π functions with normalization constant N . The operator Q ^ = P ^ 13 + P ^ 23 is written in terms of the permuation operator P ^ i j that swaps particles i and j. The hyperangular basis functions φ are determined from solutions of the equations [19]
φ s ( α ) = s 2 φ s ( α ) ,
φ s ( π / 2 ) = 0 ,
φ s ( 0 ) + 8 3 φ s ( π / 3 ) = 0 ,
which are solved by φ s ( α ) = sin ( s ( π / 2 α ) ) . The channel index s is obtained by solving the transcendental equation
8 3 sin s π 6 s cos s π 2 = 0 .
Each channel supports an associated set of hyperradial eigenfunctions F j ( s ) ( R ) satisfying
2 2 m d 2 d R 2 + 1 R d d R + 2 s 2 2 m R 2 + m ω 2 R 2 2 F j ( s ) ( R ) = E 3 b ( s , j ) F j ( s ) ( R ) ,
with quantum number j = 0 , 1 , 2 , and relative three-body eigenenergy E 3 b ( s , j ) .
The transcendental equation (Equation (6)) has both real and imaginary solutions. For s 2 > 0 , Equation (7) has the solution
F j ( s ) ( R ) = e R 2 / 2 a ho 2 L j ( s ) ( R 2 / a ho 2 ) R a ho s ,
where L j ( s ) is a generalized Laguerre polynomial of degree j, and the associated eigenenergies are E 3 b ( s , j ) = ω ( s + 2 j + 1 ) . We label the channels with s 2 > 0 as universal because they depend solely on the trapping parameters and not on microscopic details of the interaction.
The lone imaginary root of Equation (6) s = i s 0 denotes the Efimov channel, where s 0 1.00624 is Efimov’s constant [9]. Equation (7) must be supplemented with the additional boundary condition F j ( s 0 ) ( R ) sin [ s 0 ln ( R / R t ) ] to preserve the Hermiticity of the problem [19,20], where R t is a three-body parameter setting the phase of the log-periodic oscillation. The hyperradial solutions in this channel are
F j ( s 0 ) ( R ) = a ho R W E 3 b ( s 0 , j ) / 2 ω , i s 0 / 2 ( R 2 / a ho 2 ) ,
where W is a Whittaker function [21]. The eigenenergy spectrum in the Efimov channel is obtained by solving
arg Γ 1 + i s 0 E 3 b ( s 0 , j ) / ω 2 + s 0 ln R t a ho = arg Γ 1 + i s 0 ,
which is understood mod π and was first obtained in Ref. [22]. To make the connection between the three-body parameters R t and κ * , we take the j limit of Equation (10), using the relations [21]
arg Γ 1 + i s 0 E 3 b ( s 0 , j ) / ω 2 = s 0 2 ψ 1 E 3 b ( s 0 , j ) / ω 2 , s 0 2 ln E 3 b ( s 0 , j ) / ω 2 ,
where ψ in the above equation is the digamma function. Inserting this limiting result into Equation (10) gives the free-space result for the binding energies of Efimov trimers
E 3 b ( s 0 , j ) = j 2 2 m R t 2 e 2 s 0 ( arg Γ ( 1 + i s 0 ) π j ) .
Setting j = 0 in the above equation gives R t = 2 exp ( Im ln [ Γ ( 1 + i s 0 ) ] / s 0 ) / κ * in terms of the wave number of the j = 0 trimer E 3 b ( s 0 , 0 ) = 2 κ * 2 / m . We now choose the three-body parameter R t such that E 3 b ( s 0 , 0 ) matches the ground-state trimer energy at unitarity for neutral atoms studied in Ref. [23]. In that work, the universality of κ * with the van der Waals length was established, finding κ * = 0.226 / r vdW . Therefore, we make the arbitrary distinction that j = 0 is the “ground-state” Efimov trimer in the sense of Ref. [23] and j < 0 comprises the remainder of the Thomas collapse [24]. Taking the j limit of the Whittaker function in Equation (9) gives
F j s 0 ( R ) = j K i s 0 ( 2 R | E 3 b ( s 0 , j ) | / ω / a ho ) ,
where K i s 0 is a modified Bessel function of imaginary argument, which is the free-space bound-state wave function for an Efimov trimer [9]. The Efimov channel depends on the trapping parameters and also on R t , which is sensitive to the microscopic length scale r vdW . The Efimov channel is therefore not universal. Furthermore, although R t depends on r vdW , the model remains zero-range involving only contact potentials. We contrast this with the approach of Ref. [18], which used more realistic finite-ranged potentials.
For completeness, we discuss also the j + limit of Equation (10), where Stirling’s formula gives
arg Γ 1 + i s 0 E 3 b ( s 0 , j ) / ω 2 s 0 2 ln ( E 3 b ( s 0 , j ) / ω ) π E 3 b ( s 0 , j ) 2 ω .
Therefore, Equation (10) in this limit becomes
E 3 b ( s 0 , j ) / ω = 2 j + 1 π s 0 ln ( R t 2 / a ho 2 ) 2 arg Γ ( 1 + i s 0 ) + s 0 π ln ( E 3 b ( s 0 , j ) / ω ) ,
which is in general difficult to solve. However, Equation (15) is written in a form “ f ( x ) = a + g ( f ( x ) ) ” that can be iterated. In the E 3 b ( s 0 , j ) limit, this produces corrections in powers of j as
E 3 b ( s 0 , j ) / ω = j + 2 j + 1 π s 0 ln ( j R t 2 / a ho 2 ) 2 arg Γ ( 1 + i s 0 ) + O ( j 1 ln j ) ,
which was first obtained in Ref. [25].
The normalization constants N s , j are defined as
N s , j 2 2 3 3 1 F j ( s ) | F j ( s ) ϕ s | ϕ s ,
where the hyperangular inner product is defined in terms of the hyperangular solid angle as
ϕ s | ϕ s d Ω | ϕ s ( Ω ) | 2 ,
= 12 π s sin s * π 2 cos s π 2 s π 2 sin s π 2 4 π 3 3 cos s π 6 ,
d Ω 2 d α sin 2 2 α 0 π / 2 d ρ ^ d r ^ .
The result in Equation (18) was first calculated in Ref. [25] using a change of variables trick due to Efimov in Ref. [26]. The hyperradial inner product is defined as
F j ( s ) | F j ( s ) 0 d R R | F j ( s ) ( R ) | 2 ,
F j ( s ) | F j ( s ) = a ho 2 Γ s + 1 + j 2 j ! ( s 2 > 0 ) ,
F j ( s 0 ) | F j ( s 0 ) = a ho 2 π · Im ψ 1 E 3 b ( s 0 , j ) / ω + i s 0 2 sinh ( s 0 π ) · Γ 1 E 3 b ( s 0 , j ) / ω + i s 0 2 2 ,
F j ( s 0 ) | F j ( s 0 ) = j a ho 2 π s 0 4 | E 3 b ( s 0 , j ) / ω | sinh ( π s 0 ) ,
where the relevant integrals can be found tabulated in Ref. [27]. The factor ( 2 / 3 ) 3 in Equation (17) is due to the Jacobian (D) of the transformation D ( r 1 , r 2 , r 3 ) / D ( C , R , Ω ) = ( 3 / 2 ) 3 , required so that Ψ | Ψ = Ψ CM | Ψ CM ψ s , j | ψ s , j = 1 .

3. Three-Body Contact C 3

Any one of the universal relations involving the contacts can be used as a starting point; however, we choose to obtain them through the limiting behavior of few-body correlation functions [2,3]. We begin by deriving the extensive three-body contacts C 3 for each trapped three-body eigenstate.
The non-normalized triplet correlation function is defined in second-quantization as [28]
G ( 3 ) ( r 1 , r 2 , r 3 ) ψ ^ ( r 1 ) ψ ^ ( r 2 ) ψ ^ ( r 3 ) ψ ^ ( r 3 ) ψ ^ ( r 2 ) ψ ^ ( r 1 ) ,
in terms of the bosonic field operators ψ ^ ( r ) = k a ^ k e i k · r / V and system volume V. For three bosons in vacuum, this is equivalent in first quantization to
G ( 3 ) ( r 1 , r 2 , r 3 ) = 3 ! | Ψ ( r 1 , r 2 , r 3 ) | 2 .
Following [2], we integrate out the center of mass dependence, and take the limit R 0 at fixed Ω
3 ! d 3 C Ψ ( r 1 , r 2 , r 3 ) = 3 ! | ψ s , j ( R , Ω ) | 2 , = R 0 | ψ s c ( R , Ω ) | 2 8 3 s 0 2 C 3 ( s , j ) ,
for the jth eigenstate in the Efimov channel. In the above equation, we have indexed the extensive three-body contact specific to the eigenstate ( s , j ) as C ( s , j ) . In Equation (27), ψ s c ( R , Ω ) is zero-energy three-body scattering state [2]
ψ s c ( R , Ω ) = 1 R 2 sin s 0 ln R R t ϕ s 0 ( Ω ) ϕ s 0 | ϕ s 0 .
Integrating both sides of Equation (27) over d Ω , and taking advantage of the orthogonality of the hyperangular eigenfunctions so that only s = i s 0 needing to be considered yields
s 0 2 4 F j ( s 0 ) | F j ( s 0 ) F j ( s 0 ) ( R ) sin [ s 0 ln R / R t ] 2 = R 0 C 3 ( s 0 , j ) .
From the asymptotic behavior of the Whittaker functions [21], the R 0 limit in the above equation can be taken with the simple result
a ho 2 C 3 ( s 0 j ) = s 0 Im ψ 1 E 3 b ( s 0 , j ) / ω + i s 0 2 .
In the above equation, C 3 is made dimensionless by rescaling in powers of the oscillator length. In the following two subsections, we derive first the | E 3 b | / ω limiting forms of Equation (30) and then show some results at intermediate energies.

3.1. Limiting Cases

In this subsection, we derive the | E 3 b | / ω limiting forms of Equation (30).
Case I ( E 3 b / ω ) : In this limit, the digamma function has the asymptotic form [21]
ψ ( z ) = ln z + O ( z 1 ) ( z with | arg z | < π ) .
Using the property ln z = ln | z | + i arg z , we find
Im ψ 1 E 3 b ( s 0 , j ) / ω + i s 0 2 = arg 1 E 3 b ( s 0 , j ) / ω + i s 0 , = arctan s 0 1 E 3 b ( s 0 , j ) / ω . s 0 | E 3 b ( s 0 , j ) | / ω .
We therefore obtain the result for
C 3 ( s 0 , j ) m | E 3 b ( s 0 , j ) | 2 , κ * 2 e 2 π j / s 0 , ( E 3 b ( s 0 , j ) / ω ) ,
where we have used the result of Equation (12). The above result matches the operational definition κ * κ * E 3 b ( s , j ) | a = 2 2 C 3 ( s , j ) / m in the free-space limit [2,3].
Case II ( E 3 b / ω + ) : In this limit, the asymptotic expansion of the digamma function (Equation (31)) should be used carefully because we are taking the limit | arg z | π . Instead, we use the reflection formula for the digamma function [21] to obtain
Im ψ 1 E 3 b ( s 0 , j ) / ω + i s 0 2 = π Im cot π 2 ( 1 + i s 0 E 3 b ( s 0 , j ) / ω ) + Im ψ 1 + E 3 b ( s 0 , j ) / ω + i s 0 2 .
From the discussion of Case I, we know that the final term on the right in the above equation will vanish in the limit E 3 b / ω + . Therefore, if the term involving the imaginary part of the cotangent is nonvanishing, it is the only term that needs to be considered. Following this reasoning and using properties of cotangent for imaginary argument [21], we find that
π Im cot π 2 ( 1 + i s 0 E 3 b ( s 0 , j ) / ω ) = π sinh ( π s 0 ) cosh ( π s 0 ) + cos ( π E 3 b ( s 0 , j ) / ω ) ) ,
which is nonvanishing. We now take the limit ( E 3 b / ω + ) and use Equation (16) to find
a ho 2 C 3 ( s 0 , j ) s 0 π cosh ( π s 0 ) + cos ( s 0 ln ( j R t 2 / a ho 2 ) 2 arg Γ ( 1 + i s 0 ) ) sinh ( π s 0 ) , s 0 π 1 + 0.085 cos ( s 0 ln ( j R t 2 / a ho 2 ) 2 arg Γ ( 1 + i s 0 ) ) . ( E 3 b / ω )
The leading order result a ho 2 C 3 ( s 0 , j ) s 0 / π was obtained previously in Ref. [25] for an analogous expression for the decay rate Γ in this limit. Having outlined the limiting behavior of C 3 ( s 0 , j ) in Section 3.1, we now discuss its behavior in general and over a range of a ho κ * as shown in Figure 1. We find that C 3 ( s 0 , j ) matches the limiting forms (Equations (33) and (36)) by | E 3 b ( s 0 , j ) | / ω 10 2 as shown in Figure 1a. At intermediate energies, as a ho decreases, Efimov trimers whose free-space energies are comparable to ω are shifted above the zero-point energy of the trap as studied in Ref. [29]. This process repeats log-periodically as a ho κ * 0 . A single log-period is shown in Figure 1a, where the shift of the j = 2 Efimov trimer is clearly visible. For | E 3 b | / ω < 1 , we find that a ho 2 C 3 as a function of E 3 b / ω interpolates smoothly between the limiting behaviors as shown in Figure 1b.

3.2. Results

In general, we find that C 3 ( s 0 , j ) for a particular Efimov trimer is increased compared to the free space limit (see Equation (33)). We echo the conclusion of Ref. [18] that the restriction of a trimer to a reduced region of space should correspond to an increase in the probability to find three bosons at short distances as encoded in C 3 ( s 0 , j ) .
Additionally, the phase of the log-periodic oscillation of C 3 ( s 0 , j ) for E 3 b ( j ) / ω 1 appears to adjust. To see this behavior as a function of a ho κ * , we rewrite Equation (36) as
a ho 2 C 3 ( s 0 , j ) s 0 π 1 + 0.085 cos ( s 0 ln ( j ) 2 s 0 ln ( κ * a ho ) + Λ ) , ( E 3 b / ω )
where the phase factor Λ 0.697 . The factor 2 s 0 ln ( κ * a ho ) = 2 s 0 ln ( κ * a ho e π / s 0 ) mod 2 π , which gives the log-periodic behavior seen in Figure 1.

4. Two-Body Contact C 2

In this section, we derive the extensive two-body contacts C 2 for each trapped three-body eigenstate. Analogous to the treatment in Section 3, we begin from the definition of the non-normalized pair correlation function [28]
G ( 2 ) ( r 1 , r 2 ) ψ ^ ( r 1 ) ψ ^ ( r 2 ) ψ ^ ( r 2 ) ψ ^ ( r 1 ) .
For three bosons in vaccum, this is equivalent in first quantization to
G ( 2 ) ( r 1 , r 2 ) = 3 ! d 3 r 3 | Ψ ( r 1 , r 2 , r 3 ) | 2 .
We then integrate over the center of mass coordinate c ( r 1 + r 2 ) / 2 , and do a change of variables { c , r 3 } { C , r 3 , 12 } where r 3 , 12 r 3 c . The Jacobian for this transformation is unity, and we obtain
d 3 c G ( 2 ) ( r 1 , r 2 ) = 3 ! d 3 r 3 , 12 | ψ ( r , r 3 , 12 ) | 2 .
In the limit | r | 0 , the left-hand side of the above equation is related to C 2 [2], and therefore we find the following relation between the relative three-body wave function and the extensive two-body contact
3 ! d 3 r 3 , 12 | ψ ( r , r 3 , 12 ) | 2 = | r | 0 C 2 ( 4 π ) 2 r 2 .
In order to take the above limit, we rewrite the three-body relative wave function as
Ψ ( r , r 3 , 12 ) = χ 0 ( r , ρ ) r ρ + χ 0 ( r ( 2 ) , ρ ( 2 ) ) r ( 2 ) ρ ( 2 ) + χ 0 ( r ( 3 ) , ρ ( 3 ) ) r ( 3 ) ρ ( 3 ) ,
where the functions χ 0 have the important property of being finite for | r | 0 and vanishing as | ρ | 0 [30]. The vector superscripts indicate different Jacobi trees interrelated through the kinematic rotations
r ( j ) = r ( i ) cos γ ( k ) + ρ ( i ) sin γ ( k ) , ρ ( j ) = r ( i ) sin γ ( k ) ρ ( i ) cos γ ( k ) , γ ( k ) = ϵ i j k π / 3 ,
where we have implicitly defined r r ( 1 ) , ρ ρ ( 1 ) . With some rearrangement and the help of the kinematic rotations, we take the limit in Equation (41) to obtain
C 2 = 3 ! ( 4 π ) 2 lim | r | 0 d 3 r 3 , 12 r 2 χ 0 ( r , ρ ) r ρ + χ 0 ( r ( 2 ) , ρ ( 2 ) ) r ( 2 ) ρ ( 2 ) + χ 0 ( r ( 3 ) , ρ ( 3 ) ) r ( 3 ) ρ ( 3 ) 2 , = 3 ! ( 4 π ) 2 d 3 r 3 , 12 χ 0 ( 0 , ρ ) ρ 2 .
Written in terms of the trapped eigenstates outlined in Section 2, this becomes
C 2 ( s , j ) = 72 π 2 | N s , j | 2 | sin ( s π / 2 ) | 2 A j , j s , s ,
where we have indexed the extensive two-body contact specific to the eigenstate ( s , j ) as C 2 ( s , j ) . The notation A j , j s , s is shorthand for diagonal components of the array
A j , j s , s 3 2 d z F j s ( z ) [ F j ( s ) ( z ) ] * .
It is possible to obtain analytic expressions for the array A j , j s , s using the trapped hyperradial eigenstates outlined in Section 2. We are concerned however only with the components with j = j and s = s required to evaluate C 2 ( s , j ) via Equation (45). For s = i s 0 , Equation (46) is an integral over Whittaker functions that can be found in a table [27] with result
A j , j s 0 , s 0 = 3 2 a ho Re Γ [ 1 / 2 ] Γ [ 1 / 2 + i s 0 ] Γ [ i s 0 ] Γ 1 i s 0 E 3 b ( s 0 , j ) / ω 2 Γ 2 + i s 0 E 3 b ( s 0 , j ) / ω 2 × 3 F 2 1 2 , 1 2 + i s 0 , 1 E 3 b ( s 0 , j ) / ω + i s 0 2 ; 1 + i s 0 , 2 E 3 b ( s 0 , j ) / ω + i s 0 2 ; 1 ,
where 3 F 2 is a generalized hypergeometric function. We note that an equivalent expression to Equation (47) was first obtained in Ref. [25] for the loss-rate Γ . For channels with s 2 > 0 , Equation (46) is an integral over generalized Laguerre polynomials that can be recast using the recurrence relation [21]
L j ( s ) ( z ) = l = 0 j ( 2 s 1 ) / 2 + j l j l L l 1 / 2 ( z ) ,
where the · is the generalized binomial coefficient. Using this relation, the array A j , j s , s can be found tabulated in Ref. [27], and we find that
A j , j s , s = 1 2 3 2 a ho 2 s 1 2 ! 1 / 2 j s + j 1 / 2 j 3 F 2 j , j s , 1 / 2 ; j s + 1 / 2 , j + 1 / 2 ; 1 .
In the following subsection, we derive some limiting cases of A j , j s , s , and then discuss results for C 2 ( s , j ) in general.

4.1. Limiting Cases

Due to the 3 F 2 function in Equations (47) and (49), it is unclear how to derive limiting expressions for C 2 ( s , j ) in the limit j + [31]. It is, however, possible to derive expressions in the opposite limits, which includes the j limit in the Efimov channel and the case j = 0 for the universal channels.
For s = i s 0 in the j limit, the 3 F 2 function can be expanded up to order ( E 3 b ( s 0 , j ) ) m where m is any nonnegative integer [31]. Alternatively, using the asymptotic form of the hyperradial wave function in this limit (Equation (13)), the relevant integral for A j , j s 0 , s 0 can be found tabulated in Ref. [27], and we obtain
A j , j s 0 , s 0 = j a ho 3 π 2 8 | E 3 b ( s 0 , j ) | / ω cosh π s 0 ,
which gives the following expression for C 2 ( s , j )
a ho C 2 ( s 0 , j ) = j 144 π 2 | E 3 b ( s 0 , j ) | / ω sinh ( π s 0 / 2 ) tanh ( π s 0 ) 9 π s 0 sinh ( π s 0 / 2 ) 8 3 π cosh ( π s 0 / 6 ) + 18 cosh ( π s 0 / 2 ) , C 2 ( s 0 , j ) 53.097 κ * e π j / s 0 . ( E 3 b / ω )
This matches the result in Ref. [32] for C 2 for an Efimov trimer obtained via an alternative approach by considering the high-momentum tail of the single-particle momentum distribution.
For s 2 > 0 with j = 0 , the array A 0 , 0 s , s has the form
A 0 , 0 s , s = a ho 1 2 3 2 ( s 1 / 2 ) ! ,
and we obtain the following expression
a ho C 2 ( s , 0 ) = 144 2 π Γ ( s + 1 / 2 ) Γ ( s ) sin ( π s / 2 ) 8 3 π cos ( π s / 6 ) 18 cos ( π s / 2 ) + 9 π s sin ( π s / 2 ) ,
= s + 16 2 s 1 / 2 ,
which vanishes in the limit s + .

4.2. Results

Having outlined some of the limiting behaviors of C 2 ( s , j ) in Section 4.1, we now analyze its behavior in general over a range of universal channels and for various values of a ho κ * in the Efimov channel as shown in Figure 2. First, we discuss results for the universal channels shown in Figure 2a,b. In Figure 2a, we see that our asymptotic prediction for s + with j = 0 given by Equation (53) matches our calculation of a ho C 2 as s n 1 (solid black line). At larger j, the behavior of a ho C 2 is shown in Figure 2b, where we see a gradual decay. The j + asymptotic form of a ho C 2 was not derived in Section 4.1, although we estimate its form by fitting a ho C 2 to a power law j α . We find that, for extremely large energies E 3 b / ω > 10 6 , this exponent approaches α 0.5 in both universal and Efimov channels. This power law can be motivated by approximating 1 / 2 j s j s in the argument of the 3 F 2 function in Equation (49). The 3 F 2 function then reduces to 2 F 1 [31] and the j + limit can be taken to obtain a j 1 / 2 power law scaling for a ho C 2 .
Our results for a ho C 2 in the Efimov channel are shown in Figure 2c,d for both positive and negative energies starting at j = 0 . We see that the E 3 b / ω asymptotic form given in Equation (51) (dot-dashed black line) matches our calculation of a ho C 2 as shown in Figure 2c. At positive energies, a gradual decay of a ho C 2 is evident as in the universal case (Figure 2c), and the log-periodic oscillations that were prominent in C 3 in Figure 1 are absent. For varying a ho κ * , our results follow the same trends. As j + , we observe numerically the same j 1 / 2 power discussed for the universal channels. Additionally, in Figure 2d, we show the behavior of a ho C 2 s 0 , j as a function of E 3 b / ω in the region | E 3 b | / ω < 1 . Here, we find, as in Figure 1b, a smooth interpolation between limiting behaviors.
In general, we find that, over the same energy ranges, a ho C 2 is larger in the Efimov channel than in the universal channels. Physically, we understand this as a larger probability for finding correlated pairs of atoms at short-distances, which is measured by C 2 [4]. This should be larger in the absence of a repulsive hyperradial barrier for s 2 > 0 (see Equation (7)). As s increases, this barrier becomes more repulsive (see Equation (7)), and it is therefore less likely to find correlated pairs at short-distances. This is reflected in Figure 2a,b. When comparing to the free-space result in Equation (51) for an Efimov trimer, we find that the trap leads to a relative increase of C 2 . Tighter traps increasingly confine the atoms, and therefore it makes sense that the probability to find correlated pairs should increase accordingly.

5. Discussion

In this work, analytic solutions of the trapped unitary three-boson problem from Ref. [19] are utilized to derive analytic expressions for the two and three-body contacts for each trapped three-body eigenstate. These results can be used along with the set of universal relations associated with the contacts [2,3] to predict, for instance, the high-frequency tail of the rf transition rate, the virial theorem, loss-rates, and total energy. We derive the large energy asymptotics of the contacts and compare them with existing results in the literature when possible, finding good agreement. At intermediate energies, the only results in the literature that we are aware of are from Ref. [18], which were obtained using more realistic Gaussian potentials with a finite-range. We take a wave function based approach to obtain the contacts from the short-distance behavior of the pair and triplet correlation functions. In Ref. [18], however, the contacts were obtained from their operational definitions [2,3]
a a E 3 b ( s , j ) | κ * = 2 8 π m a C 2 ( s , j ) ,
κ * κ * E 3 b ( s , j ) | a = 2 2 m C 3 ( s , j ) ,
by taking the derivatives of three-body eigenenergies with variations of a and κ * . This amounted to changing the underlying parameters of the finite-range Gaussian potentials and approximating the derivatives via finite differencing. In Table 1, we compare against results in that work, finding good agreement for the eigenenergies and contacts to within a few percent or less. In Ref. [18], the convergence of eigenenergies was estimated as 10 3 E ho or better. The eigenenergies calculated in the present work can be calculated from Equation (10) to arbitrary precision. We attribute the disagreement therefore to the absence of finite-range effects in the zero-range model in the present work.
Finally, we note that the method outlined in these notes to obtain C 2 and C 3 through the limiting behavior of the correlation functions applies also away from unitarity at finite scattering lengths provided the zero-range approximation can be made [33]. Additionally, the derivation of C 2 from a three-body model presented in Section 4 may be generalized in the spirit of Ref. [17] to study the dynamics of the two-body contact following an interaction quench. We address this latter problem elsewhere [34] and leave the former as the subject of future study.

Funding

This research was funded by the Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO) Grant No. 680-47-623.

Acknowledgments

We acknowledge discussions with John Corson, José D’Incao, and Servaas Kokkelmans.

Conflicts of Interest

The author declares no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Chin, C.; Grimm, R.; Julienne, P.; Tiesinga, E. Feshbach resonances in ultracold gases. Rev. Mod. Phys. 2010, 82, 1225–1286. [Google Scholar] [CrossRef]
  2. Werner, F.; Castin, Y. General relations for quantum gases in two and three dimensions. II. Bosons and mixtures. Phys. Rev. A 2012, 86, 053633. [Google Scholar] [CrossRef]
  3. Braaten, E.; Kang, D.; Platter, L. Universal Relations for Identical Bosons from Three-Body Physics. Phys. Rev. Lett. 2011, 106, 153005. [Google Scholar] [CrossRef] [PubMed]
  4. Braaten, E.; Kang, D.; Platter, L. Universal relations for a strongly interacting Fermi gas near a Feshbach resonance. Phys. Rev. A 2008, 78, 053606. [Google Scholar] [CrossRef]
  5. Tan, S. Energetics of a strongly correlated Fermi gas. Ann. Phys. 2008, 323, 2952–2970. [Google Scholar] [CrossRef]
  6. Tan, S. Large momentum part of a strongly correlated Fermi gas. Ann. Phys. 2008, 323, 2971–2986. [Google Scholar] [CrossRef]
  7. Tan, S. Generalized virial theorem and pressure relation for a strongly correlated Fermi gas. Ann. Phys. 2008, 323, 2987–2990. [Google Scholar] [CrossRef]
  8. Zwerger, W. The BCS-BEC Crossover and the Unitary Fermi Gas; Springer: New York, NY, USA, 2011; Volume 836. [Google Scholar]
  9. Efimov, V. Weakly-bound states of three resonantly-interacting particles. Sov. J. Nucl. Phys 1971, 12, 101. [Google Scholar]
  10. 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]
  11. Laurent, S.; Pierce, M.; Delehaye, M.; Yefsah, T.; Chevy, F.; Salomon, C. Connecting Few-Body Inelastic Decay to Quantum Correlations in a Many-Body System: A Weakly Coupled Impurity in a Resonant Fermi Gas. Phys. Rev. Lett. 2017, 118, 103403. [Google Scholar] [CrossRef] [PubMed]
  12. Wild, R.J.; Makotyn, P.; Pino, J.M.; Cornell, E.A.; Jin, D.S. Measurements of Tan’s Contact in an Atomic Bose-Einstein Condensate. Phys. Rev. Lett. 2012, 108, 145305. [Google Scholar] [CrossRef] [PubMed]
  13. Fletcher, R.J.; Lopes, R.; Man, J.; Navon, N.; Smith, R.P.; Zwierlein, M.W.; Hadzibabic, Z. Two- and three-body contacts in the unitary Bose gas. Science 2017, 355, 377–380. [Google Scholar] [CrossRef] [PubMed]
  14. Makotyn, P.; Klauss, C.E.; Goldberger, D.L.; Cornell, E.; Jin, D.S. Universal dynamics of a degenerate unitary Bose gas. Nat. Phys. 2014, 10, 116–119. [Google Scholar] [CrossRef]
  15. Smith, D.H.; Braaten, E.; Kang, D.; Platter, L. Two-Body and Three-Body Contacts for Identical Bosons near Unitarity. Phys. Rev. Lett. 2014, 112, 110402. [Google Scholar] [CrossRef] [PubMed]
  16. Eigen, C.; Glidden, J.A.P.; Lopes, R.; Cornell, E.A.; Smith, R.P.; Hadzibabic, Z. Universal prethermal dynamics of Bose gases quenched to unitarity. Nature 2018, 556, 221–224. [Google Scholar] [CrossRef] [PubMed]
  17. Colussi, V.E.; Corson, J.P.; D’Incao, J.P. Dynamics of Three-Body Correlations in Quenched Unitary Bose Gases. Phys. Rev. Lett. 2018, 120, 100401. [Google Scholar] [CrossRef] [PubMed]
  18. Blume, D.; Sze, M.W.C.; Bohn, J.L. Harmonically trapped four-boson system. Phys. Rev. A 2018, 97, 033621. [Google Scholar] [CrossRef]
  19. Werner, F.; Castin, Y. Unitary Quantum Three-Body Problem in a Harmonic Trap. Phys. Rev. Lett. 2006, 97, 150401. [Google Scholar] [CrossRef] [PubMed]
  20. Danilov, G.S. On the three-body problem with short-range forces. Sov. Phys. JETP 1961, 13, 349. [Google Scholar]
  21. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions: With Formulas, Graphs, and Mathematical Tables; Courier Corporation: Chelmsford, MA, USA, 1964; Volume 55. [Google Scholar]
  22. Jonsell, S.; Heiselberg, H.; Pethick, C.J. Universal Behavior of the Energy of Trapped Few-Boson Systems with Large Scattering Length. Phys. Rev. Lett. 2002, 89, 250401. [Google Scholar] [CrossRef] [PubMed]
  23. Wang, J.; D’Incao, J.P.; Esry, B.D.; Greene, C.H. Origin of the Three-Body Parameter Universality in Efimov Physics. Phys. Rev. Lett. 2012, 108, 263001. [Google Scholar] [CrossRef] [PubMed]
  24. 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]
  25. Werner, F. Atomes Froids Piégés en Interaction Résonnante: Gaz Unitaire et Problème à Trois Corps. Ph.D. Thesis, Université Pierre et Marie Curie, Paris, France, 2008. [Google Scholar]
  26. Efimov, V. Effective interaction of three resonantly interacting particles and the force range. Phys. Rev. C 1993, 47, 1876–1884. [Google Scholar] [CrossRef]
  27. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series, and Products; Academic Press: Cambridge, MA, USA, 2007; p. 822. [Google Scholar]
  28. Cohen-Tannoudji, C.; Guéry-Odelin, D. Advances in Atomic Physics: An Overview; World Scientific: Singapore, 2011. [Google Scholar]
  29. D’Incao, J.P.; Wang, J.; Colussi, V.E. Efimov Physics in Quenched Unitary Bose Gases. Phys. Rev. Lett. 2018, 121, 023401. [Google Scholar] [CrossRef] [PubMed]
  30. Naidon, P.; Endo, S. Efimov physics: A review. Rep. Prog. Phys. 2017, 80, 056001. [Google Scholar] [CrossRef] [PubMed]
  31. Olver, F.W. NIST Handbook of Mathematical Functions Hardback and CD-ROM; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  32. Castin, Y.; Werner, F. Single-particle momentum distribution of an Efimov trimer. Phys. Rev. A 2011, 83, 063614. [Google Scholar] [CrossRef]
  33. Portegies, J.; Kokkelmans, S. Efimov Trimers in a Harmonic Potential. Few-Body Syst. 2011, 51, 219. [Google Scholar] [CrossRef]
  34. Colussi, V.E.; van Zwol, B.E.; D’Incao, J.P.; Kokkelmans, S.J.J.M.F. Bunching, clustering, and the buildup of few-body correlations in a quenched unitary Bose gas. arXiv, 2019; arXiv:1901.10355. [Google Scholar]
Figure 1. Behavior of a ho 2 C 3 for both negative and positive eigenenergies in the Efimov channel. (a) Results for four values of a ho κ * chosen evenly spaced over a single log-period as indicated in the key. The j = 0 indicates the ground-state trimer in the free-space limit in the sense of Equation (12) and are contained in the grey shaded region. The black dashed-dotted line is the asymptotic result for a ho 2 C 3 in the limit E 3 b / ω given by Equation (33). The solid black lines are the asymptotic results for a ho 2 C 3 in the limit E 3 b / ω given by Equation (36) that oscillate around a ho 2 C 3 s 0 / π indicated by the black dashed line; (b) the transition of a ho 2 C 3 given by Equation (30) as a function of E 3 b / ω from positive to negative energies. The black dashed-dotted line is the same as in (a).
Figure 1. Behavior of a ho 2 C 3 for both negative and positive eigenenergies in the Efimov channel. (a) Results for four values of a ho κ * chosen evenly spaced over a single log-period as indicated in the key. The j = 0 indicates the ground-state trimer in the free-space limit in the sense of Equation (12) and are contained in the grey shaded region. The black dashed-dotted line is the asymptotic result for a ho 2 C 3 in the limit E 3 b / ω given by Equation (33). The solid black lines are the asymptotic results for a ho 2 C 3 in the limit E 3 b / ω given by Equation (36) that oscillate around a ho 2 C 3 s 0 / π indicated by the black dashed line; (b) the transition of a ho 2 C 3 given by Equation (30) as a function of E 3 b / ω from positive to negative energies. The black dashed-dotted line is the same as in (a).
Atoms 07 00019 g001
Figure 2. Behavior of a ho C 2 for both negative and positive eigenenergies in the universal and Efimov channels. (a) The behavior of a ho C 2 ( s , 0 ) for s 2 > 0 compared against the s + asymptotic result (black solid line) Equation (54); (b) results for a ho C 2 ( s , j ) for s 2 > 0 over a range of energies in the three universal channels s 1 4.465 , s 10 22.932 , and s 20 42.967 ; (c) results for a ho C 2 ( s 0 , j ) for four values of a ho κ * chosen evenly spaced over a single log-period (same as in Figure 1) as indicated in the key. The black dashed-dotted line is the asymptotic result for a ho C 2 in the limit E 3 b / ω given by Equation (51). The grey shaded region indicates the j = 0 eigenstates; (d) the transition of a ho C 2 s 0 , j as a function of E 3 b / ω from positive to negative energies. The black dashed-dotted line is the same as in (c).
Figure 2. Behavior of a ho C 2 for both negative and positive eigenenergies in the universal and Efimov channels. (a) The behavior of a ho C 2 ( s , 0 ) for s 2 > 0 compared against the s + asymptotic result (black solid line) Equation (54); (b) results for a ho C 2 ( s , j ) for s 2 > 0 over a range of energies in the three universal channels s 1 4.465 , s 10 22.932 , and s 20 42.967 ; (c) results for a ho C 2 ( s 0 , j ) for four values of a ho κ * chosen evenly spaced over a single log-period (same as in Figure 1) as indicated in the key. The black dashed-dotted line is the asymptotic result for a ho C 2 in the limit E 3 b / ω given by Equation (51). The grey shaded region indicates the j = 0 eigenstates; (d) the transition of a ho C 2 s 0 , j as a function of E 3 b / ω from positive to negative energies. The black dashed-dotted line is the same as in (c).
Atoms 07 00019 g002
Table 1. Comparison of the results for E 3 b / ω , a ho C 2 , and a ho 2 C 3 in this work with the results of Ref. [18] for the energetically lowest eight eigenstates. To obtain these results, we have chosen an oscillator length a ho κ * 0.43 .
Table 1. Comparison of the results for E 3 b / ω , a ho C 2 , and a ho 2 C 3 in this work with the results of Ref. [18] for the energetically lowest eight eigenstates. To obtain these results, we have chosen an oscillator length a ho κ * 0.43 .
sj E 3 b / E ho [18] E 3 b / E ho a ho C 2 [18] a ho C 2 a ho 2 C 3 [18] a ho 2 C 3
i s 0 00.50000.500035.635.70.4780.469
i s 0 12.92172.901922.422.50.3340.327
i s 0 25.08955.058718.819.00.3230.313
s 1 05.51875.465321.621.6 9 × 10 6 0
i s 0 37.20057.159717.317.00.3260.310
s 1 17.52587.465320.520.6 3 × 10 5 0
s 2 07.84527.81847.157.17 2 × 10 7 0
i s 0 49.28749.235115.415.70.3360.311

Share and Cite

MDPI and ACS Style

Colussi, V.E. Two-Body and Three-Body Contacts for Three Bosons in the Unitary Regime: Analytic Expressions and Limiting Forms. Atoms 2019, 7, 19. https://doi.org/10.3390/atoms7010019

AMA Style

Colussi VE. Two-Body and Three-Body Contacts for Three Bosons in the Unitary Regime: Analytic Expressions and Limiting Forms. Atoms. 2019; 7(1):19. https://doi.org/10.3390/atoms7010019

Chicago/Turabian Style

Colussi, Victor E. 2019. "Two-Body and Three-Body Contacts for Three Bosons in the Unitary Regime: Analytic Expressions and Limiting Forms" Atoms 7, no. 1: 19. https://doi.org/10.3390/atoms7010019

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