Classical Trajectory of Molecules in Electromagnetic Field: A Handy Method to Simulate Molecular Vibrational Spectra

Within harmonic approximations, molecular vibrational spectra are simulated in a standard way through force field diagonalization and following transformation of Cartesian to normal-mode tensor derivatives. This may become tedious for large systems of many thousands of atoms and also not necessary because of a limited resolution required to interpret an experiment. We developed an algorithm based on the real-time real-field molecular dynamics, effectively at zero temperature, invoked in a molecule by the electromagnetic field of light. The algorithm is simple to implement and suitable for parallel computing, and it can be potentially extended to more complicated molecular-light interaction modes. It circumvents the diagonalization and is suitable to model vibrational optical activity (vibrational circular dichroism and, to a lesser extent, Raman optical activity). For large molecules, it becomes faster than diagonalization, but it also enables the assignment of vibrational spectral bands to local molecular motions.