The Diode Current Formula (Shockley Equation) Defined

When you are designing a clipping circuit, sizing a rectifier, or debugging a temperature-sensitive analog sensor, the rule-of-thumb "a silicon diode drops 0.7V" is not enough. You need the exact exponential relationship between voltage and current. The diode current formula, formally known as the Shockley diode equation, models the ideal behavior of a p-n junction under forward and reverse bias.

The core formula is:

ID = IS (eVD / (n VT) - 1)

Here is the exact definition of every symbol, including the physical constants and typical bench values you will encounter in real datasheets and SPICE simulations.

Shockley Equation Symbol Reference
Symbol Parameter Unit Typical Value / Notes
ID Diode forward current Amperes (A) mA for signal diodes; A for power rectifiers
IS Reverse saturation current Amperes (A) 10-15 to 10-9 A (highly temp-dependent)
VD Voltage across the diode junction Volts (V) 0.5V to 0.8V for Si; 0.2V to 0.4V for Schottky
n Ideality factor (emission coefficient) Unitless 1.0 (ideal) to 2.0 (real silicon, e.g., 1.75 for 1N4148)
VT Thermal voltage (kT/q) Volts (V) ~0.02585 V at 300K (26.85°C)

Thermal voltage (VT) is derived from fundamental physics: VT = kT/q, where k is the Boltzmann constant (1.3806 × 10-23 J/K), T is absolute temperature in Kelvin, and q is the elementary charge (1.602 × 10-19 C). At a standard bench temperature of 300K, VT is exactly 25.85 mV.

Rearranged Forms: Solving for Voltage and Saturation Current

On the bench, you rarely know VD and need to find ID. Usually, you know the current your circuit is drawing (e.g., 15 mA through a biasing network) and need to calculate the exact voltage drop to set your transistor base bias. You must rearrange the formula.

Solving for Junction Voltage (VD)

By isolating the exponential term and taking the natural logarithm (ln) of both sides, we get:

VD = n VT ln(ID / IS + 1)

Practical approximation: In forward bias, ID is orders of magnitude larger than IS. Therefore, the "+ 1" becomes mathematically irrelevant. The working formula for forward bias simplifies to:

VD ≈ n VT ln(ID / IS)

Solving for Saturation Current (IS)

If you are characterizing an unknown diode on a curve tracer and need to extract its SPICE model parameters, you solve for IS:

IS = ID / (eVD / (n VT) - 1)

Assumptions, Limits, and Unit Traps

The Shockley equation is an elegant model, but it is not a universal law. It relies on specific physical assumptions that break down at the extremes of real-world operation.

  • Low-Level Injection: The formula assumes the injected minority carrier density is much lower than the majority carrier density. At high currents (typically >100 mA for small signal diodes), high-level injection occurs, and the current grows linearly rather than exponentially.
  • Ignores Bulk Series Resistance (RS): Real diodes have physical semiconductor bulk and lead resistance. At high currents, the voltage drop across RS dominates. A complete SPICE model adds ID * RS to the VD term.
  • No Breakdown Modeling: The formula predicts a tiny, constant reverse current (-IS) as VD goes negative. It completely fails to model Zener or Avalanche breakdown.
⚠️ Callout: The Unit Traps That Break Your Math

The most common reason hobbyists and students get "Math Error" or impossibly huge numbers on their calculators is unit mismatch in the exponent.

  • Trap 1: Thermal Voltage. VT is ~26 mV. If you type 26 into your calculator instead of 0.02585, your exponent becomes tiny, and your calculated current drops to near zero.
  • Trap 2: Saturation Current. Datasheets often list IS or reverse leakage in pA or nA. You must convert to base Amperes. A 1N4148 SPICE model uses IS = 2.52 nA. You must enter 2.52e-9, not 2.52.
  • Trap 3: Natural Log vs Base-10 Log. Ensure your calculator is using ln (base e), not log (base 10), when using the rearranged VD formula.

Worked Example 1: Forward Bias Signal Diode Calculation

Scenario: You are designing a precision analog clipper using a generic small-signal silicon diode. You apply exactly 0.65 V across the diode junction. What is the resulting forward current?

Given Parameters:

  • VD = 0.65 V
  • IS = 1.5 × 10-12 A (1.5 pA, typical for a generic small-signal Si diode)
  • n = 1.3 (ideality factor)
  • VT = 0.02585 V (assuming 300K bench temperature)

Step-by-Step Solution:

  1. Calculate the thermal voltage multiplier:
    n × VT = 1.3 × 0.02585 V = 0.033605 V
  2. Calculate the exponent:
    VD / (n VT) = 0.65 / 0.033605 = 19.3423
  3. Evaluate the exponential term:
    e19.3423 = 251,358,500
  4. Subtract 1:
    251,358,500 - 1 = 251,358,499 (The -1 is negligible here, proving the forward-bias assumption).
  5. Multiply by IS:
    ID = (1.5 × 10-12 A) × 251,358,499 = 0.000377 A

Final Answer: The diode current is 0.377 mA (or 377 µA). This magnitude is highly realistic for a signal diode operating just below the traditional "0.7V knee".

Worked Example 2: Finding Junction Voltage from Known Current

Scenario: You are biasing a transistor and routing 15 mA through a standard 1N4001 rectifier used as a crude voltage drop. You need to know the exact VD to calculate the base voltage.

Given Parameters:

  • ID = 15 mA = 0.015 A
  • IS = 5.0 × 10-12 A (5 pA)
  • n = 1.2
  • VT = 0.02585 V

Step-by-Step Solution (Using Rearranged Formula):

  1. Calculate the current ratio:
    ID / IS = 0.015 A / (5.0 × 10-12 A) = 3,000,000,000 (or 3 × 109)
  2. Take the natural logarithm:
    ln(3 × 109) = 21.821
  3. Calculate the thermal multiplier:
    n × VT = 1.2 × 0.02585 V = 0.03102 V
  4. Multiply to find VD:
    VD = 0.03102 V × 21.821 = 0.6769 V

Final Answer: The junction voltage drop is 0.677 V. This confirms the bench rule-of-thumb that a silicon diode at moderate current drops roughly 0.65V to 0.7V, but gives you the exact 3-digit precision needed for analog biasing.

Decision Path: Selecting the Right Diode for Your Calculated Current

Once you use the Shockley equation to determine the expected current (ID) and your circuit dictates the reverse voltage (VR), you must select a physical component. Do not leave this to an open-ended "it depends". Use this decision matrix to terminate your design with a concrete part number.

Diode Selection Decision Tree
Calculated ID (Forward) Max VR (Reverse) Primary Constraint Concrete Part Pick
< 100 mA < 75 V Fast switching, low capacitance 1N4148 (or 1N914)
100 mA to 1.0 A Up to 1000 V General purpose, high surge tolerance 1N4007 (Use 1N4001 for <50V)
1.0 A to 3.0 A < 40 V Low forward voltage drop (efficiency) 1N5822 (Schottky)
> 3.0 A Up to 600 V High power rectification, needs heatsink MBR360 or 6A10
💡 The Default Bench Recommendation

If your calculations put you in the ambiguous middle ground (e.g., 50 mA to 200 mA) and switching speed is not critical, default to the 1N4007. It is virtually indestructible on the bench, handles up to 1000V reverse bias, and costs pennies. If you are working strictly with low-voltage logic or high-frequency RF signals (< 50 mA), default to the 1N4148 for its 4ns reverse recovery time.

For deeper reading on semiconductor physics and SPICE modeling parameters, refer to the HyperPhysics diode theory pages hosted by Georgia State University, or the comprehensive PN Junction Diode tutorials at Electronics Tutorials. Mastering the math behind the component ensures your designs work on the first PCB spin, rather than requiring iterative bench hacks.