Current Density Equations are the transport laws expressing total carrier current flow as the sum of drift (field-driven) and diffusion (concentration-gradient-driven) components — they connect the electrostatic potential and carrier density distributions solved by the Poisson and continuity equations to the actual current flowing through every point in a semiconductor device.

What Are the Current Density Equations?

• Electron Current: J_n = qnmu_nE + qD_n*(dn/dx), where the first term is drift (carriers moving in the electric field direction) and the second term is diffusion (carriers moving down the concentration gradient).
• Hole Current: J_p = qpmu_pE - qD_p*(dp/dx), with drift in the field direction and diffusion down the hole concentration gradient (note the sign difference from electrons).
• Einstein Connection: Diffusivity D and mobility mu are not independent — they are related by D = mu*kT/q, halving the number of transport parameters required and ensuring thermodynamic consistency.
• Total Current: The total electrical current density is J = J_n + J_p — both carrier types contribute to the current at every point, with their relative contributions determined by the local electric field and carrier gradients.

Why the Current Density Equations Matter

• Drift vs. Diffusion Regimes: Different device regions are dominated by different current mechanisms — the MOSFET channel above threshold is drift-dominated (field-driven at high field); the base of a bipolar transistor is diffusion-dominated; the subthreshold MOSFET channel is also diffusion-dominated. Understanding which mechanism controls current is essential for device optimization.
• I-V Characteristics: Integrating the current density equations over the device cross-section gives terminal current as a function of applied voltage — the measured I-V characteristic that defines transistor performance. Compact model equations such as BSIM are closed-form approximations to the exact current density integrals.
• Equilibrium Condition: At thermal equilibrium, J_n = J_p = 0 everywhere — drift and diffusion exactly cancel. This requires that the electric field created by band bending precisely compensates the concentration gradient at every point, a condition maintained by the Fermi level being spatially constant.
• Quasi-Fermi Level Representation: An equivalent and often more physically transparent form is J_n = qnmu_n*(dE_Fn/dx) / q, where E_Fn is the electron quasi-Fermi level — current flows whenever quasi-Fermi levels have a spatial gradient, providing an elegant graphical interpretation using band diagrams.
• High-Field Extensions: At high electric fields (above approximately 10^4 V/cm in silicon), carriers reach velocity saturation and the linear drift term mu*E must be replaced by a velocity-saturation model that caps the drift current — required for accurate short-channel transistor simulation.

How the Current Density Equations Are Used in Practice

• TCAD Implementation: The current density equations are discretized on the device mesh using the Scharfetter-Gummel scheme, which handles the exponential variation of carrier density with potential to provide stable, convergent solutions across many orders of magnitude in carrier concentration.
• Compact Model Foundation: Long-channel MOSFET current formulas (linear and saturation I-V), diode equations, and bipolar transistor gain expressions are all derived from closed-form integration of the current density equations under appropriate approximations.
• Current Flow Visualization: TCAD post-processing visualizes current flow line plots (streamlines of J_n and J_p) throughout the device, enabling identification of parasitic current paths, leakage channels, and efficiency-limiting recombination zones.

Current Density Equations are the transport laws at the heart of semiconductor device physics — expressing how both drift in electric fields and diffusion down concentration gradients contribute to current flow, they connect the electrostatics and carrier statistics solved by Poisson and continuity equations to the observable terminal currents that define device performance and are parameterized in every compact model used in circuit simulation.