乘风万里 发表于 2022-7-18 22:34
计算单原子的基态能量以及电子占据时，需要关掉对称性 ISYM = 0  ，使用Gaussian smearing ：ISMEAR = 0。 ...

非常感谢您详细的解答，对我解决占据计算问题起到很大作用，之前没有把SIGMA调的那么小，也没有设置电子自洽步数足够大，再次感谢您

ChemG 发表于 2022-7-18 16:23
原来是这样，如果是金刚石结构的，那就是Fd-3m (227号)

感谢您的解答，我现在懂了根据结构来选择空间群了

计算单原子的基态能量以及电子占据时，需要关掉对称性 ISYM = 0  ，使用Gaussian smearing ：ISMEAR = 0。以及十分小的SIGMA，SIGMA=0.001，this is necessary if atomic orbitals are almost degenerated。同时打开自旋极化ISPIN = 2，得到区分上下自旋的占据。肯定需要很多步电子自洽，NELM=300。这几个参数可以很好的算出3s2,3p2占据。其它注意事项请仔细阅读VASP算原子的官方教程：https://www.vasp.at/wiki/index.php/Calculation_of_atoms，也可以看看大师兄的教程：https://www.bigbrosci.com/2017/10/30/ex08/.

白袍将军小火炉 发表于 2022-7-18 14:59
感谢您的回答，但是这样算电子占据可能是不对的，如图1；那么选一个cif的话，算总能量除以原子数目就可以 ...

原来是这样，如果是金刚石结构的，那就是Fd-3m (227号)

原子哪有空间群，空间群源自于对称性，一个原子怎么转都对称。算原子也不用下cif，直接Si坐标 0 0 0，自己再写个晶胞就得了。