下载地址
相关软件
原子核质量表更新:new nuclear mass table(60)
软件简介
基于宏观-微观方法计算原子核的结合能以及基态形变
作者:王 宁
联系方式 Email: wangning@gxnu.edu.cn; ningmany@gmail.com
Web: http://www.imqmd.com/wangning
参考文献: [1] Ning Wang, Min Liu and Xizhen Wu, Phys. Rev. C 81 (2010) 044322
[2] Ning Wang, et al., Phys. Rev. C 82 (2010) 044304
版本号3.3
程序文件以及参数:
1)WS3.3.for : 该方法的Fortran程序; 里面镶嵌了Woods-Saxon势下单粒子能级的求解程序WSBETA; 以及最优化DSC方法来搜索能量最小值。
2)AME2003.dat : 原子核质量过剩(mass excess)的实验数据:用来和理论计算对比。
3)INPUT.WS : 用户自己定义模型参数;主要是液滴能、单粒子势,以及是否考虑镜像核约束IMIRR的控制参数。
4)EMASS.DAT : 计算结果
5)PES.DAT : 如果用户选择输出势能面E(beta2,beta4), 结果输出至该文件。
6)WS3.2.dat : 采用WS 方法对整个核素表的计算结果; 参看Phys.Rev.C81(2010)044322
7)WS3.3.dat : 采用WS*方法对整个核素表的计算结果; 参看Phys.Rev.C82(2010)044304
8)WS3.3.exe : Intel Fortran编译器产生的可执行文件。
Fortran程序在Intel Fortran编译器下测试过
运行程序后需要对一些输入条件做选择:
第一步:input KEY: 参数组选择, 0 --> WS; 1 --> WS*; 2 --> 用户自定义(设置INPUT.WS中的相应参数)
第二步:势能面PES是否输出选择, 0 --> 不输出(节省时间); 1 --> 输出在PES.dat中
第三步:只针对用户自定义参数情况(否则直接到下一步);选择要考虑的原子核形变自由度;最多考虑四个自由度 beta2,beta3,beta4,beta6 (NDIM=4);
第四步:选择要计算的原子核,质子数Z和质量数A; 这块儿可以批量计算,设置原子核质子数和质量数的范围: (Zmin, Zmax) (Amin, Amax);
如果只计算一个核,则设置Zmin=Zmax=Z, Amin=Amax=A即可。如果只计算已测得2149个原子核质量,选择Zmin<=0, Zmax=110
==========================================================
A Fortran code for nuclear mass calculation based on the macroscopic-microscopic method
Reference:
[1] Ning Wang, Min Liu and Xizhen Wu, Phys. Rev. C 81 (2010) 044322
[2] Ning Wang, et al., Phys. Rev. C 82 (2010) 044304
Author: Ning Wang
Email: wangning@gxnu.edu.cn; ningmany@gmail.com
Web: http://www.imqmd.com/wangning
version 3.3
---------------------------------------------------------------
attached files
---------------------------------------------------------------
1) WS3.3.for : Fortran code to calculate nuclear binding energy, in which the WSBETA subroutine to solve the single particle levels under axially deformed Woods-Saxon potential is used; and the DSC method to find the minimum value on the potential energy surface.
2) AME2003.dat : measured nuclear mass excesses
3) INPUT.WS : input model parameters (used only for USER case)
4) EMASS.DAT : output binding energy of nucleus and deformations at the ground state
5) PES.DAT : the potential energy surface E(beta2,beta4)
6) WS3.2.dat : results with WS parameter set; (see Phys.Rev.C81(2010)044322)
7) WS3.3.dat : results with WS* parameter set; (see Phys.Rev.C82(2010)044304)
8) WS3.3.exe : execute file from Intel Fortran compiler v10
---------------------------------------------------------------
How to use
---------------------------------------------------------------
run the execute file, and then
First :input KEY: to choose the parameter set, 0 --> WS; 1 --> WS*; 2 --> USER (to set the values in INPUT.WS)
Second :whether record the potential energy surface or not, 0 --> No (to save time); 1 --> YES (in PES.dat)
Third :only for USER case (otherwise to next step);to set the degree of freedom of nuclear deformation;maximum (NDIM=4) beta2,beta3,beta4,beta6;
Finally:set the charge number Z and mass number A of nuclei under consideration; for a selected nucleus, set (Zmin=Zmax=Z, Amin=Amax=A); to calcualte the known 2149 nuclei,set Zmin<=0, Zmax=110
---------------------------------------------------------------