NuShell@MSU
いくつかの論文中では、NuShell@MSU で shell model 計算 を行っているため、自分でも結果を再現できる。以下では、論文の値を再現する入力ファイルのサンプルを置く。NuShell@MSU のインストール方法は Programs Memo を参照。
実行方法
たとえば、下記の ne30gs_in.ans を実行するには、以下のようにすればよい。
NuShell 用のコマンドプロンプトを起動。
c:\aaa\NuShell>
- rsh ディレクトリに移動したら、計算用のディレクトリを作成。
c:\aaa\NuShell> cd rsh c:\aaa\NuShell\rsh> mkdir ne30gs c:\aaa\NuShell\rsh> cd ne30gs c:\aaa\NuShell\rsh\ne30gs>
このディレクトリに ne30gs_in.ans を保存。
NuShell の入力として ne30gs_in.ans の内容をリダイレクト。
c:\aaa\NuShell\rsh\ne30gs> shell < ne30gs_in.ans
- ne30gs.bat というバッチファイルが出来る。
- バッチファイルの実行。(拡張子 .bat は省略できる。)
c:\aaa\NuShell\rsh\ne30gs> ne30gs
- エネルギー準位や基底状態の成分を示すファイルが出力される。
エネルギー準位 : ne30gs.lpt
基底状態の成分 : l0a0eh.lpe
PRC41(1990)1147
Phys. Rev. C 41, 1147 (1990) Table IV の 30Ne の Binding energy : 90754 keV (Full(2hbarw)) を再現した。
- Model space : sdpf
- Interaction : WBMB (sdpfmw)
- Trancation : 2p2h without 0p0h
- 得られた Binding energy : -90.755 MeV
入力ファイル : ne30gs_in.ans
ne30gs lpe, 1 ! option (lpe or lan), neig (zero=10) ----------- sdpf ! model space (*.sp) name (a8) 14 ! number of valence particles y ! any restrictions (y/n) h 2,2 sdpfmw ! interaction (*.int) name (a8) 0.0, 0.0, ! min, max J (f,f) 5.0, 5.0, ! min, max T (f,f) 0.0, 0.0, ! J,T to leave out (end with zero or blank) 0 ! parity (0 for +) (1 for -) (2 for both) 28 0.0000 ! A value for sd-shell neucleons (i) st ! option -----------
出力エネルギー準位 : ne30gs.lpt
Interaction file information from ne30gs.mit spe taken from first line in *.int file Interaction spe sdpfw -3.9478 1.6466 -3.1635 0.0000 0.0000 0.0000 0.0000 sdpfmw 0.0000 0.0000 0.0000 14.0080 18.5240 13.5350 13.5230 sdpfmc 0.0000 0.0000 0.0000 -8.3630 -8.3630 -8.3630 -8.3630 sdpfcm 0.0000 0.0000 0.0000 3.0000 3.0000 3.0000 3.0000 Interaction file information from l0a0eh.nus Interaction spe-norm tbme-norm sdpfw 1.00000 1.11294 N name Njtp T E(MeV) J Ex(MeV) (* for yrast) 1 l0a0eh 1 5 -90.755 0 + 0.000 * .
基底状態の成分 : l0a0eh.lpe
NuShell@MSU V1.0 R1.000 B.A. Brown and W.D.M. Rae 2007 Based on NuShell code by W.D.M. Rae, Garsington, Oxford, UK, 2006/7 onejtih 2*J, 2*T 0 10 Lowest Energy Levels Level 1 E= -90.755 MeV Parity + 4.768 % Partition 6 4 2 2 0 0 0 1.476 % Partition 7 3 2 1 0 1 0 3.116 % Partition 7 3 2 2 0 0 0 1.067 % Partition 6 3 3 1 0 1 0 1.287 % Partition 6 3 3 2 0 0 0 6.113 % Partition 8 2 2 1 0 1 0 4.492 % Partition 7 2 3 1 0 1 0 45.163 % Partition 8 2 2 2 0 0 0 3.976 % Partition 7 2 3 2 0 0 0 2.234 % Partition 8 2 2 0 2 0 0 7.273 % Partition 8 2 2 0 0 2 0 2.801 % Partition 6 2 4 2 0 0 0 1.459 % Partition 8 2 2 0 0 0 2 Average n 7.499 2.318 2.183 1.450 0.093 0.382 0.074 Norm, total %, total n 0.9999998 99.99999 14.00000 Calculations successfully completed and data written to disk.
PRC81(2010)041302(R)
Phys. Rev. C 81, 041302(R) (2010), FIG. 1, 30Ne のエネルギー準位 の "USD" の計算を再現した。
- Model space : sd
- Interaction : USDA (usda)
- Restriction : なし
入力ファイル : ne30usda_in.ans
出力エネルギー準位 : ne30usda.lpt
First 2+ energy : 1800 MeV (論文), 1765 keV (NuShell 出力)
First 4+ energy : 2500 MeV (論文), 2516 keV (NuShell 出力)
PRC84(2011)037302
Phys. Rev. C 84, 037302 (2011), FIG. 3, 25F のエネルギー準位を再現した。
- Model space : spsdpf
- Interaction : WBP (wbp)
- Trancation : sd-shell (for positive parity), p-, sd-, pf-shell (for negative parity)
- Restriction : 0hbw (for positive parity), 1hbw (for negative parity)
入力ファイル : f25_in.ans
出力エネルギー準位 : f25.lpt
Spin-parity |
論文 |
NuShell 出力 |
5/2+ |
0 keV |
0 keV |
1/2+ |
911 keV |
911 keV |
3/2+ |
3074 keV |
3074 keV |
9/2+ |
3230 keV |
3230 keV |
5/2+ |
3757 keV |
3757 keV |
3/2+ |
- |
3988 keV |
7/2+ |
- |
4052 keV |
1/2+ |
- |
4070 keV |
1/2- |
4296 keV |
4296 keV |
5/2+ |
4800 keV |
4800 keV |
PRL102(2009)152501
Phys. Rev. Lett. 102, 152501 (2009), Table I, 23O の 5/2+ のエネルギーと C2S を再現した。
- Model space : sd
- Interaction : USDB (usdb)
- Restriction : none
23O の励起準位の計算用入力ファイル : o23_in.ans
23O の励起準位の結果 : o23.lpt
- 5/2+ のエネルギー : 2.593 MeV を再現。
23O と 24O の間の C2S 計算用入力ファイル : o23o24_in.ans
23O(1/2+) の C2S の結果: b1707b080.lsf
23O(1/2+) の C2S : 1.810 を再現。
23O(5/2+) の C2S の結果: b5707b080.lsf
23O(5/2+) の C2S : 5.665 を再現。
PRC86(2012)054604
- Model space : sd
- Interaction : WBP
- Restriction : 0hbw (for positive parity), 1hbw (for negative parity)
22C の励起準位の計算用入力ファイル : 22C_in.ans
22C の励起準位の結果 : 22C.lev
21C の励起準位の計算用入力ファイル : 21C_in.ans
21C の励起準位の結果 : 21C.lev
5/2+1 のエネルギー : 1.109 MeV を再現。
3/2+1 のエネルギー : 2.191 MeV を再現。
21C と 22C の間の C2S 計算用入力ファイル : 21C22C_in.ans
21C(1/2+1) の C2S の結果: t190lp0a0.lsf
21C(1/2+1) の C2S : 1.403 を再現。
21C(4/2+1) の C2S の結果: t590lp0a0.lsf
21C(4/2+1) の C2S : 4.212 を再現。
21C(3/2+1) の C2S の結果: t390lp0a0.lsf
21C(3/2+1) の C2S : 0.342 を再現。
F.C. Barker, Aust. J. Phys., 33(1980)177
Title: The 7Be(p,γ)8B Cross section at Low Energies
- S. Cohen と D. Kurath の相互作用
- Model space name: ppn
- Interaction name: ckipn (S. Cohen and D. Kurath, Nucl. Phys. A73, 1 (1965) and Nucl. Phys. A101, 1 (1967))
- Restriction : None
7Be と 8B の間の Spectroscopic factor 計算用入力ファイル : 7be8b_ckipn_in.ans
8B の 2+, 1+, 3+ 状態へ遷移する Spectroscopic factor の結果: a31332422.lsf, a31332222.lsf, a31332622.lsf
- Table 1 の Calculated values の CK の欄の値を再現
J = 2: (論文) 0.282 + 0.751 = 1.033, (NuShell@MSU) 0.056141 + 0.977015 = 1.033156
J = 1: (論文) 0.140 + 0.306 = 0.446, (NuShell@MSU) 0.123977 + 0.321496 = 0.445473
J = 3: (論文) 0.338, (NuShell@MSU) 0.337577
- Table 1 の Calculated values の CK の欄の値を再現
- Kumar の相互作用
- Model space name: ppn
- Interaction name: pkupn (N. Kumar, NPA225, 221 (1974))
- label.dat の pkupn の 説明には ppn pkupn a d ! (8-16)TBME Nucl. Phys. A101, 1 (1967) とあるが、間違いのようだ。label.old の説明が正しいようで、ここには、PPN PKUPN A D ! N. Kumar, NPA225, 221 (1974) と書いてある。
- Restriction : None
7Be と 8B の間の Spectroscopic factor 計算用入力ファイル : 7be8b_pkupn_in.ans
8B の 2+, 1+, 3+ 状態へ遷移する Spectroscopic factor の結果: x3133y422.lsf, x3133y222.lsf, x3133y622.lsf
- Table 1 の Calculated values の CK の欄の値を再現
J = 2: (論文) 0.250 + 0.751 = 1.001, (NuShell@MSU) 0.067299 + 0.932811 = 1.000110
J = 1: (論文) 0.206 + 0.221 = 0.427, (NuShell@MSU) 0.059451 + 0.367916 = 0.427367
J = 3: (論文) 0.308, (NuShell@MSU) 0.308232
- Table 1 の Calculated values の CK の欄の値を再現
J. A. Tostevin et al, Phys. Rev. C 70, 064602 (2004)
- Title: Correlated two-nucleon stripping reactions
- Model space : sd
- Interaction : USD (Interaction name: w)
- Restriction : None
28Mg と 26Ne の間の Two-nucleon spectroscopic amplitudes (TNA) 計算用入力ファイル : 26ne28mg_in.ans
- TABLE II の値を再現(ただし、符号の決め方がいまいちよくわからない。OXBASH を使うと、値がかなり違い、符号も違う。)
NuShell@MSU (b060aw040.tna, b460aw040.tna, b860aw040.tna)
OXBASH (b060aw040_oxbash.tna, b460aw040_oxbash.tna, b860aw040_oxbash.tna)
- TABLE II の値を再現(ただし、符号の決め方がいまいちよくわからない。OXBASH を使うと、値がかなり違い、符号も違う。)
Jπf |
E* (MeV) |
[0d3/2]2 |
[0d3/20d5/2] |
[0d5/2]2 |
[1s1/20d3/2] |
[1s1/20d5/2] |
[121/2]2 |
Paper |
|||||||
0+1 |
0.0 |
-0.30146 |
- |
-1.04685 |
- |
- |
-0.30496 |
2+1 |
2.02 |
-0.05030 |
0.37358 |
-0.63652 |
-0.06084 |
-0.13916 |
- |
4+1 |
3.50 |
- |
0.33134 |
1.59639 |
- |
- |
- |
2+2 |
3.70 |
0.04721 |
-0.07248 |
0.85297 |
0.16158 |
0.17590 |
- |
NuShell@MSU |
|||||||
0+1 |
0.0 |
0.30146 |
- |
1.04685 |
- |
- |
0.30496 |
2+1 |
2.02 |
-0.05030 |
-0.37357 |
-0.63651 |
0.06083 |
-0.13915 |
- |
4+1 |
3.50 |
- |
0.33131 |
-1.59646 |
- |
- |
- |
2+2 |
3.70 |
-0.04721 |
-0.07248 |
-0.85298 |
0.16158 |
-0.17591 |
- |
OXBASH |
|||||||
0+1 |
0.0 |
0.30169 |
- |
1.04754 |
- |
- |
0.30408 |
2+1 |
2.02 |
-0.05026 |
0.37367 |
-0.63730 |
-0.06053 |
-0.13876 |
- |
4+1 |
3.50 |
- |
-0.33199 |
-1.59754 |
- |
- |
- |
2+2 |
3.70 |
0.04693 |
-0.07256 |
0.85171 |
0.16073 |
0.17632 |
- |