Monday, April 08, 2013

Reaction rate constants for the uni-, bi-, and tri-molecular reactions

Eyring equation: k = (kBT/h) e^(-delta G/RT)

The above equation can be use for uni-molecular directly. The simplified interpretation is that the reaction rate constant is proportional to the frequency of the reaction mode (kBT/h) and the relative concentrations of the transition state and the reactant(s):
[P(TS)/1 bar] / [P(reactant)/1 bar] = e^(-delta G/RT), where delta G is the standard Gibbs energy difference between the transition state and the reactant at temperature T and the standard state.

For a bi-molecular reaction, [P(TS)/1 bar] / {[P(A)/1 bar][P(B)/1 bar]} = e^(-delta G/RT),and again, the delta G is the standard Gibbs energy difference between the transition state and the two reactants at temperature T and the standard state..

For the bi- and tri-molecular reactions, we have learned that the reaction rate constants should be s^(-1)M^(-1) and s^(-1)M^(-2). Why? And how can we still use the Eyring equation to obtain the reaction rate constant given the Gibbs energy of activation?

Note that delta G refers to the difference between the STANDARD Gibbs energy of the transitions state (TS) and that of the reactant(s). The STANDARD state for gases is 1 bar!

Therefore, the Eyring equation can still be used to calculate the rate constant k for the bi- and tri-molecular reactions as long as we use PRESSURE/(1 bar) consistently in the expression of the reaction rate.

For example, A + B -> C

Reaction rate = d[P(C)/1 bar)]/dt = -d[P(A)/1 bar)]/dt = -d[P(B)/1 bar)]/dt = k[P(A)/1 bar)][P(B)/1 bar)]

Note that [P(A)/1 bar)], [P(B)/1 bar)], [P(C)/1 bar)] are unit-less. Pressures are divided by 1 bar because the delta G are usually determined at 1 bar. GAMESS calculations also give the values of delta G at 1 bar by default.

Therefore, we can see the left-hand-side (lhs) of the equation is in s^(-1) and so is the rhs.

d[P(C)/1 bar)]/dt = k[P(A)/1 bar)][P(B)/1 bar)]

What if the users of the above equation prefer to use M (mol/L)?
 Again, we can simply derive the equation in terms of concentrations.

d[P(C)/1 bar)]/dt = k[P(A)/1 bar)][P(B)/1 bar)]

d([C]RT/1 bar)/dt = k ([A]RT/1 bar) ([B]RT/1 bar)

 d([C])/dt = k(RT/1 bar) [A] [B]

Therefore the rate constant of a bi-molecular reaction becomes the product of 3 terms:
 the rate constant in s^(-1)M^(-1) = (kBT/h) e^(-delta G/RT) (RT/1 bar)

For a tri-molecular reaction, the rate constant = (kBT/h) e^(-delta G/RT) (RT/1 bar)^2 

Note that many rate constants of bi-molecular reactions are in the units of (s, cm, molecule),
We simply need to find the value of (RT/1 bar) in (s, cm, molecule)
 
For example, at 298.15K,
(RT/1 bar) = (8.3145 * 298.15) / 100000 = 0.024790 m^3/mol = 4.1164 * 10^(-20) cm^3/molecule


T(K)
kB T/p (cm^3/molecule)
100
1.38065E-20
200
2.76130E-20
298.15
4.11641E-20
300
4.14195E-20
400
5.52260E-20
500
6.90325E-20
600
8.28390E-20
700
9.66455E-20
800
1.10452E-19
900
1.24259E-19
1000
1.38065E-19
1100
1.51872E-19
1200
1.65678E-19
1300
1.79485E-19
1400
1.93291E-19
1500
2.07098E-19
1600
2.20904E-19
1700
2.34711E-19
1800
2.48517E-19
1900
2.62324E-19
2000
2.76130E-19
2100
2.89937E-19
2200
3.03743E-19
2300
3.17550E-19
2400
3.31356E-19
2500
3.45163E-19