| 3 | 1/1 | 返回列表 |
| 查看: 1272 | 回復(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” 然后可以看到每一個原子上的受力,每個力都小于你的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ā)表 | |
|---|---|---|---|---|
|
[考研]
|
簡木ChuFront 2026-03-19 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 一志愿中南化學(xué)337求調(diào)劑 +4 | niko- 2026-03-19 | 5/250 |
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +4 | 西南交通專材355 2026-03-19 | 4/200 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 生物學(xué)調(diào)劑招人。! +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 一志愿天大材料與化工(085600)總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑,求各位老師收留 +10 | 秋有木北 2026-03-14 | 10/500 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[基金申請] 被我言中:新模板不強(qiáng)調(diào)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|