me = 9.109*10^-31; me2 = 0.01*me; k = 1.38*10^-23; 
    em = 1.602*10^-19; Ed = 1*10^-3*em; h = 6.626*10^-34; 
    NC = 2*(2*π*me2*k*4/(h^2))^(3/2); 
    ND = 10^13*10^6;

NC

4*ND/NC*Exp[Ed/k/4]

nd = Sqrt[ND*NC]*Exp[-Ed/k/4/2]

Plot[4*ND/NC*Exp[Ed/k/T], {T, 10^5, 2.7169*10^6}, 
      ImageSize -> 550]

4*ND/NC

Solve[4*ND/NC*Exp[Ed/k/T] == 1.0361, T]