第十章 分子反应动力学

理解碰撞、过渡态与活化过程,建立反应速率理论的微观图景。

6 节
智慧数字教材 |物理化学 · 预览版
✨ AI 未来课堂
CHAPTER 10

第十章 分子反应动力学

10.1

双分子反应的简单碰撞理论

阿伦尼乌斯根据实验从宏观的角度总结出了化学反应的动力学基本规律,即,指出反应的速率常数与温度的关系取决于活化能和指前因子。为什么会有这样的宏观规律存在?人们希望能从理论上或从微观的角度作出解释。更重要的是,人们希望像化学热力学能预言反应的方向和限度一样,化学动力学能够从分子水平上预言化学反应的动力学过程,直接从理论推导出化学反应的速率常数。经典的化学动力学理论包括:双分子反应的简单碰撞理论、、理论等。

是在阿伦尼乌斯关于"活化分子"和"活化能"概念的基础上,利用已经建立起来的气体分子运动论,在1918年由刘易斯建立起来的。简单碰撞理论(simple collision theory,SCT)把分子看作无结构的硬球,并假设两个反应物分子发生反应的先决条件是必须发生碰撞,但并非每一次碰撞都能导致反应发生。 在热平衡系统中,分子的平动能符合。如果互碰分子对的平动能不够大,则碰撞不会导致发生反应,碰撞后随即分离。只有那些相对平动能在分子连心线上的分量超过某一临界值的分子对,才能把平动能转化为分子内部的能量,使旧键断裂而发生原子间的重新组合。这种能导致旧键断裂的碰撞称为有效碰撞。碰撞理论认为只要知道分子的碰撞频率,再求出可导致旧键分裂的有效碰撞在总碰撞中的分数,即可求得反应速率和。

10.1.1

分子碰撞频率 ZABZ_{\mathrm{AB}}(硬球碰撞模型)

以气相双分子反应为例:

A+B⟶ 产物 \mathrm{A}+\mathrm{B} \longrightarrow \text { 产物 }

两个分子的碰撞过程实质上是在分子的作用力下发生的,两个分子先是相互接近,接近到一定距离时,它们之间开始产生斥力,斥力随分子间距离的减小而很快增大,之后分子就改变原来的方向而互相远离,于是就完成了一次碰撞过程。

研究双分子碰撞最简单的方法是固定其他所有分子的位置,只有被观察的分子 A 在运动。如图 10.1 所示,以 A 分子沿某一方向运动的距离为轴, A 分子和 B 分子的半径和 dABd_{\mathrm{AB}} 为半径的圆截面,所有运动着的 A 分子和 B分子,两者质心的投影落在半径为 dABd_{\mathrm{AB}} 的圆截面之内,都有可能发生碰

图10.1 双分子反应硬球碰撞模型
图10.1 双分子反应硬球碰撞模型

撞(图 10.2),其中 dABd_{\mathrm{AB}} 称为半径,数值上等于 A 分子和 B 分子的半径之和。 图10.2中虚线圆的面积称为(collision cross section),截面积 σ\sigma 数值上等于 πdAB2\pi d_{\mathrm{AB}}^{2} 。图10.3显示了硬球碰撞模型的剖面图,若 A 分子沿圆柱体轴向运动的平均速率为 vAv_{\mathrm{A}} ,则圆柱的体积为 πdAB2vA\pi d_{\mathrm{AB}}^{2} v_{\mathrm{A}} 。

图 10.2 分子间的碰撞和有效半径
图 10.2 分子间的碰撞和有效半径
图10.3 双分子反应硬球碰撞模型剖面图
图10.3 双分子反应硬球碰撞模型剖面图

设反应系统的体积为 VV ,体系中 A 与 B 的分子数分别为 NA,NBN_{\mathrm{A}}, N_{\mathrm{B}} ,则单位时间、单位体积内 A 分子与静止 B 分子的碰撞次数为

πdAB2vANBV\pi d_{\mathrm{AB}}^{2} v_{\mathrm{A}} \frac{N_{\mathrm{B}}}{V}

实际气相反应中分子 B 是运动的,因此以平均相对速率 vABv_{\mathrm{AB}} 代替分子 A 的平均速率 vA v_{\text {A }} 有

vAB=8RTπμ( 其中 μ 为折合质量 ,μ=MAMBMA+MB)v_{\mathrm{AB}}=\sqrt{\frac{8 R T}{\pi \mu}}\left(\text { 其中 } \mu \text { 为折合质量 }, \mu=\frac{M_{\mathrm{A}} M_{\mathrm{B}}}{M_{\mathrm{A}}+M_{\mathrm{B}}}\right)

则单位时间、单位体积内运动着的 A 与 B 分子的为

ZAB=πdAB28RTπμ(NBV)(NAV)Z_{\mathrm{AB}}=\pi d_{\mathrm{AB}}^{2} \sqrt{\frac{8 R T}{\pi \mu}}\left(\frac{N_{\mathrm{B}}}{V}\right)\left(\frac{N_{\mathrm{A}}}{V}\right)
(10.1)

其中, σ(πdAB2)\sigma\left(\pi d_{\mathrm{AB}}^{2}\right) 的量纲为 m2\mathrm{m}^{2} ,平均相对速率 8RTπμ\sqrt{\frac{8 R T}{\pi \mu}} 的量纲为 m⋅s−1\mathrm{m} \cdot \mathrm{s}^{-1} ,浓度 (NiV)\left(\frac{N_{i}}{V}\right) 的量纲为 m−3\mathrm{m}^{-3} ,因此 ZABZ_{\mathrm{AB}} 的量纲为 s−1⋅ m−3\mathrm{s}^{-1} \cdot \mathrm{~m}^{-3} 。

反应体系中所含反应物分子的数目等于阿伏伽德罗常数 LL 乘体系中该反应物的摩尔数 nn ,则该反应物在体系中的浓度为

NAV=LnAV=L[ A]NBV=LnBV=L[ B]\begin{aligned} & \frac{N_{\mathrm{A}}}{V}=\frac{L n_{\mathrm{A}}}{V}=L[\mathrm{~A}] \\ & \frac{N_{\mathrm{B}}}{V}=\frac{L n_{\mathrm{B}}}{V}=L[\mathrm{~B}] \end{aligned}

其中 [A][\mathrm{A}] 、 [B][\mathrm{B}] 为物质的量的浓度。 由此,在反应体系中,分子 A 与 B 每立方米、每秒碰撞的次数,即可表示为

ZAB=πdAB28RTπμ[ A][B]L2 m−3⋅ s−1Z_{\mathrm{AB}}=\pi d_{\mathrm{AB}}^{2} \sqrt{\frac{8 R T}{\pi \mu}}[\mathrm{~A}][\mathrm{B}] L^{2} \mathrm{~m}^{-3} \cdot \mathrm{~s}^{-1}
(10.2)

每立方米、每秒发生碰撞的摩尔组数为

ZABL=ZAB′=πdAB28RTπμ[ A][B]L mol⋅ m−3⋅ s−1\frac{Z_{\mathrm{AB}}}{L}=Z_{\mathrm{AB}}^{\prime}=\pi d_{\mathrm{AB}}^{2} \sqrt{\frac{8 R T}{\pi \mu}}[\mathrm{~A}][\mathrm{B}] L \mathrm{~mol} \cdot \mathrm{~m}^{-3} \cdot \mathrm{~s}^{-1}
(10.3)

对于同种分子之间的反应

2 A⟶ 产物 2 \mathrm{~A} \longrightarrow \text { 产物 }

dABd_{\mathrm{AB}} 是两个 A 分子的半径之和,即分子 A 的直径 dA。MAd_{\mathrm{A}}。 M_{\mathrm{A}} 为 A 分子的摩尔质量。可表示为

ZAA=22πdA28RTπMA[NAV]2=22πdA28RTπMA[ A]2L2 m−3⋅ s−1Z_{\mathrm{AA}}=\frac{\sqrt{2}}{2} \pi d_{\mathrm{A}}^{2} \sqrt{\frac{8 R T}{\pi M_{\mathrm{A}}}}\left[\frac{N_{\mathrm{A}}}{V}\right]^{2}=\frac{\sqrt{2}}{2} \pi d_{\mathrm{A}}^{2} \sqrt{\frac{8 R T}{\pi M_{\mathrm{A}}}}[\mathrm{~A}]^{2} L^{2} \mathrm{~m}^{-3} \cdot \mathrm{~s}^{-1}
(10.4)

同种分子每立方米、每秒发生碰撞的摩尔组数为

ZAAL=ZAA′=22πdA28RTπMA[ A]2L mol⋅ m−3⋅ s−1\frac{Z_{\mathrm{AA}}}{L}=Z_{\mathrm{AA}}^{\prime}=\frac{\sqrt{2}}{2} \pi d_{\mathrm{A}}^{2} \sqrt{\frac{8 R T}{\pi M_{\mathrm{A}}}}[\mathrm{~A}]^{2} L \mathrm{~mol} \cdot \mathrm{~m}^{-3} \cdot \mathrm{~s}^{-1}
(10.5)

然而,根据 ZAB,ZAAZ_{\mathrm{AB}}, Z_{\mathrm{AA}} 计算得到的 kk 值比实验值大很多,说明并不是每一次碰撞都能发生反应,即 ZAB,ZAAZ_{\mathrm{AB}}, Z_{\mathrm{AA}} 中只有一部分碰撞是能发生反应的。

10.1.2

有效碰撞分数 qq 的计算

两个分子的碰撞强度取决于它们的相对平动能。如图10.4中,法向的相对移动并不使两分子产生碰撞压力,径向(质心连线)上的相对移动才产生碰撞压力,即只有在质心连线上的碰撞能才对反应有贡献。只有那些相对平动能在分子连心线上的分量超过某一临界值的分子对,才能把平动能转化为分子内部的能量,使旧键断裂而发生原子间的重新组合。这种能导致旧键断裂的有

图 10.4 径向速率,连心线分量
图 10.4 径向速率,连心线分量

图10.4中,分子 A 和 B 的相对速率为 vABv_{\mathrm{AB}} ,相对平动能为 εr=12μvAB2\varepsilon_{\mathrm{r}}=\frac{1}{2} \mu v_{\mathrm{AB}}^{2} ,在径向 (质心连线方向)上的速率分量是 v1v_{1} ,径向相对平动能分量为

12μv12=12μvAB2cos⁡2θ=12μvAB2(1−sin⁡2θ)=12μvAB2(1−b2dAB2)=εr(1−b2dAB2)=εr′\begin{aligned} \frac{1}{2} \mu v_{1}^{2} & =\frac{1}{2} \mu v_{\mathrm{AB}}^{2} \cos ^{2} \theta \\ & =\frac{1}{2} \mu v_{\mathrm{AB}}^{2}\left(1-\sin ^{2} \theta\right) \\ & =\frac{1}{2} \mu v_{\mathrm{AB}}^{2}\left(1-\frac{b^{2}}{d_{\mathrm{AB}}^{2}}\right) \\ & =\varepsilon_{\mathrm{r}}\left(1-\frac{b^{2}}{d_{\mathrm{AB}}^{2}}\right) \\ & =\varepsilon_{\mathrm{r}}^{\prime} \end{aligned}

式中:θ\theta 是 vABv_{\mathrm{AB}} 与 dABd_{\mathrm{AB}} 的夹角;bb 称为碰撞参数,相当于通过 A 和 B 分子质心的与 vABv_{\mathrm{AB}}平行线间的垂直距离。

一个化学反应所需的最小能量,即发生反应的临界能,称为反应的阈能 εc\varepsilon_{\mathrm{c}} 。当双分子径向碰撞的能量 εr′\varepsilon_{\mathrm{r}}^{\prime} 大于 εc\varepsilon_{\mathrm{c}} 时,反应能发生;反之,则反应不能发生。即反应发生

的条件为

εr′=εr(1−b2dAB2)⩾εc\varepsilon_{\mathrm{r}}^{\prime}=\varepsilon_{\mathrm{r}}\left(1-\frac{b^{2}}{d_{\mathrm{AB}}^{2}}\right) \geqslant \varepsilon_{\mathrm{c}}
(10.6)

在双分子模型中用 EcE_{\mathrm{c}} 表示,它等于某化学反应中的反应物发生所需要达到的最低相对平动能 εc\varepsilon_{\mathrm{c}} 的 LL 倍,即

Ec=LεcE_{\mathrm{c}}=L \varepsilon_{\mathrm{c}}
(10.7)

如果假设反应速率比分子间能量传递的速率慢很多,即反应发生时分子的能量分布仍然遵守平衡时的玻尔兹曼能量分布定律,并且将分子的碰撞看作二维运动,则根据能量分布定律,能量在 EcE_{\mathrm{c}} 以上的分子数占总分子数的分数即为

q=e−Ec/(RT)q=\mathrm{e}^{-E_{\mathrm{c}} /(R T)}
(10.8)

例如,对某一反应有 Ec=186.2 kJ⋅ mol−1, T=300 K,q=3.7×10−33E_{\mathrm{c}}=186.2 \mathrm{~kJ} \cdot \mathrm{~mol}^{-1}, ~ T=300 \mathrm{~K}, q=3.7 \times 10^{-33} ,即每 2.7×10322.7 \times 10^{32} 次碰撞中有一次碰撞是。

这里需要指出的是,双分子简单中的反应阈能 EcE_{\mathrm{c}} 与阿伦尼乌斯 EaE_{\mathrm{a}} 的概念是不同的。中活化能 EaE_{\mathrm{a}} 是活化分子的平均能量和反应物分子的平均能量之差,而碰撞理论中的反应阈能 EcE_{\mathrm{c}} 是指反应物分子发生反应时所需要的最低能量。 EaE_{\mathrm{a}} 与 EcE_{\mathrm{c}} 的关系为

Ea=RT2 dln⁡k dT=RT2(12T+EcRT2)=Ec+12RTE_{\mathrm{a}}=R T^{2} \frac{\mathrm{~d} \ln k}{\mathrm{~d} T}=R T^{2}\left(\frac{1}{2 T}+\frac{E_{\mathrm{c}}}{R T^{2}}\right)=E_{\mathrm{c}}+\frac{1}{2} R T
(10.9)

在温度不是很高时,可近似认为 Ec≈EaE_{\mathrm{c}} \approx E_{\mathrm{a}} 。遗憾的是简单本身并不能预言 EcE_{\mathrm{c}} 的值,而需借助阿伦尼乌斯 EaE_{\mathrm{a}} 求算 qq ,这是简单碰撞理论的一大缺陷。

10.1.3

速率常数 k\boldsymbol{k} 的计算

根据双分子,反应速率表示为

r=ZAB′×q=ZAB′e−Ec/(RT)r=Z_{\mathrm{AB}}^{\prime} \times q=Z_{\mathrm{AB}}^{\prime} \mathrm{e}^{-E_{\mathrm{c}} /(R T)}
(10.10)

而根据的,反应速率为

r=−d[ A]dt=k[ A][B]r=-\frac{\mathrm{d}[\mathrm{~A}]}{\mathrm{d} t}=k[\mathrm{~A}][\mathrm{B}]
(10.11)

因此

r=−d[ A]dt=ZABLe−Ec/(RT)=LπdAB28RTπμe−Ec/(RT)[A][B]r=-\frac{\mathrm{d}[\mathrm{~A}]}{\mathrm{d} t}=\frac{Z_{\mathrm{AB}}}{L} \mathrm{e}^{-E_{\mathrm{c}} /(R T)}=L \pi d_{\mathrm{AB}}^{2} \sqrt{\frac{8 R T}{\pi \mu}} \mathrm{e}^{-E_{\mathrm{c}} /(R T)}[\mathrm{A}][\mathrm{B}]
(10.12)

得到为

k=LπdAB28RTπμe−Ec/(RT)=ZAB0e−Ec/(RT)k=L \pi d_{\mathrm{AB}}^{2} \sqrt{\frac{8 R T}{\pi \mu}} \mathrm{e}^{-E_{\mathrm{c}} /(R T)}=Z_{\mathrm{AB}}^{0} \mathrm{e}^{-E_{\mathrm{c}} /(R T)}
(10.13)

或 k=22LπdAA28RTπMAe−Ec/(RT)=ZAA0e−Ec/(RT)k=\frac{\sqrt{2}}{2} L \pi d_{\mathrm{AA}}^{2} \sqrt{\frac{8 R T}{\pi M_{\mathrm{A}}}} \mathrm{e}^{-E_{\mathrm{c}} /(R T)}=Z_{\mathrm{AA}}^{0} \mathrm{e}^{-E_{\mathrm{c}} /(R T)}(两个同种分子)

式(10.13)和式(10.14)可综合写为

k=Z0e−Ec/(RT)k=Z^{0} \mathrm{e}^{-E_{\mathrm{c}} /(R T)}
(10.15)

例10.1 600 K 时,反应 2NOCl⟶2NO+Cl22 \mathrm{NOCl} \longrightarrow 2 \mathrm{NO}+\mathrm{Cl}_{2} 的 kk 值为 60( mol⋅dm−3)−1⋅ s−160\left(\mathrm{~mol} \cdot \mathrm{dm}^{-3}\right)^{-1} \cdot \mathrm{~s}^{-1} ,实验为 105.5 kJ⋅ mol−1105.5 \mathrm{~kJ} \cdot \mathrm{~mol}^{-1} 。已知 NOCl 分子直径为 0.283 nm ,摩尔质量为 65.5 g⋅ mol−165.5 \mathrm{~g} \cdot \mathrm{~mol}^{-1} 。试计算反应在该温度下的速率常数。

解:

Ec=Ea−12RT=[105.5×103−(12×8.314×600)]J⋅ mol−1=103.0 kJ⋅ mol−1k=22LπdAA28RTπMe−Ec/(RT)=2×3.14×(2.83×10−10 m)2×6.022×1023 mol−1×8.314×600 J⋅ mol−13.14×65.5×10−3 kg⋅ mol−1exp⁡(−103000 J⋅ mol−18.314×600 J⋅ mol−1)=5.09×10−2( mol⋅ m−3)−1⋅ s−1=50.9( mol⋅dm−3)−1⋅ s−1\begin{aligned} E_{\mathrm{c}}= & E_{\mathrm{a}}-\frac{1}{2} R T=\left[105.5 \times 10^{3}-\left(\frac{1}{2} \times 8.314 \times 600\right)\right] \mathrm{J} \cdot \mathrm{~mol}^{-1}=103.0 \mathrm{~kJ} \cdot \mathrm{~mol}^{-1} \\ k= & \frac{\sqrt{2}}{2} L \pi d_{\mathrm{AA}}^{2} \sqrt{\frac{8 R T}{\pi M}} \mathrm{e}^{-E_{\mathrm{c}} /(R T)} \\ = & 2 \times 3.14 \times\left(2.83 \times 10^{-10} \mathrm{~m}\right)^{2} \times 6.022 \times 10^{23} \mathrm{~mol}^{-1} \times \\ & \sqrt{\frac{8.314 \times 600 \mathrm{~J} \cdot \mathrm{~mol}^{-1}}{3.14 \times 65.5 \times 10^{-3} \mathrm{~kg} \cdot \mathrm{~mol}^{-1}}} \exp \left(\frac{-103000 \mathrm{~J} \cdot \mathrm{~mol}^{-1}}{8.314 \times 600 \mathrm{~J} \cdot \mathrm{~mol}^{-1}}\right) \\ = & 5.09 \times 10^{-2}\left(\mathrm{~mol} \cdot \mathrm{~m}^{-3}\right)^{-1} \cdot \mathrm{~s}^{-1}=50.9\left(\mathrm{~mol} \cdot \mathrm{dm}^{-3}\right)^{-1} \cdot \mathrm{~s}^{-1} \end{aligned}

双分子反应的推导出的与实验测得的速率常数之间的接近程度从下面的两个例子中可以有一个概念性的了解。 (1) H2\mathrm{H}_{2} 与 Ar 的反应: k实验 =1.13×103 mol−1⋅ L⋅ s−1,kSCT=3.6×103 mol−1⋅ L⋅ s−1k_{\text {实验 }}=1.13 \times 10^{3} \mathrm{~mol}^{-1} \cdot \mathrm{~L} \cdot \mathrm{~s}^{-1}, \quad k_{\mathrm{SCT}}=3.6 \times 10^{3} \mathrm{~mol}^{-1} \cdot \mathrm{~L} \cdot \mathrm{~s}^{-1} (2) CH3\mathrm{CH}_{3}-与 H2\mathrm{H}_{2} 反应: k实验 =1.42×104 mol−1⋅ L⋅ s−1,kSCT=6.02×106 mol−1⋅ L⋅ s−1k_{\text {实验 }}=1.42 \times 10^{4} \mathrm{~mol}^{-1} \cdot \mathrm{~L} \cdot \mathrm{~s}^{-1}, \quad k_{\mathrm{SCT}}=6.02 \times 10^{6} \mathrm{~mol}^{-1} \cdot \mathrm{~L} \cdot \mathrm{~s}^{-1} 就大量反应的总体结果而言,对于结构极其简单的反应 k理论 ≈k实际 k_{\text {理论 }} \approx k_{\text {实际 }} ,二者较接近;而对于大多数反应 k理论 ≫k实际。这是由于双分子碰撞理论采用了过于简化的"硬球模 k_{\text {理论 }} \gg k_{\text {实际。这是由于双分子碰撞理论采用了过于简化的"硬球模 }}型",即把发生碰撞的分子看成刚性硬球。而实际上,分子是有内部结构的。为弥补 kk的偏差,碰撞理论引入了概率因子 PP 。

10.1.4

概率因子 PP

引入概率因子后,在双分子反应的的公式中须乘 PP ,即

k=PZ0e−Ec/(RT)k=P Z^{0} \mathrm{e}^{-E_{\mathrm{c}} /(R T)}
(10.16)

前面我们讨论了双分子中阈能与阿伦尼乌斯的关系,将式(10.9)代入式(10.16)得到

kSCT=PZ0e−Ec/(RT)=PZ0e1/2e−Ea/(RT)k_{\mathrm{SCT}}=P Z^{0} \mathrm{e}^{-E_{\mathrm{c}} /(R T)}=P Z^{0} \mathrm{e}^{1 / 2} \mathrm{e}^{-E_{\mathrm{a}} /(R T)}
(10.17)

从式(10.17)及 k=Ae−Ea/(RT)k=A \mathrm{e}^{-E_{\mathrm{a}} /(R T)} 可看出,阿伦尼乌斯公式中 AA 与双分子速率公式中有效碰撞频率是相互对应的,即

A=PZ0e1/2A=P Z^{0} \mathrm{e}^{1 / 2}
(10.18)

但是,简单本身不能求算 PP 值的大小。因此,概率因子 PP 是一个经验性的校正系数。概率因子的取值可以从 0∼10 \sim 1 ,包括了使真正数降低的各种因素造成的影响。比如某些基团的屏蔽效应使分子之间能量传递有过程,分子内部能量传递或再分配需要一定时间等因素。

双分子的成功之处在于揭示了反应速率的本质,提供了关于化学反应过程的简明而直观的图像。在运用分子运动论知识推算碰撞数的过程中,自然得出了与一致的结果。双分子碰撞理论对阿伦尼乌斯经验公式作出了理论解释,

预言阿氏活化能 EaE_{\mathrm{a}} 稍与 TT 有关,即 Ea=Ec+12RTE_{\mathrm{a}}=E_{\mathrm{c}}+\frac{1}{2} R T 。所以当反应过程中温度变化范围较大时,以 ln⁡k∼1/T\ln k \sim 1 / T 作图不能得到完好的直线。通过 kSCT k_{\text {SCT }} 与 kAk_{\mathrm{A}} 的比较,揭示了 AA 的物理含义,即指前因子 AA 与以反应物浓度为单位值时的有效有关。

双分子的不足之处显而易见,它不能完全脱离动力学实验数据而预言,即碰撞理论是一个半经验的理论。SCT 本身未能提出计算阈能 EcE_{\mathrm{c}} 的方法,需借助通过实验测定的反应活化能而求得;SCT 本身也无法给出概率因子 PP 的数值,也需要先由实验求得,然后再给出解释。双分子碰撞理论的成果离人们预言速率常数 kk 的期望还有很大差距,但它是最早提出的速率理论,它所提出的一些概念至今仍十分有用,它为我们描绘了一幅虽然粗糙但十分明确的反应图像。

10.1 双分子反应的简单碰撞理论