计算化学公社

标题: vasp运行报错ZPOTRF解决方案 [打印本页]

作者
Author: sky 时间: 2018-9-26 17:52
标题: vasp运行报错ZPOTRF解决方案
本帖最后由 sky 于 2018-9-27 10:50 编辑

官方论坛里主要有两个解决方案：
1. 输入的结构可能有问题，比如某两个原子距离太近
2. 改算法，有人反映IALGO=48就能运行

但我最近碰到的错误用以上方法却无法解决，最后发现是编译的问题。
在天河二号上，vasp5.3以及后面的版本无法使用自旋计算（ISPIN=2）,否则报错ZPOTRF：
Error EDDDAV: Call to ZHEGV failed. Returncode = 25 2 48
或者是报错：
BRMIX: very serious problems
the old and the new charge density differ
old charge density: 1292.00000 new 1292.06824
0.826E+02

该问题几年前就有人反映过：
http://muchong.com/t-9732230-1-authorid-2613852

我这次发现vasp5.2版本不受影响，但5.3和5.4版本都会报错。如果要继续使用5.3以上版本，编译的时候需要把所有优化选项改为O1。搜索makefile文件，把O2、O3(大写字母O)统统改为O1就行，亲测有效，但计算的速度会比原来慢一倍左右。
(管理员回复的解决方案为更改makefile中的FREE参数，由原来的
FREE = -free -names lowercase
改成
FREE = -free -heap-arrays)
亲测该方法无效，我看过管理员的makefile文件，用的依然是O1优化。

另外我很想知道是只有天河2号才有这问题，还是其它计算集群都有问题？希望能有人反馈一下你们编译和计算的情况，谢谢！

作者
Author: comchem 时间: 2023-10-25 23:02
我也遇到这个问题了，应该是编译的问题，现在还没有找到解决办法。

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