2.1.6Band Theory & Carrier Physics

Carrier concentration equations (n, p, ni)

1,782 words8 min readdifficulty · medium1 backlinks

1. From first principles: what are we integrating?

WHY the limits? Electrons only conduct if they sit above EcE_c (conduction band). Holes are missing electrons below EvE_v (valence band). So each integral covers only its own band.


2. Deriving the density of states gc(E)g_c(E)

Near the band edge the electron behaves like a free particle with an effective mass mnm_n^*: EEc=2k22mnE - E_c = \frac{\hbar^2 k^2}{2 m_n^*}

Why? The parabolic approximation: the bottom of a band always looks like a parabola, and curvature \to effective mass.

Count kk-states in a sphere of radius kk (spin factor 2, volume VV): N(k)=2V(2π)343πk3N(k) = 2 \cdot \frac{V}{(2\pi)^3}\cdot \frac{4}{3}\pi k^3

Convert to energy using k=2mn(EEc)/k=\sqrt{2m_n^*(E-E_c)}/\hbar and differentiate gc(E)=1VdNdEg_c(E)=\frac{1}{V}\frac{dN}{dE}:

gc(E)=12π2(2mn2)3/2EEc\boxed{g_c(E) = \frac{1}{2\pi^2}\left(\frac{2m_n^*}{\hbar^2}\right)^{3/2}\sqrt{E-E_c}}


3. The Boltzmann approximation (WHY it simplifies everything)

For a non-degenerate semiconductor, EFE_F sits deep in the gap, so for EEcE \ge E_c we have EEFkTE-E_F \gg kT: f(E)=11+e(EEF)/kTe(EEF)/kTf(E)=\frac{1}{1+e^{(E-E_F)/kT}} \approx e^{-(E-E_F)/kT}


4. The master result: nn and NcN_c

Plug gcg_c and Boltzmann ff into the integral: n=Ec12π2(2mn2)3/2EEc  e(EEF)/kTdEn = \int_{E_c}^{\infty}\frac{1}{2\pi^2}\left(\frac{2m_n^*}{\hbar^2}\right)^{3/2}\sqrt{E-E_c}\;e^{-(E-E_F)/kT}dE

Substitute x=(EEc)/kTx=(E-E_c)/kT, use 0xexdx=π2\int_0^\infty \sqrt{x}\,e^{-x}dx=\tfrac{\sqrt\pi}{2}:

n=Nce(EcEF)/kT,Nc=2(2πmnkTh2)3/2\boxed{n = N_c\, e^{-(E_c-E_F)/kT}}, \qquad N_c = 2\left(\frac{2\pi m_n^* kT}{h^2}\right)^{3/2}

By the identical argument for holes (using 1fe(EFEv)/kT1-f\approx e^{-(E_F-E_v)/kT}):

p=Nve(EFEv)/kT,Nv=2(2πmpkTh2)3/2\boxed{p = N_v\, e^{-(E_F-E_v)/kT}}, \qquad N_v = 2\left(\frac{2\pi m_p^* kT}{h^2}\right)^{3/2}

Figure — Carrier concentration equations (n, p, ni)

5. The intrinsic concentration nin_i and the mass-action law

Multiply nn and pp — the EFE_F terms cancel: np=NcNve(EcEv)/kT=NcNveEg/kTnp = N_c N_v\, e^{-(E_c-E_v)/kT} = N_c N_v\, e^{-E_g/kT}

This is independent of EFE_F → true for any doping. In an intrinsic crystal n=pnin=p\equiv n_i, so:

ni=NcNv  eEg/2kTnp=ni2\boxed{n_i = \sqrt{N_c N_v}\; e^{-E_g/2kT}} \quad\Rightarrow\quad \boxed{np = n_i^2}

Location of intrinsic level EiE_i (set n=pn=p): Ei=Ec+Ev2+kT2ln ⁣NvNcE_i = \frac{E_c+E_v}{2} + \frac{kT}{2}\ln\!\frac{N_v}{N_c} So EiE_i sits near midgap, nudged by the small ln(Nv/Nc)\ln(N_v/N_c) term.

Rewriting n,pn,p around EiE_i (very handy form): n=nie(EFEi)/kT,p=nie(EiEF)/kTn = n_i\, e^{(E_F-E_i)/kT}, \qquad p = n_i\, e^{(E_i-E_F)/kT}


6. Worked examples



7. Feynman + Mnemonic

Recall Explain to a 12-year-old (click to reveal)

Imagine a theater. The balcony (conduction band) has empty seats; the main floor (valence band) is packed. Heat gives a few floor-people enough energy to jump to the balcony — each leaves an empty seat (a "hole") downstairs. NcN_c is how many balcony seats there really are, and the exponential eEg/2kTe^{-E_g/2kT} is how likely someone can afford the jump. The taller the jump (EgE_g), the fewer make it; the hotter the room (TT), the more do. And magically, (people upstairs) × (empty seats downstairs) always equals a fixed number ni2n_i^2.


8. Forecast-then-Verify checkpoints


Flashcards

What integral defines nn?
n=Ecgc(E)f(E)dEn=\int_{E_c}^{\infty} g_c(E)f(E)\,dE — density of states × occupancy over the conduction band.
Why does f(E)e(EEF)/kTf(E)\to e^{-(E-E_F)/kT}?
For EEFkTE-E_F\gg kT (non-degenerate), the "+1" in Fermi–Dirac is negligible → Boltzmann tail.
State nn in terms of NcN_c.
n=Nce(EcEF)/kTn=N_c\,e^{-(E_c-E_F)/kT}.
What is NcN_c physically?
Effective density of states at the conduction edge, Nc=2(2πmnkT/h2)3/2N_c=2(2\pi m_n^*kT/h^2)^{3/2}.
Derive the mass-action law.
np=NcNve(EcEv)/kT=ni2np=N_cN_v e^{-(E_c-E_v)/kT}=n_i^2; EFE_F cancels, so it holds for any doping.
Formula for nin_i.
ni=NcNveEg/2kTn_i=\sqrt{N_cN_v}\,e^{-E_g/2kT}.
Why Eg/2E_g/2 not EgE_g in nin_i?
Because ni=npn_i=\sqrt{np}, and npeEg/kTnp\propto e^{-E_g/kT}; the square root halves the exponent.
Compact forms around EiE_i.
n=nie(EFEi)/kTn=n_i e^{(E_F-E_i)/kT}, p=nie(EiEF)/kTp=n_i e^{(E_i-E_F)/kT}.
In n-type with NDniN_D\gg n_i, find pp.
nNDn\approx N_D, then p=ni2/NDp=n_i^2/N_D.
Does doping change nin_i?
No — nin_i depends only on Nc,Nv,Eg,TN_c,N_v,E_g,T.
Where does EiE_i sit?
Near midgap: Ei=Ec+Ev2+kT2ln(Nv/Nc)E_i=\frac{E_c+E_v}{2}+\frac{kT}{2}\ln(N_v/N_c).
Why does gcEEcg_c\propto\sqrt{E-E_c}?
Counting kk-states on spheres; energy k2\propto k^2 gives density EEc\propto\sqrt{E-E_c}.

Connections

  • Fermi-Dirac distribution — the occupancy factor f(E)f(E) we approximated.
  • Density of states — origin of gc(E)EEcg_c(E)\propto\sqrt{E-E_c}.
  • Effective mass — sets mn,mpm_n^*, m_p^* inside Nc,NvN_c, N_v.
  • Doping and charge neutrality — combines with np=ni2np=n_i^2 to solve for n,pn,p.
  • Fermi level position — computed from nn vs NcN_c.
  • Intrinsic vs extrinsic semiconductors — regimes where these equations apply.
  • Band gap Eg — controls the exponential in nin_i.

Concept Map

seats × occupancy

occupancy

integrate over band

integrate over band

derives

non-degenerate case

requires Ec-EF ≳ 3kT

yields

compact form

if violated

n·p product

Density of states g(E)

Carrier count n, p

Fermi-Dirac f(E)

Electron conc n

Hole conc p uses 1-f

Parabolic band + effective mass

Boltzmann approx

Non-degenerate check

Effective DOS Nc, Nv

n = Nc·exp(-(Ec-EF)/kT)

Use Fermi-Dirac integral F1/2

Intrinsic conc ni

Hinglish (regional understanding)

Intuition Hinglish mein samjho

Dekho, semiconductor mein hume count karna hota hai ki kitne mobile electrons conduction band mein hain (nn) aur kitne holes valence band mein hain (pp). Directly gin nahi sakte, isliye do cheezein multiply karte hain: kitni seats available hain (density of states g(E)g(E)) aur har seat filled hone ka chance kitna hai (Fermi-Dirac f(E)f(E)). Inko band ke upar integrate karo, bas — carrier concentration mil gaya.

Non-degenerate case mein f(E)f(E) ka "+1" ignore karke Boltzmann tail e(EEF)/kTe^{-(E-E_F)/kT} use karte hain, aur square-root DOS integrate karne par simple result aata hai: n=Nce(EcEF)/kTn=N_c e^{-(E_c-E_F)/kT} aur p=Nve(EFEv)/kTp=N_v e^{-(E_F-E_v)/kT}. Yahan Nc,NvN_c, N_v ko "effective number of states" samjho — jaise poora band ek hi level pe collapse ho gaya ho.

Sabse important jugaad: nn aur pp ko multiply karo to EFE_F cancel ho jaata hai, aur milta hai np=NcNveEg/kT=ni2np=N_cN_v e^{-E_g/kT}=n_i^2. Yeh mass-action law hai. Intrinsic mein n=p=ni=NcNveEg/2kTn=p=n_i=\sqrt{N_cN_v}\,e^{-E_g/2kT}. Yaad rakhna Eg/2E_g/2 aata hai kyunki square root hai. Ek aur baat: nin_i sirf material aur temperature pe depend karta hai — doping se nin_i kabhi nahi badalta, sirf nn aur pp ka balance shift hota hai.

Practical trick: n-type mein NDniN_D\gg n_i ho to nNDn\approx N_D, phir p=ni2/NDp=n_i^2/N_D se minority carrier nikaal lo. Bas yehi 20% concepts 80% numericals solve kar dete hain, isliye in teen equations ko rat mat lo — derive karke feel lo.

Go deeper — visual, from zero

Test yourself — Band Theory & Carrier Physics

Connections