GaAs/AlGaAsヘテロ接合太陽電池のシミュレーション

1. はじめに：GaAs太陽電池とその応用

ガリウムヒ素（GaAs）系太陽電池は、変換効率の高さ、軽量性、信頼性がコストよりも重要視される用途に用いられる。GaAsは直接遷移型の半導体であり、そのバンドギャップは約1.42 eVであるため、わずか数マイクロメートルの厚さの材料内で太陽光を効率的に吸収することができる。この特性から、高性能な太陽光発電デバイスに最適な材料と言える。

GaAs系太陽電池は、表面面積が限られており放射線耐性が求められる衛星や宇宙機器において広く利用されている（??）。また、太陽光を小型の半導体セル上に集光させる集光型太陽光発電や、高高度飛行機やレーザー電力変換器といった特殊用途にも用いられる。現代の宇宙用太陽電池では、より高い変換効率を得るためにGaAsと他のIII–V族半導体を組み合わせた多接合構造が採用されることが多い。

このチュートリアルでは、GaAs/AlGaAsヘテロ接合太陽電池のシミュレーションを行う。より広いバンドギャップを持つ薄いAlGaAsウィンドウ層をGaAs接合の上に配置することで、不要な光吸収や表面再結合による損失を抑制する。OghmaNanoを用いて、光の伝達行列法による計算とドリフト拡散モデリングを組み合わせ、このデバイスがいかにして太陽光を吸収し、光生成キャリアを輸送して電力を生み出すかを解析する。

2. GaAs/AlGaAs太陽電池の作成とその仕組みの理解

GaAs/Al0.3Ga0.7Asヘテロ接合太陽電池のシミュレーションを行う。この太陽電池は、GaAsのp-n接合の上に配置された薄くてバンドギャップの大きいAlGaAsウィンドウ層から構成される。比較的単純な構造ではあるが、このデバイスは光吸収、キャリア生成、ドリフト拡散輸送、再結合、ならびに半導体ヘテロ接合がキャリア収集に与える影響といったIII-V族半導体太陽電池の重要な側面をよく示している。

まずOghmaNanoを起動し、ファイルリボンから「新規シミュレーション」をクリックします。これによりデバイスライブラリが開きます（??）。次に「GaAsデモ」をダブルクリックして利用可能なGaAsシミュレーション一覧を表示します（??）。「GaAs/Al0.3Ga0.7Asヘテロ接合太陽電池」を選択し、保存先フォルダを指定して「保存」をクリックしてください。

読み込み完了後、OghmaNanoはデバイス構造タブ（??）にてデバイス構造を表示します。緑色の矢印は上方からの光照射を示しており、色付きの領域が4つの半導体層を表しています。この構造や材料特性は、レイヤー編集器および電気的パラメータ編集器を用いて変更可能です。

このデバイスは、光が入射する表面側から裏面電極に向かって、4つの半導体層で構成されている：

GaAs吸収層は、入射した光子を電子–正孔対に変換する役割を担う。GaAsの直接遷移型バンドギャップは約1.42 eVであり、この値は高い光吸収係数をもたらすため、比較的薄い半導体層内でも太陽光を効率的に吸収できる。p型発電層とn型吸収層が主な光電変換接合を形成しており、この接合を越えて拡散したキャリアによって空乏層および内部電界が生じ、光生成された電子と正孔がそれぞれ対応する電極へと輸送される。Al0.3Ga0.7Asからなるウィンドウ層は、光照射面にヘテロ接合を形成する。アルミニウムの添加によりGaAsよりもバンドギャップが拡大するため、下層のGaAsで吸収可能な光子に対してこのウィンドウ層はより透明性を発揮する。また電子親和力およびバンドギャップの差異により界面において伝導帯および価電子帯に不連続性が生じ、これが少数キャリアの表面への移動を抑制し再結合損失を低減する効果がある。ただし過大なバンドオフセットは逆にキャリアの抽出を妨げる可能性もある。高濃度にドーピングされたp+ウィンドウ層およびn+裏面接触層は、低抵抗な電気的接続の実現にも寄与している。裏面側では、低濃度ドーピングの吸収層と高濃度ドーピングの接触層との境界部に新たな電界が形成され、これが裏面接触からの少数キャリアである正孔の移動を抑制する。以上のように、各層の厚さ、ドーピング濃度、およびバンド配列の組み合わせが、吸収された光子が電流へと変換される効率を決定するのである。

メインウィンドウの「デバイス構造」タブ左側にある「レイヤーエディター」ボタンをクリックすることで、デバイス構造を確認・変更できます。このボタンは、多層構造のデバイスを表す3本の色付き水平バーで示されています。クリックするとレイヤーエディター（??）が開き、各層の厚み、光学特性を持つ材料、およびその種類を変更可能です。太陽電池においては特に層の厚みが重要であり、これが光の吸収量や光生成されたキャリアが電極に到達するまでに移動しなければならない距離を決定するからです。

電気関連リボン（??）を使用すると、ドーピング／イオン編集器および勾配／ひずみ編集器にアクセスできる。ドーピング編集器（??）では、デバイス全体におけるドナーおよびアクセプターの濃度を制御する。これらの値が各層における空乏層の位置や幅、内部電界の強さ、ならびに電気伝導度を決定する。したがって、ドーピング条件の変更はキャリアの収集効率、再結合の様相、そして最終的なJ–V特性に大きな影響を与えることになる。

グラディエント／ひずみエディタ（??）では、材料組成の変数である\(x\)、\(y\)、\(z\)を制御する。これらの値は各層の始端および終端において個別に指定可能であり、均一組成の半導体構造だけでなく組成が傾斜した構造もシミュレーションできる。AlxGa1-xAsのような合金においては、アルミニウム含有率がバンドギャップやGaAsとのバンドアライメントに影響を与える。したがって合金組成を変更することで、ヘテロ接合におけるバンドオフセットがキャリア輸送および再結合に与える影響を調査することが可能となる。各組成変数の意味は、選択された半導体材料モデルによって定義される。

3. シミュレーションの実行と結果の確認

GaAs/AlGaAs太陽電池のシミュレーションを実行するには、青色のボタンをクリックしてください シミュレーションを実行 ボタンが存在する ファイル リボン形式、または押下操作 F9。開く 端子 タブをクリックして計算状況を確認してください??）。初めのうちは、端末出力が複雑に見え、画面上に多数の数値情報が表示されるかもしれません。しかし、各項目が何を表しているかを理解すれば、その出力内容は比較的容易に解釈できるようになります。

シミュレーションの第一段階として光学計算が実行される。この様子は??の上部に表示される色付きテキストで確認できる。OghmaNanoは入射太陽光スペクトル内の各波長における光吸収および光生成の様子を計算する。この計算が完了すると、プログラムは電流–電圧（J–V）スイープ処理を開始し、太陽電池に一連の電圧を印加してそれに伴う電流密度を算出する。

J–V計算における各ラインは、それぞれ1つの印加電圧に対応している。top=で示される値は上部電極に印加される電圧であり、btm=は0 Vに固定されている下部電極の電圧である。計算された電流密度は、印加電圧の隣にA/m²単位で表示される。f()=は電気回路ソルバーの数値残差であり、計算結果が半導体方程式をどの程度正確に満たしているかを示す指標である。この残差値が小さいほど、数値計算の収束性が良好であることを意味する。

印加電圧が増加するにつれて、光電流の大きさは徐々に減少していく。図??において、電圧スイープは開回路電圧VOCに近づいている。約0.92 Vから0.94 Vの範囲で電流密度の符号が負から正へと変化しており、これは正味の電流がゼロを通過したことを示している。これが開回路状態であり、この時点では順方向ダイオード電流が光生成電流と正確に釣り合っている。シミュレーションはVOCを超えても継続され、J–V特性曲線の順バイアス領域を解析した後、指定された最大電圧である1.25 Vに達した時点で終了する。

各印加電圧において、OghmaNanoは計算処理を行う ポアソンの方程式 それと共に ドリフト拡散方程式およびキャリア連続の式これらは、デバイス全体における静電ポテンシャル、電子・正孔濃度、およびキャリア電流の変化の様子を記述するものである。

\[ \nabla\cdot(\varepsilon\nabla\varphi) = -q\,(p-n+N_D^+-N_A^-) \] \[ \nabla\cdot\mathbf{J}_n = q(R-G), \qquad \nabla\cdot\mathbf{J}_p = -q(R-G) \]

ここで\(G\)は光学モデルから算出されるキャリア生成率であり、\(R\)はキャリア再結合を表す。この再結合には放射再結合、ショックレー・リード・ホール再結合、オージェ再結合などが含まれる可能性がある。生成、再結合、およびキャリア抽出のバランスが太陽電池の電気的出力を決定する。より詳細な背景情報についてはドリフト拡散モデリングの理論を参照されたい。

計算が完了したら、Outputタブ（??）を選択してください。OghmaNanoはシミュレーション結果をシミュレーションディレクトリに書き出します。その中にはJ–V特性曲線が記載されたjv.csvや、主要な太陽電池性能指標が記載されたsim_info.datも含まれます。optical_outputおよびsnapshotsディレクトリには光学的・電気的特性に関するさらなる情報が格納されています。詳細については出力ファイルの説明書をご参照ください。

jv.csvをダブルクリックすると、照射下でのJ–V特性曲線が表示される（??）。逆バイアスおよび小さな順バイアス下では電流値はほぼ一定のままであるが、印加電圧が開回路電圧に近づくにつれて急激に増加する。次にsim_info.datをダブルクリックすると、算出されたデバイスの性能指標を確認できる（??）。

GaAs/AlGaAs太陽電池のシミュレーションで得られた主要な結果を以下にまとめる。

短絡電流密度JSCとは、印加電圧がゼロの状態における光電流の値である。負の符号は、太陽電池から流れ出る電流がダイオードの順方向電流とは逆向きであるという規約を反映している。その値は約26.9 mA cm⁻²であり、入射した太陽光の大部分が電荷キャリアに変換されていることを示している。

開回路電圧であるVOCとは、回路内の正味電流がゼロとなる状態における電圧のことです。この状態では光生成によって生じたキャリアが再結合と釣り合っており、デバイスは外部回路へ電流を供給できなくなります。計算されたVOCの値は0.922 Vであり、これは内部電位差の値である1.387 Vよりも低くなっています。この差は再結合現象や照射下で維持される有限のキャリア密度の影響を反映しています。詳細については内部電位差の計算をご参照ください。

フィルファクターは、J–V特性曲線が理想的な長方形形状にどれだけ近づいているかを示す指標であり、一方で電力変換効率は入射光エネルギーのうち電気エネルギーに変換される割合を表します。

\[ \mathrm{FF}=\frac{|J_{mp}|V_{mp}}{|J_{SC}|V_{OC}}, \qquad \mathrm{PCE}=\frac{|J_{mp}|V_{mp}}{P_{in}} \]

AM1.5G照射条件下で入射光強度が1000 W/m²の場合、計算された最大出力密度は200.48 W/m²となり、これは 電力変換効率は20.05%であるフィルファクターが約80.8%であることから、シミュレーション対象のデバイスは動作電圧範囲の大部分において十分な光電流を維持していることがわかる。最大出力点は約0.80 V付近に存在するが、その正確な位置は設定された電圧分解能に依存する。 J–Vエディタ.

シミュレーションにより得られた効率値は約 20% これは、GaAs系太陽電池で達成可能な高い光電変換性能を示している。GaAsの直接遷移型バンドギャップにより強力な光吸収が可能となり、またAlGaAsウィンドウ層によって表面での光損失や再結合損失の低減が実現される。さらに層の厚み、ドーピング濃度、および光学設計を最適化することで性能向上が見込まれる。例えば、反射防止コーティングを施すことでデバイスによる入射太陽光の吸収率を高めることが可能である。 反射防止コーティングのチュートリアル.

5. 光学シミュレーションを確認する

メインウィンドウ上部にある「光学」タブをクリックすると、光学関連のリボンが開きます（??）。なお、光強度（Suns）は1.0に設定されています。「伝達行列」をクリックすると、光学シミュレーションエディタが起動します。

エディタのツールバーでは、光生成プロファイルの算出方法を選択できる。具体的には伝達行列法、有限差分TMM法、またはより単純なプロファイル（指数関数型、平坦型、ファイルからの読み込み、QEスペクトルからの生成、一定値）から選択可能である。この例では、選択されているモデルが有限差分TMM法である。伝達行列法ではデバイスを複数の薄膜の積層体として扱い、各界面における反射や干渉、各層内での光吸収を考慮しながら各層内を伝播する光の場を計算する。電気モデルを再実行せずに光場を再計算するには「光シミュレーション実行」ボタンをクリックすればよい。背景情報としては伝達行列法とは何かおよびOghmaNanoにおける伝達行列法を参照されたい。??の「光子分布」タブでは波長および深さごとの光子密度が示され、??の「吸収された光子分布」タブでは光子がどの位置で吸収されたかが表示される。また「生成率」「反射光」「透過光」の各タブでは積分された生成プロファイルならびに反射・透過スペクトルが得られる。

これらのマップから明らかな特徴がいくつかある：

• 短波長の光は表面付近で吸収される。GaAsは青色光を非常に強く吸収するため、こうした光子は最初の数百ナノメートル以内、すなわち発電層およびp-n接合付近で吸収される。
• 波長の長い光はより深くまで侵入し、n型GaAs吸収層全体で吸収される。このため、吸収層は数マイクロメートルの厚さが必要となるのである。
• 870 nmを超える波長の光子は、デバイス全体を通過して吸収されることなく透過する。このような光子のエネルギーはGaAsのバンドギャップ値よりも低く、キャリアを生成することができない。
• 特定の波長域（例えば760 nm付近）における暗い水平線は、デバイス自体の特性によるものではなく、大気吸収によって生じたAM1.5Gスペクトルの減少部分に起因するものである。

GaAs層内で吸収されたすべての光子は電子・正孔対を生成する。この吸収された光子の密度が、ドリフト拡散方程式における生成率 \(G\) となる。AlGaAs製のウィンドウ層はGaAsよりもバンドギャップが広いため、両者のバンドギャップの間のエネルギーを持つ光子はこの層を透過し、その下にあるGaAs層で吸収される。こうして生成されたキャリアは効率よく回収される。一方、ウィンドウ層自体または未保護表面に極めて近いGaAs層で吸収された光子は、再結合によって失われる可能性がはるかに高くなる。

6. 電気的物質パラメータ

メインウィンドウに戻り、選択してください デバイス構造 タブを選択してクリックする 電気的特性値 左側のツールバーから選択すると、これが開く 電気的パラメータ編集機能 (??各半導体層ごとにタブが存在する。スクリーンショットにはアルミニウムが示されている_0.3ガリウム_0.7ウィンドウ層としては、一般的にAlが選択される。_xガリウム_1-xOghmaNanoにおけるLua材料モデルの表示例 コードエディタアルミニウムの含有率 \(x\) これは以下を通じて制御される 勾配／ひずみ編集機能 (??ツールバーでは電気物性の制御も行え、ドリフト拡散現象や平衡状態の設定が可能である。 ショックレー・リード・ホール（SRH） トラップ機能は有効化されているが、動的トラップ、オージェ再結合、および励起子輸送機能は無効化されている。後者のモデルは特に無秩序構造を持つ半導体デバイスにおいて有用だが、本シミュレーションでは必要とされない。

OghmaNanoでは、材料特性を直接指定することも、内蔵のmaterials databaseを用いて計算することも可能である。GaAs、AlGaAs、InGaAsといった特性がよく把握されている無機半導体に対しては、このデータベースがバンドギャップ、電子親和力、carrier mobility、再結合といった特性値を、組成や温度の関数として算出するためのモデルを提供している。これらのモデルはLua scriptsとして実装されており、利用可能な材料はLua material model libraryに一覧表示される。Scriptボタンを押下するとシミュレータは選択されたLua材料モデルを使用するが、押下しない場合はParametersタブに直接入力された値が用いられる。このより簡易な手法は、材料パラメータが個別に変化しうる無秩序な構造を持つ半導体や特性が十分に把握されていない半導体の解析に特に有用である。

これら2つの手法は組み合わせることも可能であり、その例がパラメータタブ（??）に示されている。ここでは、キャリア移動度、有効状態密度、バンドギャップ、電子親和力といった特性値は直接入力されるが、SRHトラップパラメータなどはLuaモデルから供給され、「値はスクリプト参照」と記載される。この動作を制御するには、Scriptボタンの横にある矢印をクリックし「材料データベースを上書き」を選択する（??）。チェックが入っているパラメータはテキストボックス内の値が使用され、チェックがないパラメータは材料データベースの値が用いられる。この機能は、ヘテロ接合におけるバンドオフセットの解析や、基礎となる材料モデルを変更せずに個々のパラメータを調整する際に有用である。同メニューには「材料スクリプトのインポート」機能もあり、データベースから別のLuaモデルを選択することも可能だ。

この例では、y方向における電子および正孔の移動度、バンドギャップ、電子親和力、有効状態密度、相対誘電率が上書き設定される一方で、再結合パラメータは引き続きLuaモデルから取得される。ヘテロ接合におけるバンド整合性を調査する際に有用であり、バンドギャップや電子親和力を変更することで、共通の材料データベースを修正することなく伝導帯および価電子帯のオフセットを変化させることが可能となる。また、この上書き機能はパラメータスイープや実験データとの比較にも利用できる。

スクリプトメニューには、データベースから別の材料モデルを選択するための「材料スクリプトのインポート」オプションも用意されている。詳細については、Luaを用いた電気的パラメータの定義、材料モデルライブラリ、および材料データベースの説明書を参照されたい。ユーザー独自の材料定義はoghma_local内に保存される。

7. エネルギー帯および擬フェルミ準位

半導体デバイスを理解するための最も有用なツールの一つがスナップショットビューアである。このビューアを使用すると、シミュレーションの各段階においてエネルギーバンド、キャリア密度、電界といった物理量をプロットすることが可能となる。例えばJ–Vスイープの際には、印加電圧の変化に伴うこれらの物理量の変動が確認でき、過渡シミュレーションでは時間経過に伴う変化を観察できる。ビューアを開くには、Outputタブに移動してsnapshotsフォルダをダブルクリックする。ウィンドウ下部にあるスライダーを操作することで、シミュレーションの各ステップ間を移動できる。詳細についてはスナップショットビューアの説明書を参照されたい。

興味深い出発点として、暗状態および照明下におけるデバイスを比較することが挙げられる。暗状態でのシミュレーションを行うには、Opticalリボンを選択し、Light intensity (Suns)の値をゼロに設定して(??)、再度シミュレーションを実行する。得られた約0 Vにおけるエネルギーバンド図が??に示されている。電子および正孔のフェルミ準位が一致しており、デバイス全体にわたって完全に平坦である点に注目されたい。これは熱平衡状態の特徴であり、半導体のバンドが大きく湾曲しているにもかかわらず、デバイス内を流れる正味の電流は存在しない。

さて、照明強度を1サンに戻してシミュレーションを再実行してください。得られた約0 Vにおけるバンド構造が??に示されています。伝導帯および価電子帯の形状は暗状態時とほぼ同様ですが、電子および正孔の擬フェルミ準位はもはや一致していません。これは、光生成された電子と正孔がデバイス内で非平衡状態のキャリア分布を形成したためです。

バンド構造図を詳細に見てみると、最も顕著な特徴は約0.25µmから0.55µmの範囲における著しいバンドの曲がり具合である。これはGaAsのp-n接合における空乏層に相当する。n型吸収層のドーピング濃度がp型発電層よりも低いため、空乏層は主に吸収層側へと広がっている。この結果生じる内部電界は、光生成された電子・正孔対の分離を促進し、電子をn型接触部へ、正孔をp型接触部へとそれぞれ輸送する役割を果たしている。

照射された表面側、すなわちAlGaAsウィンドウ層とGaAs発電層が接する部分においても興味深い現象が見られる。伝導帯における不連続性は、両材料の電子親和力およびバンドギャップの違いに起因する。このヘテロ接合によってエネルギー障壁が形成され、p型発電層内の少数キャリアである電子が表面まで到達して再結合するのを抑制する効果がある。またデバイスの裏面付近でも、高濃度にドーピングされたn+接触層への移行に対応するごく小さな変化が確認できる。

電子と正孔の擬フェルミ準位間の隔たりは、半導体が平衡状態からどの程度離れているかを示す指標となる。縮退していないキャリア統計においては、局所的な電子および正孔の数密度はこの準位間隔と以下の関係にある。

\[ np = n_i^2 \exp\left(\frac{E_{F,n}-E_{F,p}}{k_B T}\right) \]

照射光の作用により、光生成によってキャリア密度が増加し、\(E_{F,n}\)と\(E_{F,p}\)の間に分離が生じる。この分離は短絡状態、すなわち外部から印加される電圧がゼロの場合でもデバイス内部に存在し得る。開回路状態では、端子電圧は両電極界面で形成される擬フェルミ準位の分離度によって決定され、\(qV_{OC}\)は理想的な選択的接触面におけるその分離度を表している。再結合現象はこの分離度の大きさを制限するため、開回路電圧の決定において極めて重要な役割を果たす。スナップショットスライダーを用いて電圧掃引を行うことで、太陽電池が開回路電圧に近づく際にエネルギーバンドおよび擬フェルミ準位がどのように変化していくかを観察することが可能である。

9. さらなる検討事項

10. まとめおよび学習内容

このチュートリアルでは、広バンドギャップを持つAlGaAsウィンドウ層、p型GaAs発電層、n型GaAs吸収層、および高濃度ドープされた裏面電極から構成されるGaAs/Al0.3Ga0.7Asヘテロ接合太陽電池をシミュレーションした。1太陽光照射下において、シミュレーション結果は短絡電流密度26.91 mA cm⁻²、開放電圧0.922 V、フィルファクター80.8%、および電力変換効率20.05%を示した。これらの結果は、半導体ヘテロ構造の光学的・電気的特性がいかにしてその太陽光発電性能を決定づけるかを示している。

光生成分布やエネルギーバンド図を検討することで、太陽光がGaAs層内でどのように吸収されるか、内部電界が光生成キャリアをどのように分離するか、またAlGaAs/GaAsヘテロ接合がキャリア輸送や再結合に影響を与えるバンドオフセットをもたらすかを明らかにした。暗状態および照明下におけるデバイスの比較からは、擬フェルミ準位の分裂の重要性ならびにそれが開回路電圧とどのような関係にあるかも確認できた。

最後に、OghmaNanoがいかにして光学的な伝達行列計算とドリフト拡散モデリングを組み合わせているか、またそのLua材料データベースおよびパラメータ上書き機能によって半導体の物性がどのように制御可能かについても検討した。これらの機能が連携することで、層の厚さ、ドーピング、合金組成、および材料特性がデバイス性能に与える影響を調査するための枠組みが提供される。同様の手法は他のIII–V族太陽電池や半導体ヘテロ構造にも応用可能である。