Maxwell–Boltzmann Statistics: Non-Degenerate Carrier Transport, Thermal Velocity Distributions, and Semiconductor Kinetic Theory

Executive Overview

Maxwell–Boltzmann (MB) statistics describes the thermodynamic energy distribution and velocity kinetics of non-interacting, classical particles in thermal equilibrium. In solid-state physics and semiconductor device engineering, Maxwell–Boltzmann statistics serves as the non-degenerate approximation to quantum Fermi–Dirac statistics, valid when carrier concentrations are well below the quantum effective density of states ($n \ll N_c$, $p \ll N_v$) and the Fermi level $E_F$ lies deep within the bandgap ($E_c - E_F \ge 3 k_B T$). Under these conditions, quantum state filling effects and Pauli exclusion can be neglected, allowing electron and hole dynamics to be modeled via classical kinetic theory. Maxwell–Boltzmann kinetics underpins core semiconductor relationships, including the Law of Mass Action ($n \cdot p = n_i^2$), the Einstein Relation between carrier diffusion and mobility ($D/\mu = k_B T / q$), standard drift-diffusion transport equations, thermionic emission over Schottky barriers, and classical hot-carrier velocity distributions. This article provides a comprehensive theoretical derivation, velocity distribution kinetics, validity boundary conditions, Python numerical scripts, and semiconductor engineering applications.

Theoretical Derivation & Non-Degenerate Limit

Classical Kinetic Postulates

Maxwell–Boltzmann statistics applies to systems where: 1. Particles are distinguishable or sufficiently sparse that quantum mechanical wave packet overlap is negligible ($\lambda_{\text{thermal}} \ll d_{\text{interparticle}}$). 2. Any number of particles can occupy a single state ($n_i \ge 0$). 3. The average occupation number of any state is much less than unity ($f(E) \ll 1$).

Classical Limit of Fermi–Dirac Statistics

The exact quantum Fermi–Dirac distribution for electrons in the conduction band is:

When the Fermi level $E_F$ is situated deep inside the bandgap such that $(E_c - E_F) \ge 3 k_B T$ (approx. $\ge 0.077\text{ eV}$ at $300\text{ K}$), the exponential term for any conduction band state $E \ge E_c$ satisfies:

Neglecting the $+1$ in the denominator yields the Maxwell–Boltzmann distribution function:

Where $A = \exp(E_F / k_B T)$ is the fugacity coefficient.

Maxwellian Velocity & Energy Distributions

In a non-degenerate 3D semiconductor crystal with isotropic parabolic bands ($E(v) = \frac{1}{2} m_n^* v^2$), the probability density function $P(v) dv$ for an electron possessing speed between $v$ and $v + dv$ is obtained by integrating over velocity space angles in spherical coordinates:

Characteristic Velocity Quantities

From the velocity distribution $P(v)$, three fundamental velocity metrics are derived:

1. Most Probable Velocity ($v_p$): The peak of $P(v)$, found by setting $\frac{dP(v)}{dv} = 0$:

2. Mean / Average Thermal Velocity ($\langle v \rangle$ or $v_{\text{th}}$):

For electrons in silicon at $300\text{ K}$ ($m_n^* \approx 0.26 m_0$): $$v_{\text{th}} \approx \sqrt{\frac{8 \times (1.38 \times 10^{-23}) \times 300}{\pi \times (0.26 \times 9.11 \times 10^{-31})}} \approx 2.3 \times 10^7\text{ cm/s} \quad (2.3 \times 10^5\text{ m/s})$$

3. Root-Mean-Square (RMS) Velocity ($v_{\text{rms}}$):

The average kinetic energy of a classical non-degenerate carrier gas is directly proportional to temperature:

   P(v) Probability Density
    ^
    |      *  (v_p = sqrt(2 k_B T / m*))
    |    *   *
    |   *     *  <v> = sqrt(8 k_B T / pi m*)
    |  *        *   v_rms = sqrt(3 k_B T / m*)
    | *          *
    +----------------------------------------> Speed v

Carrier Density & Mass Action Law

Conduction Band Electron Concentration ($n$)

Integrating $f_{\text{MB}}(E)$ multiplied by the 3D density of states $N_c(E) = \frac{1}{2\pi^2} \left( \frac{2m_n^*}{\hbar^2} \right)^{3/2} \sqrt{E - E_c}$:

Where $N_c$ is the effective density of states in the conduction band:

Valence Band Hole Concentration ($p$)

Similarly, hole concentration $p$ in a non-degenerate valence band is:

Where $N_v = 2 \left( \frac{2\pi m_p^* k_B T}{h^2} \right)^{3/2}$.

The Law of Mass Action

Multiplying $n$ and $p$ eliminates the Fermi level $E_F$, yielding the fundamental thermodynamic relation for non-degenerate semiconductors in thermal equilibrium:

Where $n_i(T)$ is the intrinsic carrier concentration:

For silicon at $300\text{ K}$ ($E_g = 1.12\text{ eV}$): $n_i \approx 1.0 \times 10^{10}\text{ cm}^{-3}$. The Law of Mass Action holds independently of doping concentrations ($N_D, N_A$), provided the semiconductor remains non-degenerate.

Transport Consequences: The Einstein Relation

In non-degenerate semiconductor transport, carrier drift under an electric field $\mathbf{E}$ is balanced by carrier diffusion down a concentration gradient $\nabla n$.

Mathematical Derivation of $D/\mu$

In equilibrium under a potential gradient $\phi(x)$ (where $\mathbf{E} = -d\phi/dx$), the electron Fermi level remains spatially constant ($dE_F/dx = 0$). The conduction band edge varies as $E_c(x) = E_{c0} - q\phi(x)$.

The non-degenerate electron concentration is:

Taking the spatial gradient of $n(x)$:

Setting total electron current density $J_n = J_{\text{drift}} + J_{\text{diff}} = 0$:

Dividing by $q n \mathbf{E}$ yields the celebrated Einstein Relation for Non-Degenerate Carriers:

Where $V_t$ is the thermal voltage ($25.85\text{ mV}$ at $300\text{ K}$).

(Note: In degenerate semiconductors, the Einstein relation generalizes to $\frac{D_n}{\mu_n} = \frac{k_B T}{q} \frac{F_{1/2}(\eta_c)}{F_{-1/2}(\eta_c)}$, demonstrating that MB kinetics underestimates diffusion at high carrier densities).

Limits of Validity & Breakdown Regimes

Maxwell–Boltzmann statistics fails when system parameters cross quantum or non-equilibrium thresholds:

1. High Doping / Degeneracy Threshold ($n \ge 0.1 N_c$):

• In Si ($N_c = 2.86 \times 10^{19}\text{ cm}^{-3}$), when $N_D > 3 \times 10^{18}\text{ cm}^{-3}$, $E_F$ approaches within $2 k_B T$ of $E_c$.
• The MB approximation underestimates carrier density $n$ for a given $E_F$, requiring full Fermi–Dirac integrals.

2. Low-Temperature Carrier Freeze-Out ($T < 100\text{ K}$):

• At cryogenic temperatures, thermal energy $k_B T$ is smaller than donor/acceptor ionization energies ($E_d \approx 45\text{ meV}$ for P in Si). Dopants fail to ionize, invalidating simple MB thermal activation models.

3. High Field / Hot-Carrier Transport ($\mathbf{E} > 10^4\text{ V/cm}$):

• Under strong electric fields in sub-10 nm channel regions, carriers gain kinetic energy faster than they can relax via optical phonon emission.
• The electron energy distribution function (EEDF) becomes non-Maxwellian, developing a high-energy tail characterized by an elevated carrier temperature $T_e > T_{\text{lattice}}$, requiring numerical Boltzmann Transport Equation (BTE) or Monte Carlo solvers.

Quantitative Python Script: Velocity Distribution & Validity Check

The following Python script computes the 3D Maxwellian speed distribution for electrons in Si, GaAs, and InGaAs, and calculates the percentage error of the MB approximation vs exact FD statistics as a function of doping.

import numpy as np
import matplotlib.pyplot as plt

# Constants
k_B = 1.380649e-23       # J/K
k_B_eV = 8.617333e-5     # eV/K
m_0 = 9.1093837e-31      # kg
q = 1.602176634e-19      # C
T = 300.0                # K

# Effective masses (m* / m0)
m_star_Si = 0.26
m_star_GaAs = 0.067
m_star_InGaAs = 0.041

def maxwell_speed_pdf(v, m_eff):
    """3D Maxwell-Boltzmann speed PDF P(v)."""
    m = m_eff * m_0
    factor = 4.0 * np.pi * (m / (2.0 * np.pi * k_B * T))**(1.5)
    return factor * (v**2) * np.exp(-m * (v**2) / (2.0 * k_B * T))

v_grid = np.linspace(0, 1e6, 1000)  # m/s

# Calculate Characteristic Velocities for Silicon
m_Si = m_star_Si * m_0
v_p_Si = np.sqrt(2.0 * k_B * T / m_Si)
v_th_Si = np.sqrt(8.0 * k_B * T / (np.pi * m_Si))
v_rms_Si = np.sqrt(3.0 * k_B * T / m_Si)

print("==================================================================")
print("MAXWELL-BOLTZMANN ELECTRON THERMAL SPEEDS IN SILICON (T = 300 K)")
print("==================================================================")
print(f"Most Probable Velocity (v_p)   : {v_p_Si*1e-2:10.2f} cm/s ({v_p_Si:8.1f} m/s)")
print(f"Average Thermal Velocity (<v>) : {v_th_Si*1e-2:10.2f} cm/s ({v_th_Si:8.1f} m/s)")
print(f"RMS Velocity (v_rms)           : {v_rms_Si*1e-2:10.2f} cm/s ({v_rms_Si:8.1f} m/s)")
print("==================================================================")

# MB vs FD Validity Check
N_c_Si = 2.86e19  # cm^-3
n_ratios = np.logspace(-3, 1, 5)  # n / N_c from 0.001 to 10

print("\n==================================================================")
print("MB APPROXIMATION ACCURACY VS CARRIER DENSITY RATIO (n / N_c)")
print("==================================================================")
for r in n_ratios:
    n_conc = r * N_c_Si
    # MB assumption: eta_MB = ln(r)
    eta_MB = np.log(r)
    # Joyce-Dixon FD: eta_FD = ln(r) + r/sqrt(8)
    eta_FD = np.log(r) + r / np.sqrt(8.0)
    error_meV = (eta_FD - eta_MB) * (k_B_eV * T) * 1000.0
    status = "NON-DEGENERATE (MB Valid)" if r < 0.1 else "DEGENERATE (FD Required)"
    print(f"n/N_c = {r:6.3f} | n = {n_conc:8.2e} cm^-3 | Error: {error_meV:6.1f} meV | Status: {status}")
print("==================================================================")

Semiconductor Engineering Applications

1. Thermionic Emission in Schottky Barrier Diodes: The current density $J_{\text{SBD}}$ flowing over a metal-semiconductor Schottky barrier of height $\Phi_{Bn}$ is derived by integrating the MB velocity flux of carriers with kinetic energy exceeding $q\Phi_{Bn}$: $$J_{\text{SBD}} = A^ T^2 \exp\left( -\frac{q\Phi_{Bn}}{k_B T} \right) \left[ \exp\left( \frac{q V_a}{\eta k_B T} \right) - 1 \right]$$ Where $A^ = \frac{4\pi q m_n^* k_B^2}{h^3}$ is the effective Richardson constant ($112\text{ A/cm}^2\text{K}^2$ for n-type Si).

2. Subthreshold Leakage ($I_{\text{off}}$) in MOSFETs: In weak inversion ($V_{gs} < V_t$), the channel surface potential $\psi_s$ is controlled electrostatically by the gate. The subthreshold current is dominated by diffusion of non-degenerate electrons obeying MB statistics, producing the exponential subthreshold swing $S$: $$I_{\text{sub}} \propto \exp\left( \frac{q\psi_s}{k_B T} \right) \implies S = \frac{\partial V_{gs}}{\partial \log_{10} I_d} = \ln(10) \frac{k_B T}{q} \left( 1 + \frac{C_d}{C_{\text{ox}}} \right) \approx 60\text{ mV/dec at } 300\text{ K}$$

3. TCAD Drift-Diffusion Transport Solvers: Standard commercial TCAD tools (e.g., Synopsys Sentaurus, Silvaco Atlas) utilize MB carrier statistics as the baseline computationally efficient model for low-to-moderate doping regions, switching dynamically to Fermi–Dirac integrals in heavily doped source/drain regions.

References

1. Lundstrom, M. (2000). Fundamentals of Carrier Transport (2nd ed.). Cambridge University Press. 2. Sze, S. M., & Ng, K. K. (2006). Physics of Semiconductor Devices (3rd ed.). John Wiley & Sons. 3. Vasileska, D., Choudhury, S. M., & Goodnick, S. M. (2010). Computational Electronics: Semiclassical and Quantum Device Modeling and Simulation. CRC Press. 4. Markowich, P. A., Ringhofer, C. A., & Schmeiser, C. (1990). Semiconductor Equations. Springer-Verlag.