tsunoda _1

Size: px

Start display at page:

Download "tsunoda _1"

れんかありの
4 years ago
Views:

1 モンテカルロ殻模型によるベータ崩壊の研究東京大学原子核科学研究センター (CNS) 角田佑介

2 モンテカルロ殻模型による研究対象の核種 Hg r プロセス Sn Sm Zr モンテカルロ殻模型 (MCSM) により従来の殻模型計算の手法では困難な核種も計算可能に理研仁科センターアーカイブより引用

原子核の形状 : spherical oblate prolate 一部の核種では異なる形状を持つ固有状態が近いエネルギーに現れる ( 変形共存 ) Motivation 66 変形共存の領域 96 Zr Hg Sm Sn oblate Q(fm 2 ) 0 + 2 0 + 1 0+ 3 spherical 0 + 4 Q(fm 2

3 原子核の形状 : spherical oblate prolate 一部の核種では異なる形状を持つ固有状態が近いエネルギーに現れる ( 変形共存 ) Motivation 66 変形共存の領域 96 Zr Hg Sm Sn oblate Q(fm 2 ) spherical Q(fm 2 ) prolate S. Leoni et al., PRL 118, (2017) Heyde et al., Rev. Mod. Phys. 83, 1467 (2011) 陽子中性子の一方が閉殻に近い領域で変形共存が見られる殻構造との関係 ( 陽子数 Z=28) 同位体近傍の原子核形状とベータ崩壊の関係を主に議論

4 Zr (Z=40) Zr 同位体 (Togashi et al.) プレスリリース原子核での形の量子相転移とスーパーコンピュータ京によるシミュレーション 2016/10/18 Zr Sn Sm Hg プレスリリースより引用

5 Sn (Z=50) Sn 同位体 (Togashi et al.) プレスリリース原子核形状の 2 次相転移をスパコンシミュレーションで発見 2018/08/11 Zr Sn Sm Hg プレスリリースより引用

6 Sm (Z=62) Sm 同位体 (YT et al.) 中性子数 N=82( 魔法数 ) から中性子が増えると球形から変形度が大きくなる Zr Sn Sm Hg Sm 励起エネルギー ( 実験値 ) 4 + : 青 2 + : 赤

7 Hg (Z=80) Hg 同位体 (YT et al.) プレスリリース水銀原子核はハムレット 2018/10/02 実験との共著論文が Nature Physics 誌に掲載 B. A. Marsh et al., Nature Physics 14, 1163 (2018) Sn Zr Sm Hg 1 次相転移が 6 回連続して起きていると考えられるプレスリリースより引用

8 原子核の量子構造に関する殻模型の概要原子核の中の陽子や中性子は入れ物 ( ポテンシャル ) の中で軌道運動量子論ではとびとびのエネルギー ( 図中で水平方向の線 ) それらの軌道に下から順番に陽子や中性子を入れていく上の方の軌道では異なったパターンで入ったものが混在する ( 図では 1 パターンのみが示されている ) 中性子各軌道に入れる粒子数には上限がありその結果ジルコニウム 90 では左の図のようになっている陽子シェード部分は完全に詰まっている 90 Zr の場合記者会見の発表スライドより引用

9 モンテカルロ殻模型 (MCSM) 重い核種を計算したり原子核の変形を十分に扱うには広い模型空間が必要だが非常に広い模型空間では直接対角化は不可能モンテカルロ殻模型では MCSM 基底による小さなハミルトニアン行列を対角化 T. Otsuka et al., PPNP47, 319 (2001) 補助場 MC+ 変分的手法で固有エネルギーを最小化固有状態 Slater 行列式 1 基底ではHFに相当角運動量パリティ射影 MCSM 基底から原子核形状の情報を抽出 (T-plot)

10 Monte Carlo shell model (MCSM) Advanced MCSM dimension Conventional We can perform MCSM calculations in large model spaces. The largest case corresponds to ~10 42 Dimension of Hamiltonian matrix for many-body states to be diagonalized in the conventional shell-model calculation year

重要さ射影後の基底と波動関数との overlap probability T-plot of 0 + 1 state of 68

11 MCSM による原子核形状の解析 (T-plot) 殻模型計算の相互作用を用いて Constrained HF により Potential energy surface (PES) を計算点の位置 : 変形度射影前の MCSM 基底の四重極変形点の面積 : 重要さ射影後の基底と波動関数との overlap probability T-plot of state of 68 (Z=28, N=40) oblate triaxial 角運動量パリティ射影 spherical prolate MCSM 波動関数 MCSM 基底 Slater 行列式

3/2 0f 7/2 core: occupied 68 This model space is wide enough to discuss how magic numbers 28,

12 Monte Carlo Shell model (MCSM) calculation for isotopes with a focus on their shapes unoccupied pfg9d5 model space partially occupied d 5/2 0g 9/2 1p 1/2, 0f 5/2 1p 3/2 0f 7/2 core: occupied 68 This model space is wide enough to discuss how magic numbers 28, 50 and semi-magic number 40 are retained/broken Interaction: A3DA interaction is used with minor corrections

13 66 (Z=28) の変形共存 66 (Z=28,N=38) で変形共存が予言される実験では 0 + 2,3,4 状態から球形の状態へと遷移する半減期を測定し MCSM 計算と比較 oblate Q(fm 2 ) Q(fm 2 ) prolate spherical 実験グループとの共著論文 : S. Leoni et al., PRL 118, (2017)

14 Comparison to experiments for 66 S. Leoni et al., PRL 118, (2017)

15 ベータ崩壊 Fermi 遷移 : ΔL=0, ΔJ=0, Δπ=+ Gamow-Teller (GT) 遷移 : ΔL=0, ΔJ=1, Δπ=+ 第一禁止遷移 : ΔL=1, ΔJ=0,1,2, Δπ=- B GT = & ( ) (* ( +, -. / 01 GT 遷移のみを考えるとき ft = f: phase space factor t: 半減期 (9 : /9 < ) >?(@A)

16 A=66 の GT 遷移 B(GT) の大きい (log ft の小さい ) 遷移について実験値を再現している

17 で示した球形の状態間の遷移が起こりやすい

Shape Evolution of isotopes 0 + 1 states of 68-78 68 70 72 Energy of prolate state comes

18 Shape Evolution of isotopes states of Energy of prolate state comes down Barrier becomes low γ-soft deformation (strong fluctuation in the γ direction)

19 Extension of model space for region unoccupied pfg9d5 model space partially occupied unoccupied pfsdg model space 0g 7/2, 1d 3/2 2s 1/2 partially occupied 50 1d 5/2 50 1d 5/ g 9/2 1p 1/2, 0f 5/2 1p 3/2 0f 7/ g 9/2 1p 1/2, 0f 5/2 1p 3/2 0f 7/ core: occupied core: occupied

20 78 in pfg9d5 and pfsdg model spaces pfg9d5 pfsdg pfg9d pfsdg 0 + 1

21 78 の実験と計算の比較 R. Taniuchi et al., Nature 569, 53(2019) 変形状態様々な理論計算

22 Summary モンテカルロ殻模型 (MCSM) により中重核から重い核までを計算 MCSM の特性を用いた T-plot による原子核形状の解析幅広い質量領域の核種を MCSM で計算し原子核形状を中心とする性質を議論原子核形状とベータ崩壊の遷移強度の関係 Hg など重い原子核の計算 r プロセス核の計算も可能に

スライド 1

スライド 1 原子核物理学 8. 中性子過剰核平成 20 年度原子核物理学 1 核図表存在すると予想される核種約 6000 114 存在が確認された核種約 3000 184 安定な核種 207 ほとんどの不安定核は中性子過剰核 82 126 50 82 r-process (rapid neutron capture) 8 20 8 28 20 28 50 70 超新星爆発における爆発的元素合成中性子捕獲反応とベータ崩壊の競争