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

24小时热门版块排行榜    

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

飞翔……

金虫 (正式写手)

[求助] 静态自洽计算的能量 与 结构优化时的总能 不一致已有1人参与

研究的金属氧化物体系
1. 采进行结构优化，计算收敛后，提取energy  without entropy的最后一项作为体系的总能；
2. 修改NSW, IBRION后进行静态自洽计算，得出的总能和上面计算的相差很大，但结构并没有发生变化；
很困惑为什么会出现这种情况，是什么原因造成的？
INCAR在附件中，vasp采用的是5.2的版本。
走过路的虫友都来瞧一瞧，看一看，帮忙出言献策！
回复此楼

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

soft study

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2015-11-17 11:13:20
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

飞翔……

金虫 (正式写手)
引用回帖:
2楼: Originally posted by 未岸亡灵 at 2015-11-17 13:30:39
检查OSZICAR，看电子步是否收敛。有可能你自洽的时候电子步并未收敛

谢谢你的回答！
检查了静态自洽时候的电子步，确实未收敛，请问有什么办法可以解决吗？
（另外，我也检查了结构优化和静态自洽的总能相一致的情况，发现静态自洽时，电子步均未收敛，所以我在想是不是其他原因导致的）
赞一下 回复此楼
3楼2015-11-17 14:37:12
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

未岸亡灵

铁杆木虫 (著名写手)

暂时离开小木虫

【答案】应助回帖

★
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-11-21 08:10:22
检查OSZICAR，看电子步是否收敛。有可能你自洽的时候电子步并未收敛
赞一下(1人) 回复此楼
2楼2015-11-17 13:30:39
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

未岸亡灵

铁杆木虫 (著名写手)

暂时离开小木虫

【答案】应助回帖

★
liliangfang: 金币+1, 谢谢交流 2015-11-18 08:01:33
引用回帖:
3楼: Originally posted by 飞翔…… at 2015-11-17 14:37:12
谢谢你的回答！
检查了静态自洽时候的电子步，确实未收敛，请问有什么办法可以解决吗？
（另外，我也检查了结构优化和静态自洽的总能相一致的情况，发现静态自洽时，电子步均未收敛，所以我在想是不是其他原因导 ...

解决办法自行百度，有很多谈如何收敛的。
NELM默认值是60，有些计算大约是能量不断递减的，故到最后虽然没有收敛，但是能量值和优化时得到差异不大也很正常，因为还是有可能再跑几步就达到你的收敛标准了。
而你现在的情况，很可能是遇到的是能量震荡的过程，并且震荡幅度还比较大，OSZICAR里dE能量一会儿加一会儿减的，自然无法收敛，即使设置再多步，也不可能收敛。这时候就要考虑换算法什么的了。

当然，不收敛的结果都是不可信的，因为算出来的本身就不是个绝对能量，能量收敛到的精度还不一样，自然在做差的时候误差会很大
赞一下(1人) 回复此楼

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

飞翔……
4楼2015-11-17 17:51:15
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖

飞翔……

金虫 (正式写手)
送红花一朵
引用回帖:
4楼: Originally posted by 未岸亡灵 at 2015-11-17 17:51:15
解决办法自行百度，有很多谈如何收敛的。
NELM默认值是60，有些计算大约是能量不断递减的，故到最后虽然没有收敛，但是能量值和优化时得到差异不大也很正常，因为还是有可能再跑几步就达到你的收敛标准了。
而你 ...

非常感谢你抽空回复！
查了如何收敛，谈的都是结构优化这一步，但是结构优化这一部，我看了我计算的结构，虽然第一个离子步中，电子步没有收敛，但是第二个离子步之后，电子就能收敛至标准内。（不知道这一步是否对结果又影响？有人说没有，我的理解也是没有）
现在的情况是结构优化的时候可以收敛，此时的电子步和离子步都没有很大的震荡，但若将此时的CONTCAR进行静态计算，60个电子步内却又不能收敛，所以很困惑。（优化的时候已经收敛的结构，进行静态自洽计算时却又不收敛）
赞一下 回复此楼
5楼2015-11-17 20:25:35
 已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答
信息提示
关闭
请填处理意见
关闭