Improving anharmonic infrared spectra using semiclassically prepared molecular dynamics simulations

Classical molecular dynamics is a convenient method for computing anharmonic infrared spectra of polyatomic molecules and condensed phase systems. However it does not perform well for predicting accurate intensities and it lacks nuclear quantization, two deficiencies that are usually accounted for by empirical scaling factors. In this paper we show on the examples of the trans isomer of nitrous acid and naphthalene that both issues can be alleviated by preparing the initial conditions according to semiclassical quantization based on a normal mode representation. The method correctly reproduces fundamental frequencies obtained with quantum mechanical methods. At increasing temperatures, the effective frequencies are found to follow the same trends as path-integral based methods. In the low-temperature limit, the band intensities predicted by the method are also found to agree with quantum mechanical considerations.