应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » HSE计算的一些问题，希望大家帮我解决谢谢
52/1返回列表
查看: 3745  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wangyan6930

木虫 (著名写手)

[求助] HSE计算的一些问题，希望大家帮我解决谢谢 已有4人参与

经过最近的多次测试，兄弟想请教如下几个问题，我也是做了多次测试很出的结果，已经成功的计算了Si的能带，并和实验值符合的很好，但是测试新的例子的时候出下了以下几个问题，请大家帮忙：
1、关于布里渊区K点的选择，我现在做的测试是NaNbO3的能带计算，空间群是221PM-3M对称性，高对称点是X（0.5，0，0） R（0.5，0.5，0.5） M （0.5，0.5，0），G（0，0，0）,我选择的路径是X-R-M-G-R(CASTEP默认路径）和R-X-G-M-R-G分别进行了计算，想请教一下KPOINTS的路径怎么样选取更加合理？
2、我HSE计算能带的文件里面我同时处理了态密度，发现从DOS上读取的带隙和能带图上读取的带隙有一些差距，这个是由于DOS计算的时候加展宽的原因吗？还是同一个计算结果里面同时选择的DOS和能带本来就不合适？也就是HSE的DOS和能带应该分开来算？3
HSE计算出来的带隙约为2.5左边，跟实验值的3.19有些差距，我选择了我的两次不同计算过程的文件给您发过去，请您帮我查看一下错误可以吗？积分路径是RXGMRG 计算出来的能带结果是2.5左右,能不能麻烦帮我看下问题出在什么地方，附件是我选的两次不同计算路径的参数，麻烦帮我看一下，谢谢！

我把我计算的参数贴出来请大家给我指出问题，积分路径R-X-G-M-R-G，主要问题是计算出来的能带和实验值有差距，反复测试了数十次和查资料实在找不到问题，请大家帮忙，谢谢！
PBE计算

INCAR

ENCUT = 520
EDIFF = 1e-5
EDIFFG = 2e-4
ISMEAR=0
SIGMA=0.1
IBRION=-1
NSW=0
ISIF=2
POTIM=0.1
LWAVE=T
LCHARG=T
#NPAR=1

KPOINTS
Automatic mesh
0
Monkhorst Pack
  5   5   5
  0.  0.  0.

POSCAR
Na Nb O3  
    1.00000000000000
     4.0096004013579449    0.0000000011715256    0.0000000011715255
     0.0000000723430511    4.0096004013579449    0.0000000011715256
     0.0000000723430511    0.0000000723430511    4.0096004013579449
   Na   Nb   O
     1     1     3
Direct
  0.0000000000000000  0.0000000000000000  0.0000000000000000
  0.5000000000000000  0.5000000000000000  0.5000000000000000
  0.0000000000000000  0.5000000000000000  0.5000000000000000
  0.5000000000000000  0.5000000000000000 -0.0000000000000000
  0.5000000000000000  0.0000000000000000  0.5000000000000000

HSE自洽计算
INCAR

ENCUT = 520
EDIFF = 1e-5
EDIFFG = 2e-4
ISMEAR=0
SIGMA=0.1
IBRION=-1
NSW=0
ISIF=2
POTIM=0.1
LWAVE=T
LCHARG=T
#NPAR=1
LHFCALC=.TRUE.
HFSCREEN=0.2
ALGO = D
TIME = 0.4

KPOINTS
Automatic mesh
0
Monkhorst Pack
  5   5   5
  0.  0.  0.

HSE能带计算
INCAR
ENCUT = 520
EDIFF = 1e-5
EDIFFG = 2e-4
ISMEAR=0
SIGMA=0.1
IBRION=-1
NSW=0
ISIF=2
POTIM=0.1
LWAVE=T
LCHARG=T
#NPAR=1
LHFCALC=.TRUE.
HFSCREEN=0.2

KPOINTS

Automatically generated mesh
     60
Reciprocal lattice
    0.00000000000000    0.00000000000000    0.00000000000000             1
    0.20000000000000    0.00000000000000    0.00000000000000             6
    0.40000000000000    0.00000000000000    0.00000000000000             6
    0.20000000000000    0.20000000000000    0.00000000000000            12
    0.40000000000000    0.20000000000000    0.00000000000000            24
    0.40000000000000    0.40000000000000    0.00000000000000            12
    0.20000000000000    0.20000000000000    0.20000000000000             8
    0.40000000000000    0.20000000000000    0.20000000000000            24
    0.40000000000000    0.40000000000000    0.20000000000000            24
    0.40000000000000    0.40000000000000    0.40000000000000             8
   0.50000000  0.50000000  0.50000000       0
   0.50000000  0.44444444  0.44444444       0
   0.50000000  0.38888889  0.38888889       0
   0.50000000  0.33333333  0.33333333       0
   0.50000000  0.27777778  0.27777778       0
   0.50000000  0.22222222  0.22222222       0
   0.50000000  0.16666667  0.16666667       0
   0.50000000  0.11111111  0.11111111       0
   0.50000000  0.05555556  0.05555556       0
   0.50000000  0.00000000  0.00000000       0
   0.50000000  0.00000000  0.00000000       0
   0.44444444  0.00000000  0.00000000       0
   0.38888889  0.00000000  0.00000000       0
   0.33333333  0.00000000  0.00000000       0
   0.27777778  0.00000000  0.00000000       0
   0.22222222  0.00000000  0.00000000       0
   0.16666667  0.00000000  0.00000000       0
   0.11111111  0.00000000  0.00000000       0
   0.05555556  0.00000000  0.00000000       0
   0.00000000  0.00000000  0.00000000       0
   0.00000000  0.00000000  0.00000000       0
   0.05555556  0.05555556  0.00000000       0
   0.11111111  0.11111111  0.00000000       0
   0.16666667  0.16666667  0.00000000       0
   0.22222222  0.22222222  0.00000000       0
   0.27777778  0.27777778  0.00000000       0
   0.33333333  0.33333333  0.00000000       0
   0.38888889  0.38888889  0.00000000       0
   0.44444444  0.44444444  0.00000000       0
   0.50000000  0.50000000  0.00000000       0
   0.50000000  0.50000000  0.00000000       0
   0.50000000  0.50000000  0.05555556       0
   0.50000000  0.50000000  0.11111111       0
   0.50000000  0.50000000  0.16666667       0
   0.50000000  0.50000000  0.22222222       0
   0.50000000  0.50000000  0.27777778       0
   0.50000000  0.50000000  0.33333333       0
   0.50000000  0.50000000  0.38888889       0
   0.50000000  0.50000000  0.44444444       0
   0.50000000  0.50000000  0.50000000       0
   0.50000000  0.50000000  0.50000000       0
   0.44444444  0.44444444  0.44444444       0
   0.38888889  0.38888889  0.38888889       0
   0.33333333  0.33333333  0.33333333       0
   0.27777778  0.27777778  0.27777778       0
   0.22222222  0.22222222  0.22222222       0
   0.16666667  0.16666667  0.16666667       0
   0.11111111  0.11111111  0.11111111       0
   0.05555556  0.05555556  0.05555556       0
   0.00000000  0.00000000  0.00000000       0
回复此楼

» 猜你喜欢
1楼 2015-01-06 23:23:42
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wangyan6930

木虫 (著名写手)
送红花一朵
引用回帖:
11楼: Originally posted by 未岸亡灵 at 2015-01-21 07:47:56
恭喜你...

谢谢

[ 发自小木虫客户端 ]
回复此楼
14楼2015-01-21 08:02:15
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

chuanghua304

禁虫 (职业作家)
★ ★ ★ ★ ★
感谢参与，应助指数 +1
wangyan6930: 金币+5, ★有帮助, 楼上有没有用VASP做HSE计算能带，我就是想讨论一下我碰到的问题，我的KPOINTS是4*4*4的，确实不是很；晶格常数是计算所得，跟实验值接近；ISMEAR=-5对能带计算有没有影响？HSE计算好像大家使用的ISMEAR=0 来做HSE计算的，我可以和楼主交流一下不？ 2015-01-07 09:54:16
本帖内容被屏蔽
2楼2015-01-07 08:29:31
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lyd924

金虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
ljw4010: 金币+2, 谢谢交流！ 2015-01-07 21:23:39
wangyan6930: 金币+5, ★有帮助 2015-01-21 01:13:53
能隙应该主要从能带来判断吧，态密度展宽与smear有关，会造成能隙与能带图中的有点差异
赞一下 回复此楼
3楼2015-01-07 11:47:14
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wangyan6930

木虫 (著名写手)
引用回帖:
3楼: Originally posted by lyd924 at 2015-01-07 11:47:14
能隙应该主要从能带来判断吧，态密度展宽与smear有关，会造成能隙与能带图中的有点差异

楼主你好，我就是利用脚本程序处理出来的能带来看带隙的，不是从DOS上看的，看出来还是有点小，能不能从参数上分析一下我的有什问题。
赞一下 回复此楼
4楼2015-01-07 14:39:36
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答
信息提示
关闭
请填处理意见
关闭