RRKM

The existence of the polyad number as a bottleneck to energy flow on short time scales is potentially important for efforts to control molecnlar reactivity rising advanced laser techniqnes, discussed below in section Al.2.20. Efforts at control seek to intervene in the molecnlar dynamics to prevent the effects of widespread vibrational energy flow, the presence of which is one of the key assumptions of Rice-Ramsperger-Kassel-Marcns (RRKM) and other theories of reaction dynamics [6]. 

RRKM fit to microcanonical rate constants of isolated tran.s-stilbene and the solid curve a fit that uses a reaction barrier height reduced by solute-solvent interaction [46],... 

Khundkar L R, Marcus R A and Zewail A H 1983 Unimolecular reactions at low energies and RRKM-behaviour isomerization and dissociation J. Phys. Chem. 87 2473-6... 

A3.12 Statistical mechanical description of chemical kinetics RRKM... 

In the statistical description of ununolecular kinetics, known as Rice-Ramsperger-Kassel-Marcus (RRKM) theory [4,7,8], it is assumed that complete IVR occurs on a timescale much shorter than that for the unimolecular reaction [9]. Furdiemiore, to identify states of the system as those for the reactant, a dividing surface [10], called a transition state, is placed at the potential energy barrier region of the potential energy surface. The assumption implicit m RRKM theory is described in the next section. 

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 / =... 

The rapid IVR assumption of RRKM theory means that a microcanonical ensemble is maintained as the A molecnles decompose so that, at any time t, k(E) is given by... 

A RRKM nnimolecnlar system obeys the ergodic principle of statistical mechanics[H],... 

Figure A3,12.2(a) illnstrates the lifetime distribution of RRKM theory and shows random transitions among all states at some energy high enongh for eventual reaction (toward the right). In reality, transitions between quantum states (though coupled) are not equally probable some are more likely than others. Therefore, transitions between states mnst be snfficiently rapid and disorderly for the RRKM assumption to be mimicked, as qualitatively depicted in figure A3.12.2(b). The situation depicted in these figures, where a microcanonical ensemble exists at t = 0 and rapid IVR maintains its existence during the decomposition, is called intrinsic RRKM behaviour [9].

Figure A3.12.2. Relation of state oeeupation (sehematieally shown at eonstant energy) to lifetime distribution for the RRKM theory and for various aetiial situations. Dashed eiirves in lifetime distributions for (d) and (e) indieate RRKM behaviour, (a) RRKM model, (b) Physieal eounterpart of RRKM model, (e) Collisional state seleetion. (d) Chemieal aetivation. (e) Intrinsieally non-RRKM. (Adapted from [9].)...

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. 

The above describes the fundamental assumption of RRKM theory regarding the intramolecular dynamics of A. The RRKM expression for k E) is now derived. 

The tenn (E-E ) is tire sum of states at the transition state for energies from 0 to E-E. Equation (A3.12.15) is the RRKM expression for the imimolecular rate constant. 

Figure A3.12.3. Hannonic RRKM imimolecular rate constants for C2Hj H+C2H4 dissociation classical state counting (solid curve), quantal state counting (dashed curve). Rate constant is in units of s and energy in kcal moK (Adapted from [17].)...

RRKM theory allows some modes to be uncoupled and not exchange energy with the remaining modes [16]. In quantum RRKM theory, these uncoupled modes are not active, but are adiabatic and stay in fixed quantum states n during the reaction. For this situation, equation (A3.12.15) becomes... 

The RRKM rate constant is often expressed as an average classical flux tlirough the transition state [18,19 and 20]. To show that this is the case, first recall that the density of states p( ) for the reactant may be expressed as... 

The RRKM rate constant written this way is seen to be an average flux tlirough the transition state. 

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. 

As a result of possible recrossings of the transition state, the classical RRKM lc(E) is an upper bound to the correct classical microcanonical rate constant. The transition state should serve as a bottleneck between reactants and products, and in variational RRKM theory [22] the position of the transition state along q is varied to minimize k E). This minimum k E) is expected to be the closest to the truth. The quantity actually minimized is N (E - E ) in equation (A3.12.15). so the operational equation in variational RRKM theory is... 

Variational RRKM theory is particularly important for imimolecular dissociation reactions, in which vibrational modes of the reactant molecule become translations and rotations in the products [22]. For CH —> CHg+H dissociation there are tlnee vibrational modes of this type, i.e. the C—H stretch which is the reaction coordinate and the two degenerate H—CH bends, which first transfomi from high-frequency to low-frequency vibrations and then hindered rotors as the H—C bond ruptures. These latter two degrees of freedom are called transitional modes [24,25]. C2Hg 2CH3 dissociation has five transitional modes, i.e. two pairs of degenerate CH rocking/rotational motions and the CH torsion. 

Variational RRKM calculations, as described above, show that a imimolecular dissociation reaction may have two variational transition states [32, 31, 34, 31 and 36], i.e. one that is a tight vibrator type and another that is a loose rotator type. Wliether a particular reaction has both of these variational transition states, at a particular energy, depends on the properties of the reaction s potential energy surface [33, 34 and 31]- For many dissociation reactions there is only one variational transition state, which smoothly changes from a loose rotator type to a tight vibrator type as the energy is increased [26],... 

The classical mechanical RRKM k(E) takes a very simple fonn, if the internal degrees of freedom for the reactant and transition state are assumed to be hamionic oscillators. The classical sum of states for s harmonic oscillators is [16]... 

Only in the high-energy limit does classical statistical mechanics give accurate values for the sum and density of states tenns in equation (A3.12.15) [3,14]. Thus, to detennine an accurate RRKM lc(E) for the general case, quantum statistical mechanics must be used. Since it is difficult to make anliannonic corrections, both the molecule and transition state are often assumed to be a collection of hannonic oscillators for calculating the... 

Regardless of the nature of the intramolecular dynamics of the reactant A, there are two constants of the motion in a nnimolecular reaction, i.e. the energy E and the total angular momentum j. The latter ensures the rotational quantum number J is fixed during the nnimolecular reaction and the quantum RRKM rate constant is specified as k E, J). 

For a RRKM calculation without any approximations, the complete vibrational/rotational Flamiltonian for the imimolecular system is used to calculate the reactant density and transition state s sum of states. No approximations are made regarding the coupling between vibration and rotation. Flowever, for many molecules the exact nature of the coupling between vibration and rotation is uncertain, particularly at high energies, and a model in which rotation and vibration are assumed separable is widely used to calculate the quantum RRKM k(E,J) [4,16]. To illustrate this model, first consider a linear polyatomic molecule which decomposes via a linear transition state. The rotational energy for tire reactant is assumed to be that for a rigid rotor, i.e. 

It is straightforward to introduce active and adiabatic treatments of K into the widely used RRKM model which represents vibration and rotation as separable and the rotations as rigid rotors [41,42]. Eor a synnnetric top, tlie rotational energy is given by... 

If K is adiabatic, a molecule containing total vibrational-rotational energy E and, in a particular J, K level, has a vibrational density of states p[E - EjiJ,K). Similarly, the transition state s sum of states for the same E,J, and Kis [ -Eq-Ef(J,K)]. The RRKM rate constant for the Kadiabatic model is... 

In these models the treatment of K is the same for the molecule and transition state. It is worthwhile noting that mixed mode RRKM models are possible in which K is treated differently in the molecule and transition state [39],... 

---