Self-Consistent Green Function Method in Nuclear Matter

Symmetric nuclear matter is studied within the Brueckner-Hartree-Fock (BHF) approach and is extending to the self-consistent Green’s function (SCGF) approach. Both approximations are based on realistic nucleon-nucleon interaction; that is, CD-Bonn potential is chosen.The single-particle energy and the equation of state (EOS) are studied.The Fermi energy at the saturation point fulfills theHugenholtz-VanHove theorem. In comparison to the BHF approach, the binding energy is reduced and the EOS is stiffer. Both the SCGF and BHF approaches do not reproduce the correct saturation point. A simple contact interaction should be added to SCGF and BHF approaches to reproduce the empirical saturation point.


Introduction
The correct treatment of short-range correlations when performing nuclear matter calculations using the basic nucleonnucleon interaction (NN) is of great importance [1].One of these correlations has been studied using Brueckner-type resummation of ladder diagrams.This resummation allows to rewrite the ground-state energy of nuclear matter, using as an effective interaction the -matrix, which takes care of the short-range repulsive core in the nucleon-nucleon interaction [2,3].Calculations using realistic interactions lead to results, which lie along a line (the Coester line) shifted with respect to the phenomenological saturation point ( 0 ≈ 0.16 fm −3 , / ≈ −16 MeV) [4].The remaining discrepancy can be attributed to relativistic effects and three-body forces contributions [5].
Self-consistent approaches based on the in-medium matrix approximation for nuclear matter have been studied [6].In this way, a spectral function for nucleons in nuclear matter including two-particle correlations is obtained.The ladder diagrams involved in the calculation of the in-medium -matrix include also hole-hole (ℎℎ) propagation.The matrix approximation takes into account some of the higherorder hole line contributions as compared to the -matrix approach.It would be instructive to study the saturation properties of nuclear matter for the self-consistent -matrix approximation with realistic interactions.
Strictly speaking, however, the Bethe-Brandow-Petschek theorem [7] only defines the energy variable to be used in the calculation of self-energy or single-particle potential for the hole states.The choice for the propagator of the particle states is not defined on this level of the hole-line expansion and, therefore, has been discussed in a controversial way.The conventional choice has been to ignore self-energy contributions for the particle states completely and approximate the energies by the kinetic energy only.This conventional choice for the single-particle spectrum, however, is not very appealing as it leads to a gap at the Fermi surface; the propagator for single-particle states with momenta below the Fermi momentum   is described in terms of a bound singleparticle energy while the corresponding spectrum for the particle states starts at the kinetic energy for the momentum   .
Jeukenne et al. [8] argued that it would be more natural to choose the propagator according to the Green function method, that is, define the single-particle propagator with a single-particle energy which includes the real part of the self-energy as a single-particle potential for particle and hole states.This leads to a spectrum that is continuous at the Fermi momentum, which provided the name "continuous choice" for this approach.This continuous choice leads to an enhancement of correlation effects in the medium and tends to predict larger binding energies for nuclear matter than the conventional choice.
In the present work, we want to explore some features of the single-particle energy and the corresponding Green function for infinite symmetric nuclear matter at zero temperature.As a starting point we consider the Brueckner-Hartree-Foch (BHF) approximation employing a self-consistent continuous choice spectrum with an exact treatment of Pauli's operator [9,10].The contribution of hole-hole ladder terms are then added in a perturbative approach [11].Also, the binding energy of nuclear matter is studied.

The Model
2.1.BHF Approximation.Starting from realistic nucleonnucleon (NN) interaction, we have to use more advanced many-body approximations like the BHF which have the capability to account for the effects of correlations, which are due to the strong tensor and short-range components of such realistic NN interaction.The single-particle energy of a particle in the BHF approximation corresponds to the Hartree-Fock expression using the -matrix for the effective interaction.This means that the self-energy of a nucleon in nuclear matter with momentum  is given by (note that spin and isospin quantum numbers are suppressed) [12,13] with the occupation probability of a free Fermi gas with a Fermi momentum   : The matrix elements in (1) denote antisymmetrized matrix elements of the Brueckner -matrix that are determined by solving the Bethe-Goldstone equation for a given realistic NN interaction : The single-particle energies   of the intermediate states should be the corresponding BHF single-particle energies which are defined in terms of the real part of the BHF selfenergy of (1) by with a starting energy parameter Ω =  +   in the Bethe-Goldstone equation ( 3).
The Pauli operator ( 1 ,  2 ) restricts the intermediate states to particle states with momenta  1 ,  2 , which are above the corresponding Fermi momentum.However, the singleparticle spectrum is often parameterized in the form of an effective mass so that a so-called angle-averaged propagator can be defined, which reduces the Bethe-Goldstone equation to an integral equation in one dimension.Hassaneen et al. [4] have studied the properties of nuclear matter in the framework of BHF approach with an angle-average approximation to the Pauli operator using five different two-body NN potentials.It is found that BHF theory is able to produce saturation due to short-range correlations but fails to reproduce the empirical region of saturation quantitatively ( 0 = 0.16 fm −3 ; / = −16 MeV).The exact Pauli operator has been treated in [9,10].

Self-Consistent
Green's Function.One of the drawbacks of the BHF approximation is the fact that it does not provide results for the equation of state, which are consistent from the point of view of thermodynamics.As an example we mention that BHF results do not fulfill, for example, the Hugenholtz-Van Hove theorem [14].This is due to the fact that the BHF approximation does not consider the propagation of particle and hole states on equal footing.An extension of the BHF approximation, which obeys this symmetry, is the self-consistent Green's function (SCGF) method using the so-called -matrix approximation.During the last years techniques have been developed, which allow to evaluate the solution of the SCGF equations for microscopic NN interactions [15][16][17][18].Those calculations demonstrate that for the case of realistic NN interactions, the contribution of particleparticle ladders dominates the contribution of corresponding hole-hole propagation terms.This justifies the use of the BHF approximation and a procedure, which goes beyond BHF and accounts for hole-hole terms, in a perturbative way [11][12][13].
A perturbative treatment of hole-hole propagation avoids the appearance of the so-called pairing instabilities.A survey of pairing properties of nuclear matter media is presented in [19].This leads to a modification of the self-energy in the BHF approximation by adding a hole-hole term of the form [11] Using this expression of the self-energy induced by ℎℎ contributions, we are not following the rules of holeline or Bethe-Brueckner-Goldstone expansion.In fact, the Bethe-Brandow-Petschek theorem [7] does not apply to the self-energy of hole states if this correction is used.More information about SCGF theory is reviewed by Dickhoff and Barbieri [20].

Physics Research International 3
The quasiparticle energy for the extended self-energy can be defined as The spectral functions for hole and particle strength,  ℎ (, ) and   (, ), are obtained from the real and imaginary parts of the self-energy where the plus-minus sign on the left-hand side of this equation refers to the case of hole (ℎ,  <   ) and particle states (,  >   ), respectively.
In the case of two-body interactions, the hole spectral function gives access, through the Koltun sum rule [21], to the binding energy per particle: with deg denoting the degeneracy of the single-particle level, which is 4 for nuclear matter.The kinetic energy per nucleon is and the potential energy per particle has the form V/() = /() − /().Introducing the removal energy   () one can express the total binding energy per nucleon as In Brueckner-Hartree-Fock,   () =  2 /2 + (); that is, it coincides with the quasiparticle energy, and () denotes the step function, () = (  − ); therefore, Figure 1: The single-particle potential or real part of self-energy as a function of momentum.The calculations are done for symmetric nuclear matter with a Fermi momentum   = 1.333 fm −1 using CD-Bonn potential.The solid line represents the SCGF approximation, while the dashed one shows the BHF approach.

Results and Discussion
In Figures 1 and 2, we display the effect of hole-hole ladder contribution on the predicted self-energy or single-particle (s.p.) potential and single-particle energy or the quasiparticle energy ( qp ).Results for the BHF approximation of s.p. potential and s.p. energy are shown with dotted line whereas the solid line represents the SCGF one.The calculations are done at the empirical value for the saturation density of symmetric nuclear matter ( = 0.16 fm −3 ) using CD-Bonn potential [22].The corrections of the single-particle energies due to the Δ ∑ 2ℎ1 term are larger for momenta below the Fermi momentum and tend to zero for momenta above   .This can be seen already from Figure 1 and is shown explicitly in Figure 2, where the single-particle potential, that is, the quasiparticle energy minus the kinetic energy, is compared to the corresponding value obtained in the BHF approximation.
One finds that the single-particle potential derived from the quasiparticle energy tends to a constant for momenta below   .This means that the effective mass that describes the momentum dependence of the quasiparticle energy  qp is essentially equal to the bare mass for all momenta below   .At first sight, one may be tempted to consider this quasiparticle spectrum also in the energy denominators of the Bethe-Goldstone equation and the 2h1p correction term of (6).
From Figure 1, one can see also that the BHF self-energy does not have a simple parabolic shape as a function of the momentum.There is a characteristic dip in the selfenergy always occurring at momenta slightly above the Fermi momentum of the kind of nucleons under consideration.The momentum dependence of the single-particle potential in a homogeneous infinite system is a sign of the nonlocality of the single-particle potential and is typically characterized in terms of an effective mass m * , which can be used to parameterize the momentum dependence of the singleparticle energies according to (5).
In Figure 3, the binding energy per particle of symmetric nuclear matter is plotted versus the density.The results of the full SCGF calculation, computed at a temperature of  = 0 MeV, are indicated by solid line passing through solid circles in Figure 3.The BHF results are indicated by the triangles that are connected by a solid line.The empirical saturation point is given by the big square.The hole-hole contributions to the nucleon self-energy, which are included in the SCGF approach, yield a repulsive contribution the energy per nucleon, which increases with increasing nuclear density.This shifts the saturation point to a lower density and binding energy per nucleon so that also the saturation point obtained with SCGF is within the Coester band.Due to this effect, the saturation density in the SCGF approximation,  sat = 0.211 fm −3 and  sat = −11.34MeV, is considerably smaller than the saturation density in the BHF approach that has the value  sat = 0.391 fm −3 and  sat = −23.97MeV.After all the saturation density is reduced by almost 50%; however, both values are not near the empirical one ( sat ≈ 0.16 fm −3 ,  sat ≈ −16 MeV).
A nice property of the SCGF approach is its thermodynamic consistency.In particular, it fulfills the Hugenholtz-Van Hove theorem, which states that whenever a further nucleon is added to the nuclear system in equilibrium, that is, at the saturation point, the energy gain is exactly the binding energy of a nucleon [14], the chemical potential  =  sat .In Figure 3, the chemical potentials are given by the dotted lines.The violation of the theorem in the zero temperature continuous choice BHF calculation amounts to about 20 MeV.
In order to obtain the experimental saturation value in a symmetric nuclear matter, one has to add threebody interaction terms or a density-dependent two-nucleon  interaction.Therefore, it is quite natural to supplement the effective interaction or the self-energy of BHF and SCGF calculations by a simple contact interaction, which we have chosen following the notation of the Skyrme interaction to be of the form where  is the matter density and  0 ,  3 , and  are parameters.For a fixed value of  (typically  = 0.5), we have fitted  0 and  3 in such a way that BHF and SCGF calculations plus the contact term (ct) of ( 14) yield the empirical saturation point for symmetric nuclear matter; that is,  corr / = ct.The results for these fitting parameters are listed in Table 1, and the corresponding energy versus density curves are displayed in Figure 4.For all cases, the fit yields an attractive two-body contact interaction and a repulsive  3 term.The results are compared to those obtained by Baldo and Shaban [23] and supplemented by a microscopic three-body force and relativistic BHF.One can see from Figure 4, that the present calculations plus the contact term give stiffer EOS like relativistic BHF especially at high density [24].

Conclusion
The single-particle energy is calculated self-consistently using BHF and SCGF approximations.The contribution of the ℎℎ terms leads to a repulsive contribution to the single-particle energy which decreases with momentum.This implies that Figure 4: / in MeV for symmetric nuclear matter as a function of density using CD-Bonn potential.The results are obtained using both BHF and SCGF approaches with and without an additional contact interaction of the form displayed in (14), in comparison with microscopic three-body force by Baldo and Shaban [23] and relativistic BHF [24] by Müther et al.
the equation of state derived within the SCGF approach tends to be stiffer than the corresponding one evaluated within the BHF approximation.Also, the inclusion of ℎℎ propagation in the SCGF approach leads to a somewhat stiffer chemical potential than in BHF approximation.The Hugenholtz-Van Hove theorem is well fulfilled in the SCGF approximation while it is badly violated in the BHF approach.Therefore, the present approximation should avoid the so-called pairing instability which occurs in the SCGF approach using the Galitskii-Feynman propagator.A more systematic study of symmetric nuclear should be useful, which goes beyond the present model and treats the propagation of hole-hole scattering terms in a nonperturbative way [15][16][17][18]25].
Both of the present methods fail to reproduce the empirical saturation point.Thus, these microscopic approaches are supplemented by a density-dependent contact term to obtain the empirical saturation point of symmetric nuclear matter.

Figure 2 :
Figure2: The single-particle energy determined from the SCGF approximation, solid line, and BHF approach, dashed line.The results are obtained for symmetric nuclear matter at a Fermi momentum   = 1.333 fm −1 using CD-Bonn potential.

Figure 3 :
Figure 3: The binding energy per nucleon / (black curves) and the chemical potential  (red curves) as a function of the density .The box indicates the empirical saturation point.See text for details.

Table 1 :
The values of the parameters  0 and  3 that are used to fit the saturation point.