The Shockley Equation Diode: Formula and Worked Examples

The behavior of a semiconductor p-n junction is not linear; it is an exponential relationship governed by quantum mechanics and thermodynamics. For hobbyists and engineers alike, the Shockley diode equation is the foundational mathematical model that predicts how current flows through a diode based on the applied voltage and junction temperature. While circuit simulators like SPICE use heavily modified versions of this math to account for parasitic effects, the core equation remains the bedrock of semiconductor theory. Below, we break down the formula, define every variable, and walk through bench-realistic worked examples.

The Core Equation Diode Model and Symbol Definitions

The standard Shockley equation diode model defines the forward current ($I_D$) as a function of the voltage dropped across the junction ($V_D$). The formula is expressed as:

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

To use this formula accurately, you must understand the physical meaning and standard units of every symbol. The thermal voltage ($V_T$) is itself a derived value: $V_T = (k · T) / q$.

Rearranged Forms: Solving for Voltage, Ideality, and Temperature

On the bench, you rarely know $I_D$ and need to find it. More often, you are measuring a known current and need to back-calculate the junction voltage, or you are characterizing an unknown diode to find its ideality factor. Here are the algebraically rearranged forms of the equation diode model:

• Solving for Junction Voltage ($V_D$):
  VD = n · VT · ln( (ID / IS) + 1 )
• Solving for Ideality Factor ($n$):
  n = VD / ( VT · ln( (ID / IS) + 1 ) )
• Solving for Saturation Current ($I_S$):
  IS = ID / ( eVD / (n · VT) - 1 )

Assumptions, Limitations, and Unit Traps

The Shockley equation is an idealized model. According to HyperPhysics semiconductor references, it applies strictly under specific conditions and fails if you ignore its physical boundaries.

When it applies and its assumptions:

• Low-Level Injection: The injected minority carrier concentration must be much less than the majority carrier concentration.
• No Bulk Resistance: The equation models only the p-n junction. It ignores the series resistance ($R_S$) of the semiconductor bulk material and wire leads.
• Steady State: It assumes DC conditions, ignoring junction capacitance and reverse recovery time.

Realistic Answer Magnitudes:

If your calculation yields a forward voltage ($V_D$) of 2.5V for a standard silicon diode at 1A, your math is wrong. Realistic magnitudes for silicon are 0.5V to 1.1V. For $I_S$, expect extremely small numbers ($10^{-12}$ A to $10^{-8}$ A). If you calculate an $I_S$ in the milliamp range, you have a unit error.

Unit Mistakes That Break the Math:

1. Using Celsius instead of Kelvin: Plugging 25 into $T$ instead of 298.15 will collapse your thermal voltage and result in wildly exponential, incorrect currents.
2. Mixing mA and A in $I_S$: Datasheets and SPICE models sometimes list $I_S$ in nanoamps or picoamps. You must convert to base Amperes before plugging it into the exponent, or your $I_D$ result will be off by orders of magnitude.
3. Forgetting the $-1$ in Reverse Bias: If $V_D$ is negative, the exponent approaches 0, and $e^0 = 1$. The $1 - 1$ yields $I_D = -I_S$. If you drop the $-1$ for 'simplicity', you will calculate 0A reverse current instead of the actual leakage current.

Worked Example 1: Finding Forward Current in a 1N4148

Scenario: You have a 1N4148 small-signal switching diode on your breadboard. Your multimeter measures exactly 0.650 V across the diode terminals. What is the forward current?

Given Values (derived from typical 1N4148 SPICE parameters at 300 K):

• $V_D = 0.650$ V
• $I_S = 2.52 \times 10^{-9}$ A (2.52 nA)
• $n = 1.752$
• $T = 300$ K $\rightarrow V_T = 0.02585$ V

Step-by-Step Solution:

1. Calculate the denominator of the exponent:
   $n \cdot V_T = 1.752 \cdot 0.02585\text{ V} = 0.045289\text{ V}$
2. Calculate the full exponent:
   $V_D / (n \cdot V_T) = 0.650\text{ V} / 0.045289\text{ V} = 14.352$
3. Evaluate the exponential term:
   $e^{14.352} \approx 1,710,206$
4. Subtract 1 (negligible here, but shown for completeness):
   $1,710,206 - 1 = 1,710,205$
5. Multiply by $I_S$ to find $I_D$:
   $I_D = (2.52 \times 10^{-9}\text{ A}) \cdot 1,710,205 = 0.004309\text{ A}$

Final Answer: $I_D = 4.31$ mA. This is a highly realistic operating point for a 1N4148 in a standard 5V logic pull-up circuit.

Worked Example 2: Calculating Junction Voltage for a 1N4007 Rectifier

Scenario: You are designing a 1A power supply and using a 1N4007 rectifier. You need to know the exact junction voltage drop to calculate power dissipation ($P = V_D \cdot I_D$). We will use the rearranged formula to solve for $V_D$.

Given Values (typical 1N4007 parameters at 300 K):

• $I_D = 1.0$ A
• $I_S = 2.7 \times 10^{-8}$ A (27 nA — note this is higher than the 1N4148 due to the massive physical junction area required to handle 1A)
• $n = 2.0$ (power rectifiers are heavily recombination-dominated)
• $V_T = 0.02585$ V

Step-by-Step Solution:

1. Set up the rearranged equation:
   $V_D = n \cdot V_T \cdot \ln( (I_D / I_S) + 1 )$
2. Calculate the current ratio:
   $I_D / I_S = 1.0\text{ A} / (2.7 \times 10^{-8}\text{ A}) = 37,037,037$
3. Add 1 and take the natural log:
   $\ln(37,037,038) \approx 17.427$
4. Multiply by the thermal voltage term:
   $n \cdot V_T = 2.0 \cdot 0.02585\text{ V} = 0.0517\text{ V}$
5. Calculate final $V_D$:
   $V_D = 0.0517\text{ V} \cdot 17.427 = 0.9009\text{ V}$

Final Answer: $V_D \approx 0.90$ V. This aligns perfectly with the 1N4007 datasheet, which specifies a typical forward voltage of 0.9V to 1.0V at 1.0A. Power dissipation is therefore $0.90\text{ V} \cdot 1.0\text{ A} = 0.90\text{ W}$.

Decision Tree: Picking the Right Diode for Your Operating Point

The equation diode model proves that $I_S$ and $n$ vary wildly depending on the physical construction of the semiconductor. You cannot use a small-signal diode for a power supply, nor can you use a power rectifier for high-frequency switching. Use this decision path to terminate your design process with a concrete part number.

By matching your target $I_D$ and acceptable $V_D$ to the physical parameters defined in the Shockley equation, you eliminate the guesswork from component selection and ensure your circuit operates within safe thermal limits.