Modeling of the excited states of multichromophoric systems is crucial for the understanding of photosynthesis functioning. The excitonic Hamiltonian method is widely used for such calculations. Excited states of the combined system are constructed from the wave functions of individual chromophores while their interactions are described by excitonic couplings. In the current study we enhance earlier proposed non-orthogonal product approach to incorporate dynamic correlation effects accounted for by the multireference perturbation theory. We discuss the problems of constructing the excitonic Hamiltonian including charge-transfer states for the molecular systems where the overlap contribution to the excitonic couplings is non-negligible. The benchmark calculations were performed for a model system. It was shown that the overlap component of the excitonic coupling is of great importance. The enhanced method provides an accurate description of the excited states energies and other properties.
The new version has been updated according to reviewers' suggestions. More points have been added on the graphs, additional explanations have been added in several sections.