Indeed it is better to have the factor
1/(4 * pi * eps_0)
for the reduced dipole moment
eps_0 being the vaccum permittivity, thus
1/(4 * pi * eps_0) = 1e-7 * c^2
with c the celerity of light
reprivate A and B method in thermal cond,
delete non defined function name in Stockmayer potential