24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1999  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

superdonnie

铜虫 (小有名气)

[求助] VASP结构优化的收敛问题 已有3人参与

各位大侠大家好,

我使用vasp做hcp铁的结构优化,发现力的收敛存在问题,希望大侠们能帮我看看是怎么回事。

我做结构优化时INCAR中使用了EDIFF=1E-8; EDIFFG=-1E3,从OUTCAR中提取的力的情况如下:

POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      0.00000      0.00000      0.00000         0.000000      0.000000      0.000000
      0.00000      1.41908      2.00688         0.000000      0.000000      0.000000
-----------------------------------------------------------------------------------
--
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      0.00000      0.00000      0.00000         0.000000      0.000000      0.000000
      0.00000      1.44122      1.94569         0.000000      0.000000      0.000000
-----------------------------------------------------------------------------------
--
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      0.00000      0.00000      0.00000         0.000000      0.000000      0.000000
      0.00000      1.44245      1.94239         0.000000      0.000000      0.000000
-----------------------------------------------------------------------------------
--
POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      0.00000      0.00000      0.00000         0.000000      0.000000      0.000000
      0.00000      1.44283      1.94136         0.000000      0.000000      0.000000
-----------------------------------------------------------------------------------

我发现在结构优化的过程中所有的离子构型的力都是0,按照之前ediffg=-1E3的标准,计算早在初始离子构型就已经达到了收敛标准,那么为什么还会继续进行这么多步呢?是否有什么地方我可能给漏掉了?

还望各位大侠们不吝赐教。
回复此楼

» 猜你喜欢

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

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

superdonnie

铜虫 (小有名气)

引用回帖:
5楼: Originally posted by xrdcvw73 at 2015-08-26 21:50:18
个人觉得精度过高
一般觉得1E-4   -1E-2 就ok
很多教程都是这个精度,如果你需要提高精度,可以尝试多步骤优化,先打到较低精度,然后一步一步提高

我使用的IBRION=2,也就是比较robust的conjugate gradient算法,根据输出可以达到1E-8的总能量精度。而力的精度即使使用-1E-2,也会得出同样的结果,也就是明明看起来一开始就达到了收敛水平,但还是又优化了两次结构。谷歌搜索也没发现有人遇见过同样的问题,真是头疼啊
7楼2015-08-26 21:56:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

bird007

荣誉版主 (职业作家)

踩单车的季节

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
superdonnie: 金币+2, 有帮助, 谢谢交流~ 2015-08-27 09:29:18
liliangfang: 金币+1, 谢谢交流 2015-08-29 17:37:26
EDIFF=1E-8
这个精度 太高了吧
一日,老师在课堂上教偶认字——帅.偶百思不得其解,正在这时,同桌mm递来一面镜子,偶照之,顿时恍然大悟......[此广告位招商!][此广告位招商!]
2楼2015-08-26 19:43:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Litter_horsr

银虫 (小有名气)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
superdonnie: 金币+2, 有帮助, 谢谢交流~ 2015-08-27 09:29:37
liliangfang: 金币+1, 谢谢交流 2015-08-29 17:37:32
1E-8这个估计有点悬。
-1E3这个写法估计也不对
想好好学习的,一般妒能
3楼2015-08-26 21:24:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

superdonnie

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by bird007 at 2015-08-26 19:43:24
EDIFF=1E-8
这个精度 太高了吧

因为要做声子计算,所以按照计算声子的phonopy程序对VASP输入的要求将其精度取得较高,从输出来看计算的能量是达到了这样的精度的。感谢大侠参与讨论
4楼2015-08-26 21:46:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 279学硕食品专业求调剂院校 20+7 孤独的狼爱吃羊 2026-04-12 29/1450 2026-04-16 09:00 by screening
[考研] 294求调剂 +10 淡然654321 2026-04-15 10/500 2026-04-16 08:58 by screening
[考研] 327求调剂 +25 Xxjc1107. 2026-04-13 27/1350 2026-04-15 23:22 by Equinoxhua
[考研] 322求调剂 +7 123安康 2026-04-12 14/700 2026-04-15 21:48 by noqvsozv
[考研] 一志愿华中农业071010,320求调剂 +11 困困困困坤坤 2026-04-14 12/600 2026-04-15 18:58 by zhuwenxu
[考研] 生物学调剂 +9 纸扇zhishan 2026-04-13 9/450 2026-04-15 18:28 by AN流800
[考研] 291求调剂 +7 关忆北. 2026-04-14 7/350 2026-04-15 14:39 by 黑科技花岗岩
[考研] 332求调剂 +15 蕉蕉123 2026-04-10 15/750 2026-04-13 23:12 by pies112
[考研] 0856专硕求调剂 希望是a区院校 +24 好好休息好不好 2026-04-09 27/1350 2026-04-13 22:22 by pies112
[考研] 297工科,求调剂? +13 河南农业大学-能 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
[考研] 生物学调剂 +11 小冉要努力 2026-04-10 13/650 2026-04-13 11:46 by 电化学及催化
[考研] 0854调剂 +10 长弓傲 2026-04-11 11/550 2026-04-13 10:38 by wp06
[考研] 电气工程专硕320求调剂 +5 小麻子111 2026-04-10 5/250 2026-04-12 10:47 by zhouyuwinner
[考研] 326求调剂 +6 Shansyn 2026-04-10 6/300 2026-04-12 09:46 by hammer3
[考研] 22408调剂315分 +3 zhuangyan123 2026-04-09 3/150 2026-04-12 00:25 by 蓝云思雨
[考研] 284求调剂 +11 archer.. 2026-04-09 12/600 2026-04-11 20:23 by 蓝云思雨
[考研] 270求调剂 +14 杨乐369 2026-04-11 14/700 2026-04-11 20:16 by 蓝云思雨
[考研] 一志愿985机械学硕380求调剂 +5 关关雎鸠10 2026-04-11 5/250 2026-04-11 10:10 by 知念。A
[考研] 考研调剂 +26 硕星赴 2026-04-09 27/1350 2026-04-10 22:24 by 猪会飞
[考研] 085601初试330分找调剂 +10 流心奶黄包l 2026-04-09 10/500 2026-04-10 08:14 by Sammy2
信息提示
请填处理意见