8.4 Marcus theory of electron transfer

Fundamentals of Marcus Theory

Marcus theory provides a quantitative framework for predicting the rates of electron transfer reactions. Developed by Rudolph Marcus in the 1950s (earning him the 1992 Nobel Prize in Chemistry), it connects thermodynamics and kinetics by showing how both the energy released in a reaction and the structural rearrangements required to accommodate that transfer jointly determine how fast the transfer occurs. This framework applies broadly, from simple inorganic redox couples to the charge-separation steps in photosynthesis.

Key Concepts

Reorganization energy ($\lambda$) is the energy the system would need to rearrange all nuclear coordinates from the reactant geometry to the product geometry without actually transferring the electron. It has two components:

• Inner-sphere ($\lambda_i$): Changes in bond lengths and angles within the donor and acceptor molecules themselves. For example, when $\text{Fe}^{2+}$ loses an electron to become $\text{Fe}^{3+}$, the metal-ligand bonds shorten.
• Outer-sphere ($\lambda_o$): Reorientation of surrounding solvent molecules in response to the new charge distribution. In polar solvents like water, this component is often dominant.

Driving force ($\Delta G°$) is the standard Gibbs free energy change for the electron transfer. A negative $\Delta G°$ means the reaction is exergonic (thermodynamically favorable). The more negative $\Delta G°$ is, the more energy is released.

Activation energy ($\Delta G^\ddagger$) is the energy barrier the system must overcome for the transfer to occur. Marcus theory's central achievement is providing an exact relationship between $\Delta G^\ddagger$, $\lambda$, and $\Delta G°$.

Potential Energy Surfaces

Marcus theory models the reactant and product states as two parabolic potential energy surfaces plotted against a generalized nuclear coordinate (sometimes called the reaction coordinate). This coordinate captures all the structural changes: bond lengths, bond angles, and solvent positions.

The parabolas have equal curvature (the harmonic oscillator approximation), and they're offset both horizontally (different equilibrium geometries) and vertically (by $\Delta G°$). The point where the two parabolas intersect is the transition state. At this geometry, the electron can transfer because the reactant and product states have the same energy, satisfying the Franck-Condon principle (nuclei don't move during the instantaneous electron hop).

The height of the intersection point above the reactant minimum gives $\Delta G^\ddagger$. Changing either $\lambda$ or $\Delta G°$ shifts the parabolas relative to each other, raising or lowering this intersection.

Mathematical Framework and Applications

The Marcus Equation

The central result of the theory is:

$\Delta G^\ddagger = \frac{(\lambda + \Delta G°)^2}{4\lambda}$

This expression for the activation barrier feeds directly into the rate constant through transition-state-like kinetics:

$k_{ET} = A \exp\!\left[-\frac{(\lambda + \Delta G°)^2}{4\lambda k_B T}\right]$

where:

• $k_{ET}$ is the electron transfer rate constant
• $A$ is the pre-exponential factor, determined primarily by the electronic coupling $H_{DA}$ between donor and acceptor (stronger orbital overlap means larger $A$)
• $k_B$ is the Boltzmann constant
• $T$ is the absolute temperature

How to read the equation: The numerator $(\lambda + \Delta G°)^2$ controls the barrier height. Because $\Delta G°$ is negative for exergonic reactions, increasing the magnitude of $\Delta G°$ first decreases the numerator (lowering the barrier) and then, once $|\Delta G°|$ exceeds $\lambda$, increases it again. This produces three distinct regimes:

1. Normal region ($|\Delta G°| < \lambda$): Making the reaction more exergonic lowers the barrier and speeds up the rate. This matches everyday chemical intuition.
2. Activationless point ($-\Delta G° = \lambda$): The barrier vanishes ($\Delta G^\ddagger = 0$), and the rate reaches its maximum. The two parabolas intersect exactly at the reactant minimum.
3. Inverted region ($|\Delta G°| > \lambda$): Making the reaction even more exergonic now raises the barrier and slows down the rate. This counterintuitive prediction was the most controversial aspect of Marcus theory and took decades to confirm experimentally.

Applications

Chemical systems:

• Self-exchange reactions like $\text{Fe}^{2+}/\text{Fe}^{3+}$ are the simplest test cases because $\Delta G° = 0$, so the barrier depends only on $\lambda$.
• The Creutz-Taube ion (a mixed-valence ruthenium dimer) provided early evidence for strong electronic coupling and delocalized electron transfer.

Biological systems:

• In photosynthetic reaction centers, the initial charge separation after chlorophyll excitation is extremely fast (picoseconds) because $-\Delta G°$ is tuned close to $\lambda$, placing the system near the activationless point. Subsequent back-electron-transfer (charge recombination) is much slower because it falls in the inverted region, which is exactly what the organism needs to preserve the separated charges.
• In the mitochondrial electron transport chain, each step involves a modest driving force and relatively small reorganization energies, keeping transfer rates fast and efficient.

Key factors that modulate $k_{ET}$:

• Donor-acceptor distance: The electronic coupling $H_{DA}$ decays exponentially with distance, roughly as $H_{DA} \propto e^{-\beta r}$, where $\beta$ depends on the intervening medium (about 1.0–1.4 Å$^{-1}$ through protein, higher through vacuum).
• Solvent polarity: Higher dielectric constant solvents increase $\lambda_o$, which shifts the activationless point to larger $|\Delta G°|$ values and changes which regime the reaction falls in.
• Temperature: Higher $T$ increases the rate in the normal region (classical barrier crossing) but the effect is more complex in the inverted region, where nuclear tunneling can become significant.

Experimental techniques:

• Transient absorption spectroscopy tracks the appearance and decay of radical ion intermediates on femtosecond-to-nanosecond timescales, directly measuring $k_{ET}$.
• Electrochemical methods (cyclic voltammetry, chronoamperometry) probe redox potentials and heterogeneous electron transfer rates at electrode surfaces, providing $\Delta G°$ and rate data in a single experiment.