Shockley Diode Interactive Calculator

Results:

The saturation current I_S represents the reverse leakage current that flows when the diode is reverse-biased by more than a few thermal voltages (typically greater than 100 mV). This current arises from thermally generated minority carriers diffusing across the junction—electrons in the p-region and holes in the n-region that are swept across by the built-in electric field. The magnitude is proportional to the intrinsic carrier concentration squared (n_i²), the junction area, and the diffusion coefficient, while inversely proportional to the doping concentration and minority carrier diffusion length. For silicon at room temperature, n_i ≈ 1×10¹⁰ cm^-3, which is tiny compared to typical doping levels of 10¹⁶-10¹⁸ cm^-3. This explains why I_S is typically in the range of 10^-15 to 10^-9 amperes for practical devices. Higher I_S values indicate either high temperature operation, low doping concentrations, poor crystal quality with high defect density, or deliberate design choices such as in Schottky diodes where I_S is several orders of magnitude higher than in p-n junction diodes.

The ideality factor quantifies deviations from the ideal Shockley model, which assumes current flow is dominated purely by minority carrier diffusion. When n = 1.0, diffusion current is the sole mechanism. However, real diodes experience recombination-generation current in the depletion region, where electron-hole pairs recombine through trap states at mid-gap energies. This recombination current has a different voltage dependence—proportional to exp(qV_D/2kT) instead of exp(qV_D/kT)—which corresponds to n = 2.0. The observed ideality factor represents a weighted combination of these mechanisms. High-quality silicon diodes at moderate to high current densities show n ≈ 1.02-1.08 because diffusion dominates. At very low currents (nanoampere range) or in diodes with high defect densities, recombination becomes significant and n approaches 1.5-2.0. Schottky diodes achieve n very close to 1.0 because they operate via majority carrier transport over a potential barrier with no minority carrier injection. Surface recombination, tunneling through thin barriers, and non-uniform current distribution across large area junctions can push n even higher than 2.0 in poorly designed devices.

The Shockley equation fails under several conditions that violate its underlying assumptions. At high forward currents, series resistance from the neutral regions, contacts, and bondwires becomes significant, adding an ohmic voltage drop IR_s that reduces the actual junction voltage below the applied terminal voltage. For currents above approximately 10 mA in small-signal diodes or 1 A in power rectifiers, this resistance effect dominates and the I-V curve becomes nearly linear rather than exponential. High-level injection occurs when the injected minority carrier density becomes comparable to the majority carrier density, changing the boundary conditions—this typically happens at current densities exceeding 100 A/cm² and requires solving the full drift-diffusion equations numerically. At very high reverse voltages approaching breakdown, tunneling mechanisms (Zener effect below 6 V, avalanche multiplication above 6 V) generate currents that increase dramatically, rendering the simple exp(V/nV_T) term meaningless. The equation also assumes steady-state DC conditions and fails for high-frequency AC signals where junction capacitance and carrier transit time effects create frequency-dependent impedance. Finally, at cryogenic temperatures below approximately 150 K, carrier freeze-out reduces the effective doping and the assumption of complete ionization breaks down.

The standard method involves measuring forward I-V characteristics at multiple current levels and extracting parameters from the linearized form of the Shockley equation. Taking the natural logarithm yields: ln(I + I_S) ≈ ln(I_S) + V_D/(nV_T). For I >> I_S, this simplifies to ln(I) ≈ ln(I_S) + V_D/(nV_T), which is a straight line when ln(I) is plotted versus V_D. The slope is 1/(nV_T), allowing extraction of n if temperature is known, and the y-intercept at V_D = 0 gives ln(I_S). In practice, measure the diode at 5-10 current points spanning 10 µA to 10 mA, use a precision source-measure unit with four-wire Kelvin connections to eliminate lead resistance, and maintain constant temperature in a controlled chamber. The slope method works best in the range where diffusion current dominates (typically 100 µA to 1 mA for small-signal diodes). At lower currents, recombination effects increase the apparent ideality factor, while at higher currents, series resistance reduces the slope. Advanced techniques involve measuring at multiple temperatures to separately extract the activation energy and prefactor of I_S, or using reverse-bias measurements where the small reverse current directly approaches I_S for large negative voltages.

The "-1" term ensures that I = 0 when V_D = 0 (zero bias condition) and that the equation correctly predicts negative current (reverse bias) when V_D < 0. It arises from the boundary condition that the excess minority carrier concentration at the depletion region edge must go to zero at equilibrium. Mathematically, the full solution to the diffusion equation yields exp(qV_D/kT) - 1, not just exp(qV_D/kT). For forward bias voltages exceeding approximately 100 mV (about 4V_T at room temperature), the exponential term is much larger than unity—exp(0.1/0.026) ≈ 48—so the "-1" contributes less than 2% error and can be dropped for simplified calculations. However, this approximation completely fails in three critical situations: reverse bias operation where exp(V_D/nV_T) approaches zero and the "-1" term dominates the current; zero-bias resistance calculations where dI/dV = (I_S/nV_T) at V_D = 0, a result that requires the "-1" term; and noise analysis where shot noise is proportional to 2qI_S at zero bias. Precision circuit analysis for applications like ultra-sensitive photodetectors or Johnson noise thermometry must retain the full equation.

Ionizing radiation creates permanent damage in semiconductor crystal lattices through displacement of silicon atoms and introduction of deep-level trap states. These defects dramatically increase recombination rates, reducing minority carrier lifetime from typical values of 1-10 microseconds down to nanoseconds after exposure to 100-1000 krad (Si) total ionizing dose. The immediate effect is an increase in reverse leakage current—I_S can increase by factors of 10³ to 10⁶ after heavy irradiation, which is catastrophic for applications requiring low leakage such as switched integrators or long-duration charge storage. The forward voltage at a given current also increases because enhanced recombination increases the effective ideality factor toward n = 2.0. Neutron radiation is particularly damaging because the recoil energy from neutron collisions creates large displacement cascades, while gamma and electron radiation primarily cause ionization damage that can partially anneal at elevated temperatures. Radiation-hardened diodes use specialized processing including neutron transmutation doping for precise defect engineering, gettering layers to trap mobile impurities, and guard ring structures to control surface leakage paths. For space applications, designers must derate diode current capability by 50-80% to account for end-of-life degradation after years of accumulated radiation dose.