Home
3803 words
19 minutes
遷移状態理論の礎:絶対反応速度論の定式化と統計力学的背景

last_modified: 2026-01-04

免責事項: 本記事は、Henry Eyringによる論文 J. Chem. Phys. 1935, 3, 107-115 を基に、生成AIによって作成された解説記事です。数式や論理の正確な内容は、必ず原著論文を参照してください。

1. 序論:化学反応速度論における「衝突」から「活性化」への転換#

1930年代初頭までの化学反応速度論において、反応速度を理論的に予測するための主要な枠組みは「衝突理論(Collision Theory)」であった。しかし、この理論は単純な気相反応においても実験値との間に無視できない乖離を示していた。本稿では、Henry Eyringが1935年に発表した論文に基づき、これらの課題を克服するために提唱された「遷移状態理論(Transition State Theory: TST)」または「絶対反応速度論(Absolute Reaction Rate Theory)」の基礎概念と数理的導出を解説する。

1.1 従来手法の課題:衝突断面積と立体因子の不確定性#

当時の標準的な手法では、二分子反応の速度は、反応する分子間の衝突回数に基づいて推定されていた 。この際、衝突頻度の算出には、運動量輸送の測定から得られる断面積(Cross-sectional area)が用いられていた 。しかし、この物理的な断面積は、実際に原子の組換え(質量の移動)を伴う「反応有効断面積」とは明確な相関を持たないという根本的な問題があった 。

具体的には、以下の要因が乖離の原因となっていた。

  1. 衝突の激しさ(エネルギー): 全ての衝突が反応を引き起こす十分なエネルギーを持っているわけではない 。
  2. 配向(Orientation): 運動量の交換には十分であっても、原子の交換に適した配向で衝突しているとは限らない 。

従来の理論では、この不一致を解消するために経験的な「立体因子(Steric factor)」または「配向因子(Orientation factor)」を導入していた 。この因子は通常 11 から 10110^{-1} の間であるが、反応によっては 10810^{-8} 程度まで小さくなる場合もあり、理論的な予測能力を著しく損なう要因となっていた 。

1.2 Eyringのアプローチ:ポテンシャル曲面と統計力学の融合#

Eyringは、原子間の力が量子力学(電子の運動と分布)によって決定される一方で、原子核自体の運動は、ひとたびその力が決定されれば、古典力学に従って記述できると仮定した 。これはBorn-Oppenheimer近似の概念に基づくものである。

彼の提案する核心は、ポテンシャルエネルギー曲面(Potential Energy Surface: PES)が既知であると仮定すれば、統計力学(または分子運動論)の手法を用いて、詳細な反応速度を第一原理的に計算できるという点にある 。このアプローチは、Herzfeld, Tolman, Fowlerらの先行研究や、PelzerとWignerによるオルソ・パラ水素変換の研究を発展させたものである 。


2. 理論的枠組み:活性錯合体の定義#

絶対反応速度論において最も重要な概念は「活性錯合体(Activated Complex)」である。これは反応物から生成物へと至る経路上の峠(Saddle point)に位置する準安定な化学種として定義される。

2.1 ポテンシャルエネルギー曲面上の運動#

原子系(反応系)の状態は、原子核間距離を座標とする多次元空間上の点として表現される 。

  • 安定な化合物: ポテンシャル曲面上の谷(極小点)に対応する 。
  • 反応: ある低い領域(反応物)から別の低い領域(生成物)への系の移動に対応する 。
  • 活性化状態: 反応経路(最低エネルギーパス)上の最も高い点、すなわち鞍点(Saddle point)に対応する 。

熱反応において、系はボルツマン因子に従い、エネルギー的に最も有利な経路(最低の峠)を経由して反応する 。

2.2 活性錯合体の性質#

活性錯合体(活性化状態にある系)は、通常の安定分子と同様の物理的性質(質量、慣性モーメント、振動数など)を持つが、一つの決定的な違いがある 。

  • 分解モード: 活性錯合体は鞍点に位置するため、通常の安定分子における振動モードのうちの1つが、負の曲率を持つ(虚振動に対応する)。これは、エネルギー障壁を越える方向への並進運動として扱われる 。

nn 原子の非線形活性錯合体の場合、自由度は以下のように配分される 。

  • 全自由度:3n3n
  • 重心の並進:3
  • 回転:3
  • 振動:3n73n - 7 (通常の分子は 3n63n-6 であるが、1つが反応座標に割り当てられるため)
  • 反応座標上の並進:1

この「反応座標上の並進」を統計力学的に処理することが、Eyringの理論の要である。なお、障壁頂上付近は平坦であるため、トンネル効果は通常無視できると仮定される 。


3. 数理的導出:絶対反応速度定数の定式化#

Eyringは、反応速度定数を、活性錯合体の濃度と、それが障壁を越えて生成物側へ移動する速度との積として定式化した 。

3.1 分配関数(Partition Functions)による平衡定数の記述#

統計力学によれば、ある状態の重み付き状態数(分配関数)は以下のように記述される。

並進の分配関数(1自由度あたり): 一般化座標 qiq_i と運動量 pip_i を用いて、

ftrans=1hdqiexp(pi22mikT)dpi=(2πmikT)1/2lihf_{trans} = \frac{1}{h} \int dq_i \int_{-\infty}^{\infty} \exp\left(\frac{-p_i^2}{2m_i kT}\right) dp_i = \frac{(2\pi m_i kT)^{1/2} l_i}{h}

ここで lil_i は長さの次元を持つ積分範囲である 。単位長さあたりでは li=1l_i=1 となる 。

振動の分配関数: 調和振動子近似を用いると、量子力学的分配関数は以下のようになる 。

fvib=(1exp(hνikT))1f_{vib} = \left(1 - \exp\left(\frac{-h\nu_i}{kT}\right)\right)^{-1}

古典極限(hνikTh\nu_i \ll kT)では kT/hνikT/h\nu_i となる 。

回転の分配関数: 非線形分子の場合:

frot=8π2(8π3ABC)1/2(kT)3/2h3f_{rot} = \frac{8\pi^2 (8\pi^3 ABC)^{1/2} (kT)^{3/2}}{h^3}

ここで A,B,CA, B, C は主慣性モーメントである 。

反応物と活性錯合体の間の平衡を考えると、活性錯合体の濃度 CC^\ddagger は、反応物の濃度と分配関数の比によって決定される。

C=(Reactants)×FFreactantsexp(E0kT)C^\ddagger = (\text{Reactants}) \times \frac{F^\ddagger}{F_{reactants}} \exp\left(\frac{-E_0}{kT}\right)

ここで E0E_0 は絶対零度における活性化エネルギーである 。

3.2 障壁通過速度と普遍頻度因子 kT/hkT/h#

ここで、活性錯合体の分配関数 FF^\ddagger から、反応座標(分解モード)に対応する並進の自由度を分離する。 反応座標上の並進の分配関数 ftransf^\ddagger_{trans} は、障壁頂上の微小な長さ δ\delta を考えると、

ftrans=(2πmkT)1/2δhf^\ddagger_{trans} = \frac{(2\pi m^* kT)^{1/2} \delta}{h}

となる。ここで mm^* は反応座標方向の有効質量である。

活性錯合体が障壁を越える平均速度 vˉ\bar{v} は、マクスウェル分布より以下のように計算される 。

vˉ=pˉm=0exp(p2/2mkT)(p/m)dpexp(p2/2mkT)dp=(kT2πm)1/2\bar{v} = \frac{\bar{p}}{m^*} = \frac{\int_0^\infty \exp(-p^2/2m^*kT) (p/m^*) dp}{\int_{-\infty}^\infty \exp(-p^2/2m^*kT) dp} = \left(\frac{kT}{2\pi m^*}\right)^{1/2}

注意すべきは、分子子の積分範囲が 00 から \infty である点である(順方向への通過のみを考えるため) 。

反応速度 RateRate は、活性錯合体の濃度(単位長さあたり)に通過速度を掛けたものである。

Rate=Cinternal×(2πmkT)1/2h×(kT2πm)1/2=Cinternal×kThRate = C^\ddagger_{internal} \times \frac{(2\pi m^* kT)^{1/2}}{h} \times \left(\frac{kT}{2\pi m^*}\right)^{1/2} = C^\ddagger_{internal} \times \frac{kT}{h}

ここで、CinternalC^\ddagger_{internal} は反応座標の並進自由度を除いた活性錯合体の濃度項を指す。 この計算により、質量 mm^* や長さ δ\delta が相殺され、非常にシンプルで普遍的な因子 kT/hkT/h が導出される 。

3.3 一般化された速度定数の式#

以上の考察より、任意の次数の反応に対する比速度定数 kk は以下の一般式で表される 。

k=κkThFFnexp(E0kT)k = \kappa \frac{kT}{h} \frac{F'^\ddagger}{F_{n}} \exp\left(\frac{-E_0}{kT}\right)

ここで:

  • κ\kappa (kappa): 透過係数(Transmission coefficient)。一度障壁を越えた系が戻らずに反応を完結する確率(通常は約1) 。
  • kT/hkT/h: 頻度因子。常温付近で約 6×1012s16 \times 10^{12} \text{s}^{-1} の値を持つ。
  • FF'^\ddagger: 反応座標の並進を除いた活性錯合体の分配関数 。
  • FnF_n: 反応物(通常状態)の分配関数 。
  • E0E_0: 活性化エネルギー 。

この式(論文中の式(10))は、反応速度が「遷移状態の熱力学的安定性(エントロピー項を含む分配関数)」と「エネルギー障壁(ボルツマン因子)」によって決定されることを示しており、従来の衝突理論の曖昧さを排除した厳密な形式となっている。


4. 具体的な反応系への適用#

Eyringは論文中で、いくつかの具体的な反応系に対してこの理論を適用し、数式を展開している。

4.1 3原子反応 (A+BCAB+CA + BC \to AB + C)#

直線型の活性錯合体 ABCA-B-C を経由する場合、速度定数 k1k_1 は以下のように具体化される [cite: 130-135]。

k1=κkThI3I2σ2σ3(m3m1m2)3/2(2πkT)3/2h3exp(E0kT)k_1 = \kappa \frac{kT}{h} \frac{\frac{I_3}{I_2} \frac{\sigma_2}{\sigma_3} \left(\frac{m_3}{m_1 m_2}\right)^{3/2} \frac{(2\pi kT)^{3/2}}{h^3} \dots}{\dots} \exp\left(\frac{-E_0}{kT}\right)

ここでは、反応物 AA(原子)、BCBC(二原子分子)、および活性錯合体 ABCABC(直線状三原子)の並進、回転、振動の各分配関数の比が詳細に計算される。 特に重要なのは、活性錯合体において、反応物の並進や回転の自由度の一部が、**「曲げ振動(Bending vibration)」**に変換される点である 。

並進の分配関数は 10810^8 程度の大きさを持つが、これが活性錯合体において「硬い」振動に変わると、その寄与は 11 に近くなる。この並進から振動への自由度の変換に伴う分配関数の減少(約 10810^{-8} のファクター)が、衝突理論における「衝突断面積」に対応する概念となる 。

4.2 4原子反応および単分子反応#

4原子反応(例:AB+CDAC+BDAB + CD \to AC + BD)や単分子分解についても同様の式が導出される 。 特に単分子反応においては、高圧極限において以下の形に帰着することが示される 。

k=κ(σσ)νiνiexp(E0kT)k = \kappa \left(\frac{\sigma}{\sigma^\ddagger}\right) \frac{\prod \nu_i}{\prod \nu_i^\ddagger} \exp\left(\frac{-E_0}{kT}\right)

この式は、PolanyiとWignerによって近似的に得られた式と一致しており、活性化エントロピーの変化が振動数の変化として現れていることを示唆している 。


5. 同位体効果の予測#

本理論の最も強力な検証材料の一つが、同位体効果(Isotope Effect)の予測である。 古典的な衝突理論では、同位体置換による反応速度の変化は、質量の変化に伴う衝突速度の違い(m1/2\sim m^{-1/2})としてのみ説明される 。例えば、水素(H)を重水素(D)に置換した場合、速度比は 21.41\sqrt{2} \approx 1.41 と予測される。

しかし、Eyringの式(式(2)や(5))は、量子化された振動エネルギー(零点エネルギー)の影響を明示的に含んでいる 。

vibrational term=1exp(hν/kT)1exp(hν/kT)\text{vibrational term} = \frac{1 - \exp(-h\nu/kT)}{1 - \exp(-h\nu^\ddagger/kT)}

活性錯合体と反応物では結合の強さが異なるため、同位体置換による零点エネルギーの変化量(12hΔν\frac{1}{2}h\Delta\nu)が異なる。これにより、活性化エネルギーの実効値が変化し、単なる質量効果(1.4倍)よりもはるかに大きな反応速度の違いが生じることが説明される 。この予測は実験事実とよく一致する。


6. 衝突理論との整合性と「立体因子」の正体#

最後に、Eyringは自身の理論がどのような条件下で従来の衝突理論(式(13))に帰着するかを示し、衝突理論における「立体因子」の物理的意味を明らかにした 。

Z=2N1N2(r1+r2)2(2πkTm1+m2m1m2)1/2Z = 2 N_1 N_2 (r_1 + r_2)^2 \left(2\pi kT \frac{m_1 + m_2}{m_1 m_2}\right)^{1/2}

Eyringの理論において、以下の条件が満たされるとき、結果は単純な衝突理論と一致する 。

  1. 反応に関与する分子の内部振動数が、活性錯合体になっても変化しない。
  2. 衝突によって失われる並進の自由度が、活性錯合体において「回転」または「極めて緩い(力の定数がゼロに近い)曲げ振動」に変換される。

逆に言えば、実験的に観測される立体因子(P<1P < 1)や異常な反応速度は、並進自由度が「硬い」振動自由度へと変換されることによるエントロピーの減少(分配関数の減少)として、統計力学的に厳密に説明できる 。例えば、重い原子を含む系で曲げ振動が硬い場合、反応速度は衝突理論の予測より 1/1001/100 程度になることが理論から直接導かれる 。


7. 結論と歴史的意義#

Henry Eyringによる1935年の本論文は、化学反応速度論におけるパラダイムシフトを成し遂げた。

  1. 第一原理的計算の可能性: ポテンシャルエネルギー曲面さえ分かれば、経験的なパラメータ(立体因子など)に頼らずに反応速度を絶対値として計算できる枠組みを提示した 。
  2. 普遍的な頻度因子: kT/hkT/h という普遍的な因子を導出し、反応速度の上限を物理定数と温度のみで規定した 。
  3. 量子効果の取り込み: 零点エネルギーの変化を通じて同位体効果を説明し、トンネル効果や非断熱遷移(本論文では除外されているが)への拡張の基礎を築いた 。

この理論は、反応物と遷移状態の構造(慣性モーメントや振動数)という、より具体的な物理量に基づいて反応速度を議論することを可能にし、その後の計算化学や反応動力学の発展の礎となった。活性錯合体は「分解しつつある分子」であるが、通常の分子と同様に統計力学的に扱えるという洞察は、化学反応の微視的理解における記念碑的な達成である。


参考文献#

  1. J. Chem. Phys. 1935, 3, 107.
遷移状態理論の礎:絶対反応速度論の定式化と統計力学的背景
https://ss0832.github.io/posts/20260103_physchem_transition_state_theory/
Author
ss0832
Published at
2026-01-04

Related Posts

遷移状態理論の系譜:熱力学・反応速度論・統計力学の融合と1935年の総合
2026-01-04
Keith J. LaidlerとM. Christine Kingによるレビュー論文『The Development of Transition-State Theory』(1983)に基づく学術的解説。1935年のEyring、Evans、Polanyiによる理論の確立に至るまでの、熱力学、気体分子運動論、統計力学という3つの異なる源流の発展過程と、その後の科学的受容について詳細に論じる。
溶液反応における遷移状態理論の展開:Evans-Polanyiによる熱力学的定式化と圧力効果の解明
2026-01-04
1935年のM. G. EvansとM. Polanyiによる論文『Some Applications of the Transition State Method...』の学術的解説。Eyringの統計力学的アプローチと対をなす、熱力学に基づいた反応速度式の導出、活性化体積の概念の導入、および高圧下での反応速度変化に関する理論的枠組みを詳細に論じる。
遷移状態理論の体系化と限界の定式化:Wignerによる1938年の総括と量子補正の導入
2026-01-04
Eugene Wignerによる1938年の論文『The Transition State Method』の包括的解説。反応速度論における遷移状態法の基礎となる3つの仮定(断熱性、古典性、平衡性)を明確化し、反応速度式を統計力学および熱力学から導出する過程、さらに量子効果(トンネル効果)による補正項の導入について、数理的背景と共に詳細に論じる。
遷移状態理論の現代的到達点:変分原理と量子効果による限界の克服
2026-01-04
Donald G. Truhlarらによる1996年の総説『Current Status of Transition-State Theory』に基づき、従来の遷移状態理論(TST)が抱える「再交差」や「トンネル効果」といった本質的な限界を、変分遷移状態理論(VTST)や多次元トンネル補正がいかにして克服したかを数理的・歴史的背景と共に解説する。
Claisen転位の反応速度論と微視的機構:気相実験と半経験的分子軌道法によるアプローチ
2026-01-02
F. W. SchulerとG. W. Murphyによる1950年の速度論的研究、およびM. J. S. DewarとC. Jieによる1989年の理論的研究に基づき、ビニルアリルエーテル類のClaisen転位における反応速度論的パラメータと、遷移状態の電子的性質(芳香族性対ビラジカル性)について詳説する。実験的活性化パラメータと計算化学による遷移状態構造の相関を論じる。
σ-ホール相互作用を駆動力とするトリクロロメチル基の還元的脱塩素化:チオフェノールによる反応機構の理論的・実験的解明
2026-01-05
Caballero-GarcíaらによるPhysical Chemistry Chemical Physics (2016) の論文に基づき、電子求引性基に結合したトリクロロメチル基における塩素原子の求電子性と、硫黄求核剤による還元メカニズムを解説する。