vasp运行报错ZPOTRF解决方案

sky

官方论坛里主要有两个解决方案：
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号才有这问题，还是其它计算集群都有问题？希望能有人反馈一下你们编译和计算的情况，谢谢！

comchem

我也遇到这个 问题了，应该是编译的问题，现在还没有找到解决办法。