vasp计算能带报错，出现BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES怎...

wmg166

vasp计算能带。先结构优化，再自洽（ISTART=0, SIGMA=0.01,ICORELEVEL=1,LVTOT=T,LVHAR=T），最后非自洽(ISTART=1,ICHARG=11)，PREC与静态一致，ENCUT不变。KPATH.in改为KPOINTS...计算能带时报错，INCAR 和  KPOINTS如最下方所示，


K-Path Generated by VASPKIT.
   20
Line-Mode
Reciprocal
   0.0000000000   0.0000000000   0.0000000000     GAMMA         
   0.0000000000   0.5000000000   0.0000000000     X              

   0.0000000000   0.5000000000   0.0000000000     X              
   0.5000000000   0.5000000000   0.0000000000     M              

   0.5000000000   0.5000000000   0.0000000000     M              
   0.0000000000   0.0000000000   0.0000000000     GAMMA         

   0.0000000000   0.0000000000   0.0000000000     GAMMA         
   0.0000000000   0.0000000000   0.5000000000     Z              

   0.0000000000   0.0000000000   0.5000000000     Z              
   0.0000000000   0.5000000000   0.5000000000     R              

   0.0000000000   0.5000000000   0.5000000000     R              
   0.5000000000   0.5000000000   0.5000000000     A              

   0.5000000000   0.5000000000   0.5000000000     A              
   0.0000000000   0.0000000000   0.5000000000     Z              

   0.0000000000   0.5000000000   0.0000000000     X              
   0.0000000000   0.5000000000   0.5000000000     R              

   0.5000000000   0.5000000000   0.0000000000     M              
   0.5000000000   0.5000000000   0.5000000000     A  




出现这样的提示怎么办：
LDA part: xc-table for Pade appr. of Perdew
found WAVECAR, reading the header
  number of k-points has changed, file:     6 present:   180
  trying to continue reading WAVECAR, but it might fail
WARNING: stress and forces are not correct
POSCAR, INCAR and KPOINTS ok, starting setup
FFT: planning ...
reading WAVECAR
the WAVECAR file was read successfully

===================================================================================
=   BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES
=   PID 19114 RUNNING AT localhost.localdomain
=   EXIT CODE: 9
=   CLEANING UP REMAINING PROCESSES
=   YOU CAN IGNORE THE BELOW CLEANUP MESSAGES



INCAR 内容 ：
Global Parameters
ISTART =  1            (Read existing wavefunction; if there)
# ISPIN =  2           (Spin polarised DFT)
ICHARG =  11         (Non-self-consistent: GGA/LDA band structures)
LREAL  = .FALSE.       (Projection operators: automatic)
# ENCUT  =  400        (Cut-off energy for plane wave basis set, in eV)
PREC   =  Normal       (Precision level)
LWAVE  = .TRUE.        (Write WAVECAR or not)
LCHARG = .TRUE.        (Write CHGCAR or not)
ADDGRID= .TRUE.        (Increase grid; helps GGA convergence)
LVTOT  = .TRUE.      (Write total electrostatic potential into LOCPOT or not)
LVHAR  = .TRUE.      (Write ionic + Hartree electrostatic potential into LOCPOT or not)
# NELECT =             (No. of electrons: charged cells; be careful)
# LPLANE = .TRUE.      (Real space distribution; supercells)
# NPAR   = 4           (Max is no. nodes; don't set for hybrids)
# NWRITE = 2           (Medium-level output)
# KPAR   = 2           (Divides k-grid into separate groups)
NGX    = 500         (FFT grid mesh density for nice charge/potential plots)
NGY    = 500         (FFT grid mesh density for nice charge/potential plots)
NGZ    = 500         (FFT grid mesh density for nice charge/potential plots)
ICORELEVEL=1
NBANDS=368

Static Calculation
ISMEAR =  0            (gaussian smearing method)
SIGMA  =  0.01         (please check the width of the smearing)
LORBIT =  11           (PAW radii for projected DOS)
NEDOS  =  2001         (DOSCAR points)
NELM   =  60           (Max electronic SCF steps)
EDIFF  =  1E-08        (SCF energy convergence; in eV)



===================================================================================
   Intel(R) MPI Library troubleshooting guide:
      https://software.intel.com/node/561764
===================================================================================

每天进步一点点

我前几天出现过这个问题，好像和内存相关。因为我的体系原子很多，我选的高对称点路径太多了，老师让我选一部分计算就行，如果实在要全部路径，可以分开计算。再就是调小撒点密度，然后增加计算节点，调整并行参数。在我的计算里我调小了路径，用了三个节点，加了参数NCORE🟰8，就好算很多了。而且这个LWAVE LCHARG如果后续不用可以关掉，不然占内存。我也是小白有的地方描述不准确，请见谅

wmg166

每天进步一点点 发表于 2024-5-22 02:04
我前几天出现过这个问题，好像和内存相关。因为我的体系原子很多，我选的高对称点路径太多了，老师让我选一 ...

我当前的KPOINTS如下，   是不是 删掉最后几行就可以调小撒点密度了 ？
K-Path Generated by VASPKIT.
   20
Line-Mode
Reciprocal
   0.0000000000   0.0000000000   0.0000000000     GAMMA         
   0.0000000000   0.5000000000   0.0000000000     X              

   0.0000000000   0.5000000000   0.0000000000     X              
   0.5000000000   0.5000000000   0.0000000000     M              

   0.5000000000   0.5000000000   0.0000000000     M              
   0.0000000000   0.0000000000   0.0000000000     GAMMA         

   0.0000000000   0.0000000000   0.0000000000     GAMMA         
   0.0000000000   0.0000000000   0.5000000000     Z              

   0.0000000000   0.0000000000   0.5000000000     Z              
   0.0000000000   0.5000000000   0.5000000000     R              

   0.0000000000   0.5000000000   0.5000000000     R              
   0.5000000000   0.5000000000   0.5000000000     A              

   0.5000000000   0.5000000000   0.5000000000     A              
   0.0000000000   0.0000000000   0.5000000000     Z              

   0.0000000000   0.5000000000   0.0000000000     X              
   0.0000000000   0.5000000000   0.5000000000     R              

   0.5000000000   0.5000000000   0.0000000000     M              
   0.5000000000   0.5000000000   0.5000000000     A              

wmg166

wmg166 发表于 2024-5-22 11:37
我当前的KPOINTS如下，   是不是 删掉最后几行就可以调小撒点密度了 ？
K-Path Generated by VASPKIT.
...

用另一台服务器计算出现这样的提示，咋办 ：
reading WAVECAR
random initialization beyond band          369
the WAVECAR file was read successfully
WARNING: dimensions on CHGCAR file are different
ERROR: charge density could not be read from file CHGCAR for ICHARG>10
ERROR: charge density could not be read from file CHGCAR for ICHARG>10
ERROR: charge density could not be read from file CHGCAR for ICHARG>10
ERROR: charge density could not be read from file CHGCAR for ICHARG>10
ERROR: charge density could not be read from file CHGCAR for ICHARG>10

每天进步一点点

wmg166 发表于 2024-5-22 11:37
我当前的KPOINTS如下，   是不是 删掉最后几行就可以调小撒点密度了 ？
K-Path Generated by VASPKIT.
...

可以试试删掉最后几行的一些高对称点，然后第二行20修改成10看看能不能算？这个新的报错我不太懂，查了一下好像是说INCAR的PREC和ENCUT必须和静态自洽的一致。

北大-陶豫

wmg166 发表于 2024-5-22 11:37
我当前的KPOINTS如下，   是不是 删掉最后几行就可以调小撒点密度了 ？
K-Path Generated by VASPKIT.
...

第二行的20表示每条能带上撒20个点，一般不需要这么多，调小这个数就可以避免爆内存的问题。虽然撒点数减少之后计算比较粗略，但可以在 ~/.vaspkit 配置文件里可以设置能带外推，另外，得到粗略的能带图后可以找出自己感兴趣的少数几条能带，用较高的分辨率重算。

wmg166

北大-陶豫 发表于 2024-5-22 15:12
第二行的20表示每条能带上撒20个点，一般不需要这么多，调小这个数就可以避免爆内存的问题。虽然撒点数减 ...

谢谢大师

wmg166

北大-陶豫 发表于 2024-5-22 15:12
第二行的20表示每条能带上撒20个点，一般不需要这么多，调小这个数就可以避免爆内存的问题。虽然撒点数减 ...

大神，KPOINTS我改成这样：
K-Path Generated by VASPKIT.
   20
Line-Mode
Reciprocal
   0.0000000000   0.0000000000   0.0000000000     GAMMA         
   0.0000000000   0.5000000000   0.0000000000     X              

   0.0000000000   0.5000000000   0.0000000000     X              
   0.5000000000   0.5000000000   0.0000000000     M              

   0.5000000000   0.5000000000   0.0000000000     M              
   0.0000000000   0.0000000000   0.0000000000     GAMMA         计算能带还是如一楼所示那样的提示，咋办，改成 ： 0.0000000000   0.0000000000   0.0000000000     GAMMA         
   0.0000000000   0.5000000000   0.0000000000     X   
好像还是同样的提示，咋办 还需要修改 High symmetry points 文件吗，High symmetry points 文件也改成：0.0000000000   0.0000000000   0.0000000000     GAMMA         
   0.0000000000   0.5000000000   0.0000000000     X       吗？

乐平

你的问题在于

found WAVECAR, reading the header
  number of k-points has changed, file:     6 present:   180

你最前一步自洽计算的时候，INCAR 里设置了
NGX    = 500         (FFT grid mesh density for nice charge/potential plots)
NGY    = 500         (FFT grid mesh density for nice charge/potential plots)
NGZ    = 500         (FFT grid mesh density for nice charge/potential plots)
ICORELEVEL=1
NBANDS=368

这些没有？
如果没有，这些就不要设置。


如果设置了，这些就不要变。


另外，你的 ENCUT 是注释状态，根本有设置…… 如果 PREC 一直保持不变（比如你这里是 Normal）倒还说得过去；如果 PREC 一会儿设置 Accuracy ，一会儿设置 Normal，那么 ENCUT 不设置的话，会随着 PREC 的不同而变化，导致你的 WAVECAR 内容会变动，到时候读取的内容也不会不一样。