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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » <求助>VASP 结构优化收敛，静态计算不收敛。（INCAR KPOINTS POSCAR附上）
51/1返回列表
查看: 4466  |  回复: 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

铜虫 (小有名气)

【答案】应助回帖

★
感谢参与，应助指数 +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的回帖
查看全部 5 个回答

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的回帖

nami纳米nm

铁虫 (小有名气)

【答案】应助回帖

建议离子弛豫把EDIFFG的值设置大一点
赞一下 回复此楼
欲戴皇冠，必承其重
5楼2018-07-07 10:22:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 265求调剂 +3 Jack?k?y 2026-03-17 3/150 2026-03-21 03:17 by JourneyLucky
[考研] 278求调剂 +6 烟火先于春 2026-03-17 6/300 2026-03-21 01:57 by JourneyLucky
[考研] 279分求调剂 一志愿211 +11 chaojifeixia 2026-03-19 12/600 2026-03-21 01:49 by 星空星月
[考研] 297求调剂 +9 戏精丹丹丹 2026-03-17 9/450 2026-03-21 01:49 by JourneyLucky
[考研] 271材料工程求调剂 +8 .6lL 2026-03-18 8/400 2026-03-21 00:58 by JourneyLucky
[考研] 一志愿西南交大，求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 260求调剂 +3 朱芷琳 2026-03-20 3/150 2026-03-20 20:35 by 学员8dgXkO
[考研] 289求调剂 +6 怀瑾握瑜l 2026-03-20 6/300 2026-03-20 20:30 by 学员8dgXkO
[考研] 一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 一志愿中国海洋大学，生物学，301分，求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 6/300 2026-03-19 17:18 by fei626-918
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 080500，材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
信息提示
关闭
请填处理意见
关闭