Full-Potential Korringa Kohn Rostoker Study of Structural and Mechanical Properties in SiGeSn Alloys

技术摘要：全势 KKR 研究 Si$_{1-(x+y)}$Ge$_y$Sn$_x$ 合金的结构与力学性质

问题陈述
对于开放式共价键晶体结构（如金刚石结构的 IV 族半导体），精确的第一性原理描述在传统的 Korringa–Kohn–Rostoker (KKR) 方法中仍然具有挑战性。这些传统方法严重依赖于对有效晶格势的球对称近似，特别是 muffin-tin 近似 (MTA) 和原子球近似 (ASA)。虽然这些近似降低了计算复杂度，但往往无法捕捉到在具有方向性键合的开放结构中至关重要的非球对称势贡献。因此，传统的 KKR 计算可能会在总能量、原子间作用力和平衡结构参数方面产生显著误差，并可能错误地预测半导体的金属特性。

方法论
作者利用 MuST 软件包中新开发的实现的全势 KKR (KKR-FP) 方法来解决这些局限性。该研究采用多重散射理论 (MST) 形式体系，通过格林函数法求解 Kohn-Sham 方程，而非显式的哈密顿量对角化。

• 势处理： 研究系统地比较了三种有效势的处理方式：MTA、ASA 以及保留了非球对称贡献的全势 (KKR-FP) 方法。
• 计算细节： 计算在局部密度近似 (LDA) 下进行交换相关处理。单点散射问题通过角动量截断 $l_{max} = 4$ 进行求解。布里渊区积分使用 $12 \times 12 \times 12$ 的 Monkhorst–Pack k 点网格，对于态密度 (DOS) 计算则使用更密的网格。
• 研究体系： 该方法以金刚石结构的元素硅 (Si)、锗 (Ge) 和 $\alpha$-锡 (Sn) 进行基准测试。随后，该框架被应用于具有特定组分的三元 Si$_{1-(x+y)}$Ge$_y$Sn$_x$ 合金：Si$_{0.5}$Ge$_{0.375}$Sn$_{0.125}$、Si$_{0.375}$Ge$_{0.25}$Sn$_{0.375}$ 以及二元 Si$_{0.5}$Sn$_{0.5}$。
• 分析： 通过将计算的总能量拟合到三阶 Birch–Murnaghan 状态方程，确定平衡晶格参数 ($a_0$) 和体积模量 ($B_0$)。通过计算混合焓 ($\Delta H_{mix}$) 来评估相对于元素参考态的合金形成能量有利程度。

主要贡献与结果

1. Si 的势近似基准测试：

   • MTA 失效： Muffin-tin 近似未能产生 Si 的明确能量极小值，总能量随晶格参数单调下降，导致其对于这种开放结构而言是不可靠的。
   • ASA 局限性： 原子球近似虽然产生了确定的极小值，但显著低估了平衡晶格参数（$a_0 = 5.03$ Å，实验值为 $5.431$ Å，误差约为 7.4%）。此外，KKR-ASA 在费米能级处错误地预测了有限的态密度，暗示了金属行为。
   • KKR-FP 的成功： 全势方法准确地重现了实验平衡晶格参数（$a_0 = 5.433$ Å，偏差 $\approx 0.04\%$）和体积模量（$B_0 = 94.9$ GPa，实验值 $\approx 97.9$ GPa）。至关重要的是，KKR-FP 恢复了 Si 的半导体特性，在费米能级处表现出清晰的能隙。

2. 向 Ge 和 $\alpha$-Sn 的扩展：

   • KKR-FP 方法被应用于 Ge 和 $\alpha$-Sn。对于 Ge，计算得到的 $a_0$ 为 $5.73$Å（对比实验值$5.652$Å），$B_0$为$70.86$GPa（对比实验值$\approx 75$GPa）。对于$\alpha$-Sn，计算得到的$a_0$为$6.7$Å（对比实验值$6.483$ Å）。尽管存在一些偏差，但这些结果处于先前理论预测的范围内，并确立了该方法在 IV 族材料中的可迁移性。

3. SiGeSn 合金的结构与力学性质：

   • 晶格参数： 计算得到的各类三元合金平衡晶格参数表现出对 Vegard 定律的正偏差。对于 Si$_{0.375}$Ge$_{0.25}$Sn$_{0.375}$，计算得到的 $a_0$ ($6.11$ Å) 比 Vegard 定律的估计值 ($6.00$ Å) 高出约 1.83%。
   • 体积模量： 观察到随着 Sn 含量的增加和组分的变化，体积模量大幅下降。$B_0$ 从 Si$_{0.5}$Ge$_{0.375}$Sn$_{0.125}$ 的 $88.18$GPa 下降到 Si$_{0.375}$Ge$_{0.25}$Sn$_{0.375}$的$58.23$ GPa，表明压缩性增加。
   • 混合焓： 所有研究的组分均得到了正的混合焓，表明相对于选定的元素参考态，同质合金形成的能量存在惩罚。具体而言，Si$_{0.5}$Ge$_{0.375}$Sn$_{0.125}$ 的 $\Delta H_{mix}$ 为 $+180.92$ meV/atom，而 Si$_{0.375}$Ge$_{0.25}$Sn$_{0.375}$ 为 $+346.32$ meV/atom。二元 Si$_{0.5}$Sn$_{0.5}$ 合金的值为 $+237.87$ meV/atom。

意义与主张
本文声称，全势处理对于准确描述像 Si 这样具有开放结构的半导体的结构和电子性质是必不可少的。该研究确立了 MuST 软件包中的 KKK-FP 实现是研究成分复杂的 IV 族合金的一个有前景的框架。结果强调，非球对称势贡献对于恢复正确的平衡晶格参数、体积模量和半导体带隙至关重要。

关于 SiGeSn 合金，作者得出结论，正的混合焓表明相对于所选的元素参考态，存在向化学偏析或相分离发展的能量驱动力。然而，他们也审慎地指出，混合焓本身并不能决定相对于所有竞争相在有限温度下的热力学稳定性，因为在有限温度下，构型熵可能会有利于混合相的存在。本研究为理解这些合金的结构行为和混合热力学提供了第一性原理基础，为解释合金稳定性的实验观察提供了视角，但并未在缺乏进一步热力学分析的情况下预测具体的相分离间隙或相边界。