We know that high doping densities cause the band gap to shrink in semiconductor. I found the general relationship between doping density and band gap change as
delta E=(-3*q^2/16*pi*epsilon)*square root(q^2*N/epsilon*k*T)---(i)
where delta E is the change of band gap energy,N is the doping density, q is the electronic charge, epsilon is the dielectric constant of the semiconductor, k is Boltzmann's constant,T is the temperature in Kelvin and * is "multiplication sign" , / is "division sign" and ^ is "to the power" sign. For silicon,epsilon = 11.9 and T=300K.
Now how can I reduce the equation (i) to
delta E= -22.5*square root (N/10^18 cm^-3) meV
for silicon.