Arrhenius Equation

The Arrhenius equation is a fundamental mathematical expression in physical chemistry, chemical kinetics, and biophysics that quantifies the temperature dependence of chemical reaction rate constants. Formulated by the Swedish physicist and chemist Svante Arrhenius in 1889 building on thermodynamic equilibrium relations established by Jacobus Henricus van 't Hoff in 1884, the equation relates the macroscopically observed rate constant k of an elementary or multistep chemical reaction to the thermodynamic temperature T (in kelvin) and an energy barrier termed the activation energy (E_a):

k = A exp(−*E*_a / RT)

where A is the pre-exponential factor (or frequency factor, having the same units as k) and R = 8.314462 J mol⁻¹ K⁻¹ is the universal molar gas constant. When expressed on a per-molecule basis rather than a per-mole basis, the equation takes the form k = A exp(−*E*_a / k_BT), where k_B = R / N_A = 1.380649 × 10⁻²³ J K⁻¹ is the Boltzmann constant. Across disciplines ranging from heterogeneous catalysis and atmospheric photochemistry to cellular enzymes and soil respiration under climate change, the Arrhenius equation links molecular collision mechanics and statistical mechanics with macroscopic reaction rates.

The Exponential Boltzmann Factor and Variables

In the Arrhenius equation k = A exp(−*E*_a / RT), the rate constant k is decomposed into the product of two distinct physical quantities:

1. The Pre-Exponential (Frequency) Factor (A): Represents the asymptotic ceiling rate constant at infinite temperature (T → ∞, where exp(−*E*_a / RT) → 1), corresponding to the frequency of molecular collisions or vibrational attempts possessing the proper steric orientation to react. Its units match k: s⁻¹ for a first-order reaction, and L mol⁻¹ s⁻¹ for a second-order reaction.
2. The Exponential Boltzmann Factor (exp(−*E*_a / RT)): A dimensionless probability factor bounded between 0 and 1 that specifies the fraction of colliding reactant molecules (or thermal states in a canonical ensemble) possessing translational and internal energy equal to or greater than the molar activation energy E_a (typically 20 to 250 kJ mol⁻¹ for covalent bond rearrangements). Because E_a / RT sits inside the negative argument of the natural exponential function exp(−*E*_a / RT), modest increases in temperature T or reductions in E_a via catalysis produce large, nonlinear increases in k.

For example, at room temperature (T = 298.15 K, RT = 2.479 kJ mol⁻¹), an uncatalyzed reaction with E_a = 80 kJ mol⁻¹ has a Boltzmann fraction of exp(−80,000 / 2,479) = exp(−32.27) ≈ 9.66 × 10⁻¹⁵. Lowering E_a to 40 kJ mol⁻¹ via an enzyme raises the fraction to exp(−40,000 / 2,479) = exp(−16.14) ≈ 9.83 × 10⁻⁸, accelerating the reaction by exp(+16.14) ≈ 1.02 × 10⁷ (ten millionfold) at constant temperature.

Logarithmic Linearization and the Two-Temperature Form

Taking the natural logarithm of both sides of k = A exp(−*E*_a / RT) converts the nonlinear curve into a linear equation suitable for linear regression of experimental kinetic data:

ln(k) = ln(A) − (E_a / R) (1 / T)

When experimental values of ln(k) are plotted on the vertical axis against reciprocal temperature (1 / T, in K⁻¹) on the horizontal axis—a graph known as an Arrhenius plot—reactions obeying simple Arrhenius behavior yield a straight line with slope = −*E*_a / R (so E_a = −*R* × slope) and y-intercept (at 1 / T = 0) = ln(A). Evaluating the linearized equation at two discrete temperatures T₁ and T₂ (with rate constants k₁ and k₂) eliminates ln(A) and yields the two-temperature Arrhenius relation:

ln(k₂ / k₁) = (E_a / R) (1 / T₁ − 1 / T₂) = [ E_a (T₂ − T₁) ] / ( R T₁ T₂ )

Differentiating ln(k) with respect to T provides the formal differential definition of the experimental Arrhenius activation energy recommended by IUPAC: E_a = R T² (d ln(k) / d*T*) = −*R* (d ln(k) / d(1/T)).

Van 't Hoff Equilibrium and Hard-Sphere Collision Theory

Arrhenius derived his rate law in 1889 by analyzing the temperature sensitivity of acid-catalyzed sucrose inversion and linking reaction kinetics to Jacobus Henricus van 't Hoff's 1884 equation for the equilibrium constant K_eq = k_f / k_r: d ln(K_eq) / d*T* = Δ*H*° / (RT²). Arrhenius proposed that reactant molecules exist in thermal equilibrium with high-energy "active" molecules governed by a Boltzmann exponential factor exp(−*E*_a / RT).

During the late 1910s, Max Trautz (1916) and William Cudmore McCullagh Lewis (1918) grounded the Arrhenius equation in the kinetic theory of gases and the Maxwell-Boltzmann distribution. For a bimolecular gas-phase reaction A + B → Products, molecules are modeled as hard spheres with collision cross-section σ_AB = π (r_A + r_B)² and reduced mass μ = m_Am_B / (m_A + m_B). Integrating the Maxwell-Boltzmann velocity distribution over all collisions whose line-of-centers translational kinetic energy exceeds the threshold energy E₀ yields the collision-theory rate constant:

k(T) = p Z_AB(T) exp(−*E*₀ / RT) = p N_A σ_AB (8 k_B T / π μ)^1/2 exp(−*E*₀ / RT)

where Z_AB(T) ∝ T^1/2 is the molar collision frequency factor and p (typically 10⁻⁶ to 1) is the dimensionless steric factor, which accounts for the requirement that polyatomic reactants must collide in a specific geometric orientation for frontier molecular orbitals to overlap productively.

Eyring-Polanyi Transition State Theory

In 1935, Henry Eyring, Michael Polanyi, and Meredith Gwynne Evans formulated Transition State Theory (TST) using quantum statistical mechanics and potential energy surfaces. TST models a chemical reaction as a trajectory along a multidimensional potential energy surface from a reactant valley across a first-order saddle point—the transition state or activated complex ([AB]^‡)—into the product valley. Assuming a quasi-equilibrium K^‡ between reactants and the activated complex and a universal vibrational crossing frequency k_BT / h (where h = 6.626 × 10⁻³⁴ J s is the Planck constant; at 300 K, k_BT / h ≈ 6.25 × 10¹² s⁻¹), the Eyring equation expresses the rate constant in terms of the standard Gibbs free energy of activation (Δ*G*^‡ = Δ*H*^‡ − *T*Δ*S*^‡):

k = κ (k_B T / h) exp(−Δ*G*^‡ / RT) = κ (k_B T / h) exp(Δ*S*^‡ / R) exp(−Δ*H*^‡ / RT)

where κ is the transmission coefficient, Δ*H*^‡ is the enthalpy of activation, and Δ*S*^‡ is the entropy of activation. Comparing the derivatives d ln(k)/d*T* of the Arrhenius and Eyring equations for a unimolecular or solution reaction yields E_a = Δ*H*^‡ + RT and A = e κ (k_B T / h) exp(Δ*S*^‡ / R), showing that the pre-exponential factor A is governed by the activation entropy Δ*S*^‡.

Enzymatic Catalysis, Q₁₀ Scaling, andQuantum Tunneling

In biochemistry and physiology, reaction rates over physiological temperature intervals (0–40 °C) are summarized using the temperature coefficient Q₁₀, defined as the factor by which a rate increases when temperature is raised by 10 K:

Q₁₀ = ( k₂ / k₁ )^{10 / (T₂ − T₁)} = exp[ 10 E_a / ( R T₁ T₂ ) ]

Near 25–30 °C, an enzymatic reaction with E_a ≈ 50 kJ mol⁻¹ exhibits Q₁₀ ≈ 1.95. Biological enzymes obey the Arrhenius equation up to an optimal temperature T_opt (37–42 °C in mesophiles; 75–95 °C in hyperthermophiles), above which thermal disruption of non-covalent bonds unfolds the tertiary protein structure (Δ*H*_denat ≈ 200–600 kJ mol⁻¹). In global climate change—pioneered by Arrhenius in his 1896 calculation of surface warming from atmospheric carbon dioxide and the greenhouse effect—microbial decomposition of soil organic carbon accelerates with temperature according to E_a ≈ 50–100 kJ mol⁻¹.

Departures from linear Arrhenius plots diagnose specialized physical regimes: (1) quantum mechanical tunneling of electrons, protons (H⁺), or hydride ions (H⁻) through narrow barriers produces convex curvature at low temperatures and large primary kinetic isotope effects (k_H / k_D > 15–80); (2) super-Arrhenius glassy dynamics near glass transitions follow the Vogel-Fulcher-Tammann law k = A exp[−*B* / (T − T₀)]; and (3) multistep reactions with an exothermic pre-equilibrium can display negative apparent activation energies.

Why does a 10 °C temperature increase roughly double the rate of many chemical and biological reactions near room temperature?

In the two-temperature Arrhenius equation ln(k₂/k₁) = E_a(T₂ − T₁) / (RT₁T₂), setting T₁ = 295 K (22 °C), T₂ = 305 K (32 °C), and k₂/k₁ = 2.0 (so ln 2 = 0.693) yields an activation energy of E_a = 0.693 × 8.314 × 295 × 305 / 10 ≈ 51.8 kJ mol⁻¹. Because many uncatalyzed solution reactions and enzyme-catalyzed metabolic steps near room temperature possess activation energies in the 45 to 60 kJ mol⁻¹ window, a 10 K rise increases the high-energy Maxwell-Boltzmann tail exp(−*E*_a/RT) by roughly a factor of two (Q₁₀ ≈ 1.8 to 2.3), even though the mean thermal speed of molecules increases by less than 2 percent.

How can a chemical reaction have a negative apparent activation energy?

Elementary single-step reactions with a physical barrier cannot have a negative activation energy, but a multistep reaction mechanism consisting of a rapid, reversible exothermic pre-equilibrium (A + B ⇌ I, with equilibrium constant K₁ and negative enthalpy Δ*H*°₁ < 0) followed by a slower rate-determining conversion of the intermediate into products (I → P, with rate constant k₂ and small barrier E_a,2) has an overall observed rate constant k_obs = K₁k₂. Because the apparent activation energy is the sum E_a,obs = Δ*H*°₁ + E_a,2, whenever the pre-equilibrium is more exothermic than the second barrier is positive (|Δ*H*°₁| > E_a,2), raising the temperature shifts the pre-equilibrium backward away from intermediate I faster than it accelerates step 2 (as in the classic termolecular oxidation 2 NO + O₂ → 2 NO₂), causing the overall reaction to slow down as temperature increases.

How does the Eyring equation improve upon the empirical Arrhenius equation?

Whereas the 1889 Arrhenius equation k = A exp(−*E*_a/RT) treats the pre-exponential factor A and activation energy E_a as empirical fitting parameters extracted from a plot of ln(k) versus 1/T, the 1935 Eyring-Polanyi equation k = κ(k_BT/h) exp(Δ*S*^‡/R) exp(−Δ*H*^‡/RT) derives the rate constant from quantum statistical mechanics on a potential energy surface. By separating the free-energy barrier Δ*G*^‡ into the enthalpy of activation (Δ*H*^‡, bond stretching and electronic reorganization) and the entropy of activation (Δ*S*^‡, conformational ordering or disordering at the saddle point), an Eyring plot of ln(k/T) versus 1/T directly measures whether a reaction is slowed by an energetic barrier or by an entropic bottleneck.

1. Arrhenius, S. (1889). "Über die Reaktionsgeschwindigkeit bei der Inversion von Rohrzucker durch Säuren." Zeitschrift für Physikalische Chemie, 4U(1), 226–248. https://doi.org/10.1515/zpch-1889-0416
2. Eyring, H. (1935). "The activated complex in chemical reactions." The Journal of Chemical Physics, 3(2), 107–115. https://doi.org/10.1063/1.1749604
3. Laidler, K. J. (1984). "The development of the Arrhenius equation." Journal of Chemical Education, 61(6), 494–498. https://doi.org/10.1021/ed061p494
4. Klinman, J. P., & Kohen, A. (2013). "Hydrogen tunneling links protein dynamics to enzyme catalysis." Annual Review of Biochemistry, 82, 471–496. https://doi.org/10.1146/annurev-biochem-051710-133623