| 9 | 1/1 | 返回列表 |
| 查看: 1700 | 回復: 8 | ||
lingch1990木蟲 (小有名氣)
|
[求助]
大家來幫我看看,實在是沒折了,能想到的都試了。。。 已有3人參與
|
|
對一個層狀結(jié)構(gòu)的原胞進行優(yōu)化,加了范德華修正之后(LVDW=.TRUE.),計算結(jié)果就出問題了,體系能量為正值,仔細檢查了各個環(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的計算,除了LVDW = .TRUE.,還要加入以下標簽 VDW_RADIUS = 30 VDW_SCALING = 0.75 VDW_D = 20.0 VDW_C6 = 可采用默認值 VDW_R0 = 可采用默認值 可參考以下鏈接 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ù),你看看上面的這個官方說明 |

木蟲 (小有名氣)

木蟲 (小有名氣)
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 305求調(diào)劑 +3 | RuiFairyrui 2026-03-28 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 材料與化工272求調(diào)劑 +9 | 阿斯蒂芬2004 2026-03-28 | 9/450 |
|
|
[考研] 070305高分子化學與物理 304分求調(diào)劑 +4 | c297914 2026-03-28 | 4/200 |
|
|
[考研] 一志愿北化085600材料專碩275|有文章專利|求調(diào)劑 +5 | Micky11223 2026-03-25 | 5/250 |
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +6 | 芋泥小鈴鐺 2026-03-28 | 6/300 |
|
|
[考研] 311求調(diào)劑 +3 | 希望上岸阿小楊 2026-03-23 | 3/150 |
|
|
[考研] 307求調(diào)劑 +8 | 超級伊昂大王 2026-03-24 | 9/450 |
|
|
[考研] 287求調(diào)劑 +10 | land xuxu 2026-03-26 | 10/500 |
|
|
[考研] 考研化學308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學大學 ,080500材料科學與工程學碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 尋找調(diào)劑 +5 | 倔強芒? 2026-03-21 | 8/400 |
|
|
[考研] 282求調(diào)劑 +3 | wcq131415 2026-03-24 | 3/150 |
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|