where ρ ρ = N V is the number density of molecules, h ( r ) = g ( r ) − − 1 , g ( r ) is the radial distribution function, u ( r ) is the direct interaction between pairs. β β = 1 k B T with T being the Thermodynamic temperature and k B the Boltzmann constant.