原子核壳模型程序 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,下次或可一试。
- 在这个脚本的基础上,可以很容易地写循环,生成若干不同核的数据。
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】凌霞软件回馈社区,博客园 & 1Panel & Halo 联合会员上线
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】博客园社区专享云产品让利特惠,阿里云新客6.5折上折
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 用 C# 插值字符串处理器写一个 sscanf
· Java 中堆内存和栈内存上的数据分布和特点
· 开发中对象命名的一点思考
· .NET Core内存结构体系(Windows环境)底层原理浅谈
· C# 深度学习:对抗生成网络(GAN)训练头像生成模型
· 为什么说在企业级应用开发中,后端往往是效率杀手?
· 本地部署DeepSeek后,没有好看的交互界面怎么行!
· 趁着过年的时候手搓了一个低代码框架
· 用 C# 插值字符串处理器写一个 sscanf
· 推荐一个DeepSeek 大模型的免费 API 项目!兼容OpenAI接口!