The diode Shockley equation is the foundational mathematical model for the non-linear current-voltage (I-V) relationship of a PN junction. If you are designing a rectifier, simulating a clipping circuit in SPICE, or just trying to understand why your silicon diode drops 0.7V at 10mA but 0.6V at 1mA, this formula is the governing physics. The direct answer to how a diode behaves in forward bias is captured by the equation: I = IS (eVD / (n · VT) - 1).
The Core Formula and Symbol Definitions
Below is the standard Shockley diode equation. Unlike ideal resistors governed by Ohm's Law, the PN junction current scales exponentially with applied voltage.
I = IS (eVD / (n · VT) - 1)
| Symbol | Parameter | Standard Unit | Typical Bench Value (Silicon) |
|---|---|---|---|
| I | Diode forward current | Amperes (A) | 1 mA to 1 A |
| IS | Reverse saturation (scale) current | Amperes (A) | 10-15 A to 10-9 A |
| VD | Voltage across the diode junction | Volts (V) | 0.5 V to 0.8 V |
| n | Ideality factor (emission coefficient) | Dimensionless | 1.0 to 2.0 |
| VT | Thermal voltage | Volts (V) | ~0.02585 V (at 300 K) |
The thermal voltage VT is derived from fundamental constants: VT = kT / q, where k is the Boltzmann constant (1.3806 × 10-23 J/K), T is the absolute junction temperature in Kelvin, and q is the elementary charge (1.602 × 10-19 C). For standard room temperature (300 K / 27°C), VT is approximately 25.85 mV. You can verify this physics baseline via Georgia State University's HyperPhysics reference.
Assumptions, Limits, and Realistic Magnitudes
The Shockley equation is an idealized model. It applies beautifully to low-level injection in the forward-bias region, but it makes several assumptions that break down on a real workbench.
When the Formula Applies
- Forward Bias Region: It accurately predicts current from the microamp range up to moderate milliamp ranges where the junction's internal resistance is negligible.
- Low-Level Injection: It assumes the injected minority carrier concentration is much smaller than the majority carrier concentration.
- Uniform Doping: It assumes an abrupt or uniformly doped PN junction.
Where the Model Fails
At high forward currents (e.g., pushing 3A through a 1N5408), the bulk semiconductor material and the physical leads introduce series resistance (RS). The actual voltage becomes Vtotal = VD + (I · RS). The pure Shockley equation will massively underpredict the voltage drop at high currents. Furthermore, in reverse bias, the equation simply predicts I ≈ -IS, completely ignoring surface leakage currents and avalanche/Zener breakdown mechanisms.
Realistic Magnitudes and Unit Traps
Warning: The Unit Mistakes That Break the Math
1. Temperature in Celsius: If you plug 25°C directly into VT = kT/q instead of converting to 298.15 K, your thermal voltage will be near zero, and your exponential term will overflow.
2. IS Scaling: Datasheets and SPICE models often list IS in picoamps (pA) or nanoamps (nA). If your IS is 2.5 nA, you must input 2.5 × 10-9 A. Forgetting the negative exponent will yield currents in the giga-ampere range.
3. Thermal Voltage in mV: VT is 25.85 mV. If you use 25.85 in the denominator instead of 0.02585 V, the exponent shrinks by a factor of 1000, and the diode appears to never turn on.
Rearranged Forms for Circuit Analysis
On the bench, you rarely know all variables. Often, you measure the current and need to find the junction voltage, or you are characterizing an unknown diode and need to extract the ideality factor. Here are the algebraically rearranged forms:
- Solving for Junction Voltage (VD):
VD = n · VT · ln( (I / IS) + 1 ) - Solving for Saturation Current (IS):
IS = I / ( eVD / (n · VT) - 1 ) - Solving for Ideality Factor (n):
n = VD / ( VT · ln( (I / IS) + 1 ) )
For practical forward-bias calculations where I is significantly larger than IS (which is almost always true when the diode is conducting), the '+ 1' inside the natural log can be safely ignored, simplifying the VD equation to: VD ≈ n · VT · ln(I / IS).
Worked Examples with Unit Tracking
Let's run two real-world scenarios, tracking every unit to ensure the dimensional analysis holds up. For both examples, we assume a standard ambient junction temperature of 300 K, yielding a thermal voltage VT = 0.02585 V.
Problem 1: Calculating Forward Current (I) for a Signal Diode
Scenario: You are biasing a 1N4148 small-signal switching diode with a forward voltage of 0.65 V. The diode has a typical saturation current IS = 1.0 × 10-12 A (1 pA) and an ideality factor n = 1.2. What is the forward current?
- Identify knowns: VD = 0.65 V, IS = 1.0 × 10-12 A, n = 1.2, VT = 0.02585 V.
- Calculate the denominator of the exponent:
n · VT = 1.2 · 0.02585 V = 0.03102 V. - Calculate the full exponent:
VD / (n · VT) = 0.65 V / 0.03102 V = 20.9542 (dimensionless). - Evaluate the exponential term:
e20.9542 = 1,259,020,000 (dimensionless). - Subtract 1 and multiply by IS:
I = (1.0 × 10-12 A) · (1,259,020,000 - 1)
I ≈ (1.0 × 10-12 A) · (1.259 × 109)
I = 1.259 × 10-3 A. - Final Answer: I = 1.26 mA. (This aligns perfectly with typical 1N4148 bench measurements at 0.65V).
Problem 2: Calculating Forward Voltage Drop (VD) for a Power Rectifier
Scenario: A 1N4007 rectifier diode is conducting 1.0 A of continuous DC current. Due to recombination effects in the wider depletion region of a power diode, the ideality factor is higher: n = 1.8. The saturation current is IS = 1.5 × 10-9 A (1.5 nA). What is the junction voltage drop?
- Identify knowns: I = 1.0 A, IS = 1.5 × 10-9 A, n = 1.8, VT = 0.02585 V.
- Use the rearranged VD formula:
VD = n · VT · ln( (I / IS) + 1 ) - Calculate the current ratio:
I / IS = 1.0 A / (1.5 × 10-9 A) = 666,666,666.67 (dimensionless). - Add 1 and take the natural logarithm:
ln(666,666,667.67) = 20.3176 (dimensionless). - Calculate the thermal multiplier:
n · VT = 1.8 · 0.02585 V = 0.04653 V. - Multiply to find VD:
VD = 0.04653 V · 20.3176 = 0.9453 V. - Final Answer: VD = 0.945 V. (Note: A multimeter reading on a physical 1N4007 at 1A will likely show ~1.05V to 1.1V. The ~0.1V difference is the I · RS bulk resistance drop that the pure Shockley equation ignores).
Frequently Asked Questions
Why does the diode Shockley equation fail at high forward currents?
The pure Shockley equation only models the PN junction itself. It assumes the P and N semiconductor regions have zero electrical resistance. In reality, the physical silicon die, the bond wires, and the external leads all possess bulk series resistance (RS). At low currents (e.g., 1 mA), the voltage drop across this resistance is negligible (microvolts). At high currents (e.g., 3 A), the I · RS drop becomes significant, causing the I-V curve to transition from exponential to linear. SPICE simulation models fix this by adding an RS parameter in series with the ideal Shockley diode model.
How do I calculate thermal voltage (VT) if my ambient temperature changes?
Thermal voltage is strictly dependent on absolute temperature. You must convert your Celsius temperature to Kelvin by adding 273.15. The formula is VT = (k · T) / q. A useful bench shortcut is that VT scales linearly at approximately 0.08617 mV per degree Kelvin. Therefore, at 20°C (293.15 K), VT is 25.26 mV. At a hot 85°C (358.15 K), VT rises to 30.86 mV. This temperature shift is exactly why diode forward voltage drops decrease by roughly 2 mV per °C rise in temperature, a critical factor in thermal runaway scenarios for paralleled diodes.
What is the difference between the ideality factor (n) for silicon vs. Schottky diodes?
The ideality factor (n) accounts for deviations from the ideal diffusion current model. For standard silicon PN junction diodes (like the 1N4148 or 1N4007), n typically ranges from 1.0 to 2.0, often sitting around 1.2 to 1.5 due to recombination currents in the depletion region. Schottky diodes (like the BAT54 or 1N5819), however, rely on majority-carrier thermionic emission over a metal-semiconductor barrier rather than minority-carrier diffusion. Because they lack the minority-carrier recombination lag, Schottky diodes exhibit an ideality factor much closer to the theoretical ideal, typically between 1.05 and 1.15. For deeper semiconductor physics models, refer to the Electronics Tutorials PN Junction documentation.
