手机版

形成能 DFT 计算

发布时间：2026-08-31 来源：科研学术网

小中大

字号：

一、项目背景与科学问题

形成能（formation energy）是材料热力学稳定性的核心判据。对化合物，形成焓定义为单位化学式从单质生成化合物时的能量变化，其正负与大小决定化合物能否稳定存在及分解倾向；对缺陷（空位、间隙、掺杂原子），缺陷形成能决定缺陷的平衡浓度与类型，直接影响半导体的载流子浓度、发光性质与催化活性。第一性原理 DFT 计算可在原子层面预测形成能，为新材料合成与缺陷工程提供热力学依据。

本项目包含两部分：一是计算 SiC、GaN、AlN、ZnO 四种半导体的形成焓并与实验对比，评估 GGA-PBE 泛函的精度；二是计算 ZnO 中氧空位的形成能随氧化学势的变化，分析缺陷浓度与生长条件的关系。

二、计算方法与理论背景

2.1 化合物形成焓

化合物 AₓBᵧ 的形成焓定义为 ΔH_f = E(AₓBᵧ) − xE(A_bulk) − yE(B_bulk)，其中 E(AₓBᵧ) 为化合物晶胞总能，E(A_bulk) 与 E(B_bulk) 为参考单质（稳定相）的每原子能量。ΔH_f 为负表示形成放热、化合物热力学稳定；绝对值越大结合越强。参考单质需取对应标准状态：Si、Ga、Al 取各自晶体，N 取 N₂ 分子，Zn 取金属，O 取 O₂ 分子。

2.2 缺陷形成能

缺陷形成能定义为 E_f = E_def − E_host + Σ nᵢμᵢ，其中 E_def 为含缺陷超胞总能，E_host 为完美超胞总能，nᵢ 为移除/添加的原子数，μᵢ 为相应元素的化学势。化学势的取值反映生长环境：富氧（O-rich）条件下 μ_O 高、氧空位形成能低；贫氧（Zn-rich）条件下 μ_O 低、氧空位形成能高。因此缺陷形成能是环境依赖的量，其范围可由单质化学势上下限约束。

2.3 超胞与收敛

缺陷计算需采用足够大的超胞以减弱缺陷-镜像相互作用，通常 3×3×2 以上；原子位置需完全弛豫（缺陷引起局域晶格畸变），能量收敛至 1 meV/原子量级。k 点与截断能需与完美超胞保持一致以保证误差抵消。

三、计算设置与模型构建

软件：VASP；泛函：GGA-PBE，PAW 赝势，截断能 500 eV。
化合物：SiC（闪锌矿）、GaN（纤锌矿）、AlN（纤锌矿）、ZnO（纤锌矿）原胞，全优化晶格与原子。
参考态：Si/Ge/Al/Zn 金属晶体、N₂/O₂ 分子（大盒子）。
缺陷：ZnO 3×3×2 超胞，移除一个 O 原子形成中性氧空位，全弛豫。
化学势扫描：μ_O 在 Zn-rich 与 O-rich 极限间变化。

四、结果分析与案例图解读

图 6 的左面板给出四种半导体形成焓的 DFT（蓝色柱）与实验（橙色柱）对比：SiC −0.62/−0.68 eV、GaN −1.04/−1.13 eV、AlN −2.94/−3.13 eV、ZnO −3.21/−3.63 eV（每化学式单位）。所有 DFT 值均为负且量级正确，说明 PBE 准确预测了化合物的热力学稳定性；同时 DFT 值系统地偏正（绝对值偏小）0.06–0.42 eV，相对偏差约 5%–12%，体现 GGA 对强离子/共价键合体系的轻微欠束缚。值得注意的是，形成焓绝对值随离子性增强而增大（SiC < GaN < AlN < ZnO），DFT 对这一趋势的复现完全正确，验证了方法对化学键类型的敏感度。

图 6 的右面板给出 ZnO 中中性氧空位的形成能随氧化学势的变化（红色曲线）。横轴 μ_O 从贫氧（Zn-rich，−5 eV）到富氧（O-rich，−1 eV）变化，形成能 E_f 随 μ_O 升高而线性下降：在 Zn-rich 极限约 4.5 eV，在 O-rich 极限约 1.2 eV。这一趋势的物理意义明确：富氧条件下氧原子容易从环境中获得，移除氧原子形成空位所需能量低，因此氧空位浓度高；反之贫氧条件下氧空位稀少。E_f 变化范围约 3 eV，说明生长气氛对 ZnO 缺陷化学有巨大影响。由 E_f 可进一步估算平衡空位浓度 n ∝ exp(−E_f/k_BT)，为实验（如 PL 光谱的绿光发射与氧空位关联）提供定量预测。

两面板结合，左面板验证了 DFT 对化合物热力学稳定性的预测能力，右面板揭示了缺陷形成能的环境依赖性，共同构成”从晶体形成到缺陷化学”的完整热力学图景，为材料合成气氛控制与缺陷工程提供定量依据。