Molecular Dynamics Simulations of Bimolecular Electron Transfer: the Distance-Dependent Electronic Coupling

Understanding the distance dependence of the parameters underpinning Marcus theory is imperative when interpreting the results of experiments on electron transfer (ET). Unfortunately, most of these parameters are difficult or impossible to access directly with experiments, necessitating the use of computer simulations to model them. In this work, we use molecular dynamics simulations in conjunction with constrained density functional theory calculations to study the distance dependence of the electronic coupling matrix element, |H_RP|, for bimolecular ET. Contrary to what is typically assumed for such intermolecular reactions, we find that the magnitude of |H_RP| does not decay exponentially with the center-of-mass separation of the reactants, r_COM. The addition of other simple measures of donor/acceptor (D/A) orientation did not improve the correlation of |H_RP| with r_COM. Using the minimum distance separation, r_min, of the reactants as the structural descriptor allowed the system to be partitioned into high-coupling/close-contact and low-coupling/non-contact regimes, but large fluctuations of |H_RP| were still found for the close-contact reactant pairs. Despite the persistent large fluctuations of |H_RP|, its mean value was found to decay piecewise exponentially with increasing r_min, which was attributed to significant changes in the average D/A pair structure. The results herein advise one to use caution when interpreting the experimental results derived from spherical reactant models of bimolecular ET.