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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » <求助>VASP 结构优化收敛，静态计算不收敛。（INCAR KPOINTS POSCAR附上）
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 4479  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料求调剂 +11 一样YWY 2026-04-05 11/550 2026-04-05 23:36 by 来看流星雨10
[考研] 机械专硕274求调剂，不挑专业学校 +6 泛泛2333 2026-04-05 7/350 2026-04-05 23:20 by chyhaha
[考研] 一志愿 南京航空航天大学 ，080500材料科学与工程学硕 +10 @taotao 2026-03-30 10/500 2026-04-05 17:57 by jj987
[考研] 材料专硕 调剂 +14 CXN123456 2026-04-03 14/700 2026-04-05 17:18 by Hdyxbekcb
[考研] 081200-11408-276学硕求调剂 +4 崔wj 2026-04-04 5/250 2026-04-05 14:06 by imissbao
[考研] 求调剂 +3 电气小神童 2026-04-04 3/150 2026-04-05 10:17 by barlinike
[考研] 323求调剂 +8 李佳乐1 2026-04-04 8/400 2026-04-04 22:26 by hemengdong
[考研] 材料调剂 +10 懒羊羊轻置玉臀 2026-04-02 11/550 2026-04-04 21:56 by laoshidan
[考研] 319求调剂 +4 星星不眨眼喽 2026-04-03 4/200 2026-04-04 16:25 by 中飞院空管学院�
[考研] 土木304求调剂 +4 兔突突突， 2026-03-31 4/200 2026-04-04 13:34 by 1753564080
[考研] 一志愿北京科技大学材料工程085601，求调剂 +17 cdyw 2026-04-02 18/900 2026-04-04 11:14 by w_xuqing
[考研] 一志愿南昌大学324求调剂 +13 hanamiko 2026-04-01 13/650 2026-04-03 18:30 by ls刘帅
[考研] 085600 295分求调剂 +19 W55j 2026-03-30 23/1150 2026-04-03 09:53 by 千千运气
[考研] 材料考研调剂 +10 Gs大王 2026-04-02 10/500 2026-04-03 09:47 by 遗忘消失的灆
[考研] 279求调剂 +5 傅文秋 2026-04-02 5/250 2026-04-02 18:10 by 笔落锦州
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 09:52 by yulian1987
[考研] 296求调剂 +4 汪！？！ 2026-03-31 7/350 2026-04-01 22:04 by 客尔美德
[考研] 310分求调剂 +4 成功上岸wang 2026-04-01 4/200 2026-04-01 20:35 by liu823948201
[考研] 267求调剂 +13 uiybh 2026-03-31 13/650 2026-04-01 10:25 by 探123
[考研] 求0861交通运输专硕or材料专硕调剂 +4 勒布朗@ 2026-03-31 4/200 2026-04-01 09:54 by 一只好果子?
信息提示
关闭
请填处理意见
关闭