形成能DFT计算在晶体结构筛选中的实战应用

形成能DFT计算看起来是个简单活——优化结构、算总能、减一下参考态能量就行——但实际做起来，坑的深度和广度都容易被低估。一个高熵合金项目，目标是从十五个候选组成里筛选出热力学最稳定的三到四个，用形成能是否足够负作为第一道筛子。项目开始时以为两周搞定，实际上花了五周才拿到可信结果。

参考态选择：最容易踩的坑

形成能定义里最容易被忽略的细节是参考态的选择。对于二元合金，参考态通常用纯元素的同素异形体。但问题来了：如果计算纯元素时用的是错误的晶格结构，形成能会系统性偏正，项目可能把实际上稳定的组分误判为不稳定。这个项目的第一轮计算就出了这个问题：铌的参考态用了错误的晶格结构，导致含铌组分的形成能DFT计算结果整体偏正0.1到0.2 eV每原子，三个实际上稳定的组分在第一轮被错误地筛掉了。

补救方法说起来简单：把每个元素参考态的所有可能晶体结构都算一遍，取能量最低的那个作为参考态。实际做起来，十个元素每个算2到3个结构，就是二十多个单点能计算。项目在这里多花了十天，但这一步骤省掉的后续麻烦远超这个时间成本。

振动修正：该不该加

形成能的温度依赖性来自振动自由能贡献。对于高温稳定的化合物，振动修正可能把形成能从负值推到正值，改变相图上的相边界。这个项目在900K工作温度下运行，振动修正要不要加，组里两派意见分歧很大。最终的决定是：对形成能绝对值最大的三个候选组成做振动修正，其他组成暂时用0K形成能筛选。

结果出来后，振动修正对形成能的影响在0.03到0.08 eV每原子之间，对于形成能是否小于负0.1 eV每原子这个判断阈值来说，影响不算致命。但有一组成分的形成能在0K是负0.09 eV每原子，加振动修正后变成负0.02 eV每原子，从勉强稳定变成实际上不稳定——这个案例后来被写进了组里的计算规范。

形成能DFT计算的价值不在于给出一个绝对值，而在于给相竞争格局一个可比较的量化标尺。参考态选择、振动修正、磁性基态——每一个细节都可能改变排序结果。