ペロブスカイト太陽電池におけるCELIV測定: パートA
1. 概要
このチュートリアルではまず、線形電圧増加による電荷抽出 （Charge Extraction by Linearly Increasing Voltage, CELIV）技術の基本概念を紹介し、 それがどのようにして薄膜半導体中のキャリア移動度と密度を調べるために使用されるかを説明します。 理論的背景を確立した後、実践的な手順へと進み、 OghmaNano内でCELIVシミュレーションを設定して実行する方法を示します。 この理論と実践的シミュレーションの組み合わせにより、 手法の明確な理解と、それを自分のデバイスに適用するためのツールの両方を得ることができます。
2. CELIVの背景
線形電圧増加による電荷抽出（CELIV）技術は、 有機太陽電池およびペロブスカイト太陽電池における 電荷キャリアの移動度と密度を研究するために用いられる実験手法です。 典型的なセットアップでは、太陽電池は信号発生器に接続され、 線形に増加する逆バイアスが印加されます。 その間、結果として生じる電流トランジェントがオシロスコープで記録されます。 電荷キャリアを生成するために、 連続照明または短いレーザーパルスが追加される場合もあります。 この構成の模式図は ?? に示されています。
電圧ランプが印加されると、測定される電流は2つの異なる寄与から構成されます。 1つはデバイスの幾何学的静電容量による矩形ベースライン、 もう1つは半導体中の移動可能な電荷キャリアの抽出による追加ピークです。 幾何学的静電容量項は常に存在しますが、 キャリア関連ピークは移動度とキャリア密度に関する直接的な情報を提供します。 これらの電圧および電流トランジェントの例は ?? に示されています。
キャリア移動度はCELIVトランジェントから直接決定できます。 抽出電流が最大値に達する時間を測定することで求められます。 均一なキャリア分布および再結合が無視できる最も単純な場合、 移動度は次のCELIV式で与えられます: \[ \mu = \frac{2 d^{2}}{3 A t_{\text{max}}^{2}} \] ここで \(d\) は活性層厚さ、 \(A\) は電圧ランプレート（\(A = \mathrm{d}V/\mathrm{d}t\)）、 \(t_{\text{max}}\) は電流ピークが発生する時間です。 この関係は実験で観測されるトランジェントを、 関心のある基本輸送パラメータと結び付けます。
3. CELIVの限界
CELIV移動度式の導出はいくつかの単純化仮定に依存しています。 第一に、輸送を支配する電荷キャリアの種類は1種類のみであると仮定します。 すなわち電圧ランプが印加されたときに移動するキャリア種は1つだけです。 第二に、キャリアは一様に分布しており、 ?? の左側に示された理想化ケースのように、 デバイスから滑らかに外へ移動すると仮定します。 第三に、抽出中の再結合過程は無視できるとされ、 顕著なキャリア損失は起こらないと仮定されます。 実際には、これらの条件が満たされることはほとんどありません。 ?? の右側に示すように、 電荷キャリアはしばしば浅い・中程度・深いトラップを占有し、 その遅れた放出により抽出トランジェントは広がったり歪んだりします。 多くの実験技術と同様に、 CELIVで抽出された移動度は真の微視的値ではなく近似値として扱うべきであり、 実際には見かけの移動度はトランジェント中に変化する場合があります （doi:10.1063/1.4818267）。 この制限は特に無秩序有機半導体において重要であり、 エネルギー無秩序とトラッピングが輸送に強く影響します。 対照的にペロブスカイト材料では、 より高い移動度と低いトラップ密度を示すことが多く、 標準的なCELIV解析はより堅牢で解釈しやすくなります。
💡 理解度チェック: 簡単な例にCELIV理論を適用してみましょう。
厚さ d = 200 nm のデバイスを考えます。
電圧ランプレートは A = 2 × 106 V/s、
電流ピークは tmax = 5 µs に観測されました。
CELIV移動度式を用いてキャリア移動度を推定してください。
答えを表示
CELIV移動度は次式で与えられます
\[
\mu = \frac{2 d^{2}}{3 A t_{\text{max}}^{2}}
\]
値を代入すると:
d = 200 × 10-9 m,
A = 2 × 106 V/s,
tmax = 5 × 10-6 s.
\[ \mu = \frac{2 (200 × 10^{-9})^{2}}{3 (2 × 10^{6}) (5 × 10^{-6})^{2}} \approx 1.1 × 10^{-4} \; \text{cm}^2/\text{Vs} \]
この簡単な例は、CELIVトランジェントのピーク時間から 移動度推定値を求められることを示しています。
4: OghmaNanoでCELIVシミュレーションを実行する
New simulation をクリックします。 すると利用可能なデバイスタイプのライブラリが開きます （??）。 Perovskite cells（赤で強調表示）をダブルクリックして ペロブスカイト例フォルダを開きます。 MAPbI₃ device、 Perovskite solar cell、 そして専用CELIVデモなどの事前設定シミュレーションが表示されます （??）。 このチュートリアルでは Perovskite solar cell – CELIV example を選択してください。 プロンプトが表示されたら、 書き込み権限のあるフォルダにシミュレーションを保存します。
Perovskite solar cell – CELIV example を選択すると、 メインシミュレーションウィンドウが開きます （??）。 計算を開始するには、 Run simulation（青い再生アイコン）をクリックするか F9 を押します。 OghmaNanoは時間依存ドリフト–拡散方程式を解き、 CELIVトランジェントを生成します。
jv.csv（JVカーブデータ）、
time_j.csv（抽出電流 vs 時間）、
time_v.csv（印加電圧 vs 時間）、
optical_output（電場分布）が含まれます。
任意のファイルをダブルクリックすると内蔵プロットツールで結果を表示できます。
ステップ5: CELIV結果を表示する
Output タブ
（??）
を開き、
time_v.csv をダブルクリックします。
これにより印加電圧プログラム —
CELIVランプ —
が
??
に示すようにプロットされます。
次に
time_j.csv をダブルクリックして
抽出電流トランジェントを表示します
（??）。
トランジェントは典型的なCELIV形状ですが、
**反転して表示**されることに気づくでしょう。
OghmaNanoの符号規約では、
デバイスから外へ流れる電流は負です。
??
のランプが印加されると、
電荷がデバイスから引き出されるため、
デバイスは負電流（抽出ピーク）を出力します。
後の反対符号ピーク（
??
では正）は、
バイアスが戻り、
トランジェントが初期状態へ戻った後に
電荷が再びデバイスへ流入する際に発生します。
?? のプロットは、 シミュレーション電流を −1 倍して 実験で一般的に用いられる表示形式に合わせたCELIVトランジェントを示しています。 この形式では、 ベースラインはデバイスの幾何学的静電容量に対応し、 鋭い正ピークは印加電圧ランプ中に移動キャリアが抽出されることを示します。 このピークの大きさとタイミングが、 キャリア移動度と密度を決定するために用いられる主要量です。 このような従来の向きでトランジェントを表示することで、 OghmaNanoシミュレーションを 公開されたCELIV測定結果と直接比較しやすくなります。 移動度を抽出するには、 電流が最大値に達する時間 （\(t_\text{max}\)） を特定し、 この値をCELIV移動度式 \(\mu = \tfrac{2d^2}{3At_\text{max}^2}\) に代入します。 ここで \(d\) は活性層厚さ、 \(A\) は電圧ランプレートです。
div style="clear: both;">🧪 計算例:
??
からCELIVピーク時間を読み取り、
移動度を推定してください。
活性層厚さ d = 600 nm、
電圧ランプレート
A = 5 × 105 V/s、
測定ピーク時間
tmax = 3 µs
と仮定します。
どの移動度が得られますか？
答えを表示
CELIV関係式 \[ \mu = \frac{2 d^{2}}{3 A t_{\text{max}}^{2}} \] を使用し、 \(d = 600\times10^{-9}\,\mathrm{m}\), \(A = 5\times10^{5}\,\mathrm{V\,s^{-1}}\), \(t_{\text{max}} = 3\times10^{-6}\,\mathrm{s}\) を代入します。
\[ \mu = \frac{2(600\times10^{-9})^{2}}{3\,(5\times10^{5})\,(3\times10^{-6})^{2}} \approx 5.3\times10^{-8}\ \mathrm{m^{2}\,V^{-1}\,s^{-1}} \] \(\mathrm{cm^{2}\,V^{-1}\,s^{-1}}\) へ変換すると （\(1\ \mathrm{m^{2}} = 10^{4}\ \mathrm{cm^{2}}\)）: \[ \mu \approx 5.3\times10^{-4}\ \mathrm{cm^{2}\,V^{-1}\,s^{-1}}. \]
この値は無秩序有機半導体で報告される典型的範囲にあります。 （ランプレートや照明強度で \({t_\text{max}}\) が変化すると、 \(\mu\) もそれに応じて変化します。）
👉 次のステップ: 次のCELIV（C-LIV）シミュレーションセクションへ進みましょう — パートB — ここでは次の内容を行います: 電圧ランプレートの調整、 照明強度の実験、 電子／正孔移動度の変更が CELIVトランジェント形状をどのように変化させるかを探索します。