VASP计算精度和是否收敛问题

KCM@SHU

（1）各位老师，我最近在自学VASP，然后我看网上的教程了解到VASP收敛主要看能量和力是否收敛，我做结构优化后，我在看输出文件OSZICAR的时候发现第一步离子步的时候电子步计算到了最大的步数还是没有收敛，但是程序也没有停下来，反而继续开始了后面的离子步，但是我看最后离子步都跑完了好像d E好像都没有达到我设置的EDIFF，然后我看OUTCAR里面的TOTAL-FORCE也没有达到EDIFFG，虽然最后显示reached required accuracy - stopping structural energy minimisation是不是其实还是没收敛呀？
（2）各位老师我想问一下，我看网上说有时候也要看d eps是否小于EDIFF这个是判断什么呢？还有像我这种程序显示收敛但是能量和力都没有达到设定的情况，是不是有办法可以像Gaussian里面一样读取最后一步作为初始步继续优化，达到收敛的概率会大一点呢？
（3）各位老师我想问一下VASP对于高精度计算的定义主要是指哪一类呢？我主要是结构优化、计算能带、PDOS、吸附能这一块，因为我看我输出文件里面warning建议对于精度不高的计算可以用LREAL= Auto

卢夫

qqqm 发表于 2024-9-10 15:49
楼主你好,请问你做结构优化用时多久呢,有多少原子和核数

就我自己的体系100个原子 64核 一般优化四五个小时左右  我一般是先gam版粗优化 然后再std版优化

shinnashiori

haonan 发表于 2025-9-12 18:36
请问解决了吗，我也遇到了这个问题

额，楼上本来就有解决办法啊，，，

haonan

目前还不会计算 发表于 2024-9-28 22:02
File "check_VASP_convergency.py", line 203
    raise FileNotFoundError(f"The file {file_path.nam ...

请问解决了吗，我也遇到了这个问题

目前还不会计算

好滴好滴，，确实是我看错了，，我现在改到了正确的了，可以了。感谢感谢

乐平

目前还不会计算 发表于 2024-9-29 16:26
确实是版本很低，我默认安装的那个Anaconda，然后版本只有2.7。
Python 2.7.3 :: Anaconda 1.4.0 (64-bi ...

…… 无语了

你是在我推荐的清华的镜像网站里下载的吗？

如果是，你不仔细看看日期吗？ 选最近的日期对应的 Anaconda3 下载啊

目前还不会计算

乐平 发表于 2024-9-29 10:34
你默认的 Python 不会是 2.7 吧……

你在命令行输入 Python 看看输出的版本号

确实是版本很低，我默认安装的那个Anaconda，然后版本只有2.7。
Python 2.7.3 :: Anaconda 1.4.0 (64-bit)

乐平

目前还不会计算 发表于 2024-9-28 22:02
File "check_VASP_convergency.py", line 203
    raise FileNotFoundError(f"The file {file_path.nam ...

你默认的 Python 不会是 2.7 吧……

你在命令行输入 Python 看看输出的版本号

最好是截图，免得乱猜也猜不到你具体的情况……

目前还不会计算

乐平 发表于 2024-9-28 20:31
因为你的 python 是纯 Python，没有其他库

建议安装 Anaconda3，几乎常用的库都包含了，你不需要自己 ...

File "check_VASP_convergency.py", line 203
    raise FileNotFoundError(f"The file {file_path.name} does not exist.")
                                                                       ^
SyntaxError: invalid syntax
我安装好了库，现在又变成这样了，不知道为啥

乐平

目前还不会计算 发表于 2024-9-28 19:56
我输入：python3 check_VASP_convergency.py 之后，就直接显示这个，这是为啥呀，大佬
Traceback (most  ...

因为你的 python 是纯 Python，没有其他库

建议安装 Anaconda3，几乎常用的库都包含了，你不需要自己额外安装其他库。

从清华的镜像网站下载就可以，链接如下（比直接从国外的网站下载要快）
https://mirrors.tuna.tsinghua.edu.cn/help/anaconda/

目前还不会计算

乐平 发表于 2024-6-21 16:31
你不需要取点，只用看趋势就可以了

如果实在不放心，可以点图上工具栏里的“放大镜”工具，然后框选最 ...

我输入：python3 check_VASP_convergency.py 之后，就直接显示这个，这是为啥呀，大佬
Traceback (most recent call last):
  File "check_VASP_convergency.py", line 9, in <module>
    import numpy as np
ModuleNotFoundError: No module named 'numpy'

KCM@SHU

qqqm 发表于 2024-9-10 15:49
楼主你好,请问你做结构优化用时多久呢,有多少原子和核数

这个不太好说，我一般一次优化5-10个小时，一般优化一次不一定成功，用的核数是28，原子大概在70-200左右

qqqm

楼主你好,请问你做结构优化用时多久呢,有多少原子和核数

乐平

KCM@SHU 发表于 2024-6-20 12:37
好的老师，我用您的代码画出来图形了，我再琢磨一下咋样可以抓点，我不太确定力有没有达到我设定值，感觉 ...

你不需要取点，只用看趋势就可以了

如果实在不放心，可以点图上工具栏里的“放大镜”工具，然后框选最大受力的末尾十来个点。框选的时候尽量靠近点，这样纵坐标会自动缩放。你会看到纵坐标的刻度变成 10^-3 左右。足以说明最大受力满足你设置的 EDIFFG = -0.02

KCM@SHU

乐平 发表于 2024-6-20 10:16
我的 Python 脚本里用了 Numpy 来求最大受力。核心的求解公式就一行

好的老师，我用您的代码画出来图形了，我再琢磨一下咋样可以抓点，我不太确定力有没有达到我设定值，感觉好像在我设定值附近波动，老师计算出来图形上面那个能量是不是就是可以用来计算吸附能的自由能呢？