第一性原理评估纤锌矿III-氮化物的热力学稳定性：交换关联泛函与自洽方法的比较研究

《Materials Today Communications》：First-principles evaluation of the thermodynamic stability of wurtzite III-nitrides: A comparative study of exchange-correlation functionals and self-consistent methods

　　

密度泛函理论（DFT）已成为预测材料稳定性和性质的重要工具。在本工作中，研究人员对纤锌矿相III-族氮化物AlN、GaN和InN的形成能进行了详细的第一性原理研究。研究人员评估了计算方法和交换关联泛函对这两种化合物能量稳定性的影响。采用了两种不同的方法求解Kohn-Sham方程：全势线性缀加平面波（FP-LAPW）方法（在WIEN2k和ELK代码中实现）和赝势+平面波（PP-PW）方法（在Quantum Espresso和VASP中实现），使用了模守恒赝势、超软赝势和投影缀加波赝势。考虑了几种交换关联（XC）泛函，包括基于GGA的泛函（PW91、PBE、PBEsol和WC）以及meta-GGA（SCAN）。结果表明，形成能受方法和XC泛函的选择显著影响。为验证InN相的动力学稳定性，计算了声子色散和振动态密度。还分析了热性质和拉曼活性声子模式，并与现有实验数据吻合良好。本研究为提高基于DFT的模拟的预测精度提供了有价值的基准，并为氮化物基光电子材料的设计提供了见解。

论文解读：第一性原理评估纤锌矿III-氮化物的热力学稳定性

研究背景与问题
III-族氮化物半导体（AlN、GaN、InN）因其可调的能带隙（从InN的0.7 eV到AlN的6.2 eV）在发光二极管、高电子迁移率晶体管和光电器件中具有广泛技术应用。然而，基于密度泛函理论（DFT）的第一性原理计算在预测这些材料的热力学稳定性时存在显著挑战。特别是，标准广义梯度近似（GGA）泛函如PBE对InN的纤锌矿（WZ）相形成能常给出正值，暗示其热力学不稳定，这与实验事实（WZ-InN为稳定相）明显矛盾。现有数据库（如Materials Project、Aflowlib、OQMD）中InN的形成能数值差异高达0.5 eV/原子，表明方法依赖性误差尚未被系统解耦。现有研究多采用单一计算实现，缺乏对两种根本不同数值方法（PP-PW与FP-LAPW）和多种交换关联（XC）泛函的统一比较。为此，研究人员系统评估了计算方法与XC泛函对WZ-III氮化物（AlN、GaN、InN）形成能的影响，并验证InN的动力学稳定性，旨在建立可靠的计算基准，指导氮化物光电器件的设计。该论文发表在《Materials Today Communications》。

研究内容与结论
研究人员利用赝势+平面波（PP-PW）方法（Quantum Espresso和VASP代码）和全势线性缀加平面波（FP-LAPW）方法（WIEN2k和ELK代码），结合多种GGA泛函（PW91、PBE、PBEsol、WC）和meta-GGA泛函（SCAN），对WZ-XN（X=Al, Ga, In）化合物的形成能进行了全面计算。结论表明：FP-LAPW方法（特别是WIEN2k）给出的形成能最接近实验值，正确预测了InN的稳定性；PBE泛函结合超软赝势（US-PP）高估形成能，对InN给出定性错误的正值；PBEsol和SCAN泛函显著改善精度。通过声子色散和振动态密度计算确认了WZ-InN的动力学稳定性，拉曼活性模式与实验吻合，热性质（热容、自由能、熵）符合预期趋势。本研究为改进DFT模拟的预测精度提供了关键基准，并揭示了方法选择对重离子体系的重要性。

主要技术方法
研究人员采用第一性原理DFT框架，分别使用Quantum Espresso（QE）和VASP进行PP-PW计算（采用模守恒NC、超软US和投影缀加波PAW赝势），以及WIEN2k和ELK进行FP-LAPW计算（设置R_MT为In/Al/Ga 2.30 Bohr、N 1.90 Bohr，R_MT×K_max=9）。所有XC泛函（PW91、PBE、PBEsol、WC、SCAN）在统一条件下实施。通过总能量随晶格常数变化确定平衡晶格常数，按公式E_F = (E_total - Σn_iμ_i)计算形成能（参考态为金属Al、Ga、In和N₂分子）。声子计算采用有限位移法（Phonopy+QE）和密度泛函微扰理论（DFPT），使用2×2×2超胞，验证WZ-InN的动力学稳定性。

研究结果
平衡晶格常数分析：PBE和PW91泛函预测InN的晶格常数比实验高约1.5%（如PBE-US给出3.595 ?，实验范围3.530–3.545 ?），而PBEsol和WC预测值落入实验范围。GaN和AlN的偏差类似，PBE呈现最大偏离（GaN偏大1.19%）。对于GaN，所有NC赝势组合在±3%范围内未找到能量最小值（标记*）。该结果通过总能量-晶格常数曲线（图2）和表1得出。

形成能分析：形成能随AlN→GaN→InN依次降低（变负），符合电负性趋势。PBE泛函在所有代码（QE、VASP）下均给出最负值（对InN甚至为正），与实验矛盾。例如，PBE-US在QE下对InN形成能为0.394 eV/原子（正数），而FP-LAPW（WIEN2k）的PBE结果均为负值，正确预测稳定性。PBEsol和SCAN泛函结果与实验吻合较好，但VASP和ELK对GaN仍低估实验值。WIEN2k的FP-LAPW方法在所有三个化合物中均给出实验范围内的形成能。图5进一步显示，PW91和PBE在QE的各赝势下均高估形成能，PBEsol改善但仍偏高，仅PBEsol-NC对InN落在实验误差带。这些差异源于In的4d半芯态与N的2s态之间的闭壳层排斥，FP-LAPW通过全电子处理避免了赝势的芯-价分割误差。该结果通过形成能二维色图（图3）和对比图（图4、5）得出。

声子色散与动力学稳定性：PBE-US下的有限位移法和DFPT均给出无虚频声子谱，确认WZ-InN动力学稳定。DFPT/PBE结果在Γ点附近出现数值假象（U形异常），而DFPT/PBEsol和Phonopy结果正常，表明该假象源于PBE功能限制和自旋-轨道耦合（SOC）缺失。声子态密度显示低频声学模主要由In贡献，高频光学模由N贡献。该结果通过图6得出。

拉曼活性模式与热性质：计算拉曼光谱显示五个峰（E₂(low) 82.43 cm^-1、A₁(TO) 419.91 cm^-1、E₁(TO) 439.14 cm^-1、E₂(high) 444.0 cm^-1、A₁(LO) 540.64 cm^-1），与实验相比低估5–9%，归因于PBE和未含SOC。热性质计算（图8）显示振动自由能随温度单调下降，熵逐步趋近饱和，恒容热容C_V遵循低温T³定律并趋近Dulong-Petit极限。该结果通过图7、8和表2得出。

总结讨论与结论翻译
讨论部分指出，精确形成能计算对设计异质结构（如HEMT、LED）至关重要，尤其InN基体系的热力学混溶极限和相分离预测需要可靠基准。FP-LAPW方法结合PBEsol或SCAN泛函优于PBE+US。高温设备的热性质（热容、自由能）对散热和结构完整性至关重要。本研究局限性包括未含SOC和混合泛函（HSE06），未来需引入准谐或分子动力学方法。结论翻译如下：本工作对纤锌矿相III-族氮化物（AlN、GaN、InN）的热力学和振动性质进行了全面的第一性原理研究。通过结合PP-PW和FP-LAPW方法及多种DFT代码与泛函（GGA和meta-GGA），系统评估了方法选择对预测形成能的影响。结果表明，WIEN2k代码实现的FP-LAPW方法给出的形成能与实验值更接近，优于PP-PW方法。这强化了WIEN2k作为验证其他DFT代码结果可靠基准的作用。值得注意的是，对于InN，泛函选择至关重要；PBEsol在标准PBE基础上提供了改进精度，但PBE特别是与超软赝势结合时，可能产生定性错误的不稳定性趋势。这些差异凸显了半局域形式对重离子体系的固有局限。为验证结果，研究人员计算了InN的声子色散和振动态密度，确认其动力学稳定性。分析的拉曼活性模式与实验光谱高度一致，进一步证明了所选计算框架的可靠性。这些发现通过识别筛选新型合金和异质结构的可靠模拟方法，指导了材料设计。准确评估合成可行性对于优化光电器件生长条件和减少缺陷至关重要。尽管本研究提供了稳健框架，但受限于不含自旋-轨道耦合和非谐效应。未来需纳入准谐或分子动力学方法以捕捉温度依赖的声子重整化和寿命，实现与实验更接近的比较。