# Analytical model of overlapping Feshbach resonances

###### Abstract

Feshbach resonances in ultracold collisions often result from an interplay between many collision channels. Simple two-channel models can be introduced to capture the basic features, but cannot fully reproduce the situation when several resonances from different closed channels contribute to the scattering process. Using the formalism of multichannel quantum defect theory we develop an analytical model of overlapping Feshbach resonances. We find a general formula for the variation of the scattering length with magnetic field in the vicinity of an arbitrary number of resonances, characterized by simple parameters. Our formula is in excellent agreement with numerical coupled channels calculations for several cases of overlapping resonances in the collisions of two Li atoms or two Cs atoms.

###### pacs:

34.10.+x,34.50.Cx,03.65.Nk## I Introduction

Magnetically tunable Feshbach resonances have become an essential tool for the control of atomic interactions in ultracold quantum gases Chin et al. (2010). This is because such resonances allow the continuous tuning of the two-body -wave scattering length, a key parameter that controls phenomena such as Bose-Einstein condensation Inouye et al. (1998); Cornish et al. (2000), the crossover physics between Bose-Einstein condensate and Bardeen-Cooper-Schriefer paring in fermonic gases of mixed species Bourdel et al. (2004); Zwierlein et al. (2004), and the formation of weakly bound Feshbach molecules Jochim et al. (2003); Lang et al. (2008). Near the pole position of a resonance at magnetic field , the scattering length takes on the following simple form as the field is varied Moerdijk et al. (1995); Timmermans et al. (1999):

(1) |

where is a constant “background” scattering length in the region of the pole, and is the width of the resonance; at the “zero crossing” where . While scattering lengths are difficult to measure accurately, except for the locations of extreme features such as poles or zero crossings, they can generally be calculated accurately for alkali metal atom species using coupled channels calculations based on models of the ground state potentials Chin et al. (2010). An accurate function is essential to interpret experiments involving three-body physics, since a key component of the theory is to know the scattering length at the laboratory field Kraemer et al. (2006); Knoop et al. (2009); Zaccanti et al. (2009); Gross et al. (2010a); Ferlaino et al. (2011); Berninger et al. (2011, 2013); Dyke et al. (2013).

While Eq. (1) is applicable to an isolated single resonance, there are many cases of multiple overlapping resonances, and a generalization to such cases is needed. We derive here an analytic representation of due to a set of overlapping resonances from different closed channels interacting with a single open channel based on the Mies version of multichannel quantum defect theory (MQDT) Mies (1984); Mies and Julienne (1984); Julienne and Mies (1989); Mies and Raoult (2000); Julienne and Gao (2006); Julienne (2009a):

(2) |

where the sum in the denominator represents a shift due to the mutual interaction of the resonances with the open scattering channel. We will show the connection between this formula and an alternative product form that is based on a square well model of multiple resonances Lange et al. (2009).

This paper will give the MQDT framework for deriving Eq. (2) and illustrate it using fits to coupled channel calculations involving overlapping resonances relevant to the Efimov physics of colliding Li atoms or Cs atoms. We will discuss how mutual interaction between the resonances affects their positions and local widths. In particular, we will show how a narrow resonance on the shoulder of a broad one can be viewed as a “local” resonance in with a modified background and width due to the presence of the other resonance.

This work is organized as follows. In Sec. II we briefly review the formulation of MQDT for van der Waals interactions. Sec. III is devoted to derivation and analysis of the formulas describing a single Feshbach resonance using two-channel MQDT theory. In Sec. IV we extend the model to the case of an arbitrary number of closed channels and resulting overlapping resonances. Sec. V contains comparison of our theory to numerical and experimental results as well as characterization of some known resonances in lithium and cesium. Conclusions are drawn in Sec. VI.

## Ii Quantum defect theory.

The collision of two cold atoms, including internal spin degrees of freedom and the effect of short-range forces, can generally be described by a multichannel matrix Schrödinger equation

(3) |

where , , and are respectively the interatomic distance, energy, reduced mass, and the unit matrix. References Köhler et al. (2006); Chin et al. (2010) discuss the various kinds of basis sets that can be used to set up to describe the interaction matrix for the collision of two ground S state alkali metal atoms. The number of scattering channels depends on the number of internal spin states of each atom and the number of orbital angular momenta needed to represent the collision. The matrix asymptotically approaches a diagonal form

(4) |

Here is the orbital angular momentum quantum number, is the coefficient of the van der Waals potential, and is the threshold energy for the th channel. Depending on collision energy , some of the channels can be open (), and some are closed (). Here we assume a single -wave open channel ().

We start our analysis by a brief review of the MQDT formalism of Mies Mies (1984); Mies and Julienne (1984); Mies and Raoult (2000), which offers some different physical insights than alternative MQDT formalisms Burke et al. (1998); Gao (2008, 2011). Because of the long-range van der Waals potential, we will like Gao Gao (2008, 2011) use as the unit of length and as the unit of energy. In MQDT one replaces by a set of reference potentials which reproduce the asymptotic form of and parameterizes the solution of the new diagonal reference problem at short distances using functions with WKB-like normalization , . These functions are connected with the long-distance scattering solutions , (for open channels) or exponentially decaying solutions (for closed channels) using the MQDT functions , and Julienne and Mies (1989); Mies and Raoult (2000); Julienne (2009a); Idziaszek et al. (2011)

(5) |

where ensures unit normalization of the bound state function; see Ref. Julienne (2009a) for a discussion of the normalization of the MQDT functions. Eigenvalues of the reference closed channel occur where . The general solution of the coupled channels problem can be written as

(6) |

where and are diagonal matrices of the reference solutions, is the quantum defect matrix which contains information about the short-range couplings, and gives the amplitudes. The short-range processes which determine are assumed to be present only at length scales . This brings an important simplification to the problem, as the energy scales associated with the short range are much bigger than , which is of the order of milikelvins. As a result, for ultracold collisions the matrix can be regarded as energy-independent. Moreover, assuming that the potential varies from its long-range form only at short distances makes it possible to use the analytic theory of van der Waals interactions for the MQDT functions Gao (1998, 2000).

The observable properties of the system are given in terms of the open channel block of the scattering matrix . Within the framework of MQDT, it can be obtained from the quantum defect matrix and the open channel quantum defect functions, as follows Mies (1984):

(7) |

where is a diagonal matrix with elements giving the phase shifts for the reference potentials, and

(8) |

The matrices of the reference channel quantum defect functions , and are diagonal. The renormalized open channel matrix is

(9) |

where the indices and respectively denote open and closed channels. Here we assume a single open channel.

## Iii MQDT description of a two-channel resonance.

We will now provide the MQDT description of a single magnetically tunable Feshbach resonance in the simplest case when only one open and one closed channel are present Julienne and Gao (2006). This is general, since Mies et al. Mies et al. (2000) showed how a problem with a single isolated resonance due to multiple closed channels can be reduced to a problem with a single effective closed channel; see also Ref. Nygaard et al. (2006). By choosing the reference potentials to reproduce the true scattering lengths of the uncoupled open and closed channels, the quantum defect matrix will contain only off-diagonal terms

(10) |

The dimensionless short range coupling parameter is assumed to be independent of and . Substituting this into equations (7)-(9), we obtain the matrix (which in the case of a single open channel is just a complex number ) in a factored form with the two factors representing the background and resonant scattering parts:

(11) |

The open reference channel phase shift as , and for convenience we choose to write the background scattering length in our length units as , using the mean scattering length of the van der Waals potential, , introduced by Gribakin and Flambaum Gribakin and Flambaum (1993).

The and functions control the threshold behavior of resonance scattering, giving the respective amplitude and phase relations between the short- and long-range reference functions. Discussion of the physical meaning of these MQDT functions in the cold collision context is given in Refs. Julienne and Mies (1989); Mies and Raoult (2000); Julienne and Gao (2006); Julienne (2009b). When the collision energy in the open channel becomes large compared to , then and . The threshold behavior for -waves as is Mies and Raoult (2000)

(12) |

The function in Eq. (11) vanishes at an eigenvalue , where the quantum number labels the vibrational level. Near such an eigenvalue we can expand the energy-dependence as Mies and Julienne (1984)

(13) |

where is the field-dependent position of the “bare” closed channel eigenvalue, is the magnetic moment difference between the bare open and closed channel states (measured in per gauss units), and is the magnetic field at which the bare bound state crosses the open channel threshold. Note that is the mean vibrational spacing near , and is the corresponding vibrational frequency. If does not represent the last bound state of the closed channel, a useful approximation is .

It is straightforward to relate the expression in Eq. (11) to conventional resonant scattering by rewriting the resonant factor as

(14) |

where

(15) |

is a constant that represents a short-range decay width. Note that the form in Eq. (14) gives the entire -dependence in the term from the expansion of , while the and MQDT functions give the near-threshold variation with energy. The two terms in the denominator proportional to represent the energy-dependent shift and width due to the threshold resonance. The numerator of the pole term in Eq. (14) represents the threshold decay width

(16) |

where is the short-range coupling between the bare open and closed channels. Thus, represents the width when the open channel scattering wave function is replaced by the wave function with short range WKB normalization. Following Ref. Mies and Julienne (1984), when the quantity can be interpreted as the short-range probability that the bound state decays into the open channel during a single vibrational cycle. Multiplying by converts the short-range probability into the proper threshold probability. Consequently, the bare bound state decay rate can be interpreted as the decay probability per cycle times the closed channel vibrational frequency (number of cycles per second).

Simple algebraic transformations of (11) give the energy-dependent scattering length Bolda et al. (2002); Krych and Idziaszek (2009); Idziaszek and Julienne (2010), defined as , in the form

(17) |

By taking the limit and making use of relations (13)-(15), we obtain the standard scattering length

(18) |

We can now define a resonance width by writing the “pole strength” in the numerator in Eq. (18) as

(19) |

Note that the so defined is the same as in Eq. (22) of Chin et al. Chin et al. (2010); the dimensionless resonance strength parameter defined by Chin et al. Chin et al. (2010) is . The (or ) parameter is very important for characterizing the properties of an isolated resonance and also may turn out to be of particular importance in the analysis of three-body recombination near the resonance Wang and Julienne (2013).

With the above definitions, and using the threshold properties in Eq. (12), the -wave scattering length is given by the familiar formula Moerdijk et al. (1995); Timmermans et al. (1999); Chin et al. (2010)

(20) |

where the width and shift are proportional to :

(21) | |||||

(22) |

The latter formula for the shift, previously given in Refs. Köhler et al. (2006); Julienne and Gao (2006); Chin et al. (2010) without derivation, shows that the shift can not be much larger in magnitude than , towards which it tends for large .

## Iv Many closed channels.

In many physical systems coupling to a single closed channel is not sufficient to describe the scattering and bound states. In fact, coupled-channels calculations often show many overlapping resonances due to couplings with several closed channels which have poles near one another as a function of Takekoshi et al. (2012); Berninger et al. (2013); Gross et al. (2011). It is noteworthy that several resonances that have been used to study exotic three-body physics and the Efimov effect occur in regions with overlapping resonances Berninger et al. (2011); Gross et al. (2011). When resonances appear near one another, the simple formula (20) fails to describe the scattering length properly. However, it is straightforward to extend the results from the previous section by adding additional closed channels to the model. We thus start from the quantum defect matrix of the form

(23) |

It is in fact usually not necessary to include couplings between the closed channel states. We can instead assume that we use a basis in which the closed channels are already diagonalised Mies (1968). Then by choosing the reference potentials to reproduce the scattering lengths we get rid of the diagonal terms in the quantum defect matrix as well. As a result of a similar procedure as before and simple algebraic transformations, we obtain the energy-dependent scattering length

(24) |

where the dependence on is contained in the expansion of each terms as in Eq. (13). The standard scattering length as is

(25) |

where the resonant pole terms,

imply the limit, the widths are defined as in Eq. (15), and , , and are defined for the open channel with background scattering length . By defining a resonance width and shift for each resonance as in Eqs. (19) and (22), we find

(26) |

The mutual influence of the resonances on one another is thus contained in the terms . We note that this influence is not due to direct coupling between the closed channels, but rather to indirect interaction via the open channel.

Several important comments are in order here. Firstly, from looking at the structure of Eq. (26) one may suppose that there may be no simple local parameter similar to any more, since all the widths are needed to fully describe each resonance. However, it is possible to rewrite Eq. (26) in the form of isolated terms with new parameters, so that it is possible to define a meaningful parameter for each resonance. One particularly interesting case is the interplay between a broad resonance and a very narrow one. This is a fairly common situation if a weak resonance of high partial wave character exists near a broad -wave resonance. Assuming two resonances with , Eq. (26) can be simplified to

(27) |

where renormalizes the width of the narrow resonance and denotes the two pole positions. The new width can be either larger or smaller than , depending on the background scattering length and the relative position of resonances. One can treat the narrow resonance as an isolated one and describe it by the parameter using the renormalized width. By rewriting Eq. (27) as

(28) |

the scattering length can be approximated as a “local” isolated narrow resonance near with

(29) |

where is the same whether the “local” background or “global” background is used.

In general, it is also possible to algebraically transform formula (26) into a product form which was previously derived using a simple model of coupled square wells Lange et al. (2009)

(30) |

The transformation can be done by noticing that both (26) and (30) may be represented in the form , where and are polynomials of th and th order in . By equating the coefficients in the polynomials one obtains a set of equations connecting both formulas, which can be solved numerically. The resonance positions are given by the zeros of the denominators in (26). However, this usually does not determine the parameters uniquely. The most “intuitive” solution gives as the distance between the pole and the nearest zero of the scattering length, but other solutions are also possible. Consequently, there is no clear interpretation of the parameters as resonance widths. Remarkably, each and resonance position is a function of all the bare widths , crossing positions and the background scattering length in Eq. (26). Thus, one can always use Eq. (30) to define a pole strength for a “local” pole as in Eq. (28) as a product of a “local” width and a “local” background scattering length,

(31) |

Both and are functions of all the other pole terms (that is, all the parameters), and only their product remains well-determined when the fitting range eliminates distant poles that affect the local region.

Another important question is how to extract multiple resonance parameters by fitting numerical coupled channel calculations of . In some cases, such as Li in the spin channel, there are only two resonances Gross et al. (2011), and fitting the formula (26) is relatively easy. In other cases, such as cesium in its spin channel, the number of resonances is very high Berninger et al. (2013), and the fitting becomes computationally costly. Excluding some resonances from the fitting will matter for the uniqueness of the fit, as Eq. (26) is nonseparable. The effect is larger if any omitted resonances are quite broad or if they lie close to the region of interest. A parameter,

(32) |

which can be numerically determined by separately fitting to Eqs. (26) and (30), can be introduced as a measure of the impact of other resonances on the th one. In any case, the product obtained from the fitting tends to be robust, even if the individual terms are not. Consequently, a meaningful pole strength can be defined for each resonance.

Finally, it should be noted that our general formula in Eq. (24) shows how to obtain the scattering properties at finite energies away from . Thus, it is more general than the expressions in Eqs.(26) or (30) and can yield effective range or other, more rigorous, finite energy corrections in the presence of single or multiple resonances.

## V Applications

Precise Feshbach spectroscopy has been performed for many alkali metal species, and we will use examples of overlapping resonances here that have been studied experimentally for ultracold Li Gross et al. (2009, 2011) and cesium Chin et al. (2004); Berninger et al. (2013). Both of these species have been used in experimental studies of exotic three-body Efimov physics Kraemer et al. (2006); Ferlaino et al. (2011); Berninger et al. (2011); Gross et al. (2009, 2010b, 2011), where it is important to understand the character of the resonances and the mapping of the experimental field to scattering length. In order to demonstrate in practice how several overlapping resonances can be described, we apply our MQDT formulas to analyze the results of numerical coupled channels calculations using full Hamiltonian models that have been calibrated to reproduce a variety of experimental data. To account for the changing spin character of the channel states at moderate magnetic fields as well as the influence of resonances which were not included in the fit, we allow a small linear variation with in as a first-order correction. We used simple least-square fitting procedures that converge slowly and have to be performed carefully for nonlinear and diverging functions.

Figures 1 and 2 give the results of our fitting. Figure 1 shows the two -wave resonances in collision of two Li atoms in the spin channel. For this system the characteristic length ( is the Bohr radius) for 1393.39 atomic units Yan et al. (1996) (). The coupled channels scattering model was derived from our very accurate model for the collision of Li atoms Zürn et al. (2013). We fit the binding energy data for the Li spin channel from Ref. Dyke et al. (2013), using a changed scattering length from that of the Li potential to account for the failure of isotopic mass scaling, and find coupled channels pole positions at G and G Gau () for two , atoms, in excellent agreement with the experimental values, 844.9(8) G and 893.7(4) G, reported by Gross et al. Gross et al. (2011). Our least squares fitting of Eq. (26) to a discrete set of coupled channels calculations tabulated on a 1 G grid between 750 G and 950 G yields the same two positions, with G, G, widths G, G and the background scattering length (the definitions ensure that and have the same sign to fulfill the requirement that be positive definite Chin et al. (2010)). These four and values found are insensitive to the fit range used. Note that differs in each case from the pole position because of the large shifts involved. Figure 1 also shows the relative deviation of the fit from the numerics, indicating that the fit quality is better than everywhere. A fit of the same quality is obtained if the product form in Eq. (30) is used. Using MHzG, we obtained the parameters for these resonances equal to for the narrow one and for the wide one. Calculating the parameter defined as in Eq. (32) gives and .

In the case of cesium, for which and atomic units Berninger et al. (2013), many overlapping resonances in various partial waves (up to the -wave) were observed, and several regions of overlapping resonances are present, some quite complex Berninger et al. (2013). We focus here on three particular regions of magnetic field for which interesting three-body features have also been reported Kraemer et al. (2006); Ferlaino et al. (2011); Berninger et al. (2011). The first is a set of four -wave resonances near G, the second involves an -wave and two narrow -wave resonances near G, and finally the third region has an extremely broad -wave resonance near G with a narrow -wave one on its shoulder not ().

Figure 2 shows a comparison of the fits obtained by doing a least-squares fit of Eq. (26) to the coupled channels function from the Supplemental Material for Ref. Berninger et al. (2013). Table I lists the fitting parameters found using the respective fitting ranges G to G, G to G, and G to G. Note that the background scattering lengths obtained from locally fitting the resonances in each region differ from the global , which is expected to be on the order of . This means that other resonances not included in the fit are near enough to affect the local background and widths. Nevertheless, the fit quality in all cases stays at the level. A comparable quality fit is obtained by fitting to the product form in Eq. (30). The resonance at G is especially interesting here, as it is broadened by the presence of two stronger ones, that is, its nearby zero-crossing is shifted G from the pole of the resonance. This is consistent with the value of 12.98, where from Eq. (32) is an analogue of in Eq. (27) for a narrow resonance affected by two broad ones.

Table I also shows the and parameters describing each of the resonances. The parameter defined for each resonance is insensitive to fitting ranges. Our respective values of of and for the G and G strong -wave resonances compare with the values and estimated by Chin et al. Chin et al. (2010) from an earlier coupled channels model fit to Eq. (1) for a single isolated resonance. The parameter for the broad resonances tend to be the order of unity. However, the narrower resonances on their shoulders show large departures of from unity, implying a renormalization of the pole strength due to interactions among the resonances. We also found that the 7 G -wave and G -wave resonances, both being narrow ones sitting on the shoulder of the nearby broad -wave one, have respective values of 1.6 and 16, differing by an order of magnitude. These values are relevant to explaining the three-body physics that has been explored in the vicinity of these resonances Ferlaino et al. (2011); Berninger et al. (2011); Wang and Julienne (2013).

When the formula in Eq. (27) was applied to each of the narrow resonances near G and G, with respective parameters of and , the fit quality was only to accurate. However, we found that by using a more complex formula that takes into account the small change in the width of the broad resonance due to its interaction with the narrow ones, we were able to obtain a fit accurate to approximately , comparable to that obtained with Eqs. (26) or (30).

[G] | [G] | ||||||

492.68 | -wave | 493.315 | 0.1518 | 12.98 | 3.91 | 56 | |

495.04 | -wave | 499.617 | 3.397 | 1.05 | 3.84 | 100 | |

501.44 | -wave | 502.186 | 2.088 | 0.093 | 3.79 | 5.4 | |

505.38 | -wave | 505.483 | 0.2389 | 0.194 | 3.85 | 1.3 | |

544.19 | -wave | 544.108 | 0.02179 | 4.95 | 1.52 | 1.0 | |

548.79 | -wave | 556.693 | 6.626 | 1.01 | 3.90 | 160 | |

554.07 | -wave | 553.793 | 0.5048 | 0.387 | 1.34 | 1.6 | |

786.17 | -wave | 884.638 | 92.26 | 0.856 | 3.66 | 1480 | |

820.32 | -wave | 819.386 | 0.4508 | 3.38 | 2.05 | 16 |

## Vi Conclusion

In conclusion, we presented a simple analytical model to describe the variation of the -wave scattering length with magnetic field when there is an arbitrary number of Feshbach resonances. We introduced simple parameters characterizing the system, analogous to the resonance width and background scattering length in the two-channel case. We discussed the non-separability of the scattering length formula for overlapping resonances and provided examples where we accurately reproduced coupled channels numerical results and found the resonance parameters. Apart from characterization of overlapping sets of resonances, our model should be quite helpful for precise mapping of the scattering length to the laboratory field, which is a critical aspect in the interpretation of experiments with three-body recombination and Efimov physics.

This work was supported by the Foundation for Polish Science International PhD Projects Programme co-financed by the EU European Regional Development Fund, by AFOSR MURI Grant No. FA9550-09-1-0617, and in part by the National Science Foundation under Grant No. NSF PHY11-25915.

## References

- Chin et al. (2010) C. Chin, R. Grimm, P. S. Julienne, and E. Tiesinga, Rev. Mod. Phys. 82, 1225 (2010).
- Inouye et al. (1998) S. Inouye, M. R. Andrews, J. Stenger, H.-J. Miesner, D. M. Stamper-Kurn, and W. Ketterle, Nature 392, 151 (1998).
- Cornish et al. (2000) S. L. Cornish, N. R. Claussen, J. L. Roberts, E. A. Cornell, and C. E. Wieman, Phys. Rev. Lett. 85, 1795 (2000).
- Bourdel et al. (2004) T. Bourdel, L. Khaykovich, J. Cubizolles, J. Zhang, F. Chevy, M. Teichmann, L. Tarruell, S. J. J. M. F. Kokkelmans, and C. Salomon, Phys. Rev. Lett. 93, 050401 (2004).
- Zwierlein et al. (2004) M. W. Zwierlein, C. A. Stan, C. H. Schunck, S. M. F. Raupach, A. J. Kerman, and W. Ketterle, Phys. Rev. Lett. 92, 120403 (2004).
- Jochim et al. (2003) S. Jochim, M. Bartenstein, A. Altmeyer, G. Hendl, S. Riedl, C. Chin, J. Hecker Denschlag, and R. Grimm, Science 302, 2101 (2003).
- Lang et al. (2008) F. Lang, K. Winkler, C. Strauss, R. Grimm, and J. Hecker Denschlag, Phys. Rev. Lett. 101, 133005 (2008).
- Moerdijk et al. (1995) A. J. Moerdijk, B. J. Verhaar, and A. Axelsson, Phys. Rev. A 51, 4852 (1995).
- Timmermans et al. (1999) E. Timmermans, P. Tommasini, M. Hussein, and A. Kerman, Phys. Rep. 315, 199 (1999).
- Kraemer et al. (2006) T. Kraemer, M. Mark, P. Waldburger, J. G. Danzl, C. Chin, B. Engeser, A. D. Lange, K. Pilch, A. Jaakkola, H.-C. Nägerl, et al., Nature 440, 315 (2006).
- Knoop et al. (2009) S. Knoop, F. Ferlaino, M. Mark, J. G. Danzl, T. Kraemer, H.-C. Nägerl, and R. Grimm, Nature Phys. 5, 227 (2009).
- Zaccanti et al. (2009) M. Zaccanti, B. Deissler, C. DâErrico, M. Fattori, M. Jona-Lasinio, S. Müller, G. Roati, M. Inguscio, and G. Modugno, Nature Physics 5, 586 (2009).
- Gross et al. (2010a) N. Gross, Z. Shotan, S. Kokkelmans, and L. Khaykovich, Phys. Rev. Lett. 105, 103203 (2010a).
- Ferlaino et al. (2011) F. Ferlaino, A. Zenesini, M. Berninger, B. Huang, H. C. Naegerl, and R. Grimm, Few Body Sys. 51, 113 (2011).
- Berninger et al. (2011) M. Berninger, A. Zenesini, B. Huang, W. Harm, H.-C. Nägerl, F. Ferlaino, R. Grimm, P. S. Julienne, and J. M. Hutson, Phys. Rev. Lett. 107, 120401 (2011).
- Berninger et al. (2013) M. Berninger, A. Zenesini, B. Huang, W. Harm, H.-C. Nägerl, F. Ferlaino, R. Grimm, P. S. Julienne, and J. M. Hutson, Phys. Rev. A 87, 032517 (2013).
- Dyke et al. (2013) P. Dyke, S. E. Pollack, and R. G. Hulet, arXiv:1302.0281 (2013).
- Mies (1984) F. H. Mies, J. Chem. Phys. 80, 2514 (1984).
- Mies and Julienne (1984) F. H. Mies and P. S. Julienne, J. Chem. Phys. 80, 2526 (1984).
- Julienne and Mies (1989) P. S. Julienne and F. H. Mies, J. Opt. Soc. Am. B 6, 2257 (1989).
- Mies and Raoult (2000) F. H. Mies and M. Raoult, Phys. Rev. A 62, 012708 (2000).
- Julienne and Gao (2006) P. S. Julienne and B. Gao, in Atomic Physics 20, edited by C. Roos, H. Häffner, and R. Blatt (AIP, Melville, New York, 2006), pp. 261–268.
- Julienne (2009a) P. S. Julienne, Chapter 6 of Cold Molecules: Theory, Experiment, Applications, ed. by R. V. Krems, W. C. Stwalley, B. Friedrich, CRC Press (arXiv:0902.1727) pp. 221–243 (2009a).
- Lange et al. (2009) A. D. Lange, K. Pilch, A. Prantner, F. Ferlaino, B. Engeser, H.-C. Nägerl, R. Grimm, and C. Chin, Phys. Rev. A 79, 013622 (2009).
- Köhler et al. (2006) T. Köhler, K. Góral, and P. S. Julienne, Rev. Mod. Phys. 78, 1311 (2006).
- Burke et al. (1998) J. P. Burke, C. H. Greene, and J. L. Bohn, Phys. Rev. Lett. 81, 3355 (1998).
- Gao (2008) B. Gao, Phys. Rev. A 78, 012702 (2008).
- Gao (2011) B. Gao, Phys. Rev. A 83, 062712 (2011).
- Idziaszek et al. (2011) Z. Idziaszek, A. Simoni, T. Calarco, and P. S. Julienne, New Journal of Physics 13, 083005 (2011).
- Gao (1998) B. Gao, Phys. Rev. A 58, 1728 (1998).
- Gao (2000) B. Gao, Phys. Rev. A 62, 050702 (2000).
- Mies et al. (2000) F. H. Mies, E. Tiesinga, and P. S. Julienne, Phys. Rev. A 61, 022721 (2000).
- Nygaard et al. (2006) N. Nygaard, B. I. Schneider, and P. S. Julienne, Phys. Rev. A 73, 042705 (2006).
- Gribakin and Flambaum (1993) G. F. Gribakin and V. V. Flambaum, Phys. Rev. A 48, 546 (1993).
- Julienne (2009b) P. S. Julienne, Faraday Discuss. 142, 361 (2009b).
- Bolda et al. (2002) E. L. Bolda, E. Tiesinga, and P. S. Julienne, Phys. Rev. A 66, 013403 (2002).
- Krych and Idziaszek (2009) M. Krych and Z. Idziaszek, Phys. Rev. A 80, 022710 (2009).
- Idziaszek and Julienne (2010) Z. Idziaszek and P. S. Julienne, Phys. Rev. Lett. 104, 113202 (2010).
- Wang and Julienne (2013) Y. Wang and P. S. Julienne, unpublished (2013).
- Takekoshi et al. (2012) T. Takekoshi, M. Debatin, R. Rameshan, F. Ferlaino, R. Grimm, H.-C. Nägerl, C. R. Le Sueur, J. M. Hutson, P. S. Julienne, S. Kotochigova, et al., Phys. Rev. A 85, 032506 (2012).
- Gross et al. (2011) N. Gross, Z. Shotan, O. Machtey, S. Kokkelmans, and L. Khaykovich, Comptes Rendus Phys. 12, 4 (2011).
- Mies (1968) F. H. Mies, Phys. Rev. 175, 164 (1968).
- Gross et al. (2009) N. Gross, Z. Shotan, S. Kokkelmans, and L. Khaykovich, Phys. Rev. Lett. 103, 163202 (2009).
- Chin et al. (2004) C. Chin, V. Vuletić, A. J. Kerman, S. Chu, E. Tiesinga, P. J. Leo, and C. J. Williams, Phys. Rev. A 70, 032701 (2004).
- Gross et al. (2010b) N. Gross, Z. Shotan, S. Kokkelmans, and L. Khaykovich, Phys. Rev. Lett. 105, 103203 (2010b).
- Yan et al. (1996) Z.-C. Yan, J. F. Babb, A. Dalgarno, and G. W. F. Drake, Phys. Rev. A 54, 2824 (1996).
- Zürn et al. (2013) G. Zürn, T. Lompe, A. N. Wenz, S. Jochim, P. S. Julienne, and J. M. Hutson, Phys. Rev. Lett. 110, 135301 (2013).
- (48) Units of gauss rather than tesla, the accepted SI unit for the magnetic field, have been used in this paper to conform to the conventional usage in this field of physics.
- (49) Following the usage of Ref. Chin et al. (2010), an -wave resonance means an -wave bound states is coupled to an -wave open channel.