Approximate selfenergy for Fermi systems with large swave scattering length: a step towards density functional theory.
Abstract
In the present work, we start from a minimal Hamiltonian for Fermi systems where the swave scattering is the only low energy constant at play. ManyBody Perturbative approach that is usually valid at rather low density is first discussed. We then use the resummation technique with the ladder approximation to obtain compact expressions for both the energy and/or the onshell selfenergy in infinite spindegenerated systems. Diagrammatic resummation technique has the advantage in general to be predictive in a region of density larger compared to manybody perturbation theory. It also leads to nondiverging limit as . Still, the obtained expressions are rather complex functional of the Fermi momentum . We introduce the full phasespace average or the partial phasespace methods respectively applied to the energy or to the selfenergy to simplify their dependences in terms of while keeping the correct limit at low density and the nondiverging property at large . Quasiparticle properties of Fermi system in various regime of density and scattering length are then illustrated. Our conclusion is that such simplified expressions where the direct link is made with the low energy constant without finetuning can provide a clear guidance to obtain density functional theory beyond the perturbative regime. However, quasiparticle properties close or near unitary cannot be reproduced unless this limit is explicitly used as a constraint. We finally discuss how such approximate treatment of quasiparticle can guide the development of density functional theory for strongly interacting Fermi systems.
1 Introduction
Strongly interacting manybody Fermi systems can sometimes be described by rather simple density functional theories (DFT). This is the case of nuclear systems where simple functionals, like those based on the Skyrme type contact interactions [1, 2, 3, 4], are nowadays widely used. With very few parameters, the functional can describe static, thermodynamical and dynamical properties very accurately in a unified framework. This is even more surprising in view of the complexity of the strong multibody interaction between nucleons. The question (a) “How such simplicity can emerge in strongly interacting Fermi systems?” is still largely open in the nuclear physics community (see discussion in [5]).
This situation is not unique in nature. Simple DFTs apply also to the case of Fermi gas at unitarity. This gases are characterized by infinite swave scattering length in the dilute regime. In this case, the energy becomes directly proportional to the free FermiGas energy. This situation can be seen as one of the simplest DFT one could ever imagine. Still, while in a DFT framework unitary gases can be described in a rather simplistic manner (see for instance [6, 7, 8]), their treatment starting from a particleparticle interactions, requires rather advanced manybody techniques like MonteCarlo (MC) methods [9, 10, 11, 12, 13, 14, 15], SelfConsistent Green Function (SCGF) [16, 17], Brueckner HartreeFock (BHF) [18], Blod Diagrammatic Monte Carlo [19, 20] eventually associated to resummation technique based on conformalBorel transformation [21]. These approaches generally rely on rather involved numerical methods and usually prevent from connecting analytically the energy with the lowenergy constants (LEC) associated to the underlying interaction. To our opinion, to reply to the question (a) it would be desirable to also give some hints on the other question (b) “Can we qualitatively or quantitatively connect the parameters of the DFTs used in strongly interacting Fermions with the lowenergy constants of the interaction?” Assuming that we can directly connect the parameters used in a functional to the LEC of the interaction, this would render the DFT completely nonempirical. This would also be at variance with the strategy used nowadays to construct a DFT. Indeed, most currently used DFTs in cold atoms or in nuclei are usually directly adjusted either on experimental datas or on pseudodatas obtained using one of the abinitio methods mentioned above. Such direct fitting procedure is very powerful because it includes automatically complex manybody correlations in the DFT. It also leads in general to a very precise description of the global properties of Fermi systems. This is for instance the case in atomic nuclei were the precision on the ground state energy is better than 2–3 for medium mass nuclei and goes down to 0.5 % for heavy systems (see for instance [22, 23]). This strategy has also some drawbacks. Among them, we usually face the difficulty that some components of the functional are not really constrained by the experiments. This is for instance the case of the density dependence of the symmetry energy in nuclei that is particularly important for the physics of exotic nuclei. Another example that was pointed out recently is the failure of empirical functionals to properly describe the low density limit of neutron matter [24]. One should mention that, with recent progresses in the nuclear interaction and in abinitio manybody techniques, there is an increasing interest in developing DFTs directly starting from a clear manybody framework. Among the recent works, we mention the Density Matrix Expansion (DME) proposed already some times ago [25, 26] that has reached now a certain level of maturity [27, 28, 29, 30, 31, 32]. Another clearly defined approach is to write the effective action and use the inversion method as proposed in [33, 34] (see also the recent interesting progress of Ref. [35]). Still, the quantitative description of strongly interacting systems beyond the low density limit and/or beyond the HartreeFock approximation is a rather difficult task.
For this reason, we explore here qualitatively how a DFT can be obtained for three dimensional infinite Fermi systems where the parameters of the DFT are directly linked to the LEC of the interaction. More specifically, we consider the simplified problem where the interaction is described by a single LEC, , and where the interaction strength can vary from the perturbative to the nonperturbative regime. Such physical situation was explored in different regimes using standard manybody techniques starting from an EffectiveFieldTheory approach [36, 37]. For instance, the low density limit was studied in Ref. [38]. This case is particularly highlighting since in this case, up to third order in perturbation, the energy can be written as a simple polynomial (and potentially polylogarithmic from fourth order) of where denotes the density. The perturbative approach breaks down when increases. In this case, DFT have also been obtained using diagrammatic resummation techniques [39, 40]. Both perturbation and resummation to obtain compact expressions of the energy in terms of will be briefly discussed here. As we will see, the bruteforce resummation however generally suffers from the lack of predictive power especially close to the unitary limit. Following Ref. [39], we show that, using a procedure called hereafter phasespace average approximation, the energy can be written as a simplified functional of that in addition improves the description of strongly interacting systems. The work of Ref. [39] was actually the starting point of several new developments in the nuclear manybody context. In Ref. [24], guided by the simplified expression of the energy a hybrid functional was proposed where some of the parameters are directly connected to . Similarly, in [8, 41], a nonempirical functional was proposed that could reproduce both cold atoms gases and neutron matter up to fm including the effective range effect. Such new functionals were also used in Ref. [41] (see Fig. 6 of this reference) to understand the quantitative values of parameters that are used in empirical functional like Skyrme DFT. It was shown that the LEC are strongly renormalized due to inmedium effects. This actually was also recently shown using BruecknerHartreeFock in Ref. [42] and was encoded in the ELYO functional through density dependent coupling constants in Ref. [43]. A review on the novel scientific activities in this field can be found in Ref. [44] (see also the recent work [45] for application to finite systems including pairing).
The thermodynamical properties of strongly interacting systems was studied in Ref. [46] using one of the functional proposed recently. While most of the observed properties of systems close to unitarity were reproduced very accurately, two difficulties have been identified. The first one is that the dynamical response function in the superfluid phase can without surprise only be achieved by introducing explicitly the pairing field in the functional. The explicit treatment of superfluidity is not the subject of the present work and we will concentrate on normal systems. A second source of difficulty is the absence of clear prescription for the effective mass in the large limit. Such effective mass and more generally quasiparticle properties are rather standard quantities helping to understand Fermi liquids. Its knowledge are of particular importance for instance to understand certain properties like the static response of neutron matter recently calculated with an abinitio theory in [47, 48]. It turns out, for instance for neutron matter, that the effective mass in neutron systems is scarcely known (see Fig. 6 of Ref. [46]) and has only been very recently estimated using AFDMC in Ref. [49] and BHF calculation [50, 51]. For this reason, we also explore the possibility to obtain selfenergies, for which direct contact with the Fermi liquid theory can be made, as functionals of in the nonperturbative regime. In order to achieve this goal, we also use resummation techniques and extend the phasespace average technique to the selfenergy. Finally, we briefly discuss how such analytical form can be useful in DFT approach.
2 DFT for dilute systems from manybody perturbation theory
We concentrate here on systems where the only low energy constant (LEC) at play is the swave scattering length . Infinite systems composed of spindegenerated particles of mass , i.e. a relevant situation for nonpolarized neutron matter and/or spin degenerated cold atoms, are investigated. Following Ref. [38] and using the EffectiveField Theory approach for homogeneous dilute Fermi gas, the swave interaction is simply written as a zerorange interaction that identifies with a constant in momentum space:
(1) 
where and are the relative momenta of the incoming and outgoing particles. The constant is linked to the scattering length through ():
(2) 
using the convention that a negative is attractive, so that the the swave scattering phase shift verifies . The model case, where the interaction is dominated by has been widely exploited in Fermi systems in the past [52, 53, 54, 55, 56] (see also [57]). The interaction (1) has a wellknown ultraviolet (UV) divergence. In the present sections, we summarize some known results for this model. Note that, the results have been obtained with proper treatment of the UV divergence using standard techniques (in particular minimal subtraction scheme of dimensional regularization). For more details see [58, 59, 60, 61, 62, 63, 64].
Hereafter, we use the twobody Hugenholtz energy diagrams convention for homogeneous dilute Fermi system based on particlehole propagator:
(3) 
The great simplification in infinite system stems from the fact that the relevant singleparticle states for the FreeGas (FG) are identified as planewaves with energy and occupation numbers at zero temperature given by where is the Fermi momentum of the system. The Feynman rules with the proper symmetry factor to estimate a contribution to the energy can be found in Ref. [38]. The diagrams contributing to the ground state energy up to third order are shown in table 1. Note that second order diagrams, i.e. composed by two vertex (black dot), and the third order diagrams are vanishing in this context.
We focus first here on the possibility to obtain the energy as a function of the quantity . The natural manybody approach if is to start from a perturbative expansion in powers of . We can then expand the energy as a series:
(4) 
where is linear in , is quadratic in , … denotes here the FreeGas (FG) energy defined as:
(5) 
where the spin degeneracy is equal to for the present case. Here we have use the fact that the Fermi momentum is linked to the density through:
(6) 
The perturbative approach leads to the socalled LeeYang expansion of the energy [52, 53, 54]. The same expansion was obtained in Ref. [38] using EFT technique. Keeping in mind that only connected diagrams contribute to the energy, the first order term in Eq. (4) is given by:
(7)  
where we recall the associated diagram. We recognize here the HartreeFock (HF) contribution at leading order in . The second order contribution stems from the direct and exchange terms due to the coupling between the 2 particle2 hole (2p2h) excited state and the uncorrelated HF ground state. It is given by:
Higher order contributions can be evaluated analytically (or numerically) in a similar way. For instance, the energy per particle at third order has been historically obtained by Efimov and Amusia [65, 66], Baker [67] and Bishop [68] and was more recently discussed in different works [40, 69, 70, 71, 72, 73, 74]. Very recently, the fourth order has also been worked out [75]. For now on, we will systematically assume that .
2.1 Link with DFT and shortcomings
Since our goal is to make connection between manybody approaches starting from the bare interaction and DFT, let us make simple early stage remarks relevant for the present discussion. Once the energy is written in terms of powers of , this energy can directly be interpreted as a functional of the local density using the equation (6). For infinite systems, the local density is just a constant. A standard technique to export a DFT in infinite system to finite systems is to use a Local Density Approximation (LDA) where the equation (4) is first transformed into an integral over space of the energy density functional:
(9) 
Here, contains the kinetic and the potential energy contributions that are both written in terms of the local density , itself linked to a local Fermi momentum . This approach, that has some connections with the ThomasFermi approximation [76, 77], leads to functionals of the local density consistently with the Hohenberg–Kohn theorem [78]. Such direct mapping has initiated several novel ways to obtain DFT in the nuclear physics context [24, 43, 45] (for a review see [44]). Similar strategy is now currently used in a different context in order to incorporate quantum corrections that might stabilize quantum droplets [79, 80].
The simple strategy discussed above to obtain a DFT for manybody systems has several shortcomings:

The perturbative approach provides a systematic and constructive approach to write the energy as a function of a polynomial of and/or of increasing orders. It however faces the difficulty that the number of diagrams to be evaluated and/or the complexity of the integrals that appear both significantly increase from one order of the perturbation to the next order.

Unless the expansion is made up to infinity in Eq. (4), the deduced energy only applies below a certain value of the density and/or . When the swave scattering length becomes large, the perturbative expansion at low density is not valid anymore. Typical examples in nature with large are unitary Fermi gas [81, 82] or nuclear neutron matter [57]. For instance, the swave scattering length for neutronneutron or protonproton interaction is , leading to a range of validity in density for the perturbative expansion. Compared to the saturation density , i.e. typical density in nuclear system, the perturbation theory is not appropriate to describe properly these systems at the relevant density scale.

A pure LDA approximation misses some important aspects that affect the system properties. An illustration is the effective mass that is standardly used in the Fermi Liquid Theory (FLT) approach. Consistently with the expansion of the energy given by Eq. (4), one might also obtain an expansion of the effective mass as:
(10) In the present case where only is considered, a correction to the bare mass starts to appear only at second order in perturbation, leading to the socalled Galitskii formula [83] (see section 4.1). Again, the expansion (10) truncated at a given order is usually restricted to the low density region.
The difficulties (i) and (ii) can eventually be solved by treating the problem numerically in a nonperturbative framework. With the increase of interest of systems with varying swave scattering length, several efficient numerical approaches have been developed. To quote some of them, we mention the Brueckner HartreeFock (BHF) [18, 57, 84, 85, 86, 87], the Selfconsistent Green function (SCGF) [16, 17, 88], the QuantumMonte Carlo (QMC) approach and/or AuxiliaryField Diffusion MonteCarlo (AFDMC) [9, 12, 13, 89, 90, 91], or the recently proposed approach based on Blod Diagrammatic Monte Carlo [19, 20, 21]. These direct numerical techniques, while very effective in some cases, do not lead to an energy written as a functional of . Their results can still serve as pseudodatas on which a DFT can be adjusted. This is particularly useful when properties of strongly correlated Fermi systems cannot be directly probed experimentally.
Here, however, we follow a different goal that is to have an analytical guidance to design a DFT for systems with large scattering length. Alternatively to a direct numerical method, the problem can be approximately solved by selecting certain classes of diagrams and by summing up these diagrams to all orders [39, 40, 57, 69]. As we will see below, this approximate treatment has also the advantage to automatically lead to analytical functions of . In the present work, our strategy is to use the resummation technique as a starting point to obtain a DFT for Fermi systems in the nonperturbative regime. Illustrations of the resummation technique applied to the energy are given in the next section. We then discuss, how a DFT can be deduced from it. As we will see below, considering directly the energy does not provides information on quasiparticle properties (item (iii)). For this reason, the methodology employed for the energy is extended to the selfenergy in section 4.
3 Diagrammatic resummation technique for the energy
Motivated by systems with anomalously large scattering lengths, several strategies have been proposed for selecting certain classes of diagrams and for providing compact expressions of the energy in terms of . In this section, we summarize some recent attempts to make resummation of diagrams using the interaction (1).
In 2000, Steele [69] laid the groundwork to describe the ground state properties of a Fermi system at low density with large scattering length based on resummation. He calculated explicitly the energy per particle up to fourth order in perturbation for ladder and ring diagrams. In this work, each contribution is written as an integral over the phasespace of a combination of the particleparticle, holehole, particlehole (resp. holeparticle) scattering loop functions. Noting that the contributions of energy diagrams composed by particleparticle loops are dominant when the order of perturbation increases, he proposed to retain only these contributions and sum them up to obtain a geometric series to be integrated in the accessible phasespace.
The resulting energy is schematically represented in Fig. 1a as an infinite sum of diagrams. Although rather tedious, the technique to obtain resummed expression is rather standard starting for instance from Greenfunction approach (see for instance [57]). After averaging over angles and using the minimal subtraction scheme of dimensional regularization [58, 59, 60, 61, 62, 63, 64], the energy can be written in the compact form (see for instance [39, 40, 69] and A for the definition of and ):
(11) 
In the seminal work of Ref. [40] based on the vacuummedium propagator consisting to rewrite the particlehole propagator (3) as:
(12) 
a larger class of diagrams were considered by the reorganization of the manybody diagrammatic calculation (see also [92] for further details). The diagrammatic representation of this resummation is shown in Fig. 1b. Contrary to the resummation of particleparticle ladder diagram at order by considering only particleparticle loops, this resummation takes into account diagrams composed by combined particleparticle and holehole loops. Paying particular attention to the combinatorial occurrence of each diagram, it was shown that the energy writes as a phasespace integration of an arctangent function [40]:
(13) 
In the following, results of the numerical integrations of Eqs. (11) and (13) will be referred as Geometric series Exact Integration (GEI) and Arctangent series Exact Integration (AEI) respectively. Both Eqs. (11) and (13) present several interesting features compared to a perturbative expansion. Firstly, they could be expanded in powers of and, noteworthy, their second order expansions match the LeeYang formula. Note that the Eq. (13) is slightly more predictive in the sense that it is valid up to third order in ^{1}^{1}1Only the ladder diagram contribution term at third order is accounted for since ring diagrams composed by particlehole loops (the nonvanishing diagram of table 1) are not taken into account in the resummation.. We show in Fig. 2a a comparison of the energy obtained by integrating numerically the two resummed expressions as a function of .
One of the motivations for the use of resummation is that, contrary to any truncation, the energy is not diverging as , i.e. in the unitary gas regime (see panel (b) of Fig. 2). This was firstly discussed in Ref. [39] for the GEI case and latter in [40] for the AEI case. As noted in these works, the ratio:
(14) 
generally refereed to as the Berstch parameter, significantly depends on the class of diagram selected for resummation. In the two cases considered here, we have:
(15) 
These values in both cases significantly differ from the experimentally observed value of the Bertsch parameter [13, 17, 93, 94]. It should be however kept in mind that the value of corresponds to the one of a superfluid unitary gas while superfluidity is not accounted for in the present resummation. Therefore, to be consistent, one should a priori compare with the value of the Bertsch parameter in normal systems. In [18], using Brueckner HartreeFock approach, a value was found, which is compatible with the experimental result of Ref. [94] giving . This value is actually consistent with the AEI case given by Eq. (13). Nevertheless, one obvious conclusion is that the choice of certain diagrams significantly affects the energy behavior as increases.
3.1 Phasespace average approximation for the resummed energy
Despite the fact that the selection of diagrams influences the results of a resummation approach, the resulting expressions of the energy in terms of is an interesting steps towards a DFT for interacting systems beyond the LeeYang formula. Still, the deduced expressions are rather complicated especially due to the necessity to perform explicit integrations on phasespace for all values of (note also that also appears in the definition of and , see A). This complexity can however be partially reduced using what we call below a PhaseSpace (PS) approximation.
The PS approximation was discussed for the Geometric series case in Ref. [39], it consists simply in replacing the numerator and denominator entering in the integral respectively by their average values integrated over the phasespace, leading in this way to a much simpler approximation. Let us introduce the notation defined as:
We see for instance that the GEI can be written as:
The PS approximation consists in replacing this expression simply assuming:
(17)  
This approximation still insures that the LeeYang expression is recovered up to second order in . Using the integrals given in A.1, we obtain:
(18) 
This compact form, called hereafter Geometric PhaseSpace (GPS), was introduced at several occasions in the nuclear physics and/or cold atom context [39, 75, 95, 96]. The energy obtained in the GPS approximation is shown in Fig. 3.
Eq. (18) could be interpreted as a minimal Padé approximation in at low density, the Padé recently shown for instance in Fig. 2 of Ref. [75]. We mention that Padé approximations Padé can be obtained (see for instance Ref. [95]) that could reproduce the development (LABEL:eq:devll) to a given desired order in for any .
One important conclusion for the present work is that the energy obtained from Eq. (18) largely extend the domain of density over which it reproduces the exact MonteCarlo result compared to the LeeYang formula, i.e. compared to the second order perturbation theory. Essentially, above in Fig. 3, the LeeYang expression deviates significantly from the exact calculation. Note that the inclusion of third order perturbation theory only slightly extend the domain of validity. On contrary, we see in Fig. 3a that the resummed formula follows closely the exact results and therefore it could be useful to obtain a compact form for a DFT beyond the perturbative regime. The GPS approximation has indeed been recently used in Ref. [24] to obtain a nuclear EDF where some of the parameters are directly connected to the physical swave scattering length, contrary to the widely used Skyrme EDF [2].
Although the main goal of the present work is to obtain DFT suitable beyond the perturbative regime, we would like to mention also that the approximated form (18) leads to a Bertsch parameter , that is closer to the one obtain at unitarity for superfluid systems [24, 39] compared to the one obtained with direct integration. It should be noted however that a relatively small difference in the value of leads to large deviations in the energy due to the fact that it is multiplied by the FreeGas energy. For this reason, it was proposed in Eq. (LABEL:eq:devll) to relax the slightly the low density constraint and adjust directly the denominator on the unitary gas in [8, 97]. Such strategy turns out to be be highly predictive for systems close and/or at unitarity [41, 46]. The functional proposed originally in [8] (assuming that the effective range cancels out, ) will be called Geometric series Unitary Limit (GUL) in the following. Note however, that here our primary goal is to see how far we can go beyond the perturbative regime without adjusting any parameter.
3.2 PhaseSpace approximation with Arctangent resummation [Ladder approximation]
Using the same approximation as above, a phasespace approximation of Eq. (13) can be obtained leading to the following compact expression, called hereafter Arctangent PhaseSpace (APS):
(19) 
Illustrations of the energy dependence obtained in the APS approximation are shown in Fig. 3. We note that the APS closely follow the GPS case at low density. This is indeed expected since both are constructed to match the same LeeYang expansion for low density Fermi gas. More surprisingly, the APS turns out to be very effective up to unitarity. It indeeds gives a Bertsch parameter equal to that is very close to . This is an interesting finding since, contrary to the GPS case where the unitary limit can only be reproduce at the price of degrading the description of the low density regime, in the APS case, both low density (second order expansion in ) and unitary limit can be very reasonably accounted for without adjusting any parameter.
For the sake of completeness and although that we do not expect to gain so much in terms of predictive power compared to the APS, we mention that similarly to the GUL case, an AUL (Arctangent Unitary Limit) approximation can be made. In that case, relaxing the constrain on the second order term in the LeeYang energy and imposing the value for the Bertsch parameter, we obtain:
(20) 
where .
The four approximations [GSP, APS, GUL, AUL] introduced here are rather simple functions of compared to the original GEI, AEI integral equations and therefore provide much simpler functionals of the density . One should mention a drawback of the phasespace approximation (see Fig. 3). By using phasespace average in the denominator, one restrict the value of that could be used. Indeed, while the integrated GEI and/or AEI integral forms can be applied from negative to positive values of around unitarity, this is not the case for the phasespace expressions where a pole appears for a certain positive value of . From now on, we will only consider the case where is negative that is also the relevant situation for neutron matter.
In summary, we have shown here that several functionals can be obtained that reproduces quite well the energy of Fermi gases at unitarity. We would like to mention that the value is the admitted value of superfluid unitary gas. It might then be surprising to reproduce this value with a functional originally motivated by the diagrammatic expansions of Ref. [40, 92] where superfluidity is not treated. It is however important to keep in mind that whatever is the motivation/strategy to produce a DFT, the only final criteria is the ability of the functional to accurately describe the groundstate energy of the system at various densities. This is actually the only purpose of a DFT constructed in the spirit of the original work of Hohenberg and Kohn who have shown the existence of a functional able to reproduce the exact energy consistently with the exact local onebody density.
In the following, however we would like to consider directly the selfenergy that is a priori clearly beyond the scope of a DFT approach. In this case, the discussion made above for the energy cannot be made and superfluidity should be explicitly introduced to describe superfluid systems. Such a treatment is beyond the scope of the present work and in the rest of the article, we will concentrate on nonsuperfluid systems.
4 Diagrammatic resummation for the selfenergy
Scientists working on the description of finite systems using DFT also often use these DFTs to get much more information than simply the energy of the system. Going from a Hohenberg–Kohn [78] to a Kohn–Sham [99] framework, a Slater determinant and/or more generally quasiparticle vacuum to treat superfluidity is first introduced to construct the normal and anomalous densities. In the nuclear physics context and although strong debates exist on the possibility to interpret physically the singleparticle and/or quasiparticle properties (see for instance [100]), the singleparticle shell evolution is standardly considered as a relevant output of the nuclear DFT [101, 3]. It is also standard to make connection with FermiLiquid theory [102]. As noted in Ref. [46], functionals that do not explicitly use the concept of effective masses and/or pairing gap will have difficulties to reproduce the excitation properties in strongly interacting Fermi systems. This is one of the main motivation of the present work. Starting from a manybody diagrammatic approach based on Green functions, a natural way to make connection with the Fermi liquid theory is to use the concept of selfenergy [57]. In the present section, we explore the possibility to start from a welldefined manybody framework based on selection of diagrams followed by resummation and see if the phasespace average approach used for the energy in previous section can be exported to the selfenergy. We will use as a starting point the resummation performed in Ref. [92] that is consistent with the GEI or AEI energy depending on the selected diagrams. In Ref. [92], “only” systems in the normal phase are considered. Similar semianalytical developments in the superfluid phase would require to use the GorkovGreen approach instead of the standard Green function framework. As far as we know, although this would be a very useful while most probably extremely involved problem, this has not been achieved yet. For this reason, we will only consider normal systems here.
4.1 Selfenergy for Fermi systems at low density
We restrict the present discussion to the onshell approximation of the selfenergy, so that the selfenergy , also called complexvalued singleparticle potential, is a function solely of the singleparticle momentum . Selfenergy are slightly more complex to obtain in the sense that more diagrams contributes to the selfenergy compared to the energy. This is illustrated in table 1 where we show the important diagrams in both cases for a contact interaction up to third order. By closing the legs, we recover the diagrams that have been estimated to obtain and previously.
Selfenergy  Energy  
Diagrams  Diagrams  
Using the perturbative approach, the proper selfenergy relevant for dilute systems can also be expanded as:
(22)  
The first two terms are wellknown and are respectively given by:
(23) 
Here, we have defined and we introduced the reduced momentum . The two functions and related respectively to the real and imaginary parts of the selfenergy were first derived by Galitskii [83] (Eqs. (34) and (35) of [83]). For the sake of completeness, they are also recalled in B. The third order contribution to the selfenergy (diagrams ( – ) in table 1) was studied within EFT in Ref. [70].
In the following, we will introduce the notation and for the real and imaginary parts of the selfenergy. These quantities enters respectively into the singleparticle (sp) energy and the lifetime of the quasiparticle (qp) respectively. We will essentially focus here on the real part of the selfenergy. The quasiparticle properties can be obtained from the behavior of the selfenergy close to () [103]. For instance, the chemical potential and the effective mass are respectively linked to the value of and its derivative at :
(24) 
Starting from the expansion (22) and making a Taylor expansion of each terms around (or ), one obtains a systematic approach to compute quasiparticle properties in powers of . For instance using the Taylor expansion:
(25) 
leads to the following expression of the chemical potential and effective mass up to second order in :
(26) 
These expressions are well known results also discussed in Ref. [83]. In particular, the latter equation is often refereed as the Galitskii mass formula.
4.2 Selfenergy resummation
The perturbative approach to the selfenergy is valid for weakly interacting systems in the low density regime. For this reason, as for the energy, some attempts to provide resummed expressions have been made in the past directly on the selfenergy. Our starting point is the resummation approach proposed in Ref. [92] for the selfenergy identified as the first functional derivative according to the occupation number of the GEI/AEI resummed energy so that the diagram selection rules are consistent with Ref. [40]. This consistency will be extremely useful when we will perform the phasespace average on the selfenergy.
Selecting the diagrams in table 1 that give those entering in the energy resummation of Fig. 1, a compact expression was obtained for the selfenergy in [92]. Such resummation can be schematically represented as in Fig. 4.
We consider here the two cases where either both particleparticle and holehole ladder diagrams are considered or only particleparticle ladder diagrams are used for the resummation. These two cases are respectively consistent with the AEI and GEI resummation for the energy. We refer to [92] for technical details (see also [104] for a complete discussion). After resummation and angle averaging approximation, the selfenergy can be written as (again with the convention ):
(27)  
Here and take different forms depending on the type of diagrams that are used. To avoid confusion between these two cases we will use the following convention: for the selfenergy obtained from the resummation of combined and ladder diagrams and for the one where only ladder diagrams are used. Similarly to the selfenergy that could be separated into a real and imaginary part, we can decompose and as:
The expressions of these functions are given in B.1. The resummed selfenergies of Ref. [92] have a number of interesting properties. First, in the low density regime, the first and second order selfenergies given by Eq (23) are properly recovered by construction. Another interesting feature of the expression given in B, is that the selfenergies do converge also to a finite result as for all . This is illustrated in figure 5 where the results of direct integration of Eq. (27) are shown as a function of at unitarity. Note that for the multidimensional numerical integration of the equations, we used the Vegas method implemented in the Cuba library [105]. The singleparticle energies defined as:
(28) 
and obtained in Fig. 5 coincides with those reported in Ref. [40] with marked bumps. These bumps seems unphysical not only because they significantly differ from the BHF calculations of Ref. [18] but also due to the presence of singleparticle energies above the Fermi energy for . We would like to mention that there is no reason that the deduced selfenergy is predictive at large scattering length due to the neglected diagrams.
Starting from the different expressions, one can deduce from it, Eq. (24), the quasiparticle properties by direct numerical integration. We show in Fig. 6 the evolution of and as a function of obtained in the two types of resummation considered here. We see that both AEI and GEI approximation significantly extend the domain of validity compared to the perturbative theory and the corresponding chemical potential extracted from them are both rather close to the BHF result up to . For the AEI case, it is remarkable to see, especially having in mind the strange behavior of Fig. 5, that the chemical potential extracted from the AEI case perfectly matches the BHF calculation for all regime of . We see however (panels (c) and (d) of Fig. 6) that the comparison is in general less favorable for the effective mass. Both approximations overpredict the effective mass compared to the BHF calculation for , even if the agreement is slightly better than the second or third order perturbation theory. We also observe that the effective mass, as well as the chemical potential, obtained with each other strongly depend on the selected diagrams.
A last important remark for the discussion below is that the expression of the resummed selfenergies obtained either from ladders or combined and ladders resummations are consistent with the GEI and AEI approximations given by Eqs. (11) and (13) respectively. Consistent means here that they respect the HugenholtzvanHove (HvH) theorem [107]. This theorem (at zero temperature) states that the singleparticle energy, given by Eq. (28), evaluated at the Fermi surface () is equal to the chemical potential of the systems. Using the thermodynamical relation between the chemical potential and the ground state energy, , where is the unit volume, the HvH theorem leads to:
(29) 
This equation gives a strong constraint between the energy and singleparticle potential at .
We have seen here that the resummation of diagrams is only a semisuccess to predict quasiparticle properties. More specifically, the effective mass deviates rather rapidly from the expected result as increases. Although their predictive power is limited, the AEI and GEI approximation can serve as a guidance to provide simplified expressions of the selfenergy that will be useful latter on the DFT context. We discuss below two approaches to obtain compact expression of the singleparticle potential.
4.3 Partial phasespace average for the selfenergy
As we have seen in section 3.1, the phasespace average, by avoiding the estimates of rather complex integrals, automatically led to simplified expressions for the resummed energy that turned out to be rather useful in practice [24, 8, 46]. The goal here is to develop an equivalent method directly at the selfenergy level.
The first difference compared to the energy is that the phasespace average should not be made on all variables because the selfenergy depends on (or ). In the following, starting from expression (27), for any function that depends on the variable , we will introduce the two averages:
(30)  
(31) 
that denotes partial phasespace (PPS) average at fixed value of respectively relevant for and . With these notations, Eq. (27) writes:
Expressions of the PPS for some functions are given in B.2.
Starting from these expressions, our goal is to provide for the selfenergy a phasespace approximation similar to the one for the energy given in section 3.1. Following the strategy we used previously, we will impose the approximate form to fulfill specific constraints:

Low density limit: We will always impose that the selfenergy matches the exact selfenergy in the low density limit up to a certain order in .

Large limit: We also seek for expressions that do not diverge in the limit .

Consistency with the HvH theorem: while it is not a priori absolutely necessary, in some cases, we will in addition impose that the selfenergy we obtained should be consistent with either the GPS or the APS energy. Again, consistency means here that the considered selfenergy and the energy obtained through phasespace average leads to the same chemical potential using Eq. (29). Note that the latter condition is more constraining than the condition (i) and (ii). In particular, since the energy already has the constraints (i) and (ii), they will be automatically fulfilled when (iii) is explicitly imposed.
In the following, we present two strategies to get approximate selfenergies using PPS. In the first strategy called hereafter simple PPS approximation, we only impose conditions (i) and (ii). While in the second strategy, that would be called consistent PPS, the form of the energy deduced will also be constrained to one of the PS approximation discussed in section 3.1 by requiring the HvH theorem to hold.
4.4 Simple PPS approximation for the selfenergy
We illustrate here a first simple strategy we can follow to impose (i) and (ii) avoiding complicated integrals. Our starting point is Eq. (27). As an illustration, we consider the specific case of the diagrams resummation and concentrate on the real part of the selfenergy with . Starting from the expressions given in B, a direct expansion of the denominator in gives:
(32) 
where the are given in terms of a PPS in the appendix by Equation (53). For instance while the expression of is given by Eq. (44). In the following, we will often use the second order approximation:
(33) 
as a reference for the low density limit. Let us assume that we seek for an approximate form of the selfenergy that matches the expansion (32) up to second order in while being nondivergent at large value of . Guided by the approximation made for the energy, one might simply use a Padé approximate form:
(34) 
Note that, guided by the APS expression, a similar expression can be obtained that approximate the AEI. From these expressions and using the expansion of given by eq. (25), we immediately see that the low density limit of the chemical potential and effective mass given by Eq. (26) are recovered. In addition starting from the expression of , the corresponding form of the energy in infinite systems can be obtained simply using the relationship:
(35) 
that can easily be obtained from the definition of as a partial derivative of the energy with respect to the particle number.
We compare in Fig. 7 the dependences of the chemical potential and effective mass obtained with the simple GPS approximation for the singleparticle potential respectively as a function of or . By construction, contrary to an approximation where the selfenergy is truncated at a given order in , simple GPS leads both to a smooth and converging behavior of these quantities up to infinite values of . These approximated selfenergies also reproduce correctly the low density limit. We see in particular that the chemical potential obtained with the simple PPS approximation follows closely the one of the original GEI result. We note also that the effective mass is more affected by the phasespace approximation. It turns out to be slightly lower compared to the GEI case and closer qualitatively to the BHF results up to . Not surprisingly, as in the original result obtained by direct integration, the large limit is not correctly accounted for. One difference we have observed however is that the dependence of given by Eq. (34) close to unitarity remains very smooth and does not presents the bumps seen in Fig. 5.
4.5 PPS with the HugenholtzvanHove theorem constraint: illustration with the GPS functional for .
The simple PPS approximation has some advantages. Among them, we note that the direct strategy used in previous section automatically leads the correct low density limit while the resulting singleparticle potential has a rather compact expression. We also saw that it gives quite reasonable behavior much beyond the perturbative regime. One drawback is that the corresponding energy obtained by direct integration through Eq. (35) turns out to be more complex than typically the GPS and/or APS functionals given by Eqs. (18) or (19). For this reason, we explored a different strategy that consists directly in imposing the constraint (iii), i.e. the energy obtained when applying Eq. (29) should match a preselected expression of the energy (the APS, GPS, GUL, AUL expressions for instance). In practice, this strategy is much less straightforward since the chemical potential is imposed whatever the value of and . It however has the direct advantage that all nice properties that were obtained at the energy level are automatically incorporated in the singleparticle potential.
Our starting point is to presuppose that we already know the expression of the energy in terms of . As an illustration, we consider below that the energy should match the GPS energy given by Eq. (18) obtained by the ladder approximation. From the imposed energy, and using Eq. (29), we obtained that the potential at should verify:
where we have used the expression of given by equation (25) directly recognized in the chemical potential. From this, we see that imposing the HvH theorem at (or ) gives us a strong guidance on the singleparticle potential for (or ). The simplest approximation that could be directly inferred from to obtain the potential consists in replacing by , i.e.:
(36) 
The present singleparticle potential presents several specific properties:

First, similarly to the simple PPS approach presented in previous section, its expansion to second order in matches the exact result for low density Fermi gas.

It also has automatically a nondivergent behavior in the limit . Due to the HvH theorem constraint, the associated limit is compatible with the value of the associated Berstch parameter, in the present illustration.

The form of the singleparticle potential turns out to be slightly more complicated than in the simple PPS approximation, Eq. (34). It is worth mentioning however that this form has strong similarities with the singleparticle potential obtained by ladder approximation (see Eq. (50)), in particular with the presence of two terms with similar dependence as in Eq. (36).

Last, we note that the HvH constraint solely does not uniquely define the form of the potential . Indeed, we can fulfill this constraint and keep all above mentioned properties using for instance the generalized formed:
(37) where the only constraint on is that . Again this flexibility should be seen as a positive point since it might be used to impose additional constraints latter on. In the following, we will generically denotes the potential given by Eq. (37) simply by and refer to it as the GPS approximation, while the case , leading to Eq. (36) is simply called GPS approximation for the singleparticle potential and is noted as . Unless specified, results presented below will be obtained in the GPS approximation.
Equivalent strategy can be followed starting from the APS approximation. Imposing the HvH theorem consistent with the APS approximation for the energy, we endup with the following expression for the singleparticle potential:
(38)  
This approximation is called APS hereafter. The only constraint on the function is now . If we assume that this function is constant for all , approximation called simply APS hereafter, we end up with the singleparticle potential:
(39)  
Again, all constraints (iiii) are respected and whatever the explicit form of , in the infinite scattering length limit, we obtain the APS value for the Bertsch parameter.
The chemical potential and effective mass dependence in obtained with the consistent GPS and APS scheme are displayed in Fig. 8 respectively as a function of or . The conclusion are similar as for the simple GPS case presented previously, i.e. the low density limit (respectively LeeYang and Galitskii formula) are properly reproduced by construction. The BHF results is reproduced qualitatively up to while the perturbative expansion breaks down around . We note however that the result of the consistent APS approximation is slightly worse compared to the original AEI as far as the chemical potential is concerned.
5 Quadratic and quartic approximation for the singleparticle potential
Our targeted goal here is to provide DFT inspired by the manybody resummation technique presented above. The clear advantage to start from the selfenergy level instead of the energy itself is that direct connections can be made between the selfenergy and the Fermi Liquid theory. This was illustrated previously with the chemical potential and the effective mass. Such quantities are also standardly obtained with Energy Density Functionals for instance used in the nuclear physics context, like Skyrme or Gogny EDFs. Empirical functionals, especially the functionals derived using Skyrme like contact interactions lead to very simple singleparticle potentials (see discussion below) with polynomial dependence in . For instance, the original parameterization proposed in [2] leads simply to quadratic dependence of the singleparticle potential in infinite matter. Novel generations of Skyrme EDF have been proposed leading to quartic or higherorder dependence in the momentum [110, 111, 112]. The justification that such simple approximation can contain important physical aspects can be found in [57]. The different singleparticle potentials obtained in previous section presents rather complex density and momentum dependence. However, starting from the PPS approximation, one might obtain a systematic polynomial expansion to a given order in . For this, we approximate the selfenergy obtained by assuming a polynomial form. We introduce the following expansion:
(40) 
This polynomial expansion, truncated at an appropriate order will not only be useful to make contact with empirical density functional theory but will also enable to obtain practical DFT for finite systems based on the present approach.
5.1 Quadratic approximation for the selfenergy
The simplest approximation that is also certainly the most highlighting one for the present discussion is to consider quadratic singleparticle potential, i.e. keeping only and in Eq. (40). Only two constraints are then needed to obtain and . One possibility is to impose that some of the quasiparticle properties are exactly recovered. For instance, imposing the chemical potential and the effective mass to be the same as the original ones obtained with one of the PPS approximation leads to:
giving finally:
(41) 
Since they are used as constraint, the present method automatically insures that the quasiparticle properties are preserved even if a simplified expansion is used for . Since the chemical potential is also constrained, due to the relationship (35), the energy of the system will also be preserved. Said differently, the GPS (resp. APS) approximation for the selfenergy combined with the polynomial approximation of will lead to the GPS (resp. APS) reference energy given by Eq. (18) (resp. (19)). This does not necessarily implies that the potential are similar, however by construction, they should become identical as becomes close to . We compare in Fig. 9 and 10 for the GPS and APS cases, the original GPS and APS potentials given by Eqs. (36) and (39) respectively, with their quadratic approximations for different values of . We clearly see in this figure that the value and the slope of at that are respectively linked to the chemical potential and the effective mass are identical. However some differences are observed as goes to zero. This also implies that some deviations occur with the secondorder expansion of the selfenergy given by Eq. (33) when the quadratic form of the potential is used in the low density limit. However, this approximation still leads to the proper behavior given by Eq. (26) in this limit.
5.2 Quartic approximation for the selfenergy
The agreement between the original potential and the polynomial expansion can be significantly improved simply by going to a quartic form of the potential, i.e. by truncating the expansion (40) to the next order. Then, the three parameters can be adjusted by adding to the chemical potential and effective mass constraint, the additional constraint that the value of at is identical to the one of the original GPS and APS approximation. We see in Fig. 9 and 10 respectively for the GPS and APS cases, that the quartic approximation significantly improves the global shape of the potential compared to the original PPS approximation.
6 Summary
In the present work, we explore the possibility to start from a welldefined manybody approach based on the diagrammatic resummation technique and obtain approximate expressions both for the energy and/or the selfenergy. While the present work is mainly focused on this topic, the ultimate goal is to use these approximate expression as a guidance for proposing new DFT for Fermi systems beyond the perturbative regime. At the selfenergy level, two simplifications are considered that might help in the DFT context. We first propose to use a partial phasespace approximation leading to simpler density dependence of the selfenergy and ultimately of the associated energy. We show that the partial phasespace average can be made either by imposing the form of the energy simultaneously or not. In the former case, the constraint on the energy is made using the HugenholtzvanHove theorem consistency. If the energy is not used as a constrained, the associated selfenergy expression turns out to be simpler. In all cases, the selfenergy of low density Fermi gas is properly recovered while a nondiverging limit is reached when contrary to truncated perturbative manybody framework. It is found in general that the combination of a resummation technique with phasespace average approximation can be used in a wider range of densities compared to manybody perturbation theory. We note however that without any adjustment, the functionals are not predictive in the unitary regime.
Besides the simplification introduced by the phasespace averages, we note in general that the resulting energy and/or selfenergy present both smoother behaviors that turn out to be qualitatively closer to the behavior calculated through Brueckner HartreeFock approach for non superfluid systems compared to a direct use of the resummation approach without PS approximation.
Guided by some phenomenological arguments commonly used for nuclei and also by the success of simple empirical functionals like the Skyrme DFT in nuclear systems, the singleparticle potential is further simplified by assuming that the singleparticle potential can be approximated by a quadratic or quartic polynomial in . This second approximation is made in such a way that quasiparticle properties are preserved for all . Again, this has the advantage that the low density limit of the chemical potential or the effective mass identifies automatically respectively with the LeeYang and Galitskii expression. We finally briefly discuss how the possibility to have analytical density dependent expressions of these quasiparticle properties can serve to design new DFTs. More generally, in view of the recent scientific emulation that followed the use of resummed formula for the energy [8, 24, 43, 44], the approximate expressions obtained in the present work can also be a strong guidance to obtain semiempirical or nonempirical functionals constrained at low density or at unitary or both.
As an illustration of this guidance, let us assume that we start from one of the parametrized forms of the chemical potential and effective mass obtained in the present work. We can then first rewrite the energy of the system in infinite system as:
Following Ref. [77], we can show that this energy can be obtained from a density functionals valid in both finite and infinite systems of the form:
(42) 
where the local density and the kinetic density entering in this equation are linked with the KohnSham states through and . One might notice in particular that in the present approach, corrections beyond the HF are automatically incorporated in both the potential and effective mass terms. This direction will be further explored in the near future.
The use of diagrammatic resummation leads to rather complex expressions in general. For this reason, we focused here the discussion on a rather simple case of a nonsuperfluid system with only one low energy constant, the swave scattering length focusing on the onshell selfenergy. In addition, an extension to include offshell effects would be a priori desirable especially to describe the Emass [113, 114]. Another natural extension of the present work is to include also the swave effective range and/or pwave scattering volume. The functional proposed in Ref. [8] and further discussed in [41, 46] already incorporate the effect of the effective range on the energy density functional. At leading order, one could simply add the present functional, however a proper treatment of the possible interference effects as well as effects beyond the leading order is needed if for instance becomes large. This case happens for instance in nuclear system at saturation density. We did not consider here the treatment of the effective range together with large scattering length, however we mention that the work of Refs. [115, 116] can be use as a starting point.
Another important extension would be to include the effect of pairing correlations. As we mentioned in the introduction, a prerequisite for the the present study is the possibility to obtain a compact expression of the energy after summing up selected diagrams to all orders in perturbations. By itself developing a perturbative approach on top of a quasiparticle vacua is possible (see for instance [117, 118]). However, even at second order, by replacing the 2p2h energy by the 4 quasiparticle, the complexity of the integrals to be made increases significantly. As far as we know, such problem has not be solved analytically. An alternative to the analytical approach using a DFT guided by the present work would be simply to add a posteriori a pairing energy to the DFT. This procedure is standardly used for nuclei and would at least allow to extend to the socalled Superfluid LDA (SLDA) approach of Ref. [7] away from unitarity.
Appendix A Useful definitions and integrals
In the main text, several quantities are written as integrals over the phasespace. For the sake of completeness, the different functions defined in the text as well as the different variable are summarized in this appendix. Note that a complete derivation of all equations can be found in Ref. [104]
Our starting point is the interaction matrix elements (1) written in momentum space as (note that, the spin is implicit and will lead to factors in the energy). Different functions appearing in the integrals for the energy and/or selfenergy are written as a function of and where and are vectors defined through:
(43) 
After proper treatment of the UV divergence and after averaging over vectors relative angles, it could be shown that the energy take the form (11) and (13), where , and are given respectively by [40, 92]:
a.1 Integrals used for phasespace average of the energy
The following integral are used to obtain the Phasespace averaged resummed expression for the energy: