首先检查结构,有没有不稳定的结构,例如悬空键、原子距离太接近(保存任务时会提示)。如果结构没有问题,进行后面的尝试。
如果是GGA,可以使用BP86,GGA-PBE的收敛性有时候不是很好。如果BP86仍然不收敛,可以考虑使用metaGGA-SCAN。
如果体系有磁性,则指定自旋极化具体数值,有助于收敛。
检查一下不收敛的计算结果ADF LOGO > Output,找到
B A N D E N E R G Y R A N G E S The energies below are in Hartree. band min(up) max(up) occup(up) min(dn) max(dn) occup(dn) 2 -0.78949 -0.78949 1.00000 -0.76259 -0.76259 1.00000 3 -0.78726 -0.78726 1.00000 -0.75745 -0.75745 1.00000 4 -0.78496 -0.78496 1.00000 -0.75116 -0.75116 1.00000 5 -0.20317 -0.20317 1.00000 -0.19882 -0.19882 1.00000 6 -0.12008 -0.12008 1.00000 -0.07778 -0.07778 0.00000 ………………………… …………………………
检查最高占据轨道和最低占据轨道的能量。例如上例中,occup(up)也就是自旋向上的电子占据数,能量为-0.12008 Hatree(注意本例是BAND计算分子,所以没有能带,所以“能带”上的最高点和最低点是同一个数值)的轨道,占据了一个电子;能量为-0.07778 Hatree的轨道为occup(dn)(自旋向下),占据了0个电子。
指定占据的设置,参考:BAND里面如何指定占据方式
SCF关键词里面,增加Mixing设置,但mixing的数值不宜过大
SCF .... Mixing 0.05 ! more conservative mixing ... End
…表示其他关键词
增加DIIS关键词:
Diis
Variant LISTi ! invoke the LISTi method
End