However, the extra valance electron moves in the semiconductor's conduction band and not free space.
Therefore, the electron's mass is the conduction band's effective mass.
Furthermore, the interactions between the electron and proton are screened by the lattice.
As a result, we need to introduce the following substitutions: $m_e \to m_e^*$, $\epsilon_0 \to \epsilon\epsilon_0$.
We thus estimate the energy of the bound state created by the impurity:
$$E = -\frac{m_e^*}{m_e \varepsilon^2} R_E = -0.01 \text{eV (in Ge)}$$
$r_B = 4$ nm (vs $r_B = 0.5$ Å in H).
The electron is very weakly bound to the impurity! At room temperature (0.026 eV), the donor electron is easily thermally excited into the conduction band.
On the other hand, we can add a group III element to reduce the average number of electrons in the system.
Group III elements lacks 1 electron and 1 proton and are therefore known as **acceptors**.
We treat the absence of an electron as a hole and the lacking proton as an effective negative charge.
As a result, we once again end up with a Hydrogen model, except this time the charges are flipped (hole circles around a negative center).
That allows us to use the previous results and to conclude that an acceptor creates a weakly bound state above the valance band.
How large can $N_D/N_A$ be? The distance between donors should be such that the states don't overlap, so the distance must be much larger than 4 nm. Therefore **maximal** concentration of donors before donor band merges with conduction band is $N_D \lesssim (1/4\textrm{nm})^3 \sim 10^{-5}\ll N_C$.
In order to model **multiple** donor/acceptor states, we assume that they are all degenerate at the binding energy.
Therefore, we model the density of states of donors/acceptors as a Dirac delta function: