The Shockley Diode Equation: Worked Examples and Real-World Limits

The Shockley diode equation (often called the ideal diode equation) is the fundamental mathematical model for the current-voltage (I-V) characteristic of a PN junction. If you need to predict how a silicon or germanium diode will behave under forward bias or slight reverse bias, this is the formula you use. Below is the exact formula, followed by a complete symbol breakdown, algebraic rearrangements for bench debugging, and fully worked examples with strict unit tracking.

The Core Formula and Symbol Definitions

The standard form of the diode equation calculates the diode current (I) based on the voltage dropped across the junction (V_D):

I = I_S · (e^{V_D / (n · V_T)} - 1)

Where V_T is the thermal voltage, defined as V_T = kT / q. To use this equation without blowing up your calculator, you must understand the physical meaning and typical magnitude of every variable. According to standard semiconductor physics references like Georgia State University's HyperPhysics, these parameters are highly temperature-dependent.

Symbol	Parameter	Typical Silicon Value	SI Unit
I	Diode forward current	1 mA to 1 A	Amperes (A)
IS	Reverse saturation (scale) current	10-12 to 10-9 A	Amperes (A)
VD	Voltage across the diode junction	0.5 V to 0.8 V	Volts (V)
n	Ideality factor (emission coefficient)	1.0 to 2.0	Dimensionless
VT	Thermal voltage (kT/q)	~25.85 mV at 300 K	Volts (V)
k	Boltzmann constant	1.3806 × 10-23	Joules/Kelvin (J/K)
T	Absolute junction temperature	298 K to 350 K	Kelvin (K)
q	Elementary charge	1.602 × 10-19	Coulombs (C)

Rearranged Forms for Circuit Debugging

On the bench, you rarely know all variables. You might measure a voltage drop and need to find the current, or measure a current and need to back-calculate the ideality factor of an unknown diode. Here are the algebraically rearranged forms:

• Solving for Junction Voltage (V_D):
  VD = n · VT · ln(I / IS + 1)
  Use when: You know the bias current and want to predict the forward voltage drop.
• Solving for Saturation Current (I_S):
  IS = I / (eVD / (n · VT) - 1)
  Use when: Characterizing a specific diode batch on a curve tracer to build a custom SPICE model.
• Solving for Ideality Factor (n):
  n = VD / (VT · ln(I / IS + 1))
  Use when: Determining if recombination current (n ≈ 2) or diffusion current (n ≈ 1) dominates the junction.

Worked Examples with Unit Tracking

The most common reason the diode equation yields nonsense (like a forward voltage of 45 volts) is a unit mismatch in the exponential term. The exponent V_D / (n · V_T) must be strictly dimensionless. Below are two step-by-step derivations.

Example 1: Calculating Forward Voltage (V_D)

Scenario: A generic silicon signal diode is conducting 15 mA at a junction temperature of 300 K. The manufacturer specifies an ideality factor n = 1.2 and a reverse saturation current I_S = 2.0 nA. Find V_D.

1. Convert all inputs to base SI units:
   I = 15 mA = 0.015 A
   I_S = 2.0 nA = 2.0 × 10^-9 A
   T = 300 K
2. Calculate Thermal Voltage (V_T):
   V_T = (k · T) / q
   V_T = (1.3806 × 10^-23 J/K · 300 K) / (1.602 × 10^-19 C)
   V_T = 0.02585 J/C. Since 1 Joule/Coulomb = 1 Volt, V_T = 0.02585 V (or 25.85 mV).
3. Apply the rearranged V_D formula:
   V_D = n · V_T · ln(I / I_S + 1)
   V_D = 1.2 · 0.02585 V · ln(0.015 A / 2.0 × 10^-9 A + 1)
4. Solve the inner term (units cancel):
   0.015 / 2.0e-9 = 7,500,000 (dimensionless)
   ln(7,500,000 + 1) ≈ 15.83
5. Final Multiplication:
   V_D = 0.03102 V · 15.83 = 0.491 V

Bench Reality Check: If you measure a real 1N4148 at 15 mA, your multimeter will likely read closer to 0.75 V. Why the discrepancy? The ideal Shockley equation ignores bulk series resistance (R_S). At 15 mA, the voltage drop across the silicon bulk and the wire leads adds roughly 0.25 V to the ideal junction voltage.

Example 2: Calculating Forward Current (I)

Scenario: A 1N4007 rectifier has 0.65 V dropped strictly across its PN junction. I_S = 50 nA, n = 1.5, and ambient temperature is 25°C (298.15 K). Find I.

1. Convert to SI and calculate V_T:
   V_D = 0.65 V
   I_S = 50 × 10^-9 A
   V_T = (1.3806e-23 · 298.15) / 1.602e-19 = 0.02569 V
2. Calculate the exponent denominator:
   n · V_T = 1.5 · 0.02569 V = 0.038535 V
3. Calculate the dimensionless exponent:
   V_D / (n · V_T) = 0.65 V / 0.038535 V = 16.867
4. Evaluate the exponential term:
   e^16.867 ≈ 21,144,534
   Subtract 1: 21,144,533
5. Multiply by I_S:
   I = 50 × 10^-9 A · 21,144,533 = 0.001057 A
   I = 1.057 mA

Assumptions, Limits, and Unit Traps

The Shockley equation is an elegant model, but it is not a universal law. As detailed in semiconductor physics literature, it relies on strict assumptions that break down in real-world power electronics.

When the Formula Applies (and When it Doesn't)

• Low-Level Injection: The equation assumes the injected minority carrier concentration is much smaller than the majority carrier concentration. At very high forward currents, high-level injection occurs, and the current grows linearly with voltage rather than exponentially.
• No Series Resistance: The formula models only the depletion region. It completely ignores the ohmic resistance of the P and N bulk regions and the metal contacts.
• DC Steady State: The equation is strictly for DC. It does not account for junction capacitance (C_J) or diffusion capacitance, which dominate during high-frequency switching or transient recovery.
• Reverse Breakdown: The formula accurately predicts the tiny leakage current (I ≈ -I_S) in reverse bias, but it completely fails to model the avalanche or Zener breakdown regions where current spikes vertically.

Unit Mistakes That Break the Math

Because the variable V_D sits inside an exponent, minor unit errors result in catastrophic calculation failures.

• The mV Trap: If you calculate V_T as 25.85 (thinking in millivolts) but input V_D as 0.7 (thinking in volts), your exponent becomes 0.7 / 25.85 = 0.027. The calculated current will be microamps instead of amps. Rule: Always convert V_T to Volts before dividing.
• The I_S Scale Trap: Datasheets rarely list I_S directly. If you extract a leakage current spec of 5 μA at 25V reverse bias and plug that directly into the equation as I_S, your forward current calculations will be off by a factor of 1,000. Real I_S is derived from curve-fitting the forward log-linear region, not from reverse leakage specs (which are dominated by surface leakage, not bulk diffusion).
• Celsius vs. Kelvin: Plugging 25 into the T variable instead of 298.15 will shrink V_T by a factor of 12, causing the exponential term to overflow any standard calculator.

Realistic Answer Magnitudes

When debugging your math, use these sanity checks for standard silicon PN junctions at room temperature:

• V_D should almost always fall between 0.4 V and 0.9 V for forward currents between 1 mA and 10 A. If you calculate 2.5 V, you forgot to subtract the 1 in the rearranged formula, or your I_S is too small.
• n should be between 1.0 and 2.0. If you back-calculate an ideality factor of 4.5, your measurement probes are picking up voltage drop across the test leads, not just the diode junction.

Frequently Asked Questions

Why does the diode equation fail at high forward currents?

At high currents (typically above 100 mA for small signal diodes), the voltage drop across the physical silicon bulk and the wire bonds becomes significant. This is modeled in SPICE as a series resistor (R_S). The true voltage across the junction is V_D(junction) = V_D(measured) - (I · R_S). Because the Shockley equation only knows about the junction, it underpredicts the total measured voltage drop at high currents. To fix this, engineers use the modified equation: V_total = n·V_T·ln(I/I_S) + I·R_S.

How do I measure the ideality factor (n) on the bench?

You cannot measure n with a single multimeter reading. You must take multiple voltage and current readings across the diode's forward bias region (where series resistance is negligible, usually between 10 μA and 1 mA). Plot the natural log of the current, ln(I), on the Y-axis against the measured voltage, V_D, on the X-axis. The slope of this linear region is equal to 1 / (n · V_T). Since you know V_T (approx 25.85 mV at room temp), you can solve for n. A slope of roughly 14.5 V^-1 indicates an ideality factor of roughly 1.8, typical for recombination-dominated current in a standard silicon rectifier.

Does the Shockley equation apply to Schottky and Zener diodes?

It applies to Schottky diodes in the forward bias region, but with vastly different parameters. Schottky junctions (metal-semiconductor) have a much larger I_S (often 10^-6 A to 10^-9 A) and a lower forward voltage drop (0.2 V to 0.4 V). The ideality factor n for Schottkys is usually very close to 1.0. However, the equation does not apply to Zener diodes operating in their reverse breakdown region. The Shockley model predicts a flat, tiny reverse current (-I_S), whereas a Zener diode exhibits a sharp, controlled avalanche or tunneling current knee that requires a completely different piecewise mathematical model.

What happens to the diode equation at extremely low temperatures?

As temperature drops toward absolute zero, the thermal voltage V_T shrinks toward zero, which mathematically forces the exponential term toward infinity for any fixed V_D. In physical reality, the equation breaks down due to 'carrier freeze-out.' At cryogenic temperatures (e.g., below 100 K), the dopant atoms lack the thermal energy to ionize, meaning there are almost no free charge carriers. I_S drops to effectively zero, and the silicon behaves more like an insulator than a semiconductor. Standard SPICE diode models will fail to converge in simulations at these temperatures without specialized cryogenic parameter sets.