The equation of diode current, universally known as the Shockley ideal diode equation, models the non-linear current-voltage (I-V) characteristic of a PN junction. Whether you are biasing a 1N4148 signal diode or calculating the forward drop of a Schottky rectifier, this formula bridges semiconductor physics and practical circuit design. Below, we break down the formula, define every variable with realistic bench magnitudes, and walk through solved problems and real-world failure scenarios.
The Core Formula and Symbol Definitions
The standard Shockley equation for diode current is expressed as:
ID = IS * (e(VD / (n * VT)) - 1)
Where the thermal voltage VT = (k * T) / q. According to the NIST physical constants, the Boltzmann constant (k) is 1.380649 × 10-23 J/K and the elementary charge (q) is 1.602176634 × 10-19 C.
| Symbol | Parameter | Unit | Realistic Bench Magnitude |
|---|---|---|---|
| ID | Diode forward current | Amperes (A) | 1 mA to 1 A (signal to power) |
| IS | Reverse saturation current | Amperes (A) | 1 pA to 100 nA (highly temp-dependent) |
| VD | Voltage across the diode | Volts (V) | 0.2 V (Schottky) to 0.7 V (Silicon PN) |
| n | Ideality factor (emission coefficient) | Dimensionless | 1.0 (ideal) to 2.0 (recombination-heavy) |
| VT | Thermal voltage | Volts (V) | ~25.85 mV at 300 K (26.85 °C) |
When the Formula Applies (and Its Assumptions)
The Shockley equation is derived under specific physical assumptions. It applies strictly to low-level injection in a uniformly doped PN junction. It assumes:
- No series bulk resistance: It ignores the voltage drop across the semiconductor bulk material and contacts (which dominates at high currents).
- No avalanche breakdown: It does not model reverse breakdown (Zener or avalanche regions).
- Uniform temperature: It assumes the junction temperature is uniform and static.
For a deeper look at the semiconductor physics governing these assumptions, HyperPhysics provides an excellent breakdown of the depletion region mechanics.
Rearranged Forms for Bench Debugging
On the bench, you rarely solve for ID directly. More often, you measure VD and ID to extract the diode's physical parameters, or you need to find the required bias voltage. Here are the algebraically rearranged forms solving for each variable:
- Solve for VD (Forward Voltage):
VD = n * VT * ln((ID / IS) + 1) - Solve for IS (Saturation Current):
IS = ID / (e(VD / (n * VT)) - 1) - Solve for n (Ideality Factor):
n = VD / (VT * ln((ID / IS) + 1)) - Solve for VT (Thermal Voltage / Junction Temp):
VT = VD / (n * ln((ID / IS) + 1))
Worked Problems with Strict Unit Tracking
Let's run two calculations. Notice how strictly tracking units—especially converting milliamps to amps and millivolts to volts—prevents catastrophic math errors.
Problem 1: Finding Forward Current in a 1N4148 Signal Diode
Given: A 1N4148 diode has a measured forward voltage VD = 0.65 V. The ambient temperature is 300 K. The datasheet-derived ideality factor at this operating point is n = 1.75, and IS = 2.5 nA.
Find: ID
Step 1: Calculate Thermal Voltage (VT)
VT = (1.3806 × 10-23 J/K * 300 K) / (1.6022 × 10-19 C) = 0.02585 V (or 25.85 mV)
Step 2: Convert all inputs to base SI units (Volts, Amperes)
VD = 0.65 V
IS = 2.5 × 10-9 A
Step 3: Calculate the exponent
Exponent = VD / (n * VT) = 0.65 / (1.75 * 0.02585) = 0.65 / 0.0452375 = 14.3686
Step 4: Solve for ID
ID = 2.5 × 10-9 * (e14.3686 - 1)
ID = 2.5 × 10-9 * (1,738,025 - 1)
ID = 2.5 × 10-9 * 1,738,024 = 0.004345 A (or 4.35 mA)
Problem 2: Finding Forward Voltage in a BAT54 Schottky Diode
Given: A BAT54 Schottky diode is passing ID = 15 mA. Room temperature is 27 °C. For Schottky junctions, n = 1.1 and IS = 50 nA.
Find: VD
Step 1: Convert Temperature to Kelvin
T = 27 °C + 273.15 = 300.15 K
Step 2: Calculate VT
VT = (1.3806 × 10-23 * 300.15) / 1.6022 × 10-19 = 0.02586 V
Step 3: Convert current to Amperes
ID = 0.015 A
IS = 50 × 10-9 A
Step 4: Apply the rearranged VD formula
VD = n * VT * ln((ID / IS) + 1)
Ratio = 0.015 / (50 × 10-9) = 300,000
VD = 1.1 * 0.02586 * ln(300,000 + 1)
VD = 0.028446 * 12.6115 = 0.358 V
Real-World Scenario: The Thermal Runaway Trap
Equations on paper behave perfectly; components on a breadboard do not. Here is a scenario where blind trust in the ideal equation of diode current leads to hardware failure.
The Setup
You are designing a 2A DC power supply rectifier. To save costs and board space, you decide to place two 1N4007 diodes (rated for 1A each) in parallel, assuming the current will split evenly at 1A per diode. You calculate the forward voltage drop using the Shockley equation at 1A, finding VD ≈ 0.8V. You assume a total power dissipation of 0.8W per diode, which seems manageable without a heatsink.
The Numbers
According to the ideal equation, if both diodes have identical IS and n, and are at the exact same temperature, they will share the 2A load perfectly. VD remains clamped at ~0.8V for both.
The Outcome
Within 30 seconds of powering on, one diode begins to smoke and eventually fails short-circuit, taking the second diode out with it. The power supply trips its overcurrent protection.
What Went Wrong
The ideal equation lacks two critical real-world parameters: bulk series resistance (RS) and thermal coupling.
- Manufacturing Tolerance: No two diodes have the exact same IS. The diode with the slightly higher IS draws slightly more current initially.
- Thermal Feedback: IS is highly temperature-dependent, roughly doubling for every 10 °C rise in junction temperature. The diode drawing more current heats up faster, which increases its IS, which causes it to draw even more current.
- Missing RS: At 1A, the voltage drop across the semiconductor bulk resistance (which the Shockley equation ignores) becomes significant. Without this natural ballasting resistance, the exponential nature of the equation causes current hogging.
The Fix: Never parallel standard PN diodes without ballast resistors (e.g., 0.1Ω in series with each anode) to force current sharing, or simply use a single diode rated for 3A, like the 1N5408.
Unit Mistakes That Break the Math
Warning: The Exponent is Unforgiving
Because the diode current relies on an exponential function, even minor unit errors result in answers that are off by orders of magnitude. Watch out for these three bench-top traps:
- Celsius vs. Kelvin: The thermal voltage formula VT = kT/q requires absolute temperature. Plugging 25 °C directly into T yields a VT of 2.15 × 10-24 V, collapsing your calculated current to virtually zero. Always add 273.15.
- Millivolts vs. Volts in the Exponent: VT is commonly memorized as "25.85". If you plug 25.85 into the denominator instead of 0.02585, your exponent shrinks by a factor of 1000, and e0.014 evaluates to ~1.01, yielding a calculated current in the picoamp range instead of milliamps.
- Dropping the "-1" in Reverse Bias: In forward bias, e(VD/nVT) is massive, so the "-1" is negligible. But if you apply a reverse voltage (e.g., VD = -5V), the exponent becomes a large negative number, making ex approach 0. If you forget the "-1", your formula calculates 0 A. With the "-1", it correctly calculates ID = -IS (the reverse leakage current).
Mastering the equation of diode current means knowing when to trust the math and when to trust your thermal camera. Use the Shockley equation for small-signal biasing and theoretical analysis, but always verify high-current designs against datasheet V-I curves that account for bulk resistance and thermal derating. For further reading on practical semiconductor modeling, the All About Circuits semiconductor textbook offers excellent supplementary circuit models.
