- 25 -
Symmetry Projected Density Functional Theory and Neutron Halo’s E.C. Lopes1 and P. Ring1 Physik Department, Technische Universit¨ at M¨ unchen, 85748 Garching, Germany we show neutron density distributions of Ne-isotopes as a function of the radius r on a logarithmic scale. Number projected results (full lines) are compared with those without number projection (dashed line). In both cases we find halo phenomena and there is only a small difference in the density distributions. Going from 30 Ne to 40 Ne we find a phase transition from zero pairing in 30 Ne to rather strong pairing in 34 Ne and all the heavier isotopes. Because of vanishing pairing there is no difference between projected and unprojected calculations in 30 Ne. 32 Ne is a transitional nucleus with week pairing. In this case we find the largest differences between the projected and the unprojected calculations. For heavier nuclei with strong pairing, simple Hartree-Bogoliubov theory without projection provides already a rather good approximation and therefore number projection is not necessary and the conclusions found in the earlier calculations are valid.
0
10
30 -2
Ne
-4
Neutronen
10 10
32
34
Ne
Neutronen
Ne
Neutronen
-6
10
Protonen
Protonen
Protonen
-3
The appearance of halo phenomena near the drip line nuclei has challenged our traditional understanding of the nuclei as an incompressible charged liquid drop and extended nuclear physics to low density and inhomogeneous system, where the coupling to the continuum has to be treated in a consistent way. Recently Relativistic Hartree Bogoliubov (RHB) theory in the continuum has been applied successfully to the description of halo phenomena in light and medium heavy nuclei [1, 2, 3]. This theory provides a self-consistent treatment of pairing correlation in the presence of the continuum and allows a microscopic description of halo phenomena in the framework of density functional theory. Essential conditions for the formation of a neutron halo have been found: (a) the Fermi surface of the neutrons has to be close to the continuum limit, (b) there has to be a large level density in the vicinity of the Fermi surface, such that pairing correlations can develop and (c) there have to be single particle states with zero or small angular momentum barrier in close neighborhood. This is the case in 11 Li [1], where the 2s1/2 level is shifted downwards and where a neutron halo as been observed, but also in Ne-nuclei [2] close to the neutron drip line, where 2p3/2 and 2p1/2 reach the continuum limit and where a neutron halo has been predicted. The same theory also predicts further halo phenomena in the Zr-region [3] with multi-particle halos containing up to 6 neutrons. One of the essential shortcoming of mean field theories is the fact, that they violate symmetries. In HartreeBogoliubov (HFB) approximation gauge symmetry of particle number conservation is broken and the wave functions are linear combinations of wave functions with different particles. In principle one has to project onto the eigenspace of good particle number and this projection has to be carried out before the variation. In heavy nuclei, where many particles contribute to the pairing correlations one has found that an explicit projection is not necessary and the BCS-approach, where the particle number is only conserved on the average presents already a very reasonable approximation. However, in light nuclei it is a priori not clear, whether such arguments are not applicable and whether the simple Hartree-Bogolubiov theory is applicable. Recently a new method has been developed to solve the number projected Hartree-Fock-Bogoliubov problem [4]. So far this method has only been applied to simple non-relativistic models. We have now implemented this method for realistic RHB theory. For this purpose a new code has been developed to solve the number projected relativistic Hartree-Bogoliubov equations in by a variation after projection. This method has been applied in several areas of the periodic system [5]. Here we show as an example calculations for the chain of Ne-isotopes, where in Ref. [2] extended halo phenomena have been found. In Fig. 1
ρ(r) [fm ]
1
-8
10 0 10
36
38
Ne
-2
10
-4
Neutronen
10
40
Ne
Neutronen
Ne
Neutronen
-6
10
Protonen
Protonen
Protonen
-8
10
0
5
10
15
0
5
10
15
0
5
10
15
r [fm]
Figure 1: The neutron density distributions as a unction of radius r for 30 Ne, 32 Ne, 34 Ne 36 Ne, 38 Ne, and 40 Ne in Projected RHB theory (solid line) compared with that in the case without projection (dashed line).
Concluding we have implemented full number projection before variation in relativistic Hartree-Bogoliubov theory. This will allow many other applications in transitional nuclei in the future.
References [1] J. Meng and P. Ring, Phys. Rev. Lett. 77, 3963(1996). [2] W. P¨ oschl, D. Vretenar, G.A. Lalazissis, and P. Ring, Phys. Rev. Lett. 79, 3841 (1997). [3] J. Meng and P. Ring, Phys. Rev. Lett. 80, 460(1998). [4] J.A. Sheikh and P. Ring, Nucl. Phys. A665, 71 (2000). [5] E. C. Lopes, PhD thesis, Technische Universit¨ at M¨ unchen (2002)
- 26 -
Renormalization group approach to neutron matter Achim Schwenka , Bengt Frimanb and Gerald E. Brownc Department of Physics, The Ohio State University, Columbus; b GSI, Darmstadt; Astronomy, SUNY Stony Brook
The main idea of the renormalization group (RG) approach to Fermi liquids [1] is to “adiabatically” include the in-medium corrections to the effective interaction by solving the RG equations in the relevant channels. In this exploratory calculation [2], we include the particle-hole (ph) channels, which play a special role in Fermi liquid theory. The main effect of scattering in the particle-particle channel, the taming of the short-range repulsion, is taken care of by using Vlow k [3] as the starting point of the RG flow. The low-lying excitations in this channel, which are responsible, e.g., for superfluidity, are not included. These are then treated explicitly when we compute the superfluid gap by employing BCS theory for the fully reducible scattering amplitude. We derive the one-loop RG equations for the quasiparticle interaction and the scattering amplitude at zero temperature. The evolution of the effective mass is included in the RG flow, as well as a simplified treatment of the renormalization of the quasiparticle strength. As we decimate down to the Fermi surface, by taking the cutoff Λ → 0, we obtain not only the forward scattering amplitude for low-lying quasiparticle-quasihole excitations, but also the amplitude for general (non-forward) scattering processes of quasiparticles on the Fermi surface. The one-loop RG equations in the ph channels are solved at zero temperature. The corresponding diagrams, which renormalize the four-point vertex, are shown in Fig. 1. Here the momenta pi of the intermediate ph pair lie in a shell Λ − dΛ ≤ |pi − kF | ≤ Λ. When only one channel is considered, the one-loop RG equation is exact in the sense that it is equivalent to the corresponding scattering equation, and when both the direct and exchange ph channels are included, the scattering amplitude remains antisymmetric under the RG flow. The RG approach includes the induced interaction. We work in the approximation that both ph momentum transfers are small compared to the Fermi momentum. The physical system considered is neutron matter, where complications of the tensor force do not enter in the S-wave.
Figure 1: The one-loop contributions to the RG equation. γ denotes the running effective four-point vertex.
∆ [MeV]
a
c
Department of Physics and
3.5 3
direct
2.5 2 1.5 1
full
0.5 0 0
0.2 0.4 0.6 0.8
1
1.2 1.4 1.6 -1
k F [fm ]
Figure 2: The 1 S0 superfluid gap versus the Fermi momentum kF . The dots denote the pairing obtained from Vlow k only, whereas the squares and the triangles are computed from the full RG solution using different approximations for the quasiparticle strength [2]. In comparison, the dashed line is obtained by solving the BCS equation with the bare NN interaction [5], while the solid line includes the ph polarization effects [4]. Using the resulting Fermi liquid parameters, we compute the 1 S0 superfluid pairing gap in weak coupling BCS theory, Fig. 2. This application probes the angular dependence of the scattering amplitude. We generally find good agreement with the results obtained in the polarization potential model by Wambach et al. [4]. The RG method is a promising tool for studying a wide range of nuclear many-body problems. For a similar analysis of symmetric or asymmetric nuclear matter, it is necessary to extend the flow equations to incorporate tensor interactions. Fairly large renormalization effects are expected and spin non-conserving interactions are generated in the medium [6]. In particular, polarization effects on the 3 P2 –3 F2 pairing of neutrons, where the tensor and spin-orbit forces play a crucial role, may have important effects on neutron star properties.
References [1] R. Shankar, Rev. Mod. Phys. 66 (1994) 129. [2] A. Schwenk, B. Friman and G.E. Brown, Nucl. Phys. A713 (2003) 191. [3] S.K. Bogner, T.T.S. Kuo and L. Coraggio, Nucl. Phys. A684 (2001) 432c; S.K. Bogner, et al., nuclth/0108041. [4] J. Wambach, T.L. Ainsworth and D. Pines, Nucl. Phys. A555 (1993) 128. [5] U. Lombardo and H.-J. Schulze, astro-ph/0012209. [6] A. Schwenk and B. Friman, to be published.
- 27 -
Effective NN-Interaction Including Radial and Tensor Correlations T. Neff and H. Feldmeier
5 PSfrag4 replacements
2 −1
[fm ]
1 0
1
and the unitary tensor correlator ½ X ³ ´¾ 3 C = exp −i ϑ(r ) (~ σ ~ r )(~ σ p ~ )+h.a. . (3) Ω ij i ij j Ωij ∼ ∼ ∼ ∼ ∼ 2 ∼ i=
[email protected]3-9ER9,;p¦fx©ea
[ jbw¢spq c Ô Ëíì î Í ¿ ½¡Á4ïïï à,áâ car]`a
pq[]¥ º¦g©Yo_oo[4o
os_bjipq[^] [4ejip¦_ipq[^]pqj [L_iac^jbpq_bpq[]£ Z\[» ac£vµKo]`mepT2_imeoµo[eojiqw¤cpC]b]µ¼½C¾¿?À¦½CÁµa
[»de[eac_i]7p¦f knmeojio]bdec_ba
[ej _bo]{H]is co~«©¡w_imeo_bï ðij ñÕa
}Cò a KÂlÃt\xoYo[xo[^soq{ ÂÃÄÅÆHÇÀx²ÈºqÉVÊtà à ´ fde[^s_ba
pq[pqfY_imeola
pq[}Lo
pxsac_w ¿ ac[d^[ea°_t] p¦f^_imeon_bm^oji ËxÌ¥}qÍocÎÈpxÏ sac_w Í ÎÈÏ pqf_imeogo
os_bjipq[^] de[^eojOÍ Zkgllu §^ jioo®
4]`mepqpTdegj[7jbo~ac[]`d^ °_tac] £e qji o]o7 óCdefdeac_b[^o¥s!_iacpLo[7
Kpqgf acÂ,_bm³£0ªo jba°_i
acmeopqj&dx/_px x oqr]¥[eo _bôTar?s £ ±dx_9a
[/_bmeolejbo~]`o[^sogpqf ]_ijbpL[e® q[eo_barsn§^o
_bm^ol/pq©ea
a°_w p¦f%_imeo«
a
qmC_z j`_ia
sco~]a
]7jioxd^so)£õ v] spL[^]`o~óCdeo[^so_bmeo sp¡pLca
[e_iac/oac[^sjbo ]bo]zde_ip [pqjtxoj7p¦f, L[eac_bd^xoµfpqj e m a
qme
wstm jiqo~7a
pq[]£ ö7knmear]lpqjiyzm ]©Yoo[]`de^pLj`_io©¡w α = 30 ÷%ø ZnspLc ©Ypqj _iacpL[zspL[L_ij s!_~£ −(dEi/ds)
0.02
ο ο ο
0.01
0
0
5
10
15
20
i
0.6
−(dEi/ds)
o
proton anti−proton
0.4
B=0 Poth ο B=6T, α=30 ο B=6T, α=60
10
0
5
10
a
qdejio spq[xa°_iacpL[^],²AÐar]g_bm^o [eLco%©Yo_oo[_bmeo, q[eo_barsl§or [^ _bmeoµacpL[}qo
pxsac_we´!£Ñl ]7oxYos!_io³_imeo«o[eoji qw¤cpC]b]za
]z& d^stm jbLojKfpqj [C_iacejip¦_ipq[^]K_im [¥fpqjg^jbpq_bpq[]£ ªo r]bps rsde _io¤_bm^oo[eojbLwcpC]b]fpqjmea
qm^cw³stm jiqo~ acpL[^]{C]`m^ pTg[a
[ ac £CvfpLj9|lÒiÃbÓ£0knm^osp¡pq
ac[^®_ba
olÔp¦fV ]`gacfÕ_ acpL[^]ggac_bm }qocpxsa°_w Í ¿Ö ×ÙØC×ÛÚvÜ Í ÎÈÏ s [©o eejipÝxac _boz©¡w ac[C_boqj _iac[e_imeo¥o[eojbLwµ
pL]i]¼½C¾¿\ÀT½¡ÁÞ Í^¿ ß Ã xpTg[_bp_imeo acpL[7}Lo
p¡sa°_w Í ¿ÈÖ à,áâ gmeojbo¼½C¾ ¿ À¦½CÁ%jbo stmeo~]Kac_i]o y}
deo Ô2ã Ô Ë äKå Í ¿È¿ÈÖ Öà,×ÙØCá×ÛÚâ0ævç ¼è>é4ê Í epqj spq[^]boji} _iac}Loo]`_ba
_boz_bmeostm j s_bojbar]_ia
sej j _bo Ô Ë ß)ë ar]® o}
d _bo _ Í ¿ÈÖ à,áâ gmeojbo_ime oxj jbo~xd^s_ba
pq[x deo vi
2
0
1/(Z τ0) in s
−1
0.2
1
0.1 0
20
40
60
80
100
Z
¯a
qd^jbo®
ùûúVütúVýeúHþvÿ>ú c ¸£el/pqjio_`_iaVo_ A£
{e _bdejio²ALL´Oº¡ô ^£ Ý ÷ £ ÷ ©ejiaco
]boOo_ ?£
{Ýu9m¡w¡]£go}>£ Ho_`_~£ ²? Lq´)e qºeL£ T ªÑ£ dea
[C_o_ ?£
{^w¡V£xZ\[C_£ ²A eÝ´Oº¡ô x£ º¦ ÷ £ Yga
sty¡[ qoA{v£k¯p¡o >oj~{nac [ !"$#%'&)(*,+-*.(/102( +435/16871/:9;¯{^ZuKpL[xf£^uOjipxs¦£xp=£ < >