| 9 | 1/1 | 返回列表 |
| 查看: 1699 | 回復(fù): 8 | ||
lingch1990木蟲 (小有名氣)
|
[求助]
大家來幫我看看,實(shí)在是沒折了,能想到的都試了。。。 已有3人參與
|
|
對一個(gè)層狀結(jié)構(gòu)的原胞進(jìn)行優(yōu)化,加了范德華修正之后(LVDW=.TRUE.),計(jì)算結(jié)果就出問題了,體系能量為正值,仔細(xì)檢查了各個(gè)環(huán)節(jié),沒發(fā)現(xiàn)原因,求助大神! 這是INCAR SYSTEM=aa ISTART = 0 ICHARG = 2 ENCUT = 350 EDIFF = 1E-5 IBRION = 2 #IALGO = 48 #POTIM = 0.25 NSW = 500 LREAL= .FALSE. EDIFFG = -2E-2 ISMEAR = 0 SIGMA = 0.1 PREC = Accurate ISIF = 2 NPAR = 4 LVDW = .TRUE. LCHARG = .FALSE. LWAVE = .FALSE. 這是out文件 N E dE d eps ncg rms rms(c) RMM: 1 0.120694905873E+03 0.12069E+03 -0.36576E+03 2460 0.691E+02 RMM: 2 0.559986339940E+02 -0.64696E+02 -0.61631E+02 2460 0.154E+02 RMM: 3 0.439712166147E+02 -0.12027E+02 -0.11463E+02 2460 0.597E+01 RMM: 4 0.390931902337E+02 -0.48780E+01 -0.40976E+01 2460 0.381E+01 RMM: 5 0.369737547884E+02 -0.21194E+01 -0.18232E+01 2460 0.244E+01 . . . RMM: 56 0.760382798675E+02 -0.14809E-02 -0.50332E-03 5177 0.450E-01 0.165E-01 RMM: 57 0.760368953504E+02 -0.13845E-02 -0.49995E-03 5172 0.441E-01 0.169E-01 RMM: 58 0.760354517223E+02 -0.14436E-02 -0.48087E-03 5197 0.440E-01 0.171E-01 RMM: 59 0.760341433073E+02 -0.13084E-02 -0.47381E-03 5172 0.430E-01 0.174E-01 RMM: 60 0.760328257954E+02 -0.13175E-02 -0.45167E-03 5175 0.428E-01 1 F= 0.75145626E+02 E0= 0.75151151E+02 d E =0.751456E+02 curvature: 0.00 expect dE= 0.000E+00 dE for cont linesearch 0.000E+00 trial: gam= 0.00000 g(F)= 0.348E+03 g(S)= 0.000E+00 ort = 0.000E+00 (trialstep = 0.100E+01) search vector abs. value= 0.348E+03 bond charge predicted N E dE d eps ncg rms rms(c) RMM: 1 -0.605881548289E+02 -0.13662E+03 -0.21796E+02 6869 0.106E+02 0.283E+01 RMM: 2 0.109067495123E+02 0.71495E+02 -0.15808E+01 5334 0.306E+01 0.188E+01 . . . RMM: 56 0.546842852482E+02 -0.23498E-04 -0.46851E-05 4265 0.428E-02 0.109E-01 RMM: 56 0.546842852482E+02 -0.23498E-04 -0.46851E-05 4265 0.428E-02 0.109E-01 RMM: 57 0.546842815707E+02 -0.36775E-05 -0.41106E-05 4209 0.398E-02 2 F= 0.53810824E+02 E0= 0.53815411E+02 d E =-.213348E+02 trial-energy change: -21.334802 1 .order -99.407368 -347.647760 148.833023 step: 0.4655(harm= 0.7002) dis= 1.34786 next Energy= 2.109306 (dE=-0.730E+02) . . . bond charge predicted N E dE d eps ncg rms rms(c) RMM: 1 0.528533676651E+02 -0.91010E-03 -0.67153E-04 4921 0.250E-01 0.442E-02 RMM: 2 0.528535264401E+02 0.15877E-03 -0.67866E-05 4758 0.714E-02 0.298E-02 RMM: 3 0.528537676079E+02 0.24117E-03 -0.36750E-05 4007 0.507E-02 0.968E-03 RMM: 4 0.528537805196E+02 0.12912E-04 -0.17478E-05 2934 0.311E-02 0.204E-03 RMM: 5 0.528537780874E+02 -0.24322E-05 -0.12223E-05 2849 0.252E-02 7 F= 0.51998445E+02 E0= 0.52002399E+02 d E =-.231472E+02 curvature: -0.06 expect dE=-0.131E-05 dE for cont linesearch -0.151E-06 trial: gam= 0.00008 g(F)= 0.212E-04 g(S)= 0.000E+00 ort =-0.292E-01 (trialstep = 0.826E+00) search vector abs. value= 0.188E-04 reached required accuracy - stopping structural energy minimisation 最后居然也能收斂,但是能量卻始終是正值,該怎么辦?! |
木蟲 (小有名氣)
|
考慮DFT-D2 VDW的計(jì)算,除了LVDW = .TRUE.,還要加入以下標(biāo)簽 VDW_RADIUS = 30 VDW_SCALING = 0.75 VDW_D = 20.0 VDW_C6 = 可采用默認(rèn)值 VDW_R0 = 可采用默認(rèn)值 可參考以下鏈接 http://cms.mpi.univie.ac.at/vasp/vasp/DFT_D2_method_Grimme.html |

木蟲 (小有名氣)
|
http://cms.mpi.univie.ac.at/vasp/vasp/DFT_D2_method_Grimme.html 我也感覺你好像是少設(shè)置了不少參數(shù),你看看上面的這個(gè)官方說明 |

木蟲 (小有名氣)

木蟲 (小有名氣)
鐵桿木蟲 (著名寫手)
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工272求調(diào)劑 +7 | 阿斯蒂芬2004 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +6 | 腳滑的守法公民 2026-03-27 | 6/300 |
|
|
[考研] 277跪求調(diào)劑 +5 | 1915668 2026-03-27 | 9/450 |
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
[考研] 求調(diào)劑 +8 | 張zz111 2026-03-27 | 9/450 |
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[考研]
|
鐘llll 2026-03-26 | 4/200 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研]
材料調(diào)劑
5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +9 | Mr. Z 2026-03-25 | 9/450 |
|
|
[考研] 生物技術(shù)與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
[考研] 289材料與化工(085600)B區(qū)求調(diào)劑 +4 | 這么名字咋樣 2026-03-22 | 5/250 |
|
|
[考研] 化工專碩求調(diào)劑 +3 | question挽風(fēng) 2026-03-24 | 3/150 |
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 6/300 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問津 2026-03-23 | 5/250 |
|
|
[考研]
|
2117205181 2026-03-21 | 8/400 |
|