VASP报错 Error EDDDAV: Call to ZHEGV failed. Returncode = 25 2 48

北大-陶豫

我最近也遇到了同样的问题，但是尝试了楼上的方法（单核运行、使用ALGO=FAST、使用ALGO=VERYFAST等）均未成功。经反复尝试发现，若计算读取了之前的旋轨耦合计算任务产生的WAVECAR，或当前计算量很小（ENCUT极低，K点数极低）时可以正常计算。我的解决方法如下：
（1）        使用极少的K点数（如1x1x1）、极低的ENCUT（如40eV），完成一次成功的旋轨耦合计算。这步计算可以用mpirun多核并行计算。
（2）        读取任务（1）产生的WAVECAR，进行正常的计算。
我猜测题主这一报错可能的原因：
若任务没有读取之前任务的WAVECAR，就会自己猜一个初猜然后进行 SCF 迭代计算，为了减少计算量，迭代的前几步是 non-self-consistent（？）的，之后再切换到正常的 SCF 迭代。对于旋轨耦合计算任务，这步“切换”需要较大量的计算资源（很吃硬盘？很吃内存？我也没搞清楚），所以只有在使用极少的K点数（如1x1x1）、极低的ENCUT（如40eV）的时候才能成功。

下面附上我旋轨耦合计算任务的 INCAR 输入文件模板。

1. Global Parameters
   
2. ISTART =  1            (Read existing wavefunction, if there)
   
3. # ISPIN  =  1            (Non-Spin polarised DFT)
   
4. # ICHARG =  11         (Non-self-consistent: GGA/LDA band structures)
   
5. LREAL  = .FALSE.       (Projection operators: automatic)
   
6. ENCUT  =  400        (Cut-off energy for plane wave basis set, in eV)
   
7. PREC   =  Accurate   #(Precision level: Normal or Accurate, set Accurate when perform structure lattice relaxation calculation)
   
8. LWAVE  = .TRUE.        (Write WAVECAR or not)
   
9. LCHARG = .TRUE.        (Write CHGCAR or not)
   
10. ADDGRID= .TRUE.        (Increase grid, helps GGA convergence)
    
11. # LVTOT  = .TRUE.      (Write total electrostatic potential into LOCPOT or not)
    
12. # LVHAR  = .TRUE.      (Write ionic + Hartree electrostatic potential into LOCPOT or not)
    
13. # NELECT =             (No. of electrons: charged cells, be careful)
    
14. # LPLANE = .TRUE.      (Real space distribution, supercells)
    
15. # NWRITE = 2           (Medium-level output)
    
16. # KPAR   = 2           (Divides k-grid into separate groups)
    
17. # NGXF    = 300        (FFT grid mesh density for nice charge/potential plots)
    
18. # NGYF    = 300        (FFT grid mesh density for nice charge/potential plots)
    
19. # NGZF    = 300        (FFT grid mesh density for nice charge/potential plots)
    
20. 
    
21. #ALGO = Fast
    
22. 
    
23. Static Calculation
    
24. ISMEAR =  0            (gaussian smearing method)
    
25. SIGMA  =  0.05         (please check the width of the smearing)
    
26. LORBIT =  11           (PAW radii for projected DOS)
    
27. NEDOS  =  2001         (DOSCAR points)
    
28. NELM   =  60           (Max electronic SCF steps)
    
29. EDIFF  =  1E-06        (SCF energy convergence, in eV)
    
30. LSORBIT = .TRUE.
    
31. 
    
32. NCORE = 4
    
33. 
    

复制代码

请大家批评指正。

yuzj41d

jiewei 发表于 2014-12-24 10:03
vasp运行结构优化时，出现错误：Error EDDDAV: Call to ZHEGV failed. Returncode =  80 5 6

两种方法解 ...

jiewei老师，我想请问您一下，我在ALGO=Fast下运行优化时，时常出现 Error EDDRMM: Call to ZHEGV failed. Returncode =    但是结构能顺利收敛 ，请问这种情况下我优化收敛得到的结果能用吗 ? 合理吗？

Cindydu123

题主的分两步走，我试了了一下，也会出现其他问题，报错信息为：WARNING: Sub-Space-Matrix is not hermitian in DAV           59
  8.637407541478907E-005
我认为是应该是第一步K点取的太疏了，解波函数有问题。个人认为此方法不适用我的体系，还是要K点和晶格空间合适为好，需另寻方法。

dou857

卡开发发 发表于 2014-12-24 11:25
这问题以前没碰到过，按照网上的资料(应该很好找，百度都能找到，别人转载来的mail list)，要么就是结构有 ...

老师您好，vasp计算HSE06时出现forrtl: severe (174): SIGSEGV, segmentation fault occurred，网上说是修改stack，但是我修改了还是出现这个问题，这个如何解决那

卡开发发

dou857 发表于 2024-3-28 22:06
老师您好，vasp计算HSE06时出现forrtl: severe (174): SIGSEGV, segmentation fault occurred，网上说是 ...

暂时没办法判断，你可以先看看OUTCAR有没有其他异常，以及尝试ulimit -s unlimited。

wth1219

除了使用ALGO=FAST和VERYFAST之外，再提供两个略带玄学的解决方案。

一是用NCORE的平方数作为核数，我在128核的节点（NCORE=16）上跑出现这个问题，但把两个节点并起来用256核跑就没有这个问题了。

二是如若在FPMD中遇到这个问题可以添加IWAVPR=11，这会影响计算速度，理论上对于计算结果应该没什么影响。

wth1219

dou857 发表于 2024-3-28 22:06
老师您好，vasp计算HSE06时出现forrtl: severe (174): SIGSEGV, segmentation fault occurred，网上说是 ...

我曾经在128核的节点（NCORE=16）上跑出现这个问题，但把两个节点并起来用256核跑就没有这个问题了。