When the molecule is initially stationary at the minimum of some state, a convenient option to generate initial geometries and velocities is to employ a Wigner PDF, assuming that the potential energy surface is harmonic around the minimum. In this case, the Wigner probability distribution for the thermal population of the vibrational states of a molecule with at N atoms and 36