Computing Anharmonic Infrared Spectra of Polycyclic Aromatic Hydrocarbons Using Machine-Learning Molecular Dynamics

We introduce a machine learning molecular dynamics (MLMD) approach to calculate the anharmonic infrared (IR) absorption spectra of polycyclic aromatic hydrocarbons (PAHs), key carriers of interstellar aromatic IR bands. This method accounts for temperature effects in a molecule-specific way and achieves accuracy comparable to conventional quantum chemical calculations at a fraction of the cost, scaling linearly with system size. We applied MLMD to calculate the anharmonic spectra of 1,704 PAHs in the NASA Ames PAH IR Spectroscopic Database with up to 216 carbon atoms at different temperatures, demonstrating its capability for high-throughput spectral calculations of large molecular systems.