В полупроводниках при температурах, отличных от абсолютного нуля, часть электронов из валентной зоны возбуждается в зону проводимости. При этом в валентной зоне остаются незаполненные состояния, которые рассматриваются как квазичастицы — дырки. Дырка обладает положительным зарядом, а её поведение в электрических и магнитных полях аналогично движению положительно заряженного носителя. Для корректного описания процессов в полупроводниках необходимо рассматривать статистику распределения дырок по энергетическим состояниям.
Электроны в полупроводнике подчиняются распределению Ферми–Дирака:
$$ f(E) = \frac{1}{e^{(E - E_F)/kT} + 1}, $$
где
Вероятность того, что электронное состояние в валентной зоне окажется свободным, определяется как
1 − f(E).
Эта величина и есть вероятность существования дырки в данном состоянии. Таким образом, распределение дырок можно записать как
$$ f_p(E) = 1 - f(E) = \frac{1}{e^{-(E - E_F)/kT} + 1}. $$
Заметим, что эта функция описывает вероятность заполненности энергетических уровней дырками, и по форме является функцией Ферми–Дирака, но с противоположным знаком в экспоненте.
Валентная зона характеризуется своей плотностью состояний. Для электронов в валентной зоне плотность состояний вблизи её верхнего края Ev можно записать как
$$ g_v(E) = \frac{1}{2\pi^2}\left(\frac{2m_p^*}{\hbar^2}\right)^{3/2} \sqrt{E_v - E}, $$
где
Так как дырка — это отсутствие электрона, то энергетический спектр дырок удобно отсчитывать от уровня Ev, где энергия дырки Ep = Ev − E.
Численная концентрация дырок определяется интегрированием произведения плотности состояний на вероятность их занятости дырками:
p = ∫−∞Evgv(E) fp(E) dE.
Для удобства интеграл переписывают через энергию дырки:
p = ∫0∞gv(Ep) fp(Ep) dEp,
где плотность состояний для дырок имеет вид
$$ g_v(E_p) = \frac{1}{2\pi^2}\left(\frac{2m_p^*}{\hbar^2}\right)^{3/2} \sqrt{E_p}. $$
В большинстве практически используемых полупроводников уровень Ферми лежит вблизи середины запрещённой зоны или ближе к зоне проводимости. Тогда условие
EF ≪ Ev − 3kT
выполняется, и можно использовать максвелловское приближение:
fp(E) ≈ e(E − EF)/kT.
Подставляя это в выражение для концентрации, получаем:
p = Nv e−(EF − Ev)/kT,
где
$$ N_v = 2\left(\frac{2\pi m_p^* kT}{h^2}\right)^{3/2} $$
— эффективная плотность состояний в валентной зоне.
p = n = ni,
где ni — собственная концентрация носителей.
p ≈ NA,
а уровень Ферми смещается ближе к валентной зоне.