The standard equation of diode current is the Shockley diode equation, which models the steady-state direct current (DC) flowing through an ideal p-n junction. The direct answer for the forward-bias current is I = IS (eVD / (n × VT) - 1). While hobbyists often approximate a silicon diode as a simple 0.7V voltage drop, bench-level circuit design, SPICE simulation, and precision analog work require the exact exponential relationship. Below, we break down the formula, rearrange it for practical bench measurements, and solve real-world problems using actual component parameters.
The Shockley Equation: Symbols and Definitions
The Shockley equation defines the non-linear relationship between the voltage applied across a diode and the resulting current. It is derived from the physics of minority carrier diffusion across the depletion region. According to All About Circuits, this model assumes low-level injection and ignores bulk semiconductor resistance.
The core formula is:
I = IS (eVD / (n × VT) - 1)
| Symbol | Parameter | Standard Unit | Typical Silicon Value |
|---|---|---|---|
| I | Diode forward current | Amperes (A) | 1 mA to 1 A |
| IS | Reverse saturation current | Amperes (A) | 10-15 A to 10-9 A |
| VD | Voltage across the diode | Volts (V) | 0.5 V to 0.8 V |
| n | Ideality factor (emission coefficient) | Dimensionless | 1.0 to 2.0 |
| VT | Thermal voltage (kT/q) | Volts (V) | ~0.02585 V at 300K |
Note on VT: Thermal voltage is derived from fundamental constants: VT = (k × T) / q, where k is the Boltzmann constant (1.3806 × 10-23 J/K), T is the absolute temperature in Kelvin, and q is the elementary charge (1.602 × 10-19 C).
Rearranged Forms: Solving for Voltage, Saturation Current, and Ideality
On the bench, you rarely know all variables upfront. When characterizing an unknown diode or reverse-engineering a SPICE model, you need to isolate specific variables. Here are the algebraically rearranged forms of the equation of diode current:
- Solve for Diode Voltage (VD):
VD = n × VT × ln((I / IS) + 1)
Use case: Calculating the exact forward voltage drop when you know the operating current and the diode's SPICE parameters. - Solve for Saturation Current (IS):
IS = I / (eVD / (n × VT) - 1)
Use case: Extracting the leakage/saturation parameter from a datasheet's V-I curve at a specific test point. - Solve for Ideality Factor (n):
n = VD / (VT × ln((I / IS) + 1))
Use case: Determining the recombination characteristics of a junction by measuring VD at a known I.
Worked Examples with Unit Tracking
Abstract formulas are useless if you drop a decimal during calculation. Below are two solved problems using real-world component parameters, with strict unit tracking to prevent magnitude errors.
Problem 1: Finding Forward Current for a 1N4148 Signal Diode
Given: A standard 1N4148 switching diode has a forward voltage VD = 0.65 V. The manufacturer SPICE model specifies IS = 2.52 nA and n = 1.75. The ambient temperature is 27°C (T = 300 K). Find the forward current I.
- Calculate Thermal Voltage (VT):
VT = (1.3806 × 10-23 J/K × 300 K) / (1.602 × 10-19 C)
VT = 0.02585 Volts (or 25.85 mV). - Convert IS to Base Units (Amperes):
IS = 2.52 nA = 2.52 × 10-9 A. - Calculate the Exponent Term:
Exponent = VD / (n × VT)
Exponent = 0.65 V / (1.75 × 0.02585 V) = 0.65 / 0.0452375 = 14.3686 (dimensionless). - Evaluate the Exponential:
e14.3686 ≈ 1,738,562. - Calculate Final Current (I):
I = 2.52 × 10-9 A × (1,738,562 - 1)
I = 2.52 × 10-9 × 1,738,561 = 0.00438 Amperes.
Answer: I = 4.38 mA.
Problem 2: Finding Forward Voltage for a 1N4007 Rectifier
Given: A 1N4007 rectifier is passing I = 1.0 A. Assume an idealized IS = 1.0 pA (10-12 A) and n = 1.5 at T = 300 K. Find the theoretical junction voltage VD.
- Identify Knowns in Base Units:
I = 1.0 A, IS = 1.0 × 10-12 A, VT = 0.02585 V. - Calculate the Current Ratio:
I / IS = 1.0 / 10-12 = 1012 (one trillion). - Apply the Natural Logarithm:
ln(1012 + 1) ≈ ln(1012) = 27.631. - Calculate VD:
VD = n × VT × 27.631
VD = 1.5 × 0.02585 V × 27.631 = 1.071 Volts.
Answer: Theoretical VD = 1.07 V.
Assumptions, Limits, and Common Unit Mistakes
To use the equation of diode current effectively, you must understand its boundaries. As noted by Electronics Tutorials, the Shockley model is an approximation that breaks down under specific physical conditions.
When the Formula Applies (and When It Doesn't)
- Applies to: Steady-state DC conditions, low-to-moderate forward bias (where current is dominated by diffusion), and constant junction temperature.
- Fails at High Forward Bias: Above ~100 mA for small signal diodes, bulk series resistance dominates. The modified SPICE equation adds an
I × RSterm to VD. - Fails at High Reverse Bias: The equation predicts a constant reverse current of -IS. In reality, avalanche or Zener breakdown occurs, causing current to spike violently, which this formula does not model.
- Ignores AC/Transient Effects: It does not account for junction capacitance (Cj) or reverse recovery time (trr).
Realistic Answer Magnitudes
If your calculator spits out a number outside these ranges, you have a math error:
- IS (Saturation Current): Should be microscopic. For silicon, expect 1 fA (10-15 A) to 10 nA (10-8 A). If you calculate an IS in the milliamp range, your inputs are wrong.
- VD (Forward Voltage): For standard silicon, 0.5V to 0.85V. For Schottky, 0.15V to 0.45V. For LEDs, 1.8V (red) to 3.3V (blue/white).
- n (Ideality Factor): Strictly between 1.0 (ideal diffusion) and 2.0 (recombination-dominated). If you calculate n = 4.5, your measurement points are flawed or the diode is damaged.
Unit Mistakes That Break the Math
- Temperature in Celsius: Plugging 27°C directly into VT = kT/q yields a thermal voltage near zero, causing the exponent to approach infinity and crashing your calculation. Always convert to Kelvin (T = °C + 273.15).
- Mixing mA and A: If IS is in nanoamps (10-9) but you input it as 10-3 (milliamps) by mistake, the ratio I/IS shrinks by a factor of a million, resulting in a drastically underestimated VD.
- VT in mV vs VD in V: If you use VD = 0.7 (Volts) but use VT = 25.85 (millivolts) without converting VT to 0.02585, your exponent will be 1000x too small, and the diode will appear to conduct almost zero current.
Frequently Asked Questions
Why does the equation of diode current fail at high forward bias?
The base Shockley equation assumes the semiconductor material outside the depletion region has zero resistance. In physical components like the 1N4007 or 1N5408, the bulk silicon and the metal-semiconductor contacts introduce a parasitic series resistance (RS), typically ranging from 0.05 Ω to 1 Ω. At low currents (e.g., 1 mA), the voltage drop across RS is negligible. At high currents (e.g., 3 A), the I × RS drop can exceed the junction voltage itself. To model this accurately, engineers use the modified equation: Vtotal = VD + (I × RS).
How do you calculate thermal voltage in the diode equation?
Thermal voltage (VT) represents the voltage equivalent of the thermal energy of charge carriers. You calculate it using the formula VT = (k × T) / q. At standard room temperature (27°C or 300.15 K), the math works out to: (1.3806 × 10-23 J/K × 300.15 K) / 1.602 × 10-19 C = 0.02586 Volts. For quick bench estimations, most engineers round this to 26 mV. Note that VT increases linearly with temperature, which is why diodes become more conductive (and drop less voltage for a given current) as they heat up.
What is a realistic ideality factor (n) for silicon vs Schottky diodes?
The ideality factor (n) accounts for deviations from the ideal p-n junction theory, primarily due to carrier recombination in the depletion region. For standard silicon signal diodes (like the 1N4148), n typically ranges from 1.5 to 1.8. For silicon power rectifiers, it sits closer to 1.2 to 1.5. Schottky diodes, which rely on a metal-semiconductor junction rather than a p-n junction, have very little depletion-region recombination; therefore, their ideality factor is much closer to the theoretical ideal, usually ranging from 1.05 to 1.2. If you are extracting n from a datasheet graph and get a value outside the 1.0 to 2.0 range, your reading of the logarithmic current axis is likely incorrect.






