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

登录注册

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

24小时热门版块排行榜

南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知

返回列表

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

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

1ST强帖: 5
应助: 46 (小学生)
贵宾: 8.842
金币: 7440.4
散金: 9654
红花: 118
沙发: 23
帖子: 5306
在线: 1961.4小时
虫号: 537452
注册: 2008-04-01
性别: GG
专业: 理论和计算化学
管辖: 仿真模拟

[交流] 【求助成功】为什么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 ]

回复此楼

» 猜你喜欢

孩子确诊有中度注意力缺陷已经有13人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复

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

为什么 Johnson wiley 中chirality的数据库有部分文章是可以免费下的呢？已经有4人回复
求助：PLGA合成中要求压力较低，1000pa，怎么控制压力不变呢？已经有10人回复
紫外分光光度计和原子吸收分光光度计单色器位置为什么不同？已经有7人回复
为什么液相色谱压力一直不稳？已经有6人回复
Required Reviews Completed 变为ender review 又变为Required Reviews Completed 已经有12人回复
做上环手术弄得大出血，医院却说正常已经有4人回复
★★★教你图片如何处理才能达到SCI投稿要求已经有61人回复
请问这个数据是否符合要求已经有5人回复
required reviews completed已经一个月已经有19人回复
等待Required Reviews Completed 已经有27人回复
SBB 一个月Required Reviews Completed 已经有12人回复

1楼2010-05-07 09:18:42

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

fanchen021

铁杆木虫 (正式写手)

应助: 7 (幼儿园)
金币: 7994.9
散金: 1844
红花: 4
帖子: 776
在线: 569小时
虫号: 773188
注册: 2009-05-18
专业: 理论和计算化学

★
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 个回答

jinhuan123

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 529.3
散金: 21
红花: 1
帖子: 106
在线: 51.7小时
虫号: 805643
注册: 2009-07-08
专业: 催化化学

★
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

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 529.3
散金: 21
红花: 1
帖子: 106
在线: 51.7小时
虫号: 805643
注册: 2009-07-08
专业: 催化化学

★
gavinliu7390(金币+1):感谢交流！ 2010-05-07 10:46:46

OUTCAR 里确实有EDIFFG，但是具体怎么体现出来，就像你列出的这个，现在还没有研究过呢，期待......

赞一下(1人)

回复此楼

3楼2010-05-07 10:10:01

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

citrine

铁杆木虫 (正式写手)

应助: 0 (幼儿园)
金币: 8823.3
散金: 1
红花: 7
帖子: 460
在线: 299.9小时
虫号: 427638
注册: 2007-07-29
专业: 凝聚态物性 II ：电子结构

★ ★
cenwanglai(金币+2):谢谢参与.按你说的意思,是要判断每一个原子的受力的哈.但是不知道为什么我的任务会在力没有达到收敛精度的情况下报告"达到收敛精度,退停止结构能量最小化" 2010-05-07 20:59:32
gavinliu7390(金币+2):谢谢交流！ 2010-05-10 08:48:43

POSITION TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
2.20903 0.00000 6.13890 0.000000 0.000000 0.055890

最后三个数是原子收到的力，显然没有达到收敛标准。
我注意到NSW=10，这个参数应该变大，比如说100

赞一下(1人)

回复此楼

5楼2010-05-07 11:27:40

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

查看全部 12 个回答