When you are staring at a solid-state physics midterm or trying to spec a custom epitaxial layer for a power electronics project, calculating the conductivity of a doped example of semiconductor material is a foundational skill. Silicon remains the most ubiquitous example of semiconductor material in commercial use, but the mathematical framework applies universally across the periodic table's Group IV and III-V compounds.
This walkthrough dissects a classic exam problem: finding the room-temperature conductivity and drift current density of an N-type silicon wafer. We will cover the governing theorems, expose the most common unit-conversion traps, and execute every algebraic step without skipping the "obvious" arithmetic that costs students points.
Reference Data: Material Properties of Common Semiconductor Examples
Before solving any drift-diffusion problem, you need reliable baseline constants. The table below provides the room-temperature (300 K) properties for four standard materials. Note that while intrinsic carrier concentration ($n_i$) values are heavily debated in modern literature (with silicon's $n_i$ often cited as $1.0 \times 10^{10}$ cm$^{-3}$ in newer texts versus the classic $1.5 \times 10^{10}$ cm$^{-3}$), we will use the modern $1.0 \times 10^{10}$ cm$^{-3}$ consensus for this walkthrough.
| Material | Bandgap ($E_g$) at 300K | Intrinsic $n_i$ (cm$^{-3}$) | Electron Mobility $\mu_n$ (cm$^2$/V·s) | Hole Mobility $\mu_p$ (cm$^2$/V·s) |
|---|---|---|---|---|
| Silicon (Si) | 1.12 eV | $1.0 \times 10^{10}$ | 1350 | 480 |
| Germanium (Ge) | 0.66 eV | $2.4 \times 10^{13}$ | 3900 | 1900 |
| Gallium Arsenide (GaAs) | 1.42 eV | $2.1 \times 10^{6}$ | 8500 | 400 |
| Silicon Carbide (4H-SiC) | 3.26 eV | $8.2 \times 10^{-9}$ | 900 | 100 |
For deeper reading on how temperature affects these mobility figures, the All About Circuits semiconductor textbook provides excellent open-source derivations of the phonon scattering models that govern these limits.
The Exam Problem: N-Type Silicon Conductivity & Drift Current
A specific example of semiconductor material—an N-type Silicon wafer—is uniformly doped with $10^{16}$ phosphorus atoms/cm$^3$. The wafer is maintained at 300 K and subjected to a uniform electric field of $50 \text{ V/cm}$. Calculate:
1. The majority and minority carrier concentrations.
2. The total electrical conductivity ($\sigma$) of the wafer.
3. The drift current density ($J$) flowing through the material.
Which theorem applies? We use the Mass Action Law ($np = n_i^2$) to find minority carriers, and the Drift Current Equation ($J = \sigma E$, where $\sigma = q(n\mu_n + p\mu_p)$) for transport.
The Trap: The most common mistake in this specific example of semiconductor calculation is unit mismatching. Students often convert concentrations to $m^{-3}$ but forget to convert mobility from $\text{cm}^2/\text{V}\cdot\text{s}$ to $\text{m}^2/\text{V}\cdot\text{s}$, or they mix V/cm with V/m. Rule of thumb: In semiconductor device physics, stick entirely to the CGS-derived standard (cm, V/cm, A/cm$^2$) unless explicitly forced into SI meters. It saves you from $10^{22}$ arithmetic errors.
Step-by-Step Algebraic Solution
We will solve this using the standard electron charge constant $q = 1.6 \times 10^{-19} \text{ C}$ and the Silicon parameters from our reference table.
-
Determine Carrier Concentrations ($n$ and $p$)
Phosphorus is a Group V donor, so the doping concentration $N_d = 10^{16} \text{ cm}^{-3}$.
Because $N_d \gg n_i$ ($10^{16} \gg 10^{10}$), the majority electron concentration is simply the doping level:
n \approx N_d = 10^{16} \text{ cm}^{-3}
To find the minority hole concentration ($p$), apply the Mass Action Law ($np = n_i^2$):
p = n_i^2 / n
p = (1.0 \times 10^{10})^2 / 10^{16}
p = 1.0 \times 10^{20} / 10^{16}
p = 1.0 \times 10^{4} \text{ cm}^{-3} -
Calculate Total Conductivity ($\sigma$)
The general conductivity formula is $\sigma = q(n\mu_n + p\mu_p)$.
Plug in the values:
\sigma = (1.6 \times 10^{-19}) \times [(10^{16} \times 1350) + (10^{4} \times 480)]
Evaluate the bracketed terms (notice how the minority hole term is mathematically negligible, but we show the work for completeness):
\text{Electron term} = 1.35 \times 10^{19}
\text{Hole term} = 4.8 \times 10^{6}
Add them together:
\sigma = (1.6 \times 10^{-19}) \times [1.35 \times 10^{19} + 0.00000048 \times 10^{19}]
\sigma \approx (1.6 \times 10^{-19}) \times (1.35 \times 10^{19})
\sigma = 1.6 \times 1.35 = 2.16 \text{ S/cm} -
Calculate Drift Current Density ($J$)
Apply Ohm's Law in point form: $J = \sigma E$.
J = 2.16 \text{ S/cm} \times 50 \text{ V/cm}
J = 108 \text{ A/cm}^2
Units: Conductivity is Siemens per cm (S/cm). Current density is (S/cm) $\times$ (V/cm) = A/cm$^2$. The units balance perfectly.
Order of Magnitude: A conductivity of 2.16 S/cm equates to a resistivity ($\rho = 1/\sigma$) of $0.46 \text{ }\Omega\cdot\text{cm}$. If you cross-reference standard semiconductor resistivity vs. doping charts, an N-type Si wafer at $10^{16} \text{ cm}^{-3}$ doping yields roughly $0.4$ to $0.5 \text{ }\Omega\cdot\text{cm}$. Our math aligns perfectly with empirical bench measurements.
Independent Verification & Common FAQs
How to verify the answer independently:
If you are on the bench and want to verify this theoretical example of semiconductor behavior, you cannot measure current density directly without knowing the cross-sectional area. Instead, use a four-point probe station to measure the sheet resistance of the wafer. By multiplying the sheet resistance by the known epitaxial layer thickness, you extract the bulk resistivity ($\rho$). Invert that value to get $\sigma$, and compare it to our calculated $2.16 \text{ S/cm}$. If your probe reads $0.46 \text{ }\Omega\cdot\text{cm}$, your theoretical doping assumptions match the physical silicon.
Frequently Asked Questions
- Why did we ignore the hole mobility in the final calculation?
We didn't ignore it algebraically; we proved it was negligible. The minority hole concentration ($10^4$) is twelve orders of magnitude smaller than the majority electron concentration ($10^{16}$). Multiplying it by $\mu_p$ yields a value that vanishes past the 8th decimal place in the final sum. - What if the temperature was 400 K instead of 300 K?
The intrinsic carrier concentration ($n_i$) scales exponentially with temperature. At 400 K, $n_i$ increases drastically, and the mobility ($\mu_n$) decreases due to increased lattice (phonon) scattering. You would need to recalculate $n_i$ using the Boltzmann distribution and apply temperature-derating curves to the mobility before running the algebra again. - Does this apply to wide-bandgap materials like SiC?
Yes, the algebraic framework is identical. However, as seen in our table, the intrinsic carrier concentration for 4H-SiC is practically zero at room temperature, and the mobility figures are highly anisotropic depending on the crystal axis relative to the electric field.
