Generalized nonequilibrium Fermi's golden rule and its semiclassical approximations for electronic transitions between multiple states

The nonequilibrium Fermi's golden rule (NE-FGR) approach is developed to simulate the electronic transitions between multiple excited states in complex condensed-phase systems described by the recently proposed multi-state harmonic (MSH) model Hamiltonian. The MSH models were constructed to faithfully capture the photoinduced charge transfer dynamics in a prototypical organic photovoltaic carotenoid-porphyrin-C-60 molecular triad dissolved in tetrahydrofuran. A general expression of the fully quantum-mechanical NE-FGR rate coefficients for transitions between all pairs of states in the MSH model is obtained. Besides, the linearized semiclassical NE-FGR formula and a series of semiclassical approximations featuring Wigner and classical nuclear sampling choices and different dynamics during the quantum coherence period for the MSH model are derived. The current approach enables all the possible population transfer pathways between the excited states of the triad, in contrast to the previous applications that only addressed the donor-to-acceptor transition. Our simulations for two triad conformations serve as a demonstration for benchmarking different NE-FGR approximations and show that the difference between all levels of approximation is small for the current system, especially at room temperature. By comparing with nonadiabatic semiclassical dynamics, we observe similar timescales for the electronic population transfer predicted by NE-FGR. It is believed that the general formulation of NE-FGR for the MSH Hamiltonian enables a variety of applications in realistic systems.