2.1.4Band Theory & Carrier Physics

Fermi level and Fermi-Dirac distribution

1,905 words9 min readdifficulty · medium6 backlinks

WHY do we even need this?

WHAT problem are we solving? We have billions of electrons and a continuum of allowed energy states (bands). We can't track each electron. We need a statistical rule: given a state at energy EE, what's the chance an electron sits there?

WHY not just use classical Boltzmann statistics? Because electrons obey the Pauli exclusion principle: at most one electron per quantum state. Classical particles are distinguishable and can share states freely. That single constraint changes the whole distribution.


Deriving the Fermi–Dirac distribution from scratch

We derive f(E)f(E) using the grand-canonical idea: a single state exchanges electrons with a big reservoir at temperature TT and chemical potential μ\mu.

Step 1 — A state can only hold 0 or 1 electron. Why this step? Pauli exclusion. Occupation number n{0,1}n \in \{0, 1\} only.

Step 2 — Grand-canonical probability of each option. The probability of a state having energy EE and nn electrons is proportional to the Gibbs factor e(Enμn)/kBTe^{-(E n - \mu n)/k_BT}. Why this step? μ\mu (chemical potential) is the energy cost/gain of adding one particle from the reservoir; the μn-\mu n term accounts for particle exchange.

  • n=0n=0: weight =e0=1= e^{0} = 1
  • n=1n=1: weight =e(Eμ)/kBT= e^{-(E-\mu)/k_BT}

Step 3 — Normalize to get the average occupation.

f(E)=n=01+1e(Eμ)/kBT1+e(Eμ)/kBTf(E) = \langle n\rangle = \frac{0\cdot 1 + 1\cdot e^{-(E-\mu)/k_BT}}{1 + e^{-(E-\mu)/k_BT}}

Why this step? Average occupation = (sum of n×n\timesweight) / (sum of weights).

Step 4 — Divide top and bottom by the numerator's exponential to get the clean form:

f(E)=11+e(Eμ)/kBT\boxed{\,f(E) = \frac{1}{1 + e^{(E-\mu)/k_BT}}\,}

At T=0T=0, μEF\mu \to E_F, and even at finite TT in semiconductors we usually call μEF\mu \approx E_F. So:

Figure — Fermi level and Fermi-Dirac distribution

HOW f(E)f(E) behaves — reading the curve

At T=0T=0 K (the step function):

  • E<EFE < E_F: e(EEF)/kBTe=0f=1e^{(E-E_F)/k_BT} \to e^{-\infty} = 0 \Rightarrow f=1 (all states below filled).
  • E>EFE > E_F: e+f=0e^{+\infty}\to\infty \Rightarrow f=0 (all states above empty).

At T>0T>0: the sharp step smears over a width of a few kBTk_BT. Some electrons get thermally kicked above EFE_F, leaving empty states (holes) below.

Boltzmann approximation (the 80/20 shortcut)

When EEFkBTE - E_F \gg k_BT (states well above EFE_F, e.g. conduction band in a semiconductor), the 11 in the denominator is negligible:

f(E)e(EEF)/kBTf(E) \approx e^{-(E-E_F)/k_BT}

Why this matters: this is the classical Boltzmann tail — it's why we can compute carrier concentrations easily in non-degenerate semiconductors. This single approximation covers ~80% of device calculations.


Worked examples


Common mistakes (Steel-manned)


Recall Feynman: explain to a 12-year-old

Imagine a stadium with seats stacked from ground floor upward, and a rule: only one person per seat. People are lazy and want the lowest seats. When it's freezing cold (T=0T=0), everyone sits as low as possible — every low seat taken, every high seat empty. The "water line" between full and empty seats is the Fermi level. When it warms up, a few people near the line get excited and jump to higher seats, leaving some low seats empty. The Fermi–Dirac curve just tells you the chance any given seat has someone in it — 1 way below the line, 0 way above, and a fuzzy 12\tfrac12 right at the line.


Active-recall flashcards

What is the Fermi–Dirac distribution formula?
f(E)=11+e(EEF)/kBTf(E)=\dfrac{1}{1+e^{(E-E_F)/k_BT}}
By definition, what is f(E)f(E) at E=EFE=E_F?
Exactly 12\tfrac12 (50%), at any temperature.
What is the Fermi level physically?
The electrochemical potential of electrons; energy where occupation probability = 1/2; constant across a system in equilibrium.
Why is Fermi–Dirac used instead of Boltzmann for electrons?
Electrons obey the Pauli exclusion principle (max one per state), so occupation caps at 1; Boltzmann has no such cap.
What does f(E)f(E) look like at T=0T=0 K?
A step function: f=1f=1 for E<EFE<E_F, f=0f=0 for E>EFE>E_F.
What is the Boltzmann approximation of f(E)f(E) and when valid?
fe(EEF)/kBTf\approx e^{-(E-E_F)/k_BT}, valid when EEFkBTE-E_F \gg k_BT (non-degenerate, far above EFE_F).
Probability a state Δ\Delta below EFE_F is EMPTY equals what?
f(EF+Δ)f(E_F+\Delta) — same as probability a state Δ\Delta above is filled (antisymmetry: f(EF+Δ)+f(EFΔ)=1f(E_F+\Delta)+f(E_F-\Delta)=1).
Value of kBTk_BT at room temperature (300 K)?
About 0.0259 eV (≈ 26 meV).
Does f(E)f(E) give the number of electrons?
No — it's a probability (0–1). Electron count needs n=g(E)f(E)dEn=\int g(E)f(E)\,dE.
Inversion: given ff, what is EEFE-E_F?
EEF=kBTln ⁣(1ff)E-E_F = k_BT\ln\!\big(\tfrac{1-f}{f}\big).

Connections

Concept Map

forces

rules out

input to

enters

normalize avg occupation

defines

is where

equals at T=0

flat in equilibrium

at T=0 becomes

smears

Pauli exclusion principle

Occupation n 0 or 1

Classical Boltzmann stats

Grand-canonical Gibbs factor

Chemical potential mu

Fermi-Dirac distribution f of E

Fermi level E_F

f equals one half

Constant across system

Step function

Temperature T

Hinglish (regional understanding)

Intuition Hinglish mein samjho

Dekho, electrons thode "antisocial" hote hain — Pauli exclusion principle kehta hai ki ek quantum state mein bas ek hi electron beth sakta hai. Isliye saare electrons sabse neeche wale level mein nahi ghus sakte; woh neeche se upar tak seats bharte jaate hain. Fermi–Dirac distribution f(E)f(E) humein batati hai ki kisi energy EE wali seat ke bharne ki probability kitni hai. Formula hai f(E)=11+e(EEF)/kBTf(E)=\frac{1}{1+e^{(E-E_F)/k_BT}}, aur EFE_F (Fermi level) woh energy hai jahan yeh probability exactly aadhi (0.5) ho jaati hai.

Zero temperature pe yeh ekdum step ban jaati hai — EFE_F ke neeche sab full (f=1), upar sab empty (f=0). Jab temperature badhta hai, yeh step "smear" ho jaata hai thoda — kuch electrons garmi se upar chhalaang maar dete hain aur neeche khaali jagah (holes) chhod dete hain. Yeh smearing sirf kuch kBTk_BT jitni chaudi hoti hai (room temp pe kBT0.026k_BT \approx 0.026 eV).

Important baat: f(E)f(E) sirf probability hai, electron count nahi. Agar us energy pe koi state hi nahi hai (jaise semiconductor ke band gap mein), toh f=0.5f=0.5 hone ke bawajood wahan koi electron nahi hoga. Actual electron number nikalne ke liye density of states g(E)g(E) se multiply karke integrate karna padta hai.

Yeh kyun important hai? Kyunki har semiconductor device — diode, transistor, solar cell — ki working depend karti hai ki EFE_F kahan hai aur bands kitne bhare hain. Doping se EFE_F shift hota hai, aur wahi decide karta hai current kaise flow karega. Ek 80/20 shortcut yaad rakho: jab EEFkBTE-E_F \gg k_BT, toh fe(EEF)/kBTf \approx e^{-(E-E_F)/k_BT} — yehi Boltzmann approximation zyaadatar calculations mein kaam aata hai.

Go deeper — visual, from zero

Test yourself — Band Theory & Carrier Physics

Connections