| 21 | 1/1 | 返回列表 |
| 查看: 4410 | 回復(fù): 20 | |||
stractor金蟲(chóng) (著名寫手)
|
[交流]
【求助】如何看力
|
|
在VASP中,優(yōu)化一般以力為收斂標(biāo)準(zhǔn),怎么看力呢? 我在OUTCAR中發(fā)現(xiàn)了很多關(guān)于力的,比如 FORCE on cell =-STRESS in cart. coord. units (eV/reduce length): Direction X Y Z XY YZ ZX -------------------------------------------------------------------------------------- Alpha Z 228.09 228.09 228.09 。。。。。。 Total -3.54 -4.54 13.70 -4.58 -0.30 0.42 in kB -5.41 -6.94 20.92 -6.99 -0.46 0.64 external pressure = 2.86 kB Pullay stress = 0.00 kB 。。。。。。 FORCES acting on ions electron-ion (+dipol) ewald-force non-local-force convergence-correction ----------------------------------------------------------------------------------------------- -.397E+00 0.126E-01 -.195E+00 0.345E+00 0.378E-02 0.191E+00 0.492E-01 0.508E-02 0.287E-02 0.149E-02 0.121E-02 -.427E-03 -.279E+00 0.448E+00 0.232E+00 0.268E+00 -.385E+00 -.192E+00 0.151E-01 -.731E-01 -.463E-01 0.142E-02 -.120E-02 -.535E-03 。。。。。。 0.720E+00 0.667E+00 0.351E+01 -.124E+00 -.874E+00 0.151E+02 0.326E-01 0.278E-01 -.409E-01 0.250E-03 0.143E-03 -.281E-03 ----------------------------------------------------------------------------------------------- 0.963E+00 -.298E+00 0.194E+02 0.185E-12 -.108E-13 -.799E-13 -.970E+00 0.301E+00 -.194E+02 0.211E-02 0.110E-02 -.395E-02 POSITION TOTAL-FORCE (eV/Angst) ----------------------------------------------------------------------------------- 10.06599 10.06719 0.01673 -0.001720 0.022612 -0.001236 10.06703 5.03749 10.34262 0.005402 -0.011136 -0.007025 。。。。。。 |
» 搶金幣啦!回帖就可以得到:
+1/474
+1/69
+2/62
+2/58
+1/38
+1/21
+1/20
+1/12
+1/12
+1/8
+1/8
+1/7
+1/4
+1/4
+1/3
+1/3
+1/3
+1/2
+1/1
+1/1
鐵桿木蟲(chóng) (正式寫手)
榮譽(yù)版主 (知名作家)
榮譽(yù)版主 (知名作家)
金蟲(chóng) (著名寫手)
|
EDIFFG=-0.03。 個(gè)人感覺(jué)是 POSITION TOTAL-FORCE (eV/Angst) ----------------------------------------------------------------------------------- 10.06491 10.06976 0.01770 0.010282 0.007834 0.002618 ...................... 不知道是不是? [ Last edited by stractor on 2010-11-25 at 11:24 ] |
鐵桿木蟲(chóng) (正式寫手)
|
不會(huì)吧,我的所有OUTCAR里都有這一項(xiàng)啊 3.51061 2.02685 5.67094 0.000000 0.000000 0.018799 5.85101 2.02685 0.00000 0.000000 0.000000 0.715101 4.68081 4.05370 1.65492 0.000000 0.000000 -0.774121 3.51061 2.02685 3.30983 0.000000 0.000000 -0.118627 5.85101 2.02685 4.90216 0.000000 0.000000 0.016457 ----------------------------------------------------------------------------------- total drift: 0.000000 -0.000001 0.012219 -------------------------------------------------------------------------------------------------------- FORCES: max atom, RMS 0.026838 0.012964 FORCE total and by dimension 0.073335 0.026838 OPT: skip step - force has converged |
榮譽(yù)版主 (知名作家)
鐵桿木蟲(chóng) (正式寫手)
榮譽(yù)版主 (知名作家)
木蟲(chóng) (正式寫手)
鐵桿木蟲(chóng) (正式寫手)
鐵桿木蟲(chóng) (正式寫手)
木蟲(chóng) (正式寫手)
|
呵呵 這個(gè)問(wèn)題我原來(lái)也是不懂 用了半年vasp也沒(méi)關(guān)心著問(wèn)題 主要功夫是看能量去了 嘿嘿 我給個(gè)正解吧(自認(rèn)為正確的 哈哈) 在OUTCAR里是能找到如下信息的(4.6版) POSITION TOTAL-FORCE (eV/Angst) ----------------------------------------------------------------------------------- 2.32720 9.15829 7.18156 0.002637 -0.004234 0.000292 8.62320 9.15829 7.18156 0.002637 -0.004234 0.000292 0.75579 11.92243 7.16200 -0.001375 0.001914 -0.000926 7.05179 11.92243 7.16200 -0.001375 0.001914 -0.000926 0.75649 17.47136 7.12124 -0.000949 0.002844 0.000262 7.05249 17.47136 7.12124 -0.000949 0.002844 0.000262 5.47484 9.15824 7.18139 0.001707 -0.003726 -0.000903 每排后3個(gè)數(shù)據(jù)就是力,一般要所有的力數(shù)據(jù)收斂到你INCAR設(shè)定的標(biāo)準(zhǔn) relax就會(huì)結(jié)束了 |
金蟲(chóng) (著名寫手)
鐵桿木蟲(chóng) (正式寫手)
|
是的。比如,我在我的一個(gè)優(yōu)化好的目錄下輸入vef.pl 后,會(huì)出現(xiàn), 0 1.59483500 -216.043454 0 1 1.50700800 -216.063328 -0.019874 2 1.39912000 -216.202084 -0.15863 3 1.35256700 -216.214994 -0.17154 4 0.97341100 -216.346551 -0.303097 5 0.90999800 -216.357609 -0.314155 6 0.74630000 -216.456104 -0.41265 7 0.70556300 -216.463665 -0.420211 8 0.60424300 -216.549739 -0.506285 9 0.56631000 -216.556446 -0.512992 10 0.34305100 -216.586833 -0.543379 11 0.31455300 -216.591458 -0.548004 12 0.55732000 -216.624470 -0.581016 13 0.52821900 -216.628386 -0.584932 14 0.28016700 -216.658036 -0.614582 15 0.25022700 -216.661423 -0.617969 16 0.27186400 -216.670846 -0.627392 17 0.22158600 -216.673514 -0.63006 18 0.15353700 -216.679686 -0.636232 19 0.12014000 -216.681291 -0.637837 20 0.20970900 -216.685369 -0.641915 21 0.13338200 -216.686432 -0.642978 22 0.08766300 -216.687047 -0.643593 23 0.06567300 -216.687499 -0.644045 24 0.07568400 -216.687493 -0.644039 25 0.04350000 -216.688126 -0.644672 26 0.07302700 -216.688421 -0.644967 27 0.04070300 -216.688689 -0.645235 28 0.06882400 -216.688672 -0.645218 29 0.05552600 -216.689060 -0.645606 30 0.05419000 -216.689335 -0.645881 31 0.04710200 -216.689607 -0.646153 32 0.03033700 -216.689594 -0.64614 33 0.03087700 -216.689797 -0.646343 34 0.03965800 -216.689809 -0.646355 35 0.04456100 -216.689733 -0.646279 36 0.02284800 -216.689808 -0.646354 [ Last edited by fanchen021 on 2010-11-26 at 14:21 ] |
鐵桿木蟲(chóng) (正式寫手)
|
其中第一列是優(yōu)化步數(shù),第二列是最大原子受力,小于收斂標(biāo)準(zhǔn)優(yōu)化就結(jié)束了,第三列是總能,第四列是當(dāng)前步總能和第一步總能的差值,可以看總能的變化趨勢(shì)。 總的來(lái)說(shuō),UT的這個(gè)腳本看力還是很方便的 |
金蟲(chóng) (小有名氣)
|
本帖內(nèi)容被屏蔽 |
| 21 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 268求調(diào)劑 +5 | 好運(yùn)連綿不絕 2026-03-12 | 5/250 |
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +5 | Liwangman 2026-03-15 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 +6 | 妮妮ninicgb 2026-03-15 | 9/450 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
[考研] 308求調(diào)劑 +5 | 是Lupa啊 2026-03-11 | 5/250 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
[考研] 307求調(diào)劑 +5 | 超級(jí)伊昂大王 2026-03-12 | 5/250 |
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
[論文投稿]
投稿問(wèn)題
5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|
|
[考博] 26申博求助 +3 | 跳躍餅干 2026-03-10 | 4/200 |
|