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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助成功】为什么total drift没有达到EDIFFG的要求,却说reached required accuracy?
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
121/2返回列表12下一页
查看: 1874  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

fanchen021

铁杆木虫 (正式写手)
★
cenwanglai(金币+3): 2010-05-10 18:11:52
qasd(金币+1):谢谢 2010-05-10 22:58:57
那你试试 grep 'FORCES: max atom, RMS' OUTCAR
赞一下(2人) 回复此楼
11楼2010-05-10 13:54:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

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的回帖

citrine

铁杆木虫 (正式写手)
★ ★
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的回帖

stractor

金虫 (著名写手)
★
cenwanglai(金币+2):呵呵,即使设置成EDIFFG=-0.001,即使实际上力没有达到要求的收敛精度,但是还是会以"reached required accuracy"结束.不知道是不是编译的原因. 2010-05-07 20:53:45
gavinliu7390(金币+1):谢谢交流！ 2010-05-10 08:48:58
EDIFFG=-0.001,收敛标准确实挺高，建议改为-0.01至-0.03.
赞一下(1人) 回复此楼
6楼2010-05-07 12:10:18
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zengxiangm

木虫 (小有名气)
老师说是每一个原子每一个方向上的受力都要小于|EDIFFG|，但我的情况也和楼主一样，希望有高人来解答
赞一下 回复此楼
7楼2010-05-07 12:18:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★
cenwanglai(金币+2):谢谢参与.grep "F" OSZICAR|tail -1给出的F就是总能? 2010-05-07 20:55:08
gavinliu7390(金币+1):感谢交流！ 2010-05-07 23:27:39
EDIFF为负数时，以原子所受的力来判断收敛；如果为正数，以前后两次总能的变化为依据来判断收敛。
赞一下(1人) 回复此楼
8楼2010-05-07 14:11:45
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
★ ★ ★ ★
cenwanglai(金币+1):谢谢参与
cenwanglai(金币+1): 设置为-0.001纯粹是跑一下看看会出现什么情况.呵呵 2010-05-08 13:43:52
gavinliu7390(金币+2):谢谢交流！ 2010-05-10 08:49:15
zzy870720z(金币+1):感谢交流 2010-05-10 08:49:44
此外，比如你优化的步数如果设置为20步，NSW=20，到第20步结束时，哪怕力没有达到标准，VASP也会reached required accuracy而结束。
你为什么要设置-0.001？这个值太小了，费时。-0.01至-0.03够了。
赞一下(3人) 回复此楼
9楼2010-05-08 13:27:57
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主
这个体系收敛没有问题了,问题是收敛的结果不对.
vef.pl  怎么运行呢?运行vasp后,直接在工作目录下的命令提示符下输入?
$vef.pl
有文件fe.dat生成,但是是空的文件.
赞一下 回复此楼
10楼2010-05-10 08:42:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cenwanglai 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭