fux 发表于 2023-10-24 15:02
您好，我的IBRION=2，您想说的是ISMEAR=0是么？我记得以前看到过气体分子/原子体系最好设置ISMEAR=0 SIGM ...

是ISMEAR=0.
对气体分子/原子体系推荐设置ISMEAR=0以及较小的SIGMA的原因在于避免轨道出现非整数占据。通常气体分子/原子结合1*1*1的gamma-centered K点计算，此时计算出的能量可直接看作是分子的能级，而无需考虑能带结构。可以在OUTCAR中搜索k-point     1 :，检查每次SCF收敛后各能级占据情况。通常正确的计算要求没有非整数占据。

ahxb 发表于 2023-10-24 13:23
正常，vasp在笛卡尔坐标下优化，本来收敛就比量子化学程序慢，还需要考虑到镜像分子间的微弱相互作用。
另 ...

您好，我的IBRION=2，您想说的是ISMEAR=0是么？我记得以前看到过气体分子/原子体系最好设置ISMEAR=0 SIGMA=0.01的说法，一直不知道为什么。对于您提到的‘检查能级占据情况’，我在静态和DOS计算中检查过，但弛豫的我看了一下不太明白，比如我之前做的一个无机晶体的计算，静态和DOS计算中DOSCAR是能够看到在电子数=NELECT附近可能存在一个gap，但是弛豫计算时候最大的电子数都不到NELECT，那弛豫过程的能级占据应该怎么检查呢？

正常，vasp在笛卡尔坐标下优化，本来收敛就比量子化学程序慢，还需要考虑到镜像分子间的微弱相互作用。
另外关于你的INCAR有个问题： 对有机分子取IBRION=0, SIGMA=0.1未必合适，需要检查能级占据情况。