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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » HSE计算的一些问题，希望大家帮我解决谢谢
141/2返回列表12下一页
查看: 3100  |  回复: 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的回帖

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的回帖

wangyan6930

木虫 (著名写手)
引用回帖:
4楼: Originally posted by wangyan6930 at 2015-01-07 14:39:36
楼主你好，我就是利用脚本程序处理出来的能带来看带隙的，不是从DOS上看的，看出来还是有点小，能不能从参数上分析一下我的有什问题。...

在朋友的指点下 我发现了自己一个参数设置的错 不知道对结果有没有影响EDIFFG=2e-4设置错了应该设置成EDIFFG=-2e-4，但是不知道对计算结果有没有影响，我会继续做测试，也请各位高手继续指出我的错误。
赞一下 回复此楼
5楼2015-01-07 15:41:43
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

legodgz

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
liliangfang: 金币+1, 谢谢交流 2015-01-15 07:52:14
wangyan6930: 金币+10, ★有帮助 2015-01-21 01:14:36
EDIFF: allowed error in total energy
我估计不会是因为这个参数的问题
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

wangyan6930
6楼2015-01-14 10:19:11
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wangyan6930

木虫 (著名写手)
送红花一朵
引用回帖:
6楼: Originally posted by legodgz at 2015-01-14 10:19:11
EDIFF: allowed error in total energy
我估计不会是因为这个参数的问题

确实不是因为这个参数 我计算完了后 带隙跟原来计算的一样 更实验差距0.7 eV. 你能不能检查下我的计算流程了计算参数有没有什么大的问题？ 当然我跟朋友讨论，他说AEXX的值影响带隙，可以调节AEXX，你有什么看法和建议给我提醒一下，谢谢
赞一下 回复此楼
7楼2015-01-14 14:05:33
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

未岸亡灵

铁杆木虫 (著名写手)

暂时离开小木虫

【答案】应助回帖

王同学呀，我们是做理论计算的，我想重点应该是先了解理论，再去做计算。
即使改完参数，你最终得到的值和实验“偶然”吻合上了，有什么用呢？
你也别多心，别人告诉你EDIFF和AEXX的问题，你或许都没去查过这两个参数到底是什么意思，你就直接改了，这样做是很浪费时间了，尤其是hse这么吃计算量的情况下。
你确实很执着，算了一个月，我的确已经很佩服你。但是，我想，我们做理论计算的，最好还是去书，至少应该去看下manual怎么讲，因为，理论和实验室相辅相成的，实验在一定程度上也不一定正确，因为才需要理论的验证。
你再想想，实验上得到的3ev的带隙，实验中用的真实模型，和你理论中人为构建的理想模型肯定是有差异的，不一定非要一定和实验完全吻合上，一个近似理论完全和实验值吻合上了，你认为这是什么一个情况呢？
毕竟量化发展不过百年，真正要做到定量还有很多艰难险阻要走，要跨越，这才是我们真正要去做的。而现在，我们可以定性地描述大致趋势，做到有据可依。我也由衷地希望有一天，量化可以做到定量，能够切实客观地反应出我们真是的世界，我们为之而努力。
说了这么多，也没有切实地帮到你，深感抱歉！
赞一下 回复此楼
8楼2015-01-20 22:36:40
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wangyan6930

木虫 (著名写手)
引用回帖:
8楼: Originally posted by 未岸亡灵 at 2015-01-20 22:36:40
王同学呀，我们是做理论计算的，我想重点应该是先了解理论，再去做计算。
即使改完参数，你最终得到的值和实验“偶然”吻合上了，有什么用呢？
你也别多心，别人告诉你EDIFF和AEXX的问题，你或许都没去查过这两个 ...

首先谢谢！你应该是我的好友吧？你怎么知道我测试了多久？ 其实很多问题我看手册论坛已经搞明白了   但是我不确定    我手册和相关网页基本看完了   我们组只有我在使用 vasp   我只想确认一下是否合理   因为我要保证最后计算的时候没有一点问题     再次我测试的例子很小   一天可以走几趟   我只用了几个核   我其他工作还在继续   不影响正常工作   当然在测试过程中我也探索出如何节约时间   设置省时的参数    虽然现在麻烦   但是今后我计算大体系可以 很快           其实懂和不懂就差一点   有时候需要大家帮忙捅破这一张纸   在这里谢谢论坛里的朋友关注

[ 发自小木虫客户端 ]
赞一下 回复此楼
9楼2015-01-21 00:34:33
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wangyan6930

木虫 (著名写手)
引用回帖:
9楼: Originally posted by wangyan6930 at 2015-01-21 00:34:33
首先谢谢！你应该是我的好友吧？你怎么知道我测试了多久？ 其实很多问题我看手册论坛已经搞明白了   但是我不确定    我手册和相关网页基本看完了   我们组只有我在使用 vasp   我只想确认一下是否合理   因为我要 ...

在这里也想告诉大家我上面的参数也可以作为一个例子来测试

[ 发自小木虫客户端 ]
赞一下 回复此楼
10楼2015-01-21 00:41:49
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangyan6930 的主题更新
141/2返回列表12下一页
信息提示
关闭
请填处理意见
关闭