24小时热门版块排行榜

返回列表

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

lywiailyw

木虫 (小有名气)

应助: 4 (幼儿园)
金币: 1980.4
散金: 14
红花: 1
帖子: 267
在线: 120.7小时
虫号: 756914
注册: 2009-04-24
性别: GG
专业: 理论和计算化学

[求助] vasp团簇结构弛豫出现如下死循环该怎么办

OSZICAR文件中重复出现如下部分：
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.627599350918E+02 -0.62760E+02 0.36453E-07 128 0.869E-08 0.126E-08
RMM: 2 -0.627599351284E+02 -0.36652E-07 0.00000E+00    0 0.000E+00 0.747E-09
RMM: 3 -0.627599351284E+02 0.18190E-10 0.00000E+00    0 0.000E+00 0.522E-09
RMM: 4 -0.627599351284E+02 -0.21828E-10 0.00000E+00    0 0.000E+00 0.261E-09
RMM: 5 -0.627599351284E+02 0.23647E-10 0.00000E+00    0 0.000E+00
296 F= -.62759935E+02 E0= -.62759935E+02  d E =-.215858E-02  mag= 19.0000
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.627599350917E+02 -0.62760E+02 0.36429E-07 128 0.769E-08 0.142E-08
RMM: 2 -0.627599351284E+02 -0.36687E-07 0.00000E+00    0 0.000E+00 0.775E-09
RMM: 3 -0.627599351284E+02 -0.15461E-10 0.00000E+00    0 0.000E+00 0.415E-09
RMM: 4 -0.627599351284E+02 0.27285E-10 0.00000E+00    0 0.000E+00 0.398E-09
RMM: 5 -0.627599351284E+02 0.13642E-10 0.00000E+00    0 0.000E+00
297 F= -.62759935E+02 E0= -.62759935E+02  d E =-.215858E-02  mag= 19.0000
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.627599350921E+02 -0.62760E+02 0.36445E-07 128 0.665E-08 0.105E-08
RMM: 2 -0.627599351284E+02 -0.36322E-07 0.00000E+00    0 0.000E+00 0.658E-09
RMM: 3 -0.627599351284E+02 0.54570E-11 0.00000E+00    0 0.000E+00 0.372E-09
RMM: 4 -0.627599351284E+02 -0.12733E-10 0.00000E+00    0 0.000E+00 0.210E-09
RMM: 5 -0.627599351284E+02 -0.26375E-10 0.00000E+00    0 0.000E+00
298 F= -.62759935E+02 E0= -.62759935E+02  d E =-.215858E-02  mag= 19.0000
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.627599350919E+02 -0.62760E+02 0.36489E-07 128 0.485E-08 0.607E-09
RMM: 2 -0.627599351284E+02 -0.36489E-07 0.00000E+00    0 0.000E+00 0.338E-09
RMM: 3 -0.627599351284E+02 0.36380E-11 0.00000E+00    0 0.000E+00 0.192E-09
RMM: 4 -0.627599351284E+02 0.36380E-11 0.00000E+00    0 0.000E+00 0.136E-09
RMM: 5 -0.627599351284E+02 -0.40018E-10 0.00000E+00    0 0.000E+00
299 F= -.62759935E+02 E0= -.62759935E+02  d E =-.215858E-02  mag= 19.0000
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.627599350920E+02 -0.62760E+02 0.36453E-07 128 0.502E-08 0.735E-09
RMM: 2 -0.627599351284E+02 -0.36433E-07 0.00000E+00    0 0.000E+00 0.428E-09
RMM: 3 -0.627599351284E+02 0.72760E-11 0.00000E+00    0 0.000E+00 0.249E-09
RMM: 4 -0.627599351284E+02 -0.27285E-11 0.00000E+00    0 0.000E+00 0.286E-09
RMM: 5 -0.627599351284E+02 0.90949E-11 0.00000E+00    0 0.000E+00
300 F= -.62759935E+02 E0= -.62759935E+02  d E =-.215858E-02  mag= 19.0000
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.627599350920E+02 -0.62760E+02 0.36456E-07 128 0.717E-08 0.869E-09
RMM: 2 -0.627599351284E+02 -0.36400E-07 0.00000E+00    0 0.000E+00 0.467E-09
RMM: 3 -0.627599351284E+02 0.12733E-10 0.00000E+00    0 0.000E+00 0.274E-09
RMM: 4 -0.627599351284E+02 -0.36380E-11 0.00000E+00    0 0.000E+00 0.168E-09
RMM: 5 -0.627599351284E+02 -0.18190E-10 0.00000E+00    0 0.000E+00
301 F= -.62759935E+02 E0= -.62759935E+02  d E =-.215858E-02  mag= 19.0000

我这个结构是在初始优化的基础上进行的更精确优化，INCAR文件一并贴上供高手分析：
SYSTEM = sys
ISTART = 0
ICHARG = 2
INIWAV = 1
PREC = Accurate
ISMEAR = 0
SIGMA = 0.01
#POTIM = 0.1
NELMIN = 5
ISIF = 2
ENCUT = 280
EDIFF = 0.0001
EDIFFG = -0.01
IBRION = 2
ALGO = Fast
NSW = 2000
LREAL = Auto
ISPIN = 2
LORBIT = 11
ISYM = 0
LWAVE = .FALSE.
LCHARG = .FALSE.

还请有经验的虫们指点为何会出现OSZICAR中所示的循环现象，谢谢！

回复此楼

» 猜你喜欢

职称评审没过，求安慰已经有49人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复
垃圾破二本职称评审标准已经有17人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
三无产品还有机会吗已经有6人回复

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

1楼 2012-06-07 21:55:36

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

lywiailyw

木虫 (小有名气)

应助: 4 (幼儿园)
金币: 1980.4
散金: 14
红花: 1
帖子: 267
在线: 120.7小时
虫号: 756914
注册: 2009-04-24
性别: GG
专业: 理论和计算化学

引用回帖:

2楼: Originally posted by mousekingadv at 2012-06-08 19:49:46
楼主把EDIFFG注释掉试试……这看起来像是EDIFF达到精度了，EDIFFG却达不到。EDIFFG=0.01很苛刻了。可能的原因是ENCUT=280太小了（瞎猜，楼主没办法可以试试）

而且，注释掉力的话，默认值就变成能量的十倍，也就是0.001了，反而更小