Andrey Shirokov (Moscow State Univ.

)
In collaboration with
Alexander Mazur (Pacific National Univ.)
Pieter Maris and James Vary (Iowa State Univ.)

INT, Seattle, June 8, 2011

* Conventional: bound state energies are associated with variational
minimum in shell model, NCSM, etc., calculations
* Is it also true for resonant states? Can we get resonance width from
such calculations?
* Resonant states: should we, probably, study excitation energies
instead? Or the energies E = EA − EA − 1 for n−(A−1) scattering (or,
generally, with A1 + A2 = A)?
* Is it important for them to be stable with respect to ħΩ or Nmax
variation?
* Would be nice to have a simple answer from conventional calculations
without doing, say, complicated NCSM−RGM calculations.
* So, what are the general properties of eigenstates in continuum
consistent with resonance at the energy Er and width Γ?
* Some observations, examples follow; this is still work in progress.
* So, what are the general properties of eigenstates in continuum
consistent with resonance at the energy Er and width Γ?
* Some observations, examples follow; this is still work in progress.
* I discuss some general properties for oscillator basis calculations;
this is the only relavance to NCSM
* O.Rubtsova, V.Kukulin, V.Pomerantsev, JETP Lett. 90,
402 (2009); Phys. Rev. C 81, 064003 (2010):
* I.M.Lifshitz (1947):
0
H =H + V
0
Ej - E j
d (E j ) =- p
Dj
0 0
D j =E j+1 - E j
* O.Rubtsova, V.Kukulin, V.Pomerantsev, JETP Lett. 90,
402 (2009); Phys. Rev. C 81, 064003 (2010):
* I.M.Lifshitz (1947):
0
H =H + V So, the phase shift at the
0 eigenenergies Ej can be
Ej - E
d(E j ) =- p j easily calculated!
Dj
0 0
D j =E j+1 - E j
* O.Rubtsova, V.Kukulin, V.Pomerantsev, JETP Lett. 90,
402 (2009); Phys. Rev. C 81, 064003 (2010):
* I.M.Lifshitz (1947):
0
H =H + V Unfortunately, this does not work:
The dimensionality of the matrix is
0 small, the average spacing
Ej - E j between the levels is not well-
d(E j ) =- p defined.
Dj One needs sometimes Dj value
below the lowest Ej0
0 0
D j =E j+1 - E j
Direct and
inverse problem
 N

å
n '=0
l
H nn' án' | l ñ=El n l , n £N
2
N
N l
GNN (E) =- å
l =0 El - E
l
SNl (q) - GNN (E)TN,N +1S N +1,l (q) SN +1,l (ql )
tan d =- l
Þ - = f (El )
CNl (q) - GNN (E)TN,N +1CN +1,l (q) E® E l CN +1,l (ql )
 arctan(SN+1,l / CN+1,l)
N+1=10, l=0
1980

1800

1620

1440

1260

1080
degrees

900

720

540

360

180

0
0 5 10 15 20
Ecm (hW )
 arctan(SN+1,l / CN+1,l)
N+1=10, l=0
1980

1800

1620

0
1440
H =H + V
1260
0
1080 Ej - E
degrees

j
900 d (E j ) =- p
720
Dj
D j =E 0j+1 - E 0j
540

360

180

0
0 5 10 15 20
Ecm (hW )
 0
H =H + V
0
Ej - E
d(E j ) =- p j

Dj
D j =E 0j+1 - E 0j
 IS based on WS (ver. 3) phase shift
N+1=7, l=2, hW =15 MeV
180
phase shifts (degrees)

90

explicit phase shift (Numerov)

J-matrix phase shift
f(E)
shifted (0.8 MeV) phase shift

0
0 5 10
Ecm (MeV)
*J-matrix inverse N-nucleus scattering analysis suggests
values for resonant and non-resonant states that should
be compared with that obtained in NCSM
 N

å l
H nn' án' | l ñ=El n l , n £N
n '=0
2
N
N l
GNN (E) =- å
l =0 El - E
N

å
2
N l =1
l =0
l
SNl (q) - GNN (E)TN,N +1SN +1,l (q)
tan d =- l
CNl (q) - GNN (E)TN,N +1C N +1,l (q)

2E
q2 =
W
 * J-matrix:
Let us try to extract resonance information from
Eλ behavior only
N

å
n'=0
l
H nn' án' | l ñ=El n l , n £N
2
N
Nl
GNN (E) =- å
l =0 El - E
l
SNl (q) - GNN (E)TN,N +1S N +1,l (q) SN +1,l (ql )
tan d =- l
Þ - = f (El )
CNl (q) - GNN (E)TN,N +1CN +1,l (q) E® El CN +1,l (ql )
 IS based on WS (ver. 3) phase shift
N+1=4, l=2, D E = 4 MeV
180
hW =10 MeV
15 MeV
20 MeV
25 MeV
30 MeV

phase shifts (degrees)

90

Eλ should increase with ħΩ

0
5 6 7
Ecm (MeV)
Within narrow resonance Eλ is
nearly ħΩ-independent

The slope of Eλ(ħΩ) depends

however on Nmax, l, Eλ value
 * Breit-Wigner: BW G /2
d =arctan +j
Er - E
SN +1,l (ql )
tan d(El ) =-
CN +1,l (ql )
Er =El + D l

* Simple approximation:
SN +1,l (qlφ=0
) G /2
tan d (El ) =- =
CN +1,l (ql ) D
d tan d (E) G/2
= 2
dE E=El D

Derivatives calculated through El (W ), E l±(W ±D (W ))

 * Breit-Wigner: G /2
d BW
=arctan +j
Er - E
SN +1,l (ql )
tan d(El ) =-
CN +1,l (ql )
Er =El + D l

Do not expect to get a reasonable

* Simple approximation:
SN +1,l (ql ) Gφ=0
/2
tan d(El ) =- = result for Er or Γ if Γ/2Δ is small!
CN +1,l (ql ) D
If |Γ/2Δ| is large, we get good results
for Er, Γ and φ.
d tan d (E) G/2
= 2
dE E =El D

Derivatives calculated through El (W ), E l±(W ±D(W ))

 WS, version 3
N+1=10, l=2, hW =20 MeV
180
phase shifts (degrees)

90

explicit phase shift (Numerov)

BW with explicit E r, G (from S-pole)
f(E) WS, version 3
J-matrix N+1=10, l=2, hW =20 MeV
BW, meth. III 8

BW with explicit E r, G
7
J-matrix
0 BW, meth.III
0 5 10 6
Ecm (MeV)

total cross section (barn)

5

0
1 2 3 4
Ecm (MeV)
 WS, version 3
N+1=10, l=2, hW =30 MeV
180
phase shifts (degrees)

90

explicit phase shift (Numerov)

BW with explicit E r, G (from S-pole)
f(E)
J-matrix phase shift
BW, meth.III WS, version 3
N+1=10, l=2, hW =30 MeV
8
0 BW with explicit E r, G
0 5 10 7
Ecm (MeV) J-matrix
BW, meth.III
6

total cross section (barn)

5

0
1 2 3 4
Ecm (MeV)
 IS based on WS (ver. 3) phase shift
N+1=7, l=2, without shift
180
phase shifts (degrees)

90

explicit ph. shift

hW =15 MeV
25 MeV IS based on WS (vers. 3)
BW, hW =15 MeV
BW, hW =25 MeV N+1=7, l=2, without shift
7

0 Inv. scat.
1 2 3 4 5
Ecm (MeV) 6 BW, hW = 15 MeV
BW, hW = 25 MeV

total cross section (barn)

4

0
1 2 3 4
Ecm (MeV)
 Eigenvalues E 0 and resonance energies E r
IS based on WS ver.3 (N+1=7, l=2, without shift)

3
E0, Er (MeV)

E0
Er [D (hW ) = 0.5 MeV]
Er [D(hW ) = 2.5 MeV]
Inv. scat. (S-matrix)
2.5 Widths
IS based on WS ver.3 (N+1=7, l=2, without shift)
0.8

10 15 20 25 30
hW (MeV)
0.7

0.6
G (MeV)

0.5 G [D (hW ) = 0.5 MeV]

G [D(hW ) = 2.5 MeV]
Inv. scat. (S-matrix)

0.4

0.3
10 15 20 25 30
hW (MeV)
* What can we do if we obtain Eλ in a non-
resonant region above the resonance?
* We can extrapolate energies to larger (finite)
Nmax value when Eλ is in the resonant region.
* Expected dependence is
This works. However this
extrapolation seems to
be unstable and
inconvenient
 More convenient is an
exponential
extrapolation.

El (N max ) =a exp(- cN max ) + El (¥ )

* We get stable
Er and Γ; Γ is
too small as
compared
with
experiment.
* I discussed general features of continuum states
obtained in many-body calculations with oscillator basis.
* The best way to compare the calculated results with
experiment is to use “experimental” phase shifts and get
Eλ consistent with scattering data using simple inverse
scattering technique.
* Studying ħΩ dependence of Eλ obtained in NCSM, one can
get resonance energy and width. However, usually an
extrapolation to a reasonable Nmax value is required.