24小时热门版块排行榜

返回列表

chord7806

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 108.8
帖子: 31
在线: 107.3小时
虫号: 724541
注册: 2009-03-17
专业: 凝聚态物性 II ：电子结构

[求助] <求助>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 ]

回复此楼

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

第一性

» 猜你喜欢

A期刊撤稿已经有3人回复
职称评审没过，求安慰已经有34人回复
垃圾破二本职称评审标准已经有17人回复
回收溶剂求助已经有6人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
申请26博士已经有5人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
求助文献已经有3人回复
投稿返修后收到这样的回复，还有希望吗已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2012-05-29 16:11:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wanghouyuan

铜虫 (小有名气)

应助: 6 (幼儿园)
金币: 224.1
红花: 1
帖子: 75
在线: 15.2小时
虫号: 1144660
注册: 2010-11-12
性别: GG
专业: 理论和计算化学

请问你optimization和static cal在INCAR里的参数设置有何区别？

赞一下

回复此楼

子弹

2楼2012-05-29 16:50:33

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

chord7806

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 108.8
帖子: 31
在线: 107.3小时
虫号: 724541
注册: 2009-03-17
专业: 凝聚态物性 II ：电子结构

引用回帖:

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