应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » <求助>VASP 结构优化收敛，静态计算不收敛。（INCAR KPOINTS POSCAR附上）
51/1返回列表
查看: 4457  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

chord7806

新虫 (初入文坛)

[求助] <求助>VASP 结构优化收敛，静态计算不收敛。（INCAR KPOINTS POSCAR附上） 已有1人参与

最近计算了一个体系，共40原子。两种磁性离子。
先进行了结构优化。得到了CONTCAR，然后cp CONTCAR POSCAR静态计算能量始终不收敛。
希望各位老师，同学，达人，帮帮忙，看一下，到底问题出在哪？不胜感激。
下面是初始优化时的INCAR KPOINTS POSCAR
部分POSCAR：

system=******
    1.000000000000000
     5.163711736218603       0.000000000000000       0.030191807469367
     0.000000000000000       5.546700000000000       0.000000000000000
     0.000000000000000       0.000000000000000      14.830600000000000
    24  4  4  8
Direct
   0.3840999901294709   0.9585000276565552   0.1205499991774557
   0.1970999985933304   0.1957000046968459   0.4712499994784595
   0.3228000104427338   0.6952999830245971   0.4703499991446734
   0.1159000098705291   0.4585000276565551   0.1294500008225444
   0.3029000014066697   0.6957000046968459   0.2787500005215405
   0.1771999895572662   0.1952999830245970   0.2796500008553267
   0.6159000098705292   0.0414999723434448   0.3794500008225443
   0.8029000014066697   0.8042999953031542   0.0287500005215406
   0.6771999895572662   0.3047000169754029   0.0296500008553267
   0.8840999901294708   0.5414999723434448   0.3705499991774557
INCAR:

System = ******
PREC = Accurate
ISTART = 0
ICHARG = 2
ISPIN = 2
LREAL = Auto

ENCUT = 500
EDIFF = 1e-6
NELM = 40
GGA = PS

NSW = 300
IBRION = 2
ISIF = 3
EDIFFG = -0.001
MAGMOM = 24*0 -3 3 3 -3 -3 3 3 -3 8*0

LDAU = TRUE
LDAUU = 0.00 3.00 3.00 0.00
LDAUJ = 0.00 0.00 0.00 0.00
LDAUTYPE = 2
LDAUL = -1 2 2 -1
POTIM = 0.2
LMAXMIX = 4

ISMEAR = 0
SIGMA = 0.01
LORBIT = 11

LCHARG = FALSE
LWAVE = FALSE

KPOINTS：
Automatic
0
Monkhorst-Pack
6  6  2
0  0  0

计算完以后有这样的提示：
73 F= -.30352737E+03 E0= -.30352737E+03  d E =-.220886E-05  mag=     0.0000
curvature:  -0.06 expect dE=-0.192E-06 dE for cont linesearch -0.395E-08
ZBRENT: bisectioning
opt :   0.4334  next Energy=  -303.527366 (dE=-0.221E-05)
bond charge predicted
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.303527365963E+03    0.40205E-06   -0.12261E-07  2024   0.516E-03    0.155E-03
DAV:   2    -0.303527366101E+03   -0.13772E-06   -0.10919E-07  2024   0.500E-02    0.155E-03
DAV:   3    -0.303527366167E+03   -0.66240E-07   -0.41386E-08  2024   0.952E-03    0.155E-03
DAV:   4    -0.303527366279E+03   -0.11196E-06   -0.75204E-08  2024   0.228E-02    0.156E-03
DAV:   5    -0.303527366865E+03   -0.58592E-06   -0.15596E-06  2024   0.303E-02
  74 F= -.30352737E+03 E0= -.30352737E+03  d E =-.315061E-05  mag=     0.0000
curvature:  -0.02 expect dE=-0.165E-07 dE for cont linesearch -0.300E-08
trial: gam= 0.00000 g(F)=  0.285E-06 g(S)=  0.390E-06 ort = 0.815E-06 (trialstep = 0.450E+00)
search vector abs. value=  0.675E-06
reached required accuracy - stopping structural energy minimisation

OUTCAR里的精度：

Total       0.01803     0.01087     0.01049     0.00000     0.00000     0.01127
  in kB       0.06868     0.04140     0.03997     0.00000     0.00000     0.04294
  external pressure =        0.05 kB  Pullay stress =        0.00 kB
和

POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      1.95504      5.28764      1.81560         0.000493     -0.000529     -0.000062
      1.04003      1.02590      7.01562        -0.000502     -0.000445      0.000609
      1.68657      3.89042      6.97157        -0.000104      0.000333     -0.000392
      0.62891      2.52622      1.91875        -0.000493     -0.000529      0.000062
      1.57671      3.78732      4.13456         0.000502     -0.000445     -0.000609
      0.93018      1.12900      4.17862         0.000104      0.000333      0.000392
      3.19679      0.22461      5.64936        -0.000192      0.000558      0.000228
      4.13049      4.50077      0.45157         0.000433     -0.000113     -0.000710
      3.48222      1.63383      0.49540         0.000351      0.000060      0.000336
      4.55507      2.98603      5.55370         0.000192      0.000558     -0.000228
然后用优化的CONTCAR复制成POSCAR做静态计算，问题出现了，能量不收敛：
下面是提示内容：

DAV:  57    -0.295043281006E+03   -0.41838E-03   -0.36422E-01  7424   0.295E+01    0.338E+00
DAV:  58    -0.294995226835E+03    0.48054E-01   -0.55161E-01  6560   0.474E+01    0.334E+00
DAV:  59    -0.294998850811E+03   -0.36240E-02   -0.69107E-01  7696   0.140E+01    0.355E+00
DAV:  60    -0.295061875171E+03   -0.63024E-01   -0.12717E-01  6760   0.109E+01
   1 F= -.29506188E+03 E0= -.29506186E+03  d E =-.325714E-04  mag=     0.0000
而且发现OUTCAR中的 external pressure变的很大：external pressure =    -3550.84 kB，与此同时，磁矩也变的不准确了。

我的问题是：
1，有人遇到过这样的情况吗：优化结构后，算单点能量不收敛的。
2，我的INCAR中有什么是不对的吗，如果需要我贴出什么其他的信息，请在本帖留言。
真心希望问题得到解决，感激不尽。。。。。

[ Last edited by chord7806 on 2012-5-29 at 16:13 ]
回复此楼

» 收录本帖的淘帖专辑推荐

第一性

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-05-29 16:11:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wanghouyuan

铜虫 (小有名气)
请问你optimization和static cal在INCAR里的参数设置有何区别？
赞一下 回复此楼
子弹
2楼2012-05-29 16:50:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chord7806

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by wanghouyuan at 2012-05-29 16:50:33
请问你optimization和static cal在INCAR里的参数设置有何区别？

这个是我静态计算的INCAR，你看一下吧：
System = *******
PREC = Accurate
ISTART = 0
ICHARG = 2
ISPIN = 2
LREAL = Auto

ENCUT = 500
EDIFF = 1e-6
NELM = 60
GGA = PS

NSW = 0
IBRION = -1
ISIF = 3
EDIFFG = -0.001
MAGMOM = 24*0 -3 3 3 -3 -3 3 3 -3 8*0

LDAU = TRUE
LDAUU = 0.00 3.00 3.00 0.00
LDAUJ = 0.00 0.00 0.00 0.00
LDAUTYPE = 2
LDAUL = -1 2 2 -1
POTIM = 0.2
LMAXMIX = 4

ISMEAR = 0
SIGMA = 0.01
LORBIT = 11

LCHARG = FALSE
LWAVE = FALSE

我就把IBRION由原来的2改成-1了，NSW有300改成0了。。。

多谢关注。。
赞一下 回复此楼
3楼2012-05-29 17:07:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wanghouyuan

铜虫 (小有名气)

【答案】应助回帖

★
感谢参与，应助指数 +1
fzx2008: 金币+1, 谢谢回帖！ 2012-05-29 18:42:14
en ,我也是新手，也有这一类问题，即力达不到收敛标准。建议ISIF改为2（既然结构 已经优化完了），提高一下SIGMA，设置参数NELMDL。我也想验证一下，奈何机器紧张啊.
赞一下 回复此楼
子弹
4楼2012-05-29 17:29:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nami纳米nm

铁虫 (小有名气)

【答案】应助回帖

建议离子弛豫把EDIFFG的值设置大一点
赞一下 回复此楼
欲戴皇冠，必承其重
5楼2018-07-07 10:22:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chord7806 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 考研化学学硕调剂，一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 286求调剂 +3 lemonzzn 2026-03-16 5/250 2026-03-16 20:43 by lemonzzn
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 255求调剂 +3 李嘉慧， 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生，2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
信息提示
关闭
请填处理意见
关闭