版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

飞翔……

金虫 (正式写手)

应助: 4 (幼儿园)
金币: 1378.8
散金: 160
红花: 6
沙发: 1
帖子: 606
在线: 269.8小时
虫号: 673866
注册: 2008-12-13
专业: 催化化学

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

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

回复此楼

1楼2015-11-17 11:13:20

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

飞翔……

金虫 (正式写手)

应助: 4 (幼儿园)
金币: 1378.8
散金: 160
红花: 6
沙发: 1
帖子: 606
在线: 269.8小时
虫号: 673866
注册: 2008-12-13
专业: 催化化学

引用回帖:

2楼: Originally posted by 未岸亡灵 at 2015-11-17 13:30:39
检查OSZICAR，看电子步是否收敛。有可能你自洽的时候电子步并未收敛

谢谢你的回答！
检查了静态自洽时候的电子步，确实未收敛，请问有什么办法可以解决吗？
（另外，我也检查了结构优化和静态自洽的总能相一致的情况，发现静态自洽时，电子步均未收敛，所以我在想是不是其他原因导致的）

赞一下

回复此楼

3楼2015-11-17 14:37:12

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

飞翔……

金虫 (正式写手)

应助: 4 (幼儿园)
金币: 1378.8
散金: 160
红花: 6
沙发: 1
帖子: 606
在线: 269.8小时
虫号: 673866
注册: 2008-12-13
专业: 催化化学

送红花一朵

引用回帖:

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

飞翔……

金虫 (正式写手)

应助: 4 (幼儿园)
金币: 1378.8
散金: 160
红花: 6
沙发: 1
帖子: 606
在线: 269.8小时
虫号: 673866
注册: 2008-12-13
专业: 催化化学

引用回帖:

6楼: Originally posted by 未岸亡灵 at 2015-11-17 22:20:04
我遇到过同样的情况。
INCAR加个参数，NELM=200，看看能不能收敛。
或者参照http://blog.sina.com.cn/s/blog_553461e901013oo3.html...

非常感谢你的帮助！
首先通过你给的办法，让我解决了这个问题；
然后，通过这两天你的帮助，让我对静态自洽有了进一步的认识；

赞一下

回复此楼

7楼2015-11-18 22:26:08

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