掺杂体系形成能如何计算？

龙腾万里

    以g-C3N4为例，A和B为掺杂元素，前者是间隙掺杂，掺杂到空腔中；后者是元素替换，替换g-C3N4中的C元素，形成能=E(ABCN)+E(C)-E(CN)-E(A)-E(B); 其中ABCN为掺杂后的物质，CN为g-C3N4，A、B、C为元素。使用程序为VASP。

问题1：此公式是否正确？答：计算形成能时一定是所有生成物总能减去所有反应物总能量，前后总原子数一定要相同。

问题2：A、B、C的能量如何计算？是计算单个原子的能量，还是寻找其最稳定结构，计算晶胞的能量，然后除以晶胞内原子数求其平均值后的能量？答：对于形成能选后者，结合能一般选前者。

问题3：对于空位缺陷，其形成能又该如何计算？答：例如计算总原子数为64的ZnO超胞中的氧空位形成能，生成物是Zn64O63和1/2O2，反应物是Zn64O64，所以Ef=E(Zn64O63)+1/2E(O2)-E(Zn64O64)，前后总原子数一定要相同。