An Efficient Numerical Algorithm for Constructing the Wigner Function of a Quantum System with a Polynomial Potential in Phase Space

When considering quantum systems in phase space, the Wigner function is used as a function of quasi-probability density. Finding the Wigner function is related to the calculation of a Fourier transform of a certain composition of wave functions of the corresponding quantum system. As a rule, knowledge of the Wigner function is not the final goal, and the calculation of mean values of different quantum characteristics of the system is required. An explicit solution of the Schrodinger equation can be obtained only for a narrow class of potentials, therefore, in the majority of cases, numerical methods must be used to find wave functions. As a result, finding the Wigner function is related to the numerical integration of grid wave functions. When considering a one-dimensional system, the calculation of N2 Fourier integrals of the grid wave function is needed. To provide the necessary accuracy for wave functions, corresponding to the higher states of a quantum system, a larger number of grid nodes is needed. The goal of this work was to construct the numerical-analytical method for finding the Wigner function, which allows the number of computational operations to be considerably reduced. The quantum systems with polynomial potentials were considered, for which the Wigner function is represented as a series in some known functions.