界面粗さ散乱が鍵：GaNテラヘルツ量子カスケードレーザ性能の決定要因

要旨

非平衡グリーン関数（NEGF）形式を用い、界面粗さ（IFR）、非意図的ドーピング（n.i.d.）、合金無秩序（AD）がGaNベーステラヘルツ量子カスケードレーザ（QCL）の光利得に及ぼす影響を調査した。シミュレーション結果、合金無秩序の影響は無視できる一方、非意図的ドーピングは1×10¹⁷ cm⁻³以下に抑制し、電子-不純物散乱と自由キャリア吸収を回避する必要がある。最も重要なのは、界面粗さ散乱が光利得劣化の主要因であり、製造時の精密制御の必要性を強調する点である。モデルは300 Kでピーク利得60 cm⁻¹を予測し、室温動作GaN THz QCLの実現可能性を示した。

はじめに

テラヘルツ（THz）スペクトル帯は、セキュリティ検査、医療イメージング、高速通信への応用が期待されるが、実用的でコンパクトな光源は依然として乏しい。量子カスケードレーザ（QCL）は固体THz光源の有力候補だ。GaAsベースQCLの最高動作温度は~200 Kだが、GaNの高縦光学（LO）フォノンエネルギー（92 meV）は室温動作と広いアクセス周波数範囲（1–15 THz）への道を開く。最近のGaN THz QCL進展では、10–13 THzでのサブバンド間（ISB）吸収実証と、理論モデルでLOフォノン相互作用が主要制限因子と指摘されている。しかし、成長関連欠陥（IFR、n.i.d.、AD）が光利得に及ぼす影響は十分に定量化されてこなかった。本研究は、これら散乱メカニズムが光利得に与える影響を定量化し、IFRが主要劣化因子であることを特定することで、既存研究を拡張する。

手法

高品質GaN THz QCLには、厚い活性領域と低転位密度が必要だが、GaN/AlGaN格子ミスマッチにより困難を伴う。エピタキシャル成長は追加散乱源を導入する：表面粗さ由来のIFR、残留不純物（特に酸素）由来のn.i.d.、Ga表面偏析とAlアダトム移動度由来のAD。Nextnano QCLスイートに実装されたNEGF形式を用い、現在GaN THz QCLで最高動作温度構成である3量子井戸共鳴フォノン設計をモデル化。1周期の層系列（nm単位）は1.6/6.2/1.6/3.4/1.0/3.4で、6.2 nm井戸が5×10¹⁷ cm⁻³でnドープ。図1に設計構造、伝導帯プロファイル、包絡関数二乗、キャリア密度を示す（バイアス-84.5 kV cm⁻¹、10 K）。

設計活性領域構造、伝導帯プロファイル、二乗包絡関数、キャリア密度。a 1周期層系列：1.6/6.2/1.6/3.4/1.0/3.4 nm。バリアはイタリック体で表示。6.2 nm厚井戸はn=5×10¹⁷ cm⁻³でnドープ。b GaN/Al₀.₁₅Ga₀.₈₅N QCLの伝導帯プロファイルと二乗包絡関数。電界-84.5 kV cm⁻¹印加。c NEGFモデルで計算したQCLのキャリア密度と伝導帯。電界-84.5 kV cm⁻¹、温度10 K。

シミュレーションでは、PAMBE成長GaN/AlGaNの現実的パラメータを使用：IFR振幅0.25 nm、相関長1 nm；n.i.d.濃度1×10¹⁷ cm⁻³；AD散乱も包含。これらは典型的実験条件を反映する。

結果と考察

図2(a)に4構成の計算光利得スペクトルを示す：フル参照モデル（IFR、n.i.d.、AD）、および各散乱機構を1つずつ除いた3モデル。10 Kで参照構造は8.7 THzでピーク利得142 cm⁻¹。n.i.d.除去でピークが127 cm⁻¹（8.46 THz）へ低下。これは上部状態電子数が追加ドナー電子除去で減少するため。しかし閾値電流密度も低下し、不純物散乱がキャリア輸送を阻害することを示唆。AD除去ではピーク利得がわずか3%増の147 cm⁻¹（同周波数）で、低Al含有（15%）設計では合金無秩序の役割は限定的。対照的に、IFR除去でピークが191 cm⁻¹（8.7 THz）へ上昇し閾値電流低下、34%の利得向上。IFRがレーザー遷移から電子を散乱させ、光利得損失を支配している。

異なる散乱機構でのシミュレーション最大光利得vs周波数と電流-電界特性。a 光利得スペクトル。b 電流-電界特性。温度：10 K。

IFR感度を探るため、相関長1 nm固定で粗さ振幅を0.5 nmと0.75 nmに変化。図3に劇的低下を示す：ピーク利得は0.5 nmで47.9 cm⁻¹、0.75 nmで負値（-8.8 cm⁻¹）となりレーザー発振を抑制。電流-電圧曲線ではIFR増大に伴い電流密度上昇と共鳴トンネリング減衰が見られ、現実的条件下ではIFR散乱が共鳴トンネリングを上回ることを確認。これら結果は先行GaAs研究と整合し、エピタキシー時IFRを0.5 nm以下に抑える必要性を強調する。

異なるIFR振幅でのシミュレーション最大光利得vs周波数と電流-電界。a 利得スペクトル。b 電流-電界。温度：10 K。

最後に、参照モデル（IFR=0.25 nm、n.i.d.、AD）の温度依存性を評価。図4に示す通り、ピーク利得は150 Kまで~142 cm⁻¹で安定、その後低下し300 Kで61 cm⁻¹。低下は熱的バックフィリングとLOフォノン散乱増強によるが、利得はGaN二重金属波導の30 cm⁻¹損失を上回り、室温動作可能性を確認。GaAs THz QCLと比較し、GaNの高ドーピング、低屈折率、薄い周期が高光利得をもたらし、10 Kで142 cm⁻¹ピークを記録。

格子温度に対する計算最大利得の特性。

結論

8.7 THz動作GaN THz QCLの包括的NEGF研究を提示。設計は10 Kで142 cm⁻¹、300 Kで60 cm⁻¹利得を生成し、室温実現可能性を確認。合金無秩序の影響は微小（<5%利得変化）。非意図的ドーピングは≤1×10¹⁷ cm⁻³に制限し、バンドミスアラインメントと過剰不純物散乱を回避。界面粗さが光利得劣化を支配；IFRを0.5 nm以下に抑えることが正利得に不可欠。これら知見は、実用的高温デバイスへ向けたGaN THz QCL製造最適化の具体的指針を提供する。

データと材料の入手可能性

本研究で生成・解析されたデータセットは、責任著者から合理的な要請に基づき入手可能。

略語

AD

合金無秩序

IFR

界面粗さ

ISB

サブバンド間

n.i.d.

非意図的ドーピング

NEGF

非平衡グリーン関数

QCL

量子カスケードレーザ