|
1、赝势截止半径交叠大小的判断
⑴ 从POTCAR中 RCORE参数读取赝势截止半径,其中单位为Bohr
如铀的PBE赝势U: RCORE = 2.800 Bohr = 1.482 Å, 记RU = 1.482 Å.
⑵ 计算高压下体系的最近邻原子之间的键长
如金属铀在500 GPa时的最近邻U-U键键长是2.389 Å,记RU-U = 2.389 Å.
⑶ 赝势截止半径交叠率的计算
金属铀在500 GPa时,最近邻U原子的半径交叠率
, 计算结果为38.8%.
可以看出在500 GPa,最近邻的U原子之间截止半径交叠明显,高压下的计算结果可能不可靠.
2、对高压计算结果的可靠性进行验证
选取不交叠的赝势或者较小交叠的赝势对计算结果进行重新评估,计算相关的性质比较差异大小,如相变压力,状态方程,弹性常数等.
⑴ 首先选取VASP自带的截止半径更小的赝势验证计算结果,GW赝势半径相对较小可以作为比较.
⑵ 选择其他计算软件的赝势,如CASTEP的超软赝势或者QE的超软赝势验证,对于金属铀,CASTEP和QE的超软赝势在500 GPa都没有半径交叠.
对于金属铀,CASTEP和QE计算得到的相变压力和状态方程,与VASP计算的结果只有较小的误差,小于2%,从而证明了VASP高压计算结果的可靠性.
3、高压下计算结果的更精确验证可以进一步考虑全电子的计算,如WIEN2K软件
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-20 01:31
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社