Extended Brueckner-Hartree-Fock theory in many body system - Importance of pion in nuclei -
description
Transcript of Extended Brueckner-Hartree-Fock theory in many body system - Importance of pion in nuclei -
11.7.22 toki@genkenpion 1
Extended Brueckner-Hartree-Fock theory in many body system
- Importance of pion in nuclei -
Hiroshi Toki (RCNP, KEK)In collaboration with
Yoko Ogawa (RCNP)
Jinniu Hu (RCNP)
Kaori Horii (RCNP)
Takayuki Myo (Osaka Inst. of Technology)
Kiyomi Ikeda (RIKEN)
11.7.22 toki@genkenpion 2
Pion is important in Nuclear Physics !
• Yukawa (1934) predicted pion as a mediator of nuclear interaction to form nucleus
• Meyer-Jansen (1949) introduced shell model ー beginning of Nuclear Physics
• Nambu (1960) introduced the chiral symmetry and its breaking produced mass and the pion as pseudo-scalar particle
11.7.22 toki@genkenpion 3
Shell model (Meyer-Jensen)
• Phenomenological
• Strong spin-orbit interaction added by hand
• Magic number
• 2,8,20,28,50,82
QuickTime˛ Ç∆TIFFÅiîÒà≥èkÅj êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
2
8
20
40
70
Harmonicoscillator
11.7.22 toki@genkenpion 4
The importance of pion is clear in deuteron
QuickTime˛ Ç∆TIFFÅiîÒà≥èkÅj êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
QuickTime˛ Ç∆TIFFÅiîÒà≥èkÅj êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
Deuteron (1+)QuickTime˛ Ç∆
TIFFÅiîÒà≥èkÅj êLí£ÉvÉçÉOÉâÉÄǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
NN interaction S=1 and L=0 or 2
€
π
QuickTime˛ Ç∆TIFFÅiîÒà≥èkÅj êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
11.7.22 toki@genkenpion 5
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
11.7.22 toki@genkenpion 6
Variational calculation of few body system with NN interaction
C. Pieper and R. B. Wiringa, Annu. Rev. Nucl. Part. Sci.51(2001)
€
ΨVπ Ψ
Ψ VNN Ψ~ 80%
€
Ψ=φ(r12)φ(r23)...φ(rij )
VMC+GFMC
VNNN
Fujita-Miyazawa
Relativistic
Pion is keyHeavy nuclei (Super model)
11.7.22 toki@genkenpion 7
Pion is important in nucleus
• 80% of attraction is due to pion
• Tensor interaction is particularly important
€
rσ 1 ⋅
r q
r σ 2 ⋅
r q =
1
3q2S12( ˆ q ) +
1
3
r σ 1 ⋅
r σ 2q
2
Pion Tensor spin-spin
€
S12( ˆ q ) = 24π Y2( ˆ q ) σ 1σ 2[ ]2[ ]0
11.7.22 toki@genkenpion 8
Halo structure in 11LiQuickTime˛ Ç∆
TIFFÅiîÒà≥èkÅj êLí£ÉvÉçÉOÉâÉÄǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
QuickTime˛ Ç∆TIFFÅiîÒà≥èkÅj êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
QuickTime˛ Ç∆TIFFÅiîÒà≥èkÅj êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
Deuteron structure in shell model is producedby 2p-2h states
QuickTime˛ Ç∆TIFFÅiîÒà≥èkÅj êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
Deuteron wave function
Myo Kato Toki Ikeda PRC(2008)
(Tanihata..)
11.7.22 toki@genkenpion 9Pairing-blocking : K.Kato,T.Yamada,K.Ikeda,PTP101(‘99)119, Masui,S.Aoyama,TM,K.Kato,K.Ikeda,NPA673('00)207. TM,S.Aoyama,K.Kato,K.Ikeda,PTP108('02)133, H.Sagawa,B.A.Brown,H.Esbensen,PLB309('93)1.
11.7.22 toki@genkenpion 10
Tensor optimized shell model ( TOSM )
Myo, Toki, Ikeda, Kato, Sugimoto, PTP 117 (2006)
0p-0h + 2p-2h
Energy variation
(size parameter)
Gc
11.7.22 toki@genkenpion 11
11Li G.S. properties (S2n=0.31 MeV)
Tensor+Pairing
Simon et al.
P(s2)
Rm
E(s2)-E(p2) 2.1 1.4 0.5 -0.1 [MeV]
Pairing interaction couples (0p)2 and (1s)2 states.
11.7.22 toki@genkenpion 12
Unitary Correlation Operator Method
H. Feldmeier, T. Neff, R. Roth, J. Schnack, NPA632(1998)61
corr. uncorr. SM, HF, FMDCΨ = ⋅Φ ←
{ }12exp( ), ( ) ( )ij r r
i j
C i g g p s r s r p<
= − = +∑
short-range correlator
( ( ))( )
( )
s R rR r
s r+
+′ =
† 1 (Unitary trans.)C C−=
rp p pΩ= +r r r
Bare Hamiltonian
† : Hermitian generatorg g=
Shift operator depending on the relative distance r
€
C+HCΦ = EΦ
€
HΨ = EΨ
(UCOM)
† ( )C r C R r+=
11.7.22 toki@genkenpion 13
4He with UCOM
He : (0 ) with -parameterA As b
11.7.22 toki@genkenpion 14
TOSM+UCOM with AV8’
T
VLS
E
VC
VT
Few bodyCalculation
(Kamda et al)
€
Ψ =C0 0 + Cα 2 p2h :αα∑
(Myo Toki Ikeda)
11.7.22 toki@genkenpion 15
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
Y0 Y2+
€
D S12 S ≠ 0
11.7.22 toki@genkenpion 16
TOSM should be used for nuclear many body problem2p-2h excitation is essential for treatment of pion
G.S.
Spin-saturated
The spin flipped states are alreadyoccupied by other nucleons.
Pauli forbidden
€
σ1σ 2[ ]2⋅Y2(r)
11.7.22 toki@genkenpion 17
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅBQuickTime˛ Ç∆
êLí£ÉvÉçÉOÉâÉÄǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB (2011)
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅBQuickTime˛ Ç∆
êLí£ÉvÉçÉOÉâÉÄǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
HF state cannot handle the tensor interaction
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
€
δΨ H Ψ
Ψ Ψ= 0
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
11.7.22 toki@genkenpion 18
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
Total energy
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
Energy variation
11.7.22 toki@genkenpion 19
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
Variation with HF state
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
We can solve finite nuclei by solving the above equation.
11.7.22 toki@genkenpion 20
EBHF equation
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
€
∂∂ψ b
* (x)0 Veff 0
€
Tψ b (x) + d3x'ψ d* (x')V (x,x ') ψ bψ d[ ]A∫
d∑ + C0
2 ∂ 0 ˆ V αμ
∂ψ b* (x)
1
E − βν ˆ H αμαμ ,βν∑ βν ˆ V 0
+ C0
20 ˆ V αμ
αμ ,βν ,α 'μ ',β 'ν '∑
1
E − α 'μ ' ˆ H αμ
∂ α 'μ ' ˜ H β 'ν '
∂ψ b* (x)
1
E − βν ˆ H β 'ν 'βν ˆ V 0 = ε bψ b (x)
Similar to BHF equation
11.7.22 toki@genkenpion 21
Feshbach projection method
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
€
P = 0 0
€
Q = 2p2h :αα∑ 2p2h :α
€
Veff = C0
2V + C0
20
αβ∑ V α
1
E − β H αβ V 0
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
11.7.22 toki@genkenpion 22
Comparison to BHF theory
€
Veff = C0
2V + C0
20
αβ∑ V α
1
E − β H αβ V 0
= C0
2V − 0
α∑ V α
1
Eα
α V 0 + 0αβ∑ V α
1
Eα
α V β1
Eβ
β V 0 + .. ⎡
⎣ ⎢
⎤
⎦ ⎥
€
E − β H α ≈ −Eα δαβ − β V α
•There appears |C0|2.•|C0|2 is not normalized to 1.
€
C0
2+ Cα
α∑
2= 1
€
G = V − VQ
eG = V − V
Q
eV + V
Q
eV
Q
eV − ...
11.7.22 toki@genkenpion 23
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
QuickTime˛ Ç∆ êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
Nucl. Phys (2011)
Direct terms only
11.7.22 toki@genkenpion 24
Relativistic chiral mean field model
The difference between 12C and 16O is 1.5 MeV/N.
The difference comes from low pion spin states (J<3).This is the Pauli blocking effect.
P3/2
P1/2
C
O
S1/2
Pion energy Pion tensor provides large attraction for 12C
OC
Pion contribution Individual contribution
O
C
11.7.22 toki@genkenpion 25
Renaissance in Nuclear Physics by pion
• We have a EBHF theory with pion• Unification of hadron and nuclear physics• Pion is the main player ー Pion renaissance• High momentum components are produced by
pion ( Tensor interaction )ー Nuclear structure renaissance
• Nuclear Physics is truly interesting!!
11.7.22 toki@genkenpion 26
QuickTime˛ Ç∆TIFFÅiîÒà≥èkÅj êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
Monden
QuickTime˛ Ç∆TIFFÅiîÒà≥èkÅj êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
11.7.22 toki@genkenpion 27
Relativistic Brueckner-Hartree-Fock theoryBrockmann-Machleidt (1990)
€
G = V + VQ
eG
Us~ -400MeVUv~ 350MeV
€
π
€
πQuickTime˛ Ç∆
TIFFÅiîÒà≥èkÅj êLí£ÉvÉçÉOÉâÉÄǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
relativity
RBHF
Non-RBHF
€
ρ
11.7.22 toki@genkenpion 28
Important experimental data
11.7.22 toki@genkenpion 29
Deeply bound pionic atom
Toki Yamazaki, PL(1988)
Prediction
Found by (d,3He) @ GSIItahashi, Hayano, Yamazaki..Z. Phys.(1996), PRL(2004)
Physics : isovector s-wave
€
b1
b1(ρ )=1− 0.37
ρ
ρ 0
€
fπ2mπ
2 = −2mq ψ ψ
€
ψ ψρ
ψ ψ=1− 0.37
ρ
ρ 0
€
b1 ∝1
fπ2
11.7.22 toki@genkenpion 30
QuickTime˛ Ç∆TIFFÅiîÒà≥èkÅj êLí£ÉvÉçÉOÉâÉÄ
ǙDZÇÃÉsÉNÉ`ÉÉÇ å©ÇÈÇΩÇflÇ…ÇÕïKóvÇ≈Ç∑ÅB
Suzuki, Hayano, Yamazaki..PRL(2004)
Optical model analysisfor the deeply bound state.
11.7.22 toki@genkenpion 31
Nuclear structure caused by pion
• 2 p− 2 h excitation is 20%• High momentum components• Low momentum components
( Shell model ) are reduced
by 20 %
11.7.22 toki@genkenpion 32
RCNP experiment (high resolution)
Y. Fujita et al.,E.Phys.J A13 (2002) 411
H. Fujita et al., PRC
€
Ψf στ Ψ i
Not simpleGiant GT
11.7.22 toki@genkenpion 33
16O (p,d) Ep = 198 MeVΘd = 10°
16O (p,d) Ep = 295 MeVΘd = 10°
16O (p,d) Ep = 392 MeVΘd = 10°
15OLevel scheme
1/2
-0.0
- -+ +Ong, Tanihata et al
11.7.22 toki@genkenpion 34
Relative Cross Section
a) J. L. Snelgrove et al., PR187, 1246(1969) b) J.K.P. Lee et al., NP A106, 357(1968) c) G.R. Smith et al., PRC30, 593(1984)
Ex[MeV] Jπ
gs 0.0 1/2-
1 5.2 1/2+, 5/2+ 2 6.2 3/2-
3 6.8 3/2+, 5/2+
4 7.3 7/2+
5 7.6 1/2+
6 8.3 3/2+
7 8.9 5/2+,1/2+,(1/2)- 8 9.5 (3/2)+,5/2- 9 10.5 (9/2+),(3/2-),(3/2)+
10 11.6 5/2-
11 11.9 5/2-, 5/2-
• 13.6 5/2+
13.7 3/2-
1 23
45
67 8 9 10
11
12gs
c)a)
Incident proton energy [MeV]
rela
tive
cros
s se
ctio
n
σ1 /σgs
σ2 /σgs
σ4
/σgsσ6
/σgsσ12 /σgs
This workb)
11.7.22 toki@genkenpion 35
Centrifugal potential ([email protected]) pushes away the L=2 wave function
The property of tensor interaction