Describing nuclear quantum effects in vibrational properties using molecular dynamics with Wigner sampling

In this work, we develop a molecular dynamics method, combining the simple Maxwell-Boltzmann statistics with the generally applicable standard Wigner sampling. The resulting implemented single-parameter procedure, termed simplified Wigner sampling, has been extensively benchmarked by calculating vibrationally averaged rotational constants and vibrational spectra of several molecules and by comparing with results from other methods. The procedure has been also combined with the Andersen thermostat for simulation of NVT-ensembles.