| 查看: 4071 | 回复: 8 | ||||
飞翔……金虫 (正式写手)
|
[求助]
静态自洽计算的能量 与 结构优化时的总能 不一致已有1人参与
|
|
研究的金属氧化物体系 1. 采进行结构优化，计算收敛后，提取energy without entropy的最后一项作为体系的总能； 2. 修改NSW, IBRION后进行静态自洽计算，得出的总能和上面计算的相差很大，但结构并没有发生变化； 很困惑为什么会出现这种情况，是什么原因造成的？ INCAR在附件中，vasp采用的是5.2的版本。 走过路的虫友都来瞧一瞧，看一看，帮忙出言献策！  |
回复此楼» 收录本帖的淘帖专辑推荐
 soft study |
» 猜你喜欢
博士读完未来一定会好吗
已经有21人回复
-
导师想让我从独立一作变成了共一第一
已经有5人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有9人回复
-
申请2026年博士
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
关于总能计算
已经有1人回复- 价带顶到macro average potential 的值
已经有1人回复
- Vasp优化结构后进行一步静态自洽，能量提升是什么原因呢？
已经有2人回复
- VASP杂化泛函计算太密度
已经有3人回复
- VASP计算中KPOINTS文件的作用？
已经有8人回复
- 求教：关于Partical charge density的做法？做过的前辈请进
已经有2人回复
- 总能问题之疑惑
已经有11人回复
- 本人初学第一性原理，一些个人总结与大家分享，欢迎大家前来批评指正，共同学习~~~~~
已经有31人回复
- vasp 关于curie temperature的计算
已经有5人回复
- 如何计算/作图H/H_2^+/H_2/Si的基态/第一激发态/第二激发态的能量/电荷密度
已经有4人回复
- 杂化泛函计算出能带需要的费米能级应该取哪一个？
已经有1人回复
- 【求助】关于费米能级
已经有13人回复
- 静态自洽计算能量时EDIFFG的问题
已经有14人回复
- 【求助完毕】CO极性分子在金属氧化物表面吸附为什么是O上C下的构型？
已经有6人回复
- 【求助】关于拟合得到体弹性模量&自旋极化&LDA+U等问题，望大家指点指点
已经有4人回复
- 【原创】关于结构弛豫设置及收敛判据选择的一点体会
已经有140人回复
- 【求助】VASP加电场的具体计算步骤是怎么样的？
已经有31人回复
- 【求助完毕】硫酸根离子的结构优化，偶极修正及势能修正。
已经有20人回复
- 【求助】J值的算法
已经有8人回复
- 宇宙的存在方式和世界的因果关系
已经有6人回复
未岸亡灵
铁杆木虫 (著名写手)
暂时离开小木虫
- 1ST强帖: 2
- 应助: 622 (博士)
- 金币: 6517.3
- 散金: 1221
- 红花: 111
- 帖子: 1345
- 在线: 517.7小时
- 虫号: 1730773
- 注册: 2012-04-01
- 性别: GG
- 专业: 理论和计算化学
2楼2015-11-17 13:30:39
飞翔……
金虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 1378.8
- 散金: 160
- 红花: 6
- 沙发: 1
- 帖子: 606
- 在线: 269.8小时
- 虫号: 673866
- 注册: 2008-12-13
- 专业: 催化化学
3楼2015-11-17 14:37:12
未岸亡灵
铁杆木虫 (著名写手)
暂时离开小木虫
- 1ST强帖: 2
- 应助: 622 (博士)
- 金币: 6517.3
- 散金: 1221
- 红花: 111
- 帖子: 1345
- 在线: 517.7小时
- 虫号: 1730773
- 注册: 2012-04-01
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★
liliangfang: 金币+1, 谢谢交流 2015-11-18 08:01:33
liliangfang: 金币+1, 谢谢交流 2015-11-18 08:01:33
|
解决办法自行百度，有很多谈如何收敛的。 NELM默认值是60，有些计算大约是能量不断递减的，故到最后虽然没有收敛，但是能量值和优化时得到差异不大也很正常，因为还是有可能再跑几步就达到你的收敛标准了。 而你现在的情况，很可能是遇到的是能量震荡的过程，并且震荡幅度还比较大，OSZICAR里dE能量一会儿加一会儿减的，自然无法收敛，即使设置再多步，也不可能收敛。这时候就要考虑换算法什么的了。 当然，不收敛的结果都是不可信的，因为算出来的本身就不是个绝对能量，能量收敛到的精度还不一样，自然在做差的时候误差会很大 |
» 本帖已获得的红花（最新10朵）
飞翔……4楼2015-11-17 17:51:15
飞翔……
金虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 1378.8
- 散金: 160
- 红花: 6
- 沙发: 1
- 帖子: 606
- 在线: 269.8小时
- 虫号: 673866
- 注册: 2008-12-13
- 专业: 催化化学
5楼2015-11-17 20:25:35
未岸亡灵
铁杆木虫 (著名写手)
暂时离开小木虫
- 1ST强帖: 2
- 应助: 622 (博士)
- 金币: 6517.3
- 散金: 1221
- 红花: 111
- 帖子: 1345
- 在线: 517.7小时
- 虫号: 1730773
- 注册: 2012-04-01
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
liliangfang: 金币+1, 谢谢交流 2015-11-18 08:01:24
飞翔……: 金币+20, ★★★★★最佳答案, 非常非常感谢！ 2015-11-18 22:26:45
liliangfang: 金币+1, 谢谢交流 2015-11-18 08:01:24
飞翔……: 金币+20, ★★★★★最佳答案, 非常非常感谢！ 2015-11-18 22:26:45
|
我遇到过同样的情况。 INCAR加个参数，NELM=200，看看能不能收敛。 或者参照http://blog.sina.com.cn/s/blog_553461e901013oo3.html |
6楼2015-11-17 22:20:04
飞翔……
金虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 1378.8
- 散金: 160
- 红花: 6
- 沙发: 1
- 帖子: 606
- 在线: 269.8小时
- 虫号: 673866
- 注册: 2008-12-13
- 专业: 催化化学
7楼2015-11-18 22:26:08
jimsmart
捐助贵宾 (职业作家)
- 应助: 109 (高中生)
- 贵宾: 0.2
- 金币: 3581.7
- 散金: 100
- 红花: 5
- 帖子: 3110
- 在线: 567.6小时
- 虫号: 57301
- 注册: 2005-01-29
- 性别: GG
- 专业: 光谱分析
8楼2015-11-19 13:58:33
飞翔……
金虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 1378.8
- 散金: 160
- 红花: 6
- 沙发: 1
- 帖子: 606
- 在线: 269.8小时
- 虫号: 673866
- 注册: 2008-12-13
- 专业: 催化化学
9楼2015-11-19 16:18:46