| 22 | 1/1 | 返回列表 |
| 查看: 6007 | 回復(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)過檢查并作為初始結(jié)構(gòu)再次優(yōu)化,得到仍是這個(gè)結(jié)構(gòu),應(yīng)該不會(huì)有問題。奇怪的是加上一定的U值后計(jì)算沒有出現(xiàn)異常,而且再回頭重新進(jìn)行+U為零的計(jì)算又能得到正常的總能了。 請問可能造成這種情況的原因有哪些? 今天再檢查發(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)上述問題。 [ Last edited by nkleof on 2010-11-16 at 09:59 ] |
量子化學(xué)基礎(chǔ) |
» 搶金幣啦!回帖就可以得到:
+2/394
+2/132
+1/81
+1/79
+1/43
+2/32
+1/18
+1/15
+1/14
+1/12
+1/12
+1/7
+1/7
+1/5
+1/4
+1/3
+1/3
+1/2
+1/2
+1/1
鐵桿木蟲 (正式寫手)
|
在靜態(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ā)表 | |
|---|---|---|---|---|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap- 2026-03-22 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 08工學(xué)調(diào)劑 +6 | 用戶573181 2026-03-20 | 10/500 |
|
|
[考研] 284求調(diào)劑 +6 | Zhao anqi 2026-03-22 | 6/300 |
|
|
[考研] 289材料與化工(085600)B區(qū)求調(diào)劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|
[考研] 306求調(diào)劑 +5 | 來好運(yùn)來來來 2026-03-22 | 5/250 |
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
[考研] 考研調(diào)劑 +4 | 來好運(yùn)來來來 2026-03-21 | 4/200 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 一志愿中國海洋大學(xué),生物學(xué),301分,求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|