kl800.com省心范文网

Hartree Fock Calculations in the Density Matrix Expansion Approach

F. Hofmann and H. Lenske

Institut f¨ ur Theoretische Physik, Universit¨ at Gie?en D-35392 Gie?en The density matrix expansion is used to derive a local energy density functional for ?nite range interactions with a realistic meson exchange structure. Exchange contributions are treated in a local momentum approximation. A generalized Slater approximation is used for the density matrix where an e?ective local Fermi momentum is chosen such that the next to leading order o?-diagonal term is canceled. Hartree–Fock equations are derived incorporating the momentum structure of the underlying ?nite range interaction. For applications a density dependent e?ective interaction is determined from a G-matrix which is renormalized such that the saturation properties of symmetric nuclear matter are reproduced. Intending applications to systems far o? stability special attention is paid to the low density regime and asymmetric nuclear matter. Results are compared to predictions obtained from Skyrme interactions. The ground state properties of stable nuclei are well reproduced without further adjustments of parameters. The potential of the approach is further exempli?ed in calculations for A=100. . .140 tin isotopes. Rather extended neutron skins are found beyond 130 Sn corresponding to solid layers of neutron matter surrounding a core of normal composition. 21.60-n,21.60Jz,21.10.Dr

arXiv:nucl-th/9705049v1 28 May 1997

I. INTRODUCTION

Studies of nuclei far o? stability have given a new impact to nuclear structure physics. Exotic nuclei allow for the ?rst time to test nuclear theory at extreme isospin and, close to the driplines, also at low density. In the past nuclear models have been derived almost exclusively from well bound and charge symmetric stable nuclei. Hence, structure investigations far o? stability are providing important extensions into hitherto unexplored regions of the nuclear landscape. It is apparent that such studies involve far reaching extrapolations of nuclear models much beyond the range where they were initially developed. The weaker binding of strongly asymmetric systems enhances interactions between bound and unbound con?gurations and the regime of open quantum systems might be approached. In extreme dripline nuclei as e.g. 8 B [1] and 11 Li [2] evidence for a strong competition of interactions between individual nucleons and mean–?eld dynamics is found. At present nuclear structure theory is making ?rst exploratory steps into the new regions far o? β -stability. Experience shows that seemingly well established models describing stable nuclei equally well are leading to strong deviations when applied to exotic systems. A particular striking feature are the variations in predictions of binding energies and other ground state properties at large charge asymmetry. The results indicate that the isovector properties of interactions are not well enough understood in order to reproduce or predict the location of driplines. A widely applied method in non–relativistic structure calculations is the Skyrme approach [3]. Although very successful in stable nuclei, it is an open question to what extent the restricted set of operators used in the Skyrme energy functional is su?ciently accounting for the conditions in the new mass regions. The transparent structure and the good applicability is a clear advantage of the Skyrme model. It is also ?exible enough for adjusting the already incorporated parameters or introducing additional ones in order to describe new data. However, a purely empirical approach might be hampered by increasing uncertainties when approaching the limits of stability. The intention of the present work is to indicate an alternative approach which is closer to a more microscopic description but still retains to a large extent the advantages of using a local energy functional as in Skyrme theory. A more fundamental description requires to use microscopically derived in–medium interactions as e.g. the Brueckner G–matrix accounting for medium e?ects on the level of ladder diagrams [4,5]. Theoretically, such a full scale many-body calculation is not feasible as a standard approach to ?nite nuclei, except for selected cases [6,7]. A successful program is to use instead a G–matrix from in?nite nuclear matter and to apply it in the local density 1

approximation (LDA) to ?nite nuclei [4]. The LDA has been used succesfully in non–relativistic [8–10] and also covariant ?eldtheoretical models [11–18]. However, the ?nite range character of a realistic G–matrix is well known to introduce complications on the numerical level by leading to a system of integro–di?erential equations [19]. A way around this problem was discussed some time ago by Negele [8] but rarely applied later. The numerical complications arise because of the non–locality of exchange contributions to the HF–?eld. For a ?nite range interaction they involve the full one–body density matrix rather than the local densities needed for the direct parts and contact interactions of Skyrme type. Negele introduced an operator-valued expansion of the non– local parts of the density matrix. The contraction with the HF–ground state con?guration leads to a systematic expansion in terms of the local one-body density. The non–localities are contained in functions where the lowest order contribution coincides with the well known Slater approximation [8,20]. This density matrix expansion (DME) allows to integrate out the non–localities of the exchange terms. The approach presented below follows in spirit closely the original work of Ref. [8,9]. In section II the basic steps of the DME are summarized. Di?erent to Negele we choose an e?ective local ”Fermi” momentum q2 F (r) such that the next to leading order term of the DME series is canceled exactly. q2 is determined by the kinetic energy density and, in addition, the Laplacian of the one–body F density. For an in?nite system the conventional expression for the Fermi momentum is recovered. The DME results are then used to de?ne a local energy density functional. Since it originates from an e?ective interaction, assumed to be given independently e.g. by Brueckner calculation, the parameters of the model are ?xed. By variation of the energy functional HF–equations are derived. With our choice for qF the kinetic energy part of the single particle Schroedinger equation receives additional contributions. Thus, we obtain consistently e?ective masses m? p,n (r) for protons and neutrons, respectively, which are completely determined by the momentum and isospin structure of the original interaction. An important extension over the Skyrme approach is that the self– consistent local momentum qF contributes in arbitrarily high order to the exchange terms which is especially important at high and low densities. Applications to in?nite symmetric and asymmetric nuclear matter are discussed in section III. We start by choosing a somewhat simpli?ed e?ective interaction. The three Yukawa parameterization of a G–matrix calculated initially for 16 O by Toki et al. [21,22] is used. This so–called M3Y interaction has been widely and successfully applied in nuclear reaction [23,24] and structure calculations [10,25,26]. The former studies, especially for nuclear reactions, indicate that this interaction has a realistic momentum structure including a long range one pion–exchange tail. It is especially trustable at low nuclear densities around 1/3 of the saturation density. A density dependence is lacking and it is known that the original M3Y parameterization does not lead to saturation of nuclear matter. As other authors before [24], we restore the density dependence by ?tting the properties of symmetric nuclear matter at saturation including the binding energy, equilibrium density, compressibility and, in addition, also the symmetry energy. This allows to ?x the density dependence of the isoscalar and isovector vertices in the spin–scalar particle–hole interaction channels. We thus have determined a renormalized density dependent three Yukawa (D3Y) interaction which we expect to resemble closely the momentum structure of a G–matrix but with an e?ective strength ?tted to nuclear matter properties [27]. The in–medium particle–particle interaction is obtained by the appropriate transformation and used in pairing calculations. At low densities the strength of the free singlet–even NN T–matrix is approached. Results for symmetric and asymmetric nuclear matter are compared to calculations with empirical Skyrme interactions. Close to saturation a good agreement is found but deviations occur away from that point. In section IV the model is applied to ?nite nuclei. The parameters are kept ?xed as obtained from the nuclear matter ?ts. Binding energies, charge densities and radii of stable nuclei (16 O,40,48 Ca,90 Zr and 208 Pb) are well described. The good agreement con?rms the approach and makes applications to unstable nuclei meaningful. HF calculations for Sn isotopes lead to a good description of the binding energies for A≥120 but a slightly too weak binding is found in the lighter nuclides. The calculations predict rather thick neutron skins for A≥130. The paper closes in section V with a summary and conclusions.

2

II. DENSITY–MATRIX EXPANSION AND HARTREE–FOCK THEORY

The total Hartree–Fock energy of a nucleus with A nucleons and the single particle states |k is given by the kinetic energy T =

k ≤A

k|

p2 |k = 2m

d3 r

?2 h τ (r) = 2m

d3 r

?2 h {τp (r) + τn (r)} , 2m

(1)

where q = n, p denotes protons and neutrons and τq (r) the kinetic energy densities. The potential energy of a two-body interaction V ? = V

k1 k2

k1 k2 |V |k1 k2 ? k1 k2 |V |k2 k1

(2)

is expressed in terms of interactions for like (q = q ′ ) and unlike (q = q ′ ) particles as de?ned in App. A ? = Vd + Ve = V

qq′ =p,n d e Vqq ′ + Vqq ′

(3)

=

qq′

d3 r1

e d d3 r2 ρq (r1 )ρq′ (r2 )Vqq ′ (r12 ) + ρq (r1 , r2 )ρq ′ (r1 , r2 )Vqq ′ (r12 )

where direct and exchange contributions are indicated separately by the index d resp. e. Here, r12 = r1 ? r2 is the relative coordinate and the density matrices are de?ned as ρq (r1 , r2 ) =

k σ

φ? k (r1 , σ, q )φk (r2 , σ, q )

(4)

ρq (r) = ρq (r, r) where φk (r2 , σ, q ) denotes single particle wave functions.

A. Density Matrix Expansion

The exact treatment of the exchange term in the Hartree-Fock functional would, after variation, lead to a coupled system of integro-di?erential equations. A numerically simpli?ed and theoretically more transparent approach is obtained using the density-matrix expansion (DME) which was originally invented by Negele and Vautherin [8]. In the DME approach one ?rst transforms to r2 and s = r2 ? r1 . Then the density center-of-mass and relative coordinates by introducing r = r1 + 2 matrix is formally expanded into a Taylor series in s. A particular simple result is obtained if the time-reversed orbitals in the nucleus are ?lled pairwise [8]. Therefore, the DME approach is most suitable for even-even nuclei. As shown in [8], averaging over spin directions and resummation leads to a series expansion where the ?rst two terms are s s ρ r + ,r? = ρSL (sqF ) ρ(r) 2 2 35 3 2 1 2 + 3 j3 (sqF ) 4 ? ρ(r) ? τ (r) + 5 qF ρ(r) + · · · 2sqF (5)

and in the leading order term the well known Slater approximation [8,20] is recovered. Here, the 3 j1 (sqF ) is given by the spherical Bessel function jn (x) of the order Slater density ρSL (sqF ) = sq F n = 1. The quality of the expansion depends on the choice of q2 F (r) which should obviously correspond to the average relative momentum between the two interacting particles. Di?erent to Ref. [8] we choose

2 qF (r) = 2 5 τ (r) ? 1 4 ? ρ(r) . 3 ρ(r)

(6)

3

which is seen to cancel exactly the next to leading order term in Eq.(5) and has the property 2 d3 r ρ(r)qF (r) = 5/3 d3 r τ (r). s s ,r? 2 = With this choice the DME reduces to the well known Slater approximation ρ r + 2 ρSL (sqF ) ρ(r) of the exchange term but now with a modi?ed Fermi momentum that accounts for surface corrections to the Fermi momentum of the Thomas Fermi model up to the second order. The advantage of this method is that the description of the exchange potential is simpli?ed and can be treated in a way similar to the nuclear matter calculation. It is easily seen that our choice of the Fermi momentum reduces in in?nite nuclear matter to the conventional expression kF = (3/2π 2 ρ)1/3 . Once the relative and the center-of-mass coordinates are separated and the integration over the relative coordinate s can be immediately performed analytically. The potential energy of the exchange part ?nally simpli?es to: Ve =

qq′ e ?qq d3 r ρq (r)ρq′ (r)V ′ (r)

(7)

In this equation the strength of the interactions are de?ned as ? e ′ (qFq , qF ′ ; r) = V qq q

e d3 sVqq ′ (s)ρSL (sqFq (r))ρSL (sqF ′ (r)) q

(8)

In this way the exchange term of a ?nite range interaction is mapped to an e?ective contact interaction with a local strength depending on the e?ective local Fermi momentum, Eq.(6), at the place of the interaction

e 3 ?e Vqq ′ (s, r) → δ (s) Vqq ′ (qFq (r), qF ′ (r)). q

(9)

We are now at the position to reformulate the theory in terms of a local energy density functional H (r) = HT (r) + HV d (r) + HV e (r) with the total energy E= d3 r H (r) = T + 1 ? V . 2 (11) (10)

Comparing this to the equations derived before one ?nds: ?2 h [τp (r) + τn (r)] 2m 1 d HV d (r) = ρq (r)Uqq ′ (r) 2 ′ HT (r) =

qq

(12) (13) (14)

HV e (r) =

1 2

e ρq (r)Uqq ′ (r) qq′

Single particle potentials U (r, ρ) have been introduced as e.g.

d Uqq ′ (r) = d ′ d3 r′ ρq′ (r′ ) Vqq ′ (r, r )

(15) (16)

e ?e 2 2 Uqq ′ (r) = ρq ′ (r)Vqq ′ (qF , qF ′ ; r) q q

B. The Hartree–Fock Equations

In our case, the total energy depends in a complex, nonlinear way on the nuclear density ρ, the kinetic energy density τ and ?2 ρ. By using the DME we are able to simplify the exchange term. Our choice of the Fermi momentum qF (ρ, τ, ?2 ρ), Eq.(6), leads to an additional functional 4

dependence on the kinetic energy and the Laplacian of the density and therefore takes into account dynamical and surface corrections. For the variation of the HF energy functional, Eq.(10), we 2 express ρ, τ and qF by φk and its derivatives: ρ(r) =

k

|φk (r)|

2

(17) (18)

2

τ (r) =

k 2 qF (φk , ?φk , ?2 φk ) =

|?φk (r)|2 51 6ρ |?φk | ?

k

1 ? 2 φ ? φk + φk ?2 φ? k 2 k

(19)

The HF variational equation d3 r δ leads after partial integration to δE =

k

H (r) ?

k

|φk (r)| ?k

2

=0

(20)

d3 r δφk

?H ?H 2 ?H ?? ? ? ? ? ?2 φ? ? φk ?k ?φ? ? ? φ k k k

(21)

For the further calculations we assume that the nucleon–nucleon interaction has an intrinsic density dependence resulting from the medium dependence of the meson–nucleon vertex that can be parameterized by vertex corrections gγ (ρ). A model will be speci?ed in the next section. Then, the partial derivative in Eq.(21) with respect to φk is written as

2 ?H ?H ?ρ ?H ?gγ (ρ) ?ρ ?H ?qF = + + 2 ? ? ? ?φk ?ρ ?φk ?gγ (ρ) ?ρ ?φk ?qF ?φ? k

(22)

The ?rst term is the ordinary HF potential, the other terms are rearrangement terms resulting from the intrinsic density dependence of the interaction. Analogously the other partial derivatives are expressed as

2 ?H ?τ ?H ?qF ?H = + 2 ? ?φ? ? ?φ? ?τ ? ?φ? ?qF k k k 2 ?H ?H ?qF = 2 ? ?2 φ? ? ?2 φ? ?qF k k

(23) (24)

Here, the second term in Eq.(23) gives rise to an e?ective mass. Using the above de?nitions and plugging in everything into Eq.(21) leads to the Schroedinger-equation ?? ?2 h HF R ? + Uq (r) + Uq (r) φkq (r) = ?kq φkq (r). 2 m? q (r) (25)

The Hartree-Fock potential is given by

HF Uq (r) = q′ d e Uqq ′ (r) + Uqq ′ (r)

(26)

and the e?ective masses are m? q (r) = 1+ m

2m 5 h ?2 6

e ?Uqq 2 ?qF q

+ 2 ?qqq 2

?U e

.

′

(27)

Fq

The e?ective mass results from the dependence of the Fermi momentum on the kinetic energy density τ in the exchange term. It is a direct consequence of our choice of the Fermi momentum and the ?nite range nature of the interaction. 5

The variation of the density dependent vertex functions gγ (ρ) and of the Fermi momentum in the Fock energy density leads to an additional rearrangement potential

R Uq (r) =

? (HV d + HV e ) ?gγ (ρ) ?gγ (ρ) ?ρ

e e ?Uqq ′ 1 2 ?Uqq +2 2 ? qF 2 q 2 ?qFq ?qFq

(28) ?

e e ?Uqq ′ 5 2 ?Uqq +2 2 ? 2 24 ?qFq ?qFq

which is of importance for the mean–?eld or single particle potential in the HF equations. It results from the density dependence of the NN interaction and is known to describe (static) polarization of the background medium [18,28,29].

III. NUCLEAR MATTER A. The Nucleon–Nucleon Interaction

A well established approach to derive in–medium interactions from the NN–interaction in free space is Brueckner theory. However, the ladder approximation, even if re?ned by going beyond the quasiparticle approximation as e.g. in Ref. [30], is generally found to miss the empirical saturation properties of nuclear matter at least when including two–body correlations only. For the purpose of the present work we therefore use a semi-phenomenological approach. We start from the M3Y G–matrix [21,22] obtained originally for 16 O. The M3Y G-matrix is parameterized by a superposition of three Yukawa functions where the ranges were chosen to represent the long-range tail of the one-pion-exchange potential (OPEP) and medium and shortrange parts. The shorter ranges correspond to boson masses m = 490 MeV and m = 780 MeV and thus mimic σ , ω and ρ meson exchange interactions. Here, we use the parameter set obtained from the Paris NN-potential [21] which is given in Tab. I. Because of the lack of a density dependence the M3Y interaction can be considered as realistic only for a small range of nuclear densities [31] being centered around one third of the saturation density. In connection with optical potential and nuclear reaction studies various authors, e.g. Kobos et al. [32] and more recently Khoa and von Oertzen [24], have noted that the ”bare” M3Y interaction leads at high densities to the collapse of nuclear matter but by supplementing a phenomenological density dependence the correct binding properties could be restored. In spirit we apply the same approach by introducing density dependent vertices ?tted to the saturation properties of in?nite symmetric nuclear matter. In the former works only the medium dependence of the isoscalar–spin scalar channel was considered [24,32]. Here, we choose a more general description and include a density dependence also in the isovector channel. We introduce density dependent isoscalar and isovector vertex functions f0,τ (ρ), respectively, multiplying the corresponding ”bare” M3Y vertices as a whole but leaving the intrinsic momentum structure untouched. Hence, we assume that the many–body and the momentum structures are separable in which we follow the observations made in determining e?ective in–medium meson–nucleon coupling constants from Brueckner calculations [14–16,33]. Only central interactions are considered. Aiming at applications to asymmetric systems it is suitable to rewrite the central singlet and triplet even and odd components (VSE , VT E , VSO , VT O ) in terms of interactions for like (proton-proton, neutron–neutron) and unlike (neutron–proton) particles. We then consider the interaction in the ”particle–hole” channel which eventually determines the Hartree–Fock mean–?eld. Multiplying the spin and isospin exchange operators onto the interaction e?ective direct (d) and exchange (e) interactions are obtained [34].With the relations given in App. A the isoscalar and isovector interactions are then obtained as V0

d/e

Vτd/e

1 d/e d/e d/e d/e (V + Vpn + Vnp + Vnn ) 4 pp 1 d/e d/e d/e d/e = (Vpp ? Vpn ? Vnp + Vnn ) 4 =

(29) (30)

and medium dependent vertex renormalizations f0,τ (ρ) are introduced 6

V0,τ (ρ) = f0,τ (ρ)V0,τ

?d/e

d/e

(31)

taken to be the same in the direct and the exchange channel. Back transformation to the singlet/triplet representation (see App. A) is used to de?ne the corresponding in–medium particle– particle interaction. By this procedure we retain the basic relations of particle-hole and particleparticle interactions. The question arises how to choose the functional form of the vertices. As a guidance we refer to Brueckner theory and express the G–matrix R in terms of the free NN scattering matrix K and subtract the Pauli–blocked contributions from inside the Fermi sphere [18,35]: R(k, k ′ ) = K (k, k ′ ) ? P d3 qK (k, q )GF (q )R(q, k ′ ) (32)

where GF denotes the two-body propagator inside the Fermi sphere. Dependencies on the two particle center–of–mass momentum and energy have been omitted. For the coupled partial wave channels Eq.(32) leads to a linear system of equations. Without attempting to solve Eq.(32) exactly an estimate of the density dependence is obtained by choosing k=kF and, in addition, also q=kF in the interactions beneath the integral. This amounts to the reasonable assumption that the momentum dependence of K and R over the Fermi sphere is much weaker than the variation of the (singular) quantity GF . The principal value integral over GF can be performed analytically and one ?nds that the medium corrections are given in leading order by the density of states at the Fermi surface, < GF >?NF =mkF /(π ? h)2 . (< GF > depends in fact also on the two–particle center–of–mass momentum and energy which we assume to be ?xed to appropriate average values). Thus, as an approximation to Eq.(32) we ?nd for the half–o? shell G–matrix R(kF , k ) ? 1 K (kF , k ) . 1 + K (kF , kF ) < GF > (33)

Next, we consider Eq.(32) for the case that a solution R0 at a Fermi momentum k0 F is known. In Eq.(33) this corresponds to insert in the denominator the subtracted Green function < GF > 0 ? < G0 F >. Since GF is a constant and K is weakly varying we may rewrite Eq.(33) as R(kF , k ) ? and a scaling factor

0 s(kF )=

1 1+

0 )R (k , k ) s(kF 0 F F

< GF >

0 s(kF )R0 (kF , k )

(34)

1 0 , k 0 ) < G0 > 1 ? R0 (kF F F

(35)

is found. These considerations suggest to use for the vertex functions (γ = 0, τ ) the ansatz ? ?

Nγ

where z=ρ/ρ0 and the base exponent is ?xed to β =1/3. The coe?cients aγ n are to be determined by the ?t to nuclear matter properties. As seen below nuclear matter allows a ?rst (Nτ =1) and third (N0 =3) order approximation for fτ and f0 , respectively. The overall scaling factors sγ renormalize the interaction strength asymptotically for ρ → 0 such that with an appropriate readjustment of the an coe?cients the nuclear matter saturation properties are una?ected. From the above discussion it is seen that the scaling factors relate the interaction vertices to the initial renormalization point γ k0 F . In the following, sγ and an are treated as empirical parameters. The ground state energy of nuclear matter is E=

α k≤kF

fγ (ρ) = sγ ?1 +

n=1

nβ ? aγ nz

(36)

? 2 k2 h 1 + 2m 2

? |αk, α′ k′ αk, α′ k′ | V

′ αα′ k≤kF k′ ≤kF

(37)

7

? includes antisymmetrization as discussed above and α represents the spin an isospin where V quantum numbers σ and q . In nuclear matter the single particle wave functions |αk are plane waves and the Fermi momenta kFq of protons and neutrons, respectively, is related to the corresponding densities ρq by kFq = (3π 2 ρq )1/3 . The total binding energy per particle in nuclear matter is obtained as: 3 τp ρp + τn ρn E (ρ, ρτ ) = A 5 ρ 1 d ?qq ?e + ρq ρq ′ V ′ (ρ) + Vqq ′ (ρ; kFq , kF ′ ) . q 2ρ ′

qq

2 h ? 2 kF

(38)

Here, ρ = ρn + ρp and ρτ = ρn ? ρp are the isoscalar and isovector density, respectively, and τq = 2m q denotes the Fermi energy of the protons and neutrons. The e?ective density dependent interactions appearing in the above expression are given in nuclear matter as the volume integral of V ?d (s, ρ) for the Hartree contributions, e.g. ? d ′ (ρ) = V qq

?d d3 s Vqq ′ (s, ρ),

(39)

while the non-local Fock contributions are averaged over the Fermi spheres, e.g. ? e ′ (ρ; kFq , kF ′ ) = V qq q

?e d3 s ρSL (skFq )ρSL (skFq′ ) Vqq ′ (s, ρ)

(40)

For Yukawa form factors the exchange integrals can be evaluated in closed form which, however, is omitted here. In symmetric nuclear matter, Eq.(38) simpli?es to [20,24] E 3 1 ? d (ρ) + V ? e (ρ, kF ) (ρ, ρτ = 0) = τ + ρ V 0 0 A 5 2 3 1 = τ + ρf0 (ρ) d3 s V0d (s) + 5 2

h ? 2 k2

(41)

e d3 s ρ2 SL (skF ) V0 (s) .

(42)

F . For given s0,τ and exponent β the remaining parameters of where kF = (3/2π 2 ρ)1/3 and τ = 2m f0 are determined by imposing the conditions for the nuclear matter binding energy, the saturation density and the compressibility at the equilibrium density ρ0 ,

E (ρ, ρτ = 0) A P = ρ2

= ? aV

ρ=ρ0

(43) (44) (45)

d (E/A) dρ

=0

ρ=ρ0

K ∞ = 9 ρ2

d2 (E/A) dρ2

ρ=ρ0

and the isovector vertex renormalization fτ is obtained from the symmetry energy ετ (ρ), de?ned by E (ρ, 0) E (ρ, ρτ ) = + A A ρτ ρ

2

ετ (ρ) +

ρτ ρ

4

ε4 (ρ) + . . .

.

(46)

At the saturation density ρ = ρ0 the ?rst term in Eq.(46) gives the binding energy per particle of symmetric nuclear matter, ?aV , and ετ is the symmetry energy coe?cient as = ρ2 0 ? 2 E (ρ0 , ρτ ) ?ρ2 τ (47)

ρτ =0

which is known to range between 28 MeV and 32 MeV [20,36]. 8

For a ?nite range interaction the above procedure for ?xing the isovector properties of asymmetric nuclear matter involves higher order contributions from the intrinsic momentum structure of the interaction. They appear because the asymmetry in the proton and neutron Fermi momenta induces an e?ective isospin symmetry breaking component via the exchange integrals extending over the two Fermi spheres. Numerically, the mixing of the isoscalar and the isovector channel is noticed immediately from the separate averaging of the exchange contributions over the proton and neutron Slater densities ρSL (kF p ) and ρSL (kF n ), Eq.(41). In the direct part or for a momentum independent contact interaction, respectively, these contributions are absent. In those cases the potential energy can always be written in the form ρ2 V0 +ρ2 τ Vτ , which separates exactly the isoscalar and the isovector part of the interaction (except for the case that V0 would depend intrinsically also on ρτ ). Skyrme interactions account to some extent for the isovector exchange e?ect by the momentum dependent components. At extreme high asymmetries, however, di?erences to a ?nite range interaction might become noticeable.

TABLE I. Parameters of the unrenormalized M3Y interaction in the odd-even (?rst block), the spin-isospin (second block), the p-n representation as de?ned (third block) and the e?ective spin-isospin representation (fourth block) as de?ned in App. A. r0 [fm] ?c m= h [MeV] r0 VSE VT E VSO VT O V0 Vσ Vτ Vστ d Vpp d Vpn e Vpp e Vpn V0d Vτd V0e Vτe 0.250 780 11466.00 13967.00 -1418.00 11345.00 11062.00 939.00 314.00 -969.00 11375.00 10748.00 -5642.00 2594.00 11062.00 314.00 -1524.00 -4118.00 0.400 490 -3556.00 -4594.00 950.00 -1900.00 -2538.00 -36.00 224.00 450.00 -2314.00 -2761.00 536.00 -1574.00 -2538.00 224.00 -519.00 1055.00 1.414 140 -10.46 -10.46 31.39 3.49 0.00 0.00 0.00 3.49 0.00 0.00 -5.23 -10.46 0.00 0.00 -7.85 2.62 Volume integral [MeV fm3 ] -980.28 -1324.02 1600.77 823.43 131.17 155.40 241.33 295.53 372.51 -110.16 -862.65 -1127.93 131.17 241.33 -995.29 132.64

9

B. Results for In?nite Nuclear Matter

The isoscalar vertex function f0 (ρ) was adjusted to the empirical values of the binding energy (aV ? 16 MeV), the saturation density (ρ0 ?0.159 fm?3 ) and the compressibility K∞ . Empirically, K∞ is not well known (values range from 200 MeV to 350 MeV), but recent studies, e.g. [36], converge to a value of about K∞ =230 MeV. For the symmetry energy coe?cient serving to ?x fτ as = ετ (ρ0 =0.159 fm?3 )=32 MeV is used. In Tab. II two sets of vertex renormalization parameters, Eq.(36) are displayed. The ?rst set was obtained with s0 = 1, i.e. leaving the free space interaction unchanged, whereas the second set with s0 = 5 renormalizes the asymptotic behavior of the interaction as well. This value of s0 was determined from the binding energies of ?nite nuclei as will be discussed in the next section.

TABLE II. Parameterization of f0 (ρ) and fτ (ρ) a0 1 s0 = 1 s0 = 5 sτ = 1 1.365 -2.127 aτ 1 .1688 a0 2 -1.839 2.032 a0 3 0.379 -0.724 ρ0 [fm?3 ] 0.159 0.159 K [MeV] 230 230 (EB /A) [MeV] -16.00 -16.00 ετ (ρ0 )[M eV ] 32.00

10

f0(ρ), s0 = 1 3 f0(ρ), s0 = 5 fτ(ρ), sτ = 1

2

f(ρ)

1 0 0,0 0,1 0,2

-3

0,3

nuclear density ρ [fm ]

FIG. 1. Density dependent isoscalar f0 (s0 = 1 dashed line, s0 = 5 solid line) and isovector fτ (dotted line) vertex functions of the D3Y interaction.

Fig. 1 shows the density dependence of the vertex correction functions f0 and fτ . In the isovector channel only a rather small deviation from unity is found, and fτ increases only slightly for higher densities. This re?ects the fact that the original interaction already has a symmetry coe?cient close to 32 MeV. In the isoscalar channel one realizes an increase of the interaction at low densities and a suppression of the interaction at higher densities. The low density behavior obviously scales with s0 . In Fig. 2 the equations of state (EoS) for the D3Y and the Skyrme SLy4 interaction [36] are compared. For s0 =1 the D3Y and the SLy4 results are in very close agreement over the full range of densities. According to the behavior of f0 , the binding energy increases at densities lower and decreases at densities higher than the saturation density for s0 >1 while retaining its value at the equilibrium density.

11

10 D3Y, s0 = 1 D3Y, s0 = 5 Skyrme SLy4 0

EB/A(ρ) [MeV]

-10

-20 0,0 0,1 0,2

-3

0,3

nuclear density ρ [fm ]

FIG. 2. Nuclear matter equation of state calculated with the density dependent D3Y interaction for the parameterizations s0 = 1 (dashed line) and s0 = 5 (solid line) and the Skyrme SLy4 interaction (dotted line).

The isovector properties of the D3Y interactions are illustrated in Fig. 3, where the EoS for asymmetric nuclear matter is displayed for di?erent asymmetry ratios Z/A = ρp /ρ. Stable nuclei typically cover the range Z/A = 0.5 . . . 0.4 (e.g. 40 Ca, 208 Pb). In exotic nuclei the new regions with Z/A = 0.4 . . . 0.3 become accessible and far of stability Z/A = 0.2 may be reached (e.g 10 He). It is remarkable that neutron matter is still bound but this occurs at extremely low densities. The reason for this behavior is the strong attraction of the D3Y with s0 = 5 at very low densities and, as such, is not very conclusive. In stable nuclei these density regions will rarely contribute to the total binding. However, a di?erent situation may be expected far o? stability, especially in halo nuclei, where the binding sensitively depends on the interactions at low density. Comparing the D3Y to the SLy4 interaction one sees that both interactions are in agreement around the saturation density for asymmetry ratios of Z/A = 0.5 . . . 0.4. The di?erences showing up at low densities and large asymmetries are mainly due to the stronger attraction of the D3Y isoscalar interaction with s0 =5.

12

30 D3Y, s0=5 SLy4 20

Z/A

neutron matter

0.2

EB/A(ρ) [MeV]

10

0 0.4 0.5 -10

-20 0,0 0,1 0,2

-3

0,3

nuclear density ρ [fm ]

FIG. 3. Equation of state for asymmetric nuclear matter calculated with the density dependent D3Y interaction for the parameterizations s0 = 5 (solid line) and the Skyrme SLy4 (dashed line). Shown is the EoS for neutron matter and for the asymmetry ratios Z/A = 0.2, 0.4, 0.5.

In Fig. 4 the EoS for di?erent asymmetry ratios and the two sets of parameters are examined. The equilibrium density ρ0 (ξ ), the energy per nucleon EB (ξ )/A, the symmetrie energy ?τ (ξ ) and the compressibility K∞ (ξ ) for ξ =Z/A are shown. A higher value of s0 is seen to lead to a lower equilibrium density for strong asymmetries. In addition the binding energy at the equilibrium density is higher due to the stronger attraction of this parameterization of the interaction at low densities. The lower left part shows the symmetry energy of the EoS in dependence of the proton excess. One realizes a decrease of the symmetry energy for strong asymmetries and a higher value of s0 corresponding to the shift of the saturation density to lower densities. The compressibility approaches zero at high asymmetries with a remarkable ?at slope.

13

0,20 s0 = 1 s0 = 5

5

0,15

0

E/A(ρ0,ξ) [MeV]

ρ0(ξ) [fm ]

-3

0,10

-5

-10 0,05 s0 = 1 s0 = 5 0,00 0,1 30 0,2 0,3 0,4 0,5 0,1 0,2 0,3 0,4 0,5 250 200 150 100 10 s0 = 1 s0 = 5 0 0,0 0,1 0,2 0,3 0,4 0,5 0,1 0,2 0,3 0,4 s0 = 1 s0 = 5 0,5 50 0 -15

K(ρ0,ξ) [MeV]

ετ(ρ0,ξ) [MeV]

20

ξ

ξ

FIG. 4. Equilibrium density ρ0 of the EoS (upper left), binding energy per nucleon EB /A at ρ0 (upper right), symmetry coe?cient ?τ at ρ0 (lower left) and the compressibility K∞ (ρ0 ) (lower right ), all plotted as a function of the neutron excess or asymmetry ratio ξ = Z/A. Results are shown for the D3Y interaction with s0 = 1 and s0 = 5.

The properties of the D3Y interaction in the particle-particle channel were investigated by considering like–particle pairing in symmetric nuclear matter. Pairing, being mediated by the singlet-even part VSE of the interaction, was examined by Hartree-Fock-Bogoliubov (HFB) calculations. The in–medium SE particle–particle interaction was obtained by back transformation from the renormalized particle–hole interactions using the representation in terms of the exchange interactions ? ?d (see App. A), leading to VSE (ρ) = 2 γ =0,τ Vγ (ρ) + Vγ?e (ρ) . As in refs. [37,38] the following set of equations was solved simultaneously: ρ(kF ) = 4 v 2 (k, kF ) = ?(k, kF ) = 1 2 d3 k 2 v (k, kF ) (2π )3 1? ?(k, kF ) ? λ (?(k, kF ) ? λ)2 + ?2 (k, kF ) = 1 ? u2 (k, kF ) (48) (49) (50) (51)

? 2 k2 h d3 k ′ kk ′ |V ? |kk ′ ? k ′ k v 2 (k ′ , kF ) +4 2m (2π )3 d3 k ′ ? ?(k, kF ) = ? k ? k |VSE |k ′ ? k ′ u(k ′ , kF )v (k ′ , kF ) (2π )3

14

0

-4

aSE [fm]

-8

-12

D3Y, s0 = 1 D3Y, s0 = 5 bare M3Y

-16 0,0 0,2 0,4 0,6 0,8 1,0 1,2

-1

1,4

1,6

Fermi momentum kF [fm ]

FIG. 5. Scattering length app SE (kF ) for relative s–wave states. Results are shown for the bare M3Y (dotted) and the D3Y with the parameterizations s0 = 1 (dashed line) and s0 = 5 (solid line).

Fig. 5 shows the calculated in–medium scattering length for like particles (pp/nn) app SE (kF ) at a given Fermi momentum kF for di?erent parameterizations of the interaction. Because of the density dependence in the isovector and in the isoscalar channel the strength of the interaction increases at low to intermediate densities. For the D3Y with s0 = 5 the pairing interaction is enhanced by a factor of about 8 at low densities over the values of the original M3Y because of the asymptotic scaling of the isoscalar interaction. This corresponds to about 90% of the s-wave strength of the free NN T-matrix with a scattering lenght of app SE = ?17.1fm and leads to higher values of the gap in nuclear matter. In Fig. 6 ?(k, kF ) at k=kF is shown. For the bare M3Y the position of the maximum gap lies at about kF ?0.8 fm?1 . About the same result is obtained in HFB calculations with the Gogny force [37,38] and for several Skyrme interactions [39]. The position of the maximum gap in the D3Y–calculations is shifted to slightly lower densities with kF ?=0.75 fm?1 resp. kF ?0.63 fm?1 for s0 =1 resp. s0 =5. One ?nds that the height of the maximum increases by lowering the position of the maximum. Our results of ?(kF , kF ) ? 2 MeV resp. ?(kF , kF ) ? 4 MeV for the M3Y and the D3Y with s0 = 1 are in agreement with the Gogny force that yields a value of ?(kF , kF ) ? 3 MeV [37,38]. The calculated maximum gap for the D3Y with s0 = 5 has a much higher value of about 9.5 MeV. The reason for this is the much stronger pairing interaction, especially at lower densities, as seen in Fig. 5. Comparing the value of the gap at the saturation density kF = 1.33f m?1 with other calculations we ?nd a value of about 0.2 MeV for the M3Y and a nearly vanishing value for the D3Y interaction. For the Gogny force one ?nds a value of about 0.6 MeV but most other calculations favor a small to vanishing gap [39–41].

15

10

D3Y, s0 = 1 D3Y, s0 = 5 Gogny D1

8

? (k=kF,kF) [MeV]

6

4

2

0 0,0 0,2 0,4 0,6 0,8 1,0 1,2

-1

1,4

1,6

Fermi momentum kF [fm ]

FIG. 6. Gap at the Fermi momentum in nuclear matter. Results for the D3Y interaction with s0 = 1 (dotted line) and s0 = 5 (solid line) and the Gogny D1 (dashed line) are shown.

IV. HARTREE–FOCK CALCULATIONS FOR FINITE NUCLEI

The density dependent D3Y interaction derived in the preceding section for nuclear matter is applied in HF–calculations for ?nite nuclei. Stable nuclei close to the N?Z line are chosen ?rst. In these calculations mainly the isoscalar properties of the interaction are investigated. A more general test, including also isovector dynamics, is obtained from calculations for even-even Sn– isotopes covering the range from the proton dripline (100 Sn) to the A=140 mass region. The calculations for stable and unstable nulcei are performed without readjustment of parameters from their values derived for nuclear matter. The general structure of the HF–equations for a density dependent NN-interaction in the DME approach was already derived in Sect. II B. In the practical applications the DME is used only for the exchange contributions to the HF potentials. Especially at higher density and thus Fermi momentum exchange is dominated by the short range pieces of the interaction for which the local momentum approximation underlying the DME is most appropriate. Di?erent from the original work of Negele and Vautherin [8] the direct, i.e. Hartree potentials are treated in full ?nite range. The D3Y density dependence was introduced as a local renormalization of the vertices which requires to replace the interaction by [18] Vγ?d (r1 , r2 ; ρ) = gγ (r1 )Vγd (r1 ? r2 )gγ (r2 ) (52)

2 with gγ = f0,τ and where gγ (r) = gγ (ρ(r)) are the in–medium isoscalar and isovector vertex renormalizations, respectively. Vγd (r1 ? r2 ) denotes the unrenormalized ”bare” interaction in these channels. This ensures a reliable treatment of the density dependence also in the surface region where the vertices vary rapidly (see Fig. 1). The direct potentials are obtained as in Eq.(15) by folding the proton and neutron densities, respectively, with the corresponding – now density d dependent – direct interactions. To accomplish this we express the interactions Vqq ′ (r1 , r2 ) by ?d Vγ (r1 , r2 ; ρ). With the DME the exchange or Fock potentials are replaced by local potentials whose strengths are de?ned in Eq.(8). In order to preserve this separation of the center–of–mass and relative coordinates r and s the interaction must also be separable in r and s. Assuming that

16

the resulting error is small compared to the error introduced by the DME we therefore replace in the exchange interaction, as a ?rst approximation, the individual coordinates of the interacting particles by their center–o?–mass coordinate. This leads to

2 Vγ?e (s, r; ρ) = gγ (r)Vγe (s) e by which we express the exchange interactions Vqq ′ (s).

(53)

A. Details Of Numerical Calculation

For our calculations we only consider the stationary ground state of spherical even-even nuclei, therefore the single particle wave functions can be separated as φkq (r) = Rαq (r) Ψljm (θ, φ) r (54)

The functions Ψljm (θ, φ) are spinor spherical harmonics [42], the index αq represents the quantum numbers αq = {n, l, j } for protons and neutrons, resp. This allows us to simplify the density and the kinetic energy density to ρq (r) = 1 4πr2 1 4π

2 2 (r) (2jαq + 1)Rα vα q q αq 2 (2jαq + 1) vα q αq

(55)

2

τq (r) =

?r

Rαq (r) r

+

lαq (lαq + 1) r2

Rαq (r) r

2

(56)

2 for a state αq are determined by pairing correlations. Since it was The occupation weights vα q found that pairing gives minor to negligible contributions to the nuclei considered here, the BCS approximation was used. A constant pairing matrix element of Gq = 23 MeV/nucleon was used and the standard set of BCS equations [20] was solved independently for protons and neutrons, respectively. For protons also the Coulomb interaction contributes to the Hartree–Fock potential U HF :

UCoul (r) =

d3 r1 ρp (r1 )

e2 |r ? r1 |

?

3 π

1/3

e2 ρp (r)1/3

(57)

The Coulomb exchange term is treated in the local density approximation rather than in the full DME because its contribution to the potential energy is small. The spin-orbit potential is treated as in [8,43], leading to a contribution to the potential energy Es.o. = d3 r Hs.o. (r) with 1 Hs.o. (r) = ? W0 2 ρ(r) ?J(r) +

q

ρq (r) ?Jq (r) .

The spin-orbit strength W0 is related to the two-body spin-orbit potential in the short–range limit by W0 = ? 2π 3 VLS (s)s4 ds (58)

Using the M3Y parameterization [21] one ?nds a value of W0 ? 105 MeV fm5 but in correspondence to most Skyrme interactions W0 = 120 MeV fm5 from reference [43] is chosen. After variation this term gives a contribution 1 s.o. Uq (r) = ? W0 [ρ(r) ?J(r) + ρq (r) ?Jq (r)] , 2 to the Hartree–Fock potential and a single particle spin-orbit potential 17 (59)

ls Uq (r) =

1 1 d W0 [ρ(r) + ρq (r)] l · s 2 r dr

(60)

where the spin density is de?ned as Jq (r) = 1 4πr3 (2jαq + 1) jαq (jαq + 1) ? lαq (lαq + 1) ?

αq

3 2 (r). Rα q 4

(61)

B. Results for Stable Nuclei

The various contributions to the nuclear part of the self–consistent DME mean–?eld in 40 Ca are shown in Fig. 7. The isoscalar Hartree potential is weak and repulsive and binding is essentially mediated by the Fock parts. The Fock potential is strongly attractive with a depth at the center of about -110 MeV. The sum of the Hartree and Fock terms leads to an average potential depth of -70 MeV. The direct and exchange rearrangement potentials reduce the depth especially close to the center. The mean–?eld agrees qualitatively with other non–relativistic HF results, e.g [43]. It also resembles closely the Schroedinger-equivalent potential from a relativistic density dependent Hartree calculation [17,18] except for a more pronounced shell structure. Of special interest is the in?uence of the rearrangement potential whose in?uence has proven to be important as found in comparable non-relativistic [43] and relativistic calculations [17,18]. Rearrangement is of importance for the single particle potential and the energy spectrum. Because the rearrangement potential is always repulsive it lowers the single particle potential and leads to a change in the single particle spectra and the charge radii. In the calculation of the total binding energy the rearrangement e?ects cancel out. Taking this into account the total binding energy is given as

A

EB =

i=1

?i ?

1 2 q=p,n

HF d3 rρq (r)Uq (r) ? q=p,n

R Eq + Epair ,

(62)

with the rearrangement energy calculated from Eq.(27) and Eq.(28):

R Eq = R d3 r ρq (r) Uq (r) +

d3 r τq (r)

?2 h 2

1 m? q (r)

?

1 m

(63)

18

20 0

UC [MeV]

-20 -40 -60 -80

U0 U0 U0

R tot

20 0 -20

UC [MeV]

-40 -60 -80 -100 -120 0 2 4 6 8 U0 U0 U0 U0

d d,R e e,R

r [fm]

FIG. 7. Single particle potential for 40 Ca calculated with the D3Y with s0 = 5. In the lower part, d e contributions to the total isoscalar potential from the direct U0 and the exchange term U0 of the normal d,R e,R Hartree-Fock potential as well as the rearrangement potentials (U0 , U0 ) are shown. In the upper tot part, the total isoscalar potential U0 as well as its contributions, the attractive HF potential U0 and the R , are shown. repulsive rearrangement potential U0

An important contribution to mean–?eld dynamics is provided by the momentum dependence of the interaction. Since the exchange parts of the D3Y interaction were found to be most important for nuclear binding the question arises whether the momentum structure is described realistically. The e?ective mass m? q , Eq.(27) is completely determined by the Fock terms, since they are the only source of non–localities in the present model. From Fig. 8 it is seen that a very reasonable 208 value of about m? Pb which closely q /m ?0.68 is obtained in the interior of a heavy nucleus like agrees with other theoretical results and determinations from proton–nucleus elastic scattering [44]. It is worthwhile to emphasize that this result was obtained without readjustments of parameters. The D3Y interaction (or better to say the original M3Y G–matrix [21]) has apparently a 19

realistic mixture of short and long range contributions accounting reasonably well for the intrinsic momentum structure of an in–medium HF interaction. The results of this section also con?rm the DME approach. Especially promising results are obtained in conjunction with the e?ective ”Fermi momentum” qF , Eq.(6), by which non–local e?ects are accessible in a self–consistent way.

1,0

16

1,0

O

0,9 0,8 0,7 0,6

40

Ca

0,9

m (r)/m

0,8 0,7 0,6 0 1,0

48

*

2

4

6

8 1,0

0

2

4

6

8

10

Ca

0,9 0,8 0,7 0,6

208

Pb

0,9

m (r)/m

0,8 0,7 0,6 0 2 4 6 8 10

*

0

2

4

6

8

10

12

r [fm]

r [fm]

FIG. 8. E?ective mass of protons (solid line) and neutrons (dashed line) in closed shell nuclei calculated with the density dependent D3Y interaction (with the parameterization s0 = 5.)

TABLE III. interaction s0 = 1 s0 = 5 exp.

Charge radii rC [fm] for closed shell nuclei for di?erent parameterizations of the D3Y

16 O 2.79 2.88 2.73 40

Ca 3.51 3.61 3.49

48

Ca 3.52 3.60 3.47

208

Pb 5.51 5.59 5.50

20

0,10 0,08

16

0,10

O

0,08 0,06 0,04 0,02 0,00 0,10

40

Ca

ρC(r) [fm ]

-3

0,06 0,04 0,02 0,00 0,10 0,08

48

0

2

4

6

0

2

4

6

8

Ca

0,08 0,06 0,04 0,02 0,00

208

Pb

ρC(r) [fm ]

-3

0,06 0,04 0,02 0,00

0

2

4

6

8

0

2

4

6

8

10

r [fm]

r [fm]

FIG. 9. Charge density distributions of closed shell nuclei from density dependent Hartree-Fock calculations using the D3Y interaction. Results for s0 = 1 (dashed line) and s0 = 5 (dotted line) are shown. Experimental densities and their uncertainties are denoted by solid lines (taken from Ref. [45]) .

Fig. 9 compares the calculated charge density distributions of the above closed–shell nuclei with experimental data from elastic electron scattering [45]. The theoretical point particle distribution was folded with a Gaussian proton form factor [43] with r2 p = 0.8f m. Tab. III shows the calculated charge radii in comparison to the experimental data (taken from [18]). The charge density in the nuclear interior is slightly overestimated for heavy nuclei like 208 Pb and slightly underestimated for the lighter nuclei. In all cases, the density at large radii is higher than the experimental data which leads to a overestimation of charge radii as can be seen in Tab. III. In general the agreement with the experimental charge radii is satisfying and is, except for 16 O, ≤ 2% for the D3Y with s0 = 1 and ≤ 4% for the D3Y with s0 = 5. In contrast to this, the binding energies calculated with s0 = 1 are strongly underestimated by more than 0.5 MeV per nucleon, whereas the D3Y with s0 = 5 reproduces the experimental data [46] very well (Tab. IV). One sees that an improvement in the binding energies leads to higher values of the charge radii. In general it is not possible to reproduce simultaneously both the binding energies and the charge radii with one parameter set. This is in agreement with calculations in the relativistic density dependent ?eld theory which typically give comparable results [13,15]. But generally the agreement with the experimental binding energies of the D3Y with s0 = 5 is surprisingly good as can be seen in Fig. 10 realizing that the number of used parameters is lower than the one in Gogny or Skyrme interactions. The observed behavior of the D3Y for di?erent parameterizations shows the importance of the density dependence of the interaction. Although both parameterizations have the same properties in the range of the saturation density in nuclear matter the results for ?nite nuclei are quite di?erent. This behavior is caused by the di?erences at low densities. In addition, the rearrangement terms enhance the e?ect of the density dependence.

21

-7,0

-7,5

EB/A [MeV]

-8,0

-8,5

D3Y, s0=1 D3Y, s0=5 exp.

-9,0 0 50 100 150 200

A

FIG. 10. Binding energy per nucleon for closed shell nuclei from density dependent Hartree-Fock calculations using the D3Y interaction. Results for s0 = 1 (triangles) and s0 = 5 (circles) are shown. Experimental binding energies are denoted by open squares (taken from Ref. [46]). Binding energies are shown for the nuclei 16 O, 40 Ca, 48 Ca, 90 Zr, 132 Sn and 208 Pb. Lines are drawn to guide the eye.

TABLE IV. Binding energy per nucleon EB /A [MeV] for closed shell nuclei for di?erent parameterizations of the D3Y interaction s0 = 1 s0 = 5 exp. O -7.07 -8.18 -7.98

16

Ca -7.74 -8.62 -8.55

40

Ca -7.91 -8.69 -8.67

48

Zr -8.00 -8.66 -8.70

90

132

Sn -7.77 -8.38 -8.36

208

Pb -7.32 -7.88 -7.87

22

C. The Ground States of Sn Isotopes

The Sn isotopes are of particular interest for nuclear structure and also astrophysical questions because of the closure of the Z=50 proton shell. The known isotopes [46] cover the range from the proton dripline at 100 Sn [47] to the doubly magic 132 Sn nucleus which is already β –unstable. Here, we are mainly interested to investigate the isovector properties of the D3Y interaction and, as a more general aspect, to test an e?ective interaction, determined from symmetric nuclear matter and stable nuclei, in regions far o? stability. No attempt is made to optimize the parameters in order to obtain a perfect ?t of data. The parameter set with s0 =5 is used. Theoretical and experimental binding energies per particle are compared in Fig. 11. The D3Y binding energies are seen to be shifted slightly to the left of the data [46] but a fair description of the overall dependence on mass and asymmetry is obtained. The strongest binding is obtained for 120 Sn while experimentally the minimum is found at approximately 115 Sn. The binding energies of the neutron–rich isotopes are reasonably well reproduced but in view of the deviations on the neutron–poor side the agreement might be fortuitous. Towards the proton dripline, i.e. N→Z, the binding energies are increasingly underestimated by up to about 0.12 MeV per nucleon close to 100 Sn. This result is to some extent unexpected because stable nuclei with N?Z were well described. However, also the binding energy of 90 Zr (Fig. 10 is underestimated the most and apparently the e?ect is enhanced in going from the N=50 to the Z=50 shell. The deviations are in a range typical for calculations with a G–matrix in ?nite nuclei. Fully phenomenological interactions like the Skyrme SLy4 force, for which results are also displayed in Fig. 11, and the Gogny D1 interaction [19] reproduce the Sn binding energies almost perfectly.

-7,9

-8,0

Sn Isotopes

-8,1

EB/A [MeV]

-8,2

-8,3

-8,4

-8,5

D3Y, s0 = 5 Skyrme SLy4 exp.

100 110 120 130 140

-8,6

A

FIG. 11. Binding energy per nucleon for the tin isotopes from density dependent Hartree-Fock calculations using the D3Y interaction with s0 = 5 (triangles) and the Skyrme SLy4 force [36] (cirles). Experimental binding energies are denoted by open squares (taken from Ref. [46]). The lines are drawn to guide the eye.

Variations of the scaling parameter s0 resulted in an overall shift of the binding energy curve as it might be expected. Improvements on the neutron–poor side are obtained on the expense of a worse description of the neutron–rich isotopes. In particular, the minimum was found to be rather insensitive on variations of s0 . For s0 =1 the whole curve is shifted up by about 0.5 MeV per nucleon with essentially the same shape. These observations lead to the conclusion that the 23

parameters a0 n themselves should be readjusted on ?nite nuclei with constraints to nuclear matter rather than retaining the nuclear matter values. Taking the point of view that the nuclear matter parameterization accounts for the static part of the interaction, where cluster contributions are included e?ectively by the ?tting procedure [27], such a re?nement could be related to non–static contributions to the in–medium interaction because of the larger polarizability of the nuclear surface. Theoretically, this would amount to go beyond the ladder approximation and to include also ring diagrams. Such contributions are obviously contained in e.g. the Skyrme parameters but missing in the D3Y interaction. Over the range of the Sn isotopes the polarization contributions can be expected to decrease towards A?132 which is known to be a rather sti? nucleus.

0,10

110

0,08

Sn Sn Sn Sn

120 132

ρn(r) [fm ]

-3

0,06

140

0,04

0,02

0,00 0 0,08

110

2

4

6

8

10

Sn Sn Sn Sn

0,06

120 132

ρC(r) [fm ]

-3

140

0,04

0,02

0,00 0 2 4 6 8 10

r [fm]

FIG. 12. Charge density distributions (lower part) and neutron density distributions (upper part) of the tin isotopes from density dependent Hartree-Fock calculations using the D3Y interaction. Results for s0 = 5 are shown.

24

Neutron density distributions and charge densities for several Sn isotopes are displayed in Fig. 12. From 110 Sn to 140 Sn the charge distributions are slightly reduced in the interior. As seen in Fig. 13 this is accompanied by a mild increase of the charge radius by about 5%. A more drastic evolution is found for the neutron densities. Beyond 132 Sn an extremely thick neutron skin is building up leading to a sudden jump in the neutron rms radii. The neutron skin thickness is more clearly visible in Fig. 14 where the di?erence of the proton and neutron rms radii is shown. In 140 Sn a solid layer of neutron matter is found extending by about 0.48 fm beyond the proton matter distribution which is more than a factor of two larger than the neutron skin in 208 Pb. The same phenomenon is predicted also by the SLy4 interaction except for a somewhat smaller neutron excess radius. The saturation of the rms values around A=132 is an indication of the double magic nature of 132 Sn. The increase is directly related to the shell structure in the heavy tin isotopes. At A=132 the 2f7/2 shell is ?lled and pairing does not contribute. At larger masses the neutron 3p–subshells become populated. Weak binding and the low angular momentum barrier allow a large extension of the valence wave functions into the exterior. However, the calculations do not lead to a ”halo”.

5,4

5,2

Sn Isotopes

<r >

2 1/2

[fm]

5,0

4,8

4,6

neutrons protons

4,4 100 110 120 130 140

A

FIG. 13. Root mean square radii of the tin isotopes from density dependent Hartree-Fock calculations using the D3Y interaction with s0 = 5 for neutrons (squares) and protons (triangles). The lines are drawn to guide the eye.

25

0,5 0,4

Sn Isotopes

[fm] - <r p> <r n>

2 1/2 2 1/2

0,3 0,2 0,1 0,0 -0,1 -0,2 100 110 120 130 140

D3Y, s0 = 5 SLy4

A

FIG. 14. Di?erence of the root mean square radii of neutrons and protons of the tin isotopes for the D3Y interaction with s0 = 5 (triangles) and the Skyrme SLy4 (cirles). The lines are drawn to guide the eye.

An interesting observation is made from Fig. 14 on the neutron–poor side. Both the D3Y and the SLy4 interactions predict a proton skin for A≤106. The origins of the proton and neutron skins are quite di?erent. Atthe neutron rich side the excess neutrons become less bound because of the increasing repulsion from the isovector potential. Since states with ? ≥1 are involved they are stored predominantly in the surface region of the nucleus. In the heavy isotopes the isovector potential acts attractively in the proton sector and thus balances to a large extent the Coulomb repulsion. Since for N→Z the isovector part of the mean–?eld becomes strongly suppressed the repulsion is no longer compensated and Coulomb e?ects become visible. Hence, the proton skin in the neutron–poor Sn nuclides is caused solely by the Coulomb interaction. The calculations were extended further into the β –unstable region. Up to A=160 the neutron dripline is still not reached. The separation energies are small but decrease rather slowly.

V. SUMMARY, DISCUSSION AND CONCLUSIONS

The density matrix expansion (DME) method was used to derive a local energy density functional for non–relativistic meson exchange interactions. The interaction energy density was separated into direct and exchange contributions. The Hartree contributions are treated exactly. From the systematic expansion of the one–body density matrix provided by the DME an average momentum density was determined such that the next to leading order terms in the DME series are canceled. In uniform systems like in?nite nuclear matter all higher order terms are canceled and the exact energy density is recovered. In ?nite nuclei the approach corresponds to a generalized Slater expansion of the density matrix. Hartree–Fock equations were derived. With our choice for q2 F a self–consistently generated e?ective mass was obtained which is completely determined by the intrinsic momentum structure of the underlying interaction. A semi–microscopic approach to the in–medium NN interaction was discussed. The M3Y G– matrix from the Paris NN potential was modi?ed by introducing density dependent vertex renormalizations in the isoscalar and isovector particle–hole channels. The functional form of the vertex functions was chosen in close analogy to the dynamical structure and the medium dependence 26

of the vertices in a Brueckner G–matrix. The parameters, however, were determined empirically by a ?t to the binding energy, equilibrium density, compressibility and the symmetry energy of symmetric nuclear matter. The momentum structure of the original G–matrix was retained. The probably ?rst attempt to parameterize a G–matrix in terms of ?xed propagators and density dependent vertices goes back to Sprung and Banerjee [48]. They approximated the momentum structure of the singlet/triplet even-odd interactions obtained from the Reid NN potential by a superposition of 5 arbitrarily chosen gaussians with channel dependent in–medium vertices. The functional form of the vertices was chosen in the same way as in section III A but using always a ?rst order polynomial. The nuclear matter binding energy curve of the full Reid G–matrix was reproduced the best for a linear dependence of the vertices on kF , corresponding to β =1/3 in Eq.(36). From the description of individual matrix elements a slight preference for a lower value β =1/6 was deduced [48]. However, nuclear matter was found to be badly underbound by about 5 MeV since only two-body correlations were taken into account. The renormalized vertices obtained in our semi–microscopic approach include e?ectively also contributions beyond the two– body cluster approximation, as was pointed out some time ago by Bethe [27]. In nuclear matter the D3Y interaction leads to an equation of state with an unusual enhancement of attraction at low densities which is not found for phenomenological forces. Unfortunately, G– matrix calculations are not available at such low densities but obviously the free NN T–matrix should be approached asymptotically. Phenomenological interactions and also the unrenormalized M3Y G–matrix do not account for this transition. Without putting too much emphasis on the details of this particular result it might be expected that such a region of enhanced attraction could indeed exist in the far tails of nuclear wave functions. Before drawing ?nal conclusions on this subject a re?ned treatment of the low density region is clearly necessary. For that purpose an extended version of the schematic model for the in–medium vertices might prove to be useful. The ground state properties of stable nuclei were reproduced reasonably well. The good description of binding energies gives con?dence to the DME approach, especially in view of the fact that an interaction from nuclear matter was used. Acknowledgement: This work was supported in part by DFG (contract Le439/4-1), GSI and BMBF.

APPENDIX A: REPRESENTATIONS OF THE NN INTERACTION

In literature several representations of the NN interaction exist which can be found in, e. g. [22]. Here we will only consider the central part of the interaction. The M3Y interaction is given in a representation specifying the two-body channel spin and the parity (SE,TE,SO,TO) which is commonly used in inelastic scattering. In our calculations we are especially interested in the p-p resp. p-n interaction for calculating the Hartree-Fock matrix elements, the (pp, pn) representation. It can be easily obtained with the transformation [49]

d Vpp = d Vpn e Vpp e Vpn

1 [2VSE 8 1 = [1VSE 8 1 = [2VSE 8 1 = [1VSE 8

d + 6VT O ] = Vnn d + 3VT E + 1VSO + 3VT O ] = Vnp e ? 6VT O ] = Vnn e + 3VT E ? 1VSO ? 3VT O ] = Vnp

(A1)

Here, the indices d resp. e stand for the direct resp. exchange part of the NN interaction. A representation especially useful if one uses a contact interaction or is interested in the direct part of the interaction is the set (V0 , Vσ , Vτ , Vστ ), which follows from the set (SE,TE,SO,TO) through 1 [3VSE + 3VT E + 1VSO + 9VT O ] 16 1 [?3VSE + 1VT E ? 1VSO + 3VT O ] Vσ = 16 V0 = 27

Vτ = Vστ

1 [1VSE ? 3VT E ? 1VSO + 3VT O ] 16 1 [?1VSE ? 1VT E + 1VSO + 1VT O ] = 16

(A2)

Our aim is to ?nd a spin averaged interaction which separates the isoscalar and the isovector interaction for both the direct and the exchange part. De?ning V0

d/e

Vτd/e

1 d/e d/e d/e d/e (V + Vpn + Vnp + Vnn ) 4 pp 1 d/e d/e d/e d/e = (Vpp ? Vpn ? Vnp + Vnn ) 4 =

(A3)

we ?nd in the (SE,TE,SO,TO) representation V0d = 1 [3VSE 16 1 [1VSE Vτd = 16 1 [3VSE V0e = 16 1 [1VSE Vτe = 16 + 3VT E + 1VSO + 9VT O ] = V0 ? 3VT E ? 1VSO + 3VT O ] = Vτ + 3VT E ? 1VSO ? 9VT O ] ? 3VT E + 1VSO ? 3VT O ] . (A4)

One sees that the de?ned direct terms of this representation are equivalent to the ones of the set (V0 , Vσ , Vτ , Vστ ) but also that there is no equivalent for the exchange terms.

[1] W. Schwab, H. Geissel, H. Lenske, K.-H. Behr, A. Br¨ unle, K. Burkhard, H. Irnich, T. Kobayashi, G. Kraus, A. Magel, G. M¨ unzenberg, F. Nickel, K. Riisager, C. Scheidenberger, B. M. Sherill, T. Suzuki, B. Voss, Z.Phys. A350 (1995) 283. [2] M. Zinser, H. Emling, H. Lenske et al., Nucl.Phys A (in print); N. Ostrowski. H. G. Bohlen, A. S. Demyanova, B. Gebauer, R. Kapakchieva, Ch. Langner, H. Lenske, M. von Lucke-Petsch, W. von Oertzen, A. A. Ogloblin, Y. E. Penionzhkevich, M. Wilpert and Th. Wilpert, Z.Phys. A343 (1992) 389; H. G. Bohlen, B. Gebauer, M. von Lucke-Petsch, W. von Oertzen, A. N. Ostrowski and H. Lenske, Z.Phys. A344 (1993) 381. [3] T. H. R. Skyrme, Nucl. Phys. 9 (1959) 615. [4] K. A. Brueckner, J.R. Buchler, S. Jorna and R.J. Lombard, Phys. Rev. 171(1968) 1188. [5] H. S. Koehler, Phys. Rep. 18C (1975) 219. [6] R. K. Tripathi, A. Faessler, and H. M¨ uther, Phys. Rev. C10 (1974) 2080. [7] H. K¨ ummel, K. H. L¨ uhrmann and J. G. Zabolitzky, Phys. Rep. 36 (1978) 1. [8] J. W. Negele and D. Vautherin, Phys. Rev. C5 (1972) 1472. [9] J. W. Negele, Rev. Mod. Phys. 54 (1982) 913. [10] F. T. Baker, L. Bombot, C. Djalali, C. Glashauser, H. Lenske, W. G. Love, M. Morlet, E. Tomasi– Gustafson, J. Van der Wiele, J. Wambach, A. Willis, Phys.Rep. (in print). [11] M. R. Anastasio, L. S. Celenza, W. S. Pong and C. M. Shakin, Phys. Rep. 100 (1983) 327. [12] C. J. Horowitz and B. D. Serot, Nucl. Phys. A464 (1987) 613. [13] B. ter Haar and R. Mal?iet, Phys. Rep. 149 (1987) 207. [14] H. M¨ uther, R. Machleidt and R. Brockmann, Phys. Lett. B202 (1988) 483. [15] R. Brockmann and R. Machleidt, Phys. Rev. C42 (1990) 1965. [16] R. Brockmann and H. Toki, Phys. Rev. Lett. 68 (1992) 3408. [17] H. Lenske and C. Fuchs, Phys. Lett. B345 (1995) 355. [18] C. Fuchs, H. Lenske and H. H. Wolter, Phys. Rev. C52 (1995) 3043. [19] J. Decharge and D. Gogny, Phys. Rev. C21 (1980) 1568. [20] M. A. Preston and R. K. Bhaduri, Structure of the Nucleus, (Addison–Wesley, Inc., 1975). [21] N. Anantaraman, H. Toki and G. Bertsch, Nucl. Phys. A398 (1983) 279. [22] G. Bertsch, J. Borysowicz, H. McMagnus and W. G. Love, Nucl. Phys. A284 (1977) 399.

28

[23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49]

G. R. Satchler, Direct Nuclear Reactions, (Oxford University Press, New York, 1983). D. T. Khoa and W. von Oertzen, Phys. Lett. B304 (1993) 8. F. J. Eckle, H. Lenske et al., Phys. Rev. C39 1662 (1989); Nucl. Phys. A506 159 (1990). H. M. Sommermann, K. F. Ratcli?, T. T. S. Kuo, Nucl.Phys A406 (1983) 109. H. A. Bethe, Phys.Rev. 167 (1968) 879. H. S. Koehler, Phys. Rev. 18C (1975) 217. B. D. Day, Rev. Mod. Phys. 39 (1967) 719. F. de Jong and H. Lenske, Phys.Rev. C, in print. G. R. Satchler and W. G. Love, Phys. Rep. 55 (1979) 183. A. M. Kobos, B. A. Brown, R. Lindsay and G. R. Satchler, Nucl. Phys. A425 (1984) 205. S. Haddad and M. Weigel, Phys. Rev. C48 (1993) 2740. G.R. Satchler, Nucl.Phys. A55 (1964) 1. A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems, (McGraw-Hill, New York, 1971). E. Chabanat, Doctorial Thesis, University of Lyon (France) 1995; E. Chabanat, P. Bonche, P. Haensel, J. Meyer, R. Schae?er, in print. H. Kucharek, P. Ring and P. Schuck, Z. Phys. A334 (1989) 119. H. Kucharek, P. Ring, P. Schuck, R. Bengtsson and M. Giro, Phys. Lett. B216 (1989) 249. S. Takahara, N. Onishi and N. Tajima, Phys. Lett. B331 (1994) 261. J. M. C. Chen, J. W. Clark, E. Krotscheck and R. A. Smith, Nucl. Phys. A451 (1986) 509. R. C. Kennedy, Phys. Rev. 144 (1966) 804. A.R. Edmonds, Angular Momentum in Quantum Mechanics, (Princeton University Press, Princeton, 1957). D. Vautherin und D. M. Brink, Phys. Rev. C5 (1972) 626. C. Mahaux and R. Sartor, Nucl. Phys. A493 (1989) 157. P. Grange and M. A. Preston, Nucl. Phys. A219 (1974) 266.; C. J. Horrowitz and B. D. Serot, ibid., A368 (1981) 503, A399 (1983) 529 G. Audi and A.H. Wapstra, Nucl.Phys. A595 (1995) 409. M. Chartier et al., Phys. Rev. Lett. 77 (1996) 2400. D. W. L. Sprung, P. K. Banerjee, Nucl.Phys. A168 (1971) 273. J. W. Negele, Phys.Rev. C1 (1970) 1260.

29

它可用非相对论的*Hartree*-*Fock*自洽场HF-SCF_物理_自然科学_专业资料。它可用非相对论的*Hartree*-*Fock*自洽场HF-SCF 第2章 结合能与化学位移 1. 2. 3. 4. ...

化学学报 ACfA CHIMICA SINICA 1996 , 54 , 979~987 *Hartree* - *Fock* 与密度...*in* *the* Hybrid *Approach* of *Hartree* - *Fock* and *Density* Functional Theoηr ...

Skyrme *Hartree*-*Fock*理论对原子核奇偶差及气泡核结构的研究_电子/电路_工程科技_专业资料。本文基于Skyrme *Hartree*-*Fock*(SHF)理论研究了原子核的奇偶差及气泡核的...

Properties of Hot Nuclear and Neutron Matter *In* a Relativistic *Hartree*-*Fock* Theory_电子/电路_工程科技_专业资料。办声 A to , e e l a r P h y s ...

用*Hartree*-*Fock*法预测取代苯亚硒酸的pK_a值_电子/电路_工程科技_专业资料。第4 4 卷第 1 l期 201 5年1 1月 当 代 化 工 Contemporary Chemical *In* ...

龙源期刊网 http://www.qikan.com.cn 用 *Hartree**Fock* 法预测取代苯亚硒酸的 pKa 值 作者:张梦阳 裴诗恩 苏凌峰 刘喜灵 徐娉婷 刘诗咏 来源:《当代化工》...

II. A 4 *HartreeFock*, MP2, and *Density* ... results from other techniques such as *matrix* ... 94 program was used *in* *the* DFT *calculations*. ...

Dirac-Brueckner-*Hartree*-*Fock* *calculations* for iso...*The* potential *matrix* elements have been adapted ... relativistic Brueckner *approach*, *density* dependent ...

{1.5}$MnO$_4$, which displays spin, charge and orbital ordering, by means of unrestricted *Hartree*-*Fock* *calculations* on *the* multiband $p$-$d$ model....

Although *Hartree*-*Fock* *calculations* often have large...0 is *the* *density* *matrix* *in* real space. A? ...*approach* is used *in* *the* CRYSTAL code for *the* ...

barriers within *the* Skyrme-*Hartree*-*Fock* theory_...*the* *density*-dependent δ -pairing interaction (DDD... they have been obtained *in* *the* *calculations* ...

Slow long-range decay of bound *Hartree*-*Fock* orbitals and enhancement of *the*...*in* *the* *density* functional *approach*. Approximate *calculations* may also lead to...

昌吉 学 院学报 2014 年第 4期 浅析《 计算材料 学》 教学 中的 *Hartree*-*Fock* 方法 以H 2 为例探讨 *H a rt r e e*-*Fo C k* 方法 朱应 涛 ...

第44 卷第 11 期 2015 年 11 月 当代化工 Contemporary Chemical Industry Vol.44,No.11 November,2015 用 *Hartree*-*Fock* 法预测取代苯亚硒酸的 pKa 值张梦阳...

原子结构的*Hartree**Fock*近似的变分法分析_IT/计算机_专业资料。本文用变分法分析了原子结构的*Hartree*-*Fock*的近似,对二电子系统进行了波函数的能量本征态的研究,并...

- Tensor interaction in Hartree--Fock calculations
- CUMULANT EXPANSION METHOD IN THE DENSITY MATRIX APPROACH TO WAVE PACKET DYNAMICS IN MOLECUL
- Point symmetries in the Hartree-Fock approach Symmetry-breaking schemes
- Density-Dependent Relativistic Hartree-Fock Approach
- Symmetry of the Atomic Electron Density in Hartree, Hartree-Fock, and Density Functional Th
- Superheav magic strctures in the relativistic Hartree-Fock-Bogoliubov approach
- Bimodal fission in the Skyrme-Hartree-Fock approach
- Spin-dependent localized Hartree-Fock density-functional approach for the accurate treatmen
- A new approach to the modelling of local defects in crystals the reduced Hartree-Fock case
- Density Matrix Approach to Local Hilbert Space Reduction