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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 明明达到了收敛标准，却还在计算
51/1返回列表
查看: 1457  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yysskk

银虫 (小有名气)

[求助] 明明达到了收敛标准，却还在计算 已有1人参与

今天在弛豫一个体系的时候，无意中发现明明早已达到收敛标准了，计算却还在继续。
这是INCAR文件：
ISTART = 0 ; ICHARG= 2
ENCUT = 400
ISMEAR = 2 ; SIGMA = 0.2
dynamic:
IBRION = 1
NSW = 60
POTIM = 0.2
EDIFFG = -0.05
LREAL = A
ALGO=Fast
NSIM=4
NCORE=24
收敛标准是力小于0.05eV
计算在第25步达到收敛标准终止，这是OUTCAR中第20步的结果：
step: 20
     9.46384       7.01935       4.17889         0.00326      -0.00108      -0.00008
    -2.71429       4.70215       4.17596         0.00494      -0.00182      -0.00012
    12.16736       2.34078       4.18754         0.00222      -0.00008      -0.00075
    -1.36790       7.01668       4.18510         0.00350      -0.00327       0.00066
    -0.05056       4.70705       4.16075        -0.00173       0.00096      -0.00206
    -0.01194       0.02207       4.18084        -0.00111       0.00130       0.00214
     1.31064       7.06307       4.17596         0.00213      -0.00165      -0.00580
     1.33486       2.31176       4.17722         0.00019      -0.00219       0.00119
     2.69068       0.00669       4.17737         0.00416      -0.00086       0.00015
     2.65993       4.66485       4.30412         0.00050      -0.00111       0.00028
     4.06083       2.32909       4.18447        -0.00073      -0.00346       0.00002
     4.06516       7.03286       4.18795        -0.00056       0.00212       0.00197
     5.43410       4.68419       4.17705         0.00031       0.00112       0.00459
     5.40923       0.01108       4.17625         0.00113      -0.00065       0.00248
     6.76245       2.33950       4.17716         0.00157      -0.00248       0.00280
     6.76425       7.02348       4.18846         0.00104       0.00034       0.00590
     8.10727       0.01938       4.17374        -0.00038      -0.00301      -0.00017
     8.11014       4.69639       4.17638         0.00089      -0.00236       0.00479
     9.47190       2.33668       4.17668        -0.00075       0.00036       0.00102
    10.81960       0.00401       4.17358        -0.00140       0.00045       0.00189
     8.10158       7.81695       2.08881        -0.00076       0.00045       0.00001
     9.45834       5.47559       2.08584         0.00150      -0.00020       0.00088
    -2.71358       7.82110       2.08338         0.00161       0.00036       0.00296
    10.81672       3.13254       2.08894         0.00151      -0.00058       0.00362
    -1.36258       5.47446       2.09050         0.00265      -0.00028       0.00764
    12.16647       0.78915       2.08912         0.00102      -0.00024       0.00154
    -0.00961       3.12356       2.07565         0.00035       0.00058       0.00685
    -0.01004       7.82529       2.09491        -0.00048       0.00152       0.00579
     1.35284       5.48351       2.08345         0.00013      -0.00046       0.00267
     1.34166       0.79285       2.09358         0.00064       0.00051       0.00303
     2.69559       3.15225       2.09265         0.00113      -0.00282       0.00124
     2.70090       7.82003       2.08123         0.00289       0.00210       0.00194
     4.05866       5.47888       2.07023         0.00368       0.00065      -0.00246
     4.04886       0.79908       2.09502         0.00125      -0.00149      -0.00024
     5.40228       7.81472       2.09643         0.00218       0.00043      -0.00091
     5.40089       3.13419       2.08366         0.00255      -0.00142      -0.00177
     6.75627       0.79039       2.08640         0.00132      -0.00163       0.00093
     6.75654       5.47721       2.10119         0.00249      -0.00157      -0.00670
     8.11037       3.13243       2.08893         0.00064      -0.00070      -0.00254
     9.46172       0.79076       2.08680         0.00058      -0.00227       0.00135
    -4.05870       7.02990       0.00000        -0.05658       0.07477      -0.00860
    -2.70580       4.68660       0.00000        -0.03048      -0.06867      -0.04895
    -1.35290       7.02990       0.00000        -0.05607       0.09248      -0.03840
    -1.35290       2.34330       0.00000        -0.04192       0.09206      -0.09107
     0.00000       0.00000       0.00000        -0.01884      -0.05779      -0.00963
     0.00000       4.68660       0.00000        -0.01839      -0.04547       0.04743
     1.35290       2.34330       0.00000        -0.01700      -0.02258      -0.02606
     1.35290       7.02990       0.00000         0.03824      -0.03998      -0.06070
     2.70580       4.68660       0.00000         0.05156       0.02615       0.31155
     2.70580       0.00000       0.00000         0.04280      -0.06079       0.01474
     4.05870       7.02990       0.00000        -0.04736       0.04942      -0.02004
     4.05870       2.34330       0.00000        -0.00997       0.02566      -0.03565
     5.41160       4.68660       0.00000         0.00947       0.01142      -0.02814
     5.41160       0.00000       0.00000        -0.06625       0.00504      -0.13090
     6.76450       2.34330       0.00000        -0.00196       0.08162       0.01979
     6.76450       7.02990       0.00000        -0.00243      -0.00834      -0.09579
     8.11740       0.00000       0.00000         0.03133      -0.04949      -0.03659
     8.11740       4.68660       0.00000         0.01062      -0.05799       0.02902
     9.47030       2.34330       0.00000        -0.05510       0.04946      -0.01725
    10.82320       0.00000       0.00000        -0.03716      -0.05629       0.07876
     1.27989       5.11303       6.94810         0.03791       0.01757       0.00890
     3.39901       4.65296       8.68506         0.01873      -0.00449       0.00277
     3.27976       3.95499       7.32438         0.01314      -0.00630       0.01211
     2.39192       4.65821       6.33101         0.02605       0.00594       0.01073
     4.03176       4.05718       9.35500         0.01971       0.00588       0.00035
     3.85154       5.64880       8.58061         0.00708      -0.00725      -0.00048
     2.41455       4.77243       9.15292         0.02689      -0.00164       0.00806
     2.84688       2.94737       7.45616         0.00482      -0.00218       0.01336
     4.26642       3.80759       6.85854         0.00994      -0.01724       0.01615
     0.69760       5.58606       6.27939         0.02434       0.01690       0.00337
     3.05540       5.70801       5.72041         0.04079      -0.02348       0.02443
  free  energy   TOTEN  =      -567.724581 eV
  energy  without entropy=     -567.781882  energy(sigma->0) =     -567.738906
只有中间Z轴坐标为0的原子是非弛豫原子，其余所有原子早已达到收敛标准，为了防止我眼花看错了，我还特意写程序帮我查找，都没有发现弛豫原子上的力超过0.05ev的。
OUTCAR文件开头的参数里也写明了   EDIFFG = -.5E-01   stopping-criterion for IOM
POSCAR里面弛豫和非弛豫原子的标记我也检查过了没问题。
想问一下是我对收敛标准理解错了么？还有其它什么限制条件？
回复此楼

» 猜你喜欢
1楼 2015-01-24 20:26:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yysskk

银虫 (小有名气)
引用回帖:
2楼: Originally posted by ggdh at 2015-01-24 21:24:21
5.41160       0.00000       0.00000        -0.06625       0.00504      -0.1309
    -4.05870       7.02990       0.00000        -0.05658       0.07477      -0.00860
    -1.35290       7.02990        ...

那是非弛豫原子，不会计入收敛标准
赞一下 回复此楼
3楼2015-01-24 21:37:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

ggdh

专家顾问 (正式写手)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
ljw4010: 金币+2, 谢谢指导解答！ 2015-01-26 21:47:00
5.41160       0.00000       0.00000        -0.06625       0.00504      -0.1309
    -4.05870       7.02990       0.00000        -0.05658       0.07477      -0.00860
    -1.35290       7.02990       0.00000        -0.05607       0.09248      -0.03840
9.47030       2.34330       0.00000        -0.05510       0.04946      -0.01725

这有一堆超过0.05的
赞一下 回复此楼
2楼2015-01-24 21:24:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ggdh

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
yysskk: 金币+50, ★★★★★最佳答案 2015-01-24 21:53:36
ljw4010: 金币+2, 谢谢指导解答！ 2015-01-26 21:47:16
引用回帖:
2楼: Originally posted by ggdh at 2015-01-24 21:24:21
5.41160       0.00000       0.00000        -0.06625       0.00504      -0.1309
    -4.05870       7.02990       0.00000        -0.05658       0.07477      -0.00860
    -1.35290       7.02990        ...

我没有仔细审题
吧所有Z！=0的行的x y z 方向上受力平方和开根号后可得
0.00343517102922
0.00526596619814
0.00234463216731
0.00483513184929
0.00285623878554
0.00273892314606
0.00639526387259
0.00249965997688
0.00425061172068
0.00124919974384
0.00353622680268
0.00294769401397
0.00473482840238
0.00280174945347
0.00405651328113
0.00600059997
0.00303865101649
0.00541348316706
0.00131624465811
0.00239470248674
0.000883289307079
0.00175054277297
0.00338870181633
0.00396495901618
0.00809138430678
0.0018626862323
0.0068834148502
0.00600540589802
0.00271245276457
0.00313856655179
0.00328129547588
0.00406518142277
0.00447398033076
0.0019596428246
0.00240112473645
0.00341347330442
0.00229438444904
0.00731812817598
0.00271130964665
0.00270403402345
0.0427209901571
0.0194588257611
0.0189473401827
0.0287926900445
0.0205713635912
0.0101449149824
0.028119838193
0.0143692170977
0.0256289621327
0.0298228519763
0.0530279115184
我想这应该是判断标准
赞一下 回复此楼
4楼2015-01-24 21:39:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyqh

铁杆木虫 (著名写手)
是不是 能量达到收敛标准，而力还没有达到收敛标准呢？
赞一下 回复此楼
雨后方能见彩虹
5楼2015-01-25 01:18:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +3 芦lty 2026-03-25 4/200 2026-03-25 23:25 by 芦lty
[考研] 调剂310 +3 温柔的晚安 2026-03-25 4/200 2026-03-25 23:16 by peike
[考研] 化学工程085602 305分求调剂 +9 RichLi_ 2026-03-25 9/450 2026-03-25 22:21 by 418490947
[考研] 292求调剂 +6 鹅鹅鹅额额额额� 2026-03-25 7/350 2026-03-25 22:04 by 无际的草原
[考研] 材料求调剂 +4 .m.. 2026-03-25 4/200 2026-03-25 21:30 by peike
[考研] 290分调剂求助 +3 吉祥止止陈 2026-03-25 3/150 2026-03-25 19:58 by barlinike
[考研] 材料与化工考研调剂 +9 孅華 2026-03-22 9/450 2026-03-25 13:09 by cmz0325
[考研] 359求调剂 +3 王了个楠 2026-03-25 3/150 2026-03-25 12:50 by Dyhoer
[考研] 材料学硕，求调剂 6+4 糖葫芦888ll 2026-03-22 9/450 2026-03-25 11:19 by greychen00
[考研] 307求调剂 +5 超级伊昂大王 2026-03-24 5/250 2026-03-24 15:46 by 星空星月
[考研] 321求调剂 +4 Ymlll 2026-03-24 4/200 2026-03-24 14:44 by sprinining
[考研] 341求调剂(一志愿湖南大学070300) +5 番茄头--- 2026-03-22 6/300 2026-03-23 23:45 by Txy@872106
[考研] 269求调剂 +4 我想读研11 2026-03-23 4/200 2026-03-23 21:25 by pswait
[考研] 一志愿重庆大学085700资源与环境，总分308求调剂 +7 墨墨漠 2026-03-23 8/400 2026-03-23 20:36 by Creta
[考研] 293求调剂 +3 涛涛Wjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[考研] 材料与化工085600，总分304，本科有两篇sci参与，求调剂 +4 幸运的酱酱 2026-03-22 5/250 2026-03-22 20:15 by edmund7
[考研] 一志愿华中农业071010，总分320求调剂 +5 困困困困坤坤 2026-03-20 6/300 2026-03-22 17:41 by hxsm
[基金申请] 学校已经提交到NSFC，还能修改吗？ 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 求调剂 +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
信息提示
关闭
请填处理意见
关闭