Systematic comparison of approximations and functionals in first-principle calculations of aluminum-based III-V ferroelectric nitrides

本研究系统评估了化学无序建模（VCA 与 SQS）及交换关联泛函（PBE、PBESol、SCAN、SCAN+rVV10）对铝基 III-V 族氮化物结构与铁电性质的影响，结果表明，SQS 方法结合 SCAN 泛函为预测这些材料的相稳定性并识别亚稳态提供了最可靠的框架。

要点

想象你是一位建筑师，正试图为一种新型超快、可存储记忆的计算机芯片设计完美的建筑材料。你拥有两种主要原料：氮化铝（一种坚固可靠的“砖块”）和第二种可混合以改变其性质的成分。你可以选择混合钪（一种重金属元素）或硼（一种微小、轻质的元素）。

目标是创造出一种能像“铁电”开关那样运作的材料——一种能够通过翻转其内部电方向来记住自身是“开”还是“关”的材料。然而，准确预测这些混合材料的行为，就像试图在混乱的风暴中预测天气一样。你需要一个计算机模型来模拟原子，但该模型本身的准确性取决于你如何设置它。

本文本质上是对不同计算机模型进行的一次大规模“压力测试”，旨在看看哪一个模型能揭示这些基于铝的氮化物材料的真相。

作者研究的两个主要问题

作者发现，获得正确答案取决于解决两个特定的谜题：

1. “拥挤的房间”与“平均人”问题（无序性）
当你将铝与钪或硼混合时，原子并不会像列队士兵那样以完美、重复的模式排列。它们杂乱无章、随机分布，就像一场拥挤的派对，每个人都在争抢空间。

• 旧方法（虚拟晶体近似）： 想象一下，通过说“平均身高是 5 英尺 9 英寸，穿着蓝色衬衫”来描述这场派对。这就是虚拟晶体近似（VCA）。它平滑了混乱。本文表明，这种方法是个糟糕的骗子；它会让材料看起来稳定，而实际上它并不稳定，反之亦然。这就像说一座由沙子和水建成的房子是坚固的，因为沙子和水的“平均”是“泥浆”。
• 新方法（特殊准随机结构）： 这就像拍摄派对的实际混乱场景，特定的人站在特定的位置。这就是特殊准随机结构（SQS）。作者发现，要得到正确答案，你必须观察原子具体、杂乱的排列方式，而不仅仅是平均值。

2. “透镜”问题（泛函）
即使你拥有了正确的杂乱排列，你仍然需要通过一个特定的数学“透镜”（称为交换 - 关联泛函）来观察它，以计算能量。作者测试了四种不同的“透镜”：PBE、PBESol、SCAN 和 SCAN+rVV10。

• 结果： 某些“透镜”（如 PBESol）模糊不清，扭曲了图像，使材料过早地显得不稳定。而其他“透镜”（如SCAN）则像高清眼镜一样，展示了材料真实的稳定性。

关于两种混合物的发现

本文揭示，掺入钪和掺入硼就像是两个完全不同的故事，尽管它们始于相同的基底材料。

故事 A：掺入钪（重金属）

• 行为： 当你加入钪时，原子倾向于彼此靠得更近。它们开始偏好一种“拥挤”的排列（称为岩盐相），而不是承载记忆开关的“宽敞”排列（纤锌矿相）。
• 惊喜： “模糊”的模型（VCA）预测这种转变会在很低的钪浓度下很快发生。但“高清”模型（SQS + SCAN）显示，该材料能保持更长时间的稳定性和实用性——直到接近 50% 的钪浓度。这与现实世界的实验观察相符。
• 转折： 存在一种奇怪的中间状态（一种五边形六方相），它充当了踏脚石。这是一个亚稳态的“休息站”，原子在最终定居到拥挤状态之前会造访此处。

故事 B：掺入硼（微小元素）

• 行为： 硼非常微小，偏好以扁平的三角形（三边）形状存在，而不是三维金字塔形状。当你加入硼时，它会迫使结构断裂并重新配置。
• 断裂： 在中等量的硼存在下，原子间的化学键实际上会断裂并重新排列。材料发生扭曲，而“记忆开关”（极化）最初实际上会增强，这是一件好事。
• 终局： 如果你加入过多的硼，材料会完全放弃三维金字塔形状，转变为扁平的层状结构（像石墨或一叠纸）。这是性格的彻底转变。

最终裁决：“黄金标准”

在测试了所有“杂乱房间”模型和“透镜”的组合后，作者得出结论，预测这些材料行为最佳的方法是结合使用：

1. SQS：以捕捉原子真实的、杂乱的随机性。
2. SCAN：使用最准确的数学“透镜”。

这为何重要？
本文并不声称要今天建造一个新的计算机芯片。相反，它提供了蓝图的蓝图。它告诉科学家：“如果你想设计一种新的铁电材料，不要使用旧的、简单的数学工具。使用这种特定的、更复杂的工具组合，否则你的预测将是错误的。”

通过使用正确的工具，他们确认了钪混合物非常稳定，对存储设备充满希望；而硼混合物则很棘手——它们可以提升性能，但前提是你必须在结构坍塌成扁平层状物之前停止添加。

简而言之：不要相信平均值；要观察混乱。不要使用模糊的透镜；要使用高清的那一个。

技术摘要

技术摘要：基于第一性原理计算的铝基 III-V 族铁电氮化物中近似方法与泛函的系统性比较

问题陈述
III-V 族氮化物半导体（特别是氮化铝钪 Al$_{1-x}$Sc$_x$N 和氮化铝硼 Al$_{1-x}$B$_x$N）中铁电性的发现，已将研究焦点从传统的钙钛矿氧化物转移开来。尽管基于密度泛函理论（DFT）的第一性原理计算在理解这些材料方面处于核心地位，但显著的方法论不确定性依然存在。先前的研究主要依赖虚拟晶体近似（VCA）来模拟化学无序，并经常使用特定的交换关联泛函（如 PBE 或 PBESol），而缺乏系统的基准测试。目前缺乏全面的研究来考察无序处理方法的选择（VCA 与特殊准随机结构 SQS）以及交换关联泛函的选取（PBE、PBESol、SCAN、SCAN+rVV10）如何定量和定性地影响这些无序三元合金预测的结构、能量和铁电性质。具体而言，尚不清楚像 SCAN 这样的 meta-GGA 泛函是否为氮化物铁电体提供了更优的预测精度，以及像 VCA 这样计算成本较低的方法是否足以捕捉复杂的相稳定性景观。

方法论
作者采用基于 48 原子超胞的统一计算框架，比较了五种竞争多晶型：极性纤锌矿（WZ）、闪锌矿（ZB）、六方（HX）、六方层状（HL）和岩盐（RS）。这种超胞方法通过变换晶格矢量，使得所有相能够在同等基础上进行描述。

该研究系统地变化了两个关键的方法论参数：

1. 无序处理：作者在 Abinit 代码中实现的虚拟晶体近似（VCA）与使用合金理论自动化工具包（ATAT）生成并在维也纳从头算模拟包（VASP）中计算的“特殊准随机结构”（SQS）进行了对比。
2. 交换关联泛函：该研究基准测试了四种泛函：PBE、PBESol、SCAN（一种 meta-GGA）以及包含范德华相互作用的 SCAN+rVV10。

计算采用投影缀加波（PAW）赝势，平面波截断能为 620 eV。作者分析了不同 Sc 和 B 浓度（$x$）下的总 Kohn-Sham 能量、晶格常数、内禀参数（$u$）、键分布、自发极化（通过贝里相位）、玻恩有效电荷以及电子带隙。

主要贡献与结果

• 无序处理的影响（VCA 与 SQS 对比）：

  • Al$_{1-x}$Sc$_x$N：VCA 严重低估了铁电 WZ 相的稳定性域。VCA 预测在约 19% Sc 含量时发生从 WZ 到非极性岩盐（RS）相的转变，而 SQS 计算预测该转变发生在接近 47.8% Sc 处。SQS 结果与高达 43% Sc 仍保持铁电性的实验报道更为一致。
  • Al$_{1-x}$B$_x$N：VCA 未能捕捉到由硼倾向于三配位引起的强烈结构畸变。SQS 揭示，在中等硼浓度（20–50%）下，会发生断键，导致畸变的 WZ 和 ZB 结构，这些结构在动力学上是稳定的，但与常规相不同。VCA 无法重现这些局部无序效应。

• 交换关联泛函的影响：

  • PBE 与 PBESol：尽管两者相似，但这些泛函得出了截然不同的结果。PBESol 显著低估了 WZ 到 RS 转变的临界 Sc 浓度（预测约为 29.1%，而 PBE 预测约为 47.8%），并错误地预测 ZB 相为高硼浓度下的基态，而 PBE 预测的是六方层状（HL）相。
  • PBE 与 SCAN：SCAN meta-GGA 泛函在大多数性质上显示出与 PBE 出色的定性和定量一致性，但提供了更一致的相能量描述。SCAN 预测 Al$_{1-x}$Sc$_x$N 中的 WZ 到 RS 转变发生在 $x \approx 0.434$，这比 PBESol 更接近实验观测值。
  • 范德华效应：在 SCAN 中包含 rVV10 修正（SCAN+rVV10）将 RS 相的稳定性移至更低的 Sc 浓度（$x \approx 0.359$），并减小了 Al$_{1-x}$B$_x$N 中 HL 相的稳定性窗口。然而，鉴于计算开销以及 SCAN 单独已能提供与实验的稳健一致性，作者建议 SCAN 是这些系统的优选泛函。

• 结构与铁电性质：

  • Al$_{1-x}$Sc$_x$N：Sc 含量增加导致 $c/a$ 比降低，自发极化减小。该系统表现出复杂的能量景观，其中 WZ、HX 和 RS 相在转变点附近近乎简并。HX 相（五配位）被确定为 Farrer 和 Bellaiche（2002）预测的低能亚稳态，出现在 WZ 和 RS 相之间。
  • Al$_{1-x}$B$_x$N：与 Sc 掺杂相反，B 掺杂导致 WZ 相迅速失稳。在中等浓度下，B 原子倾向于三配位，导致断键和结构畸变。这导致在低硼浓度下自发极化得到增强，随后结构转变为 ZB 或 HL 相。
  • 电子性质：AlN 中的 Sc 掺杂显著降低了带隙，特别是在 RS 相中，这是由于 Sc $d$ 态与 N $p$ 态的杂化引入了部分莫特（Mott）特征。相反，AlN 中的 B 掺杂通常增加带隙，但由于键的重构，其演变是非单调的。

意义与主张
本文主张，特殊准随机结构（SQS）与 SCAN meta-GGA 泛函的结合，为预测新兴氮化物基铁电材料的性质提供了最一致和可靠的理论框架。作者强调：

1. 局部无序至关重要：VCA 不足以准确模拟这些合金中的稳定性窗口和结构畸变，特别是在 Al$_{1-x}$B$_x$N 中，局部断键效应显著。
2. 泛函选择很重要：即使在相似的泛函之间（PBE 与 PBESol），定量差异也是巨大的。与杂化泛函相比，SCAN 提供了准确性与计算成本之间的更优平衡，在描述这些铁电体的结构和能量性质方面优于 GGA 泛函。
3. 独特的物理机制：该研究确立了由配位倾向驱动清晰的“结构 - 性质”关系：Sc 促进更高的配位（有利于 RS/HX 相），而 B 促进更低的配位（有利于 ZB/HL 相）。

这项工作并未提出新的实验应用，而是为下一代氮化物铁电体的准确设计与优化提供了一条方法论路线图和必要的物理洞察。它阐明了先前理论与实验之间的差异可能源于使用了不恰当的近似（VCA）或泛函（PBESol），而非底层物理的失效。