A system of fermions with a short-range interaction proportional to the scattering length is studied at finite density. At any order , we evaluate the complete contributions to the energy per particle arising from combined (multiple) particle-particle and hole-hole rescatterings in the medium. This novel result is achieved by simply decomposing the particle-hole propagator into the vacuum propagator plus a medium-insertion and correcting for certain symmetry factors in the -th power of the in-medium loop. Known results for the low-density expansion up to and including order are accurately reproduced. The emerging series in can be summed to all orders in the form of a double-integral over an arctangent function. In that representation the unitary limit can be taken and one obtains the value for the universal Bertsch parameter. We discuss also applications to the equation of state of neutron matter at low densities and mention further extensions of the resummation method.
Resummation of fermionic in-medium ladder diagrams to all orders
Physik-Department T39, Technische Universität München, D-85747 Garching, Germany
PACS: 05.30.Fk, 12.20.Ds, 21.65+f, 25.10.Cn
1 Introduction and summary
Dilute degenerate many-fermion systems with large scattering lengths are of interest e.g. for modeling the low-density behavior of nuclear or neutron star matter. Also dilute systems of ultracold atoms can nowadays be trapped. Because of the possibility to tune (magnetically) atomic interactions through so-called Feshbach resonances, ultracold fermionic gases provide an exceptionally valuable tool to explore the (non-perturbative) many-body dynamics at weak and strong coupling together with the transition from the superconducting to the Bose-Einstein condensed state. Of particular interest in this context is the so-called unitary limit in which the two-body interaction has the critical value to support a boundstate at zero energy. In that special situation the s-wave scattering length diverges, , and the (strongly) interacting many-fermion system becomes scale invariant. Its groundstate energy is then determined by a single universal number, the so-called Bertsch parameter , which measures the ratio of the energy per particle, , to that of a free Fermi gas, . Here, denotes the Fermi momentum and stands for the (heavy) fermion mass.
The calculation of is an intrinsically non-perturbative problem which has been approached in recent years by numerical quantum Monte-Carlo simulations (in a periodic box). As the state of the art, at present a value of seems to emerge from these calculations , presumably with still debatable error bars due to finite-size corrections etc. The equation of state of neutron matter at low densities has also been studied using quantum Monte-Carlo techniques [2, 3]. Due to the very large neutron-neutron scattering length, fm, neutron matter at low densities, fm, is supposed to be a fermionic gas close to the unitarity limit. The results of a variety of sophisticated many-body calculations are summarized in Figs. 3,4 of ref.  and these give indications for a value of .
On the other hand effective field theory methods have been used to rederive and systematically improve the low-density expansion for a system of fermions with short-range interactions [4, 5, 6]. In particular, effective range corrections and (logarithmic) contributions from three-particle scattering have been included. In an unpublished paper, Steele  has computed several higher order contributions proportional to and . These are generated either by multiple particle-particle and hole-hole rescatterings in the medium or they arise from particle-hole ring diagrams. Based on some numerical evidence, Steele argued that in the limit of large space-time dimensions the hole-hole contributions would be suppressed. Under this assumption the resummation of the particle-particle contributions becomes possible (in the limit ) in the form of a simple geometrical series: . (We are choosing the sign-convention such that a positive scattering length corresponds to attraction.) The Bertsch parameter following from Steele’s approximation is , surprisingly close to recent quantum Monte-Carlo results. The validity of Steele’s arguments concerning the expansion in has been critically reassessed in the more elaborate work by Schäfer et al. . There it has been shown that if the strong coupling limit is taken after the limit the universal Bertsch parameter is (see eq.(44) in ref. ).
Moreover, the particle-particle ladders (for ) have been resummed in ref.  in the form of a phase space integral over a geometrical series, , and the corresponding Bertsch parameter had the (rather small) value . In the same paper an analogous expression for the sum of all hole-hole ladders has been given (see eq.(25) in ref. ). In that case the geometrical series involves two subtractions (hole-hole ladders start to contribute at order ) and the unitary limit does not exist for this term. The same feature applies to the sum of all particle-hole ring diagrams, which have been discussed in section IV of ref. . Beyond these partial resummation results there exist still large classes of diagrams with mixed particle-particle and hole-hole ladders which have not been considered in ref.  (or elsewhere in the literature). Clearly, these should all be resummed in order to see whether they contribute in the limit to the universal Bertsch parameter . The first mixed ( and ) ladder appears at order and it has been considered by Steele . Hammer et al.  have carefully checked Steele’s calculation (for ) and they find differences for just this term. In particular, the numerical value given for it in eq.(16) of ref.  (middle term) has to be corrected and multiplied by a factor 2. Actually, Hammer et al.  find a simpler analytical representation for it in the form of a phase space integral over the squared particle-particle bubble times the hole-hole bubble, (see eq.(7) in section 3). Besides this leading term, results for mixed particle-particle and hole-hole ladders at higher orders are not known. For large their classification and combinatorial counting is already a non-trivial problem.
The purpose of the present paper is to close this gap. We will evaluate, at any order , the complete contributions to the energy per particle arising from combined (multiple) particle-particle and hole-hole rescatterings in the medium. A key to the solution of this problem is a different organization of the many-body calculation from the start. Instead of treating (propagating) particles and holes separately, we keep them together and measure the difference to the propagation in vacuum by a “medium-insertion”. The latter involves a delta-function for on-shell kinematics and a step-function which restricts momenta to the region inside the Fermi sphere. In that organizational scheme the pertinent in-medium loop (or in-medium bubble) is complex-valued. The contribution to the energy per particle at order is therefore not obtained directly from the -th power of the in-medium loop. However, after reinstalling the symmetry factors which belong to diagrams with double medium-insertions, a real-valued expression is retained for all . Known results for the low-density expansion up to and including order are accurately reproduced in our scheme. The emerging series in can even be summed to all orders in the form of a double-integral over an arctangent function. In that explicit representation the unitary limit can be taken straightforwardly and one finds the value for the universal Bertsch parameter. As an application we discuss the equation of state of neutron matter at low densities and finally we mention further possible extensions of the resummation method.
2 Preparation: In-medium propagator
We are interested in the equation of state of a (non-relativistic) many-fermion system with a short-range two-body interaction proportional to the scattering length . In perturbation theory the interaction contributions to the energy density are represented by closed multi-loop diagrams with a certain number of contact vertices. Their evaluation proceeds via Feynman rules which introduce a factor for each interaction vertex and the (non-relativistic) particle-hole propagator:
for an internal fermion-line carrying energy and momentum . The (large) fermion mass is denoted by . The identity in the second line of eq.(1) gives the separation of the in-medium propagator into the vacuum propagator and a “medium-insertion”. Note that by construction a medium-insertion puts a particle onto the mass-shell and restricts its momentum to the interior of the Fermi sphere . The Fermi momentum is related to the density by . We are considering here only the case with spin-degeneracy factor . Relevant many-body contributions to the energy density come from diagrams with at least two medium-insertions. In case of the closed ladder diagram shown in Fig. 1 this minimal pair of medium-insertions has to placed on adjacent positions of the double-ring, for the following reason. After opening one gets a planar ladder diagram whose energy denominators are all equal to differences of fermion kinetic energies. Only for this planar topology the resulting factors of (from the energy denominators) will balance the factors from the interaction vertices such that a finite result remains (at each order ) in the non-relativistic limit. The open ladder diagram shown in Fig. 1 comes in all possible variations with further medium-insertions on internal lines. Due to the momentum-independence of the contact interaction all these multi-loop diagrams factorize (successive loops have nothing in common) and they can therefore be summed together in the form of a power of the in-medium loop.
3 In-medium loop generated by a contact interaction
The basic ingredient for calculating the energy density arising from ladder diagrams with a short-range contact interaction has been identified as the in-medium loop (or the in-medium bubble). As enforced by the two minimal medium-insertions the incoming momenta and are from the region below the Fermi surface . It is convenient to introduce their half sum and half difference . In our ordering scheme the in-medium loop is built up from contributions with zero, one, and two medium-insertions: . The corresponding one-loop diagrams are shown in Fig. 2 and their evaluation is exhibited now in detail.
The contribution with zero medium-insertions is the well-known rescattering bubble in vacuum. It is normalized relative to the tree-level interaction and therefore it includes only one factor . The evaluation of the loop integral proceeds as follows:
using residue calculus for the energy integral and in the last step the rule of dimensional regularization has been applied. In other regularization schemes the occurring scale-dependent constant is absorbed into to define the renormalized (physical) scattering length [7, 8]. The resummation of infinitely many rescatterings in the vacuum in the form of a geometrical series leads to the unitarized scattering length approximation:
which implies the relation for the s-wave scattering phase shift .
Next, we turn to the term from diagrams with one medium-insertion. There are two equal contributions which differ in their representation merely by the sign of the loop momentum and together they read:
The integral over a shifted Fermi sphere is most suitably performed by using spherical coordinates with the range of the magnitude . The quantity is a directional cosine which covers the full range since the displacement vector satisfies the condition (see also the left part of Fig. 3). Carrying out this procedure one obtains the following result for the real part:
with the logarithmic function:
written in terms of the two dimensionless variables and . Since both external momenta are from inside the Fermi sphere one has the additional constraint . It is worth to note that the function in eq.(6) is equal to the sum of the particle-particle bubble and the hole-hole bubble, , with the unusual feature that the latter is also taken at momenta below the Fermi surface. For comparison the particle-particle bubble reads [7, 8]:
Finally, we come to the right diagram in Fig. 2 with two-medium insertions. Obviously, it generates a purely imaginary contribution. The total imaginary part of the in-medium loop has the following representation:
where we have suitably arranged terms with no, one, and two step-functions . The first term in eq.(8) of the form makes no contribution to the imaginary part since the corresponding phase space is completely Pauli-blocked: . This equation expresses the obvious fact that on-shell scattering of two particles from below the Fermi surface into the region above the Fermi surface is not possible due to energy conservation. Thus there remains the imaginary part due to the second term in eq.(8). After visualizing the occurring product of - and -functions one sees that the imaginary part has a nice geometrical interpretation: namely as times that part of the solid angle of a (centered) sphere of radius which lies inside the intersection region of two spheres of radius with their centers displaced by . The corresponding configuration of spheres is sketched in the right part of Fig. 3. After inclusion of the appropriate prefactor the result for the imaginary part of the in-medium loop reads:
with the (non-smooth) function:
where lies in the interval . It is interesting to observe that the diagram with two medium-insertions alone gives twice the total imaginary part. Putting the real and imaginary pieces together the complex-valued in-medium loop reads:
and if the contribution from the diagram with two medium-insertions is taken out, the imaginary part of that same expression changes sign:
Actually, this special property of the in-medium loop turns out to be crucial in order to derive the correct expression for the energy per particle from powers of the (complex-valued) in-medium loop. The derivation of the energy per particle is the topic of the next section.
4 Energy per particle
In this section we show how the contributions to the energy per particle at any order can be constructed from the in-medium loop and how the emerging series in can be summed to all orders. Consider the (open) ladder diagram with contact interactions. It is given by the -th power of the in-medium loop times a factor . Closing the two open fermion-lines introduces an integration over the allowed phase space . The emerging integrand for the energy density at order would be complex-valued and therefore it cannot yet be the correct one. The deficit of the (naive) iteration method at this intermediate stage becomes evident if one draws all diagrams with (repeated) pairs of adjacent medium-insertions. The corresponding set of diagrams at fourth order is shown in Fig. 4. These (decorated) diagrams have additional symmetry factors which are not respected by the binomial expansion of . Furthermore, the same symmetry factors correct for the overcounting of certain diagrams as introduced by the terms of the binomial series. As a result of this combinatorial analysis one has to reweight in the binomial series expansion of the -th power of coming from the diagrams with repeated double medium-insertions with the appropriate symmetry factor . This crucial amendment leads to the following summation formula:
Note that the identity introduces the binomial coefficients for the -th power of a sum, and in this way one can easily reproduce the result on the right hand side of eq.(13). In the diagrammatic representation of the energy density, is the number of different possibilities to attach double medium-insertions on a ring with segments and the additional factor comes from the rotations which transform the ring into itself.
Inspection of the right hand side of eq.(13) shows that the resulting
homogeneous polynomials of degree in and are manifestly real
for all . The first five terms read: : , : , :
, : , : .
In Fig. 4 and the appended caption it is shown explicitly how the real-valued expression results from the sum of (complex-valued) diagrams once their symmetry factors are taken into account. The same graphology gives at second order : , and at third order : .
At this point we have achieved a representation which includes, at any order , the complete contributions from all mixed particle-particle and hole-hole ladders. One can even go further and sum up the whole series of ladder diagrams to all orders. The pertinent series can be solved easily by a (complex) logarithm. After inclusion of the kinetic energy and the Fock exchange-term (which reduces the Hartree term discussed so far by a factor ) the complete expression for the energy per particle reads:
The occurring arctangent function refers to the usual branch with odd parity, , and values in the interval . Other branches of the arctangent function are excluded by the weak coupling limit , which has to give zero independent of the sign of the scattering length . Note that there occurs a discontinuity (by an amount ) when the denominator passes through zero from positive to negative values. This happens for all positive values of since the function has a logarithmic singularity at . The integral-representation for given in eq.(14) has some similarity with the expression one obtains from the resummed particle-hole ring diagrams (see eq.(29) in ref. ). In that case the arctangent function involves two subtractions (particle-hole ring diagrams start to contribute at order ) and the integral extends over a 4-momentum transfer . In the present case the reduction to a double-integral (over the quarter unit-disc) has been obtained by employing the following master formula for integrals over the interior of two Fermi spheres:
where and . Surprisingly, the imaginary part function defined in eq.(10) occurs here as the pertinent weighting function. Integrals over the outer region of two Fermi spheres can be reduced in a similar way:
5 Expansion in powers of
Several orders in the low-density expansion of the energy per particle are known . The terms stemming from particle-particle and hole-hole ladders can be used as a check of our calculation which is organized differently by not treating separately particles and holes. We find up to and including fourth order:
As it must be, the linear and quadratic coefficients agree analytically. Since only double-integrals are involved the other numerical coefficients can be obtained with high precision. One finds again good agreement with existing calculations [4, 7, 9]. The third order coefficient is: , where the numbers on the left hand side correspond to our separation and those on the right side refer to the sum of two-fold particle-particle and two-fold hole-hole rescatterings. The same comparison for the fourth order coefficient gives: , with our separation versus the sum of triple particle-particle, triple hole-hole, and combined particle-particle and hole-hole scatterings. At this point it is important to note that a factor of is missing in Steele’s numerical result for the latter contribution (see middle term in eq.(16) of ref. ). This error has been confirmed by Hammer et al.  who have carefully checked Steele’s calculation. It has become evident that the squared imaginary part of the in-medium loop includes also important many-body correlation effects. Their precise mapping into the (traditional) particle-hole counting scheme is not obvious.
As mentioned in the introduction, Steele  has suggested a resummation to all orders in form of a simple geometrical series:
Although the original arguments (via a -expansion) have been (partly) invalidated , it may serve as a useful reference, in particular since the associated Bertsch parameter comes out fairly realistic. We adapt the expansion of the energy per particle in powers of :
such that all coefficients become for Steele’s resummation result. Performing the same expansion with our complete expression for given in eq.(14) we find the following values for the first dozen expansion coefficients :
One observes that the first few coefficients stay below while the higher ones show a tendency to increase appreciably with . The apparent region of “convergence” of the power series in can therefore roughly be estimated as . Despite various attempts we have not succeeded to derive the exact value of in terms of (reasonable) mathematical constants. The value represents an extremely accurate analytical approximation, which is far more precise than all previous determinations of this coefficient [4, 7].
For large scattering lengths the low-density expansion is of very limited validity. In that case one has to take the full expression for the energy per particle as given by eq.(14). The solid line in Fig. 5 shows its ratio to the (free) Fermi gas energy as a function of the dimensionless parameter . The behavior of this ratio is demonstrated for both signs of the scattering length ( for attraction and for repulsion). Outside the region the curve is almost flat while a sharp peak develops inside with a maximum value of about at . The dashed curve in Fig. 5 corresponds to Steele’s  resummation via a geometrical series eq.(18). The behavior on the attractive side () is quite similar while on the repulsive side () the artificial pole at causes essential differences.
6 Unitary limit
The unitary limit is of special interest since in this limit the strongly interacting many-fermion system becomes scale invariant. The energy per particle is then determined by just a constant multiple of the (free) Fermi gas energy:
with the so-called Bertsch parameter. Returning to the expression for the energy per particle given in eq.(14) one sees that the unitary limit can be performed straightforwardly. The formula for calculating the Bertsch parameter reads:
The resulting numerical value is to be compared with obtained by Schäfer et al.  from the resummation of particle-particle ladders. Actually, is defined by the principal value integral and the necessity to treat the pole-singularity limits the precision in the numerical computation of this number. One observes that the additional mixed particle-particle and hole-hole ladders increase the Bertsch parameter by more than a factor . Moreover, the divergence of the subset of hole-hole ladders (for ) encountered in ref.  has disappeared. Clearly, the value as obtained here via an analytical calculation rooted in perturbation theory is still considerably larger than the result from recent quantum Monte-Carlo simulations .
In the work by Haussmann et al.  a self-consistent, thermodynamically consistent treatment of the unitary Fermi gas at finite temperatures has been presented. In their non-perturbative approach the exact one- and two-particle Green functions serve as an infinite set of variational parameters and extensive numerical work enters into the solutions of the stationarity constraints and the thermodynamical potentials . The resulting Bertsch parameter at zero temperature was found to be . This value is remarkably close to the experimental determinations by Bartenstein et al. , , and Bourdel et al. , .
7 Application to neutron matter and outlook
As an application of our analytical result eq.(14) for the complete resummation of in-medium ladder diagrams we consider the equation of state of neutron matter. Due to the very large neutron-neutron scattering length fm [13, 14] neutron matter at low densities is supposed to be a Fermi gas close to the unitary limit. Recent quantum Monte Carlo simulations [2, 3] based e.g. on the Argonne nucleon-nucleon potential give some indication for such a behavior.
In Fig. 6 we show the energy per particle of neutron matter as a function of the neutron density . The solid line results from our analytical formula eq.(14) inserting fm for the scattering length and MeV for the fermion mass. The dash-dotted line corresponds to the sophisticated many-body calculation by the Urbana group , to be considered as representative of realistic neutron matter calculations. The dashed line in Fig. 6 reproduces Steele’s suggestion eq.(18) in the form of a simple geometrical series. One observes good agreement up to rather high neutron densities of fm, where the dimensionless parameter reaches values up to . At higher neutron densities repulsive effects from three-nucleon forces (which are included in the Urbana calculation ) start to play a more significant role. The inclusion of the effective range fm for s-wave -scattering can also lead to sizeable changes in the equation of state, as demonstrated in ref.. Note however, that only the resummed particle-particle ladders have been considered in that work.
Fig. 7 shows the neutron matter equation of state in a different representation by plotting the ratio of the energy per particle to the (free) Fermi gas energy against the dimensionless parameter . The dots in this figure reproduce results of various (quantum Monte-Carlo) calculations of low-density neutron matter and have been taken over from Figs. 3,4 in ref. (see original references therein). The solid and dashed curves are the same as in the right part of Fig. 5 only continued further out in the parameter .
The complete resummation of ladder diagrams (with a short-range contact interaction proportional to the s-wave scattering length ) as achieved in this work suggests numerous possible extensions. The s-wave effective range parameter and the p-wave scattering volumes should be included via terms in the contact interaction. Asymmetries with respect to the spin and/or isospin degrees of freedom can be studied via an appropriate modification of the medium-insertion. The generalization of the resummation method to finite temperatures would be an equally interesting project. Work along these lines in progress.
I thank J.W. Holt, A. Schwenk and W. Weise for many useful discussions. I thank H.W. Hammer for communicating to me results from own unpublished work which have been very valuable. This work is partially supported by the DFG Excellence Cluster “Origin and Structure of the Universe”.
-  M. MacNeil Forbes, S. Gandolfi, A. Gezerlis, “Resonantly interacting fermions in a box”, cond-mat/1011.2197.
-  J. Carlson, J. Morales, V.R. Pandharipande, D.G. Ravenhall, Phys. Rev. C68, 025802 (2003).
-  A. Gezerlis, J. Carlson, Phys. Rev. C81, 025803 (2005).
-  H.W. Hammer, R.J. Furnstahl, Nucl. Phys. A678, 277 (2000).
-  R.J. Furnstahl, J.V. Steele, N. Tirfessa, Nucl. Phys. A671, 396 (2000).
-  R.J. Furnstahl, H.W. Hammer, Ann. Phys. (NY) 302, 206 (2002).
-  J.V. Steele, “Effective field theory power counting at finite density”, nucl-th/0010066.
-  T. Schäfer, C.W. Kao, S.R. Cotanch, Nucl. Phys. A762, 82 (2005).
-  H.W. Hammer et al., unpublished; and private communications.
-  R. Haussmann, W. Rantner, S. Cerrito, W. Zwerger, Phys. Rev. A74, 023610 (2007).
-  M. Bartenstein et al., Phys. Rev. Lett. 92, 120401 (2004).
-  T. Bourdel et al., Phys. Rev. Lett. 93, 050401 (2004).
-  D.E. Gonzales Trotter et al., Phys. Rev. C73, 034001 (2006).
-  Q. Chen et al., Phys. Rev. C77, 054002 (2008).
-  A. Akmal, V.R. Pandharipande, D.G. Ravenhall, Phys. Rev. C58, 1804 (1998).
-  A. Schwenk, C.J. Pethick, Phys. Rev. Lett. 95, 160401 (2005).