One of the currently popular sampling techniques is based on Wigner quasiprobability distributions obtained at the harmonic approximation. Let us consider the Hamiltonian Ĥ = 1 2µ p2 + µω 2 2 ξ 2 for a single vibrational mode with ξ being the displacement of the coordinate from the equilibrium position and p = -ih∂ ξ . The ground vibrational state is then given by the wavefunction