Vibronic transitions in large molecular systems: rigorous prescreening conditions for Franck-Condon factors

In this work, rigorous prescreening conditions for the calculation of Franck-Condon factors (FCFs) are derived and implemented. These factors play an important role in numerous applications including the prediction and simulation of vibronic spectra, electron transfer rates, and nonradiative transition probabilities. For larger systems it is crucial to calculate only relevant FCFs, as the computational burden becomes otherwise prohibitive due to the sheer number of Franck-Condon integrals. By exploiting rigorous prescreening criteria one can significantly reduce the computational effort and systematically refine results to the desired target accuracy. In this work, such criteria are derived via the use of sum rules obtained through a coherent state generating function for the FCFs in the harmonic oscillator approximation, following the prescription of Doktorov et al. These sum rules allow efficient and rigorous prescreening prior to the calculation of entire batches of Franck-Condon integrals, reducing the subsequent computational burden. To illustrate the benefit of employing such conditions, they are applied in this work to calculate FC profiles for vibronic spectra of formic acid, thymine, anthracene, and a polycyclic aromatic hydrocarbon derivative containing about 450 vibrational degrees of freedom. Since the prescreening step gives stringent upper and lower bounds for the loss of FC intensity in these spectra arising from the neglected FCFs, the present approach renders most previous a priori selection schemes obsolete and has the potential to complement or even replace other approximate treatments.