Wigner sampling is one of possibilities for producing initial conditions in MD simulations, that emulate the quantum harmonic ground vibrational state. In this work we combine Wigner sampling with MD simulations to provide anharmonic vibrational properties of molecules. We also propose a new simplified sampling procedure, which is based on Wigner sampling but does not require preliminary calculation of molecular Hessian.