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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » <求助>VASP 结构优化收敛，静态计算不收敛。（INCAR KPOINTS POSCAR附上）
51/1返回列表
查看: 4456  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +6 Ncdx123456 2026-03-13 7/350 2026-03-16 07:44 by L135790
[考研] 311求调剂 +3 26研0 2026-03-15 3/150 2026-03-15 09:12 by JourneyLucky
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 309求调剂 +4 花与叶@ 2026-03-10 4/200 2026-03-14 21:26 by a不易
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 一志愿浙江大学0856材料与化工求调剂 +4 yansheng@211 2026-03-09 5/250 2026-03-14 02:10 by JourneyLucky
[考研] 一志愿北京化工大学材料与化工296分求调剂 +16 稻妻小编 2026-03-09 18/900 2026-03-14 02:00 by JourneyLucky
[考研] 一志愿天津大学，英一数二305分求调剂，四六级已过 +8 小小番的茄 2026-03-09 8/400 2026-03-14 01:53 by JourneyLucky
[考研] 312求调剂 +6 陌宸希 2026-03-10 6/300 2026-03-14 00:40 by JourneyLucky
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 332求调剂 +3 zjy101327 2026-03-11 6/300 2026-03-13 22:48 by JourneyLucky
[考研] 求材料调剂 085600英一数二总分302 前三科235 精通机器学习 一志愿哈工大 +4 林yaxin 2026-03-12 4/200 2026-03-13 22:04 by 星空星月
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 0856化工原理 +6 z2839474511 2026-03-10 6/300 2026-03-13 10:41 by houyaoxu
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[基金申请] 提交后的基金本子，已让学校撤回了，可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
[考研] 085602化工求调剂 +7 董boxing 2026-03-10 7/350 2026-03-10 17:07 by BruceLiu320
信息提示
关闭
请填处理意见
关闭