An efficient approach for the calculation of Franck-Condon integrals of large molecules.

A general and efficient approach for the calculation of Franck-Condon integrals (FCIs) of large molecules is presented. In a first step, by exploiting the diagonally dominant and sparse structure of the Duschinsky matrix, a model system is constructed for which the Duschinsky matrix takes a block-diagonal form. For each of these blocks separately, the FCIs are calculated discarding all below a certain threshold. From those integrals retained the FCIs of the model system are obtained by simple multiplication. These serve as an estimate for the FCIs of the exact system which are calculated for those integrals which lie above a certain threshold. By systematically decreasing the threshold, the simulation can be reliably converged to the exact result with an arbitrary accuracy. Using this scheme, a considerable reduction of the number of FCIs which have to be calculated is achieved which leads to an improved scaling behavior of the computational effort with system size. The approach has been tested thoroughly for a set of molecules including difficult cases. For the larger systems a speedup of up to three orders of magnitude compared to an exact calculation is observed while the errors can be kept negligible. With this approach accurate calculations of FCIs are feasible also for large molecules encountered in "real-life" chemistry, especially biochemistry and material science.