In the Skyrme-HFB approximation, the total energy
of a
nucleus is given as a sum of kinetic, Skyrme, pairing and Coulomb
terms:

The derivation of the energy is explained in detail in several articles (see Refs. [3,11] for the HF energy and [5] for the HFB case). The expression of the Skyrme force is given in section 3.

In equation (17), the kinetic energy of both neutrons and protons is given by isoscalar kinetic density. The neutron and proton masses being approximated by their average value, and one has

(17) |

where the factor in parentheses takes into account the direct part of the center-of-mass correction [12]. Since we only consider even-even nuclei, the Skyrme part of the energy functional is time even, it can be written as a sum of isoscalar () and isovector () parts (see

Index stands for neutrons and protons while the absence of index indicates the total (isoscalar) density. Using the same notations for the densities, the pairing energy density reads

(19) |

In this last expression we have added a prime to the parameters of the interaction. Although the derivation of the general Skyrme energy density functional is based on a unique force, its

The energy density defined in (17) involves a Coulomb term for protons. This term contains a direct part which can be expressed using the charge density and leads to a local field after variation of the energy and an exchange part which would lead to a non local potential. Both parts deserve a special discussion. The direct part of the Coulomb energy depends on the charge density and reads

We use the point proton density to simplify this expression,

The error introduced by this approximation has been estimated by studying the next order term [15] or more recently by comparison with the exact treatment of the Coulomb energy [16]. The Slater approximation seems to have little consequence on the nuclei bulk properties although some significant effects can be expected on the position of the proton drip-line and the Coulomb displacement energy of single particle levels. Finally, the Coulomb contributions to the pairing energy and fields are not included in the