VASP结构收敛正常结束但是力不收敛问题请教

KCM@SHU · 发表于 Post on 2024-6-27 21:51:19

本帖最后由 KCM@SHU 于 2024-6-27 21:57 编辑

各位老师好，我最近接触到VASP，我想问一下各位老师在VASP结构收敛的时候程序显示达到设定收敛值，但是在检查力的收敛的时候发现没有达到设定值，我设定的是EDIFFG = -0.01，但是实际受力明显大于0.01，这种情况一般可以怎么做呢？（如下图每一步的收敛图是按照这个链接：http://bbs.keinsci.com/forum.php ... 65&fromuid=1532里面老师的方法绘制的）（1）我尝试了将上次优化结果的构型再当成初始构型进行计算结果发现也没有达到设定值；（2）我看网上说VASP显示达到设定值但是力还没收敛的情况下可以尝试将POTIM调小，我把POTIM调小到0.2后还是没有收敛。针对VASP显示达到结构收敛但是力没收敛这种情况各位老师有什么好的建议吗？

乐平 · 发表于 Post on 2024-6-27 21:55:20

你用什么判断力没达到收敛标准？

KCM@SHU · 发表于 Post on 2024-6-27 22:03:11

乐平发表于 2024-6-27 21:55
你用什么判断力没达到收敛标准？

老师我是按照您之前的那个代码绘制的图判断的，我看图里面最后一步受力明显大于0.01，刚刚图片没上传上去，我修改了一下原贴。

乐平 · 发表于 Post on 2024-6-27 22:14:55

本帖最后由乐平于 2024-6-27 22:24 编辑

KCM@SHU 发表于 2024-6-27 22:03
老师我是按照您之前的那个代码绘制的图判断的，我看图里面最后一步受力明显大于0.01，刚刚图片没上传上去 ...

也没差太多……

另外，由于你的 KPOINTS 文件名没有S （你的是 KPOINT），于是，你的 OUTCAR 里 645 行警告说没找到 KPOINTS 文件，于是自动产生了 1 1 1 k-点

还有，你的赝势为什么都选的是 _GW 结尾的啊…… 你又不是算光学性质

Available pseudopotentials - VASP Wiki
GW potentials are recommended for calculations involving unoccupied states, e.g., computing optical properties or using many-body perturbation theory.

当然，这可能不是力没收敛的原因

KCM@SHU · 发表于 Post on 2024-6-27 22:36:25

乐平发表于 2024-6-27 22:14
也没差太多……

另外，由于你的 KPOINTS 文件名没有S （你的是 KPOINT），于是，你的 OUTCAR 里 645 ...

好的老师，我看一下，我最近刚接触VASP，赝势我当时选的时候就看文件后面带个NEW我就想着最新的就选了，我试试改一下KPOINTS和换个赝势试试，我把计算精度也提高一下试试看看能不能达到收敛标准，谢谢老师

卡开发发 · 发表于 Post on 2024-6-28 11:54:23

KCM@SHU 发表于 2024-6-27 22:36
好的老师，我看一下，我最近刚接触VASP，赝势我当时选的时候就看文件后面带个NEW我就想着最新的就选了， ...

我自己写了一小段代码，看了一下你最后一步的力是7e-3左右，下面的代码我没有考虑原子固定的情况，实际上要处理固定也不麻烦，欢迎坛友有兴趣自己来扩展。

import sys
import re
from io import StringIO
import numpy as np
"""
使用方法：python3 get_max_forces.py OUTCAR
"""
class VaspData(object):
def __init__(self, outcar) -> None:
self.outcar = outcar
@property
def max_forces(self):
pattern = '(.*?)'.join(['TOTAL-FORCE \(eV/Angst\)\n ?-+\n',
'-+\n +total drift'])
string = re.findall(pattern, self.outcar, re.S)
forces = [np.loadtxt(StringIO(s)).reshape(-1, 6)[:, 3:]
for s in string]
max_forces = [np.max(np.linalg.norm(f, axis=1)) for f in forces]
return max_forces
if __name__ == '__main__':
with open(sys.argv[1]) as f:
outcar = f.read()
outputs = VaspData(outcar)
for i, fmax in enumerate(outputs.max_forces):
print('%5d %.3e' % (i, fmax))

复制代码

乐平 · 发表于 Post on 2024-6-28 17:08:51

KCM@SHU 发表于 2024-6-27 22:36
好的老师，我看一下，我最近刚接触VASP，赝势我当时选的时候就看文件后面带个NEW我就想着最新的就选了， ...

抱歉，是我的代码里有个地方写错了，现在更新了

np.amax(np.linalg.norm(data, axis=-1), axis=-1)

复制代码

之前的轴弄错了，应该是 -1 轴

http://bbs.keinsci.com/forum.php ... 293265&fromuid=1532

乐平 · 发表于 Post on 2024-6-28 17:10:05

卡开发发发表于 2024-6-28 11:54
我自己写了一小段代码，看了一下你最后一步的力是7e-3左右，下面的代码我没有考虑原子固定的情况，实际上要 ...

卡老师读取表格的代码很优雅，值得学习

KCM@SHU · 发表于 Post on 2024-6-29 21:22:56

本帖最后由 KCM@SHU 于 2024-6-29 21:24 编辑

乐平发表于 2024-6-28 17:08
抱歉，是我的代码里有个地方写错了，现在更新了

好的，谢谢老师，老师我还想问一下NCORE这个命令是可以直接加到INCAR里面的吗？我再INCAR里面写上NCORE=16的时候，还是会出现关于NCORE=1的warning

乐平 · 发表于 Post on 2024-6-29 21:24:18

KCM@SHU 发表于 2024-6-29 21:22
好的，谢谢老师，老师我还想问一下NCORE这个命令是可以直接加到INCAR里面的吗？我再INCAR里面写上NCORE=1 ...

关于并行效率，看我之前的帖子

http://bbs.keinsci.com/forum.php ... 87&fromuid=1532

KCM@SHU · 发表于 Post on 2024-6-29 21:25:04

卡开发发发表于 2024-6-28 11:54
我自己写了一小段代码，看了一下你最后一步的力是7e-3左右，下面的代码我没有考虑原子固定的情况，实际上要 ...

谢谢老师分享

KCM@SHU · 发表于 Post on 2024-6-29 21:36:25

乐平发表于 2024-6-29 21:24
关于并行效率，看我之前的帖子

http://bbs.keinsci.com/forum.php?mod=redirect&goto=findpost&ptid=3 ...

老师我是想问一下NCORE是直接写进INCORE里面就可以修改计算方式还是要在其他文件里面修改呢？

乐平 · 发表于 Post on 2024-6-29 21:45:58

本帖最后由乐平于 2024-6-29 22:04 编辑

KCM@SHU 发表于 2024-6-29 21:36
老师我是想问一下NCORE是直接写进INCORE里面就可以修改计算方式还是要在其他文件里面修改呢？

KPAR， NCORE 当然是在 INCAR 里啊
看来你对 VASP 各个输入文件的基本内容不熟悉，建议先看看下面的视频
VASP输入文件
https://www.bilibili.com/video/B ... 828517ca6bef0a8956b

另外，仔细看我上面的帖子，对你的计算速度提升很有帮助

最后，建议用 VASPKIT 来准备你的输入文件 INCAR, POTCAR, KPOINTS
下面是VASPKIT软件开发者王伟老师介绍 VASPKIT 使用的视频（从安装到使用都有详细介绍）
https://www.bilibili.com/video/B ... 828517ca6bef0a8956b

KCM@SHU · 发表于 Post on 2024-6-29 21:58:33

乐平发表于 2024-6-29 21:45
KPAR， NCORE 当然是在 INCAR 里啊

仔细看我上面的帖子，对你的计算速度提升很有帮助

好的老师我看看，那他个这个关于NCORE的warning是不管怎么设置都会出现嘛？，我理解的是这个warning是提醒我NCORE=1计算效果很慢，但是我在INCAR里面设置了INCAR=16了

乐平 · 发表于 Post on 2024-6-29 22:01:34

KCM@SHU 发表于 2024-6-29 21:58
好的老师我看看，那他个这个关于NCORE的warning是不管怎么设置都会出现嘛？，我理解的是这个warning是提 ...

你真的仔细看了我的贴子吗？

你设置了 KPAR 和 NCOER 两个参数配合了吗？

另外，之前指出了你的 KPOINT 是错的，应该是 KPOINTS ，你改名字了吗？

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[VASP] VASP结构收敛正常结束但是力不收敛问题请教

评分 Rate