1、在最开始的结构优化过程中,如体系较大,可能有磁性,可以先设置较低的精度,只优化原子位置不优化晶胞。分别算有自旋情况和和无自旋的情况下的能量,然后选低的那个。在计算中istart=1 ICHARG用默认,一般只有在非自洽计算才特别指明ICHARG的值(非自洽计算要保持电荷密度不变)
2、出现在磁性计算下,电子步出现难收敛甚至是直接报错的情况;可以参考这篇文章或者这篇
Error EDDDAV: Call to ZHEGV failed
其主要的是通过调节AMIX的值来加快收敛,一般是AMIX=0.2 BMIX=0.0001 AMIX_MAG=0.8 BMIX_MAG=0.0001 可以称为线性混合
此外可以适当的增加NBANDS的值 同过grep NBANDS OUTCAR查看有多少个能带,增加到1.3或者2倍足够了.
3、关于判断体系是半导体还是金属,有一个经验参考
在初步试探计算时 设置 ISMEAR=0 和SIGMA=0.1 然后用grep EENTRO OUTCAR 将得到的EENTR/原子数 大于1meV为金属,小于0.1为半导体,
当然最准确的还是看DOS和能带。