The new numerical version of the Wigner approach to quantum mechanics for treatment thermodynamic properties of the strongly interacting systems of particles has been developed for extreme conditions, when there are no small physical parameters and analytical approximations used in different kind of perturbation theories can not be applied. The new path integral representation of the quantum Wigner function in the phase space has been developed for canonical ensemble. Explicit analytical expression of the Wigner function has been obtained in linear and harmonic approximations. The new quantum Monte-Carlo method for calculations of average values of arbitrary quantum operators has been proposed. Preliminary calculations of the momentum distribution function of the Coulomb systems of particles have been carried out. Comparison with classical Maxwell-Boltzmann distribution shows the significant influence of quantum effects on the high energy asymptotics ("tails') of the calculated momentum distribution functions, which resulted in appearance of sharp oscillations.
In contrast to classical physics, the language of quantum mechanics involves operators and wave functions (or, more generally, density operators). However, in 1932, Wigner formulated quantum mechanics in terms of a distribution function
W(q,p), the marginals of which yield the correct quantum probabilities for
q and
p separately \cite{wigner}. Its usefulness stems from the fact that it provides a re-expression of quantum mechanics in terms of classical concepts so that quantum mechanical expectation values are now expressed as averages over phase-space distribution functions. In other words, statistical information is transferred from the density operator to a quasi-classical (distribution) function.