求助：如何判断两种IBIRION获得的弛豫结构的优劣

shinkon

如果初始结构和平衡结构已经比较接近，使用IBRION=1或2都只需20~30步即可达到收敛标准，那么是否选取任一结构进行后续计算都可以？还是必须进一步比较静态总能等数据？以前都没怎么考虑过这个问题。。

卡开发发

如果SCF设置合理，理论上无论哪种方法收敛的结果得到的梯度最大值都会落在阈值内，这样再次用另一种方法算应该第一步仍然梯度落在阈值内就结束计算了。

但实际可能发生的是：
1、停在拐点，恰好满足梯度收敛条件，严格说这种只好去计算声子验证虚频。
2、对IBRION这种没步长控制的，有一定可能步长过大直接飞过了接近的极小点到别的结构。

shinkon

卡开发发 发表于 2023-11-24 16:55
如果SCF设置合理，理论上无论哪种方法收敛的结果得到的梯度最大值都会落在阈值内，这样再次用另一种方法算 ...

我自己做了几个小测试，对我的结构来说，用IBRION=1或IBRION=2，得到的CONTCAR用vesta打开，基本看不到显著差别；再做静态计算，能量差异基本上在0.00001eV/atom这个量级，我觉得基本应该可以视作没有差异了

卡开发发

shinkon 发表于 2023-11-24 19:25
我自己做了几个小测试，对我的结构来说，用IBRION=1或IBRION=2，得到的CONTCAR用vesta打开，基本看不到显 ...

那应该问题不大。只是极端情况下会出问题，结构如果离真的极小点足够近的时候问题比较少。