24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ApolloYang

木虫 (正式写手)

应助: 15 (小学生)
金币: 2209.6
散金: 114
红花: 4
帖子: 390
在线: 236.9小时
虫号: 582657
注册: 2008-07-25
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] 求助：vasp结构优化不能收敛

优化材料是只含C H O的晶体
POSCAR
$i
8.710999488830570    0.000000000000000    0.000000000000000
0.000000000000000    5.169000148773190    0.000000000000000
0.000000000000000    0.000000000000000    10.962000846862800
H  C  N
8  20 20
Direct
0.2440000176429750  -0.1260000020265580 0.1790000051259990
0.0150000005960464 0.4250000119209290 0.1700000017881390
  -0.2440000176429750 0.1260000020265580 0.6790000051259990
  -0.0150000005960464  -0.4250000119209290 0.6700000017881390
0.7440000176429750 0.1260000020265580 0.1790000051259990
0.5150000005960465  -0.4250000119209290 0.1700000017881390
0.2559999823570250  -0.1260000020265580 0.6790000051259990
……
  ……
……

驰豫计算的INCAR
ENCUT = 400
ISTART = 0
ICHARGE = 2
NBANDS=300
ISMEAR = 0
SIGMA = 0.2
NSW = 100
IBRION = 3
ISIF = 5
IALGO=48
#ALGO= Fast
#ISYM = 0
POTIM = 0.2
EDIFF = 1E-4
EDIFFG = -1E-3
PREC = Accurate

问题出现在：
驰豫步数为100步，但是算完几圈之后，每次驰豫迭代圈数就很小，而且不收敛。
请问这样都优化是否有效？有什么解决方法：
（IBRION＝1，3等都试过，均出现这种情况）
……
……
RMM: 8 -0.371093682477E+03 -0.10907E-06 -0.36825E-05  1557 0.774E-03
4 F= -.37109368E+03 E0= -.37100276E+03  d E =-.236048E+00
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.371095287506E+03 -0.16051E-02 -0.22090E-02  2432 0.209E-01 0.128E-01
RMM: 2 -0.371095076068E+03 0.21144E-03 -0.17902E-04  2225 0.178E-02 0.141E-01
RMM: 3 -0.371094841359E+03 0.23471E-03 -0.61025E-05  1821 0.918E-03 0.557E-02
RMM: 4 -0.371094754099E+03 0.87260E-04 -0.70320E-05  1476 0.105E-02
5 F= -.37109475E+03 E0= -.37100357E+03  d E =-.107162E-02
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.371096827810E+03 -0.19865E-02 -0.91344E-02  2432 0.423E-01 0.254E-01
RMM: 2 -0.371095937496E+03 0.89031E-03 -0.74723E-04  2432 0.364E-02 0.279E-01
RMM: 3 -0.371095012793E+03 0.92470E-03 -0.23466E-04  2240 0.172E-02 0.112E-01
RMM: 4 -0.371094658377E+03 0.35442E-03 -0.30030E-04  2426 0.213E-02 0.300E-02
RMM: 5 -0.371094641050E+03 0.17327E-04 -0.15466E-05  1451 0.539E-03
6 F= -.37109464E+03 E0= -.37100313E+03  d E =0.113049E-03
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.371095354007E+03 -0.69563E-03 -0.97841E-03  2432 0.141E-01 0.862E-02
RMM: 2 -0.371095293345E+03 0.60662E-04 -0.91573E-05  2055 0.126E-02
7 F= -.37109529E+03 E0= -.37100411E+03  d E =-.652295E-03
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.370481562311E+03 -0.69726E-06 -0.74509E-05  1330 0.179E-02 0.875E-03
RMM: 2 -0.370481560887E+03 0.14238E-05 -0.66862E-06 896 0.735E-03
8F= -.37048156E+03 E0= -.67048156E+03  d E =-.719179E+00
   N    E                   dE          d eps    ncg    rms       rms(c)
RMM: 1 -0.370481561444E+03 0.86728E-06 -0.31566E-06 813 0.553E-03 0.318E-03
RMM: 2 -0.370481561069E+03 0.37484E-06 -0.50675E-07 480 0.372E-03
  9 F= -0.37048156E+03 E0= -.67048156E+03  d E =-.719179E+00
…………
如这次OSZCAR到第6次驰豫时候，迭代圈数就只有两步了，算到100步也未收敛

[ Last edited by ApolloYang on 2011-10-10 at 16:32 ]

回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

1楼 2011-10-10 16:08:32

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

rxd337

木虫 (正式写手)

养殖场经理

应助: 13 (小学生)
金币: 2475.2
红花: 1
帖子: 372
在线: 117.9小时
虫号: 256580
注册: 2006-06-03
性别: GG
专业: 材料物理化学

【答案】应助回帖

★
uuv2010(金币+1): 多谢提示 2011-10-13 09:14:23

我觉得问题可能在于两个方面：
首先是你设置计算参数的技术问题，例如IBRION=3，IALGO=48等，这些计算设置可能是对某些体系有更快的计算，但是vasp的每个tag都是有默认的设置的，这些默认的设置往往是比较稳健的。我曾经调试IALGO算法时发现48在计算大体系时可能会有些问题，你的体系不小接近50个原子呢。所以我建议最好还是用vasp默认的设置来计算：

ENCUT = 400
ISTART = 0
ICHARG= 2
NBANDS=300
ISMEAR = 0
SIGMA = 0.2
NSW = 100
IBRION = 2
ISIF = 5
#IALGO=48
#ALGO= Fast
#ISYM = 0
POTIM = 0.2
#EDIFF = 1E-4
#EDIFFG = -1E-3
PREC = Accurate
如果还有问题，那么就可能不是计算参数的问题，可能在于你的晶胞不够“优良”，vasp确实在能势面上找不到最优结构，那就需要核查下晶胞的结构。