Unimolecular RRKM Theory

Introduction

A hot molecule created by a collision or light absorption does not necessarily have a single thermal temperature. Its total internal energy may be approximately known, and it may dissociate or isomerise before fully equilibrating with a bath. RRKM theory describes the unimolecular reaction rate as a function of that energy by comparing the number of accessible transition-state states with the density of states in the reactant. It connects molecular vibrations and rotations to pressure-dependent gas kinetics.

Core explanation

RRKM stands for Rice–Ramsperger–Kassel–Marcus. It refines earlier statistical unimolecular rate ideas by counting quantum states of reactant and activated complex, including vibrational and rotational structure and zero-point thresholds. In a common simplified expression, k(E) = N‡(E−E0) / [hρ(E)] , where E is reactant internal energy, E0 is the threshold, N‡ is the number of activated-complex states available up to E−E0 excluding the reaction-coordinate mode, ρ(E) is the reactant density of states near E, and h is Planck's constant. Symmetry, degeneracy and angular-momentum restrictions can modify the expression for a specific system.

The numerator counts ways to reach the dividing surface at the available energy. The denominator counts how many reactant states share the same energy interval. A larger N‡ tends to increase the rate, while a larger reactant density spreads population over more internal states and tends to decrease the rate for fixed numerator. The formula's units work because ρ has units of states per energy, so hρ has units of time when energy is measured consistently.

Below the threshold E0, the simplest over-barrier RRKM rate is zero because no activated-complex state is accessible. This does not rule out tunnelling or other quantum pathways. Just above threshold, the number of accessible transition-state states may rise slowly. At higher energy, more vibrational and rotational combinations open and k(E) generally changes. A thermal rate constant can be obtained by averaging energy-resolved rates over an appropriate energy distribution, but that averaging requires a defined bath and population model.

The central statistical assumption is that energy redistributes among accessible intramolecular modes sufficiently rapidly that the reactant behaves as an equilibrated microcanonical ensemble before it reacts. If energy remains localised in one bond or a trajectory escapes quickly along a special direction, the observed behavior can deviate from RRKM. Such nonstatistical dynamics are not proved by one failed fit; uncertainties in surface barriers, state counts, collisions and competing channels must also be checked.

For multiple unimolecular channels from the same energized species, each channel has its own threshold and transition-state state sum. At a fixed E, its RRKM rate competes with the others, giving a simple microcanonical branching model k1(E)/[k1(E)+k2(E)] if no other processes intervene. In a gas bath, collisions can remove or add internal energy before reaction; then a master equation couples energy grains and channels. Pressure-dependent falloff arises from the competition between chemical reaction and collisional redistribution.

A reliable RRKM calculation requires stationary structures and frequencies for reactant and each transition state, correct zero-point thresholds, rotational and symmetry treatment, and a model of anharmonic or hindered internal rotations when relevant. Low-frequency torsions can alter state counts substantially. As with other surface-based methods, the electronic structure error in E0 may dominate the predicted rate, particularly near threshold.

Step-by-step reasoning

Identify a unimolecular channel and its reactant minimum and verified saddle. Choose an internal energy E relative to the reactant ground-state convention. Calculate the zero-point-corrected threshold E0, reactant density ρ(E) and transition-state cumulative state count N‡(E−E0), omitting the unstable reaction mode. Form their ratio divided by h and include justified degeneracy factors. Compare k(E) with energy-randomisation and collisional timescales before interpreting measured kinetics.

Visual explanation

Draw a reactant well containing many horizontal quantum-state lines near energy E and a saddle threshold E0 with fewer accessible transverse states above it. Label the reactant line density ρ(E) and the activated-complex cumulative count N‡(E−E0). Beside the diagram show k(E) starting near zero at threshold and increasing as more transition-state states become accessible.

Real-world analogy

An energized molecule has many ways to distribute its internal energy. RRKM compares the number of “gate-compatible” arrangements with all arrangements available in the reactant well. If energy mixes rapidly among them, the fraction reaching the gate helps set the rate. The analogy should not suggest random conscious choices: state counts are quantum statistical quantities and dynamics can violate the mixing assumption.

Real-world example

In mass spectrometry, a molecular ion formed with excess internal energy may fragment through several channels. An energy-resolved RRKM model can estimate how each channel's rate grows above its threshold, helping interpret fragment yields. The observed signal also depends on how the ion was energized, the instrument time window and collisions. A fragment not seen experimentally may have a low rate within that window rather than an impossible pathway.

Why?

Why use a microcanonical model? An isolated energized molecule can have a specified total energy without a thermal bath temperature. Why is the transition-state reaction mode excluded? It describes passage through the dividing surface, not a bound vibration counted among transverse states. Why does ρ(E) matter? More reactant states share population at that energy. Why check energy mixing? The statistical state ratio assumes access to those states before reaction.

Common misconception

RRKM is not simply the Eyring equation with E substituted for T. It is energy resolved and uses sums and densities of quantum states. Another misconception is that above-threshold energy must immediately cause fragmentation; the rate depends on how many activated-complex states are available relative to the reactant state density and on competing processes.

Worked example

Question: At one internal energy, channel A has an accessible transition-state sum NA‡ = 200 and channel B has NB‡ = 50. Both leave the same reactant, so they share ρ(E) and h. Ignoring other factors, what branching fraction is predicted for A?

Reasoning: RRKM rates are proportional to their respective state sums when their denominators are identical. Thus kA:kB = 200:50 = 4:1. The fraction through A is 200/(200+50) = 0.80. This assumes both thresholds and state counts were evaluated at the same E and that no collision or other channel intervenes.

Answer: The simple fixed-energy RRKM branching fraction for A is 0.80.

Quick check

1. What is the major statistical assumption behind a basic RRKM rate expression? Answer: Internal energy redistributes rapidly among accessible reactant states before unimolecular reaction.

Exam focus

Write k(E) as a transition-state sum of states divided by h times reactant density of states, and define every symbol. Distinguish microcanonical energy from canonical temperature. State the zero-point-corrected threshold and the need to exclude the unstable reaction mode from bound-state counting.

Advanced insight

Angular momentum is conserved in isolated gas-phase reactions, so a fully resolved theory may use k(E,J) rather than k(E) alone. Loose dissociation transition states can require variational placement of the dividing surface rather than a tight saddle. Anharmonic torsions and mode coupling affect both ρ and N‡. When reaction competes with collisions, master equations propagate populations across energy grains and can reveal nonexponential or pressure-dependent behavior absent from one k(E) value.

Summary

RRKM theory gives an energy-resolved statistical rate for unimolecular reactions. Its numerator counts accessible transition-state states above a threshold; its denominator counts reactant states per energy. The model assumes intramolecular energy randomisation before reaction. Threshold accuracy, state counting, competing channels and collisions determine whether its predictions apply to a real experiment.

Practice questions

1. What do N‡ and ρ represent in the simplified RRKM expression? Answer: N‡ is the cumulative activated-complex state count; ρ is the reactant density of states at the specified energy.

2. Why does the ordinary over-barrier rate vanish below E0 in the basic formula? Answer: No activated-complex state is energetically accessible below the threshold.

3. What can cause observed dynamics to deviate from RRKM assumptions? Answer: Slow intramolecular energy redistribution, prompt directed trajectories or inaccurate barriers and state counts.

4. How do collisions enter a pressure-dependent unimolecular problem? Answer: They transfer internal energy between molecules and couple energy-resolved reaction rates through a master equation.

Sources: IUPAC Gold Book, RRKM theory; Journal of Chemical Theory and Computation, microcanonical tunnelling and RRKM; ACS Physical Chemistry Au, statistical and nonstatistical reactivity.