24小时热门版块排行榜    

查看: 4519  |  回复: 14
本帖产生 1 个 1ST强帖 ,点击这里进行查看

缺心眼桌

木虫 (正式写手)

[求助] 静态自洽计算能量时EDIFFG的问题

大家好,我计算一个体系能量时,静态自洽计算的设置是
ENCUT=400
ISTART=0
ICHARG=2
ISMEAR=0
SIGMA=0.08
NSW=0
EDIFF=1E-5
EDIFFG=-0.05请问有什么问题吗?我的EDIFFG是负数,力的收敛标准,我原以为力的收敛标准只在优化结构的时候用,按理说静态自洽计算时不用设置这一项吧,那这里会对我的结果有什么影响吗?谢谢大家!
回复此楼

» 收录本帖的淘帖专辑推荐

VASP

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buct2010

金虫 (正式写手)

【答案】应助回帖

★
youzhizhe(金币+1, 1ST强帖+1): 谢谢交流。 2011-05-04 14:20:03
引用回帖:
Originally posted by 缺心眼桌 at 2011-05-04 09:48:46:
大家好,我计算一个体系能量时,静态自洽计算的设置是
ENCUT=400
ISTART=0
ICHARG=2
ISMEAR=0
SIGMA=0.08
NSW=0
EDIFF=1E-5
EDIFFG=-0.05请问有什么问题吗?我的EDIFFG是负数,力的收敛标准,我原 ...

你计算体系能量直接在OSIZCAR中找到就可以了。静态计算是为了得到CHGCAR,它只进行一次迭代计算。我觉得你的精度应该设置高一点,EDIFFG=-0.05偏大,-0.0001,-0.001都可以(因为你的SIGMA=0.08和EDIFF=1E-5的精度都是高的)。收敛标准是能量还是力我觉得影响不大,因为得到CHGCAR可以得到DOS和能带结构,它们都是定性分析的(定量的禁带宽度VASP的误差是很大的)。
2楼2011-05-04 11:10:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

缺心眼桌

木虫 (正式写手)

引用回帖:
Originally posted by buct2010 at 2011-05-04 11:10:49:
你计算体系能量直接在OSIZCAR中找到就可以了。静态计算是为了得到CHGCAR,它只进行一次迭代计算。我觉得你的精度应该设置高一点,EDIFFG=-0.05偏大,-0.0001,-0.001都可以(因为你的SIGMA=0.08和EDIFF=1E-5的 ...

我的体系很大,所以设置了EDIFFG=-0.05(参考了体系差不多的文献),我也只用了一个k点的,这是我当时优化结构时设置的,静态自洽计算时没改。手册里说如果需要准确的能量也必须再做一次静态自洽计算的,再就是我接下来也需要电荷密度和能带。我想确定的是,优化结构和静态自洽计算时,如果标准既有EDIFF,也有EDIFFG,那程序会按照哪一个来结束?谢谢!
3楼2011-05-04 15:51:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buct2010

金虫 (正式写手)

【答案】应助回帖

★
bingmou(金币+1): 感谢指出! 2011-05-04 21:55:55
缺心眼桌(金币+1): 谢谢说明 2011-05-04 22:18:17
引用回帖:
Originally posted by 缺心眼桌 at 2011-05-04 15:51:22:
我的体系很大,所以设置了EDIFFG=-0.05(参考了体系差不多的文献),我也只用了一个k点的,这是我当时优化结构时设置的,静态自洽计算时没改。手册里说如果需要准确的能量也必须再做一次静态自洽计算的,再就是 ...

EDIFFG=-0.05,一个k点的设置是很粗糙的,你应该设置精确度高一点再弛豫一遍,试想你优化的构型都很粗糙怎么能指望一次静态计算得到可靠的能量结果呢?最后收敛肯定取决于收敛条件最苛刻的那个设置。
4楼2011-05-04 21:54:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

缺心眼桌

木虫 (正式写手)

引用回帖:
Originally posted by buct2010 at 2011-05-04 21:54:59:
EDIFFG=-0.05,一个k点的设置是很粗糙的,你应该设置精确度高一点再弛豫一遍,试想你优化的构型都很粗糙怎么能指望一次静态计算得到可靠的能量结果呢?最后收敛肯定取决于收敛条件最苛刻的那个设置。

谢谢。不过我的体系实在是很大,k点多了优化不动。而且我看用pwscf的经常就是用一个点优化结构啊。
静态计算的时候,也就是说EDIFFG还是要设置的吗?
5楼2011-05-04 22:20:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

缺心眼桌

木虫 (正式写手)

引用回帖:
Originally posted by buct2010 at 2011-05-04 21:54:59:
EDIFFG=-0.05,一个k点的设置是很粗糙的,你应该设置精确度高一点再弛豫一遍,试想你优化的构型都很粗糙怎么能指望一次静态计算得到可靠的能量结果呢?最后收敛肯定取决于收敛条件最苛刻的那个设置。

静态计算的时候,EDIFFG的设置有没有必要呢?设置与否,对能量结果有何影响?
6楼2011-05-04 22:23:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buct2010

金虫 (正式写手)

【答案】应助回帖

缺心眼桌(金币+1): 2011-05-05 08:42:00
引用回帖:
Originally posted by 缺心眼桌 at 2011-05-04 22:23:35:
静态计算的时候,EDIFFG的设置有没有必要呢?设置与否,对能量结果有何影响?

EDIFFG当然可设可不设,不设就是默认呗,还有静态计算就是为了得CHARGCAR,能量是优化得到的啊。
7楼2011-05-04 23:17:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

缺心眼桌

木虫 (正式写手)

引用回帖:
Originally posted by buct2010 at 2011-05-04 23:17:01:
EDIFFG当然可设可不设,不设就是默认呗,还有静态计算就是为了得CHARGCAR,能量是优化得到的啊。

我知道不设就是默认值,我的意思是,静态自洽计算仅仅是电子步的迭代而已,就一个离子步,原子都不动的啊,那设置EDIFFG为负数,还有意义吗?
另外,手册里说能量是要优化完后做一次静态计算取的,虽说差别可能不大,但这是正规做法。
8楼2011-05-05 08:48:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buct2010

金虫 (正式写手)

【答案】应助回帖

★ ★
youzhizhe(金币+2): 谢谢提示。 2011-05-05 10:12:59
缺心眼桌(金币+1): 谢谢。你的精度太高了,我的体系这样的话算不动。我试着再增加下吧 2011-05-05 10:35:53
引用回帖:
Originally posted by 缺心眼桌 at 2011-05-05 08:48:23:
我知道不设就是默认值,我的意思是,静态自洽计算仅仅是电子步的迭代而已,就一个离子步,原子都不动的啊,那设置EDIFFG为负数,还有意义吗?
另外,手册里说能量是要优化完后做一次静态计算取的,虽说差别可 ...

你说的对:静态自洽计算仅仅是电子步的迭代而已,就一个离子步(影响很小,或者说很鸡肋),原子都不动。所以你不能指望提高静态计算的精度来得到正确的能量,可见优化时的精度是很重要的,我才建议你提高优化的精度来得到正确的能量。静态计算设置的EDIFFG应该和优化时的一样,反正就一步,静态计算最重要的是得到CHARGCAR。
总之:优化的精度高,能量越准确。
看我算的体系:
优化:E0= -.46325923E+03
静态:E0= -.46325932E+03
  INCAR 文件:
   ISTART=0
   ICHARG=2
   ENCUT=400
   EDIFF=0.1E-06
   EDIFFG=0.1E-03
   NSW=200   
   IBRION=2   
   POTIM=0.1
   PREC=Accurate

» 本帖已获得的红花(最新10朵)

9楼2011-05-05 10:07:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

缺心眼桌

木虫 (正式写手)

引用回帖:
Originally posted by buct2010 at 2011-05-05 10:07:37:
你说的对:静态自洽计算仅仅是电子步的迭代而已,就一个离子步(影响很小,或者说很鸡肋),原子都不动。所以你不能指望提高静态计算的精度来得到正确的能量,可见优化时的精度是很重要的,我才建议你提高优化 ...

我怎么觉得那么纠结呢,EDIFFG明明是原子弛豫的标准,静态自洽计算的时候原子都不动,那用这个参数来干什么呀?
10楼2011-05-05 10:45:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 缺心眼桌 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 延安大学回应“副教授留绝笔信”:确已离世,正核实举报内容 +3 yexuqing 2026-11-04 5/250 2026-11-07 13:22 by 苏东坡二世
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 16:32 by 5sJFDKp2kwuI
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ELBZ6WwEaRj6 2026-11-04 14/700 2026-11-06 16:32 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-06 13:47 by BkA9Bi5cMhkD
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:36 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-05 10/500 2026-11-06 12:32 by BkA9Bi5cMhkD
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 04:47 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
信息提示
请填处理意见