计算化学公社

标题: VASP催化体系中一维材料大体系频率计算电子步收敛问题求助 [打印本页]

作者
Author: lkjmod 时间: 2024-11-6 19:48
标题: VASP催化体系中一维材料大体系频率计算电子步收敛问题求助
各位大佬，我在计算吸附氧原子，或者吸附羟基的频率计算时遇到了电子步不收敛的问题，INCAR和OSZICAR截图如下，该体系是个400多个原子胞大小为40×40×23的一维材料吸附氧或羟基，但频率计算时上千个电子步其误差一直卡在E-2下不去，ALGO = N或者F都尝试了，都不能解决问题，还麻烦各位大佬提提解决的建议

作者
Author: 白皛白 时间: 2024-11-6 20:55
给出的信息太少，很难说是哪里出了问题，请提供POSCAR。一个比较通用，且成功率较高的解决方案就是在单点计算或者计算频率的时候读取结构优化的WAVECAR。

作者
Author: lkjmod 时间: 2024-11-7 09:15
这是计算的POSCAR，方便给位参考一下

作者
Author: lkjmod 时间: 2024-11-7 09:17

白皛白发表于 2024-11-6 20:55
给出的信息太少，很难说是哪里出了问题，请提供POSCAR。一个比较通用，且成功率较高的解决方案就是在单点计 ...

大佬您好，我已经贴出POSCAR了您看看还需要什么信息吗，您提到的方法我也先尝试一下

作者
Author: 白皛白 时间: 2024-11-7 16:40
本帖最后由白皛白于 2024-11-7 16:44 编辑

lkjmod 发表于 2024-11-7 09:17
大佬您好，我已经贴出POSCAR了您看看还需要什么信息吗，您提到的方法我也先尝试一下

你好，我个人觉得有以下几个点可以注意一下，包括在结构优化的时候。如果只是解决频率计算SCF不收敛的话读取WAVECAR应该也足以解决问题了。
1，对于Mo元素的赝势文件可以看一下是不是用的pv或sv赝势，不同版本的VASPKIT默认的赝势是不同。
2，设置AMIX = 0.2 BMIX = 0.0001，AMIX可以多尝试几个，0.2收敛不了就试试0.15、0.1这样，这个参数可以帮助收敛。
3，一般来说碳化钼和Fe-N4这样的结构是有磁性的，自旋极化的问题有必要考虑。对于碳化钼来说可能要考虑DFT+U。
4，从晶格常数来看，感觉似乎在结构优化的时候没有优化晶格常数，对于碳纳米管来说优化轴向晶格常数是有必要的。
5，这个结构里面碳化钼和碳纳米管可能是不匹配的，打个比方就是这个碳纳米管可能容纳不下30个Mo原子的长度。强行装在同一个周期性盒子里可能导致这两个结构被拉伸或者压缩，这种结构上的不合理可能也是导致频率计算收敛不了的原因，感觉把这两个结构拆分开结构优化（优化晶格常数），然后再组装起来一起结构优化（优化晶格常数）可能更可靠一点。

作者
Author: lkjmod 时间: 2024-11-7 21:34

白皛白发表于 2024-11-7 16:40
你好，我个人觉得有以下几个点可以注意一下，包括在结构优化的时候。如果只是解决频率计算SCF不收敛的话 ...

感谢大佬的详细回复，学到了很多，我按照您的建议多测一下，看来要学习的东西还有很多

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