| 查看: 2674 | 回复: 4 | |||
[求助]
vasp 离子步只跑了3步，dE还很大，就说收敛了，是怎么回事？？
|
|
我的结构是表面上有个cluster，120个原子。我在前面已经优化过了结构，收敛了，结构比较合理。 这次计算只是去掉了一个氧原子，想算一下空位形成能。 vasp 离子步只跑了3步，d E =-.576572E+02 还很大，就说收敛了，是怎么回事？？ EDIFF、NELM、EDIFFG都是默认值，POTIM=0.1. 我看电子步跑了60步到最大值了，要不要把NELM加大呢？ 下面简单列了OSZICAR的结果： entering main loop N E dE d eps ncg rms rms(c) RMM: 1 0.856523697897E+04 0.85652E+04 -0.23656E+05 1032 0.985E+02 RMM: 2 0.303777438339E+04 -0.55275E+04 -0.63732E+04 1032 0.282E+02 RMM: 3 0.120743650014E+04 -0.18303E+04 -0.20187E+04 1032 0.199E+02 RMM: 4 0.232904630700E+03 -0.97453E+03 -0.80517E+03 1032 0.189E+02 RMM: 5 -0.301928474056E+03 -0.53483E+03 -0.39350E+03 1032 0.116E+02 ........... RMM: 58 -0.613962913659E+03 -0.87356E+01 -0.51924E+00 2073 0.441E+00 0.439E+01 RMM: 59 -0.615169002585E+03 -0.12061E+01 -0.87625E-01 2080 0.201E+00 0.433E+01 RMM: 60 -0.615934407958E+03 -0.76541E+00 -0.26866E-01 2103 0.119E+00 1 F= -.61593441E+03 E0= -.61595251E+03 d E =-.615934E+03 mag= 96.9828 curvature: 0.00 expect dE= 0.000E+00 dE for cont linesearch 0.000E+00 trial: gam= 0.00000 g(F)= 0.382E+00 g(S)= 0.000E+00 ort = 0.000E+00 (trialstep = 0.100E+01) search vector abs. value= 0.382E+00 bond charge predicted N E dE d eps ncg rms rms(c) RMM: 1 -0.616414208634E+03 -0.12452E+01 -0.20533E+01 2115 0.606E+00 0.436E+01 RMM: 2 -0.855334771199E+03 -0.23892E+03 -0.18604E+02 2106 0.223E+01 0.137E+02 RMM: 3 -0.632885173476E+03 0.22245E+03 -0.10007E+02 2064 0.215E+01 0.770E+01 RMM: 4 -0.613976729445E+03 0.18908E+02 -0.34720E+01 2069 0.120E+01 0.354E+01 RMM: 5 -0.619880022167E+03 -0.59033E+01 -0.11977E+01 2074 0.761E+00 0.289E+01 ....... RMM: 57 -0.673116148001E+03 -0.16096E-02 -0.27028E-04 1340 0.304E-02 0.145E+00 RMM: 58 -0.673114612545E+03 0.15355E-02 -0.94439E-05 1220 0.179E-02 0.148E+00 RMM: 59 -0.673112177097E+03 0.24354E-02 -0.69843E-05 1211 0.168E-02 0.150E+00 RMM: 60 -0.673112422360E+03 -0.24526E-03 -0.66699E-05 1173 0.139E-02 2 F= -.67311242E+03 E0= -.67309912E+03 d E =-.571780E+02 mag= 82.5543 trial-energy change: -57.178014 1 .order -0.248124 -0.381829 -0.114420 step: 1.0003(harm= 1.4279) dis= 0.01946 next Energy= -673.112442 (dE=-0.572E+02) bond charge predicted N E dE d eps ncg rms rms(c) RMM: 1 -0.673111494709E+03 0.68239E-03 -0.28328E-05 1092 0.102E-02 0.151E+00 RMM: 2 -0.673243049842E+03 -0.13156E+00 -0.60539E-02 2120 0.487E-01 0.216E+00 RMM: 3 -0.673276965381E+03 -0.33916E-01 -0.37351E-02 2074 0.372E-01 0.240E+00 RMM: 4 -0.673283230064E+03 -0.62647E-02 -0.19439E-02 2070 0.291E-01 0.147E+00 RMM: 5 -0.673283122864E+03 0.10720E-03 -0.50074E-03 2043 0.132E-01 0.169E+00 ......... RMM: 56 -0.673574115277E+03 0.15949E-02 -0.14751E-03 1828 0.683E-02 0.145E+00 RMM: 57 -0.673567578727E+03 0.65366E-02 -0.59063E-04 1601 0.499E-02 0.147E+00 RMM: 58 -0.673565550345E+03 0.20284E-02 -0.23339E-04 1335 0.314E-02 0.147E+00 RMM: 59 -0.673564938652E+03 0.61169E-03 -0.18123E-04 1280 0.248E-02 0.145E+00 RMM: 60 -0.673591572160E+03 -0.26634E-01 -0.33278E-03 1974 0.118E-01 3 F= -.67359157E+03 E0= -.67356894E+03 d E =-.576572E+02 mag= 84.2310 curvature: -0.00 expect dE=-0.209E-03 dE for cont linesearch -0.389E-04 trial: gam= 0.18927 g(F)= 0.188E+00 g(S)= 0.000E+00 ort = 0.115E+00 (trialstep = 0.100E+01) search vector abs. value= 0.245E+00 reached required accuracy - stopping structural energy minimisation |
回复此楼» 猜你喜欢
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
zuocuiping
木虫 (职业作家)
- 应助: 14 (小学生)
- 金币: 4975.9
- 散金: 966
- 红花: 60
- 帖子: 4045
- 在线: 1600.3小时
- 虫号: 890422
- 注册: 2009-11-01
- 性别: GG
- 专业: 半导体光电子器件
2楼2019-05-15 17:08:44
|
IALGO = 48 ENCUT = 400.00 SYSTEM = MgCuFe2O4-Au PREC = accurate ISPIN = 2 NPAR = 2 ISTART=0 ICHARG=2 #DFT+U for Mg Cu Fe O Au LDAU=T LDAUL=-1 -1 2 -1 -1 LDAUU= 0 0 4 0 0 LDAUJ= 0 0 0 0 0 Ionic Relaxation NSW = 250 number of steps for IOM ISIF = 2 IBRION = 2 conjugate-gradient algorithm POTIM = 0.1 #EDIFFG = -0.05 DOS related values: ISMEAR = 1 SIGMA = 0.2 SGI Origin optimized values: LWAVE = .FALSE. LCHARG = .FALSE. LVTOT = .FALSE. LREAL=.FALSE. |
3楼2019-05-15 18:11:50
4楼2019-05-15 18:51:04
zuocuiping
木虫 (职业作家)
- 应助: 14 (小学生)
- 金币: 4975.9
- 散金: 966
- 红花: 60
- 帖子: 4045
- 在线: 1600.3小时
- 虫号: 890422
- 注册: 2009-11-01
- 性别: GG
- 专业: 半导体光电子器件
5楼2019-05-15 21:08:54
投票:
