| 5 | 1/1 | 返回列表 |
| 查看: 1279 | 回復(fù): 2 | ||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | ||
qinxuming新蟲(chóng) (小有名氣)
|
[求助]
vasp幾何優(yōu)化離子收斂的問(wèn)題 已有1人參與
|
|
|
我做的幾何優(yōu)化計(jì)算,我認(rèn)為離子步?jīng)]有達(dá)到收斂標(biāo)準(zhǔn)啊,輸出結(jié)果怎么說(shuō)精度已經(jīng)被滿(mǎn)足了呢?達(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 |
新蟲(chóng) (小有名氣)
金蟲(chóng) (小有名氣)
|
因?yàn)槟阍O(shè)置的ionic relaxation的criteria是energy啊,可以看到你最后給出的結(jié)果d Energy值是小于4E-5的。 還有一般ionic relaxation的criteria一般都設(shè)置為force condition吧。一般EDIFFG=-0.01 eV/A 就夠了。之后如果你不放心的話(huà)在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 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 346求調(diào)劑 一志愿070303有機(jī)化學(xué) +4 | 蘿卜燉青菜 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6 | 王一一依依 2026-03-30 | 8/400 |
|
|
[考研] 318求調(diào)劑 +7 | 陳晨79 2026-03-30 | 7/350 |
|
|
[考研] 一志愿北京工業(yè)大學(xué),324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +8 | 柳溪邊 2026-03-26 | 8/400 |
|
|
[考研] 299求調(diào)劑 +10 | 15188958825 2026-03-25 | 10/500 |
|
|
[考研] 349求調(diào)劑 +6 | 李木子啊哈哈 2026-03-25 | 6/300 |
|
|
[考研] 2026年華南師范大學(xué)歡迎化學(xué),化工,生物,生醫(yī)工等專(zhuān)業(yè)優(yōu)秀學(xué)子加入! +3 | llss0711 2026-03-28 | 6/300 |
|
|
[考研] 305求調(diào)劑 +8 | RuiFairyrui 2026-03-28 | 8/400 |
|
|
[考研] 數(shù)一英一271專(zhuān)碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 311(085601)求調(diào)劑 +4 | liziyeyeye 2026-03-28 | 4/200 |
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
[考研] 266求調(diào)劑 +11 | 陽(yáng)陽(yáng)哇塞 2026-03-27 | 12/600 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫(kù) 2026-03-26 | 5/250 |
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
[考研] 342求調(diào)劑 +3 | 加油a李zs 2026-03-26 | 3/150 |
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚(yú)愛(ài)有機(jī) 2026-03-25 | 4/200 |
|
|
[考研] 打過(guò)很多競(jìng)賽,085406控制工程300分,求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] 26考研-291分-廈門(mén)大學(xué)(085601)-柔性電子學(xué)院材料工程專(zhuān)業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|