Shockley方程式 $I = I_0(e^{V/nV_T}-1)$ を使ってダイオードの電流-電圧特性を計算。材料・温度・理想係数をリアルタイムで変えて半導体の動作原理を体験しよう。
k = 1.381×10⁻²³ J/K、q = 1.602×10⁻¹⁹ C。nは理想係数（1=拡散支配、2=再結合支配）。
PN接合ダイオードの電流-電圧特性を記述する基本的な式がショックレー方程式です。これは拡散電流モデルに基づいています。
$$I = I_0\left(e^{V/(nV_T)}- 1\right)$$$I$: ダイオードを流れる電流 (A)
$I_0$: 逆方向飽和電流 (A)。材料や温度で決まる定数。
$V$: ダイオードにかかる電圧 (V)。順方向を正とする。
$n$: 理想係数 (1〜2)。1は理想ダイオード（拡散支配）、2は再結合電流の影響が強い場合。
$V_T$: 熱電圧 (V)。次の式で定義されます。
熱電圧 $V_T$ は絶対温度 $T$ に比例する物理定数で、電子の熱エネルギーを電圧の次元で表したものです。
$$V_T = \frac{kT}{q}$$$k$: ボルツマン定数 (1.381×10⁻²³ J/K)
$T$: 絶対温度 (K)
$q$: 素電荷 (1.602×10⁻¹⁹ C)
室温 (300K) では約 26mV となります。この値が、電流が指数関数的に増加し始める電圧の目安になります。
整流回路（AC/DC変換）：交流を直流に変換する最も基本的な用途です。家庭用電源アダプターや電源回路の入力段で、シリコンダイオードがブリッジ整流子として使われています。シミュレーターで見られる順方向電圧（シリコンで約0.7V）は、この時の損失電圧となります。
電圧リファレンス・温度センサー：ダイオードの順方向電圧が温度に依存する特性を利用します。例えば、定電流を流した時の電圧変化を測ることで、精度の高い温度測定が可能です。シミュレーターで温度を変えた時のグラフのシフトが、この原理です。
高周波・無線通信（ゲルマニウムダイオード）：ゲルマニウムは順方向電圧が低いため、微小な高周波信号の検波（復調）に適しています。ラジオの鉱石検波器や一部のRF回路では、この特性が活かされています。
太陽電池：太陽電池の動作原理は、光によって生成されたキャリアがPN接合を駆動する「光起電力効果」ですが、その内部モデルはダイオード特性が基礎となっています。暗状態の太陽電池のI-V特性は、まさにこのシミュレーターのグラフそのものと言えます。
このシミュレーターを使う上で、特に実務設計を意識した注意点がいくつかあります。まず、「順方向電圧Vfは常に0.7V（シリコン）と思い込む」こと。これは大きな誤解です。データシートの0.7Vは、特定の電流値（例えば10mAや100mA）での「代表値」に過ぎません。シミュレーターで電流スケールを対数から線形に切り替え、1mAと100mAでの電圧降下を比べてみてください。電流が10倍になっても、電圧は約60mV（n=1, 室温の場合）しか増えませんが、絶対値は大きく変わります。低電流回路では0.4V程度、大電流整流では1Vを超えることも珍しくありません。
次に、逆方向飽和電流I₀の扱い。シミュレーターでは固定値ですが、実際の部品では製造バラツキが大きく、データシートには「最大値」しか載っていないことがほとんどです。例えば、あるシリコンダイオードのI₀が25°Cで1nA（典型）と書かれていても、最大値は50nAかもしれず、高温ではこの値がさらに1000倍以上に膨れ上がります。逆方向のリーク電流を厳密に考慮する必要がある高インピーダンス回路では、このバラツキが設計の成否を分けます。
最後に、ショックレー方程式の「適用限界」を理解すること。この式は中電流領域の挙動を非常によく説明しますが、現実のダイオードにはこのモデルに含まれない効果がたくさんあります。順方向の大電流領域では、半導体内部の抵抗（体抵抗）による電圧降下が無視できなくなり、グラフの傾きが急になります。逆に、極めて低い順方向電圧や逆方向領域では、表面リークや再結合の影響が支配的で、シミュレーションと実測が乖離します。あくまで「第一近似の美しいモデル」として使い、必要に応じてより高次のモデル（SPICEモデルなど）に進むことが大切です。
このI-V特性の理解は、PN接合を基礎とするあらゆる半導体デバイスの解析に直結します。最も直接的なのは、バイポーラトランジスタ(BJT)です。BJTのベース-エミッタ間はまさにPN接合ダイオードそのもので、その特性が増幅動作の根幹を成します。シミュレーターで温度を上げるとグラフが左にズレる現象は、BJT回路の動作点が温度で変動する「熱暴走」の原因の一つとして、実際のアンプ設計で常に考慮されます。
また、太陽電池（光起電力デバイス）の動作原理も、この方程式を拡張することで理解できます。太陽電池は、光によって生成されたキャリアがPN接合を駆動する「発電モード」のダイオードとみなせます。その出力特性は、暗状態のI-V曲線（まさにこのシミュレーターのグラフ）を下方向にシフトさせた形でモデル化されます。ここで学ぶ逆方向飽和電流I₀は、太陽電池の性能指標の一つである「ダイオード品質因子」と深く関係しています。
さらに、現代のマイクロエレクトロニクスの中核であるCMOS集積回路においても、寄生ダイオードの挙動は重要です。例えば、ICの入力保護回路や、ラッチアップ現象の原因となる寄生サイリスタ構造の中には、PN接合が関与しています。これらの不要な現象をシミュレーションで予測・排除するためには、基礎となるダイオード特性の理解が不可欠です。
ショックレー方程式の美しさに触れたら、次はその「導出過程」を追ってみることを強くお勧めします。教科書を開き、「少数キャリアの連続の方程式」と「ボルツマン関係式」というキーワードを追いかけてください。これらを組み合わせることで、なぜ電流が指数関数の形で表されるのか、その物理的イメージ（例えば、電圧が熱電圧V_Tの数倍かかると、境界のキャリア濃度が指数関数的に変化する）が腹落ちします。数学的には、微分方程式の境界値問題を解くことになります。
実用的な次のステップは、SPICEなどの回路シミュレーターで実際のダイオードモデルを触ってみることです。SPICEのダイオードモデルは、ショックレー方程式をベースにしつつ、体抵抗Rsや遷移キャパシタンスCj0など、先述した「適用限界」を補う多数のパラメータを追加したものです。例えば、.model文の中の「N」（理想係数）、「Is」（飽和電流）、「Rs」などをいじりながら過渡解析を実行すると、スイッチング時のリカバリー特性や電圧オーバーシュートがどう変わるかを体感でき、理論と実践が結びつきます。
最後に、材料を「シリコン」「ゲルマニウム」から一歩進めて、化合物半導体（ガリウムヒ素GaAsやシリコンカーバイドSiC）の世界に目を向けてみましょう。SiCのバンドギャップはシリコンの約3倍も大きいため、高温・高耐圧動作が可能で、電気自動車のインバーターなどに不可欠です。バンドギャップが変わるとI₀や順方向電圧がどう変わり、それが電力損失や発熱にどう影響するか。このツールで学んだ基礎概念は、こうした最先端のパワーデバイスを理解するための強固な土台となります。