Effective method to compute Franck-Condon integrals for optical spectra of large molecules in solution.

The authors present a new method for the computation of vibrationally resolved optical spectra of large molecules, including the Duschinsky [Acta Physicochim. URSS 7, 551 (1937)] rotation of the normal modes. The method automatically selects the relevant vibronic contributions to the spectrum, independent of their frequency, and it is able to provide fully converged spectra with a quite modest computational time, both in vacuo and in condensed phase. Starting from the rigorous time-dependent expression they discuss indeed in which limits the spectrum of a molecule embedded in a solvent, described as a polarizable continuum, can be computed in a time-independent formalism, defining both nonequilibrium and equilibrium limits. In these cases the polarizable continuum model provides a suitable description of the solvent field. By computing the absorption spectra of anthracene in gas phase and of coumarin C153 in gas phase and cyclohexane, and the phosphorescence spectrum of the unsubstituted coumarin in ethanol they show that the method is fast and efficient.