RRK and RRKM Theory in Outline

Energy randomisation among vibrational modes and fall-off curves

Lesson 3125 of 4,500 · Kinetics and Reaction Dynamics

Learning objectives

Introduction

The simple Lindemann-Hinshelwood model treats every energised A molecule as if it were the same. Real molecules carry a distribution of internal energies, and a molecule only a little above a reaction threshold may behave very differently from one far above it. Rice–Ramsperger–Kassel (RRK) and Rice–Ramsperger–Kassel–Marcus (RRKM) theories refine the model by asking how internal energy is distributed among molecular states and how rapidly a molecule at a given energy can reach products.

Core explanation

Let E be the internal energy of an isolated energised reactant and E₀ the minimum energy required for a specified reaction channel under the model. The central quantity is an energy-resolved or microcanonical rate k(E). Molecules below threshold have no classical over-barrier decomposition through that channel, while molecules above threshold can have very different rates depending on how much energy is accessible in the reactive motion. A thermal rate constant averages contributions from many E values; a fall-off curve further reflects how collisions populate and depopulate those energy levels.

RRK uses a deliberately simple energy-sharing picture. Imagine s equivalent vibrational modes rapidly exchanging energy. A reaction occurs when at least E₀ becomes concentrated in a particular reactive mode. Under a classical continuous-energy approximation, a common illustrative expression is k RRK(E) = ν(1 − E₀/E)^(s−1) for E ≥ E₀ and zero below threshold. Here ν is an attempted reaction frequency. As E approaches E₀ from above, the fraction becomes small; as E grows much larger than E₀, the rate approaches ν. The exponent reflects how many nonreactive modes compete for energy. This expression is a model, not a universal law for all molecules.

RRKM replaces the equivalent-mode picture with quantum state counting based on molecular frequencies and a transition-state structure. A standard schematic form is k RRKM(E) = N‡(E−E₀)/(hρ(E)), where ρ(E) is the reactant density of states and N‡ is the sum of transition-state states available with energy up to E−E₀, excluding motion along the reactive coordinate. Units make sense: N‡ is a count, h has energy-times-time units, and ρ has states per energy, so the ratio has reciprocal-time units. The IUPAC RRKM definition describes the dependence on active-state sums and reactant-state density. The IUPAC RRK definition emphasises energy distribution and localisation of critical energy.

Both approaches rely on a statistical assumption: intramolecular energy redistributes among relevant modes rapidly enough that the population of states at a given E can be treated statistically before decomposition. If energy remains trapped in a prepared vibration, two molecules with the same total E may react at different rates depending on how they were excited. Then a single k(E) may fail to capture the dynamics. Direct trajectory and spectroscopic experiments can test this assumption. A primary computational study of microcanonical rates explicitly discusses the rapid-energy-redistribution assumption behind RRKM treatments.

How does this connect to pressure? Collisions with a bath gas can add or remove internal energy. At low pressure, few collisions occur before a molecule reacts or cools, so the energy distribution differs from a fully thermalised high-pressure distribution. At high pressure, frequent energy-transfer collisions approach a thermal population and the rate approaches its high-pressure limit for a given temperature. Between the limits, a master equation can track populations in many energy bins, using k(E) for decomposition and collision probabilities for transfers between bins. The resulting fall-off curve is often broader than the single-A Lindemann formula because molecules across many internal energies participate.

This theory also clarifies why molecular size matters. A large molecule has many vibrational states that can store energy. At fixed total excess energy, a greater reactant density of states may dilute energy among many possibilities, while the transition-state state count can also grow. The net effect depends on both numerator and denominator of the RRKM expression; “more modes always slows reaction” is not a reliable universal rule. Competing product channels each have their own threshold and transition-state state sum, so branching ratios can change with E and with pressure.

RRKM is not an exact dynamical solution. It usually assumes rapid redistribution and an appropriate transition-state dividing surface; tunnelling, recrossing or nonstatistical dynamics may require corrections. The purpose of this page is to understand its architecture: energy-specific state counting combined with collisional energy redistribution to predict thermal and pressure-dependent behaviour. Detailed calculations need accurate molecular frequencies, barrier energies, angular momentum treatment and collision-energy-transfer models.

Step-by-step reasoning

1. Identify the reacting energised molecule, its internal energy E and threshold E₀ for the channel of interest. 2. Ask whether intramolecular vibrational redistribution is fast relative to product formation. 3. For a rough RRK picture, estimate how often enough energy reaches a reactive mode. 4. For RRKM, compare the number of accessible transition-state states with the reactant density of states at E. 5. Include collisions that move molecules among energy levels when predicting pressure dependence. 6. Compare the predicted fall-off and product branching with experiments, checking for nonstatistical behaviour.

Visual explanation

Draw a ladder of internal-energy bins for A with many closely spaced levels. Put a reaction threshold E₀ across the ladder. From levels above E₀, draw product arrows of different thickness, indicating that k(E) varies with energy. Add upward and downward arrows between bins for bath-gas collisions. Beside the ladder, sketch a high-pressure thermal distribution and a different low-pressure distribution. This shows why one A box cannot represent the full population.

Real-world analogy

A team needs enough resources concentrated at one workstation to complete a task. If resources are spread among many workstations, total resources alone do not guarantee completion at the required station. RRK's energy-sharing model uses a similar counting idea for energy localisation. The analogy stops there: molecular energy occupies quantum states and is redistributed by physical coupling, not deliberate decisions.

Real-world example

In combustion chemistry, an energised intermediate may either fragment or be stabilised by collisions with surrounding gas. At low pressure, it can fragment before many deactivating collisions; at high pressure, stabilisation may compete more effectively. An RRKM calculation supplies channel-specific k(E), and a master equation combines those rates with the bath-gas energy-transfer model. Comparing measured pressure-dependent product yields can test whether the chosen barrier energies and collision efficiencies are plausible.

Why?

Why does RRKM divide by the reactant density of states? At fixed E, the energised reactant can occupy many internal states. The numerator counts the transition-state routes available to pass toward products, while the denominator represents how many reactant states share the population. The ratio expresses reactive opportunity per occupied reactant state at that energy, with h setting the time scale in the statistical formula.

Common misconception

“RRKM says total energy above the barrier guarantees rapid reaction.” The energy must be distributed among states able to reach the reaction channel, and the rate can still be low just above threshold. Another mistake is treating the simple RRK formula as a precision RRKM calculation. RRK uses idealised equivalent modes; RRKM counts molecule-specific states. Both rely on statistical energy redistribution, which can fail when mode-specific dynamics persists.

Worked example

Use the illustrative RRK expression with s = 4 modes, ν = 1.0 × 10¹² s⁻¹, threshold E₀ = 50 kJ mol⁻¹ and internal energy E = 100 kJ mol⁻¹. Then 1−E₀/E = 0.50, so k(E) = 1.0 × 10¹²(0.50)³ = 1.25 × 10¹¹ s⁻¹. At E = 75 kJ mol⁻¹, the factor is 1−50/75 = 1/3 and k(E) ≈ 1.0 × 10¹²/27 = 3.70 × 10¹⁰ s⁻¹. The same molecule reacts faster in this model at higher internal energy. These values are illustrative microcanonical rates, not measured thermal rates averaged over a gas distribution.

Quick check

1. What does ρ(E) represent in the RRKM expression? Answer: The energised reactant's density of internal states near energy E, measured as states per energy interval.

Exam focus

Distinguish a thermal k(T,p) from a microcanonical k(E). Explain RRK as approximate energy localisation among modes and RRKM as transition-state state-sum divided by reactant state density. State the fast intramolecular-redistribution assumption and identify collisions as the bridge from energy-specific rates to pressure fall-off. Do not try to infer a full RRKM result from barrier height alone.

Advanced insight

Angular momentum conservation can restrict which transition-state states are accessible from a given reactant state, so detailed RRKM calculations may use k(E,J) rather than k(E) alone. Quantum tunnelling can add subthreshold transmission, and multiple wells can exchange population before product formation. If a molecule is prepared in one particular vibrational mode and reacts before its energy randomises, mode-specific trajectories may show non-RRKM lifetimes. Such findings do not invalidate statistical rate theory generally; they identify where its equilibration assumptions are too strong for the experiment.

Summary

RRK and RRKM theories describe unimolecular reaction rates as functions of internal energy. RRK estimates how likely energy is to concentrate in a reactive mode; RRKM counts accessible transition-state states relative to reactant states. Both assume sufficiently rapid intramolecular redistribution for statistical treatment. Combined with bath-gas energy transfer, energy-specific rates explain pressure-dependent fall-off more realistically than a single energised-A model.

Practice questions

1. In the simple RRK expression, what is k(E) when E < E₀ under the classical threshold assumption? Answer: Zero for that over-barrier channel, because the molecule lacks the model's required energy.

2. If ρ(E) rises while N‡(E−E₀) and other factors stay fixed, what happens to the RRKM rate? Answer: It decreases because more reactant states share the same accessible transition-state opportunities.

3. Why can two molecules with the same total E react at different rates when statistical assumptions fail? Answer: Their energy can occupy different modes and may not redistribute before reaction, so access to the reactive coordinate differs.

4. What two ingredients are combined in a pressure-dependent master-equation treatment? Answer: Energy-specific reaction rates such as k(E) and collision-driven transfer probabilities between internal-energy bins.

5. Why can total pressure affect a thermal rate even if k(E) for an isolated molecule is unchanged? Answer: Collisions alter how the molecular population is distributed among internal energies, changing the weighted average of the energy-specific rates.