pnjunction.py ダウンロード/コピー

pnjunction.py をダウンロード

pnjunction.py
pnjunction.py
  1import numpy as np
  2import math
  3import argparse
  4import matplotlib.pyplot as plt
  5
  6"""
  7pnjunction.py
  8
  9PN接合ダイオードの電流-電圧 (IV) 特性を計算し、プロットするスクリプト。
 10
 11このスクリプトは、ユーザーが指定した半導体材料のパラメータと環境条件（温度）に基づき、
 12PN接合ダイオードの重要な電気的特性である逆方向飽和電流 (Js) と内部電位 (Vbi) を計算します。
 13その後、指定された印加電圧範囲におけるダイオードの電流-電圧 (IV) 特性を、
 14ニュートン法を用いた数値計算によりシミュレーションします。
 15計算されたIV特性はmatplotlibライブラリを使用してグラフとして視覚的に表示され、
 16主要な計算結果はコンソールに出力されます。
 17
 18関連リンク:
 19:doc:`pnjunction_usage`
 20"""
 21
 22# 物理定数
 23KB = 1.380649e-23
 24QE = 1.60217663e-19
 25
 26def calculate_material_params(args):
 27    """
 28    PN接合ダイオードの材料パラメータ（逆方向飽和電流、内部電位）を計算します。
 29
 30    半導体の特性（キャリア有効質量、ドーピング濃度、移動度、寿命など）と温度に基づき、
 31    伝導帯・価電子帯有効状態密度、フェルミ準位、内部電位 (Vbi)、キャリア拡散係数、拡散長を計算し、
 32    最終的にダイオードの逆方向飽和電流 Js を決定します。
 33
 34    :param args: コマンドライン引数を格納したargparse.Namespaceオブジェクト。
 35                 以下の属性を持つことが期待されます。
 36                 :param args.temp: float 温度 (K)
 37                 :param args.men: float 電子の有効質量比 (無次元)
 38                 :param args.mhp: float 正孔の有効質量比 (無次元)
 39                 :param args.ecn: float N側伝導帯端エネルギー (eV)
 40                 :param args.ndn: float N側ドーピング濃度 (/cm^3)
 41                 :param args.evp: float P側価電子帯端エネルギー (eV)
 42                 :param args.nap: float P側ドーピング濃度 (/cm^3)
 43                 :param args.mun: float 電子移動度 (cm^2/Vs)
 44                 :param args.taun: float 電子寿命 (s)
 45                 :param args.mup: float 正孔移動度 (cm^2/Vs)
 46                 :param args.taup: float 正孔寿命 (s)
 47    :returns:
 48        :returns js: float 逆方向飽和電流 (A/cm^2)
 49        :returns v_bi: float 内部電位 (V)
 50    """
 51    T = args.temp
 52    nc = 2.51e19 * (args.men**1.5) * ((T/300)**1.5)
 53    nv = 2.51e19 * (args.mhp**1.5) * ((T/300)**1.5)
 54    
 55    efn = args.ecn - (KB * T / QE) * math.log(args.ndn / nc)
 56    efp = args.evp + (KB * T / QE) * math.log(args.nap / nv)
 57    v_bi = efp - efn
 58
 59    dn = args.mun * (KB * T / QE)
 60    dp = args.mup * (KB * T / QE)
 61    ln, lp = math.sqrt(dn * args.taun), math.sqrt(dp * args.taup)
 62
 63    term_exp = math.exp(-QE * v_bi / (KB * T))
 64    js = QE * (dn / (ln)) * args.ndn * term_exp + QE * (dp / (lp)) * args.nap * term_exp
 65    return js, v_bi
 66
 67def calculate_diode_current(v_target, temp, js, n_diode, rs, area, v_initial):
 68    """
 69    指定された印加電圧におけるPN接合ダイオードの電流をニュートン法を用いて計算します。
 70
 71    理想ダイオードモデルに直列抵抗 (Rs) を加味した方程式 (V_total = V_diode + I_diode * R_s) を
 72    ニュートン法を用いて数値的に解きます。これにより、目標の印加電圧 V_target における
 73    ダイオードにかかる電圧 V_diode と、それに対応する電流 I_diode を高精度で求めます。
 74    ニュートン法の初期値として前回の計算結果を使用することで、収束性を向上させています。
 75
 76    :param v_target: float 目標とする印加電圧 (V)。
 77    :param temp: float 温度 (K)。
 78    :param js: float 逆方向飽和電流 (A/cm^2)。
 79    :param n_diode: float ダイオード理想係数 (無次元)。
 80    :param rs: float 直列抵抗 (ohm)。
 81    :param area: float ダイオードの面積 (cm^2)。
 82    :param v_initial: float ニュートン法の初期値として使用するダイオード電圧 (V)。
 83                      通常は前回の計算結果や予測値が用いられます。
 84    :returns:
 85        :returns v_diode: float ダイオードにかかる最終的な電圧 (V)。
 86        :returns current: float ダイオードに流れる最終的な電流 (A)。
 87    """
 88    if v_target == 0: return 0.0, 0.0
 89    v_diode = v_initial
 90    dv, eps, max_iter = 1e-4, 1e-7, 100
 91    vt = (n_diode * KB * temp) / QE
 92
 93    for _ in range(max_iter):
 94        get_i = lambda vd: js * (math.exp(vd / vt) - 1.0) * area
 95        im, ip = get_i(v_diode - dv), get_i(v_diode + dv)
 96        v_total = 0.5 * ((v_diode - dv + im * rs) + (v_diode + dv + ip * rs))
 97        
 98        if abs(v_total - v_target) < eps:
 99            return v_diode, 0.5 * (im + ip)
100        
101        dv_dv = ((v_diode + dv + ip * rs) - (v_diode - dv + im * rs)) / (2.0 * dv)
102        v_diode += (v_target - v_total) / dv_dv
103    return v_diode, js * (math.exp(v_diode / vt) - 1.0) * area
104
105def main():
106    """
107    PN接合ダイオードのIV特性シミュレーションのメイン処理を実行します。
108
109    コマンドライン引数から半導体材料のパラメータ、シミュレーションの電圧範囲、
110    プロットの表示オプションなどを受け取ります。
111    `calculate_material_params` 関数でダイオードの基本パラメータ (Js, Vbi) を計算し、
112    その後、`calculate_diode_current` 関数を繰り返し呼び出して、
113    指定された電圧範囲におけるダイオード電流を計算します。
114    計算された電流-電圧特性はmatplotlibでプロットされ、コンソールに主要な計算結果が出力されます。
115
116    :returns: なし
117    :rtype: None
118    """
119    parser = argparse.ArgumentParser(description='PN Junction IV Characteristics Simulator')
120    parser.add_argument('--temp', type=float, default=300.0, help='Temperature [K]')
121    parser.add_argument('--area', type=float, default=0.01, help='Diode Area [cm^2]')
122    parser.add_argument('--rs', type=float, default=10.0, help='Series Resistance [ohm]')
123    parser.add_argument('--ndn', type=float, default=1.0e16, help='N-side doping concentration [/cm^3]')
124    parser.add_argument('--ecn', type=float, default=4.05, help='N-side conduction band edge energy [eV]')
125    parser.add_argument('--men', type=float, default=0.19, help='Effective mass ratio of electron (N-side)')
126    parser.add_argument('--mun', type=float, default=1500.0, help='Electron mobility (N-side) [cm^2/Vs]')
127    parser.add_argument('--taun', type=float, default=1e-5, help='Electron lifetime (N-side) [s]')
128    parser.add_argument('--nap', type=float, default=1.0e16, help='P-side doping concentration [/cm^3]')
129    parser.add_argument('--evp', type=float, default=5.17, help='P-side valence band edge energy [eV]')
130    parser.add_argument('--mhp', type=float, default=0.16, help='Effective mass ratio of hole (P-side)')
131    parser.add_argument('--mup', type=float, default=500.0, help='Hole mobility (P-side) [cm^2/Vs]')
132    parser.add_argument('--taup', type=float, default=1e-5, help='Hole lifetime (P-side) [s]')
133    parser.add_argument('--v0', type=float, default=-1.0, help='Start voltage for simulation [V]')
134    parser.add_argument('--v1', type=float, default=1.0, help='End voltage for simulation [V]')
135    parser.add_argument('--step', type=float, default=0.02, help='Voltage step for simulation [V]')
136    parser.add_argument('--noplot', action='store_true', help='Disable plot display')
137    args = parser.parse_args()
138
139    js, v_bi = calculate_material_params(args)
140    
141    v_apps = np.arange(args.v0, args.v1 + args.step, args.step)
142    currents = []
143    
144    v_prev = args.v0 * 0.2
145    for v_app in v_apps:
146        _, i = calculate_diode_current(v_app, args.temp, js, 1.0, args.rs, args.area, v_prev)
147        currents.append(i)
148        v_prev = _
149
150    abs_currents = np.abs(currents)
151
152    # グラフ描画
153    if not args.noplot:
154        plt.figure(figsize=(8, 6))
155        plt.semilogy(v_apps, abs_currents, 'b-', label=f'Rs = {args.rs} ohm')
156        plt.title("PN Junction IV Characteristics")
157        plt.xlabel("Applied Voltage (V)")
158        plt.ylabel("|Current| (A)")
159        plt.grid(True, which="both", ls="-", alpha=0.5)
160        plt.legend()
161        
162        # 0Vの線を表示
163        plt.axvline(0, color='black', lw=1)
164        
165        print(f"Calculated Vbi: {v_bi:.4f} V")
166        print(f"Calculated Js : {js:.6e} A/cm2")
167        plt.show()
168
169if __name__ == "__main__":
170    main()