用VASP优化结构遇到的问题

吴双

想请教一下各位老师，我在用VASP优化一个结构，初始的时候EDIFF=1e-4 优化了80多个小时没优化出来，就停止了任务，把CONCAR保存成POSCAR,EDIFF=1e-3,在LOG文件出现了这个提示，想请问下大家还要继续扩大收敛线还是应该怎么做呢

夙雪夜归

也许应该说一下自己算的什么体系，大概多少原子，用了多少核心进行计算。这样大家也许能帮你看看，因为现有的信息还是挺有限的。

ps:电子步收敛限VASP默认的是1e-4,vasp手册里写的精度高的话推荐设置1e-6,甚至1e-7，而文献里我看到的很多是使用1e-5，所以你再降低收敛限的话应当是不合理的

吴双

夙雪夜归 发表于 2021-10-14 09:11
也许应该说一下自己算的什么体系，大概多少原子，用了多少核心进行计算。这样大家也许能帮你看看，因为现有 ...

我算的是一个Lindqvist型多酸W6O19因为他自身是-2价，在计算时加了两个Na,使其电荷正确分布，在这个体系上锚定金属原子Zn，一共28个原子，用的八核计算

夙雪夜归

吴双 发表于 2021-10-14 09:19
我算的是一个Lindqvist型多酸W6O19因为他自身是-2价，在计算时加了两个Na,使其电荷正确分布，在这个体系 ...

你这是团簇、单胞、还是表面模型，盒子有多大；
看你贴的图电子步跑了100步后停了，也没有达到收敛限，我怀疑你结构本身可能就有问题，所以不太好收敛，可以仔细检查一下，或者贴张图上来

吴双

夙雪夜归 发表于 2021-10-14 09:44
你这是团簇、单胞、还是表面模型，盒子有多大；
看你贴的图电子步跑了100步后停了，也没有达到收敛限， ...

嗯嗯 老师我有一个问题想请教您一下，VASP优化结构的时候，怎么看优化是否正常呢，我现在看QSZICAR发现有的电子步收敛了，但是后面就又不收敛了

吴双

夙雪夜归 发表于 2021-10-14 09:44
你这是团簇、单胞、还是表面模型，盒子有多大；
看你贴的图电子步跑了100步后停了，也没有达到收敛限， ...

老师 我这是一个单胞 用的20 20 20 的盒子

夙雪夜归

吴双 发表于 2021-10-14 10:00
老师 我这是一个单胞 用的20 20 20 的盒子

不是老师，只是说一下我的理解，仅供参考：
1，如果是一个单胞的话，直接往优化好的单胞里加原子合适吗，你锚定Zn是为了做表面催化？
2.按我的理解，优化过程是否正常，其实就是看电子步和离子步是否达到收敛，还有你的结构在优化的过程中有没有散掉，或者被优化到一个明显不正确的状态；电子步和离子步是否收敛就是看是否在设定的步数内达到了收敛限，这方面你可以去看一下vsap手册上的NELM，NSW,EDIFF,EDIFFG条目的说明

吴双

夙雪夜归 发表于 2021-10-14 10:57
不是老师，只是说一下我的理解，仅供参考：
1，如果是一个单胞的话，直接往优化好的单胞里加原子合适吗 ...

要做一个单原子催化剂，比较内聚能和结合能看是否可以在单胞表面形成单原子

夙雪夜归

吴双 发表于 2021-10-14 11:02
要做一个单原子催化剂，比较内聚能和结合能看是否可以在单胞表面形成单原子

所以你优化的到底是单胞还是表面，做催化应该是表面吧，用单胞显然不合适了。

ymeng

我觉得没有必要把EDIFF=-3，-4都已经很宽松了。遇到这个报错，根据我个人经验按照提示把CONTCAR转化为POSCAR基本没成功过，先检查一下结构合理性

ymeng

我觉得没有必要把EDIFF=-3，-4都已经很宽松了。遇到这个报错，根据我个人经验按照提示把CONTCAR转化为POSCAR基本没成功过，先检查一下结构合理性

wsz

使用Davidson算法能收敛吗（ALGO   = Normal）

吴双

夙雪夜归 发表于 2021-10-14 11:44
所以你优化的到底是单胞还是表面，做催化应该是表面吧，用单胞显然不合适了。

是一个单胞结构，在上面锚定金属形成单原子催化剂，用于催化反应，因为我优化的单胞属于一个阴离子整体成-2价，加了两个Na使其电中性，这个电中性结构可以优化出来，但是在加金属比如Zn,Cu结构就不收敛，但是观察结构的话结构是合理的

吴双

wsz 发表于 2021-10-15 10:35
使用Davidson算法能收敛吗（ALGO   = Normal）

我用的ALGO=Accurate,结构和电子不收敛

吴双

ymeng 发表于 2021-10-15 10:08
我觉得没有必要把EDIFF=-3，-4都已经很宽松了。遇到这个报错，根据我个人经验按照提示把CONTCAR转化为POSCA ...

嗯嗯