原子核壳模型程序 BigStick 的用法

BigStick 公开版本可在 https://github.com/cwjsdsu/BigstickPublick.git 下载。其中 docs/ 与 util/docs/ 下都有详细的文档说明，这里只是记下简单运用的步骤。
下面的脚本是根据 util/example/ 里面的脚本修改的，可以自动算能谱和跃迁：

#!/bin/bash

nucleus=Ti46
Np=2
Nn=4

echo ' '
echo ' Sample script to compute B(E2) for Ca48 using GX1A interaction '
echo ' '
echo ' running tropic1b to generate E2 operator '
echo ' '
./tropic1b.x << input
fp0    ! BIGSTICK-format .sps file needed for code
2      ! menu option 2: E2
2      ! proton-neutron
y      ! include cross-shell (irrelevant here)
1.86441      ! value of osc length = (A)^1/6 fm
1.5  0.5   ! effective charge for proton, neutron
fpe2pn  ! output file name
n      ! in single-column xpn format? n = double column pns format
input

echo ' '
echo ' * * * * * * * * * * * * * * * * '
echo ' * * * * * * * * * * * * * * * * '
echo ' creating density matrices '
echo ' * * * * * * * * * * * * * * * *'
echo ' * * * * * * * * * * * * * * * * '


./bigstick.x << INPUT
dp     ! menu choice -- proton neutron densities 
$nucleus                                                   
fp    !  name of .sp file  (in this case, NuShell format) 
           $Np           $Nn  ! # of valence protons, neutrons 
           0        ! 2 x Jz of systems 
upn     ! signal to read NuShell-format unnormalized interaction
gx1apn
end
ld    ! Lanczos menu option 
          5         300      ! # states to keep, max # iterations
INPUT

echo ' * * * * * * * '
echo ' * * * * * * * * * * * * * * * * '
echo ' * * * * * * * * * * * * * * * * '
echo 'finished with run '
echo ' * * * * * * * * * * * * * * * * '
echo ' * * * * * * * * * * * * * * * * '
echo ' * * * * * * * * '

./genstrength.x << input
fpe2pn        ! name of .opme 
1             ! no additional scaling
$nucleus      ! name of reference parent state file
3             ! # of parent states
$nucleus      ! name of reference daughter state file
3             ! # of daughter states
$nucleus      ! name of density file
n             ! do not read in additional density files
$nucleus.e2        ! name of output
input

echo ' '
echo " B(E2)s written to "$nucleus".e2.str"
echo ' '

• 其中用到的可执行文件 bigstick.x, tropic1b.x, genstrength.x 分别在 src/, util/trans/ 中编译生成，然后拷贝到一个文件夹即可。
• 其中用到输入文件 fp.sp, gx1apn.int (upn格式)
• 我还不是很清楚关键字　input, INPUT 可以不可以都写成 input，下次或可一试。
• 在这个脚本的基础上，可以很容易地写循环，生成若干不同核的数据。