Abstract

We consider finite-range effects when the scattering length goes to zero near a magnetically controlled Feshbach resonance. The traditional effective-range expansion is badly behaved at this point, and we therefore introduce an effective potential that reproduces the full T-matrix. To lowest order the effective potential goes as momentum squared times a factor that is well defined as the scattering length goes to zero. The potential turns out to be proportional to the background scattering length squared times the background effective range for the resonance. We proceed to estimate the applicability and relative importance of this potential for Bose-Einstein condensates and for two-component Fermi gases where the attractive nature of the effective potential can lead to collapse above a critical particle number or induce instability toward pairing and superfluidity. For broad Feshbach resonances the higher order effect is completely negligible. However, for narrow resonances in tightly confined samples signatures might be experimentally accessible. This could be relevant for suboptical wavelength microstructured traps at the interface of cold atoms and solid-state surfaces.

1. Introduction

Cold atomic gases have enjoyed many great successes since the first realizations of Bose-Einstein condensates in the mid nineties [1]. Ensembles of ultracold atomic gases can be manipulated in magnetic or optical trap geometries and in lattice setups, effectively mimicking the structure of real materials and teaching us about their properties. In particular, extreme control can be exercised over the atom-atom interactions through the use of Feshbach resonance [2]. Tuning the system into the regime of resonant two-body interactions provides a controlled way of studying strongly correlated dynamics which is believed to be crucial for material properties such as high-temperature superconductivity or giant magnetoresistance.

Recently there has been extended interest in weakly interacting Bose-Einstein condensates for use as an atomic interferometer [3] and also to probe magnetic dipolar interactions in condensates [4]. This work was based on 39K atoms where a broad Feshbach resonance exists at a magnetic field strength of 𝐡0=402.4G [5] which allows a large tunability of the atomic interaction in experiments [6]. Similar tunability has also been reported in a condensate of 7Li [7]. The atomic interaction can be reduced by tuning the scattering length, π‘Ž, to zero, also known as zero crossing. In a Gross-Pitaevskii mean-field picture we can thus neglect the usual nonlinear term proportional to π‘Ž. The question is then what other interactions are relevant. As shown in [4], the magnetic dipole will contribute here.

In the Gross-Pitaevskii picture we might also ask whether higher order terms in the interaction can contribute around zero crossing. Recently it was shown that effective-range corrections can in fact influence the stability of condensates around zero crossing [8–10]. The Feshbach resonances used thus far in experiments have typically been very broad, and as a result the effective range, π‘Ÿπ‘’, will be small, rendering the higher order terms negligible. However, around narrow resonances this is not necessarily the case and finite-range corrections are not necessarily negligible.

For the two-component Fermi gas, there has been increased interest in producing a cold atom analog of the celebrated Stoner model of ferromagnetism [11] which applies to repulsively interacting fermions. Theoretical proposals indicate that this should be possible [12–19], and an MIT experiment subsequently announced indications of the ferromagnetic transition [20]. The results caused controversy since the spin domains were not resolved [21–24]. A later experiment in the same group did not find evidence of the ferromagnetic transition [25]. However, these studies consider broad Feshbach resonances, and the situation with narrow resonances is less clear. One can imagine that finite-range corrections could play a role in driving the phase transition. In fact, a recent experiment in Innsbruck [26] has found increased lifetimes of the repulsive gas in the strongly imbalanced case, providing hope that decay into molecules can be controlled and ferromagnetism can be studied.

The systematic inclusion of finite-range effects through derivative terms in zero-range models was begun in the study of nuclear matter decades ago [27, 28]. Later on the intricacies of the cut-off problems that arise in this respect were considered by many authors both for the relativistic and nonrelativistic case (see [29] for discussion and references). In the context of cold atoms and Feshbach resonances, we need to use a two-channel model [30] in order to take the lowest order finite-range term into account. Similar models were already introduced in [31] and denoted resonance models (see f.x. [32] for a comprehensive review of scattering models for ultracold atoms). We note that whereas resonance models treat the closed-channel molecular state as a point boson the model of [30] treats the molecule more naturally as a composite object of two atoms. In the end the parameters of the two models turn out to be similarly related to the physical parameters of Feshbach resonances (see for instance the discussion of resonance models in [32]).

In Figure 1 we show calculations of scattering length and effective range for the Feshbach resonance at 𝐡=202.1 G in 40K in both a coupled-channel model [33] and in the zero-range model discussed here. We see the effective range being roughly constant at resonance and then start to diverge at zero crossing. The zero-range model provides a good approximation to the full calculations and for many-body purposes it is preferable due to its simplicity.

Whereas the earlier work of [31] considered the regime close to the resonance, we will be exclusively concerned with zero crossing. To our knowledge the intricacies of this region have not been addressed in the literature in the context of Feshbach resonances. Around zero crossing the Feshbach model turns out to have a badly behaved effective-range expansion. The parameters obtained from the effective-range expansion should therefore be used with extreme caution as the series is divergent at this point. However, as we show in this paper, the finite-range corrections obtained from the full 𝑇-matrix at low momenta via an effective potential turn out to be the same as one would naively expect based on the effective-range expansion. After introducing the effective potential we consider its applicability and importance in the case of Bose-Einstein condensates and for two-component Fermi gases where the attractive nature of the effective interaction at zero crossing could lead to collapse above a certain critical particle number or to pairing instability and superfluidity. In general, we find that tight external confinement is a necessary condition for the higher order effects to dominate the magnetic dipole interaction and be experimentally observable.

2. Two-Channel Model

We consider a two-channel 𝑠-wave Feshbach model with zero-range interactions [30] for which the on-shell open-open channel 𝑇-matrix as a function of magnetic field, 𝐡, is π‘‡π‘œπ‘œ(𝐡)=4πœ‹β„2ξ€Έπ‘Ž/π‘šπ‘π‘”ξ€·ξ€·β„1+Ξ”πœ‡Ξ”π΅/2π‘ž2ξ€·/π‘šβˆ’Ξ”πœ‡π΅βˆ’π΅0ξ€Έξ€Έξ€Έβˆ’1+π‘–π‘Žπ‘π‘”π‘ž,(1) where Ξ”πœ‡ is the difference between the magnetic moments in the open and closed channel, π‘ž is the relative momentum of the atoms of mass π‘š, π‘Žπ‘π‘” is the scattering length away from the resonance at magnetic field 𝐡0, and Δ𝐡 is the width of the resonance. We can compare this to the standard vacuum expression for the 𝑇-matrix in terms of the phase-shift given by π‘‡πœˆ=ξ€·4πœ‹β„2ξ€Έπ‘Ž/π‘šβˆ’π‘žπ‘Žcot.𝛿(π‘ž)+π‘–π‘Žπ‘ž(2)

Typically, one has the low-energy expression βˆ’π‘žcot𝛿(π‘ž)β†’βˆ’1/π‘Ž which implies that π‘‡πœˆβŸΆξ€·4πœ‹β„2ξ€Έπ‘Ž/π‘š1+π‘–π‘Žπ‘žβŸΆ0.(3) However, as we now discuss, for the realistic two-channel 𝑇-matrix for Feshbach resonances, the quantity βˆ’π‘žcot𝛿(π‘ž) is not well defined, and the conclusion that the 𝑇-matrix vanishes at zero crossing is only true for zero momentum, π‘ž=0, as we now discuss.

From (1) and (2) we obtain the relation for the phase-shift: π‘žcot𝛿(π‘ž)=βˆ’1π‘Žπ‘π‘”ξƒ©1+Ξ”πœ‡Ξ”π΅ξ€·β„2π‘ž2ξ€·/π‘šβˆ’Ξ”πœ‡π΅βˆ’π΅0ξƒͺξ€Έξ€Έβˆ’1.(4) We now expand the right-hand side in powers of π‘ž as is usually done in an effective-range expansion. This yields π‘žcot𝛿(π‘ž)=βˆ’1π‘Ž+(𝐡)βˆžξ“π‘›=1βˆ’1π‘Žπ‘π‘”ξ‚Έβˆ’π‘Žπ‘π‘”π‘Ÿπ‘’02ξ‚Ήπ‘›ξ‚Έπ‘Žπ‘π‘”ξ‚Ήπ‘Ž(𝐡)βˆ’1𝑛+1π‘ž2𝑛,(5) where π‘Ž(𝐡)=π‘Žπ‘π‘”(1βˆ’Ξ”π΅/(π΅βˆ’π΅0)) is the common parametrization from single-channel models and π‘Ÿπ‘’0=βˆ’2ℏ2/(π‘šΞ”π΅Ξ”πœ‡π‘Žπ‘π‘”) is the background value of the effective range around the resonance. From (5) we can now read off all coefficients in an effective-range expansion with their full 𝐡-field dependence. For instance, the effective range is given simply by π‘Ÿπ‘’=π‘Ÿπ‘’0[(π‘Žπ‘π‘”/π‘Ž)βˆ’1]2, which is divergent when π‘Ž(𝐡)β†’0. We also clearly see that all the other coefficients are divergent in that limit. This is signaled also before doing the full expansion in π‘ž as the first term in (5) diverges at zero crossing. However, in effective potentials derived from the 𝑇-matrix these problems are not transparent as the lowest order coefficient is proportional to π‘Ž(𝐡) (see (12)). Below we will discuss what kind of constraints this introduces on the applicability of the effective-range expansion near zero crossing. We note that similar issues were briefly discussed in a different context in [34] where an equivalent to (7) below was obtained.

Let us first consider the low-π‘ž limit and compare the full 𝑇-matrix with the effective-range expansion as zero crossing is approached. Taking the low-π‘ž limit of (4) at zero crossing where Δ𝐡/(π΅βˆ’π΅0)=1, we find π‘žcot𝛿(π‘ž)βŸΆβˆ’1π‘Žπ‘π‘”βˆ’Ξ”πœ‡Ξ”π΅β„2π‘ž2,/π‘š(6) which diverges as π‘žβˆ’2. Therefore the coefficients of the expansion in (5) must necessarily diverge in order to retain any hope of describing the low-π‘ž behavior. Furthermore, since the expansion is an alternating series and therefore slowly converged, we also conclude that many terms must be retained for a fair approximation at very small but nonzero π‘ž. The same conclusion can be reached by considering the radius of convergence of (5), which we find by locating the pole in (4) at ℏ2π‘ž2/π‘š=Ξ”πœ‡(π΅βˆ’π΅0βˆ’Ξ”π΅). This radius indeed goes to zero at zero crossing. We are thus forced to conclude that the effective-range expansion breaks down near zero crossing.

2.1. Effective Potential at Zero Crossing

Since the effective-range expansion is insufficient we consider the full 𝑇-matrix in the low-π‘ž limit at zero crossing. To lowest order we have π‘‡π‘œπ‘œξ€·π΅=𝐡0ξ€Έ+Δ𝐡=βˆ’4πœ‹β„2π‘Žπ‘π‘”π‘šβ„2π‘ž2ξ€·π‘žπ‘šΞ”πœ‡Ξ”π΅+𝑂4ξ€Έ.(7) Using the expression for π‘Ÿπ‘’0, this can be written 4πœ‹β„2π‘šπ‘Ž2π‘π‘”π‘Ÿπ‘’02π‘ž2.(8) Knowing the 𝑇-matrix at low π‘ž we can now proceed to find an effective low-π‘ž potential through the Lippmann-Schwinger equation: 𝑉=π‘‡βˆ’π‘‡πΊ0𝑉,(9) where 𝐺0=(πΈβˆ’π»0+𝑖𝛿)βˆ’1 is the free space Green’s function [35]. This equation can be solved for 𝑇(π‘ž,π‘žβ€²)βˆπ‘ž2+π‘žβ€²2 (the symmetrized version of the full 𝑇-matrix) in an explicit cut-off approach [29, 35] and then be expanded to order π‘ž2 for consistence with the input 𝑇-matrix. In the long-wavelength limit we can take the cut-off to zero [35] and for the on-shell effective potential we then obtain the obvious answer: 𝑉(π‘ž)=4πœ‹β„2π‘šπ‘Ž2π‘π‘”π‘Ÿπ‘’02π‘ž2(10) in momentum space. The effective potential in real-space is now easily found by canonical substitution (πͺβ†’βˆ’π‘–βˆ‡) and appropriate symmetrization [36]. We have 𝑉(𝐫)=βˆ’4πœ‹β„2π‘šπ‘Ž2π‘π‘”π‘Ÿπ‘’0212ξƒ¬β†βˆ‡2𝐫𝛿(𝐫)+𝛿(𝐫)β†’βˆ‡2𝐫.(11) Notice that the Lippmann-Schwinger approach is nonperturbative as opposed to the perturbative energy shift method [36, 37].

2.2. Comparison to Effective-Range Expansion and Energy-Shift Method

Away from zero crossing one can easily relate the effective-range expansion to an effective potential through the perturbative energy shift method [18, 25, 26]. To second order the 𝑠-wave effective potential is 𝑉(𝐫)=4πœ‹β„2π‘Žπ‘šξƒ¬π‘”π›Ώ(𝐫)+22ξƒ©β†βˆ‡2𝐫𝛿(𝐫)+𝛿(𝐫)β†’βˆ‡2𝐫,ξƒͺξƒ­(12) where the first term is the effective interaction usually employed in mean-field theories of cold atoms [35]. In terms of π‘Ž and π‘Ÿπ‘’, we have 𝑔2=π‘Ž2/3βˆ’π‘Žπ‘Ÿπ‘’/2 [36, 37] with the field-dependent π‘Ž=π‘Ž(𝐡) and π‘Ÿπ‘’=π‘Ÿπ‘’(𝐡).

At zero-crossing the first term in (12) vanishes and one might expect the second term to vanish as well. However, in the naive effective-range expansion of the two-channel model discussed above we saw that π‘Ÿπ‘’ diverges as π‘Žβˆ’2 and we therefore have limπ‘Žβ†’0π‘Žπ‘”2π‘Ž=βˆ’2π‘π‘”π‘Ÿπ‘’02.(13) In particular, if we for a moment ignore π‘ž4 terms in the effective-range expansion, we recover exactly the same effective potential as in (11) at zero crossing. The finite limiting result in (13) shows that the potential in (12) is well defined as π‘Žβ†’0, provided that appropriate regularization and renormalization are performed. Equation (12) thus applies equally well at resonance (π‘Žβ†’βˆž) where the gradient terms are small and at zero crossing where the lowest order delta function term is unimportant. It is thus a well-defined effective potential over the entire range of a Feshbach resonance.

We therefore see that even though the effective-range expansion has divergent coefficients at zero crossing, the lowest order does in fact give the same effective potential as the full 𝑇-matrix if we apply it naively. The effective-range expansion should thus be viewed as an asymptotic series. However, we cannot use the effective-range expansion to estimate the validity of the second-order effective potential since the radius of convergence goes to zero at zero crossing as discussed above.

The two-channel model in (1) compares well with a coupled-channel calculation [33] as shown in Figure 1. It also compares well to other scattering models [38, 39] that include finite-range effects. In fact, the model used here compares well with the analytical models of [38] when π‘Ž(𝐡) and π‘Ÿπ‘’(𝐡) have the field-dependence introduced above. This can be seen for instance in Figure 12 of [38], although a difference is that our π‘Ž(𝐡) and π‘Ÿπ‘’(𝐡) are parametrization and not taken from coupled-channels values as in [38] (our Figure 1 quantifies the difference which is largest on π‘Ÿπ‘’(𝐡)). However, here we are concerned with the behavior when π‘Ž(𝐡)β†’0 in the context of Feshbach resonances which is not addressed in [38, 39]. We note that the resonance models of [31] and the two-channel and resonance models in [32] are very similar to the model employed here, but again those references do not consider the specific problems arising when π‘Ž(𝐡)β†’0. In addition and in contrast to previous discussions, here we construct appropriate zero-range pseudo potentials that work around zero crossing.

3. Relation to Experiments

Above we only retained terms of order π‘ž2 in the full 𝑇-matrix. We now estimate the energy regime in which this expression is valid. Demanding that the π‘ž4 term be smaller than the π‘ž2 term gives the criterion: ℏ2π‘ž2π‘šβ‰ͺℏ2π‘š||π‘Žπ‘π‘”π‘Ÿπ‘’0||.(14) We relate this condition to recent experiments with bosonic condensates of 39K working around zero crossing [3]. The resonance used there is very broad (Δ𝐡=βˆ’52G) with π‘Žπ‘π‘”=βˆ’29π‘Ž0 and π‘Ÿπ‘’0=βˆ’58π‘Ž0 (π‘Ž0 is the Bohr radius). The right-hand side of (14) is 2.3β‹…10βˆ’7 eV, corresponding to a temperature of about 3 mK. Since the experiments are performed at much lower temperatures the approximation above is certainly valid. However, as π‘Žπ‘π‘” and particularly π‘Ÿπ‘’0 are small, the front factor in (11) is also small. The relevant scale of comparison is the outer trap parameter 𝑏 [9] which is typically of order 1 μm, yielding a vanishing ratio |π‘Ž2π‘π‘”π‘Ÿπ‘’0|/𝑏3∼10βˆ’9. For broad Feshbach resonances the higher order interactions can thus be safely ignored. For very narrow resonances the situation potentially changes as π‘Ÿπ‘’0 can be very large and make the potential in (11) important. As an example, we consider the narrow resonance in 39K at 𝐡0=25.85G with Δ𝐡=0.47G, π‘Žπ‘π‘”=βˆ’33π‘Ž0, and π‘Ÿπ‘’0=βˆ’5687π‘Ž0 [5]. The right-hand side of (14) is now 2β‹…10βˆ’9 eV, corresponding to 24 μK. This is again much higher than experimental temperatures. A more careful argument can be made from the energy per particle of the noncondensed cloud. Ignoring the trap, we have 𝐸/𝑁=0.770π‘˜π΅π‘‡π‘(𝑇/𝑇𝑐)5/2 (𝑇𝑐 is the critical temperature) [35]. For a sample of 3β‹…104 a critical temperature of 100 nK was reported in [6]. Using this 𝑇𝑐 we find that 𝑇β‰ͺ900 nK for (14) to hold. Again this is within the experimental regime. The effective potential approach should therefore be applicable around zero crossing for narrow resonances. However, even with this narrow resonance we find |π‘Ž2π‘π‘”π‘Ÿπ‘’0|/𝑏3∼10βˆ’7, and the effect is still completely negligible.

In order to increase the relevance of the higher order term, we now consider some very narrow resonances that have been found in 87Rb. In particular, the resonance at 𝐡0=9.13G [40] which was recently utilized in nonlinear atom interferometry [41]. We have Δ𝐡=0.015G, π‘Žπ‘π‘”=99.8π‘Ž0, and Ξ”πœ‡=2.00πœ‡π΅ [42], which gives π‘Ÿπ‘’0=βˆ’19.8β‹…103π‘Ž0 and a ratio |π‘Ž2π‘π‘”π‘Ÿπ‘’0|/𝑏3=2.92β‹…10βˆ’5(1πœ‡m/𝑏)3. A trap length of π‘βˆΌ0.5 μm as used in [41] would thus yield 10βˆ’4 and demonstrates that higher order corrections can safely be neglected. For a ratio of 1 we need π‘βˆΌ0.03 μm which is unrealistically small in current traps or optical lattices. However, a resonance of width Δ𝐡=0.0004G is known in the same system at 𝐡0=406.2G [43] with π‘Žπ‘π‘”=100π‘Ž0 and Ξ”πœ‡=2.01πœ‡π΅ [42]. In this case we find π‘Ÿπ‘’0=βˆ’7.4β‹…105π‘Ž0 and a much more favorable ratio of |π‘Ž2π‘π‘”π‘Ÿπ‘’0|/𝑏3=0.001(1πœ‡m/𝑏)3. Here we see that a ratio of 1 is achieved already for π‘βˆΌ0.1πœ‡m which not far off from tight traps or optical lattice dimensions. In terms of temperature we still have to be in the ultralow regime of 𝑇≲30nK according to (14) for the latter resonance.

Consider now a fermionic two-component system where 𝑠-wave interactions are dominant. Since we have π‘Ÿπ‘’0<0 for all Feshbach resonances [42], the effective potential in (10) is attractive, and the system could potentially be unstable toward a paired state or become unstable to collapse above a critical particle number. For simplicity we will use the semiclassical Thomas-Fermi approach to describe a gas with equal population of the two components and estimate the critical particle number. Assuming an isotropic trapping potential with length scale βˆšπ‘=ℏ/π‘šπœ” where πœ” is the trap frequency, the ground-state density, 𝜌(𝐱), can be found by minimization and satisfies ξ‚Έπœ‡βˆ’1β„πœ”2𝐱𝑏2ξ‚Ή=12ξ€·π‘˜πΉξ€Έ(𝐱)𝑏2βˆ’4π›Όξ€·π‘˜30πœ‹πΉξ€Έ(𝐱)𝑏5,(15) where 𝜌(𝐱)=π‘˜πΉ(𝐱)/6πœ‹2 and 𝛼=π‘Ž2𝑏𝑔|π‘Ÿπ‘’0|/𝑏3. The maximum allowed momentum and chemical potential, πœ‡, is found by solving for the turning point of the right-hand side of (15) which gives π‘˜max𝑏=3πœ‹ξ‚„2𝛼1/3,πœ‡max=3ξ€·π‘˜10β„πœ”max𝑏2.(16) We can now compare this π‘˜max to the value obtained from the noninteracting density within the Thomas-Fermi approximation at the center of the trap. In terms of the number of particles in each component, 𝑁, at the center of the trap we have π‘˜πΉ(0)π‘β‰ˆ1.906𝑁1/6 [35]. By equating these two expression we obtain an estimate for the critical number of particles, 𝑁max. Inserting the relevant units, we have 𝑁max=2β‹…1025ξ‚΅π‘Ž0π‘Žπ‘π‘”ξ‚Ά4ξ‚΅π‘Ž0π‘Ÿπ‘’0ξ‚Ά2𝑏1πœ‡mξ‚Ά6,(17) where π‘Ž0 is the Bohr radius. We note that the scaling 𝑁maxβˆπ›Όβˆ’2 can also be obtained by considering the point at which the monopole mode becomes unstable.

Typical numbers for common fermionic species 6Li or 40K in the lowest hyperfine states [42] lead to 𝑁max∼1012 for 𝑏=1 μm. This is of course a huge number and experiments are well within this limit. Even if one reduced the trap length by a factor of ten and made the presumably unrealistic assumption that the particle number remains the same we still have 𝑁β‰ͺ𝑁max. The reason is that the 𝑠-wave Feshbach resonances utilized in the two-component gases are generally broad in order to study the universal regime. If we consider the narrow resonance at 𝐡0=543.25G in 6Li [44] with Δ𝐡=0.1G, π‘Žπ‘π‘”=60π‘Ž0, and Ξ”πœ‡=2.00πœ‡π΅ [42], we have 𝑁max∼2β‹…1013(𝑏/1πœ‡m)6. This is somewhat better but we still need π‘βˆΌ0.06 μm to get to an experimentally relevant 𝑁max∼106. We have to conclude that higher order 𝑠-wave interactions are highly unlikely to be observable through monopole instabilities. In light of this it seems better to consider 𝑝-wave resonances which are much more narrow in general. However, also here extremely small trap sizes appear necessary [45].

The instability toward Cooper pairing around zero crossing can also be estimated in simple terms. In general the critical temperature is π‘‡π‘βˆΌπ‘‡πΉexp(βˆ’1/𝑁0|π‘ˆ|), where 𝑁0=π‘šπ‘˜πΉ(0)/2πœ‹2ℏ2 is the density of states at the Fermi energy in the trap center and π‘ˆ<0 is a measure of the attraction. For the latter we use the effective potential in momentum space from (10) and make the assumption that π‘žβˆΌπ‘˜πΉ(0). Using the expression for π‘˜πΉ(0) in terms of 𝑁 above, we find 1𝑁0||π‘ˆ||=1.5β‹…1012βˆšπ‘ξ‚΅π‘1πœ‡mξ‚Ά3ξ‚΅π‘Ž0π‘Žπ‘π‘”ξ‚Ά2π‘Ž0||π‘Ÿπ‘’0||.(18) For broad resonances in 6Li or 40K this exponent is of order 103 and 𝑇𝑐 is thus vanishingly small. However, the scaling with trap size can help and if we imagine reducing to 𝑏=0.1 μm, we find 𝑇𝑐≲0.5𝑇𝐹 for 𝑁=106 atoms. For the narrow resonance in 6Li discussed above, we find that π‘‡π‘βˆΌ0.5𝑇𝐹 with 𝑁=106 can be achieved for π‘βˆΌ0.5 μm and π‘‡π‘βˆΌ0.1𝑇𝐹 for 𝑁=105. Thus there may be a possibility to reach the pairing instability near zero crossing if high particle numbers can be cooled in tight traps and narrow resonances are used.

While the suboptical wavelength trapping sizes needed for the above effects to be large are not achievable with typical optical or magnetic traps or optical lattice setups, they could potentially be reached via hybrid setups where atoms are trapped near a surface. Inspired by surface plasmon subwavelength optics [46], nanoscale trapping for neutral atoms has been studied [47, 48], and micropotential traps with width less than 100 nanometer (<0.1 μm) are within reach [49]. In these very tightly confined systems, it is very likely that finite-range effects could be enhanced. Devices that provide an interface between atoms and solid-state systems are under intense study at the moment, and our considerations here imply that finite-range corrections should be considered when the scattering length is tuned close to zero.

3.1. Dipole-Dipole Interactions

The discussion above ignores the dipole-dipole interaction discussed in the introduction which will compete against the higher order effective potential from the Feshbach resonance. A simple estimate can be made along the lines of the discussion in [35]. The external trapping potential is the characteristic scale of spatial variations and we thus find a ratio, π‘Ÿ, of magnetic dipole-dipole, π‘ˆmd, to higher-order 𝑠-wave zero-range interaction strength, π‘ˆ2, which can be written as π‘ˆπ‘Ÿ=mdπ‘ˆ2=π‘Ž0𝑏2π‘Ž2𝑏𝑔||π‘Ÿπ‘’0||𝑏=35.71πœ‡mξ‚Ή2ξ‚Έ100π‘Ž0π‘Žπ‘π‘”ξ‚Ή21000π‘Ž0||π‘Ÿπ‘’0||.(19) For π‘Ÿ<1 the higher order interaction term will therefore dominate the magnetic dipole term. For the case of narrow resonances in 87Rb discussed above we find π‘ŸβˆΌ0.11(𝑏/1πœ‡m)2 for the resonance at 𝐡0=9.13G and π‘ŸβˆΌ0.05(𝑏/1πœ‡m)2 for the one at 𝐡0=406.2G. For the narrow resonance in 6Li at 𝐡0=543.25G we find π‘ŸβˆΌ1.4(𝑏/1πœ‡m)2. These ratios clearly indicate that magnetic dipole-dipole interactions can be suppressed relative to higher order zero-range terms for narrow Feshbach resonances and standard trap sizes. This dominance becomes even stronger for the tight traps needed for the realization of the effects discussed above, and we thus conclude that interference of the magnetic dipole-dipole term is not a major concern.

4. Conclusions and Outlook

In this paper we have discussed the effective potential around a Feshbach resonances as the scattering length is tuned to zero and finite-range corrections become important. We showed that the effective-range expansion is badly behaved and the effective potential most be defined from the 𝑇-matrix. We have demonstrated that the low momenta effective potential obtained from the full 𝑇-matrix agrees with one obtained naively from the effective-range expansion when the scattering length goes to zero. Thus even though the effective-range expansion has divergent coefficients at zero-crossing the first terms of the associated effective potential yield consistent results. We then estimated the effects of the terms on different condensates. Since the effective potential at zero crossing is attractive it may induce various instabilities which we considered for the case of a two-component Fermi gas under harmonic confinement.

For the broad Feshbach resonances used in current experiments the effective potential discussed here are negligible and the dipole-dipole interaction dominates completely at zero crossing. However, for narrow resonances in very tightly confined systems some of the effects might be detectable. In particular, future generations of microtraps with suboptical wavelength trap sizes using surface plasmons could be small enough to make finite-range effects important. The competing dipole interaction is small for narrow resonances in tight confinement. However, it is conceivable that effects of spherically symmetric higher order terms could be separated from dipolar effects which change with system geometry [4].

Small trapped Fermi systems have recently become an experimental reality with particle numbers ranging from two to ten [50]. For two atomic fermions with different internal states, the system turns out to be well described by the analytic zero-range model of Busch et al. [51–56], and similarly for three fermions [57, 58]. Effective-range corrections to these results have also been studied [59–62]. Mesoscopic Fermi systems (less than about 50 particles) have been studied in harmonic traps using a number of numerical methods, [63–79] with particular emphasis on the unitary regime where the scattering length diverges. It would be interesting to investigate the situation also around zero-crossing of a narrow resonance where the effective range is sizable. A preliminary study along this line for three bosons is discussed in [80].

Another interesting direction of future work is the study of the contact introduced by Tan [81–90] to describe the universal behavior of strongly interacting quantum gases at a broad resonance where the range corrections are negligible, for instance through the tail of the momentum distribution which is predicted to behave as 𝐢/π‘˜4, where 𝐢 is the contact and π‘˜ the momentum of a single particle. The relations found by Tan [81, 82] have subsequently been confirmed experimentally in three dimensions [91–93]. While the contact originally pertains to two-body correlations, signatures of three-body physics in momentum distributions have also been studied both theoretically [86, 94–98] and experimentally [99]. While a few studies have considered the universal behavior when including the effective range term [100, 101], it would be very interesting to consider the regime around zero crossing for a narrow resonance where the background effective range parameter.

Acknowledgments

The author would like to thank Martin Thøgersen for very fruitful collaborations. Correspondence with Georg Bruun about two-channel models is highly appreciated. The author is grateful to Nicolai Nygaard for discussions and for producing Figure 1. The author acknowledges the hospitality of the Niels Bohr Institute, Blegdamsvej 17, 2100 Copenhagen Ø, Denmark. This work was supported by the Villum Kann Rasmussen foundation.