InGaNBi四元系合金におけるバンドギャップエンジニアリング：組成依存の構造的、電子的、歪み特性

概要

InGaNBi合金は、可視から中赤外領域までバンドギャップを調整可能な多機能プラットフォームとして次世代フォトニックデバイスに有望である。第一原理密度汎関数理論を用い、インジウム（In）とビスマス（Bi）濃度が、GaN基板上での成長時の格子定数、バンドギャップエネルギー、スピン軌道分裂、歪みにどう影響するかを体系的に調査した。格子定数はInとBi増加にほぼ直線的に増大し、Bi最大9.375%、In最大50%でバンドギャップを3.273 eVから0.651 eV（0.38–1.9 µm）まで調整可能。計算されたスピン軌道分裂は、Bi 3.125%、6.25%、9.375%でそれぞれ0.220 eV、0.360 eV、0.600 eVに達する。In/Bi比を調整することで、GaN上で許容可能な歪みを維持しつつ合金を設計でき、高効率LEDとレーザーダイオードへの道を開く。

はじめに

ウルツ鉱型InₓGa₁₋ₓN合金とInGaN/GaN量子井戸（QW）は、調整可能なバンドギャップとGaN基板との互換性により、高効率発光ダイオード（LED）とレーザーダイオード（LD）への応用で強い関心を集めてきた。しかし、[0001]配向InₓGa₁₋ₓN/GaN QWは、双軸圧縮歪みによる強い内部電場に悩まされ、電子-正孔対の発光エネルギーと振動子強度を低下させる。また、幾何学的欠陥（積層欠陥、スレッディング転位）の高密度が非放射再結合中心となり、電子漏れとオーゲ再結合がInGaN LEDの悪名高い効率ドロープに寄与する。

赤外領域では、GaAsへのビスマス（Bi）添加がバンドギャップエネルギー（E_g）低減とスピン軌道（SO）分裂増大に有効であり、オーゲ再結合抑制につながる戦略として浮上。Biは最大のV族元素で、III-Vビスマイド合金のバンド構造に大きな影響を与える。AlNBi、GaNBi、GaSbBi、InPBi、InSbBiなどで実証されている。これらの三元系では、Bi導入がBi不純物状態と重い/軽い正孔帯の相互作用を通じて価電子帯（VB）を摂動し、Biの大きな原子サイズが歪み関連バンドギャップ変調を誘起する。より最近では、GaAsNBi、InGaAsBi、GaAsPBiなど四元系ビスマイド系が、バンドギャップ、SO分裂、歪みを同時エンジニアリングできる拡張された調整可能性で注目されている。Biを他のIII族またはV族元素と組み合わせることで、光電子デバイスの機能設計空間を拡大できる。

InGaNBiが可視と中赤外応用を橋渡しする可能性を鑑み、第一原理計算でInとBi濃度が合金の構造的、電子的、歪み特性をどう形成するかを解明する。In ≤ 50%、Bi ≤ 9.375%に濃度を制限し、大きな格子ミスマッチと限られたBi固溶度による結晶品質低下を回避。本稿は以下の構成：第2節で計算手法、第3節で格子定数、バンドギャップ、SO分裂、歪みの結果と考察、第4節で要点を総括。

方法

全計算をウィーンアブイニシオシミュレーションパッケージ（VASP）実装の密度汎関数理論（DFT）で実施。構造緩和にプロジェクタ拡張波（PAW）法、交換相関にPerdew-Burke-Ernzerhof（PBE）一般化勾配近似（GGA）を使用。PBEの系統的バンドギャップ過小評価を補正するため、修正ベッケジョンソン（MBJLDA）交換ポテンシャルと局所密度近似（LDA）相関を組み合わせて使用。Biは強いスピン軌道結合（SOC）を示すため、全電子構造計算でSOCを明示的に含めた。

価電子構成をIn: 4d¹⁰5s²5p¹、Ga: 3d¹⁰4s²4p¹、N: 2s²2p³、Bi: 5d¹⁰6s²6p³に設定。原子力を< 0.02 eV/Å、全エネルギーを10⁻⁴ eVまで収束。平面波カットオフ450 eV、4 × 4 × 4モンクホルストパックk点メッシュを用い、正確な全エネルギーと電子特性を確保。

結果と考察

構造的特性

4 × 2 × 2スーパーセル（64原子）で、0 ≤ x ≤ 0.09375、0 ≤ y ≤ 0.5の36組成のIn_yGa₁₋yN₁₋xBi_xをモデル化。各ケースでInとBi原子をできるだけ均一に分布させクラスタリング効果を最小化。三元系In_yGa₁₋yNと四元系In_yGa₁₋yN₁₋xBi_xの計算格子パラメータを図1にプロット。純粋GaNでa = 3.211 Å、c = 5.235 Åを得、先行理論（a = 3.155–3.22 Å、c = 5.144–5.24 Å）と実験値（a = 3.19 Å、c = 5.19 Å）と優れた一致。In_yGa₁₋yNの格子定数はy増加にほぼ直線的に増大し、y = 0.25でa = 3.304 Å、c = 5.365 Å、y = 0.5でa = 3.397 Å、c = 5.509 Åに達し、先行研究と整合。四元系合金では、InとBi含有量増加に伴いaとcがほぼ直線的に上昇。これはIn（1.42 Å）とBi（1.48 Å）のイオン半径がGa（1.22 Å）とN（0.71 Å）より大きいため。

代表組成In₀.₂₅Ga₀.₇₅N₀.₉₃₇₅Bi₀.₀₆₂₅の結合長解析（図2）で、Ga-N結合が2.009 Åに短縮、In-N結合が2.195 Åに伸長、Ga-BiとIn-Bi結合がそれぞれ2.592 Åと2.704 Åに達することを確認。これらの傾向はInとBi原子が局所格子を拡大することを裏付け、共有結合半径の大きさと整合。このような格子歪みと電気陰性度差は電子的・光学的挙動に影響を与えると予想される。

電子的特性

バンドギャップエネルギーをMBJLDA+SOC法で計算。この手法はIII-V合金で信頼性が実証済み。図3にIn_yGa₁₋yNのIn含有量依存E_gを、実験と他理論ベンチマーク（HSE06、mBJ、LMTO-CPA-MBJ）と共にプロット。計算GaNバンドギャップ（3.273 eV）は実験値（3.40–3.50 eV）と先行計算（3.261–3.33 eV）とよく一致。yが0から0.5に増加すると、バンドギャップは3.273 eVから1.546 eVへ滑らかに減少し、理論と実験両方と一致。

四元系In_yGa₁₋yN₁₋xBi_x合金に拡張すると、図4の等高線マップでInとBi含有量の両方でバンドギャップが顕著に非線形減少することを示す。Bi最大9.375%、In最大50%でバンドギャップは3.273から0.651 eVまで変化し、波長0.38–1.9 µmに対応、可視から中赤外全域をカバー。重要なのは、Bi導入がE_gを低下させるだけでなくスピン軌道分裂Δ_SOを増大させ、Bi 3.125%、6.25%、9.375%でそれぞれ0.220 eV、0.360 eV、0.600 eVに達すること。これらの値はIn分率にほぼ依存せず、BiがΔ_SO増大の主要ドライバーであることを示唆。より高組成合金でΔ_SO > E_gを達成できれば、高効率LEDとLDの重要なボトルネックであるオーゲ再結合を抑制可能。

GaN、In₀.₂₅GaN、In₀.₂₅GaN₀.₉₆₈₇₅Bi₀.₀₃₁₂₅の投影バンド構造と全状態密度（TDOS）（図5）で、InとBiの異なる役割が明らかに。In導入は伝導帯底（CBM）と価電子帯頂（VBM）の両方を低エネルギー側へシフトさせ、ギャップを効果的に狭める。対照的にBiはフェルミ準位直下に欠陥帯を導入し、VB端と強く混成してVBMをさらに押し上げる。TDOSに-1.0から-0.5 eVの間で顕著なBi由来ピークが現れ、キャリアダイナミクスを変調する局在状態形成を確認。

GaN上でのInGaNBi歪み

実用デバイス統合には、GaN上に成長したIn_yGa₁₋yN₁₋xBi_x層の歪み状態を管理可能な範囲内に保つ必要がある。図6にInとBi分率関数の双軸歪みをマッピング。最大検討組成（y = 0.5、x = 0.09375）で合金は8.5%の圧縮歪みを経験—大きいが適切なバッファ層で許容可能。より重要なのは、In分率≤0.0625、Bi分率≤0.028で歪みが1%未満に留まり、控えめな合金化で準格子整合系が得られ、高品質エピタキシャル成長に適することを示す。

結論

第一原理DFT計算により、InGaNBi合金が可視から中赤外スペクトルまで光電子応用向けに高度に調整可能なプラットフォームを提供することを実証。格子定数はInとBi含有量にほぼ直線的に増大し、InとBiの大きな原子サイズを反映。バンドギャップはBi最大9.375%、In最大50%で3.273 eVから0.651 eVまで設計可能、スピン軌道分裂は最高Bi分率で0.600 eVに達し、Δ_SO > E_g時にオーゲ再結合抑制可能。バンド構造解析でInは主にCBMとVBMをシフト、BiはVBMを押し上げる局在欠陥状態を導入することを示す。歪み解析で、適度なIn（≤6.25%）とBi（≤2.8%）濃度では合金がGaNとほぼ格子整合し、高品質エピタキシャル成長を可能にする。これらの知見は、調整可能な発光波長と低減された非放射損失を持つ高効率InGaNBi LEDとレーザーダイオード設計の基盤を築く。