计算单个原子能量时如何选择空间群

白袍将军小火炉

在VASP计算中需要计算原子的能量（例如Si原子），但是下载cif时，有多种不同空间群的Si原子，请问根据什么来选择用哪种空间群的Si原子？

ChemG

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

白袍将军小火炉

ChemG 发表于 2022-7-18 14:28
原子哪有空间群，空间群源自于对称性，一个原子怎么转都对称。算原子也不用下cif，直接Si坐标 0 0 0，自己 ...

感谢您的回答，但是这样算电子占据可能是不对的，如图1；那么选一个cif的话，算总能量除以原子数目就可以了，我想用后者的思路，图1是一个原子算的电子占据情况，对于硅原子而言，有分数占据应该是不合理的，图2感觉占据态更贴近物理实际，所以我想了一阵，不知道选择空间群的依据是什么，因为硅原子的cif有很多不同种类的cif。谢谢您，希望您赐教！

ChemG

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

原来是这样，如果是金刚石结构的，那就是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/.

白袍将军小火炉

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

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

白袍将军小火炉

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

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