Ab initio molecular dynamics simulation of the aqueous Ru2+/Ru3+ redox reaction: the Marcus perspective

A well-behaved (low spin) transition metal aqua ion, Ru(aq)(2+), is used as a model system in an ab initio molecular dynamics study of a redox half reaction to which the Marcus theory of electron transfer is assumed to apply. Using constraint methods, we show that aqueous Ru(2+) can be reversibly transformed to Ru(3+) under the control of the classical solvent electrostatic potential as order parameter. The mean force is found to vary linearly with the order parameter in accordance with the Marcus theory. As can be expected for a half reaction, the slope in the oxidized and reduced states are asymmetric differing by approximately a factor of two. As a further test, we verify that the corresponding quadratic potential of mean force is in excellent agreement with the free energy profile obtained from the Gaussian distribution of potential fluctuations sampled from free (unconstrained) runs of the reduced and oxidized system. Similar to experimental electrochemical methods, our simulation scheme enables us to manipulate the effective thermodynamic driving force and align the free energy minima of product and reactant state. The activation energy and reaction entropy computed under these conditions are discussed and analyzed from the Marcus perspective.