PallasInstitute的个人博客分享 http://blog.sciencenet.cn/u/PallasInstitute

博文

VASP中PBE赝势在高压下交叠明显的处理方法

已有 2786 次阅读 2022-5-5 16:25 |系统分类:科研笔记

1、赝势截止半径交叠大小的判断

POTCAR RCORE参数读取赝势截止半径,其中单位为Bohr

 如铀的PBE赝势U RCORE = 2.800 Bohr = 1.482 Å R= 1.482 Å.

计算高压下体系的最近邻原子之间的键长

   如金属铀在500 GPa时的最近邻U-U键键长是2.389 Å,记RU-U = 2.389 Å.

赝势截止半径交叠率的计算

金属铀在500 GPa时,最近邻U原子的半径交叠率

image.png

, 计算结果为38.8%.

可以看出在500 GPa,最近邻的U原子之间截止半径交叠明显,高压下的计算结果可能不可靠.

 

2、对高压计算结果的可靠性进行验证

选取不交叠的赝势或者较小交叠的赝势对计算结果进行重新评估,计算相关的性质比较差异大小,如相变压力,状态方程,弹性常数等.

首先选取VASP自带的截止半径更小的赝势验证计算结果,GW赝势半径相对较小可以作为比较.

选择其他计算软件的赝势,如CASTEP的超软赝势或者QE的超软赝势验证,对于金属铀,CASTEPQE的超软赝势在500 GPa没有半径交叠.

对于金属铀,CASTEPQE计算得到的相变压力和状态方程,与VASP计算的结果只有较小的误差,小于2%,从而证明了VASP高压计算结果的可靠性.

 

3、高压下计算结果的更精确验证可以进一步考虑全电子的计算,如WIEN2K软件




https://wap.sciencenet.cn/blog-3469498-1337167.html

上一篇:[转载]ubuntu桌面卡死解决
下一篇:[转载]一些常见的材料学数据库
收藏 IP: 101.206.170.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-5-16 01:21

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部