Thanks for your reply. I’m still getting some weird results. No sure how to set units to get a numerical value. Can you provide me any clue? Thanks once again
T = Quantity[190, "Kelvins"];
kb = Quantity[1, "BoltzmannConstant"];
n = Quantity[10^11, "Meters"^3];
\[Epsilon] = Quantity[8.854 *10^-12, "Farads"/"Meters"];
e = Quantity[1.60 *10^-19, "Coulombs"];
N[LD = ((kb T \[Epsilon])/(e^2 n))^(1/2)]
N[ND = 4 \[Pi]/(3 n (LD)^3)]