计算化学公社

标题: 用VASP优化结构遇到的问题 [打印本页]

作者
Author: 吴双 时间: 2021-10-13 18:14
标题: 用VASP优化结构遇到的问题
(, 下载次数 Times of downloads: 40)
想请教一下各位老师，我在用VASP优化一个结构，初始的时候EDIFF=1e-4 优化了80多个小时没优化出来，就停止了任务，把CONCAR保存成POSCAR,EDIFF=1e-3,在LOG文件出现了这个提示，想请问下大家还要继续扩大收敛线还是应该怎么做呢

作者
Author: 夙雪夜归 时间: 2021-10-14 09:11
也许应该说一下自己算的什么体系，大概多少原子，用了多少核心进行计算。这样大家也许能帮你看看，因为现有的信息还是挺有限的。

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

作者
Author: 吴双 时间: 2021-10-14 09:19

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

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

作者
Author: 夙雪夜归 时间: 2021-10-14 09:44

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

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

作者
Author: 吴双 时间: 2021-10-14 09:59

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

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

作者
Author: 吴双 时间: 2021-10-14 10:00

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

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

作者
Author: 夙雪夜归 时间: 2021-10-14 10:57

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

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

作者
Author: 吴双 时间: 2021-10-14 11:02

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

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

作者
Author: 夙雪夜归 时间: 2021-10-14 11:44

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

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

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

作者
Author: wsz 时间: 2021-10-15 10:35
使用Davidson算法能收敛吗（ALGO = Normal）

作者
Author: 吴双 时间: 2021-10-15 18:47

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

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

作者
Author: 吴双 时间: 2021-10-15 18:48

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

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

作者
Author: 吴双 时间: 2021-10-15 18:49

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

嗯嗯

作者
Author: 吴双 时间: 2021-10-15 18:54

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

不好意思刚才打错了，我用的是ALGO=FAST 没有收敛我用您说的这个尝试一下

作者
Author: 夙雪夜归 时间: 2021-10-15 19:01

吴双发表于 2021-10-15 18:47
是一个单胞结构，在上面锚定金属形成单原子催化剂，用于催化反应，因为我优化的单胞属于一个阴离子整体成 ...

所以你的意思是没有真空层吗，这样的话该如何进行催化反应的探究呢

作者
Author: 吴双 时间: 2021-10-15 19:07
是一个这样的结构，建立的20 20 20 的笼子 K点是111

作者
Author: 吴双 时间: 2021-10-15 19:11

夙雪夜归发表于 2021-10-15 19:01
所以你的意思是没有真空层吗，这样的话该如何进行催化反应的探究呢

您好在我上传的结构中整体结构是W6O19整体是-2价，加了两个个Na（紫色的原子），使其电中性，灰色的是锚定的过渡金属Zn,这个结构优化过程中电子步有过收敛但是收敛几次后，又开始不收敛，想跟您请教下这个应该怎么解决

作者
Author: archer 时间: 2021-10-15 21:57
修改IBRION

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