氧原子能量的正确计算

1、Ex8 氧原子能量的正确计算前面的扩展练习中，要求大家结合自己的化学知识，判断计算结果是否正确。相比很多同学已经知道了氧原子含有磁性，在计算用需要考虑进去。在本节练习中，大师兄教给你如何分子氧原子在轨道中的占据情况，判断计算是否合理，并正确进行氧原子能量的计算。氧原子的电子构型是 1s22s22p4 轨道占据情况如图：所以氧原子含有两个单电子，具有磁性，且磁矩为 2 B。看一下前面的计算结果：OSZICAR ：能量信息：电子占据的情况：（OAR 中查看）Band1 含有两个电子，对应的 2s2；Band 2-4 为三个简并态，每个能带上有 1.33 个电子，对应的 2p4 。但在这里，电子的占

2、据是不正确的：因为 2px 轨道上含有两个电子，而不是把 2px 上的一个电子平均分配到3 个p 轨道上；此外，通过这个电子的占据情况，不能正确得到氧原子的磁矩信息。那么怎么解决这个问题呢？分两步解决问题：第一步：由于氧原子具有磁性，自旋极化计算要进行考虑，在 VASP计算中，需要在 INCAR 中添加一项： ISPIN = 2ISPIN 的取值有两个，1 和 2 ：1 代表不考虑自旋极化，是 VASP 的默认值；2 代表打开自旋极化该选项，计算中考虑。自旋极化什么时候该加，什么时候用默认值呢？要计算的体系具有磁性时，必须要打开自旋极化。一般来说，单原子的计算，O2 分子（基态为三重态），以及

3、含 Fe， Co，Ni 的体系需要考虑自旋极化。添加该选项后，重新计算，结果如下：INCAROSZICAR可以看出，加了自旋极化后，电子迭代步数增加了，这是因为体系中的电子被强制分成了 alpha 和 beta 两种，并分别进行计算，因此增加了计算量。体系的磁矩大小在 OSZICAR 中查看，对应右下角的mag=该项，从这里可以看到，磁矩信息与已知的 2B 吻合，这一点是合理的，正确的。OUITCAR 中的能量信息：考虑自旋极化后，体系的能量降低了。重点通过查看体系的电子占据情况，判断该计算是否合理：从图中可以看出，结果由两个 spin 组成： spin component 1 和2：spin

4、 component 1种含有 4 个电子，另一个种含有 2 个电子：参考该图，可以看出 spin component 1 为自旋向上的 alpha 电子，另一个是自选向下的 beta 电子。但是：1）ponent 1 中，能带 2-4 中 3 个 alpha 电子的能量是一样的，这不正确，因为其中一个 2px 轨道中的 alpha 电子已经成对，能量要比 py 和 pz 低；描述失败！2）ponent 2 中，另一个 beta 电子还是平均分配在了三个能带上面。电子占据依然不合理。上面的结果是由于体系的高对称性导致的简并所造成的。 8x 8 x 8 的立方体格子，在这里具有高阶的点群对称性：

5、 O_h，见OAR在 GGA 泛函中，为了获得体系更低的能量，对于原子来说，通常会采用一种 symmetry broken solution 的处理方法。但是在 VASP 计算中，体系的对称性则是通过晶胞来获得，即把这个晶胞当成一个原子来处理，因此需要手动改变晶胞的形状来消除对称性造成的简并结果。（此处的解释可能有些牵强，有的可以查找与 symmetrybroken solution 相关的文献）INCAR， KPOS， POTCAR 等均保持不变。修改CAR，如下：通过改变格子在三个方向的大小，降低对称性。大师兄还尝试了三斜的格子，结果和这个是一致的。（有的可以尝试下）下面看一下改变晶胞为

6、7.5 x 8.0 x 8.9 后的计算结果：OSZICAROSZICAR 中电子迭代步数又增加了，这是因为对称性的降低增加了计算量。磁矩为 2B，这一点跟前面一样，正确！体系的对称性：从 O_h 降低到了 D_2h能量信息：跟前面能量相比，可以看出，对称性降低后，体系的能量进一步降低了。电子占据情况：分析：1） 在 alpha 电子（spin component 1）中，能带 1 是 2s 中的alpha 电子；能带 2-3 为两个简并轨道，对应 py 和 pz 电子，能带 4 位px 中的 alpha 电子。这里能带 4 的能量应该比 2 和 3 要低，但结果恰恰相反，描述不合理；2） 在

7、 beta 电子（spin component 2）中，能带 1 是 2s 的beta 电子，能带 2 为px 中的 beta 电子，此时该 beta 电子占据了一个轨道，而不是分布在三个 p 轨道上，结果是合理的。计算到这里，已经是 VASP 官网对于单原子能量最为完整和准确的计算了。虽然上面 px 的 alpha 电子比 py 和 pz 能量高，至少电子的轨道电子占据情况是正确的，且磁矩和简并的错误已经消除。由此可见， VASP 对于单原子的电子占据情况（至少对于 O 原子来说），虽然取得了部分理想结果，但还是有的地方。此外，体系的对称性可以通过 ISYM 来控制（扩展阅读 ISIM）,I

8、SIM0 的时候，不考虑对称性进行计算，大师兄也尝试过了将 ISYM=0 应用在 8x8x8 （test3）以及 7.5x8.0 x8.9 （test4） 的格子里进行计算。对于电子占据情况来说，均得到了与前面 7.5x8.0 x8.9 一致的结果。但是使用 7.5x8.0 x8.9所得到的能量更低（test4 的能量比 test3 低）。使用 7.5x8.0 x8.9 格子时，ISYM=0 （test4）和采用默认值，得到的能量结果是相同的，因此，改变晶胞的对称性在计算单原子能量的时候是必须的，单独通过 ISIM 这一项不考虑对称性是远远不够的。扩展阅读：1）ISIM：2）重现本节，并自己认

9、真分析结果，查看电子占据情况，体系的磁性，以及能量信息。总结：1）知道什么时候使用自旋极化计算，怎么进行自旋极化计算；2）知道为什么自旋极化后计算量增加了；3）知道在哪里查看体系的磁性信息；4）知道为什么对称性降低后计算量增加了；5）会分析单原子中电子轨道的占据情况；判断结果是否合理；6）知道通过改变晶格大小调节体系的对称性从而消除简并对能量的影响；7）知道ISYM 这一个参数，在计算单原子能量中，ISYM 的作用远远不及改变晶格大小；8）知道vasp 也存在自己局限性，不要完全相信。PS：如果喜欢大师兄的文章，欢迎关注 ， ，转发。此外， 专注于科研思维的碰撞与科研生活的 ，本书中已经详细解释或者指明的易出错部分，不建议在群中继续 ，请大家认真学习并主动积极地去思考和练习。历史练习Learn-VASP-The-Hard-Way（ex7） Learn-VASP-The-Hard-Way（ex6） Learn-VASP-The-Hard-Way（ex5） Learn-VASP-The-Hard-Way（ex4） Learn-VASP-The-Hard-Way（ex3） Learn-VASP-Th