同一种材料，换一组参考能量，形成能就可能差出不少。问题往往出在元素参考态、原子数换算或计算设置上。VASP输出的是所建模型的总能量，形成能需要把目标材料与对应参考体系的能量相减。因此，开始计算前应先确定研究对象、参考反应和结果单位。

　　一、VASP怎么计算材料形成能

　　以元素单质为参考时，形成能表示这些元素生成目标材料的能量变化。计算对象可以使用原胞、常规晶胞或超胞，但相减时必须按实际原子数配平。

　　1、写清形成能表达式

　　①确认POSCAR中的元素种类及数量，写出目标材料由元素参考态生成的反应。

　　②用目标晶胞总能量减去各元素“原子数×每原子参考能量”，得到该晶胞的形成能。

　　③根据用途换算为每化学式或每原子结果，并标明eV/化学式或eV/atom。

　　例如，MgO晶胞含4个Mg和4个O，采用金属Mg与O₂分子作为参考，则晶胞形成能为E(Mg₄O₄)－4μMg－2E(O₂)。其中μMg是金属Mg的每原子能量。计算结果除以4得到每化学式形成能，除以8才是每原子形成能。

　　2、优化结构并提取静态能量

　　①分别准备目标材料和参考体系的POSCAR、INCAR、POTCAR、KPOINTS，检查元素顺序与赝势对应关系。

　　②完成结构优化，确认电子与离子收敛；将CONTCAR作为后续静态计算的结构，设置NSW为0。

　　③在静态计算中使用经过测试的截断能、k点密度和电子收敛条件，检查任务是否正常结束。

　　④从最终收敛结果提取能量，记录模型原子数和取值字段，再代入形成能表达式。

　　对于接近零温的固体能量比较，使用展宽方法时可关注OUTCAR中的energy(sigma->0)，并测试SIGMA与k点收敛。TOTEN包含展宽对应的电子熵项，不能不看设置就与其他能量混用。孤立原子或分子的占据情况还需单独检查，确认非物理展宽熵已足够小。

　　二、VASP形成能计算需要哪些参考能量

　　普通化合物的元素形成能，至少需要目标材料总能量及所有组成元素的参考能量。具体选什么结构，取决于所定义的参考状态。

　　1、固态元素使用对应单质相

　　金属元素通常采用所选温压条件下的稳定单质结构，其他固态元素也要确认同素异形体。例如碳采用石墨还是金刚石，会改变参考零点，报告时应明确说明。

　　①建立相应单质结构，检查晶相和磁性，优化后完成静态计算。

　　②将晶胞总能量除以其中的原子数，得到该元素的每原子参考能量。

　　单质晶胞含几个原子并不要求与化合物相同。需要统一的是换算后的计量关系，以及能量计算方法的可比性。

　　2、气态元素使用分子参考

　　①对以O₂、N₂或H₂为参考的体系，分别建立分子模型，设置足够的真空间隔。

　　②优化分子键长，并检查自旋设置；O₂应考虑三重态，不能直接按非磁性分子处理。

　　③增大盒子验证周期镜像影响，孤立分子通常使用Γ点计算。

　　④将分子总能量除以所含原子数，例如氧的每原子参考能量为E(O₂)/2。

　　采用孤立O原子代替O₂，会改变反应定义，不能继续把结果视为同一个氧气参考下的形成能。以自由原子为参考的能量差，通常用于讨论内聚或结合过程，还需明确符号约定。

　　3、有限温压条件需要化学势

　　零温静态计算得到的O₂能量，并不等于任意温度和氧分压下的氧化学势。若研究富氧、贫氧环境或合成窗口，需要进一步引入相应的温压贡献。

　　①确定温度、气体分压及热力学数据来源，按所用模型处理零点能、热修正和熵贡献。

　　②用目标相与竞争相的稳定性条件限制各元素化学势，明确参考零点与偏移量。

　　4、数据库参考能量需要核对计算方法

　　数据库中的总能量、修正后能量和形成能属于不同数据项。尤其是GGA与GGA+U混合结果，往往经过专门的兼容性处理。

　　①核对泛函、赝势版本、U值、磁性与能量修正方案。

　　②若无法建立兼容关系，重新计算所需参考体系；采用数据库方案时，应按其要求处理整组能量。

　　三、VASP形成能结果如何检查

　　形成能数值异常时，先回到计量与参考态，再检查计算精度。盲目增加计算步数，通常解决不了参考零点选错的问题。

　　1、复核归一化和收敛误差

　　①对照POSCAR重新统计原子数，检查分子能量是否已经换算为每原子值，以及是否重复除以化学式数量。

　　②逐步提高ENCUT、加密k点，比较形成能变化，确定精度是否满足研究要求。

　　③对可能有磁性的体系测试合理磁序，核对单质和化合物是否收敛到适当状态。

　　各体系可以使用不同的k点网格，但需要达到相近的收敛精度。含局域电子的材料还应检查U值及参考方案，不能默认给所有体系套同一个U值就能消除误差。

　　2、判断形成能能说明什么

　　形成能为负，只表示相对所选元素参考态更有利。材料仍可能分解为其他化合物，需要收集竞争相能量，构建相图或比较分解反应。与实验生成焓比较时，还应核对温度、压力、热修正及分子参考的系统误差，不能把数值差异全部归因于结构优化不足。

　　总结

　　VASP形成能计算需要一组能够相互比较的总能量。目标结构、元素单质和气体分子各自算对之后，还要检查原子数、参考态和单位是否对应。零温能量、有限温压化学势与数据库修正值也应区分清楚，否则公式代入无误，结果仍可能失去原本的物理含义。希望这些方法能帮助大家把参考能量的来源交代完整，在分析材料稳定性时更有依据。如果您还有相关需求，也欢迎联系我们进一步交流咨询。