Brownian simulation of electron transfer dynamics

Several algorithms for the numerical simulation of electron transfer reactions in condensed media within the framework of the generalized Zusman model are considered. The model accounts for the population of a number of electronic and vibrational states of the donor-acceptor pair as well as the complex dynamics of solvent relaxation. The proposed algorithms are based on the calculations of the particle's survival function along the Brownian trajectory in configuration space. The method is shown to have certain advantages compared to the conventional grid methods and is expected to be the most efficient one for multi-mode/multi-level problems. The validity conditions for the method are studied and some test computations for the reaction dynamics are performed in a number of special cases. The results obtained are compared with the corresponding analytic estimates. The work was supported by the Russian Foundation for Basiс Research (project N 08-03-00534).