Electric-field effects on defect migration energetics in GaN

技術要約：GaNにおける欠陥移動エネルギーへの電界効果

問題提起
窒化ガリウム（GaN）およびAlGaNベースのデバイスの信頼性、特に高出力および放射線環境下においては、点欠陥の移動と再構成によって損なわれる。照射による欠陥形成は十分に研究されているが、動作電界と温度の複合的な影響下における既存の欠陥の挙動については、理解が進んでいない。高電子移動度トランジスタ（HEMT）やパワーデバイスにおいて、動作電界はしばしば破壊限界（〜3 MV/cm）に達し、閾値電圧の不安定性や性能ドリフトを招く欠陥輸送を駆動する可能性がある。既存の計算手法にはトレードオフが存在する。密度汎関数理論（DFT）は高い精度を提供するが、小さなスーパーセルに制限される。一方で、標準的な原子間ポテンシャルを用いた古典的分子動力学（MD）は、外部電界下での欠陥移動をシミュレートするために必要な動的な電荷移動や静電結合をモデル化する能力を欠いている。

手法
これらの制限に対処するため、著者らはGaNに特化したReaxFF反応力場を開発した。その手法は主に以下の3つの段階で構成される。

1. 力場の開発： PBE汎関数フレームワークを用いたDFT計算（VASP）により、包括的なトレーニングデータセットを生成した。このデータセットには、様々なGaおよびGaN相（ウルツ鉱、岩塩、ジンクブレンドなど）の状態方程式（EOS）、N2分子プロファイル、および中性点欠陥（空孔、アンチサイト、フレンケル対、格子間原子、ショットキー欠陥を含む）のスーパーセルの全エネルギーが含まれている。
2. 検証： ReaxFFパラメータをこのDFTデータセットに対して最適化した。開発された力場は、格子定数、結合エネルギー、体積弾性率、および欠陥形成・移動エネルギーを、先行研究のDFTおよび実験結果と比較することで検証された。
3. 電界シミュレーション： 検証済みのRe型ReaxFFモデルをLAMMPS内で使用し、0から0.04 V/Å（0–40 MV/cm）の範囲の外電界下における固有点欠陥（$V_{Ga}$、$V_N$、$Ga_{i1,oct}$、$N_{i2,split}$）の移動障壁を決定するために、ぬげ弾性バンド法（CI-NEB）を実施した。シミュレーションでは、動的な電荷平衡（QEq）と局所的な格子歪みを考慮している。

主な貢献と結果

• 力場の検証： 開発されたReaxFFモデルは、DFTの構造特性（例：格子定数の誤差0.3%以内、体積弾性率の誤差1.5%以内）および欠陥エネルギーを正確に再現している。特定の構成において絶対的な欠陥形成エネルギーに偏差（最大26%）は見られるものの、モデルは文献値と比較して相対的な傾向と移動障壁を高い忠実度で捉えている。
• 非線形電界応答： 単純な双極子－電界近似（$\Delta E_B = -\mu_{eff} \cdot E$）とは異に、本研究は電界が移動障壁を異方的かつ非線形に変化させることを明らかにしている。この応答は、特定の欠陥の種類およびその移動メカニズムに強く依存する。
  • 空孔： 窒素空孔（$V_N$）の移動は、格子歪みや電荷変化が最小限であるため、剛体双極子モデルと一致する、電界に対してほぼ線形で対称な応答を示す。ガリウム空孔（$V_{Ga}$）は軽微な非対称性を示し、前方障壁が増加する一方で後方障壁は減少するが、その傾きは不均一である。これは、Gaの大きな原子サイズが局所的な歪みと電荷再分配を誘起することに起因する。
  • 格子間原子： 格子間原子（$Ga_{i1,oct}$ および $N_{i2,split}$）は顕著な非線形性と非対称性を示す。$Ga_{i1,oct}$ については、電界がレート決定ステップとなる鞍点構成をシフトさせるほど強力であり、どのエネルギーピークが障壁を支配するかというクロスオーバーを引き起こす。これらの欠陥は複数の原子の協同運動を伴うため、複雑な部分電荷の再分配と歪み応答が生じ、線形予測から大きく逸脱する。
• 偏差のメカニズム： 線形挙動からの逸脱は、反応座標に沿った電界誘起の部分電荷再分配と局所的な格子歪みから生じる。有効双極子モーメントは固定されておらず、欠陥が移動するにつれて連続的に進化する。これは反応力場によって自然に捉えられる現象であるが、固定双極子モデルでは見落とされる。

意義と主張
本論文は、電界は単に欠陥の運動を線形にバイアスするだけでなく、電荷－格子結合を通じて欠陥の移動景観を複雑に修飾し得ることを主張している。

• 速度論的影響： 高い障壁（>2 eV）を持つ空孔に対する電界誘起の障壁減少は、デバイスの時間スケールにおいて拡散速度を大きく変えるほどではないが、障壁が1 eV付近の格子間原子については速度論的に重要である。約0.05 eVの減少は、室温において移動速度を1桁加速させる可能性があり、高電界領域における欠陥再分布を加速させる可能性がある。
• デバイスの信頼性： これらの知見は、高電界領域（例：HEMTのゲートエッジ）において、格子間原子が大幅に移動しやすくなり、トラップ生成や閾値電圧ドリフトに寄与することを示唆している。応答の方向的な非対称性は、欠陥が特定のデバイス領域へと正味のドリフトを起こすことを意味しており、長期的な信頼性の問題を引き起こす可能性がある。
• 手法の範囲： 著者らは、QEq電荷平衡スキームの古典的な性質上、得られた結果は中性欠陥に関するものであると述べている。実際の動作条件下での高度に荷電した欠陥（例：$V_{Ga}^{3-}$）は、さらに強い電界応答を示す可能性があることを認めており、現在の知見は電界効果の下限を示していることを示唆している。それにもかかわらず、本研究は、電界と欠陥の結合が異方的かつ欠陥特異的であり、現実的な動作条件下での正確な予測には反応性シミュレーションが必要であることを確立している。