PN 接合ダイオード IV 特性シミュレータ
プログラムの動作
pnjunction.py は、PN接合ダイオードの電流-電圧（IV）特性をシミュレートし、その結果をグラフで表示するPythonプログラムです。
このプログラムの主な機能は以下の通りです。
与えられた材料パラメータ（キャリア有効質量、移動度、寿命、ドナー/アクセプター濃度、バンド端エネルギーなど）と動作温度に基づいて、PN接合の内部電位 (\(V_{bi}\)) と逆飽和電流密度 (\(J_S\)) を計算します。
直列抵抗 (\(R_S\)) を考慮したショックレーダイオード方程式を用いて、指定された電圧範囲におけるダイオード電流を数値的に計算します。
計算された電流-電圧特性を半対数グラフとしてプロットし、ダイオードの電気的挙動を視覚的に表示します。
すべての物理的および設定パラメータはコマンドライン引数として指定可能であり、ユーザーはさまざまな条件でのダイオード特性を簡単に探索できます。
このプログラムは、半導体デバイスの基本的な電気的特性、特にPN接合ダイオードの順方向・逆方向バイアス時の挙動をシミュレーションを通じて理解することを目的としています。
原理
pnjunction.py は、PN接合ダイオードの挙動を記述する半導体物理学の基本原理に基づいています。
物理定数:
ボルツマン定数: \(k_B = 1.380649 \times 10^{-23} \, \text{J/K}\) (
KB)素電荷: \(q = 1.60217663 \times 10^{-19} \, \text{C}\) (
QE)
有効状態密度: 伝導帯有効状態密度 \(N_C\) と価電子帯有効状態密度 \(N_V\) は、温度 \(T\) と有効質量 \(m^*\) に依存します。 $\(N_C = 2.51 \times 10^{19} \cdot (m_{en}^{1.5}) \cdot \left(\frac{T}{300}\right)^{1.5}\)\( \)\(N_V = 2.51 \times 10^{19} \cdot (m_{hp}^{1.5}) \cdot \left(\frac{T}{300}\right)^{1.5}\)\( ここで、\)m_{en}\( は電子有効質量比、\)m_{hp}$ は正孔有効質量比です。
フェルミ準位: n型半導体の伝導帯端 \(E_{Cn}\) とドナー濃度 \(N_D\)、およびp型半導体の価電子帯端 \(E_{Vp}\) とアクセプター濃度 \(N_A\) を用いて、フェルミ準位 \(E_{Fn}\) および \(E_{Fp}\) を計算します。 $\(E_{Fn} = E_{Cn} - \frac{k_B T}{q} \ln\left(\frac{N_D}{N_C}\right)\)\( \)\(E_{Fp} = E_{Vp} + \frac{k_B T}{q} \ln\left(\frac{N_A}{N_V}\right)\)$
内部電位 (Built-in Potential): PN接合の内部電位 \(V_{bi}\) は、接合が熱平衡状態にあるときのフェルミ準位の差に等しいです。 $\(V_{bi} = E_{Fp} - E_{Fn}\)$
拡散係数と拡散長: アインシュタインの関係式により、キャリア移動度 \(\mu\) から拡散係数 \(D\) を求めます。 $\(D_n = \mu_n \frac{k_B T}{q}\)\( \)\(D_p = \mu_p \frac{k_B T}{q}\)\( 拡散長 \)L\( は、拡散係数 \)D\( とキャリア寿命 \)\tau\( から計算されます。 \)\(L_n = \sqrt{D_n \tau_n}\)\( \)\(L_p = \sqrt{D_p \tau_p}\)$
逆飽和電流密度 (Saturation Current Density): 理想的なPN接合における逆飽和電流密度 \(J_S\) は、少数キャリアの拡散によって生じます。 $\(J_S = q \left( \frac{D_n}{L_n} N_D e^{-\frac{q V_{bi}}{k_B T}} + \frac{D_p}{L_p} N_A e^{-\frac{q V_{bi}}{k_B T}} \right)\)\( ここで、\)N_D\( はn型側のドナー濃度、\)N_A$ はp型側のアクセプター濃度です。
ダイオード電流の計算: 直列抵抗 \(R_S\) を考慮したショックレーのダイオード方程式を用いて、ダイオードに流れる電流 \(I_D\) を計算します。 $\(I_D = A \cdot J_S \left( e^{V_D / (n \cdot V_T)} - 1 \right)\)\( ここで、\)A\( は接合面積、\)n\( は理想係数（本プログラムでは \)1.0\( を使用）、\)V_T = k_B T / q\( は熱電圧、\)V_D\( はダイオード両端にかかる電圧です。 外部印加電圧 \)V_{app}\( は、ダイオード電圧 \)V_D\( と直列抵抗 \)R_S\( による電圧降下の和として表されます。 \)\(V_{app} = V_D + I_D \cdot R_S\)\( この方程式は \)V_D\( について直接解くことが難しいため、プログラムではニュートン法のような数値解法（`calculate_diode_current` 関数内で `dv_dv` を用いた更新式）を使用して、与えられた \)V_{app}\( に対応する \)V_D\( と \)I_D$ を反復的に計算します。
必要な非標準ライブラリとインストール方法
本プログラムは、以下の非標準ライブラリを使用しています。
numpy: 数値計算（特に配列操作）に使用されます。
matplotlib: グラフの描画に使用されます。
これらのライブラリは、Pythonのパッケージマネージャー pip を使ってインストールできます。コマンドラインで以下のコマンドを実行してください。
pip install numpy matplotlib
必要な入力ファイル
このプログラムは、外部の入力ファイルを必要としません。 すべての必要なパラメータは、プログラム実行時にコマンドライン引数として直接指定されます。
生成される出力ファイル
このプログラムは、いかなるファイルも生成しません。 計算結果は以下の形式で出力されます。
標準出力:
計算された内部電位 (\(V_{bi}\)) の値。
計算された逆飽和電流密度 (\(J_S\)) の値。
グラフィカルウィンドウ:
印加電圧と電流の絶対値の半対数プロット（IV特性曲線）。このウィンドウは、
--noplotオプションが指定されていない場合に表示されます。
コマンドラインでの使用例 (Usage)
pnjunction.py は、以下の書式でコマンドラインから実行します。
python pnjunction.py [OPTIONS]
利用可能なオプション（引数）は以下の通りです。
オプション |
型 |
デフォルト値 |
説明 |
|---|---|---|---|
|
float |
|
温度 [K] |
|
float |
|
接合面積 [cm\(^2\)] |
|
float |
|
直列抵抗 [ohm] |
|
float |
|
n型側ドナー濃度 [cm\(^{-3}\)] |
|
float |
|
n型側伝導帯底エネルギー [eV] |
|
float |
|
n型側電子有効質量比 |
|
float |
|
n型側電子移動度 [cm\(^2\)/Vs] |
|
float |
|
n型側電子寿命 [s] |
|
float |
|
p型側アクセプター濃度 [cm\(^{-3}\)] |
|
float |
|
p型側価電子帯頂上エネルギー [eV] |
|
float |
|
p型側正孔有効質量比 |
|
float |
|
p型側正孔移動度 [cm\(^2\)/Vs] |
|
float |
|
p型側正孔寿命 [s] |
|
float |
|
印加電圧開始点 [V] |
|
float |
|
印加電圧終了点 [V] |
|
float |
|
印加電圧ステップサイズ [V] |
|
action |
|
プロット表示を無効にする |
コマンドラインでの具体的な使用例
1. デフォルトパラメータでの実行
最も基本的な実行方法です。すべてのパラメータはデフォルト値が使用されます。
python pnjunction.py
実行結果の説明: 標準出力には、計算された内部電位 (\(V_{bi}\)) と逆飽和電流密度 (\(J_S\)) が表示されます。
Calculated Vbi: 0.8143 V
Calculated Js : 1.348624e-12 A/cm2
また、PN Junction IV Characteristics というタイトルのグラフウィンドウが開き、-1.0V から 1.0V までの印加電圧に対するダイオード電流の絶対値（半対数スケール）がプロットされます。デフォルトの直列抵抗 \(R_S=10.0 \, \text{ohm}\) の曲線が表示されます。順方向バイアスでは、ある閾値電圧を超えると電流が指数関数的に増加し、高電圧側では直列抵抗によって電流の増加が制限される様子が確認できます。逆方向バイアスでは、電流はほぼ一定の低い値（飽和電流）を示します。
2. 特定のパラメータを変更して実行
温度を350Kに上げ、直列抵抗を5.0 ohm、接合面積を0.05 cm\(^2\) に変更して実行します。
python pnjunction.py --temp 350 --rs 5.0 --area 0.05
実行結果の説明: 標準出力の \(V_{bi}\) と \(J_S\) の値は、温度上昇と面積増加の影響で変化します。
Calculated Vbi: 0.7712 V
Calculated Js : 1.139454e-10 A/cm2
\(V_{bi}\) はわずかに減少し、\(J_S\) は温度と面積の増加により顕著に増加していることがわかります。 グラフウィンドウに表示される曲線は、これらのパラメータ変更を反映します。温度上昇により逆飽和電流が増加し、ダイオードの「オン」になる電圧（順方向電圧降下）がわずかに減少する傾向が見られます。直列抵抗が小さくなったため、高電流領域での電圧降下が少なくなり、同じ電圧でより多くの電流が流れるようになります。面積の増加も電流値全体を増加させます。
3. プロットを無効にして数値結果のみを出力
グラフ表示をせずに、計算された \(V_{bi}\) と \(J_S\) の値のみを確認したい場合に便利です。
python pnjunction.py --noplot
実行結果の説明: グラフウィンドウは表示されず、標準出力に以下の情報のみが表示されます。
Calculated Vbi: 0.8143 V
Calculated Js : 1.348624e-12 A/cm2
これは、グラフが不要な場合や、他のプログラムで出力を処理したい場合に役立ちます。