Abstract
A third order expansion for Wigner-Kirkwood commutation function, a complex function in classical phase space that accounts for the Heisenberg uncertainty relation, is approximated and integrated over momentum to give a real function in position configuration space. Metropolis Monte Carlo computer simulation results are given for liquid Lennard-Jones ⁴He below 10\,K.