Infrared absorption line shapes in the classical limit: a comparison of the classical dipole and fluctuating...

Infrared spectroscopy is a versatile technique for probing the structure and dynamics of condensed-phase systems. Simulating infrared absorption spectra with molecular dynamics (MD) offers a powerful means to establish a molecular-level interpretation of experimental results, as well as a basis for the parametrization of more accurate simulation force-fields. Two distinct methods for the calculation of infrared absorption line shapes of high-frequency (Planck's omega/k(B)T>>1) vibrational probes from MD simulations are examined: The classical dipole approximation (CDA) and the fluctuating frequency approximation (FFA). Although these two formalisms result in expressions for the infrared absorption line shape that appear very different, both approximations are shown to yield identical results for the infrared line shape of a harmonic system in the condensed-phase. The equivalence of the FFA and CDA is also demonstrated in the case where the transition dipole of the oscillator fluctuates in response to the environment (i.e., where the Condon approximation has been relaxed). Finally we examine the effects of solute anharmonicity and demonstrate that the CDA and FFA are not equivalent in general, and the magnitude of the deviations increases with anharmonicity. We conclude that the calculation of infrared absorption line shapes via the CDA is a promising alternative to the FFA approach in cases where it may be difficult or undesirable to employ the latter, particularly when the effects of anharmonicity are small.