Finite temperature infrared spectroscopy of polycyclic aromatic hydrocarbon molecules: Path-integral molecular dynamics

Abstract

The vibrational spectra of the naphthalene, pyrene, and coronene molecules have been computed in the 0-3500 cm(-1) infrared range using classical and quantum molecular dynamics simulations based on a dedicated tight-binding potential energy surface. The ring-polymer molecular dynamics (RPMD) and partially adiabatic centroid molecular dynamics (CMD) methods have been employed to account for quantum nuclear effects. The contributions of quantum delocalization to the line shift and broadening are significant in the entire spectral range and of comparable magnitude as pure thermal effects. While the two methods generally produce similar results, the CMD method may converge slower at low temperature with increasing Trotter discretization number. However, and contrary to the CMD method, the RPMD approach suffers from serious resonance problems at high frequencies and low temperatures.