Continuum QRPA calculation with the Skyrme effective force

1
Continuum QRPA calculation with the Skyrme
effective force

• Niigata University
• Kazuhito Mizuyama
• Masayuki Matsuo Yasuyoshi Serizawa

2
Mean field theory and Physics of the neutron rich
region
(http//www.rarf.riken.go.jp/RIBF/nuclearchart-e.h
tm)
Halo, Skin, Soft dipole, Pigmy etc

• Shallow fermi level
• Many valence neutrons

Continuum QRPA can treat both the continuum and
pairing.
3
Continuum QRPA
Continuum QRPA
Exact HFB green function
Halo, Skin, Soft dipole, Pigmy etc
satisfies correct asymptotic
outgoing boundary conditions
4
-- Theoretical defects --
Continuum QRPA

• Self consistency is broken.
• Sum rule is not conserved.

Skyrme HFB Landau-Migdal approx. cQRPA

• Serizawa, in the previous talk.
• E. Khan, et.al. PRC66(2002)024309

To improve these defects,
Skyrme HFB Landau-Migdal approx. cQRPA
Continuum QRPA calc. with the velocity dependent
terms of the Skyrme interaction.
5
Extension of continuum QRPA with Skyrme force
p-h interaction
p-p interaction
Linear response equation of continuum QRPA
6
Numerical check
Rbox15fm lcut4 Ecut60MeV
p-h channel
SLy4, (SkM, SGII, SIII, only for 22O)
p-p interaction
Smoothing constant e1.0MeV
V0315MeV fm-3
?0 0.32 fm-3 ( Mix pairing (Dobaczewski et al.
nucl-th/0203060,0203038))
External field
(E1-dipole)
7
Check(1) Renormalization factor
The renormalization factor becomes close to 1 in
the full calculation.
?
?
Consistent result with M. Yamagami, proceedings
of the International Symposium on FRONTIERS OF
COLLECTIVE MOTIONS (CM2002)_at_Aizu, Japan 6-9 Nov.
2002.
Residual interaction with velocity dependent
terms
?
?
8
Check(2)
16O E1 dipole excitation
Smoothing constant e1.0MeV
16O
E dB(E1)/dE e2fm2
SLy4
GDR energy position agree with
SIII
Ex MeV
(cf. K.F.LIU, N.V.GIAI, Phys.Lett.Vol.65,23(1976))
9
Check(3)
20O E1 dipole excitation
Landau-Migdal
1,?,?
20O
lt?ngt 1.93 MeV
dB(E1)/dE e2fm2/MeV
SLy4
Ex MeV
M. Yamagami, proceedings of the International
Symposium on FRONTIERS OF COLLECTIVE
MOTIONS (CM2002)_at_Aizu, Japan 6-9 Nov. 2002.
10
--Summary--

• We formulated the continuum QRPA based on the
  Skyrme energy functional with keeping the
  velocity dependent terms.
• We constructed its new numerical code.
• We checked the renormalization factors for some
  kinds of parameters of the Skyrme interaction in
  oxygen isotopes.
• We also calculated E1 strength functions for
  oxygen isotopes and compared with others work.