Donor charge-state degeneracy ($n_{\mathrm{full}}/n_{\mathrm{empty}}=r{e}^{\beta (\mu -{ϵ}_D)}$)

Let an isolated donor impurity admit an empty ionized state with statistical weight $g_0$ and a singly occupied neutral state with statistical weight $g_1$, with double occupation excluded. Its grand canonical ensemble weights are $g_0$ and $g_1{e}^{-\beta ({ϵ}_D-\mu )}$, so $r=g_1/g_0$ multiplies the occupied-to-empty ratio. A nondegenerate empty state and two allowed electron spin orientations give $r=2$; a simplified single-state occupation model uses $r=1$. With negligible electron holes and uncompensated donor impurities, charge neutrality and nondegenerate conduction occupation imply $n^2/(N_D-n)=(N_c/r)e^{-E_i/(k_{B}T)}$, where $E_i$ is the donor binding energy and $N_c$ the effective density of states.