| 查看: 1396 | 回复: 6 | |||||||||||||
[求助]
明明达到了收敛标准，却还在计算已有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里面弛豫和非弛豫原子的标记我也检查过了没问题。 想问一下是我对收敛标准理解错了么？还有其它什么限制条件？ |
回复此楼» 猜你喜欢
谈谈两天一夜的“延安行”
已经有15人回复
-
EST投稿状态问题
已经有6人回复
-
职称评审没过，求安慰
已经有15人回复
-
垃圾破二本职称评审标准
已经有11人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有16人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
ggdh
专家顾问 (正式写手)
-
专家经验: +69 - 应助: 52 (初中生)
- 金币: 1747.4
- 散金: 10
- 红花: 18
- 帖子: 697
- 在线: 206.1小时
- 虫号: 96968
- 注册: 2005-11-08
- 性别: GG
- 专业: 化学生物学与生物有机化学
- 管辖: 量子化学
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
ljw4010: 金币+2, 谢谢指导解答！ 2015-01-26 21:47:00
感谢参与，应助指数 +1
ljw4010: 金币+2, 谢谢指导解答！ 2015-01-26 21:47:00
ggdh
专家顾问 (正式写手)
-
专家经验: +69 - 应助: 52 (初中生)
- 金币: 1747.4
- 散金: 10
- 红花: 18
- 帖子: 697
- 在线: 206.1小时
- 虫号: 96968
- 注册: 2005-11-08
- 性别: GG
- 专业: 化学生物学与生物有机化学
- 管辖: 量子化学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
yysskk: 金币+50, ★★★★★最佳答案 2015-01-24 21:53:36
ljw4010: 金币+2, 谢谢指导解答！ 2015-01-26 21:47:16
yysskk: 金币+50, ★★★★★最佳答案 2015-01-24 21:53:36
ljw4010: 金币+2, 谢谢指导解答！ 2015-01-26 21:47:16
zyqh
铁杆木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 6649.8
- 散金: 1454
- 红花: 22
- 帖子: 2365
- 在线: 593.4小时
- 虫号: 2762045
- 注册: 2013-10-29
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
锐利的碎片
木虫 (正式写手)
star watcher
- 1ST强帖: 1
- 应助: 136 (高中生)
- 金币: 3637.1
- 散金: 252
- 红花: 22
- 帖子: 988
- 在线: 1224.9小时
- 虫号: 961933
- 注册: 2010-03-05
- 专业: 凝聚态物性 II ：电子结构