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

24小时热门版块排行榜    

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

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的回帖
查看全部 5 个回答

wanghouyuan

铜虫 (小有名气)
请问你optimization和static cal在INCAR里的参数设置有何区别？
赞一下 回复此楼
子弹
2楼2012-05-29 16:50:33
已阅   回复此楼   关注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的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +8 Ldrop2023 2026-08-13 8/400 2026-08-16 11:28 by jiege890620
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 25/1250 2026-08-16 10:53 by hmhminy
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
关闭
请填处理意见
关闭