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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 静态自洽计算能量时EDIFFG的问题
51/1返回列表
查看: 3524  |  回复: 14
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 1ST强帖 ，点击这里进行查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

缺心眼桌

木虫 (正式写手)

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

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

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

VASP

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2011-05-04 09:48:46
已阅   回复此楼   关注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
赞一下(1人) 回复此楼

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

1051358187
9楼2011-05-05 10:07:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

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的误差是很大的）。
赞一下(1人) 回复此楼
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点的设置是很粗糙的，你应该设置精确度高一点再弛豫一遍，试想你优化的构型都很粗糙怎么能指望一次静态计算得到可靠的能量结果呢？最后收敛肯定取决于收敛条件最苛刻的那个设置。
赞一下(1人) 回复此楼
4楼2011-05-04 21:54:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭