Microscopic based density matrix treatments of electron-transfer reactionsin condensed phases

Several non-phenomenological density matrix treatments of electron-transfer (ET) reactions in condensed phases are developed and examined. The methods consider the donor and acceptor system (the solute) under the influence of the surrounding fluctuating solvent. The main emphasis is placed on semicl assical methods, where the starting point is the Hamiltonian of the quantum mechanical electronic states of the solute. The diagonal elements of the H amiltonian include the fluctuations of the solute electronic energies as a result of the interaction between the solute and the field from the classic ally moving solvent molecules. The fluctuating Hamiltonian is used to const ruct a Liouville equation, which is treated by three approaches. The first method is based on a direct numerical integration of the relevant Liouville equation. The second involves the use of a second-order Liouville equation , and the third involves the use of a Redfield type equation. The methods a re examined by simulating electron transfer between two sodium-like atoms t hat are held at a 4 Angstrom separation in water. The simulations generate the fluctuations of the electronic energies of the states that are involved in the electron-transfer process. The fluctuating energies are then used i n evaluating the rate constant of the reaction as a function of its assumed free energies. The results of the three approaches are similar to the corr esponding results obtained from the Marcus equation. However, the Redfield equation converges much more quickly than the direct Liouville equation and its second-order version. The problems associated with the semiclassical t reatments are briefly considered, emphasizing the approximation involved in treating the solvent motion classically. Some of these problems can be ove rcome by a previously developed density matrix approach(1) that uses classi cal simulations to evaluate the Franck-Condon factors of the solvent vibron ic states. This vibronic density matrix treatment is briefly described and used in simulating an electron-transfer reaction in the reaction center fro m Rps. viridis.