| 查看: 1718 | 回复: 8 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
lingch1990木虫 (小有名气)
|
[求助]
大家来帮我看看，实在是没折了，能想到的都试了。。。 已有3人参与
|
||
|
对一个层状结构的原胞进行优化，加了范德华修正之后（LVDW=.TRUE.），计算结果就出问题了，体系能量为正值，仔细检查了各个环节，没发现原因，求助大神！ 这是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 最后居然也能收敛，但是能量却始终是正值，该怎么办？！ |
回复此楼» 猜你喜欢
313求调剂
已经有4人回复
-
311求调剂一志愿合肥工业大学
已经有6人回复
-
0703化学调剂，求导师收
已经有11人回复
-
284求调剂
已经有5人回复
-
求调剂
已经有4人回复
-
求调剂
已经有4人回复
-
总分322求生物学/生化与分子/生物信息学相关调剂
已经有7人回复
-
070300化学求调剂
已经有6人回复
-
070300一志愿211，312分求调剂院校
已经有7人回复
-
085600 295分求调剂
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 正大天晴达沙替尼、伊马替尼火箭神速”之我见 转载
已经有43人回复
- 【转载】西游记科研版
已经有17人回复
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
lingch1990: 金币+5, ★有帮助 2014-07-01 11:08:00
感谢参与，应助指数 +1
lingch1990: 金币+5, ★有帮助 2014-07-01 11:08:00
|
考虑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 |
3楼2014-07-01 10:13:14
lingch1990
木虫 (小有名气)
- 应助: 15 (小学生)
- 金币: 4313.7
- 散金: 258
- 红花: 4
- 帖子: 209
- 在线: 439.9小时
- 虫号: 1900956
- 注册: 2012-07-19
- 性别: GG
- 专业: 理论和计算化学
4楼2014-07-01 11:08:41
【答案】应助回帖
感谢参与，应助指数 +1
|
http://cms.mpi.univie.ac.at/vasp/vasp/DFT_D2_method_Grimme.html 我也感觉你好像是少设置了不少参数，你看看上面的这个官方说明 |
5楼2014-07-01 12:29:20
lingch1990
木虫 (小有名气)
- 应助: 15 (小学生)
- 金币: 4313.7
- 散金: 258
- 红花: 4
- 帖子: 209
- 在线: 439.9小时
- 虫号: 1900956
- 注册: 2012-07-19
- 性别: GG
- 专业: 理论和计算化学
6楼2014-07-01 12:50:23
