MS计算形成能：Materials Studio形成能计算详解

形成能的物理意义

形成能（Formation Energy）是衡量材料热力学稳定性的关键物理量，表示从组成元素的稳定单质状态形成化合物所需的能量变化。形成能越负，说明化合物相对于单质越稳定，越有可能被实验合成。

在材料设计中，形成能常用于：

• 筛选稳定的新材料候选
• 评估掺杂和缺陷的形成难易程度
• 比较不同晶相的相对稳定性
• 构建相图和稳定性相图

形成能计算公式

化合物形成能

对于化合物 A_xB_y，形成能定义为：

ΔH_f = E(A_xB_y) – x·E(A) – y·E(B)

其中：

• E(A_xB_y) 为化合物的总能量（每原子或每公式单位）
• E(A)、E(B) 为对应元素在稳定单质状态下的总能量
• x、y 为化学计量比

形成能判断标准：

• ΔH_f < 0：化合物比单质稳定，可能存在
• ΔH_f > 0：化合物不稳定，可能分解为单质
• |ΔH_f| 越大，热力学稳定性越高

缺陷形成能

对于含缺陷的体系，缺陷形成能定义为：

ΔH_f(defect) = E(defect) – E(perfect) – Σ n_i·μ_i + q·E_F

其中：

• E(defect)：含缺陷超胞的总能量
• E(perfect)：完美超胞的总能量
• n_i：加入(i<0)或移除(i>0)的原子数
• μ_i：对应原子的化学势
• q：缺陷电荷态
• E_F：费米能级（相对于VBM）

MS/CASTEP形成能计算步骤

1. 计算化合物总能量

CASTEP参数设置：

任务：Energy（单点能计算）
泛函：GGA-PBE
赝势：OTFG ultrasoft
截断能：500 eV（需做收敛测试）
k点：自动生成，间距0.04 Å⁻¹
SCF收敛：1.0×10⁻⁶ eV/atom

确保结构已充分优化后再做单点能计算。

2. 计算单质总能量

需要计算化合物中每个元素在稳定单质状态下的能量：

常见单质参考结构：

注意事项：

• 单质计算必须使用与化合物相同的泛函、截断能和赝势
• 对于气体分子（O₂, N₂, H₂），需要使用大盒子（≥10 Å）和Γ点
• 对于磁性单质（Fe, Co, Ni），必须开启自旋极化

3. 单质计算的具体设置

固体单质（以Fe为例）：

结构：BCC Fe（a=2.87 Å）
自旋极化：开启（Fe磁性）
初始磁矩：2.2 μB
k点：12×12×12

气体分子（以O₂为例）：

结构：O₂分子放在10×10×10 Å盒子中
自旋极化：开启（O₂三重态）
初始磁矩：2.0 μB
k点：Γ点（1×1×1）

4. 形成能计算

收集所有能量数据后，代入公式计算：

示例：TiO₂形成能计算

形成能 = E(TiO₂) – E(Ti) – E(O₂) = -99.68 – (-7.83) – (-8.70) = -83.15 eV/f.u. = -9.42 eV/atom

实验值约 -9.40 eV/atom，计算精度良好。

不同类型形成能计算

1. 二元化合物

AB型化合物： ΔH_f = E(AB) – E(A) – E(B)

A_xB_y型化合物： ΔH_f = [E(A_xB_y) – x·E(A) – y·E(B)] / (x+y)

2. 三元化合物

ABO₃型钙钛矿： ΔH_f = E(ABO₃) – E(A) – E(B) – 3/2·E(O₂)

注意O的能量使用O₂分子的能量除以2得到单个O原子的参考能量。

3. 固溶体形成能

对于合金或固溶体 A_{1-x}B_x：

ΔH_mix = E(A_{1-x}B_x) – (1-x)·E(A) – x·E(B)

需要构建特殊准随机结构（SQS）来模拟无序固溶体。

4. 表面形成能

表面能计算： γ = [E(slab) – n·E(bulk)] / (2·A)

其中：

• E(slab)：slab模型总能量
• n·E(bulk)：等量体相原子能量
• 2A：上下两个表面的总面积
• 除以2因为slab有两个表面

形成能相图

凸包图（Convex Hull）

凸包图是分析多组分体系稳定性的重要工具：

1. 计算不同组分比例下所有可能化合物的形成能
2. 在组分-形成能平面上绘制数据点
3. 连接最稳定的点形成凸包（convex hull）
4. 凸包上的点：稳定相
5. 凸包上方的点：不稳定相（会分解为凸包上的相邻相）

凸包上方能量（energy above hull）： ΔE_hull = E(compound) – E_hull(同一组分)

ΔE_hull < 0.05 eV/atom 通常被认为可能稳定（考虑计算误差）。

缺陷形成能计算

化学势范围确定

缺陷形成能中的化学势μ_i需要在合理范围内取值：

1. 上限：单质化学势（富集条件）
2. 下限：由化合物形成能约束
3. 约束条件：Σ n_i·μ_i = ΔH_f(compound)

示例：TiO₂中化学势范围

Ti富集条件：μ_Ti = μ_Ti(bulk), μ_O = μ_O(O₂) + ΔH_f(TiO₂) O富集条件：μ_O = μ_O(O₂), μ_Ti = μ_Ti(bulk) + ΔH_f(TiO₂)

费米能级范围

对于带电缺陷，费米能级E_F在带隙内变化：

• VBM处：E_F = 0
• CBM处：E_F = E_gap

不同电荷态的缺陷形成能随费米能级线性变化，斜率为q。

常见问题与解决方案

问题1：形成能与实验值偏差大

可能原因：

1. 截断能不足：增大ENCUT到500 eV以上
2. k点不够密：减小k点间距到0.03 Å⁻¹
3. 结构未充分优化：先做高质量结构优化
4. 泛函选择：PBE可能不够精确，考虑HSE06
5. 范德华力：对分子晶体或层状材料加DFT-D校正

问题2：气体分子计算收敛困难

O₂、N₂等开壳层分子的计算容易出问题：

解决方案：

# CASTEP设置
Spin polarized: Yes
Initial spin: 2.0（O₂三重态）
# 使用更大盒子（15 Å）
# 降低smearing值

问题3：含磁性的形成能计算

磁性化合物的形成能计算需要特别注意：

1. 化合物和单质都必须正确设置磁性
2. 比较不同磁结构的形成能
3. 注意磁矩对总能量的影响
4. DFT+U计算中U值需一致

问题4：形成能为正是否一定不稳定

形成能略正（<0.1 eV/atom）的化合物可能在特定条件下亚稳存在：

• 高压条件下可能稳定
• 动力学稳定（声子谱无虚频）但热力学不稳定
• 需要结合动力学稳定性判断

计算精度优化建议

收敛测试清单

泛函选择建议

实操案例：SrTiO₃形成能

以SrTiO₃钙钛矿为例：

1. 计算组分能量：
   • SrTiO₃：PBE/500eV，结构优化后E = -89.45 eV
   • Sr(FCC)：E = -1.72 eV/atom
   • Ti(HCP)：E = -7.83 eV/atom
   • O₂分子：E = -8.70 eV
2. 形成能计算： ΔH_f = -89.45 – (-1.72) – (-7.83) – 1.5×(-8.70) = -89.45 + 1.72 + 7.83 + 13.05 = -16.85 eV/f.u. = -6.72 eV/atom
3. 对比实验： 实验值约 -16.70 eV/f.u.，偏差 < 1%。

总结

形成能计算是材料热力学稳定性分析的核心方法。通过MS/CASTEP进行高质量的第一性原理计算，可以准确预测化合物的形成能，为新材料设计提供可靠的理论指导。在实际计算中，需要注意单质参考态的正确选择、计算参数的一致性以及磁性体系的特殊处理，以确保形成能计算的准确性和可比性。

我们提供专业的Materials Studio形成能计算服务，涵盖各类化合物、固溶体、缺陷体系的热力学稳定性分析，支持从参数测试到凸包图构建的完整计算流程。