24小时热门版块排行榜    

查看: 4573  |  回复: 8

飞翔……

金虫 (正式写手)

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

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

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

soft study

» 猜你喜欢

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

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

未岸亡灵

铁杆木虫 (著名写手)

暂时离开小木虫


【答案】应助回帖


感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-11-21 08:10:22
检查OSZICAR,看电子步是否收敛。有可能你自洽的时候电子步并未收敛
2楼2015-11-17 13:30:39
已阅   关注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的回帖

未岸亡灵

铁杆木虫 (著名写手)

暂时离开小木虫


【答案】应助回帖


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

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

当然,不收敛的结果都是不可信的,因为算出来的本身就不是个绝对能量,能量收敛到的精度还不一样,自然在做差的时候误差会很大

» 本帖已获得的红花(最新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的回帖

未岸亡灵

铁杆木虫 (著名写手)

暂时离开小木虫


【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
liliangfang: 金币+1, 谢谢交流 2015-11-18 08:01:24
飞翔……: 金币+20, ★★★★★最佳答案, 非常非常感谢! 2015-11-18 22:26:45
引用回帖:
5楼: Originally posted by 飞翔…… at 2015-11-17 20:25:35
非常感谢你抽空回复!
查了如何收敛,谈的都是结构优化这一步,但是结构优化这一部,我看了我计算的结构,虽然第一个离子步中,电子步没有收敛,但是第二个离子步之后,电子就能收敛至标准内。(不知道这一步是否 ...

我遇到过同样的情况。
INCAR加个参数,NELM=200,看看能不能收敛。
或者参照http://blog.sina.com.cn/s/blog_553461e901013oo3.html
6楼2015-11-17 22:20:04
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

飞翔……

金虫 (正式写手)

引用回帖:
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的回帖

jimsmart

捐助贵宾 (职业作家)

我建议你向售后工程师索要一个能够一致的数据案例, 再观察.重复他的计算.找感觉.
8楼2015-11-19 13:58:33
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

飞翔……

金虫 (正式写手)

引用回帖:
8楼: Originally posted by jimsmart at 2015-11-19 13:58:33
我建议你向售后工程师索要一个能够一致的数据案例, 再观察.重复他的计算.找感觉.

谢谢你的参与!
你说的是服务器的售后工程师吗?他只对计算机系统很了解,维护计算机的正常运转,但对于我们计算的东西基本上不清楚。
9楼2015-11-19 16:18:46
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 飞翔…… 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:27 by oEVWOejN9taj
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:22 by oEVWOejN9taj
[基金申请] fileCode有新解读? +9 Tide man 2026-08-08 17/850 2026-08-08 14:56 by archvillain
[基金申请] 这样的filecode谁见过 +5 布布和一二 2026-08-08 8/400 2026-08-08 14:54 by archvillain
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
信息提示
请填处理意见