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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助成功】为什么total drift没有达到EDIFFG的要求,却说reached required accuracy?
52/1返回列表
查看: 2420  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

[交流] 【求助成功】为什么total drift没有达到EDIFFG的要求,却说reached required accuracy? 已有8人参与

EDIFFG=-0.001,呵呵,比较猛吧.
OUTCAR后面报告"reached required accuracy-stoping structural energy minimisation".

但是TOTAL-FORCE和total drift都没有达到EDIFFG的要求.PREC=Accurate,关闭实空间后,每个原子每个方向上的力和total drift还是没有达到要求,计算就退出来了,运行离子步小于NSW.

这是什么原因呢?编译错了吗?

POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      2.20903      0.00000      6.13890         0.000000      0.000000      0.055890
      0.00000      0.00000      0.00000         0.000000      0.000000      0.401501
      1.10452      1.91308      3.08095         0.000000      0.000000     -0.139539
      2.20903      0.00000      1.56202         0.000000      0.000000      0.013240
      0.00000      0.00000      4.65030         0.000000      0.000000      0.044110
      1.10452      1.91308      7.62370         0.000000      0.000000      0.036490
      5.52258     -1.91308      6.13890         0.000000      0.000000      0.055890
      3.31355     -1.91308      0.00000         0.000000      0.000000      0.401501
      4.41806      0.00000      3.08095         0.000000      0.000000     -0.139539
      5.52258     -1.91308      1.56202         0.000000      0.000000      0.013240
      3.31355     -1.91308      4.65030         0.000000      0.000000      0.044110
      4.41806      0.00000      7.62370         0.000000      0.000000      0.036490
      2.20903      3.82615      6.13890         0.000000      0.000000      0.055890
      0.00000      3.82615      0.00000         0.000000      0.000000      0.401501
      1.10452      5.73923      3.08095         0.000000      0.000000     -0.139539
      2.20903      3.82615      1.56202         0.000000      0.000000      0.013240
      0.00000      3.82615      4.65030         0.000000      0.000000      0.044110
      1.10452      5.73923      7.62370         0.000000      0.000000      0.036490
      5.52258      1.91308      6.13890         0.000000      0.000000      0.055890
      3.31355      1.91308      0.00000         0.000000      0.000000      0.401501
      4.41806      3.82615      3.08095         0.000000      0.000000     -0.139539
      5.52258      1.91308      1.56202         0.000000      0.000000      0.013240
      3.31355      1.91308      4.65030         0.000000      0.000000      0.044110
      4.41806      3.82615      7.62370         0.000000      0.000000      0.036490
      0.00000      0.00000      6.92021         0.000000      0.000000     -0.138665
      2.20903      0.00000      3.86577         0.000000      0.000000     -0.147815
      1.10452      1.91308      0.78101         0.000000      0.000000     -0.125211
      3.31355     -1.91308      6.92021         0.000000      0.000000     -0.138665
      5.52258     -1.91308      3.86577         0.000000      0.000000     -0.147815
      4.41806      0.00000      0.78101         0.000000      0.000000     -0.125211
      0.00000      3.82615      6.92021         0.000000      0.000000     -0.138665
      2.20903      3.82615      3.86577         0.000000      0.000000     -0.147815
      1.10452      5.73923      0.78101         0.000000      0.000000     -0.125211
      3.31355      1.91308      6.92021         0.000000      0.000000     -0.138665
      5.52258      1.91308      3.86577         0.000000      0.000000     -0.147815
      4.41806      3.82615      0.78101         0.000000      0.000000     -0.125211
-----------------------------------------------------------------------------------
    total drift:                                0.000296     -0.000609     -0.021872

INCAR:

ENCUT = 400
EDIFF = 1E-4
EDIFFG = -0.001
IBRION = 2
ISMEAR = 0      
SIGMA = 0.05   
NSW = 10
ALGO = Fast
#LREAR = Auto  
POTIM = 0.2


呵呵,经过一段时间观察发现:EDIFFG<0是对自由原子(T T T)受力来做判断的,对固定原子(F F F)不做为判断依据.只要自由原子达到力的收敛标准就可以了.

[ Last edited by cenwanglai on 2010-5-21 at 13:37 ]
回复此楼

» 猜你喜欢
1楼 2010-05-07 09:18:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sungjen

银虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
7楼: Originally posted by zengxiangm at 2010-05-07 12:18:40
老师说是每一个原子每一个方向上的受力都要小于|EDIFFG|，但我的情况也和楼主一样，希望有高人来解答

请问你是最后怎么解决的 这个问题
赞一下 回复此楼
希望有一天成为大神
12楼2014-10-04 11:33:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

jinhuan123

金虫 (小有名气)
★
cenwanglai(金币+1):你用什么?呵呵 2010-05-07 10:05:10
gavinliu7390(金币+1):谢谢交流！ 2010-05-10 08:47:44
这个貌似没有用到过哦
赞一下(1人) 回复此楼
2楼2010-05-07 09:50:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jinhuan123

金虫 (小有名气)
★
gavinliu7390(金币+1):感谢交流！ 2010-05-07 10:46:46
OUTCAR 里确实有EDIFFG，但是具体怎么体现出来，就像你列出的这个，现在还没有研究过呢，期待......
赞一下(1人) 回复此楼
3楼2010-05-07 10:10:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanchen021

铁杆木虫 (正式写手)
★
cenwanglai(金币+2):呵呵,怎么运行呢?我在运行结束后含有OUTCAR的目录下运行$vef.pl,生成的fe.dat是空的呢. 2010-05-07 20:56:45
gavinliu7390(金币+1):谢谢交流！ 2010-05-10 08:48:24
用vef.pl  看一下啊
赞一下(1人) 回复此楼
4楼2010-05-07 10:11:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +17 木水思豆 2026-08-25 20/1000 2026-08-26 01:26 by chengyan1220
[基金申请] 今天务委会开完了，明天出结果吗 +18 angus9576 2026-08-25 22/1100 2026-08-26 00:41 by merchancy
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 今天基金会出结果吗？20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] filecode，4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤？
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭