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)]
Quantity[8.10637*10^8, (
Sqrt["BoltzmannConstant"] Sqrt["Farads"] Sqrt["Kelvins"])/(
"Coulombs" ("Meters")^2)]
Quantity[
7.86337*10^-38, (("Coulombs")^3 ("Meters")^3)/((
"BoltzmannConstant")^(3/2) ("Farads")^(3/2) ("Kelvins")^(3/2))]