Unit 5: Temperature Dependence, Arrhenius Theory & Composite Kinetics

Determining the empirical order of a reaction is the vital first step toward elucidating its molecular mechanism. Four classical experimental protocols are employed.

1. The Method of Initial Rates

By measuring the initial rate $r_0 = (d[P]/dt)_{t \to 0}$ before back-reactions, product inhibition, or substantial reactant depletion occur:

Varying $[A]_0$ while holding $[B]_0$ constant:

2. The Isolation Method (Flooding)

All reactants except one are supplied in massive stoichiometric excess ($[B]_0, [C]_0 \gg [A]_0$). Their concentrations remain virtually constant throughout the reaction:

where $k_{\text{eff}} = k [B]_0^\beta [C]_0^\gamma$. The partial order $\alpha$ is then determined directly using standard integrated rate plots for species $A$. Repeating by flooding different components yields all partial orders.

3. The Fractional-Life Method

Measuring the time $t_{f}$ required for concentration to decrease by a fixed fraction $f$ (e.g., $f = 1/2$ for half-life, $f = 3/4$):

Plotting $\ln t_{1/2}$ versus $\ln[A]_0$ yields a straight line with slope $1 - n$.

4. Differential (van 't Hoff) Method

Directly taking natural logarithms of the differential rate equation:

A plot of $\ln r$ versus $\ln[A]$ gives slope $\alpha$.

Thermodynamic Activation Parameters Matrix for Prototypical Reactions

Applying Eyring Transition State Theory to rate constants across varied chemical mechanisms demonstrates how activation enthalpy ($\Delta H^\ddagger$), activation entropy ($\Delta S^\ddagger$), and activation free energy ($\Delta G^\ddagger$) govern reaction velocity at $298.15\text{ K}$:

| Reaction Mechanism | Prototypical Example | $k(298\text{ K})$ | $E_a$ ($\text{kJ/mol}$) | $\Delta H^\ddagger$ ($\text{kJ/mol}$) | $\Delta S^\ddagger$ ($\text{J}/(\text{mol}\cdot\text{K})$) | $\Delta G^\ddagger$ ($\text{kJ/mol}$) | |---|---|---|---|---|---|---| | Radical Recombination | $2 CH_3^\bullet \longrightarrow C_2H_6$ | $2.5 \times 10^{10}\text{ M}^{-1}\text{s}^{-1}$ | $0.0$ | $-2.5$ | $-52.0$ | $+13.0$ | | Bimolecular Radical Transfer| $OH^\bullet + CH_4 \longrightarrow H_2O + CH_3^\bullet$ | $6.4 \times 10^6\text{ M}^{-1}\text{s}^{-1}$ | $+18.5$ | $+16.0$ | $-65.0$ | $+35.4$ | | Unimolecular Isomerization | $\text{Cyclopropane} \longrightarrow \text{Propene}$ | $1.2 \times 10^{-15}\text{ s}^{-1}$ | $+272.0$ | $+269.5$ | $+40.0$ | $+257.6$ | | Bimolecular Gas Diels-Alder | $1,3\text{-Butadiene} + \text{Ethene} \longrightarrow \text{Cyclohexene}$ | $3.2 \times 10^{-18}\text{ M}^{-1}\text{s}^{-1}$| $+115.0$ | $+112.5$ | $-142.0$ | $+154.8$ | | Alkaline Ester Hydrolysis | $EtOAc + OH^- \longrightarrow AcO^- + EtOH$ | $0.11\text{ M}^{-1}\text{s}^{-1}$ | $+47.0$ | $+44.5$ | $-110.0$ | $+77.3$ | | Acid Sucrose Inversion | $\text{Sucrose} + H_3O^+ \longrightarrow \text{Glc} + \text{Fru}$ | $1.8 \times 10^{-4}\text{ M}^{-1}\text{s}^{-1}$ | $+108.0$ | $+105.5$ | $+32.0$ | $+96.0$ | | Enzyme Turnover ($k_{\text{cat}}$)| Carbonic Anhydrase ($CO_2 + H_2O$)| $1.0 \times 10^6\text{ s}^{-1}$ | $+38.0$ | $+35.5$ | $-10.0$ | $+38.5$ |

Chemical reaction rates are extraordinarily sensitive to temperature. As an approximate historical rule of thumb (the van 't Hoff rule), the rate of a typical homogeneous chemical reaction approximately doubles or triples for every $10^\circ\text{C}$ temperature rise ($Q_{10} \approx 2\text{--}3$).

Contrast with Kinetic Molecular Theory

From the kinetic theory of gases, the average molecular speed scales as:

and the binary collision frequency scales as:

For a $10^\circ\text{C}$ rise from $300\text{ K}$ to $310\text{ K}$:

The collision frequency increases by a mere $1.65\%$.

This colossal discrepancy between a $1.65\%$ increase in molecular collisions and a $200\%\text{--}300\%$ increase in chemical reaction rate proved conclusively that only an extraordinarily small, highly energetic fraction of molecular collisions possess sufficient energy to overcome a critical barrier and undergo chemical rearrangement.

In 1889, the Swedish chemist Svante Arrhenius synthesized J.H. van 't Hoff's thermodynamic equilibrium equation with Boltzmann's energy distribution to propose the foundational law of chemical kinetics.

Analogy to the van 't Hoff Isochore

van 't Hoff's thermodynamic relation for the temperature dependence of the equilibrium constant $K = k_1 / k_{-1}$ is:

Arrhenius proposed that each individual rate constant satisfies a similar differential equation:

where $E_a$ is the empirical activation energy (SI unit: $\text{J/mol}$ or $\text{kJ/mol}$).

Integrated Forms of the Arrhenius Equation

Assuming $E_a$ is independent of temperature over moderate intervals, indefinite integration yields:

where:

• $A$ is the pre-exponential factor (or frequency factor), possessing the same units as the rate constant $k$.
• $\exp(-E_a / R T)$ is the Boltzmann fraction of collisions with energy exceeding $E_a$.

Two-Temperature Comparative Form

Integrating between temperatures $T_1$ and $T_2$:

Arrhenius Plot Linearization

Plotting $\ln k$ versus $1/T$ (in $\text{K}^{-1}$):

• Slope $= -\frac{E_a}{R}$
• $y$-intercept $= \ln A$

The parameters $A$ and $E_a$ provide complementary microscopic insights into the reaction coordinate.

The Activation Energy ($E_a$)

The activation energy $E_a$ represents the minimum threshold energy colliding reactant molecules must possess along the line-of-centers (reaction coordinate) to induce bond distortion, overcome Coulombic electron repulsion, and reach the top of the potential energy barrier (the activated transition state complex).

Key energetic features:

• Reactions with low activation energies ($E_a < 20\text{ kJ/mol}$, such as radical-radical recombinations) proceed extremely rapidly and exhibit weak temperature sensitivity.
• Reactions with high activation energies ($E_a > 150\text{ kJ/mol}$, such as thermal cracking of alkanes) are slow at ambient temperature and exhibit extreme exponential acceleration with temperature.
• $E_a$ is always positive for elementary thermal reactions. (Certain composite reactions exhibit apparent negative activation energies, as analyzed in Section 5.6).

The Pre-Exponential Factor ($A$)

The pre-exponential factor $A$ quantifies the frequency of collisions with favorable spatial geometry:

where:

• $Z_0$ is the total binary collision frequency at unit concentration.
• $P$ is the steric factor (orientation probability), which accounts for the requirement that molecules collide with specific mutual orientations.

For simple gas-phase atom-atom or spherical collisions, $P \sim 1$ and $A \sim 10^{11}\text{ M}^{-1}\text{s}^{-1}$. For complex polyatomic molecules requiring precise alignment of functional groups, $P$ can be as small as $10^{-4}$ to $10^{-8}$, dramatically reducing the effective rate constant.

In 1920, Richard Chace Tolman provided the rigorous statistical mechanical definition of activation energy, establishing that $E_a$ is the difference between the average energy of reactive collisions and the average energy of all collisions.

Tolman's Theorem

Let $\sigma_R(E)$ be the microscopic reaction cross-section for collisions with center-of-mass collision energy $E$. The macroscopic bimolecular rate constant is given by the statistical ensemble thermal average:

where $\mu = \frac{m_A m_B}{m_A + m_B}$ is the reduced mass.

Taking the logarithmic temperature derivative:

Evaluating the derivative inside the integral yields Tolman's Principle:

where:

• $\langle E_R \rangle$ is the average energy of all collisions that successfully result in chemical reaction.
• $\langle E_{\text{all}} \rangle$ is the average energy of all collisions in the thermal ensemble ($= \frac{3}{2} R T$ for translational motion).

Tolman's theorem proves that the empirical Arrhenius activation energy $E_a$ is precisely the energy excess that reactive molecular encounters carry above the thermal average of the bulk reactant population.

University Honors Research Monograph: Marcus Electron Transfer Theory & The Inverted Region

Rudolph A. Marcus (Nobel Prize in Chemistry 1992) revolutionized chemical kinetics by developing the microscopic theory of outer-sphere electron transfer ($D + A \longrightarrow D^+ + A^-$):

• Parabolic Potential Energy Surfaces: The reactant ($D\cdots A$) and product ($D^+\cdots A^-$) states are modeled as multi-dimensional harmonic oscillator potential wells along a collective solvent polarization and nuclear reorganization coordinate $q$.
• Free Energy Barrier Formulation: The Gibbs activation free energy is related to the standard driving force $\Delta G^\circ$ and total reorganization energy $\lambda$:

where $\lambda = \lambda_{\text{in}} + \lambda_{\text{out}}$ accounts for internal bond length distortions ($\lambda_{\text{in}}$) and dielectric solvent dipole reorientations ($\lambda_{\text{out}}$).

• The Three Kinetic Regimes:

1. Normal Regime ($-\Delta G^\circ < \lambda$): Increasing thermodynamic driving force ($-\Delta G^\circ$) lowers $\Delta G^\ddagger$, accelerating electron transfer ($\ln k \propto -\Delta G^\circ$).

2. Barrierless Regime ($-\Delta G^\circ = \lambda$): $\Delta G^\ddagger = 0$, achieving the maximum possible activationless transfer rate.

3. The Marcus Inverted Region ($-\Delta G^\circ > \lambda$): Counter-intuitively, making the reaction even more thermodynamically favorable causes $\Delta G^\ddagger$ to increase, slowing down the rate!

• Experimental Proof: Gerhard Closs and John Miller (1984) confirmed the inverted region using rigid steroid spacer molecules, a discovery fundamental to photosynthesis and solar cell design.

When a chemical transformation proceeds via a multi-step composite mechanism, the overall observed rate constant $k_{\text{obs}}$ is an algebraic combination of elementary rate constants, and the apparent activation energy $E_{a, \text{app}}$ is a composite sum.

Pre-Equilibrium Followed by Slow Elementary Step

Consider the mechanism:

The observed rate law is:

Thus, the overall observed rate constant is:

Taking natural logarithms:

Differentiating with respect to temperature:

Multiplying by $R T^2$:

where $\Delta H_1^\circ = E_{a, 1} - E_{a, -1}$ is the standard enthalpy of formation of the intermediate pre-equilibrium complex.

Apparent Negative Activation Energies

If the initial pre-equilibrium step is strongly exothermic ($\Delta H_1^\circ < 0$) and its magnitude exceeds the activation energy of the decomposition step:

When $E_{a, \text{app}} < 0$, the reaction rate slows down as temperature increases! Famous example: gas-phase termolecular oxidation of nitric oxide ($2 NO + O_2 \longrightarrow 2 NO_2$), where pre-equilibrium dimerization $2 NO \rightleftharpoons N_2O_2$ is exothermic ($\Delta H^\circ \approx -10.5\text{ kJ/mol}$), yielding $E_{a, \text{app}} \approx -4\text{ kJ/mol}$.

While the classical Arrhenius equation is exceptionally robust over moderate temperature spans, significant non-linear curvature in $\ln k$ versus $1/T$ plots emerges across broad temperature regimes.

1. Temperature-Dependent Pre-Exponential Factor (Modified Arrhenius Equation)

Collision theory and Transition State Theory demonstrate that the pre-exponential factor is not strictly constant, but varies mildly with temperature:

where:

• $m = 1/2$ from hard-sphere collision theory ($A \propto \bar{v} \propto T^{1/2}$).
• $m = 1$ from classical Eyring transition state theory ($k_B T / h$).
• $m$ can be negative or higher integer depending on partition functions.

The true activation energy becomes temperature-dependent:

2. Low-Temperature Quantum Mechanical Tunneling

For reactions involving the transfer of light particles (protons $H^+$, hydrogen atoms $H^\bullet$, or electrons $e^-$), the de Broglie wavelength $\lambda_{\text{dB}} = h / \sqrt{2 m E}$ is comparable to the activation barrier width ($d \approx 0.5\text{--}1.0\text{ Å}$).

At low temperatures ($T < 200\text{ K}$), particles tunnel through the barrier rather than climbing over it. As $T \to 0\text{ K}$, the reaction rate flattens to a temperature-independent non-zero constant ($k \to k_{\text{tunnel}}$), causing the Arrhenius plot to curve upward with $E_a \to 0$.

3. Super-Arrhenius Behavior in Glassy and Viscous Media

In supercooled liquids and polymer glass transitions, relaxation rates decelerate far faster than Arrhenius predictions, described by the Vogel-Fulcher-Tammann (VFT) equation:

reflecting cooperative molecular rearrangements as free volume collapses near the ideal glass transition temperature $T_0$.

Investigating reaction dynamics under extreme hydrostatic pressures and across broad temperature ranges allows dissection of Transition State volumes, solvation rearrangements, and the physical reality of isokinetic relationships.

1. High-Pressure Chemical Kinetics & Activation Volume ($\Delta V^\ddagger$)

Hydrostatic pressure $P$ alters chemical equilibrium and reaction rate constants through mechanical work terms $P \Delta V$. From transition state thermodynamics:

where $\Delta V^\ddagger = V^\ddagger - V_R$ is the volume of activation, representing the difference in partial molar volume between the Transition State and the reactant state.

Structural and Solvational Components of $\Delta V^\ddagger$

The activation volume decomposes into two distinct physical contributions:

1. Intrinsic Volume Change ($\Delta V_{\text{intr}}^\ddagger$):

• For bond cleavage (dissociative mechanisms, unimolecular homolysis, $S_N1$): Bonds stretch to form the transition state, resulting in a volume expansion: $\Delta V_{\text{intr}}^\ddagger > 0$ ($+5 \text{ to } +20\text{ cm}^3/\text{mol}$).
• For bond formation (associative mechanisms, cycloadditions, $S_N2$): Molecules come together to form new covalent bonds, resulting in volume contraction: $\Delta V_{\text{intr}}^\ddagger < 0$ ($-10 \text{ to } -35\text{ cm}^3/\text{mol}$). In Diels-Alder reactions, $\Delta V^\ddagger \approx -30 \text{ to } -45\text{ cm}^3/\text{mol}$, causing dramatic rate accelerations ($> 1000\times$) at $P = 10\text{ kbar}$ ($1\text{ GPa}$).

2. Solvational Volume Change ($\Delta V_{\text{solv}}^\ddagger$, Electrostriction):

• When neutral reactants develop charge in the transition state (e.g., Menshutkin reaction: $R_3N + R'X \longrightarrow [R_3N^{\delta+}\cdots R'\cdots X^{\delta-}]^\ddagger$): Intense electric fields orient and pack dipolar solvent molecules tightly around developing charges. This electrostriction effect causes huge contraction: $\Delta V_{\text{solv}}^\ddagger \approx -20 \text{ to } -50\text{ cm}^3/\text{mol}$, accelerating the reaction under pressure.

2. Enthalpy-Entropy Compensation & The Isokinetic Temperature

When kinetic parameters are measured across a series of homologous reactions (varying catalyst ligands, substituents, or solvent mixtures), a linear correlation between activation enthalpy $\Delta H^\ddagger$ and activation entropy $\Delta S^\ddagger$ is frequently observed:

where $\beta$ is known as the isokinetic temperature (or isoequilibrium temperature).

Physical Reality vs. Statistical Artifact (Exner Analysis)

At $T = \beta$:

All reactions in the series are predicted to proceed at precisely identical rates at temperature $\beta$. However, because $\Delta H^\ddagger$ and $\Delta S^\ddagger$ are mathematically derived from the slope and intercept of the same Arrhenius/Eyring plot ($\ln(k/T)$ vs $1/T$), experimental errors in the slope propagate directly into the intercept, generating an apparent linear compensation artifact:

O. Exner demonstrated that genuine physical isokinetic behavior must be verified by plotting rate constants measured at two distinct temperatures directly against each other:

If the slope $b \ne 1$ and $b \ne T_1/T_2$ with statistical significance, a true common transition-state mechanism with identical solvation reorganization coordinates exists, and the true isokinetic temperature is calculated from: