# 8.2.7: Marcus Theory

In classical transition-state theory, the expression for the rate constant of a bimolecular reaction in solution is

$$k = (\kappa \nu_{n})^{\frac{-\Delta G^{\ast}}{RT}}, \tag{6.16}$$

where $$\nu_{n}$$, the nuclear frequency factor, is approximately 1011 M-1s-1 for small molecules, and $$\Delta$$G* is the Gibbs-free-energy difference between the activated complex and the precursor complex. This theoretical framework provides the starting point for classical electron-transfer theory. Usually the transmission coefficient $$\kappa$$ is initially assumed to be unity. Thus, the problem of calculating the rate constant involves the calculation of $$\Delta$$G*, which Marcus partitioned into several parameters:

$$\Delta G^{\ast} = w^{r} + \left(\dfrac{\lambda}{4}\right) \left(1 + \dfrac{\Delta G^{o\; \prime}}{\lambda}\right)^{2}, \tag{6.17}$$

$$\Delta G^{o\; \prime} = \Delta G^{o} + w^{p} - w^{r} \ldotp \tag{6.18}$$

Here wr is the electrostatic work involved in bringing the reactants to the mean reactant separation distance in the activated complex, and wp is the analogous work term for dissociation of the products. These terms vanish in situations where one of the reactants (or products) is uncharged. $$\Delta$$G° is the Gibbs-free-energy change when the two reactants and products are an infinite distance apart, and $$\Delta$$G°' is the free energy of the reaction when the reactants are a distance r apart in the medium; $$\Delta$$G° is the standard free energy of the reaction, obtainable from electrochemical measurements (the quantity - $$\Delta$$G° is called the driving force of the reaction).

The reorganization energy $$\lambda$$ is a parameter that contains both inner-sphere ($$\lambda_{i}$$) and outer-sphere ($$\lambda_{o}$$) components; $$\lambda = \lambda_{i} + \lambda_{o}$$. The inner-sphere reorganization energy is the free-energy change associated with changes in the bond lengths and angles of the reactants. The $$\lambda_{i}$$ term can be evaluated within the simple harmonic-oscillator approximation:

$$\lambda_{i} = \left(\dfrac{1}{2}\right) \sum_{j} k_{j} (\Delta x_{j})^{2}, \tag{6.19}$$

where kj values are normal-mode force constants, and the $$\Delta$$xj values are differences in equilibrium bond lengths between the reduced and oxidized forms of a redox center.

The outer-sphere reorganization energy reflects changes in the polarization of solvent molecules during electron transfer:

$$\lambda_{o} = e^{2} \bigg[\left(\dfrac{1}{2r_{A}}\right) + \left(\dfrac{1}{2r_{B}}\right) - \left(\dfrac{1}{d}\right) \bigg] \bigg[\left(\dfrac{1}{D_{op}}\right) - \left(\dfrac{1}{D_{s}}\right) \bigg] ; \tag{6.20}$$

d is the distance between centers in the activated complex, generally taken to be the sum of the reactant radii rA and rB; Dop is the optical dielectric constant of the medium (or, equivalently, the square of the refractive index); and Ds is the static dielectric constant. This simple model for the effect of solvent reorganization assumes that the reactants are spherical, and that the solvent behaves as a dielectric continuum. (Sometimes the latter approximation is so rough that there is no correspondence between theory and experiment.)

Variations in $$lambda$$ can have enormous effects on electron-transfer rates. Some of the possible variations are apparent from inspection of Equation (6.20). First, $$\lambda_{o}$$ decreases with increasing reactant size. Second, the dependence of the reaction rate on separation distance attributable to $$\lambda_{o}$$occurs via the $$\frac{1}{d}$$ term. Third, $$\lambda\_{o}$$ decreases markedly as the solvent polarity decreases. For nonpolar solvents, Ds $$\simeq$$ Dop $$\simeq$$ 1.5 to 4.0. It is significant to note that protein interiors are estimated to have Ds $$\simeq$$ 4, whereas, Ds $$\simeq$$ 78 for water. An important conclusion is that metalloproteins that contain buried redox cofactors need not experience large outer-sphere reorganization energies.

The key result of Marcus theory is that the free energy of activation displays a quadratic dependence on $$\Delta$$G° and $$\lambda$$ (ignoring work terms). Hence, the reaction rate may be written as

$$k_{et} = (\nu_{n} \kappa)^{\frac{-(\lambda + \Delta G^{o})^{2}}{4 \lambda RT}}\ldotp \tag{6.21}$$

For intramolecular reactions, the nuclear frequency factor ($$\nu_{n}$$) is ~1013 s-1. One of the most striking predictions of Marcus theory follows from this equation: as the driving force of the reaction increases, the reaction rate increases, reaching a maximum at - $$\Delta$$G° = $$\lambda$$; when - $$\Delta$$G° is greater than $$\lambda$$, the rate decreases as the driving force increases (Figure 6.23). Figure 6.23 - Plot of log ket as a function of increasing driving force (- $$\Delta$$G°). Three $$\Delta$$G° regions are indicated: normal ( - $$\Delta$$G° < $$\lambda)$$; activationless ( - $$\Delta$$G° = $$\lambda$$); and inverted (- $$\Delta$$G° > $$\lambda$$). The corresponding two-well (ER, EP), diagrams also are shown. The dashed curve (kdiff) is for a bimolecular reaction. The predicted behavior of a unimolecular reaction (kuni) is given by the solid curve; here the rate could be as high as 1013 s-1, because it is not masked by diffusional processes.

Two free-energy regions, depending on the relative magnitudes of - $$\Delta$$G° and $$\lambda$$, are thus distinguished. The normal free-energy region is defined by - $$\Delta$$G° < $$\lambda$$A. In this region, $$\Delta$$G* decreases if - $$\Delta$$G° increases or if $$\lambda$$ decreases. If - $$\Delta$$G° = $$\lambda$$, there is no free-energy barrier to the reaction. In the inverted region, defined by - $$\Delta$$G° > $$\lambda$$, $$\Delta$$G* increases if $$\lambda$$ decreases or if - $$\Delta$$G° increases. Another widely used result of Marcus theory deals with the extraction of useful kinetic relationships for cross reactions from parameters for self-exchange reactions. Consider the cross reaction, Equation (6.22), for which the rate

$$A_{1}(ox) + A_{2}(red) \rightarrow A_{1}(red) + A_{2}(ox) \tag{6.22}$$

and equilibrium constants are k12 and K12, respectively. Two self-exchange reactions are pertinent here:

$$A_{1}(ox) + A_{1}(red) \rightarrow A_{1}(red) + A_{1}(ox) \tag{6.23a}$$

$$A_{2}(ox) + A_{2}(red) \rightarrow A_{2}(red) + A_{2}(ox) \tag{6.23b}$$

These reactions are characterized by rate constants k11 and k22, respectively. The reorganization energy ($$\lambda_{12}$$for the cross reaction can be approximated as the mean of the reorganization energies for the relevant self-exchange reactions:

$$\lambda_{12} = \frac{1}{2} (\lambda_{11} + \lambda_{22}) \tag{6.24}$$

Substitution of Equation (6.24) into Equation (6.17) leads to the relation

$$\Delta G_{12}^{\ast} = \frac{1}{2}(\Delta G_{11}^{\ast} + \Delta G_{22}^{\ast}) + \frac{1}{2}\Delta G_{12}^{\ast}(1 + \alpha), \tag{6.25a}$$

where

$$\alpha = \frac{\Delta G_{12}^{\ast}}{4(\Delta G_{11}^{\ast} + \Delta G_{22}^{\ast})}\ldotp \tag{6.25b}$$

When the self-exchange rates k11 are corrected for work terms or when the latter nearly cancel, the cross-reaction rate k12 is given by the Marcus cross relation,

$$k_{12} = (k_{11}k_{22}K_{12}f_{12})^{\frac{1}{2}}, \tag{6.26a}$$

where

$$ln f_{12} = \frac{(ln K_{12})^{2}}{4\; ln \left(\dfrac{k_{11}k_{22}}{\nu_{n}^{2}}\right)}\ldotp \tag{6.26b}$$

This relation has been used to predict and interpret both self-exchange and cross-reaction rates (or even K12, depending on which of the quantities have been measured experimentally. Alternatively, one could study a series of closely related electron-transfer reactions (to maintain a nearly constant $$\lambda$$12) as a function of $$\Delta$$G12; a plot of In k12 vs. In K12 is predicted to be linear, with slope 0.5 and intercept 0.5 In (k11k22). The Marcus prediction (for the normal free-energy region) amounts to a linear free-energy relation (LFER) for outer-sphere electron transfer.

## Cross Reactions of Blue Copper Proteins

Given the measured self-exchange rate constant for stellacyanin (k11 1.2 x 105 M 1s-1), the Marcus cross relation (Equation 6.26a) can be used to calculate the reaction rates for the reduction of CuII-stellacyanin by Fe(EDTA)2- and the oxidation of Cul-stellacyanin by Co(phen)33+. E°(Cu2+/+) for stellacyanin is 0.18 V vs. NHE, and the reduction potentials and self-exchange rate constants for the inorganic reagents are given in Table 6.3.66,67 For relatively small $$\Delta$$E° values, f12 is ~1; here a convenient form of the Marcus cross relation is log k12 = 0.5[log k11 + log k22 + 16.9$$\Delta$$E12°]' Calculations with k11, k22, and $$\Delta$$E12° from experiments give k12 values that accord quite closely with the measured rate constants.

$$Cu^{II}St + Fe(EDTA)^{2-} \rightarrow Cu^{I}St + Fe(EDTA)^{-}$$

$$k_{12}(calc.) = 2.9 \times 10^{5} M^{-1}s^{-1} \qquad (\Delta E_{12}^{o} = 0.06 V)$$

$$k_{12}(obs.) = 4.3 \times 10^{5} M^{-1}s^{-1} \qquad \qquad \qquad \qquad \qquad$$

$$Cu^{I}St + Co(phen)_{3}^{3+} \rightarrow Cu^{II}St + Co(phen)_{3}^{2+}$$

$$k_{12}(calc.) = 1.4 \times 10^{5} M^{-1}s^{-1} \qquad (\Delta E_{12}^{o} = 0.19 V)$$

$$k_{12}(obs.) = 1.8 \times 10^{5} M^{-1}s^{-1} \qquad \qquad \qquad \qquad \qquad$$

### Table 6.3 - Reduction potentials and self-exchange rate constants for inorganic reagents.

Reagent E°(V vs. NHE) k22 (M-1s-1)
Fe(EDTA)-/2- 0.12 6.9 x 104
Co(phen)33+/2+ 0.37 9.8 x 101

The success of the Marcus cross relation with stellacyanin indicates that the copper site in the protein is accessible to inorganic reagents. The rate constants for the reactions of other blue copper proteins with inorganic redox agents show deviations from cross-relation predictions (Table 6.4).68 These deviations suggest the following order of surface accessibility of blue copper sites: stellacyanin > plastocyanin > azurin. Rate constants for protein-protein electron transfers also have been subjected to cross-relation analysis.69

### Table 6.4 - Reactions of blue copper proteins with inorganic reagents.

a) M-1s-1
Protein Reagent k12 (obs.)a $$\Delta$$E12° k11 (obs.)a k11 (calc.)a
Stellacyanin Fe(EDTA)2- 4.3 x 105 0.064 1.2 x 105 2.3 x 105
Co(phen)33+ 1.8 x 105 0.186 1.2 x 105 1.6 x 105
Ru(NH3)5py3+ 1.94 x 105 0.069 1.2 x 105 3.3 x 105
Plastocyanin Fe(EDTA)2- 1.72 x 105 0.235 ~103 - 104 7.3 x 101
Co(phen)33+ 1.2 x 103 0.009 ~103 - 104 1.1 x 104
Ru(NH3)5py3+ 3.88 x 103 -0.100 ~103 - 104 4.9 x 104
Azurin Fe(EDTA)2- 1.39 x 103 0.184 2.4 x 106 2.8 x 10-2
Co(phen)33+ 2.82 x 103 0.064 2.4 x 106 7.0 x 103
Ru(NH3)5py3+ 1.36 x 103 -0.058 2.4 x 106 1.1 x 103