| 7 | 1/1 | 返回列表 |
| 查看: 1448 | 回復(fù): 6 | |||||
[求助]
明明達到了收斂標準,卻還在計算 已有1人參與
|
|
今天在弛豫一個體系的時候,無意中發(fā)現(xiàn)明明早已達到收斂標準了,計算卻還在繼續(xù)。 這是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步的結(jié)果: 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的原子是非弛豫原子,其余所有原子早已達到收斂標準,為了防止我眼花看錯了,我還特意寫程序幫我查找,都沒有發(fā)現(xiàn)弛豫原子上的力超過0.05ev的。 OUTCAR文件開頭的參數(shù)里也寫明了 EDIFFG = -.5E-01 stopping-criterion for IOM POSCAR里面弛豫和非弛豫原子的標記我也檢查過了沒問題。 想問一下是我對收斂標準理解錯了么?還有其它什么限制條件? |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +69 |
|
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的 |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +69 |
我沒有仔細審題吧所有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 我想這應(yīng)該是判斷標準 |

木蟲 (正式寫手)
star watcher
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 281求調(diào)劑(0805) +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機用戶 2026-03-17 | 9/450 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 一志愿,福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
[考研] 機械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|