Limited Diffusion of Silicon in GaN: A DFT Study Supported by Experimental Evidence

本研究将第一性原理密度泛函理论计算与超高压退火实验相结合，证明由于激活能垒过高，硅在氮化镓中的扩散极为有限，从而证实了该材料在先进电子应用中实现精确掺杂的稳定性。

要点

想象一下，氮化镓（GaN） 是一座为电子未来而建的、高科技且超耐用的城市。它是驱动我们明亮 LED 灯和高速互联网连接的材料。为了让这座城市运转，工程师需要在特定的街区加入被称为硅（Si） 原子的“市民”。这些硅原子充当电流载体（施主），使器件得以开启。

研究人员提出的核心问题是：一旦我们将这些硅市民安置在它们的家中，它们是会安守原位，还是会四处游荡？

在许多材料中，原子就像不安分的游客；如果你加热它们，它们就开始收拾行李，搬往新的地方。这种“游荡”（扩散）对电子器件是有害的，因为它会模糊芯片不同部分之间精确的界限。研究团队想要知道，氮化镓中的硅是“宅男”还是“旅行者”。

以下是他们发现的内容，以简单的方式解释：

1. “空座”理论（原子如何移动）

要在晶体城市中从一个位置移动到另一个位置，原子通常需要旁边有一个空座（空位）供其跳入。

• 研究内容： 科学家们利用强大的计算机模拟（如同超精准的视频游戏），观察一个硅原子试图跳入空座的过程。
• 结果： 他们发现，硅原子为了完成那次跳跃所必须攀登的“台阶”高得惊人。
  • 横向移动（沿着城市街道）需要攀登一面 3.2 eV 的高墙。
  • 垂直移动（向上或向下）需要攀登一面 3.8 eV 的高墙。
  • 斜向穿越城市则更加困难，需要攀登一面 10 eV 的高墙。

类比： 想象试图将一块巨石推上山。即使你给巨石一个巨大的推力（将材料加热到极端温度），它也几乎纹丝不动，因为山实在太陡峭了。

2. “直接交换”与“群体舞蹈”的失败

研究人员还检查了硅是否可以通过与邻居直接交换位置，或者通过三个原子同时进行的复杂“群体舞蹈”来移动。

• 结果： 这些方法甚至更加不可能。所需的能量就像试图跳过摩天大楼（超过 12 eV）。
• 结论： 硅被困住了。除非它找到一个非常特定的空座，否则它根本不会移动；而即使找到了，攀登的过程也太过陡峭。

3. “极端高温”测试（实验）

计算机模型固然很好，但团队需要现实世界的证据。他们取来实际的氮化镓晶体，将硅植入其中，然后对其进行超高压退火（UHPA） 处理。

• 设置： 这就像把晶体放入一个既是高压锅又是熔炉的装置中。他们将晶体加热到超过 1300°C（比披萨烤箱更热），并在 30 分钟到 3 小时的时间内，施加巨大的压力（1 GPa）进行挤压。
• 测试： 他们使用一种特殊的显微镜（SIMS）拍摄了硅原子位置的“前后”对比照片。
• 结果： 硅原子纹丝不动。“前”和“后”的照片看起来完全一样。即使经过高温烘烤和高压挤压，硅原子依然停留在被安置的确切位置。

4. 这为何重要

该论文得出结论：氮化镓中的硅是一个极其忠诚的市民。

• 不游荡： 与某些其他材料不同，那些材料中的原子在受热时会变得不安并模糊界限，而氮化镓中的硅则安守原位。
• 精确性： 这意味着工程师可以在其电子器件中创建非常清晰、精确的边界，而无需担心制造过程中的热量会模糊设计。
• 一致性： 无论晶体是在蓝宝石基底上还是在氮化镓基底上生长，也无论硅是轻微植入还是大量植入，硅就是拒绝移动。

简而言之：
研究人员证明了，氮化镓中的硅就像飓风中的一座石像。无论施加多高的温度或多大的压力，它都始终停留在它该在的位置。这使得氮化镓成为构建下一代快速、强大且精确电子器件的完美、稳定的基础。

技术摘要

技术摘要：氮化镓中硅扩散的局限性：基于第一性原理计算与实验证据的 DFT 研究

问题陈述
硅（Si）是制造 n 型氮化镓（GaN）的主要有意施主，而 GaN 是用于高功率和高频光电子及电子器件的关键材料。虽然 Si 掺杂通常通过外延生长过程中的原位掺入或离子注入来实现，但对掺杂区域横向空间分布的精确控制仍然是一个挑战。尽管离子注入可提供精确的掺杂，但它需要高温退火以修复晶格损伤并激活掺杂剂。该过程中的一个关键问题是 Si 扩散的潜在风险，这可能会模糊先进器件架构所必需的锐利掺杂分布。

现有文献关于 GaN 中 Si 扩散的报道存在冲突。一些研究表明，在高温下存在可测量的扩散，且激活能较低（例如 0.89–1.55 eV）；而另一些研究则表明扩散可忽略不计，或由特定缺陷对介导的激活焓极高（例如约 10 eV）。此外，Si 在超高压退火（UHPA）条件下的行为——该条件用于在超过典型生长极限的温度下稳定 GaN——尚需澄清。本研究旨在通过第一性原理计算研究 GaN 中 Si 的扩散机制，并利用实验 UHPA 数据对其进行验证，从而解决这些差异。

方法论
本研究采用结合密度泛函理论（DFT）计算与实验表征的双重方法。

• 理论框架：

  • 计算设置： 计算使用 SIESTA 软件包在广义梯度近似（GGA）下完成，采用 PBE 势。使用了由 324 个原子组成的超胞模型（3a√3 × 3a√3 × 3c GaN 原胞）。
  • 缺陷建模： 研究模拟了 Ga 位上的替位 Si 原子（$Si_{Ga}$）和 Ga 空位（$V_{Ga}$）。为了反映实验掺杂条件，模拟了中性情形和 n 型情形（三负电荷 $V^{3-}_{Ga}$）。
  • 扩散机制： 使用 nudged elastic band (NEB) 方法确定 Si 迁移的最小能量路径（MEP）。分析了三个晶体学方向：a 方向 [11$\bar{2}$0]、c 方向 [0001] 和 m 方向 [$\bar{1}$100]。研究了空位介导机制以及直接交换（D-Ex）或环迁移机制。
  • 热力学： 进行了声子计算以确定振动自由能贡献（$\Delta F_{vib}$），从而基于过渡态理论计算温度依赖的有效能垒和扩散系数（$D$）。

• 实验验证：

  • 样品： 使用金属有机气相外延（MOVPE）和氢化物气相外延（HVPE）在各种衬底（Ammono-GaN 和蓝宝石）上制备了 Si 注入的 GaN 样品，这些样品具有不同的位错密度（TDD）。
  • 工艺： 样品在 1 GPa 的 N$_2$ 压力下，于高达 1300°C（特定持续时间下更高）的温度下进行 UHPA。
  • 表征： 使用 X 射线衍射（XRD）评估晶体恢复情况，同时利用二次离子质谱（SIMS）测量退火前后的 Si 浓度分布，以检测任何扩散现象。

主要结果

1. 能垒与机制：

   • 空位介导扩散： 主要的扩散机制是空位介导的。计算出的能垒很高：对于 n 型系统，沿 a 方向约为 3.2 eV，沿 c 方向约为 3.8 eV。m 方向的能垒异常高（约 10 eV），使得沿该轴的直接长程迁移极不可能。
   • 替代机制： 发现没有空位的直接交换和环迁移机制具有极高的能垒（11.8–14.8 eV），使其成为不太可能的路径。
   • 温度依赖性： 虽然振动自由能略微降低了有效能垒（对于中性系统，在 1800 K 时降低了约 0.25 eV），但能垒仍然很高。在 n 型系统中，降低幅度可忽略不计。
   • 扩散系数： 即使在 1800 K 时，计算出的扩散系数（$D$）也未超过 $10^{-9}$ cm$^2$/s。当考虑实际空位浓度（通常 $<10^{18}$ cm$^{-3}$）时，有效扩散速率预计至少低三个数量级。

2. 实验观察：

   • 结构恢复： UHPA 成功去除了注入引起的损伤，并在 MOVPE 和 HVPE 样品中恢复了高质量的晶体结构。
   • Si 分布稳定性： SIMS 深度分布显示，无论生长方法（MOVPE 与 HVPE）、衬底类型、TDD（范围从 $5 \times 10^4$ 到 $10^8$ cm$^{-2}$）或注入参数（能量高达 300 keV，剂量高达 $1 \times 10^{15}$ cm$^{-2}$）如何，UHPA 后 Si 浓度均未检测到展宽或偏移。

意义与主张
本文声称解决了以往文献中关于 GaN 中 Si 扩散的不一致之处。通过将严谨的 DFT 计算与高压实验验证相结合，作者证明了：

• 由于能垒较高，特别是对于空位介导的过程，GaN 中的 Si 扩散在本质上是受限的。
• 扩散具有各向异性，横向（a 方向）扩散快于垂直（c 方向）扩散，但在标准和极端加工条件下，两者均可忽略不计。
• 即使在 UHPA 条件下（温度 >1100°C，高压），Si 分布仍保持稳定，这与类似范围内存在可测量扩散的报道相矛盾。
• 位错密度、生长方法和衬底类型等因素并未显著增强 Si 的迁移率。

研究结论指出，对于需要精确 n 型掺杂的器件制造应用，GaN 中的 Si 掺杂是稳定的，因为掺杂区域因扩散而展宽的风险极小。这为理解 GaN 中 Si 的行为提供了全面的理论和实验基础，支持了可靠 GaN 基电子和光电子器件的开发。