计算化学公社

标题: VASP报错 Error EDDDAV: Call to ZHEGV failed. Returncode = 25 2 48 [打印本页]

作者
Author: lei234 时间: 2014-12-24 09:34
标题: VASP报错 Error EDDDAV: Call to ZHEGV failed. Returncode = 25 2 48
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
请高手指点，谢谢~

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

两种方法解决：
（1）POSCAR中晶格占位有问题，认真核对一下，这种错误可能性比较大！
（2）加ALGO=VERY_FASE，使用RMM方法收敛，可能还会报类似的错误，但能继续计算！

作者
Author: 卡开发发 时间: 2014-12-24 11:25
这问题以前没碰到过，按照网上的资料(应该很好找，百度都能找到，别人转载来的mail list)，要么就是结构有问题，要么就是lapack相关的数学库或是编译的模块出了问题。

如果原子过近不妨按照楼上的说法ALGO=V试试先跑一段时间或者重新建模型试试，如果是后一种问题只好重新编译。

作者
Author: lei234 时间: 2014-12-24 12:32

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

两种方法解 ...

确实是POSCAR的坐标有问题，不过奇怪的是我没有改POSCAR,只是加了ALGO=Very_Fast,也正常结束了，但是结果应该不能用吧？

作者
Author: jiewei 时间: 2014-12-24 13:16

lei234 发表于 2014-12-24 12:32
确实是POSCAR的坐标有问题，不过奇怪的是我没有改POSCAR,只是加了ALGO=Very_Fast,也正常结束了，但是结果 ...

可以用，你建模的时候原子挨得太近，一般情况就会出现不收敛，如果你用粗糙点的方法它就过去了。
这就是为什么我们都是建议模型建好后，先用小方法或者说是力场预优化，或者用MD跑几圈，提取结构。这样就不会出现你的错误了。

作者
Author: Jack 时间: 2017-10-1 13:00

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

两种方法解 ...

请问怎么判断占位有问题，是不是看原子之间的距离是否正常？谢谢！

作者
Author: zjyx 时间: 2017-10-2 03:19
本帖最后由 zjyx 于 2017-10-2 03:23 编辑

以前也出现过这种问题,不过我的那个貌似是对称性分析上出问题了.还有一次在改了ALGO=Normal后就好了

作者
Author: mylifecg 时间: 2017-10-2 08:59

jiewei 发表于 2014-12-24 13:16
可以用，你建模的时候原子挨得太近，一般情况就会出现不收敛，如果你用粗糙点的方法它就过去了。
这就 ...

您好，请问这种warning得到的结果可靠吗？我在超算上运行同样的结构就不会出现类似的warning，但是在我的服务器上就会出现。

作者
Author: xuemei 时间: 2018-1-9 10:35

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

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

作者
Author: xuemei 时间: 2018-1-9 10:35

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

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

作者
Author: Jack 时间: 2018-1-9 18:22

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

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

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

作者
Author: 不取俗名 时间: 2018-2-20 13:33
发现大家都在讨论原子占位的问题，我自己说一个我碰见过这样错误的原因吧。
例证：优化原胞，INCAR精度很高，K点999
用的核很多，结果报如上错误；
更正：用一个节点跑，成功解决

作者
Author: helpme 时间: 2020-11-5 10:47
ALGO=Very_Fast 跑完之后，需要把它去掉再跑一下吗？

作者
Author: DMYYCG 时间: 2020-12-1 15:01

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

没有必要重新跑

作者
Author: 北大-陶豫 时间: 2022-5-17 08:07
本帖最后由北大-陶豫于 2022-5-17 08:13 编辑

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

Global Parameters
ISTART = 1 (Read existing wavefunction, if there)
# ISPIN = 1 (Non-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 = Accurate #(Precision level: Normal or Accurate, set Accurate when perform structure lattice relaxation calculation)
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)
# NWRITE = 2 (Medium-level output)
# KPAR = 2 (Divides k-grid into separate groups)
# NGXF = 300 (FFT grid mesh density for nice charge/potential plots)
# NGYF = 300 (FFT grid mesh density for nice charge/potential plots)
# NGZF = 300 (FFT grid mesh density for nice charge/potential plots)
#ALGO = Fast
Static Calculation
ISMEAR = 0 (gaussian smearing method)
SIGMA = 0.05 (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-06 (SCF energy convergence, in eV)
LSORBIT = .TRUE.
NCORE = 4

复制代码

请大家批评指正。

作者
Author: yuzj41d 时间: 2022-5-17 16:22

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

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

作者
Author: dou857 时间: 2024-3-28 22:06

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

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

作者
Author: 卡开发发 时间: 2024-3-29 09:16

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

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

作者
Author: wth1219 时间: 2024-9-23 17:29
除了使用ALGO=FAST和VERYFAST之外，再提供两个略带玄学的解决方案。

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

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

作者
Author: wth1219 时间: 2024-9-23 17:29

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

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

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

Powered by Discuz! X3.3