Level 3 — ProductionBand Theory & Carrier Physics

Band Theory & Carrier Physics

45 minutes60 marksprintable — key stays hidden on paper

Chapter: 2.1 Band Theory & Carrier Physics Level: 3 — Production (from-scratch derivations, code-from-memory, explain-out-loud) Time limit: 45 minutes Total marks: 60

Use kB=8.617×105 eV/Kk_B = 8.617\times10^{-5}\ \text{eV/K}, T=300 KT = 300\ \text{K} unless told otherwise so kBT0.02585 eVk_BT \approx 0.02585\ \text{eV}. Show all working.


Question 1 — Derive the mass action law (10 marks)

Starting from the equilibrium carrier concentration equations

n=NCe(ECEF)/kBT,p=NVe(EFEV)/kBTn = N_C\,e^{-(E_C - E_F)/k_BT}, \qquad p = N_V\,e^{-(E_F - E_V)/k_BT}

(a) Derive the mass action law np=ni2np = n_i^2 and give the explicit expression for nin_i in terms of NC,NV,EgN_C, N_V, E_g. (4)

(b) Explain out loud (in words) why the product npnp is independent of the Fermi level position, and what physically pins nin_i. (3)

(c) For an intrinsic semiconductor derive the expression for the intrinsic Fermi level EiE_i and state the condition under which EiE_i sits exactly at midgap. (3)


Question 2 — Band gaps and conductivity classification (8 marks)

(a) Contrast conductor, semiconductor, and insulator using band-diagram reasoning and typical band-gap magnitudes. Give one numeric example each. (4)

(b) A material's intrinsic carrier concentration doubling behaviour is governed by niT3/2eEg/2kBTn_i \propto T^{3/2}e^{-E_g/2k_BT}. Explain why insulators carry negligible current at room temperature while semiconductors do not, referencing the exponential term. (2)

(c) Direct vs indirect band gap: explain the difference in terms of crystal momentum kk, and state why this matters for light emission. (2)


Question 3 — Carrier concentration computation (12 marks)

Silicon at 300 K300\ \text{K} has ni=1.5×1010 cm3n_i = 1.5\times10^{10}\ \text{cm}^{-3}. It is doped with donor concentration ND=1×1016 cm3N_D = 1\times10^{16}\ \text{cm}^{-3} (fully ionised).

(a) Compute the majority electron concentration nn and minority hole concentration pp. (4)

(b) By how many orders of magnitude do majority and minority carriers differ? Comment on which dominates conduction. (3)

(c) The sample is now compensated by adding acceptors NA=6×1015 cm3N_A = 6\times10^{15}\ \text{cm}^{-3}. Recompute nn and pp using charge neutrality. (5)


Question 4 — Fermi–Dirac and Einstein relation (12 marks)

(a) Write the Fermi–Dirac distribution f(E)f(E). Show that for EEFkBTE - E_F \gg k_BT it reduces to the Boltzmann approximation, and state the value of f(E)f(E) at E=EFE = E_F. (4)

(b) Derive the Einstein relation Dμ=kBTq\dfrac{D}{\mu} = \dfrac{k_BT}{q} from the condition that in thermal equilibrium the total current (drift + diffusion) is zero in a non-uniformly doped semiconductor. Assume Boltzmann statistics. (6)

(c) Given electron mobility μn=1350 cm2/V⋅s\mu_n = 1350\ \text{cm}^2/\text{V·s} at 300 K300\ \text{K}, compute the diffusion coefficient DnD_n. (2)


Question 5 — Drift, diffusion, and code from memory (10 marks)

(a) Write the total electron current density equation combining drift and diffusion terms (1-D). Define each symbol. (3)

(b) Write, from memory, a short Python function (numpy allowed) that computes intrinsic carrier concentration given NCN_C, NVN_V, EgE_g, and TT. (4)

(c) Explain out loud what recombination and generation are, and how they balance in equilibrium. (3)


Question 6 — Temperature dependence (8 marks)

(a) Sketch (describe) the three regimes of majority carrier concentration vs temperature in a doped semiconductor: freeze-out, extrinsic/saturation, and intrinsic. (4)

(b) For silicon, using nieEg/2kBTn_i \propto e^{-E_g/2k_BT} with Eg=1.12 eVE_g = 1.12\ \text{eV}, compute the ratio ni(400)/ni(300)n_i(400)/n_i(300) (ignore the T3/2T^{3/2} prefactor). (4)

Answer keyMark scheme & solutions

Q1 (10)

(a) Multiply: np=NCNVe(ECEF)/kBTe(EFEV)/kBT=NCNVe(ECEV)/kBTnp = N_C N_V\, e^{-(E_C-E_F)/k_BT}\,e^{-(E_F-E_V)/k_BT} = N_C N_V\,e^{-(E_C-E_V)/k_BT} The EFE_F terms cancel, leaving np=NCNVeEg/kBTnp = N_C N_V e^{-E_g/k_BT} with Eg=ECEVE_g = E_C - E_V. (2) Since this is independent of EFE_F and equals ni2n_i^2 (definition, intrinsic where n=p=nin=p=n_i): ni=NCNVeEg/2kBTn_i = \sqrt{N_C N_V}\,e^{-E_g/2k_BT} (2)

(b) EFE_F appears with ++ in nn and - in pp; raising EFE_F increases nn but decreases pp by the same exponential factor, so the product cancels the EFE_F dependence (1). npnp therefore depends only on band structure (NC,NV,EgN_C,N_V,E_g) and TT (1). nin_i is pinned by the band gap and temperature — it is a material/temperature constant, not a doping quantity (1). (3)

(c) Set n=pn=p: NCe(ECEi)/kBT=NVe(EiEV)/kBTN_C e^{-(E_C-E_i)/k_BT} = N_V e^{-(E_i-E_V)/k_BT}. Take logs and solve: Ei=EC+EV2+kBT2ln ⁣NVNCE_i = \frac{E_C+E_V}{2} + \frac{k_BT}{2}\ln\!\frac{N_V}{N_C} (2) EiE_i is at midgap exactly when NC=NVN_C = N_V (log term vanishes). (1)


Q2 (8)

(a) Conductor: conduction and valence bands overlap (or partially filled band), Eg0E_g \approx 0, e.g. copper (1). Semiconductor: small gap, Eg0.5E_g \approx 0.52 eV2\ \text{eV}, e.g. Si 1.12 eV1.12\ \text{eV} (1). Insulator: large gap Eg>5 eVE_g > 5\ \text{eV}, e.g. diamond 5.5 eV\sim5.5\ \text{eV} or SiO₂ 9 eV\sim9\ \text{eV} (1). Conduction requires electrons in conduction band / holes in valence band; large gap ⇒ negligible thermal excitation (1). (4)

(b) The factor eEg/2kBTe^{-E_g/2k_BT} is exponentially small when Eg2kBTE_g \gg 2k_BT. At 300 K, 2kBT0.052 eV2k_BT\approx0.052\ \text{eV}; for insulators Eg/2kBTE_g/2k_BT is huge, giving essentially zero free carriers; for semiconductors the exponent is modest so appreciable nin_i exists. (2)

(c) Direct gap: conduction band minimum and valence band maximum at the same kk — electron–hole recombination emits a photon directly (momentum conserved). Indirect gap: minima at different kk — recombination needs a phonon to conserve momentum, so it is inefficient for light emission. This is why LEDs/lasers use direct-gap materials (GaAs) not Si. (2)


Q3 (12)

(a) nND=1×1016 cm3n \approx N_D = 1\times10^{16}\ \text{cm}^{-3} (donor >> nin_i) (2) p=ni2/n=(1.5×1010)2/1016=2.25×1020/1016=2.25×104 cm3p = n_i^2/n = (1.5\times10^{10})^2/10^{16} = 2.25\times10^{20}/10^{16} = 2.25\times10^{4}\ \text{cm}^{-3} (2)

(b) Ratio =1016/2.25×1044.4×1011= 10^{16}/2.25\times10^{4} \approx 4.4\times10^{11}, about 11–12 orders of magnitude (1). Electrons (majority) dominate conduction overwhelmingly (2). (3)

(c) Net doping NDNA=10166×1015=4×1015 cm3N_D - N_A = 10^{16} - 6\times10^{15} = 4\times10^{15}\ \text{cm}^{-3} (n-type) (2). n4×1015 cm3n \approx 4\times10^{15}\ \text{cm}^{-3} (1). p=ni2/n=2.25×1020/4×1015=5.625×104 cm3p = n_i^2/n = 2.25\times10^{20}/4\times10^{15} = 5.625\times10^{4}\ \text{cm}^{-3} (2). (5)


Q4 (12)

(a) f(E)=11+e(EEF)/kBTf(E) = \dfrac{1}{1+e^{(E-E_F)/k_BT}} (1) For EEFkBTE-E_F \gg k_BT, the exponential dominates: f(E)e(EEF)/kBTf(E)\approx e^{-(E-E_F)/k_BT} (Boltzmann tail) (2). At E=EFE = E_F: exponent =0=0, so f=1/(1+1)=1/2f = 1/(1+1) = 1/2 (1). (4)

(b) Equilibrium electron current Jn=qμnnE+qDndndx=0J_n = q\mu_n n E + qD_n \dfrac{dn}{dx} = 0 (1). With built-in field E=1qdEF...E = -\dfrac{1}{q}\dfrac{dE_F... }{} — use n=nie(EFEi)/kBTn = n_i e^{(E_F-E_i)/k_BT}; in equilibrium EFE_F flat, non-uniform EiE_i. Field from E=dϕdxE = -\dfrac{d\phi}{dx}, and dndx=nkBTd(EFEi)dx\dfrac{dn}{dx} = \dfrac{n}{k_BT}\dfrac{d(E_F-E_i)}{dx}. With qE=dEidxqE = \dfrac{dE_i}{dx} (band bending) (2): qμnnE=qDnnkBTqE    μn=qDnkBTq\mu_n n E = qD_n\frac{n}{k_BT}\,qE \;\Rightarrow\; \mu_n = \frac{qD_n}{k_BT} Hence Dnμn=kBTq\dfrac{D_n}{\mu_n} = \dfrac{k_BT}{q} (3). (6)

(c) Dn=μnkBTq=1350×0.02585=34.9 cm2/sD_n = \mu_n \dfrac{k_BT}{q} = 1350 \times 0.02585 = 34.9\ \text{cm}^2/\text{s} (since kBT/q=0.02585 Vk_BT/q = 0.02585\ \text{V}). (2)


Q5 (10)

(a) Jn=qμnnE+qDndndxJ_n = q\mu_n n E + qD_n\dfrac{dn}{dx}. qq electron charge magnitude, μn\mu_n mobility, nn electron density, EE field (drift term), DnD_n diffusion coefficient, dn/dxdn/dx concentration gradient (diffusion term). (3)

(b) Sample answer (award for correct exponential form, sqrt, and units handling):

import numpy as np
def n_intrinsic(Nc, Nv, Eg, T):
    kB = 8.617e-5        # eV/K
    return np.sqrt(Nc*Nv)*np.exp(-Eg/(2*kB*T))

(4) — 2 for sqrt(Nc*Nv), 2 for exp(-Eg/(2*kB*T)).

(c) Generation: thermal/optical energy creates electron–hole pairs (excites electrons across gap) (1). Recombination: an electron drops from CB to VB annihilating a hole, releasing energy (photon/phonon) (1). In equilibrium generation rate = recombination rate, keeping np=ni2np = n_i^2 constant (1). (3)


Q6 (8)

(a) Freeze-out (low T): dopants not fully ionised, nn rises with TT (1). Extrinsic/saturation (mid T): all dopants ionised, nNDn \approx N_D roughly flat (1). Intrinsic (high T): nin_i exceeds NDN_D, nn rises steeply/exponentially as intrinsic generation dominates (1). Correct ordering & shape (1). (4)

(b) ni(400)ni(300)=exp ⁣[Eg2kB(14001300)]\dfrac{n_i(400)}{n_i(300)} = \exp\!\left[-\dfrac{E_g}{2k_B}\left(\dfrac1{400}-\dfrac1{300}\right)\right]. Eg2kB=1.122×8.617×105=6499 K-\dfrac{E_g}{2k_B} = -\dfrac{1.12}{2\times8.617\times10^{-5}} = -6499\ \text{K} (1). (14001300)=8.333×104 K1\left(\dfrac1{400}-\dfrac1{300}\right) = -8.333\times10^{-4}\ \text{K}^{-1} (1). Exponent =6499×(8.333×104)=5.416= -6499 \times (-8.333\times10^{-4}) = 5.416 (1). Ratio =e5.416225= e^{5.416} \approx 225 (1). (4)


[
{"claim":"Q3a minority hole p = 2.25e4 cm^-3","code":"ni=1.5e10; ND=1e16; p=ni**2/ND; result = abs(p-2.25e4)<1"},
{"claim":"Q3c compensated: n=4e15, p=5.625e4","code":"ni=1.5e10; net=1e16-6e15; n=net; p=ni**2/n; result = abs(n-4e15)<1e6 and abs(p-5.625e4)<1"},
{"claim":"Q4c D_n = 34.9 cm^2/s","code":"mu=1350; Vt=0.02585; Dn=mu*Vt; result = abs(Dn-34.9)<0.2"},
{"claim":"Q6b ni ratio ~225","code":"import sympy as sp; Eg=sp.Rational(112,100); kB=8.617e-5; exp_arg=-Eg/(2*kB)*(sp.Rational(1,400)-sp.Rational(1,300)); ratio=sp.exp(exp_arg); result = abs(float(ratio)-225)<8"}
]