| 22 | 1/1 | 返回列表 |
| 查看: 5991 | 回復(fù): 21 | ||||
[交流]
【求助】vasp靜態(tài)計(jì)算能量值錯(cuò)誤
|
||||
|
優(yōu)化了一個(gè)結(jié)構(gòu)后,進(jìn)行靜態(tài)計(jì)算,+U設(shè)置為零,發(fā)現(xiàn)電子步迭代時(shí)的能量值都非常大,如下: free energy TOTEN = -386068.93443124 eV free energy TOTEN = -468284.83293180 eV free energy TOTEN = -483184.31035624 eV free energy TOTEN = -488841.30556393 eV free energy TOTEN = -491852.20379074 eV free energy TOTEN = -278816.17117876 eV free energy TOTEN = -10667.33547639 eV free energy TOTEN = 36658.42027096 eV free energy TOTEN = 3020.73440924 eV free energy TOTEN = -86463.86284949 eV free energy TOTEN = -120769.97367841 eV free energy TOTEN = -142957.25416115 eV free energy TOTEN = -175335.77020944 eV free energy TOTEN = -171597.62982132 eV free energy TOTEN = -165579.41016696 eV free energy TOTEN = -170144.20096520 eV free energy TOTEN = -144730.46667812 eV free energy TOTEN = -200129.89254844 eV free energy TOTEN = -136862.06349103 eV 優(yōu)化的結(jié)構(gòu)經(jīng)過(guò)檢查并作為初始結(jié)構(gòu)再次優(yōu)化,得到仍是這個(gè)結(jié)構(gòu),應(yīng)該不會(huì)有問(wèn)題。奇怪的是加上一定的U值后計(jì)算沒有出現(xiàn)異常,而且再回頭重新進(jìn)行+U為零的計(jì)算又能得到正常的總能了。 請(qǐng)問(wèn)可能造成這種情況的原因有哪些? 今天再檢查發(fā)現(xiàn)從以上電子步以后迭代不再進(jìn)行下去了,出現(xiàn)很多警告信息: WARNING in EDDRMM: call to ZHEGV failed 求解! 發(fā)現(xiàn)以優(yōu)化結(jié)果的CHGCAR作為初始電荷密度進(jìn)行靜態(tài)計(jì)算后,不會(huì)出現(xiàn)上述問(wèn)題。 [ Last edited by nkleof on 2010-11-16 at 09:59 ] |
量子化學(xué)基礎(chǔ) |
» 搶金幣啦!回帖就可以得到:
+3/1802
+5/1680
+1/461
+5/175
+2/146
+1/81
+1/79
+1/76
+2/64
+1/63
+1/38
+2/34
+1/33
+1/16
+1/13
+1/12
+2/8
+2/8
+1/5
+1/4
鐵桿木蟲 (正式寫手)
|
在靜態(tài)計(jì)算時(shí)+U,優(yōu)化時(shí)沒有+U。 具體設(shè)置如下: PREC = High LREAL = A ALGO = Fast NPAR = 4 NSIM = 8 ISTART = 0 ICHARG = 0 LMAXMIX = 4 NELMIN = 4 LORBIT = 12 EDIFF = 1E-5 EMAX = 20 EMIN = -20 NEDOS = 2000 ISMEAR = -5 ISPIN = 2 MAGMOM = 16*0 34*4 72*0 LDAU = .TRUE. LDAUTYPE = 2 LDAUL = 2 2 -1 LDAUU = 0 0 0 LDAUJ = 0 0 0 LDAUPRINT = 2 U值在0~2變化。 |
|
本帖內(nèi)容被屏蔽 |
|
本帖內(nèi)容被屏蔽 |
|
一般要查看能量我都是找OUTCAR的,OSZICAR的數(shù)值和OUTCAR應(yīng)該是一樣,可能有小差別,我覺得還是應(yīng)該以O(shè)UCAR為準(zhǔn)。 OUTCAR里面會(huì)記錄每一個(gè)電子步優(yōu)化的能量結(jié)果,就是我前面列出的free energy TOTEN,達(dá)到收斂后會(huì)有另外一個(gè)free energy TOTEN,以及energy without entropy,是最終的基態(tài)能量,在OUTCAR文件的最后面部分。 |
榮譽(yù)版主 (文壇精英)
木蟲 (職業(yè)作家)
| 22 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +6 | 十三加油 2026-03-21 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考研] 求調(diào)劑 +7 | Auroracx 2026-03-22 | 7/350 |
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 4/200 |
|
|
[考研] 廣西大學(xué)材料導(dǎo)師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
[考研] 299求調(diào)劑 +5 | shxchem 2026-03-20 | 7/350 |
|
|
[基金申請(qǐng)]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
[考研] 一志愿中國(guó)石油大學(xué)(華東) 本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-17 | 3/150 |
|
|
[考研] 一志愿華南師大 070300(化學(xué))304分求調(diào)劑 +3 | 0703武芊慧雪304 2026-03-18 | 3/150 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 6/300 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生(接收計(jì)算機(jī)專業(yè)調(diào)劑) +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專利,省、國(guó)級(jí)證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|