Approaching black-box calculations of pump-probe fragmentation dynamics of polyatomic molecules

A general framework for the simulation of ultrafast pump-probe time resolved experiments based on Born-Oppenheimer molecular dynamics (BOMD) is presented. Interaction of the molecular species with a laser is treated by a simple maximum entropy distribution of the excited state occupancies. The latter decay of the electronic excitation into the vibrations is based on an on-the-fly estimation of the rate of the internal conversion, while the energy is distributed in a thermostat-like fashion. The approach was tested by reproducing the results of previous femtosecond studies on ethylene, naphthalene and new results for phenanthrene