Unimolecular rate theory, RRKM

By choosing the initial conditions for an ensemble of trajectories to represent a quantum mechanical state, trajectories may be used to investigate state-specific dynamics and some of the early studies actually probed the possibility of state specificity in unimolecular decay [330]. However, an initial condition studied by many classical trajectory simulations, but not realized in any experiment is that of a micro-canonical ensemble [331] which assumes each state of the energized reactant is populated statistically with an equal probability. The classical dynamics of this ensemble is of fundamental interest, because RRKM unimolecular rate theory assumes this ensemble is maintained for the reactant [6,332] as it decomposes. As a result, RRKM theory rules-out the possibility of state-specific unimolecular decomposition. The relationship between the classical dynamics of a micro-canonical ensemble and RRKM theory is the first topic considered here. 

State specific experiments can now test unimolecular rate theories by probing microcanonical rate coefficients. Moore and coworkers [45] have studied the dissociation of ketene close to the reaction threshold in an attempt to test RRKM theory. 

Experimental studies have had an enormous impact on the development of unimolecular rate theory. A set of classical thermal unimolecular dissociation reactions by Rabinovitch and co-workers [6-10], and chemical activation experiments by Rabinovitch and others [11-14], illustrated that the separability and symmetry of normal modes assumed by Slater theory is inconsistent with experiments. Eor many molecules and experimental conditions, RRKM theory is a substantially more accurate model for the unimolecular rate constant. Chemical activation experiments at high pressures [15,16] also provided information regarding the rate of vibrational energy flow within molecules. Experiments [17,18] for which molecules are vibrationally excited by overtone excitation of a local mode (e.g. C-H or O-H bond) gave results consistent with the chemical activation experiments and in overall good agreement with RRKM theory [19]. 

Significance of Thermochemical Excitation Energy Distributions. Two fundamentally different approaches have been used for the elucidation of primary product excitation distributions for hot atom activated species. Many workers have utilized the RRKM unimolecular rate theory to invert pressure falloff data measured for the products from energetic H-for-H (19,64,65,74,83,84,85) or F-for-X (19,27) substitution reactions. More complete citations to the early literature are given elsewhere (57,58,59,77,80). This inversion procedure is computationally straightforward. However, it involves several a priori untested and... 

RRKM theory is a microcanonical transition state theory and as such, it gives the connection between statistical unimolecular rate theory and the transition state theory of thermal chemical reaction rates. Isomerization or dissociation of an energized molecule A is assumed in RRKM theory to occur via the mechanism... 

Detailed Cross-sections and Rates.—The RRKM version of transition-state theory for unimolecular reactions, as developed 25 years ago and sununarized in its useful practical form in recent books, has continued to find wide applications in unimolecular rate theory. As has been pointed out by Marcus in the 1973 Faraday Discussion on molecular beams, it is both a weakness and a strength of transition-state theory that it does not make very detailed statements on specific cross-sections and rates. With such information becoming accessible experimentally, more detailed statistical dynamical theories were to come. We have now four such detailed statistical approaches ... 

RRKM theory assumes a microcanonical ensemble of A vibrational/rotational states within the energy interval E E + dE, so that each of these states is populated statistically with an equal probability [4]. This assumption of a microcanonical distribution means that the unimolecular rate constant for A only depends on energy, and not on the maimer in which A is energized. If N(0) is the number of A molecules excited at / =... 

A situation that arises from the intramolecular dynamics of A and completely distinct from apparent non-RRKM behaviour is intrinsic non-RRKM behaviour [9], By this, it is meant that A has a non-random P(t) even if the internal vibrational states of A are prepared randomly. This situation arises when transitions between individual molecular vibrational/rotational states are slower than transitions leading to products. As a result, the vibrational states do not have equal dissociation probabilities. In tenns of classical phase space dynamics, slow transitions between the states occur when the reactant phase space is metrically decomposable [13,14] on the timescale of the imimolecular reaction and there is at least one bottleneck [9] in the molecular phase space other than the one defining the transition state. An intrinsic non-RRKM molecule decays non-exponentially with a time-dependent unimolecular rate constant or exponentially with a rate constant different from that of RRKM theory. 

In deriving the RRKM rate constant in section A3.12.3.1. it is assumed that the rate at which reactant molecules cross the transition state, in the direction of products, is the same rate at which the reactants fonn products. Thus, if any of the trajectories which cross the transition state in the product direction return to the reactant phase space, i.e. recross the transition state, the actual unimolecular rate constant will be smaller than that predicted by RRKM theory. This one-way crossing of the transition state, witii no recrossmg, is a fiindamental assumption of transition state theory [21]. Because it is incorporated in RRKM theory, this theory is also known as microcanonical transition state theory. 

Of course, in a thermal reaction, molecules of the reactant do not all have the same energy, and so application of RRKM theory to the evaluation of the overall unimolecular rate constant, k m, requires that one specify the distribution of energies. This distribution is usually derived from the Lindemann-Hinshelwood model, in which molecules A become activated to vibrationally and rotationally excited states A by collision with some other molecules in the system, M. In this picture, collisions between M and A are assumed to transfer energy in the other direction, that is, returning A to A ... 

RRKM theory, an approach to the calculation of the rate constant of indirect reactions that, essentially, is equivalent to transition-state theory. The reaction coordinate is identified as being the coordinate associated with the decay of an activated complex. It is a statistical theory based on the assumption that every state, within a narrow energy range of the activated complex, is populated with the same probability prior to the unimolecular reaction. The microcanonical rate constant k(E) is given by an expression that contains the ratio of the sum of states for the activated complex (with the reaction coordinate omitted) and the total density of states of the reactant. The canonical k(T) unimolecular rate constant is given by an expression that is similar to the transition-state theory expression of bimolecular reactions. 

Second, we calculate the unimolecular rate constant at the internal energy E via the RRKM theory. We use Eq. (7.54), where the rotational energy is neglected and where the sum and density of vibrational states are evaluated classically. Thus at E = 184 kJ/mol we get... 

The unimolecular rate constant k(E) is described within the framework of RRKM theory. In the following, we neglect the rotational energy in HCN as well as in the activated complex. The classical barrier height is Ec = 1.51 eV. 

Several conclusions can be drawn from Eqs. (76) and (77). First, the influence of fluctuations is the largest when the number of open channels u is of the order of unity, because then the distribution Q k) is the broadest. Second, the effect of a broad distribution of widths is to decrease the observed pressure dependent rate constant as compared to the delta function-like distribution, assumed by statistical theories [288]. The reason is that broad distributions favor small decay rates and the overall dissociation slows down. This trend, pronounced in the fall-of region, was clearly seen in a recent study of thermal rate constants in the unimolecular dissociation of HOCl [399]. The extremely broad distribution of resonances in HOCl caused a decrease by a factor of two in the pressure-dependent rate, as compared to the RRKM predictions. The best chances to see the influence of the quantum mechanical fluctuations on unimolecular rate constants certainly have studies performed close to the dissociation threshold, i.e. at low collision temperatures, because there the distribution of rates is the broadest. 

For some systems quasiperiodic (or nearly quasiperiodic) motion exists above the unimolecular threshold, and intrinsic non-RRKM lifetime distributions result. This type of behaviour has been found for Hamiltonians with low unimolecular thresholds, widely separated frequencies and/or disparate masses [12,62,65]. Thus, classical trajectory simulations performed for realistic Hamiltonians predict that, for some molecules, the unimolecular rate constant may be strongly sensitive to the modes excited in the molecule, in agreement with the Slater theory. This property is called mode specificity and is discussed in the next section. 

---