Locked History Actions

NuShell@MSU

NuShell@MSU

  • いくつかの論文中では、NuShell@MSU で shell model 計算 を行っているため、自分でも結果を再現できる。以下では、論文の値を再現する入力ファイルのサンプルを置く。NuShell@MSU のインストール方法は Programs Memo を参照。

実行方法

  • たとえば、下記の ne30gs_in.ans を実行するには、以下のようにすればよい。

  1. NuShell 用のコマンドプロンプトを起動。

    • c:\aaa\NuShell>
  2. rsh ディレクトリに移動したら、計算用のディレクトリを作成。
    • c:\aaa\NuShell> cd rsh
      c:\aaa\NuShell\rsh> mkdir ne30gs
      c:\aaa\NuShell\rsh> cd ne30gs
      c:\aaa\NuShell\rsh\ne30gs>
  3. このディレクトリに ne30gs_in.ans を保存。

  4. NuShell の入力として ne30gs_in.ans の内容をリダイレクト。

    • c:\aaa\NuShell\rsh\ne30gs> shell < ne30gs_in.ans
  5. ne30gs.bat というバッチファイルが出来る。
  6. バッチファイルの実行。(拡張子 .bat は省略できる。)
    • c:\aaa\NuShell\rsh\ne30gs> ne30gs
  7. エネルギー準位や基底状態の成分を示すファイルが出力される。

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

  • 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

J. A. Tostevin et al, Phys. Rev. C 70, 064602 (2004)

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

-