Gas phase infrared spectra via the phase integration quasi-classical method

We apply the recently developed phase integration method (PIM) (Monteferrante et al. Mol Phys. 2011;109:3015-3027) to the calculation of infrared spectra of gas phase molecules. The PIM combines a generalised Monte Carlo sampling of the exact thermal density of the system with classical molecular dynamics to obtain approximate time quantum correlation functions. To describe the molecules, we adopt very simple analytical potentials that have, however, proved interesting, and surprisingly challenging, benchmarks for approximate quantum dynamical schemes. We show that, in contrast with two other commonly applied methods, our spectra do not exhibit spurious features or unphysical shifts depending on the temperature. Identifying the positions of the peaks requires only a few tens of trajectories, while an accurate evaluation of the relative intensities of the peaks is computationally more demanding.