The LST employed in Ref. [3] was based on
HO densities corrected in the
asymptotic region by the contribution from the lowest-energy
quasiparticle. Since a common LST has to be carried out
for both
neutrons and protons, for each nucleus one is forced to make a
decision whether the LST is to be based on neutron or proton
density. In Ref. [3] we used a prescription (referred to as LAM)
that the
neutron densities were used for and
vice versa. In this work, we use the condition
, where
is the point where the neutron or proton logarithmic
density has a minimum as a function of
. In practice, the above
condition, dubbed RHO, does not depend on whether the neutron or proton
is considered.
![]() ![]() |
In Fig. 1 (left panel) we show the differences in
obtained in HFBTHO
by using the LST condition employing the Fermi energies
(LAM) [3] or the densities (RHO). One can see
that in the majority of neutron-rich nuclei both prescriptions lead
to identical results. However, in many proton-rich nuclei the new prescription
decreases binding up to about 500keV, and for some
medium-mass proton-rich nuclei the RHO method decreases
binding by up to 100keV. This latter
effect is due to a better description of asymptotics in the pairing
channel, which leads to extended pairing fields and reduced
pairing energies [10].
The right panel of Fig. 1 shows differences in
obtained in THO and HO bases. In most nuclei, by
using the THO basis, one obtains a small energy gain of up to 10keV.
This grows to
500keV for the very neutron-rich systems.
Again, in lighter nuclei, a better asymptotics may lead to a reduced
binding. In fact, our results show that improvements in density
profiles at large distances cannot be treated variationally.
First,
is quite insensitive to the precise description of nucleonic
densities in outer nuclear regions. Second, due to the pairing-space cutoff,
the pairing energy is not reacting variationally
on the improvement of the wave function.
Figures 2 and 3 present HFBTHO results obtained with the SLy4 and SkP Skyrme forces. It is obvious that without further improvements these traditional Skyrme forces describe nuclear masses rather poorly. The rms deviations between calculated and measured masses are as large as 3.14MeV for SkP and 5.10MeV for SLy4, respectively, as compared to about 0.70MeV deviations obtained for forces fitted specifically to masses (see Ref. [11] for a review). Moreover, pronounced kinks obtained at magic numbers suggest that the quality of the description of (semi)magic and open-shell systems is not the same. This may point to a need to systematically include dynamical zero-point corrections [12]. Work in this direction is in progress.
![]() ![]() |
![]() ![]() |
This work was supported in part by the Polish Committee for Scientific Research (KBN); by the Foundation for Polish Science (FNP); by the U.S. Department of Energy under Contract Nos. DE-FG02-96ER40963 (University of Tennessee), DE-AC05-00OR22725 with UT-Battelle, LLC (Oak Ridge National Laboratory), and DE-FG05-87ER40361 (Joint Institute for Heavy Ion Research); and by the National Nuclear Security Administration under the Stewardship Science Academic Alliances program through DOE Research Grant DE-FG03-03NA00083.