2.4 Equilibrium carrier concentrations
(a) Electrons in the conduction band
The number of conduction electrons per unit volume with energies between \(E\) and \(E+dE\) is \(dn = g_c(E)f(E)\,dE\), so that the equilibrium electron concentration is
\[n_0 = \int_{E_C}^{\infty} g_c(E)\,f(E)\,dE ,\]
where the upper limit may be taken as infinity because \(f(E)\) falls off so rapidly that states high in the band contribute nothing. In a non-degenerate semiconductor the Fermi level lies in the gap at least a few \(k_BT\) below \(E_C\), so that \(E-E_F\gg k_BT\) for every \(E\ge E_C\) and the Fermi function reduces to the Boltzmann factor \(f(E)\simeq e^{-(E-E_F)/k_BT}\). Hence
\[n_0 = C_n\int_{E_C}^{\infty}(E-E_C)^{1/2}\,e^{-(E-E_F)/k_BT}\,dE .\]
We write \(E-E_F=(E-E_C)+(E_C-E_F)\) and substitute \(x=(E-E_C)/k_BT\), so that \(dE = k_BT\,dx\):
\[n_0 = C_n\,(k_BT)^{3/2}\,e^{-(E_C-E_F)/k_BT}\int_0^\infty x^{1/2}e^{-x}\,dx = C_n\,(k_BT)^{3/2}\,\frac{\sqrt\pi}{2}\,e^{-(E_C-E_F)/k_BT},\]
since \(\int_0^\infty x^{1/2}e^{-x}dx=\Gamma(3/2)=\sqrt\pi/2\). Inserting \(C_n = 4\pi(2m_e^*)^{3/2}/h^3\), the prefactor becomes
\[\frac{4\pi(2m_e^*)^{3/2}}{h^3}\,(k_BT)^{3/2}\,\frac{\sqrt\pi}{2} = \frac{2\,\pi^{3/2}(2m_e^*k_BT)^{3/2}}{h^3} = 2\left(\frac{2\pi m_e^*k_BT}{h^2}\right)^{3/2}.\]
We therefore obtain
\[\boxed{n_0 = N_C\,e^{-(E_C-E_F)/k_BT}},\qquad \boxed{N_C = 2\left(\frac{2\pi m_e^*k_BT}{h^2}\right)^{3/2}}.\]
The quantity \(N_C\), which varies as \(T^{3/2}\), is called the effective density of states at the conduction-band edge. The result shows that the conduction band behaves as if all its states, \(N_C\) per unit volume, were concentrated at the single energy \(E_C\) and occupied with the Boltzmann probability \(e^{-(E_C-E_F)/k_BT}\).
(b) Holes in the valence band
The probability that a state is empty, that is, occupied by a hole, is
\[1-f(E) = \frac{e^{(E-E_F)/k_BT}}{1+e^{(E-E_F)/k_BT}} \simeq e^{-(E_F-E)/k_BT},\]
the approximation holding when \(E_F-E\gg k_BT\) for all \(E\le E_V\). With \(g_v(E)=C_p(E_V-E)^{1/2}\) the hole concentration is
\[p_0 = C_p\int_{-\infty}^{E_V}(E_V-E)^{1/2}\,e^{-(E_F-E)/k_BT}\,dE .\]
Writing \(E_F-E = (E_F-E_V)+(E_V-E)\) and putting \(x=(E_V-E)/k_BT\), the integral is the same as before, and
\[\boxed{p_0 = N_V\,e^{-(E_F-E_V)/k_BT}},\qquad \boxed{N_V = 2\left(\frac{2\pi m_h^*k_BT}{h^2}\right)^{3/2}},\]
where \(N_V\) is the effective density of states at the valence-band edge.
(c) The Fermi level in an intrinsic semiconductor
An intrinsic crystal must be electrically neutral, and since its carriers are created in pairs, \(n_0=p_0=n_i\). Denoting the intrinsic Fermi level by \(E_{Fi}\), we have
\[N_C\,e^{-(E_C-E_{Fi})/k_BT} = N_V\,e^{-(E_{Fi}-E_V)/k_BT}.\]
Taking logarithms,
\[\ln\frac{N_C}{N_V} = \frac{(E_C-E_{Fi})-(E_{Fi}-E_V)}{k_BT} = \frac{E_C+E_V-2E_{Fi}}{k_BT},\]
so that
\[E_{Fi} = \frac{E_C+E_V}{2} - \frac{k_BT}{2}\ln\frac{N_C}{N_V}.\]
Since \(N_C/N_V = (m_e^*/m_h^*)^{3/2}\), this may also be written
\[\boxed{E_{Fi} = \frac{E_C+E_V}{2} + \frac34\,k_BT\ln\frac{m_h^*}{m_e^*}}.\]
If the effective masses of electrons and holes were equal, \(N_C=N_V\) and \(E_{Fi}=(E_C+E_V)/2\): the intrinsic Fermi level would lie exactly at the centre of the forbidden gap. In practice the correction term is only a small fraction of \(E_g\) at room temperature, so the intrinsic Fermi level lies very close to mid-gap, slightly above it when \(m_h^*>m_e^*\).
(d) The intrinsic carrier concentration
Multiplying the expressions for \(n_0\) and \(p_0\), the Fermi level cancels:
\[n_0p_0 = \left(N_C\,e^{-(E_C-E_F)/k_BT}\right)\left(N_V\,e^{-(E_F-E_V)/k_BT}\right) = N_CN_V\,e^{-(E_C-E_V)/k_BT} = N_CN_V\,e^{-E_g/k_BT}.\]
The product \(n_0p_0\) is therefore a constant for a given material at a given temperature, however the material is doped. In the intrinsic case \(n_0=p_0=n_i\), so that
\[n_i^2 = N_CN_V\,e^{-E_g/k_BT},\qquad \boxed{n_i = \sqrt{N_CN_V}\;e^{-E_g/2k_BT}}.\]
Because \(N_CN_V\propto T^3\), the temperature dependence is
\[\boxed{n_i = B\,T^{3/2}\,e^{-E_g/2k_BT}},\]
where \(B\) is a constant characteristic of the material, \(E_g\) is the band gap in eV and \(k_B = 8.62\times10^{-5}\) eV/K. Equally, from the expression for \(n_0\) applied to intrinsic material,
\[n_i = p_i = N_C\,e^{-(E_C-E_{Fi})/k_BT}.\]
(e) Carrier concentrations referred to \(n_i\); the Fermi level in extrinsic material
Dividing \(n_0 = N_C e^{-(E_C-E_F)/k_BT}\) by \(n_i = N_C e^{-(E_C-E_{Fi})/k_BT}\), and similarly for holes, gives a more convenient pair of expressions:
\[\boxed{n_0 = n_i\,e^{(E_F-E_{Fi})/k_BT},\qquad p_0 = n_i\,e^{(E_{Fi}-E_F)/k_BT}}.\]
These give the position of the Fermi level in an extrinsic semiconductor at once. For n-type material
\[\boxed{E_F = E_{Fi} + k_BT\ln\frac{n}{n_i}},\]
and for p-type material
\[\boxed{E_F = E_{Fi} - k_BT\ln\frac{p}{n_i}}.\]
Donor doping raises the Fermi level above \(E_{Fi}\) towards \(E_C\), and acceptor doping lowers it towards \(E_V\), by an amount that grows only logarithmically with the doping.
Exam tip. Measured from the conduction-band edge, \(E_C-E_F = k_BT\ln(N_C/N_D)\) in n-type material. Increasing \(N_D\) by a factor \(a\) therefore moves \(E_F\) up by \(k_BT\ln a\), independently of the starting doping.
Example 2.1. The intrinsic carrier concentration of silicon (\(E_g = 1.12\) eV) at 300 K is \(1.5\times10^{10}\,\text{cm}^{-3}\). Calculate \(n_i\) for silicon at 400 K and at 500 K. (Take \(k_B = 8.62\times10^{-5}\) eV/K.)
Solution. Since \(n_i\propto T^{3/2}\exp(-E_g/2k_BT)\), the ratio of the concentrations at two temperatures is
\[\frac{n_i(T)}{n_i(300\,\text{K})} = \left(\frac{T}{300}\right)^{3/2}\exp\left[-\frac{E_g}{2k_B}\left(\frac1T-\frac1{300}\right)\right].\]
At 400 K,
\[\frac{n_i(400\,\text{K})}{n_i(300\,\text{K})} = \left(\frac{400}{300}\right)^{3/2}\exp\left[-\frac{1.12}{2\times8.62\times10^{-5}}\left(\frac1{400}-\frac1{300}\right)\right] = 1.540\times224.5 = 345.6,\]
so that
\[n_i(400\,\text{K}) = 1.5\times10^{10}\times345.6 = 5.18\times10^{12}\,\text{cm}^{-3}.\]
Similarly, at 500 K,
\[\begin{aligned}n_i(500\,\text{K}) &= 1.5\times10^{10}\times\left(\frac{500}{300}\right)^{3/2}\exp\left[-\frac{1.12}{2\times8.62\times10^{-5}}\left(\frac1{500}-\frac1{300}\right)\right]\\ &= 1.5\times10^{10}\times2.152\times5779 = 1.87\times10^{14}\,\text{cm}^{-3}.\end{aligned}\]
A rise of 200 K thus increases \(n_i\) by four orders of magnitude, almost entirely through the exponential factor.