Si、Ge、Si/Ge超格子における低エネルギー放射応答のアビニシオシミュレーション

概要

アビニシオ分子動力学を用いて、低エネルギーイオン照射がバルクシリコン（Si）、バルクゲルマニウム（Ge）、Si/Ge超格子（SL）に及ぼす影響を調べました。研究の結果、SL界面のGe原子はバルクGeと同様の変位抵抗を保持する一方、界面近傍のSi原子はバルクSiより高い変位閾値を示し、放射線耐性が向上することが明らかになりました。欠陥生成メカニズムはバルクとSL構造で異なり、SLではより複雑な欠陥クラスターが生成されます。欠陥形成エネルギーと移動障壁の計算から、SLにおける点欠陥は形成しにくく、空孔の移動性が低いことが示されました。これらの知見は、Si/Ge SLが放射線環境下で動作する電子・光電子デバイスに有望であることを示唆しています。

背景

近年、Si/Ge超格子は調整可能なバンド構造と向上したキャリア移動度により、次世代電子・光電子デバイスの汎用プラットフォームとして台頭しています。用途は光通信用高速フォトダイオード、マイクロエレクトロニクス、太陽電池、宇宙搭載部品に及びます。しかし、宇宙空間での高エネルギーイオン曝露はこれらの材料の光学的・電子的特性を変化させ、デバイス劣化を招きます。したがって、極限条件下でのSi/Ge SLの放射応答を理解することは不可欠です。

実験的研究では、Si/Ge SLのバルクSiに対する優れた放射線耐性が記録されています。例えば、単層Ge埋め込みSLへの電子照射ではバルクシリコンに比べ光発光劣化が減少し、Ge量子ドットin SL構造への2MeV陽子照射では例外的な耐性が示されました。Si/Ge多層膜上のGe量子井戸デバイスや熱電測定でも同様の傾向が観察され、イオン照射量によりフォノン散乱と電子輸送が強化され性能指標が向上しました。これらの知見にもかかわらず、Si/Ge SLにおける欠陥形成と進化の包括的な原子論的理解は依然として不足しています。

アビニシオ分子動力学（AIMD）は、第一原理電子構造計算から導かれる原子間ポテンシャルを用いて、原子スケールでの放射線損傷を探る強力なツールです。古典的MDとは異なり、AIMDは電荷移動と電子遮蔽効果を捉え、閾値変位エネルギー（E_d）と欠陥構造の正確な決定を可能にします。本研究では、バルクSi、バルクGe、Si₂/Ge₂ SLにおける低エネルギー変位イベントをAIMDで比較し、それらの異なる放射応答のメカニズムを特定し、欠陥形成と移行のエネルギーを定量化しました。

方法

変位カスケードはSIESTAコードでシミュレートしました。ノルム保存トロイエ・マルティンス擬ポテンシャルでイオン-電子相互作用を記述し、セペルリ・アルダー・パラメータ化による局所密度近似（LDA）で交換相関効果を処理しました。価電子状態は単一ζプラス偏極（SZP）基底セットで展開し、1×1×1 k点メッシュと60Ryカットオフを使用。SLのモデル超電池は、2層のSiと2層のGeが交互に配列した288原子で構成（Si₂/Ge₂）。図1にバルクSiとSi/Ge SLの幾何構造を示します。

一次ノックオン原子（PKA）を選択し、運動エネルギーを与えてリコイルイベントを開始。PKAが格子サイトに戻った場合、5eV刻みでより高いエネルギーでシミュレーションを繰り返し、PKAが変位した時点で閾値を0.5eV精度で絞り込みました。各元素について、バルクとSLで4-5方位の結晶方向を検証。全ランでNVEアンサンブルを使用し、数値不安定性を避けるため1.2psに制限しました。

結果と考察

バルクシリコンとゲルマニウムにおける変位イベント

AIMD計算により、バルクSiの格子定数5.50Å（実験値5.43Å、理論値5.48Åと一致）、バルクGeで5.71Å（実験値5.77Å、理論値5.65Åと整合）を得ました。表1に閾値変位エネルギー（E_d）と結果の欠陥構造をまとめ、図2-3にSiとGeの損傷構造を示します。

バルクSiでは、計算E_d値は実験測定値に近く：[001]で21eV、[110]で約47.6eV、[111]で約12.9eV。全ケースで主要欠陥タイプはフレンケル対（空孔+インタスチシャル）です。[110]方向のE_d（47eV）は古典的MD値（24eV）を上回り、AIMDにおける電荷移動効果の包含を反映。変位メカニズムは方向により異なり：[111]ではPKAがインタスチシャルサイトへ散乱し変位原子が格子サイトへ復帰、[001]ではより複雑な交換系列で2空孔2インタスチシャル生成、[110]では多重衝突でフレンケル対クラスター形成。これらはバルクSiで[110]方向が最も変位に強いことを確認します。

バルクGeでは実験・理論と整合するE_d値：[001]で約18eV、[111]と[1 1 1]で9.5eV。[111]方向では単純なインタスチシャル形成、[110]ではPKAが隣接Ge原子と位置交換し1空孔1インタスチシャル生成。Siと同様に[110]方向が最も変位困難。バルクSiとGeの全欠陥構造はフレンケル対が支配的です。

Si/Ge超格子における変位イベント

Si/Ge界面に隣接するPKA原子を選択しSi₂/Ge₂ SLを調査。表2にE_d値と関連欠陥を記載、図4-5に結果構造を示します。

SLにおいて、Si原子は[1 1 1]方向で高耐性（E_d=10eV）を示し、Si[001]とSi[00 1̅]ではそれぞれ46.5eVと42.5eVを要し複雑な反サイトクラスターを生成。Ge原子は[111]と[1 1 1̅]でバルクGeと類似の低E_d（約16eV）を保持するが、[001]（約17eV）と[00 1̅]（約18eV）でより高い閾値を示します。重要なことに、SL中のGe[110]はE_dが低下（8.5eV）し、混合欠陥クラスター（空孔、反サイト、インタスチシャル）を生成。これらの知見は、SLがSiの放射線耐性を向上させ、Geについてもある程度向上させることを確認し、従来の実験観察と一致します。

バルクSi、Ge、Si/Ge超格子における欠陥形成エネルギーと移動障壁

64原子超電池、6×6×6 k点サンプリング、500eVカットオフを用い、バルクSi、バルクGe、SLにおける空孔、インタスチシャル、反サイトの欠陥形成エネルギーを計算。結果（表3）は、バルク材料ではV_Si（3.60eV）とV_Ge（2.23eV）が最も有利な欠陥である一方、SLではほとんどの欠陥でわずかに高い形成エネルギーを示し、V_Ge（2.73eV）とV_Si（2.85eV）が依然として最低値です。

SL界面に隣接するV_GeとV_Siの移動障壁をnudged elastic band法で評価（表4）。空孔はバルクとSLの両方で[111]方向に最も移動しやすいが、障壁はバルクで一貫して低い（例：V_Si [111]方向：バルク0.11eV vs SL 0.12eV）。対照的に[100]と[110]方向ではわずかな差。これらはSLで空孔移動性が抑制され、照射下での欠陥クラスタリングと膨張を減らす可能性を示唆します。

結論

AIMD調査により、バルクSiとGeの閾値変位エネルギーは強い異方性を持ち、[110]方向が最も耐性が高いことが示されました。Si/Ge SLでは、Si原子が[111]と[001]方向で耐性向上を示し、Ge原子は[110]で低下した値を除きバルク類似の閾値を維持。欠陥形成エネルギーはSLで点欠陥の生成が困難であることを確認し、移動障壁は特に最も有利な[111]方向で空孔移動性が低下することを示します。これらの複合効果がSi/Ge超格子の優れた放射線耐性を支え、過酷な放射線環境下で動作するデバイスへの適用を魅力的にしています。