79第9章

Kramers-Grote-Hynes 理論
9.1 遷移状態理論を超えて第 3 章で扱った遷移状態理論の近似は著しい理想化である. 多自由度系, 特に凝縮系では, 反応座標方向の運動は他の自由度からの摩擦を受け, 分割面近傍で逆反応方向に引き戻されることもある. これを再交差 (recrossing) と呼ぶ. 再交差する軌道の例を図 XXに示した. (a) は, 遷移状態理論においては反応性であると見なされてしまう. (b) は確かに反応性ではあるが, 遷移状態理論では二重に勘定されてしまう. 交差の回数が増えた場合も同様に考察すれば, 遷移状態理論は一般に反応速度を過大評価することが分かる.
遷移状態理論の反応速度と真の反応速度との比k exact = κ · k TST(9.1)

として定義される補正因子 κ は, 透過係数 (transmission coefficient) と呼ばれる. 本章では, 一般化 Langevin 方程式に基づいて透過係数を見積もる Grote-Hynes (GH) 理論の概略を紹介する. また, 旧来の Kramers 理論 (10.6 節で扱う) を, GH 理論の Langevin
極限 (摩擦の記憶効果を無視する極限) として議論する.

9.2

Grote-Hynes 理論

遷移状態近傍のポテンシャル面は, 反応座標 s に沿って逆放物線で近似されるとする.
2ωb,nas22
ここでは ωb,na を非断熱的障壁周波数 (non-adiabatic barrier frequency) と呼ぶ (下V (s) ∼=−記注意参照). 式 (8.5) と同様の調和熱浴モデル Hamilton 関数を考える.
2X
p2s ωb,naH=−s2 +22i

 2
Xpi
ωi2 2+xi + sci xi22i

80

第 9 章 Kramers-Grote-Hynes 理論■注意 ωb,na を非断熱的と呼ぶ意味は, 以下に示す断熱 (adiabatic) 極限との対比による.
上の Hamilton 関数で, 熱浴座標 xi のポテンシャルエネルギーに相当する最後の 2 項は次のように書き換えられる.
"#2X ωi2 cic2i 2xi + 2 s −s(9.2)
2ωi2ωi2i
よって, s の値を固定したとき, ポテンシャルの最小点は xi = −ci s/ωi2 で与えられる. xi が常にこの値をとると仮定することは, 熱浴が s の運動に断熱的に追随して s の各点で最小エネルギーを辿ると仮定することに相当する. これは, 電子が核の変位に直ちに適合するとする断熱近似に類似している. この断熱極限では上式 (9.2) の最終項が残り, ポテンシャル V (s) を実効的に

1−22ωb,na+X c2iiωi2!

1 2
1 2 2s2 = − (ωb,na+ ζ(0))s2 = − ωb,eqs22のように変化させることになる. ここで, 式 (8.10) を用いた. また, 最後の等号で断熱的障壁周波数を22ωb,eq= ωb,na+ ζ(0)(9.3)

で定義した. (添字 eq は, 熱浴が常に平衡 (equilibrium) にあることを意味する. ) 以下, 簡単のため ωb,eq を単に ωb と書く. 式 (8.12) で述べたように, ζ(0) の付加項はポテンシャルの谷底の周波数を減少させる一方, 上記のように山頂 (障壁) の周波数を増大させる.

第 8 章と同様, 反応座標 s の動力学は一般化 Langevin 方程式

s̈(t) = ωb2 s(t) −Z t0dτ ζ(t − τ )ṡ(τ ) + R(t)(9.4)

により記述される. これの Laplace 変換は,λ2 s̃(λ) − λs(0) − ṡ(0) = ωb2 s̃(λ) − ζ̃(λ)(λs̃(λ) − s(0)) + R̃(λ)
並べ替えると,s̃(λ) =(λ + ζ̃(λ))s(0) + ṡ(0) + R̃(λ)λ2 − ωb2 + λζ̃(λ)(9.5)

8.8.3 節で見たように, Laplace 逆変換は, 次のように書かれる.
Xs(t) =Resλ=λi eλt se(λ)
i∈poles

Res は留数, 和は s̃(λ) の極について取る. すなわち,λ2 − ωb2 + λζ̃(λ) = 0により決まる s̃(λ) の極が, 実時間の動力学を決定する.

(9.6)

81

9.3 Langevin 極限 (Kramers 理論)
導出の詳細∗ は省略するが, 式 (9.6) の解を λr とすれば, 式 (9.1) のように反応速度を補正する透過係数が

κGH =λrωb=
ωbλr + ζ̃(λr )(9.7)となることを Grote と Hynes は示した. これの物理的な意味については, その Langevin
極限を次節で見た後に再び議論する.

9.3

Langevin 極限 (Kramers 理論)

歴史を遡り, 1941 年の Kramers 理論について調べる. ただし, ここでは元論文∗ の導出ではなく, Grote-Hynes 理論の Langevin 極限を考える. 後の 10.6 節で同じ問題を別の観点から取り上げる. 第 8 章で見たように, 式 (9.4) の Langevin 極限は

s̈(t) = ωb2 s(t) − ζ ṡ(τ ) + R(t)Z ∞ζ=ζ(τ )dτ = ζ̃(λ = 0)(9.8)

λ2 + ζλ − ωb2 = 0(9.10)(9.9)

0であり, Grote-Hynes 方程式 (9.6) は,という簡単な 2 次方程式になる. λ > 0 の解をとって†

λr =−ζ +を得る. これは, 摩擦の強い極限 ζ  ωb で,p
ζ 2 + 4ωb22λr '(9.11)

ωb2ζとなる. したがって, 透過係数の式 (9.7) は

κKR =ωbζ(9.12)となる. 一方, 弱い摩擦の極限 ζ  ωb では, 式 (9.11) は

λr ' ωb∗

導出については, 原論文 R. F. Grote and J. T. Hynes, J. Chem. Phys. 77, 3736 (1982) の他, B.
J. Gertner, K. Wilson, and J. T. Hynes, J. Chem. Phys. 90, 3537 (1989), Appendix.
∗H. A. Kramers, Physica 7, 284 (1941).
詳 細 な 講 義 ノ ー ト が 以 下 で 入 手 可 能.
http://www.chem.vu.nl/~zwan/lectures/kramers.pdf†
式 (9.10) は正と負の 2 解を持つ. 負の解は (Laplace 逆変換から見られるように) 速やかに減衰する項を与える. GH 理論では正の解を λr に採る.

82

第 9 章 Kramers-Grote-Hynes 理論と近似されるので, 式 (9.7) はκ'1を与える. すなわち, 摩擦による TST への補正は小さくなる. このように, 摩擦の強さ ζの関数としての透過係数 κ は, 小さな摩擦の極限では 1 であり, 摩擦が強くなるにつれてζ の逆数に比例して減衰する.
■補足 元々の Kramers 理論は, 分布関数の Fokker-Planck 方程式‡ から導かれた. それによると, 上記とは異なり, 摩擦の小さい極限では反応速度は摩擦係数に比例することが示された. これは,障壁を越えるためには反応座標が周りの媒質からエネルギーを貰う必要があることを反映している.
反応座標にエネルギーを与えるのはランダムな外力だが, その起源は摩擦力と同じく溶媒または媒質の熱運動であるから, 溶媒との結合が大き過ぎない場合には, 反応はそれによって促進される. しかし, ある臨界点を越えると摩擦は反応を抑制し始める. よって, 摩擦係数の関数としての反応速度は, その転回点で極大を持つ. これは, Kramers 転回 (turn-over) と呼ばれる.

9.4

Grote-Hynes 理論再考

上の結果をもとに, GH 理論の式 (9.7) の振舞いを考察する. 摩擦核の Laplace 変換

ζ̃(λ) は, 一般に λ が大きくなるにつれ減少する. これは, 反応座標の運動が速くなると媒質から受ける摩擦は小さくなることを表す. 定性的には, 速い運動には媒質が追随できない, もしくは媒質が摩擦を及ぼす前にすり抜けてしまうという描像である. よって, 障壁周
波数 ωb が大きいときは摩擦が小さくなり, κGH ' 1 となる. これに対し, ωb が小さくなるにつれて摩擦の影響は強く出始め, κGH < 1 となる.

■補 足 GH 方 程 式 (9.6) の 左 辺 を f (λ) と お く. f (λ) = 0 の 解 を λ > 0 で 探 す. ま ず,f (0) = −ωb2 < 0. また, 一般に λ > 0 において ζ̃(λ) > 0 としてよい§ ので, f (ωb ) = ωb ζ̃(ωb ) > 0.
よって f (λ) = 0 の解は, 0 < λ < ωb にある. この範囲で f (λ) は単調で滑らかとし, そのグラフの曲線を C とする. 上述のように, ωb が大きいとき ζ̃(ωb ) は小さいので, 曲線 C のグラフは下方へシフトし, 横軸との交点は λr ' ωb に見出される. よって, κGH = λr /ωb ' 1. 逆に, ωb が小さくなるにつれ ζ̃(ωb ) が大きくなるので, 曲線 C が上方へシフトし, 交点 λr は左にずれて κGH < 1 を与える. このように, GH 理論は上段で直感的に述べた定性的振舞いを定量的に記述している.

‡§

Fokker-Planck 方程式については, 10.3 節で扱う.
ζ(t) はランダム力の時間相関関数を温度で割ったもの (揺動散逸定理) なので, ζ(t) の Fourier 変換はランダム力のパワースペクトルに比例し (Wiener-Khinchin の定理) その係数は正である. よって, ζ(t)を Fourier-cosine 級数で表したときの係数は正であり, その Laplace 変換も λ > 0 で正となる. 例えば,K. Ando, J. Chem. Phys. 101, 2850 (1994).