计算化学公社

标题: vasp当ISIF=3时无法收敛 [打印本页]

作者
Author: IanTsian 时间: 2021-4-29 20:48
标题: vasp当ISIF=3时无法收敛
本帖最后由 IanTsian 于 2021-4-29 22:43 编辑

如题，使用vasp取ISIF=3优化二维材料晶格常数时，固定Z轴，然而在优化过程中某些离子受力无法降下去，因此无法收敛，请问一下有什么改进方法（ISIF=2时是没问题的）。INCAR文件和OPTCELL文件如图。

作者
Author: sobereva 时间: 2021-4-29 21:58
贴图方式不对，其他人看不到。重新编辑帖子并仔细看置顶的新社员必读贴了解怎么正确贴图，此问题在这里还特意强调了：http://bbs.keinsci.com/thread-18961-1-1.html

作者
Author: IanTsian 时间: 2021-4-29 22:44

sobereva 发表于 2021-4-29 21:58
贴图方式不对，其他人看不到。重新编辑帖子并仔细看置顶的新社员必读贴了解怎么正确贴图，此问题在这里还特 ...

抱歉抱歉，初次发帖，已经改正了。

作者
Author: agent99 时间: 2021-4-30 12:42
改回ISIF=2，然后手动扫描X和Y？可能会比较费时

作者
Author: IanTsian 时间: 2021-4-30 13:49

agent99 发表于 2021-4-30 12:42
改回ISIF=2，然后手动扫描X和Y？可能会比较费时

ISIF=2手动扫描是没问题的，使用ISIF=3主要是想偷个懒。

作者
Author: MisakaFaction 时间: 2021-4-30 16:41
如果受力是在低位反复震荡，可以用IBRION=1续算试试
如果一开始就有问题，可能模型的初始结构不太合理，用ISIF=2优化后的结果试试
直接用SCAN收敛有抽风的概率，请重视手册的建议如下图：
(, 下载次数 Times of downloads: 46)

作者
Author: 卡开发发 时间: 2021-5-1 09:26

IanTsian 发表于 2021-4-30 13:49
ISIF=2手动扫描是没问题的，使用ISIF=3主要是想偷个懒。

1、手动扫描会偏慢，如果沿着一个轴向可以用无导数的线搜索方法直接找极小，这个需要自己写程序。
2、OPTCELL本质上只是解决了stress上的问题，force上面还得自己想办法，我觉得可尝试更大的真空层、尽可能通过修正消除层间的偶极矩。
3、楼上提到的SCAN涉及到的问题是SCF的，在这不牵扯到这个问题（除非你的SCF确实没收敛）。但确实要小心的是，SCAN所需要的密度网格（ENAUG或者也可以通过PREC）要更高，建议做好测试，否则力会产生一些不必要的噪声引起构型优化震荡。按道理一个构型优化过程应该是由寻找下降方向和优化最优步长两个过程来完成的，IBRION=1没有搜索最有步长这个过程，如不是构型接近平衡构型不太建议使用，而IBRION=2在每个CG步骤都要使用线搜索来找最优步长，理论上来说稳定性会更好。

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