下載地址
相關(guān)軟件
原子核質(zhì)量表更新:new nuclear mass table(26)
軟件簡(jiǎn)介
基于宏觀(guān)-微觀(guān)方法計算原子核的結合能以及基態(tài)形變
作者:王 寧
聯(lián)系方式 Email: wangning@gxnu.edu.cn; ningmany@gmail.com
Web: http://www.hfljqc.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; 以及最優(yōu)化DSC方法來(lái)搜索能量最小值。
2)AME2003.dat : 原子核質(zhì)量過(guò)剩(mass excess)的實(shí)驗數據:用來(lái)和理論計算對比。
3)INPUT.WS : 用戶(hù)自己定義模型參數;主要是液滴能、單粒子勢,以及是否考慮鏡像核約束IMIRR的控制參數。
4)EMASS.DAT : 計算結果
5)PES.DAT : 如果用戶(hù)選擇輸出勢能面E(beta2,beta4), 結果輸出至該文件。
6)WS3.2.dat : 采用WS 方法對整個(gè)核素表的計算結果; 參看Phys.Rev.C81(2010)044322
7)WS3.3.dat : 采用WS*方法對整個(gè)核素表的計算結果; 參看Phys.Rev.C82(2010)044304
8)WS3.3.exe : Intel Fortran編譯器產(chǎn)生的可執行文件。
Fortran程序在Intel Fortran編譯器下測試過(guò)
運行程序后需要對一些輸入條件做選擇:
第一步:input KEY: 參數組選擇, 0 --> WS; 1 --> WS*; 2 --> 用戶(hù)自定義(設置INPUT.WS中的相應參數)
第二步:勢能面PES是否輸出選擇, 0 --> 不輸出(節省時(shí)間); 1 --> 輸出在PES.dat中
第三步:只針對用戶(hù)自定義參數情況(否則直接到下一步);選擇要考慮的原子核形變自由度;最多考慮四個(gè)自由度 beta2,beta3,beta4,beta6 (NDIM=4);
第四步:選擇要計算的原子核,質(zhì)子數Z和質(zhì)量數A; 這塊兒可以批量計算,設置原子核質(zhì)子數和質(zhì)量數的范圍: (Zmin, Zmax) (Amin, Amax);
如果只計算一個(gè)核,則設置Zmin=Zmax=Z, Amin=Amax=A即可。如果只計算已測得2149個(gè)原子核質(zhì)量,選擇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.hfljqc.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
---------------------------------------------------------------