Marcus Theory of Electron Transfer — Including the Inverted Region
Almost every kinetics rule a student learns says the same thing: make a reaction more exergonic and it goes faster. Marcus theory of electron transfer says that is true only up to a point, and that past that point extra driving force makes the reaction slower. That prediction — the inverted region — was so counterintuitive when first made that it took years of experimental work before it was accepted, and the theory earned Rudolph A. Marcus the Nobel Prize in Chemistry in 1992. It is standard CSIR-NET and GATE physical–inorganic material, and the core equation is small enough to work with by hand.
Why an electron transfer needs an activation barrier at all
An electron is light and moves far faster than nuclei. So at the instant the electron jumps, the atomic positions and the surrounding solvent molecules cannot move — the Franck–Condon principle applies here just as it does in spectroscopy. That creates a problem. Before transfer, the solvent is organised around the charge distribution of the reactants; after transfer, it needs to be organised around the products. If the electron simply jumped, the system would land in a badly solvated, high-energy arrangement, and energy would not be conserved.
Nature's solution is that thermal fluctuations must first distort the reactants and their solvent shell into a geometry where the reactant and product energy surfaces have the same energy. Only then can the electron move without any energy change. Reaching that crossing point costs energy, and that cost is the activation barrier.
The Marcus equation
- ΔG‡ — the free energy of activation for the electron transfer step.
- ΔG° — the driving force, i.e. the standard free energy change of the transfer. It is negative for an exergonic reaction, and getting this sign right is the whole game.
- λ (lambda) — the reorganisation energy: the energy it would cost to distort the reactants and their solvent shell into the product geometry without transferring the electron. Always positive.
λ splits into two physically distinct pieces:
- λinner — bond lengths and angles inside the reacting species. A metal ion whose M–L bonds change a lot on reduction has a large inner-sphere term. This is why [Co(NH3)6]3+/2+ self-exchange is slow: the change between low-spin d6 and high-spin d7 alters the bond lengths substantially.
- λouter — reorientation of the solvent. The standard dielectric-continuum estimate is λout = (Δe)²[1/(2a1) + 1/(2a2) − 1/R] [1/εop − 1/εs], where a1 and a2 are the radii of the two reactants, R their separation, and εop and εs the optical and static dielectric constants of the solvent. Because εop < εs, the bracket is positive, and polar solvents give larger λout.
Once you have ΔG‡, the rate follows the usual exponential form:
In the fuller semiclassical treatment the prefactor A contains the electronic coupling |HAB|² between donor and acceptor, which falls off roughly exponentially with distance. That is why long-range biological electron transfer chains place their redox centres close together.
Worked example 1 — a normal-region calculation
Problem: a self-assembled donor–acceptor pair has λ = 1.00 eV and ΔG° = −0.40 eV. Find ΔG‡ in kJ mol−1 and the Boltzmann factor at 298 K.
Step 1 — the barrier.
λ + ΔG° = 1.00 + (−0.40) = 0.60 eV
(0.60)² = 0.36
4λ = 4 × 1.00 = 4.00
ΔG‡ = 0.36 ÷ 4.00 = 0.090 eV
Step 2 — convert to kJ mol−1 using 1 eV = 96.485 kJ mol−1:
0.090 × 96.485 = 8.68 kJ mol−1
Step 3 — the exponential factor.
RT = 8.314 × 298 = 2478 J mol−1 = 2.478 kJ mol−1
ΔG‡/RT = 8.68 ÷ 2.478 = 3.503
exp(−3.503) = 0.0301
Step 4 — a rate estimate. Taking a typical prefactor of
A = 1 × 1013 s−1 (an assumed value, chosen to show the size of the
result rather than to describe any particular system):
k = 1 × 1013 × 0.0301 = 3.0 × 1011 s−1
— a sub-picosecond process.
Worked example 2 — the three regimes side by side
Keep λ fixed at 1.00 eV and change only the driving force.
(a) ΔG° = −0.40 eV (normal region, −ΔG° < λ)
ΔG‡ = (1.00 − 0.40)² ÷ 4.00 = 0.36 ÷ 4.00 = 0.090 eV
(b) ΔG° = −1.00 eV (activationless, −ΔG° = λ)
ΔG‡ = (1.00 − 1.00)² ÷ 4.00 = 0 ÷ 4.00 = 0 eV
The barrier vanishes completely and the rate is at its maximum. This is the peak of the
Marcus curve.
(c) ΔG° = −1.80 eV (inverted region, −ΔG° > λ)
λ + ΔG° = 1.00 − 1.80 = −0.80; squaring removes the sign: (−0.80)² = 0.64
ΔG‡ = 0.64 ÷ 4.00 = 0.16 eV = 15.44 kJ mol−1
Comparing (a) and (c). Case (c) has 1.40 eV more driving force than case
(a), yet its barrier is higher. The rate ratio at 298 K:
ΔΔG‡ = 15.44 − 8.68 = 6.76 kJ mol−1
k(c)/k(a) = exp(−6.76 ÷ 2.478) = exp(−2.728) = 0.0654
So case (c) is 1 ÷ 0.0654 = about 15 times slower despite being far more
exergonic.
The squared term is doing all the work. Because the numerator is (λ + ΔG°)², the barrier is a parabola in ΔG° with its minimum at ΔG° = −λ. Increase the driving force past that minimum and you simply climb the other side of the parabola.
| Regime | Condition | Effect of more driving force | Where it matters |
|---|---|---|---|
| Normal | −ΔG° < λ | Rate increases (the intuitive behaviour) | Most ordinary redox reactions in solution |
| Activationless | −ΔG° = λ | ΔG‡ = 0; maximum rate | The optimum for a designed electron-transfer step |
| Inverted | −ΔG° > λ | Rate decreases | Suppressing wasteful back-reactions, e.g. charge recombination in photochemical systems |
The inverted region is not a curiosity. In a light-driven charge separation, the forward step is usually arranged to sit near the top of the curve while the unwanted recombination step, which is far more exergonic, falls into the inverted region and is therefore slow. That asymmetry is what allows a separated charge pair to survive long enough to do chemistry.
Self-exchange and the cross relation
For a self-exchange reaction — an ion transferring an electron to its own oxidised form, such as Fe2+ + Fe3+ — the products are identical to the reactants, so ΔG° = 0 exactly and the equation collapses to:
With λ = 1.00 eV, ΔG‡ = 1.00 ÷ 4 = 0.25 eV = 0.25 × 96.485 = 24.1 kJ mol−1. Note that a self-exchange reaction has a real barrier even though nothing appears to change chemically — the barrier is entirely the cost of reorganising the nuclei and solvent.
Self-exchange rate constants are experimentally accessible (by isotopic labelling or NMR line broadening), which makes them useful reference data. The Marcus cross relation then predicts the rate of a reaction between two different couples from their two self-exchange rates and the equilibrium constant:
where k11 and k22 are the self-exchange rate constants, K12 is the equilibrium constant for the cross reaction (obtainable from the two standard potentials), and f is a correction factor that is close to 1 when K12 is not extreme. Setting f ≈ 1 gives the simple working form used in most exam questions. The cross relation is one of the few genuinely predictive results in solution kinetics, and its successful test against many outer-sphere couples is a large part of why the theory was accepted.
Outer sphere versus inner sphere
Marcus theory as presented here describes outer-sphere transfer, where the coordination shells of both partners stay intact and the electron tunnels across. In an inner-sphere mechanism the two metals share a bridging ligand and the electron travels through it; the rate is then governed largely by the substitution chemistry of forming that bridge. Applying the Marcus equation to a clearly inner-sphere, bridged-intermediate reaction is a category error, and exam questions sometimes test exactly that distinction.
Common mistakes that cost marks
- Getting the sign of ΔG° wrong. Exergonic means negative. Substituting +0.40 instead of −0.40 in worked example 1 gives (1.40)²/4 = 0.49 eV instead of 0.090 eV — more than five times too large.
- Confusing λ with ΔG‡. They are equal only in the special case λ = 4ΔG‡ at zero driving force. λ is a property of the reorganisation, not the barrier itself.
- Assuming more exergonic is always faster. That is exactly what the inverted region contradicts, and it is the single most tested idea in this topic.
- Mixing eV and kJ mol−1 inside one calculation. Convert once, at the end, using 1 eV = 96.485 kJ mol−1.
- Forgetting that λ depends on the solvent. The same donor–acceptor pair can sit in the normal region in one solvent and in the inverted region in another, because changing the solvent changes λout.
- Using RT in J while ΔG‡ is in kJ. RT = 2478 J mol−1 = 2.478 kJ mol−1 at 298 K; a factor of 1000 here destroys the exponential.
- Applying the theory to inner-sphere or bond-breaking reactions without saying so. The simple parabolic model assumes the electron moves without bonds breaking.
Where this appears in exams
| Exam | Typical use |
|---|---|
| CSIR-NET | Identifying the regime from λ and ΔG°; inverted region reasoning; the cross relation |
| GATE (Chemistry) | Numerical answer type: ΔG‡ from λ and ΔG°, or λ from a measured barrier |
| IIT-JAM | Outer-sphere versus inner-sphere mechanisms at a qualitative level |
| M.Sc. / research work | Photoinduced charge separation, redox catalysis, electrode kinetics |
Check the current official syllabus and notification for your paper — the depth expected on electron-transfer theory varies between examinations.
Finish the calculation. Marcus gives you the barrier; converting that barrier into a rate constant, or comparing two barriers at a given temperature, is an Arrhenius-type step. Enter ΔG‡ as the activation energy along with your prefactor and temperature, and the Arrhenius calculator returns k — which is how the 3.0 × 1011 s−1 in worked example 1 was obtained.
Open the Arrhenius Equation Calculator →Preparing for CSIR-NET, GATE, IIT-JAM or CUET-PG? ABC Chemistry runs dedicated competitive-exam batches at the coaching centre and as live online classes for students across India — details at abcchemistry.in.