计算化学公社

标题: vasp计算能带报错，出现BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES怎... [打印本页]

作者
Author: wmg166 时间: 2024-5-19 20:46
标题: vasp计算能带报错，出现BAD TERMINATION OF ONE OF YOUR APPLICATION PROCESSES怎...
本帖最后由 wmg166 于 2024-7-6 16:03 编辑

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
===================================================================================

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

作者
Author: wmg166 时间: 2024-5-22 11:37

每天进步一点点发表于 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

作者
Author: wmg166 时间: 2024-5-22 11:39

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

作者
Author: 每天进步一点点 时间: 2024-5-22 12:19

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

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

作者
Author: 北大-陶豫 时间: 2024-5-22 15:12

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

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

作者
Author: wmg166 时间: 2024-5-23 08:05

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

谢谢大师

作者
Author: wmg166 时间: 2024-7-6 16:06

北大-陶豫发表于 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    吗？

作者
Author: 乐平 时间: 2024-7-6 16:33
本帖最后由乐平于 2024-7-6 16:38 编辑

你的问题在于

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 内容会变动，到时候读取的内容也不会不一样。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)