| 3 | 1/1 | 返回列表 |
| 查看: 1271 | 回復(fù): 2 | ||
qinxuming新蟲 (小有名氣)
|
[求助]
vasp幾何優(yōu)化離子收斂的問題 已有1人參與
|
|
我做的幾何優(yōu)化計(jì)算,我認(rèn)為離子步?jīng)]有達(dá)到收斂標(biāo)準(zhǔn)啊,輸出結(jié)果怎么說精度已經(jīng)被滿足了呢?達(dá)到標(biāo)準(zhǔn)的含義是什么? 我的INCAR是 SYSTEM = crystle in a box ISMEAR=-5 ISPIN=1 PREC=Accurate ENCUT=700 ISYM=0 NELMIN = 8 # do a minimum of ten electronic steps NELM=100 EDIFF = 4E-6 # high accuracy for electronic groundstate NSW =150 # 150 ionic steps should do IBRION = 2 # ISIF=3 #all are optimized NPAR=4 那么離子的精度默認(rèn)應(yīng)該為EDIFF = 4E-5 但是我最后的OUTCAR是這樣: FREE ENERGIE OF THE ION-ELECTRON SYSTEM (eV) --------------------------------------------------- free energy TOTEN = -36.36260684 eV energy without entropy= -36.36260684 energy(sigma->0) = -36.36260684 d Force = 0.2275733E-03[-0.108E-04, 0.466E-03] d Energy = 0.3155840E-03-0.880E-04 d Force =-0.9954701E+00[-0.997E+00,-0.994E+00] d Ewald = 0.3477743E+01-0.447E+01 reached required accuracy - stopping structural energy minimisation OUTCAR前面是這樣: Electronic Relaxation 1 ENCUT = 700.0 eV 51.45 Ry 7.17 a.u. 16.85 16.85 43.15*2*pi/ulx,y,z ENINI = 700.0 initial cutoff ENAUG = 531.4 eV augmentation charge cutoff NELM = 100; NELMIN= 8; NELMDL= -5 # of ELM steps EDIFF = 0.4E-05 stopping-criterion for ELM LREAL = F real-space projection NLSPLINE = F spline interpolate recip. space projectors LCOMPAT= F compatible to vasp.4.4 GGA_COMPAT = T GGA compatible to vasp.4.4-vasp.4.6 LMAXPAW = -100 max onsite density LMAXMIX = 2 max onsite mixed and CHGCAR VOSKOWN= 0 Vosko Wilk Nusair interpolation ROPT = 0.00000 0.00000 Ionic relaxation EDIFFG = 0.4E-04 stopping-criterion for IOM NSW = 150 number of steps for IOM NBLOCK = 1; KBLOCK = 150 inner block; outer block IBRION = 2 ionic relax: 0-MD 1-quasi-New 2-CG NFREE = 1 steps in history (QN), initial steepest desc. (CG) ISIF = 3 stress and relaxation IWAVPR = 11 prediction: 0-non 1-charg 2-wave 3-comb ISYM = 0 0-nonsym 1-usesym 2-fastsym LCORR = T Harris-Foulkes like correction to forces |
金蟲 (小有名氣)
|
因?yàn)槟阍O(shè)置的ionic relaxation的criteria是energy啊,可以看到你最后給出的結(jié)果d Energy值是小于4E-5的。 還有一般ionic relaxation的criteria一般都設(shè)置為force condition吧。一般EDIFFG=-0.01 eV/A 就夠了。之后如果你不放心的話在OUTCAR搜索優(yōu)化最后一步的“TOTAL-FORCE” 然后可以看到每一個(gè)原子上的受力,每個(gè)力都小于你的EDIFFG就可以了。 例如:所有原子受力都小于0.01 eV/A POSITION TOTAL-FORCE (eV/Angst) ----------------------------------------------------------------------------------- -0.87110 7.28732 6.55372 0.003283 -0.000509 0.002108 -1.18676 3.72336 6.55369 -0.001249 0.001833 0.002096 2.05755 5.23204 6.55361 -0.000977 -0.002545 0.002264 |
新蟲 (小有名氣)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4 | 脫穎而出 2026-03-16 | 12/600 |
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|