计算化学公社

标题: 求助：VASP计算结构弛豫不收敛 [打印本页]

作者
Author: 木小虫427 时间: 2019-10-10 11:44
标题: 求助：VASP计算结构弛豫不收敛
本帖最后由木小虫427 于 2019-10-17 15:21 编辑

本人用VASP进行结构弛豫，进行到第二次弛豫的时候，后面一直在震荡：
(, 下载次数 Times of downloads: 42)
INCAR设置如下：
PREC=Accurate
ENCUT=400
ISTART=0
ICHARG=2
ISPIN=1
EDIFF=1E-5
GGA=PE
ISMEAR=0; SIGMA=0.1
EDIFFG=-0.01
NSW=100
IBRION=1
ISIF=3
POTIM=0.1
按照网上经验试过调低POTIM或者降低收敛精度，可是一直调成POTIM=0.01，EDIFFG=-0.02才行，
一步离子步就收敛了，但是看能量感觉并没有到应该的精度。
(, 下载次数 Times of downloads: 36)
也试过调高k点密度，从2*2*2 M改为6*6*6，结果电子步进行变得极慢，非常浪费时间。
请问各位前辈还有什么可以让收敛不震荡的经验分享吗？
还有想问一句，如果不做磁性计算的话，是不是不用设置MAGMOM的值？
谢谢各位前辈的帮助！

作者
Author: 卡开发发 时间: 2019-10-10 11:52
1、做不做自旋极化取决于体系本身是否有磁性，有磁性的体系就应该做自旋极化的计算。
2、不知道是否做了收敛性测试，足够高的计算精度才能够让结构优化能够做到较严格的收敛标准。
3、结构不太好的情况下尽量少用IBRION=1，IBRION=2每一步会通过线搜索确定步长，原则是会更稳定。

作者
Author: 木小虫427 时间: 2019-10-11 22:28

卡开发发发表于 2019-10-10 11:52
1、做不做自旋极化取决于体系本身是否有磁性，有磁性的体系就应该做自旋极化的计算。
2、不知道是否做了收 ...

谢谢您的解答！
请问收敛性测试是指E-Cut值的优化和K点优化吗?
还有看到一些经验都写IBRION=2发生震荡时最好换成1尝试，请问您知道为什么结构不太好的时候IBRION=2会更稳定吗？

作者
Author: 卡开发发 时间: 2019-10-12 02:06

木小虫427 发表于 2019-10-11 22:28
谢谢您的解答！
请问收敛性测试是指E-Cut值的优化和K点优化吗?
还有看到一些经验都写IBRION=2发生震荡 ...

1、是。
2、具体得看情况,。如果体系的结构特别离谱其实建议用IBRION=3；如果结构一般建议用IBRION=2，VASP的共轭梯度方法每一步会采用修正版的Brent算法来找最佳步长；如果结构的力很小，接近于势能面极小点，此时IBRION=1速度会快一些。而IBRION=1这种方法用在较远离极小点的结构上因为没有步长搜索的过程，一定可能会造成步长过长或过短。VASP里面尽管共轭梯度法每一步要搜索步长显得很慢，但原则上稳定性更好。

作者
Author: 飞翔的猪 时间: 2019-10-12 09:45
我一般参考： http://bbs.keinsci.com/thread-13065-1-1.html

作者
Author: 木小虫427 时间: 2019-10-14 09:22

卡开发发发表于 2019-10-12 02:06
1、是。
2、具体得看情况,。如果体系的结构特别离谱其实建议用IBRION=3；如果结构一般建议用IBRION=2，V ...

这样啊，明白了，谢谢您的解答~

作者
Author: 木小虫427 时间: 2019-10-14 09:22

飞翔的猪发表于 2019-10-12 09:45
我一般参考： http://bbs.keinsci.com/thread-13065-1-1.html

谢谢您的分享~

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