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

lei234

VASP报错，错误信息：Error EDDDAV: Call to ZHEGV failed. Returncode =  25 2  48
INCAR文件如下：
SYSTEM   = InAs
PREC     = accurate
EDIFF    = 1.0E-5
EDIFFG   = 1.0E-4
ENCUT    = 360.0 eV
! self-consistent calculation
ISTART   = 0
ICHARG   = 2
ALGO= Fast
POTIM    = 0.1
! Relaxation of lattice
NSW      = 0
IBRION   = -1
ISMEAR   = -5; SIGMA = 0.1
! local potenial
!LVTOT    = F
LWAVE = F
LCHARG = F
请高手指点，谢谢~

wth1219

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

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

wth1219

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

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

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

卡开发发

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

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

dou857

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

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

Cindydu123

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

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 =    但是结构能顺利收敛 ，请问这种情况下我优化收敛得到的结果能用吗 ? 合理吗？

北大-陶豫

我最近也遇到了同样的问题，但是尝试了楼上的方法（单核运行、使用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. 
    

复制代码

请大家批评指正。

DMYYCG

helpme 发表于 2020-11-5 10:47
ALGO=Very_Fast 跑完之后，需要把它去掉再跑一下吗？

没有必要重新跑

helpme

ALGO=Very_Fast 跑完之后，需要把它去掉再跑一下吗？

不取俗名

发现大家都在讨论原子占位的问题，我自己说一个我碰见过这样错误的原因吧。
例证：优化原胞，INCAR精度很高，K点999
用的核很多，结果报如上错误；
更正：用一个节点跑，成功解决

小梅的龙猫

同意楼上说的原子距离问题，我前两天建了一个参杂缺陷，可能是原子位子不对，优化精度设置较高，导致一直难以收敛，出现和楼主一样的错误，我设置了IALGO以后跑完了，不知道对不对，欢迎前辈们指导一二。

Jack

xuemei 发表于 2018-1-9 10:35
我也有同样的问题，麻烦请问你知道怎么判断占位有问题了吗?

我觉得是说原子之间距离有问题，不合理。你可以检查一下

xuemei

Jack 发表于 2017-10-1 13:00
请问怎么判断占位有问题，是不是看原子之间的距离是否正常？谢谢！

我也有同样的问题，麻烦请问你知道怎么判断占位有问题了吗?

xuemei

zjyx 发表于 2017-10-2 03:19
以前也出现过这种问题,不过我的那个貌似是对称性分析上出问题了.还有一次在改了ALGO=Normal后就好了

你好，你后来还遇到过类似问题吗？我改了ALGO以后，问题还是没有解决