Kinetic Isotope Effect

The Kinetic Isotope Effect (KIE) is the change in the rate of a chemical reaction when an atom in the reactant is replaced by one of its isotopes. For reactions involving the breaking of a bond to an isotopic atom, like a C-H bond replaced by a C-D bond, the KIE is primarily caused by the difference in zero-point energies (ZPE) of the vibrating bonds in the reactant molecule.

Zero-Point Energy Difference

The zero-point energy of a harmonic oscillator (a simple model for a vibrating bond) is given by the formula:

\( ZPE = \frac{1}{2}h\nu \)

where \(h\) is Planck's constant and \(\nu\) is the vibrational frequency.

For a C-H bond versus a C-D bond, the vibrational frequency \(\nu\) is different because of the different masses of Hydrogen (\(m_H\)) and Deuterium (\(m_D\)). The frequency is inversely proportional to the square root of the reduced mass (\(\mu\)) of the vibrating system (\( \nu \propto 1/\sqrt{\mu} \)). For a bond between carbon and hydrogen/deuterium, the reduced mass is approximately \(\mu \approx \frac{m_C m_X}{m_C + m_X}\), where X is H or D.

A common approximation for the ratio of stretching frequencies is \(\frac{\nu_H}{\nu_D} \approx \sqrt{\frac{m_D}{m_H}}\). Since the mass of deuterium (\(m_D \approx 2\)) is approximately twice the mass of hydrogen (\(m_H \approx 1\)), we get:

\( \frac{\nu_H}{\nu_D} \approx \sqrt{\frac{2}{1}} = \sqrt{2} \)

Given the stretching frequency for the C-H bond (\(\nu_H\)) is 3000 cm\(^{-1}\), the approximate frequency for the C-D bond (\(\nu_D\)) is:

\( \nu_D \approx \frac{\nu_H}{\sqrt{2}} = \frac{3000}{\sqrt{2}} \text{ cm}^{-1} \)

The difference in zero-point energy between the C-H and C-D bonds in the reactant is:

\( \Delta ZPE = ZPE_H - ZPE_D = \frac{1}{2}h\nu_H - \frac{1}{2}h\nu_D = \frac{1}{2}h(\nu_H - \nu_D) \)

Since the frequency is given in cm\(^{-1}\), we convert it to energy units by multiplying by the speed of light \(c\):

\( \Delta ZPE = \frac{1}{2}hc(\nu_H - \nu_D)_{\text{cm}^{-1}} \)

Isotope Effect Calculation

For a reaction like bond homolysis, where the bond is completely broken in the transition state, the transition state has zero bond order and effectively zero ZPE associated with the bond stretch. The primary KIE (\(k_H/k_D\)) is then primarily determined by the difference in reactant ZPEs:

\( \frac{k_H}{k_D} \approx e^{\frac{\Delta ZPE}{k_B T}} \)

where \(k_B\) is the Boltzmann constant and \(T\) is the temperature.

Let's calculate the exponent \(\frac{\Delta ZPE}{k_B T}\):

• Given:
• \( \nu_H = 3000 \text{ cm}^{-1} \)
• \( \nu_D \approx \frac{3000}{\sqrt{2}} \text{ cm}^{-1} \)
• \( h = 6.63 \times 10^{-34} \text{ Js} \)
• \( c = 3 \times 10^{10} \text{ cm/s} \)
• \( k_B = 1.38 \times 10^{-23} \text{ J/K} \)
• \( T = 298 \text{ K} \)

First, calculate the difference in frequencies in cm\(^{-1}\):

\( (\nu_H - \nu_D)_{\text{cm}^{-1}} = 3000 - \frac{3000}{\sqrt{2}} = 3000 \left(1 - \frac{1}{\sqrt{2}}\right) \approx 3000 (1 - 0.7071) \approx 3000 \times 0.2929 \approx 878.7 \text{ cm}^{-1} \)

Now, calculate \(\Delta ZPE\):

\( \Delta ZPE = \frac{1}{2}hc(\nu_H - \nu_D)_{\text{cm}^{-1}} \)

\( \Delta ZPE = \frac{1}{2} \times (6.63 \times 10^{-34} \text{ Js}) \times (3 \times 10^{10} \text{ cm/s}) \times (878.7 \text{ cm}^{-1}) \)

\( \Delta ZPE = 0.5 \times 6.63 \times 3 \times 878.7 \times 10^{-34} \times 10^{10} \text{ J} \)

\( \Delta ZPE \approx 8729 \times 10^{-24} \text{ J} = 8.729 \times 10^{-21} \text{ J} \)

Next, calculate \(k_B T\):

\( k_B T = (1.38 \times 10^{-23} \text{ J/K}) \times (298 \text{ K}) \)

\( k_B T \approx 411.24 \times 10^{-23} \text{ J} = 4.1124 \times 10^{-21} \text{ J} \)

Finally, calculate the exponent \(\frac{\Delta ZPE}{k_B T}\):

\( \frac{\Delta ZPE}{k_B T} \approx \frac{8.729 \times 10^{-21} \text{ J}}{4.1124 \times 10^{-21} \text{ J}} \approx 2.12 \)

This value is approximately equal to 2.

Therefore, the expected isotope effect \(k_H/k_D\) is:

\( \frac{k_H}{k_D} \approx e^{\frac{\Delta ZPE}{k_B T}} \approx e^2 \)

Expected Isotope Effect Value

Based on the calculation using the simplified frequency approximation and the given constants, the exponent is approximately 2. The kinetic isotope effect \(k_H/k_D\) for full bond homolysis is therefore expected to be approximately \(e^2\).