The authors extend their recent method for the computation of vibrationally resolved optical spectra of large molecules, including both the Duschinsky rotation and the effect of finite temperature in the framework of the Franck–Condon (FC) approximation, to deal with the more general case of the Herzberg–Teller (HT) model, where also the linear dependence of the transition dipole moment on the nuclear coordinates is taken into account. This generalization allows us to investigate weak and vibronically allowed transitions by far extending the range of application of the method. The calculation of the spectra of sizable molecules is computationally demanding because of the huge number of final vibrational states that must be taken into account, and the inclusion of HT terms further increases the computational burden. The method presented here automatically selects the relevant vibronic contributions to the spectrum, independent of their frequency, and it is able to provide fully converged spectra with a modest computational requirement. The effectiveness of the method is illustrated by computing the HT absorption and fluorescence Qx spectra of free-base porphyrin both at T=0K and at room temperature, performing for the first time an exact treatment of vibrations in harmonic approximation. Qx spectra are compared to experiments and FC/HT interferences are analyzed in detail.

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.