Theoretical Infrared Spectra: Quantitative Similarity Measures and Force Fields

Infrared spectroscopy can provide significant insight into the structures and dynamics of molecules of all sizes. The information that is contained in the spectrum is, however, often not easily extracted without the aid of theoretical calculations or simulations. We present here the calculation of the infrared spectra of a database of 703 gas phase compounds with four different force fields (CGenFF, GAFF-BCC, GAFF-ESP, and OPLS) using normal-mode analysis. Modern force fields increasingly use virtual sites to describe, e.g., lone-pair electrons or the o -holes on halogen atoms. This requires some adaptation of code to perform normal-mode analysis of such compounds, the implementation of which into the GROMACS software is briefly described as well. For the quantitative comparison of the obtained spectra with experimental reference data, we discuss the application of two different statistical correlation coefficients, Pearson and Spearman. The advantages and drawbacks of the different methods of comparison are discussed, and we find that both methods of comparison give the same overall picture, showing that present force field methods cannot match the performance of quantum chemical methods for the calculation of infrared spectra.