CI-NEB 计算过渡态的交上算不到一分钟自动错误断掉

chenxiankai

算一个分子在五层金的表面上扩散


我用了80个节点1920个核交的作业有四个images  换过节点数还是会自动断掉


下面是 outcar 的输出

k-point  1 :   0.0000 0.0000 0.0000  plane waves:  277047

maximum and minimum number of plane-waves per node :       159      136

maximum number of plane-waves:    277047
maximum index in each direction:
   IXMAX=   48   IYMAX=   28   IZMAX=   47
   IXMIN=  -48   IYMIN=  -28   IZMIN=  -47

NGX is ok and might be reduce to 194
NGY is ok and might be reduce to 114
NGZ is ok and might be reduce to 190

real space projection operators:
  total allocation   :     232666.38 KBytes
  max/ min on nodes  :       1351.52          0.00


parallel 3D FFT for wavefunctions:
    minimum data exchange during FFTs selected (reduces bandwidth)
parallel 3D FFT for charge:
    minimum data exchange during FFTs selected (reduces bandwidth)

下面是我的INCAR 和 KOINTS 文件

ISTART = 0
NWRITE = 2
ICHARG = 2

### Electronic Optimization SCF ###

  ENCUT = 400.0
  ISPIN = 2
   ISYM = 0
ISMEAR = 1
  SIGMA = 0.2
  EDIFF = 1.0E-4
   LDIAG = .TRUE.

### Ions relaxtion ###

    NSW = 600
IBRION = 3
  SMASS = -2
  POTIM = 0
EDIFFG = -0.07
   ISIF = 2

LORBIT = 11
   EMIN = -20.0
   EMAX = 10.0
  NEDOS = 3000

### Others ###
   GGA  = PBE
VOSKOWN = 1
   NPAR = 1
   IVDW = 11

   ALGO = Fast
   PREC = Accurate
  LREAL = A

### TS flags ###
ICHAIN = 0
IMAGES = 4
SPRING = -5
LCLIMB = .TRUE.
LTANGENTOLD = .FALSE.
LDNEB = .FALSE.
LNEBCELL = .FALSE.
IOPT = 3  ### or 7
### Write flags ###

  LWAVE = .FALSE.
LCHARG = .FALSE.
LPLANE = .TRUE.



KPOINTS


Automatic k-mesh
0
Gamma ! Gamma-centered grid
1 1 1 ! subdivisions N_1, N_2 and N_3 along recipr. vectors





哪位大侠能帮我一下 十分的感谢！！！！

夜夜衣

1、首先我先吐槽一下你的计算资源怎么会这么丰富，用80个节点将近2000核心来跑CINEB，真土豪。
2、我也是刚学vasp，感觉你的INCAR里面的参数也设置的太多了，感觉很多用默认的就完了，你可以参考一些顶级期刊的和你的课题类似的文献，一般会在si里面给出他的计算细节，电子和离子弛豫的参数可以参考他的做。
3.你用PBE计算应该是不需要设置VOSKOWN= 1，
4.不要把网上见到的所谓教程上的参数全部搬到你的INCAR。

chenxiankai

希望大家帮帮忙实质性的解决下这个问题。

夜夜衣

chenxiankai 发表于 2018-2-18 07:53
希望大家帮帮忙实质性的解决下这个问题。

这是你自己的课题，也肯定没什么人做过，熟悉你的课题的也不一定能看到你的帖子，难道你做这个课题，相关的计算的文献你一篇都没有看过吗？顶级期刊的si里面的细节都给的很丰富，先依葫芦画瓢都不会吗，你说额实质性解决难道是要别人帮你把INCAR写好了你直接拿来算就可以了吗~

hakuna

1.对照CNEB的说明，简化你的INCAR，把那些个诸如算DOS等等的参数注释掉，不明白的参数就用默认值
2.别拿资源不当回事，虽然不知道你算多大的体系，没必要这么奢侈
3.如果没有跑过CNEB，最好先跑个简单的体系，体会一下，再跑大的
4.多大的表面让你只设个Gamma点？
5.你提供的信息无法判断错误在何处

sobereva

对6L的不和谐发言进行了屏蔽。夜夜衣和hakuna都已经明显提供了显著有助于"实质性解决问题"的回复

在本论坛提问前注意阅读

在网上求助计算化学问题时的注意事项
http://bbs.keinsci.com/thread-57-1-1.html