化学反応のダイナミックスとカオス| 「遷移状態」概念再考 |神戸大学 理学部

小松崎 民樹 1(2003 年 8 月 6-10 日 第 43 回分子科学若手の会夏の学校予稿)

化学反応を制御するためには、反応の「静」と「動」の両面を理解する必要がある。「静」の側面とは物質設計に基礎を置く化学反応系のポテンシャルエネルギー曲面の設計であり、「動」の側面とは化学反応の速度論に基づいた動的な設計である。電子状態理論が円熟期を迎えつつあるのに対し、化学反応理論は未開拓な部分が多く残っているが、近年、遷移状態分光などの実験技術ならびにカオス力学系の進展により、「反応の「前」と「後」の中間領域である遷移状態において系はどのような動力学的特性を有しているか？」が次第に解明され、従来の化学反応理論を見直す必要性が問われている。すなわち、従来の反応理論では、反応座標はポテンシャルエネルギー面の幾何学情報から導かれ、運動量には陽に依存しないものとして捉えられてきたが、「なにを反応に直接携わる自由度、反応自由度、として捉えるべきなのか」という基本が改めて問われているのである。元来、Eyring, Wigner らによって導入された遷移状態とは系が反応する際に反応系から生成系へ至る過程において一回だけ交差する（一般には相空間の） 超曲面（反応分割面とも呼ぶ）" を指し、ポテンシャルエネルギー面の鞍部点（サドル）を意味するものではない。しかしながら、（反応のプロセスで一回しか通過しない）反応分割面は実存し得るのだろうか？よく知られる（Kramers-）Grote-Hynes 理論は、point-of-no-return としての遷移状態が真に存在するか否かを問う代りに、（与えられた反応分割面に対して）しばしば観測される再交差現象を（一般化）ランジュバン方程式におけるミクロな 摩擦 として解釈している。本稿では、 再交差を与えない反応分割面（遷移状態）は存在しえるのか、原点に立ち返って、なぜ系は始原系から生成系へ至りえるのか（＝なぜ反応が生起するのか） という問いに関する答えを真摯に求めたい。

目次

1 はじめに
2 化学反応を理解する上で必要なカオスの予備的基礎知識

2.1 作用─角変数 : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :2.2 正準変換摂動理論 : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :2.2.1 相空間のなかの一般化された基準振動概念 : : : : : : : : : : : : : : : : : :1 E-mail: tamiki@kobe-u.ac.jp

1

2
4

4
9
9

2.2.2 リ─正準変換摂動理論 : : : : : : : : : : : : : : : : : : : : : : : : : : : : : 11
2.2.3 小さい分母の問題 : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : 23
2.2.4 サドル領域における反応方向の作用不変性 : : : : : : : : : : : : : : : : : : 25
2.3 ポアンカレ断面 : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : 29
2.4 ハミルトン写像 : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : 36
2.5 セパラトリックス、安定・不安定不変多様体、ヘテロ（ホモ）クリニック交差 : : 43
2.5.1 ホモクリニック点の存在と 1 価の解析的積分の非存在 : : : : : : : : : : : : 49
2.5.2 相空間上の移送-大域的な遷移状態- : : : : : : : : : : : : : : : : : : : : : : 50
2.6 法双曲的不変多様体 : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : 51

3 鞍部領域の相空間構造

53

4 大域的な相空間構造

67

5 これからの化学反応動力学理論の課題
6 付録

68
73

7 おわりに
8 参考文献
9 問題集

83
84
85

3.1 遷移状態分光および原子クラスターの \相転移" に見た鞍部領域の可積分性 : : : : 53
3.2 遷移状態理論と Grote-Hynes 理論の等価性：反応における「系」とはなにか : : : 54
3.3 高次元相空間上の遷移状態：カオスの海に埋もれた規則的経路 : : : : : : : : : : : 55
4.1 Davis-Gray のセパラトリックス遷移状態理論 : : : : : : : : : : : : : : : : : : : :4.2 ホモクリニック交差の構造不安定性 : : : : : : : : : : : : : : : : : : : : : : : : : :4.3 「カオスのトポロジーが変わる」分岐現象 : : : : : : : : : : : : : : : : : : : : : :4.4 法双曲的不変多様体とその安定・不安定不変多様体の可視化 : : : : : : : : : : : :6.1 (39)-(40) 式の証明 : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :6.2 リー変換 : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : : :6.3 リ─正準変換摂動計算 (LCPT) の応用例 : : : : : : : : : : : : : : : : : : : : : : :6.4 相空間ヤコビ行列の行列要素 : : : : : : : : : : : : : : : : : : : : : : : : : : : : :

67
67
67
67

73
74
80
83

1 はじめに化学反応とは、系が環境（熱浴）から何らかの形でエネルギーを獲得して励起し、ある安定状態（反応系）から別の安定状態（生成系）へ変化するプロセス全般を指し、必ずしも化学結合の解
離生成を包含する必要はなく、相転移なども広い意味で化学反応と云える。化学反応理論は「大2

域的なエネルギー超曲面が与えられているものとした場合に」系が如何に速く、ある安定状態から別の安定状態へ移行するかを問うものであり、1930 年代後半に導入された \遷移状態" 概念は化学反応理論に於いて多大な成果を収めてきた。すなわち、遷移状態とは系が反応する際に反応系から生成系へ至る過程において少なくとも一回通過しなければならない「反応系と生成系とを\"分割する」相空間における 超曲面（これを反応分割面と呼ぶ） を意味し2 、遷移状態概念の導入は複雑な化学反応の速さを「系が如何に反応系から遷移状態へ到達するのか」という大域的性
質（〜遷移状態における状態密度分布）と「如何に遷移状態を通過するのか」という反応ボトル

ネック3 近傍の局所的性質（〜再交差の有無、通過の速度）という２つの視点から論じることを可Eyring ならびに Wigner に始まる遷移状態理論 (1935 年,1938 年)、ミクロカノニカルな遷移状態理論に当たる Rice-Ramsperger-Kassel-Marcus (RRKM) 理論 (1928 年,1952 年)、Keckによる変分的遷移状態理論 (1967 年) および Kramers-Grote-Hynes(KGH) 理論（1940 年,1980 年）能にした。

は、化学反応のダイナミックスの研究において、これまで中心的な役割を果たしてきた。いずれ\も系は与えられたエネルギーの下、「遷移状態」を探し当てるまでに、許される相空間を エル

ゴ−ド的" に経巡る `局所平衡仮説' を基本的な出発点としている4 。また、ボトルネック近傍の局所的な動力学に関しては、遷移状態理論は（その改良版も含めて）一旦、遷移状態を通過した

KGH 理論は再交差は系が反応座標方向に沿って進行することを阻害する \微視レベルでの摩擦（＝分子摩擦）" として働く、という描像に基づいている。一般に KGH 理論は摩擦が無視できる場合において遷移状態理論を包含する。1980 年代後半からの化学反応に関連する実験・理論の飛躍的進展により「ある短い時間滞在
系は再びその遷移状態へ戻ることはない（＝非再交差仮定）、する、反応の「前」と「後」の中間領域にある、遷移状態において化学反応系はどのような動力学的特性を有しているか？」が次第に解明され5 、従来の化学反応理論を抜本的に見直す必要性が問われ始めた。この結果、化学反応の統計性と力学的決定性の「中間」が化学反応ダイナミッ
クスを理解する上で非常に重要であることが次第に分かってきた。
化学反応理論において、「系」とは、通常、化学反応を起こす分子群そのものを指すのではな\ "く、「反応に直接携わる自由度とその自由度と強く結合している（ 熱 として切り離すことが困難な）自由度から構成される要素」を指す。殆どの化学反応理論では、反応自由度（＝反応座標）はポテンシャルエネルギー超曲面の幾何学的情報のみから導かれ、運動量にはあらわに依存しないものと仮定されており、前述の一般化

Langevin 方程式に基づく Kramers-Grote-Hynes 理論もこうした配位空間情報だけを含んだ反応座標を取り扱っている。しかしながら、本稿で解説する
近年の理論・実験研究の成果は「なにを反応に直接携わる自由度、反応自由度、として捉えるべきなのか？」という根幹を真剣に問う必要性を唱えている。

2 慣例的には鞍点でのヘシアン行列 @ 2 V (q) （V (q); qi はポテンシャルエネルギーと座標) の固有値のうち、負の固@qi @qj

有値を伴う固有ベクトルを零とする配位空間上の多次元面ないしは単に鞍点近傍の領域を指す場合が多い。
3 「ボトルネック」も殆んど「遷移状態」と同義に使われている場合が多いが、ここでは反応の「前」と「後」を繋ぐ相空間上の局所領域を指す。
4 唯一、理論の低摩擦極限の表式がエネルギーを供給する媒質と系の相互作用が非常に小さい場合を取り扱っていて、系内部の作用空間（相互作用のない零次では作用は保存する）上での拡散を反応の律速段階と考えている。
5年のノーベル化学賞はこの分野の実験技術に関する功績を讃えている。

Kramers

http://www.nobel.se/chemistry/laureates/1999/index.html 19993

2 化学反応を理解する上で必要なカオスの予備的基礎知識
n 自由度系の古典軌道は t = t0 における正準運動量 pi (t) とその共役な座標 qi (t) に対する初期条件 pi (t0 ); qi (t0 )(i = 1; : : : ; n) が与えられれば、ハミルトンの正準（運動）方程式

により一義的に決定され、この方程式は

@H
dpi
=
dt
@qi

(1)

dqi @H=dt @pi(2)2n 次元の相空間上の決して交わることのない6 軌道の一群を構成する。（自明な）運動の恒量である全エネルギーが存在するとき、任意の保存力学系の2n 次元の相空間上の (2n 1) 次元の超曲面に束縛される。また、エネルギー以外に m 個の運動の恒量が存在する系においては (2n m 1) 次元の超曲面に束縛されることにな
る。n 粒子系の場合、系全体の並進運動量と共役な座標の各々3 つおよび全角運動量の 3 つの計 9
つが全エネルギー以外に自明な運動の恒量であり、n 粒子系の運動は少なくとも (6n 10) 次元相

軌道はすべて空間に束縛されている。加えて、非自明な運動の恒量が存在する場合、系が動き得る相空間の次元は更に減少し、一般に非自明な運動の恒量が存在するか否かは系のもつ全エネルギーの値に依存する7 。

殆どの化学反応理論では、反応座標はポテンシャルエネルギー面の幾何学的情報のみから導か
れ、運動量にはあらわに依存しないものと殆んど無条件に仮定されている。しかしながら、力学の立場から考えるならば、運動量と座標空間から成る

2n 次元の相空間を舞台に議論するほうがむしろ自然である。この章では、化学反応のダイナミックスを理解する上で、どのような場合に運動量空間の情報が必要であり、また無視し得るのかを考察するうえで、必要なカオス力学系の基礎知識8 を学ぶ。

2.1 作用─角変数2

簡単のために、 つの独立な調和振動子系を例に考える。

qi

pm q ; p
i i

i

2

2

2 2

( 2pmi + mi !2i qi )(3)

1 (p2 + !2q2) + 1 (p2 + !2q2)
2 1 11 2 2 22

(4)

H (p; q) =ppm と置き換えると、

X

i=1

i

i

i

H (p; q) =6 力学系の解の一意性から、交差することは「初期条件１つに対して、未来が複数存在する」ことを意味するため、

2n 次元相空間上の軌道は決して交わらない。時間に依存する外力が存在する非自律系の場合には同時刻でなければ2n 次元相空間上の軌道は交差し得る。
7 例えば、T.L. Beck, D.M. Leitner, and R.S. Berry, J. Chem. Phys. 89 1681 (1988) の Ar3 の異性化ダイナミッ

クスの相空間次元の解析など。
8 カオス力学系の諸概念は化学反応を理解する上で非常に有用であり、ここに挙げる基礎知識はそのうちのほんの少しであり、興味のある方は参考文献に挙げた本を読んで新しい化学反応動力学理論の可能性を発掘してほしい。

4

= 1 と置くことができる9 。qi および pi の一般解は
qi (t) = Ai cos(!i t + i ); pi (t) = Ai !i sin(!i t + i )となり、一般性を失わずに mi(5),となり（Ai ; !i および i はそれぞれ i 番目の自由度の振幅、振動数および初期位相）、新たな運動の恒量 Ei

1 p2(t) + !2q2(t) = 1 p2(0) + !2q2(0) = constant
(6)
i i
i i
2 i
2 i
が存在することが分かる (i = 1; 2)。すなわち、運動は 4 次元相空間上の 2 次元10 曲面上に束縛されていることになる。では、その 2 次元曲面を系はどのように運動しているのであろうか？そ
Ei =









の性質は振動数比 !1 =!2 、特に有理数比か無理数比か、に大きく依存している。有理数比の場合、

q1 q2 面に投影した軌道は、Lissajous 図形のような閉曲線を形成している（図 1 左参照）。一方、
無理数比の場合（図 1 右参照）、q1 q2 面に投影した軌道は周期的な閉曲線を成さず、2 つの振幅から定まる長方形のなかをエルゴード的に埋め尽くすことが分かる（各自由度は周期的であるた(quasi-periodic) と呼ぶ）。
時間に対して不変な運動の恒量が 2 つ存在するということは p1 ; q1 ; p2 ; q2 は互いに独立ではな
く、2 つの独立変数でそれらを記述することができることを意味し、(p; q) 座標系は \冗長" であり、ハミルトニアンをもっと簡単に書き下すことができる。\（時間に対して不変な）運動の恒量
め、準周期的が存在する場合、その運動の恒量を新しい運動量とする正準座標系から運動を眺めると自由な回"転運動として焼き直すことができ、相空間構造に対するはっきりした描像を得ることができる。
以下にそれを説明する：

Poisson 括弧 fg：@u @v
fu; vg  ( @@

ちなみに、正準変換とは、

X

i

i

i

@v @u
);
@i @i(7)（u; v は任意の微分可能な関数であり、i および i は任意の正準な運動量とそれに共役な座標であPoisson 括弧の関係式
fqi; pj g = Æij ; fqi; qj g = fpi; pj g = 0(8)（Æij は Kronecker のデルタ）を保存する座標変換 (p; q) ! (p; q) を指す。Poisson 括弧の関係式がる）に対して、満たされていればその変数の組はから

dpi@H (p; q) dqi @H (p; q)
=
;
= @p
dt
@qi
dt
i

(9)

dpi
@ H (p; q) dqi @ H (p; q)=;= @ p
dt
@ qi
dt
i

(10)

のように、ハミルトンの正準方程式の形が保存し、どちらの正準座標系で運動を記述しても等価である。また、これらの座標系は

pi = pi (p; q); qi = qi (p; q)

9 本稿では以下すべて mi = 1 で記述する
10 E = P2 Ei なので、E 一定のもとで E1 ; E2 は独立ではない。恒量の数は 2 つ。
i=1

5

(11)

pi = pi (p; q); qi = qi (p ; q)(12)のように互換性をもち、正準変換は、例えば、

@F
;
@qi
@F
pi =
;
@ qi

(13)

pi =

H (p ; q; t) = H (p; q; t) +(14)@
F (q; q; t);
@t

(15)

を満たす関数（母関数）F などによりおこなうことができる11 。
0.3

0.25
0.20

0.2

0.15
0.10

0.1

q2

q2

0.05
0

0

-0.05
-0.1

-0.10
-0.15

-0.2

-0.20
-0.3

-0.25
-0.3

図

1: q1

-0.2

-0.1

0

q1

0.1

0.2

-0.3

0.3

-0.2

-0.1

0

q1

0.1

0.2

0.3

q2 面に投影した振動数比 !1 =!2 = 2 の軌道（左）、振動数比 !1 =!2 =  の軌道（右）2 

さて、 つの運動の恒量 1 ; 2 を運動量 p1 ; p2 とする正準座標系が存在するものとしよう。 1 ; 2
は時間に対して不変であるから、新しいハミルトニアン H は「共役な座標 q1 ; q2 に依存しない」

だけの関数でなければならない。

0 = ddti = ddtpi =
また、

2

 

@ H (p )@ qi

dqi @ H (p)
= @ p  !i(p) = !i( ) = constant;dtiであり、上式の第 項めは時間に不変な運動量

p

(16)
(17)

p の関数であることから !i は一定となり、新しい座標系では運動は自由運動に焼き直されることがわかる。

pi = constant; qi = !i t + i
この（ただし 1 は初期位相）(18)(p; q) が正準座標系であるためには Poisson 括弧関係式が満足されていることが示されればよい。

11 詳細は解析力学の本を参照のこと

6

2 次元の独立した調和振動子系をもうすこし考察してみよう。（6）式の不変量 Ei は実は作用変数 Ji と呼ばれる量に比例する。作用変数 Ji は

Ji 

1
2

I

(19)

pi dqi25と定義される（ここで積分は一周期についてとる）。 次元調和振動子系の例に対して（ ）式を代入すると、Ji はJi =A2i !i

Ei
=
2 !i

(20)

Ei =

(21)

となり、この場合の作用変数は時間に対して不変であり、ハミルトニアン H は作用変数にしか依存しないことが分かる

:

H=



X

i

X

!iJi :

作用変数に共役な正準変数を i とおくと、

5di @H (J)= @J = !i = 一定
dti
i(t) = !it + i



2 (

(22)
(23)

となって、（ ）式から明らかなように、変数 i は 2! 毎に相同な運動を与え、 n n は任意の整数 の不定性をもつ性質を有することから角変数と呼ばれる。

)

作用─角変数をi(p; q) で書くと、1 p2 + !2q2 ;  (p; q) = tan 1( !iqi )Ji (p; q) =2! i i i i
p


i



i

Poisson 括弧関係式@Jj @ i) = Æij ; fi; j g = fJi; Jj g = 0
@q @p

(24)

となり、作用─角変数が正準であることは@Jjfi; Jj g  ( @@qi @p
k
k
k
X

k

k

(25)

を確かに満たしていることから確認することができる。
調和振動子系に限らず、自由度の数 N と等しい数の独立な運動の恒量が存在する場合、ハミル
トニアンは作用だけの関数として表現することができる。

H (p; q) = H (J);di @H (J)= @J = !i(J)
dt
i

(26)
(27)

また、!i は運動不変量である作用の関数であるからその値も時間に対して不変である。換言すると、準周期的な規則運動を示すあらゆる力学系の「正体」は自由な回転子のあつまりとみなすことができる。このような系を可積分系と呼ぶ。化学で取り扱う系のほとんどすべては可積分系ではないが、基準モードで系を議論することは暗黙理に可積分であることを想定している。また、\"相空間上のある特定の領域および特定の自由度空間において 可積分的 であることは大いにあり得る12 。

12 例えば、後述するように、非線形性が強く殆んどすべての軌道はカオス的であるエネルギー下においても、ポテ
ンシャル曲面上のサドル領域の運動がそれに相当する。

KAM トーラスもしかりであるがここでは立ち入らない。

7

2 次元調和振動子系は相空間上でどのような構造をとっているであろうか？図 2 に、
!1 = 1、!2 = 2 の振動数が有理数比の軌道を q1 p1 面および q2 p2 面に投影したものを示す。
Stokes の定理13 から作用 Ji は
1 p dq =Ji =dpidqi(28)
2 i iH (p;q)Eであることから、図 2 の閉曲線で囲まれた面積を指していることが分かる。エネルギー一定のもとでは、初期値 fp1 (0); p2 (0); q1 (0); q2 (0)g が与えられた瞬間にある決まった一定値をとる。作用─角座標系で 4 次元相空間上に束縛された 2 次元曲面を眺めると、
トーラス

I

Z Z

H (p; q) = H (J1 ; J2 ) = E(29)

J1 ; J2 = 一定

(30)
(31)

1 = !1 t + 1 ; 2 = !2t + 2から、図 3 のように、半径 J1 および J2 で角座標 1 および 2 のちょうど \ドーナツ" のような形を思い浮かべることができる。振動数比 !1 =!2 が 2 の場合、軌道が半径 J2 の円を 1 周回る間にちょうど半径 J1 の円を 2 周回り、ドーナツの表面を閉曲線を形成する。振動数比が無理数の場合、軌道は閉曲することなく、「ドーナツの表面上を」隈無くエルゴード的に動き廻る。

図2: !1 = 1; !2 = 2 の軌道を q1

p1 面および q2

p2 面に投影したもの

n 次元可積分系（（26）-（27）式）の場合は振動数 !i が一般に作用 J の関数であり、半径 Jj などに依存して振動数が変化する。つまり、同じエネルギーであっても初期条件に依存して、あるP
トーラス（＝ n 次元のドーナツ）は閉曲軌道によって覆われ（ ni=1 ki !i(J) = 0（ki はゼロでない任意の整数））、別のトーラスは準周期的な軌道で隈無く埋め尽くされる（任意の振動数対は無理数比）ことが想像できるであろう。!i(J) は J に対する連続関数であり、有理数比を与える実数の組は実数軸上に無数（可算無限個）に存在するが、無理数比を与える組のほうが「より」無数（非可算無限個）に存在するため、一般には準周期的な軌道で埋め尽くされるトーラスが支配的に分布する。

13 R R

S



@p
@y

@q
@x



dxdy = C (pdx + qdy )
H

8

図

3: 2 次元トーラスと !1=!2 = 2 の軌道

2.2 正準変換摂動理論

2.1 章で可積分系の代表例である調和振動子系において、観測する視点を変えることによって、\見通しの良い" ダイナミックス描像を得ることができる事を見た。しかしながら、我々が取り扱う系は一般に非線形な非可積分系である。本章では可積分系を零次とする（自律系／非自律系を問わない）任意の非線型ハミルトニアンに対して、できる限り作用が保存するように非線型摂動項を取り込みながら正準変換を実行し新しい座標系へ変換する正準変換摂動理論（

canonical

perturbation theory, CPT）を紹介する。実際の手順に入るまえに、CPT が与える新しい座標描
像はどのようなものなのかを説明する。

2.2.1 相空間のなかの一般化された基準振動概念
正準変換摂動理論は問題とする m 次元ハミルトニアン H

J のみで記述されることを仮定する。
n Hn (p; q) = H0 (J) + n Hn (J; );(32)

展開可能で、かつ零次ハミルトニアン H0 は作用変数

H (p; q) =X

n=0(p; q) が摂動パラメータ  で巾級数で

X

n=1



正準変換 W は新しいハミルトニアン H の角変数依存性を最小にする（すなわち、新しい作用変数J が出来る限り保存する）ように新しい正準座標系を決定する。
H (p; q) W! H (p; q) = H (J) = n H n (J):
X

n=0

9

(33)



仮に H が、少なくとも、ある有限の摂動次数の範囲において角変数 に依存しない場合、k 番目のモードに対する新しい作用─角変数の運動方程式は次のように表すことができる。すなわち、@ H (J)dJk _=Jk = k = 0;
dt
@

(34)

Jk = constant;

(35)


_ k = @ H@ J(J)  !k (J) = constant;kk = ! k (J)t + k :

および

(36)
(37)
ここで、 k は k 番目のモードに対する初期位相である。仮りに無摂動ハミルトニアン H0 (p; q)
が n 次元調和振動子系であれば、

1 (p2 + !2q2) = M ! J
(38)
k k
k
k k
k=1
k 2
 に従う新しい座標 q(p; q) な
（ !k は k 番目のモードに基準振動数である）、(34)-(37) 式より、H
らびに共役な運動量 p(p; q) に対する運動方程式を容易に導くことができる。
d2 qk (p; q)+ ! 2q (p; q) = 0;(39)
H0 (p; q) =

M
X

X

dt2

k k

!k dqk (p; q):! k dt(40)

ei! (J)t + e i! (J)t!k ei! (J)t!k e i! (J)t

(41)
(42)

pk (p; q) =! k (= ! k (J) = ! k (p; q)) は J が一定値であるため時間に依存しない（6.1 章参照）。一般解はqk (p; q) =pk (p; q) =
（

および





k

k

k

k

(41)-(42) 式は、\ハは qk および pk の初期条件によって決まる任意の定数である）。

"

ミルトニアン H が存在する場合 、オリジナルの座標系において複雑に振る舞う（ように見える）
運動は、皆、エネルギーが同じであっても、初期条件によって振動数が変化する相空間上の分離された周期軌道として解釈できることを意味している。しかしながら、現実の化学反応系において、式(39) と (40) を満たすハミルトニアン \H " はまず存在しない。その理由は次々節で述べるが、それにも拘らず、相空間上に一般化された基準座標概念は統計性と力学的決定性の中間が本来問われるべき化学反応において非常に重要な役割を果たすことを後章で明らかにする。14

14 Goldstein などの代表的な力学の教科書は基準座標についてはかなり後の章で触れているにすぎないし、基準座標の概念自体は力学の立場からは自明なことで興味の対象にはならない。しかしながら、一般に強い非線形力学系である分子の大半が室温下で（その分子を構成する）原子団特有のグループ振動数を保有することは分子個性に由来す
る力学構造に依拠するものでそれほど自明なことではない。相空間上に一般化された基準座標・振動数概念はこれまでの基準座標概念には含まれていなかった非線形性およびダイナミックス情報を含んでおり、従来の基準振動解析を越えた新しい分子理論に発展する可能性を大いに秘めている（と思う）。

10

2.2.2 リ─正準変換摂動理論
LCPT）は、古い正準
変数 (p; q) と新しい正準変数 (p; q) の両変数から成る母関数 [(13)-(15) 式] に基づく Poincare-Von
Zeipel の手法を凌駕する画期的な方法を構築した。この方法は、新旧の正準変数から成る \やっかいな" な関数反転操作が一切なく、すべての項は単純な Poisson 括弧を繰り返すことによって書き下すことができる。\ハミルトニアン（リー母関数）" W に従う時間  に沿った相変数15 z の仮
天文物理学者である堀源一郎により開発されたリ─正準変換摂動理論（

dz= fz; W (z; )g  LW zd

想的な時間発展の、Poisson 括弧 fg の繰り込みによる形式解z() = exp

Z LW (0 )d0 z(0)  T z(0)

(43)
(44)

を（通常）リ─変換と呼ぶ。
あらゆるすべてのハミルトン系の時間発展がある時刻から別の時刻への、ハミルトンの正準方(44) の関数形で表させる任意の変換
は z(0) が正準変数ならば、正準変換であることが容易に証明できる（その逆も）（6.2 章の「リー変換」参照のこと）。z(0) の任意の関数 f に対するリー変換、すなわち T f 、は z(0) および  の関数として表される複雑な非線形関数 g を与えるが、その関数「値」は f (z()) に等しい。すな
程式の「形」を保存する「正準変換」と見做せるように、式

わち、

f (z()) = T f (z(0)) = g (z(0); ) :(45)

また、逆リー変換 T 1 を上式に施すとT 1 f (z()) = f (z(0)) = h (z(); )(46)(z())、は z() および  の関数として表される複雑な非線形関
数 h を与えるが、その関数「値」は f (z(0)) に等しい（6.2 章の付録「リー変換」を参照のこと）。
となり、逆リー変換、すなわち T 1 f

さて、リー変換を用いて次の正準変換を考えよう16 。

H (z()) = H (z(0))(47)

すなわち、リー変換を用いた新旧ハミルトニアンのお互いの関係は、W に従う仮想的な時間発展"0 の点 z(0) で評価された「旧い」ハミルトニアン H と \時刻" 経過した新しい は等しいと解釈される。さて、（46）式と（47）
点 z() で評価された「新しい」ハミルトニアン H
 とおくと、式を比べると、f ! H および h ! H
T 1 H (z()) = H (z(0)) = H (z())(48)\に沿って、 時刻

15 正準座標およびその共役な運動量をまとめて表記する。
16 ここでは、簡単のために自律系を例に説明するが、非自律系の場合は全エネルギーと時間を相変数に加える \拡"

張された相空間 を用意することでまったく同じ枠組みのなかで取り扱うことができる。

11

z (p; q)、時刻 0 での z を (p; q) と表すと、（48）式はH (p; q) = T 1 H (p; q) = H (p; q)(49)が得られる。例えば、 時刻  での をとなり、自律系の正準変換における新旧ハミルトニアンを関係づける良く知られる式H (p; q) = H (p; q)に対応することが分かる。

(50)

,

さて、新旧のハミルトニアン H H ならびにリー母関数 W が  の巾級数で展開できるとしよう。
すなわち、

H

=

H

=

W

=

LW

=

T

=

T 1

=

ここで、

n Hn ;

(51)

n H n ;(52)

n Wn+1 ;(53)

n Ln+1 ;(54)

n Tn ;

(55)

n Tn 1 ;(56)

Ln = fWn ; g

(57)

X

n=0
X

n=0
X

n=0
X

n=0
X

n=0
X

n=0

, T および T 1 も巾で展開できる。
最初に Tn 、Tn 1 および Ln に関する便利な関係式を導出しておく。(54) および (55) 式を (318) 式に、(54) および (56) 式を (319) 式にそれぞれ代入し、 の各べきの係数を等しくおけば、Tn および Tn 1 (n > 0) に対する漸化式
n 1
1
T =
T L ;
(58)
であり、W が  の巾で展開できるとしたので、自動的に LW

X

n

T 1
n

n m=0 m n m=

1 n 1LX

n m=0を得ることができる。ここで T0 および T0 1 は n m Tm1(59)= 0 における T および T 1 であるから、T0 =T0 1 = 1 である。例えば、 の 3 次まで導出すると

=
T2 =
T3 =
T1

L1 ;

1 L + 1 L2 ;2 2 2 1
1 L + 1 L L + 1 L L 1 L3 ;3 3 6 2 1 3 1 2 6 1
12

(60)
(61)
(62)

= L1 ;
1
1
T2 1 = L2 + L21 ;
2
2
1
1
1
1
T3 1 = L3 + L1 L2 + L2 L1 + L31
3
6
3
6

T1 1

となる（T および L は一般に可換ではないことに注意。Li Lj

(63)
(64)
(65)

6= Lj Li）。

H 、T 1 および H の  による巾展開の式、（52）、（56）および（51）式を H = T 1 H の両辺に
 n を関係付けることがで代入し  の各べきの係数を等しくおき、（59）式を用いると Hn および H

2きる。例えば、 次まで求めると

0 : H 0
1 : H 1= H0 ;= T0 1H1 + T1 1H0 = H1 + L1H0 ;1;= H1 + dWd2 : H 2 = T0 1 H2 + T1 1 H1 + T2 1 H0 ;= H2 + 21 L1 (H 1 + H1) + 12 L2H0 ;2= H2 + 21 fW1 ; H 1 + H1g + 21 dWdとなる。ここで、

Ln H 0 = Ln H0 = fWn; H0 g 

dWn
d

(66)
(67)
(68)
(69)
(70)
(71)
(72)

を用いた。 は非摂動ハミルトニアン H0 に従う時間である。

51）式以後の式では敢えて正準変数の引数（すなわち、(p; q) ないしは (p; q)）の表記を省略した。Poisson 括弧式を用いた計算は正準変数の組の種類によらずに不変であるため、リー変換（

正準変換摂動論においては、正準変数の組の種類は意味がなく17 引数程度の意味しか持たない。

(66)-(71) 式における引数を (p; q) とおこう。(66) 式から、ただちに O(0 ) では
H 0 (p; q) は H0 の引数を (p; q) から (p; q) に換えたものにすぎないことがわかる。一方、0 次以上
 i と Wi が順次常の各べきのオーダー O (n ) (n  1) では、１つの方程式に対し 2 つの未知関数 H
 1 と W1、O(2 ) における H 2 と W2 である。すなわち、１つに現れる。例えば、O (1 ) における Hの条件から 2 つの未知関数を決定しなければならず、そのためにはなんらかの原理が必要となる。
 (p; q) が新しい角変数  を含まない
平均化原理 そこで、変換される新しいハミルトニアン H
ここで、便宜上という要請を満たす18 ように決定することにしよう。もしこれが可能であれば、

dJk@ H (J)=k = 0dt@(73)

 (p; q) が産み出されることを意味する。これは運動をトーラス上の自より、新しい運動の恒量 Jk

由運動に焼き直すことに対応し、微分方程式の階数の低下をもたらすダイナミックスに対して見通しの良い「視座」を創出したことになる。

17 (p; q)、(p; q)、(J; )、(J;  ) のどれを使っているかに依らずに成り立つ。
18 少なくとも「ある有限次」の摂動計算に渡って。

13

問題とするハミルトニアンを、H0 に対する作用変数

J のみで表される非摂動ハミルトニアン

H0 と作用変数と角変数の両方に依存する摂動項 Hn (n  1) の和であるとする。すなわち、
H=

X

n=0

n Hn = H0 (J) +

さらに、各オーダー n において Hn が新しい角変数

Xn=1

n Hn (J; );(74)

 に依らないように、もうひとつの未知関数

Wn を決めよう。まず、 O(1 ) において H 1 と W1 は以下のようにして決定する。
H 1

1;= H1 + dWd1;= hH1i + fH1g + dW
d

(75)
(76)

ここで

1 H (p( 0); q( 0 ))d 0;hH1i = lim
!1  0 1
fH1g = H1 hH1i
Z 

平均化された部分
振動する部分k d@H0 (J )=f= constant  !k (J)k ; H0 g =d@ Jk k = !k (J) + k ( k : 初期位相)

である。

(77)
(78)
(79)
(80)

 と線形関係にある。それゆえ H1 の（ に対して）平 に依存しない部分に対応し、振動部分 fH1g は  に依存する部分に対均化された部分 hH1 i は  i が  に依存しないように決定したいならば、（76）式の第 2 項と第 3 項が打ち消しあ応する。Hから、H0 に従う時間  は新しい角変数

うように W1 を決めてやればよい。

H 1
W1

= hH1i;=fH1gd:
Z

(81)
(82)



2

O() における付加条件であり、元々W1は任意であったのでこのように取り扱っても一般性は失われない。同様にして、 O (2 ) において
これが つの未知関数 W1 、H1 を決定するにあたっての

も、W2 を以下のように決めれば H2 も に依存しないように決定することができる。

= hH2 + 21 fW1 ; H 1 + H1gi;
W2 =f2H2 + fW1; H1 + H 1ggd:
H 2

Z

(83)
(84)

ある有限次 i まで摂動展開が可能であるとしよう。そのとき、その有限次で打ち切った新しい J のみの関数であるから

ハミルトニアン H は

H (J) =i
X
k=0

k H k (J) =

14

i
X
k=0

k Hk (J; )(85)  
dH 0
0 = ddH = fH;
H0 g = fH 0 ; H g =(86)dt
 に従う運動において H 0 は運動の恒量である。
より、H\ハミルトニアン H を  のすべての次数に渡って時間  に依存しない形式に変換できるならば"Hは H と厳密に等しいことになる。

H (J; ) = H (J)(=

1

X

i

i H i (J)):

(87)

dH 0
(88)
dt
 0) が存在していることがより、元々の H に従う運動においても「非自明な」運動の恒量 H0 (= H
このとき

0 = dH= fH; H 0g = fH 0; H g =d

証明される。
リー変換に基づく関数の変換 リー変換に基づく正準変換摂動理論のもつ最大の利点のひとつ

 (p; q) のみならず、任意の物理量 f(p; q) が時間発展演算子 T を用いると容易に元の座標系 (p; q) で書き下せることである。付録 2 に示されるように、任意の微分可能な (p; q) の関数 f 、f (p; q)、に対して、元の (p; q) および  で表現される「関数値が同じ」新しい関数 f を見出すことができる。

は、一旦、W が各次数を通して求まれば、Hf (p; q) = T f (p; q) = f(p; q; ):(58) 式より  の各オーダー毎、関数 fn(p; q)

それゆえ、少なくともある有限次まで Wn が求まれば、

(p; q; ) を解析的に計算することができる。
f(p; q; ) = n fn (p; q) = n Tn f (p; q):に対する新しい関数 f

例えば、

X

X

n=0

n=0

(89)

O(2 ) まで求めると、
0 : f0 (p; q) = f (p; q)
1 : f1 (p; q) = L1 (p; q)f (p; q) =

(90)
(91)
(92)
(93)

fW1 (p; q); f (p; q)g12
2 : f2 (p; q) =2 L2 + L1 f
= 21 (fW1 ; fW1; f gg fW2; f g)となる。ちなみに、この場合は時間発展演算子 T の引数として \p" と \q" を用いたが、例えば(103) 式においては \p" と \q" を用いた。同様に (293)、(301) 式におけるリー母関数 W の引数は各々(p; q) および (p; q) を便宜上用いて説明した。このように、リー変換において引数自体は本

来意味がないことに注意しよう19 。

19 筆者自身もあちらこちらの論文で（名前自体は実体がないので）一貫性なく色々な引数を使っている。

15 (p; q)、振動数 ! k (p; q)、運動量 pk 、座標 qk などを元の座標系 (p; q) で例えば、新しい作用 Jk
表すためには

p2 + !k2qk2
Jk (p; q) = T Jk (p; q) = T k
2!k ;@ H (J)! k (p; q) = T;@Jkpk (p; q) = T pk ;qk (p; q) = T qk :!

を計算すれば良い。ここで、簡単のため、関数 f に対しリー変換を施した

(94)
(95)
(96)
(97)
(p; q) および  の関数

(p; q) と表記する。変換する関数 f として H0 の作用変数 Jk (p; q)
p2 + !k2qk2 1
Jk (p; q) = k(98)
2!k = 2 E=H (p;q) pk dqk
を用いる際は変換された関数を Jk (p; q) などで表す。実際（6.3 章の応用例参照）、pk (p; q) および qk (p; q) の i 次までの級数和 pikth (p; q) および qkith (p; q) は

をf

I

0

i
X

pikth (p; q) =

n=0
i
X

qkith (p; q) =

n=0

n

X

n

X

j

cnj p2s +1qt ;(99)

c0nj p2s qt ;

(100)

nj

0

nj

j

nj

0

nj

のようにそれぞれ書き下すことができる。ここで、
M

M
X

p2s +1  Y psl
l=1
nj

qt

nj



M
Y
l=1

njl

l=1

qltnjlおよび、各係数は元のハミルトニアンおよび!

snjl = j2snj + 1j ;

M
X
l=1

tnjl = jtnj j

!

(101)
(102)

CPT の次数 i に依存する。例えば、cnj および c0nj は

 (p; q) および qkith(p; q) における n 次の j 番目の項に掛る実数係数を、snjl および tnjl
は pikth (p; q) の n 次成分の j 番目の項の pl 、ql に掛る任意の非負整数 ( 0) である。ここで大事なことは一般に i  1 においては pikth (p; q) および qkith (p; q) は確かに時間反転対称性を保持していること、pikth (p; q) および qkith (p; q) における旧変数 pk および qk の寄与は必ずしも大きくなく、殆それぞれ pikth

んどすべてのモードが寄与していることである20 。

 (p; q) および qk0th (p; q) は元の座標系 pk および qk と一致し、同様に Jk0th (p; q) および! k0th (p; q) は H0 に対する Jk (p; q) および !k と一致する。換言すれば、摂動展開を 0 次で打ち切った物理量は皆非摂動ハミルトニアン H0 に由来する物理量であり、(p0th; q0th) は H0 = E におけ
ただし p0thk

る（準）周期軌道（＝基準座標）に相当する。一方、（少なくとも）ある有限次 i まで摂動展開が20 6 原 子ア ル ゴン ク ラ スタ ー の鞍 点 近傍 に おけ る p1 (p; q) お よび q1 (p; q) を 2 次 まで 求 めた と きの 表 式

105 (1999) 10838 の on-line supplement に公表されているので参照されたい。

J.Chem.Phys.

16 (J) は O((i+1) ) の誤差以内でH に等しく21 、(pith ; q ith ) は H = E における相空間上の（準）周期軌道を表し、それは (p; q) の可能であるとき、その有限次で打ち切った新しいハミルトニアン H関数として表すことができる。さらに無限次まで摂動展開が可能であるとき、新しいハミルトニ

アン H は元のハミルトニアン H と厳密に等しく、(p; q) は H (= H ) = E における（隠れた）相空間上の（準）周期軌道を表している。

Deprit の漸化式 Deprit により最初導出された便利な漸化式を紹介する。以下の正準変数を一切使わない導出は

Cary22 の総説によるものである。

 = T 1H に対して、 T を前から掛けて  に対する微分をとる。すなわち、
dT dH dH
H +T= d :
d
d

(103)

d

d

(104)

mTn 1m Hm

(105)

まず、 H

(318) 式を使って dT=d を取って、 T 1 を前から掛けると、dH
= T 1 dH
LW H +
を得ることができる。LW 、H 、T 1 および H の  についての巾展開の式を両辺に代入して、n(n > 0) のべきの係数を等しくおくと、
nX1
m=0

となる。

Ln m H m + nH n =

n
X
m=1

72 式より最初の和の最初の項を
dW
Ln H 0 = Ln H0 = fWn ; H0 g = n

(106)

d

2

と書き換え、 番目の和の最後の項を別途取り出しな漸化式

H n = Hn +1 n 1 (LX

Pの部分をまとめると、n > + mTn 1m Hm) + 1 dWn :n m Hmn m=10 に対する有用
(107)

n d

23 。
が導出される（非自律系では n1 @W@t が右辺に追加される ）nDragt と Finn の手法 Dragt と Finn は、ハミルトニアンが最低次を 2 次とする (p; q) の巾級


数で展開できるクラスに対して、高次計算になるほど H を算出するのに要する項の数が堀の手法よりも顕著に少なくなる方法を見出した。彼らは時間発展演算子 T をT () = e L = e (L +L + L ) ! T 0 () = e L e  L e  L   
1

2

2

3

1

T0 = 1 +

1

X

n=1

2

2

3

-Deprit
(108)

3

(109)

n Tn021 2.2.3 章「小さな分母の問題」で明らかになるように、形式上 O((i+1) ) であるが、現実系においてはこの \誤差"は一般に発散し展開級数は収束しない。

22 J.R. Cary, Physics Reports 79 130-160(1981).
23 A. J. Lichtenberg and M. A. Lieberman,

Verlag, New York, 1992).

Regular17

andChaotic

Dynamics

2ndEd.

(Springer-ここで   Lmn
(110)n
m1 !m2 !    mn !
m +2m ++nm =n（和は m1 +2m2 +    + nmn = n を満たすすべての正の整数 (m1 ; m2 ;    ; mn ) の組に対してとる）のように予め繰り込みリー変換を逐次行う形式に再定義した。Li (= fWi ; g) は (i + 2) 次の同次多
T0 =

m1 m2
m1 +m2 ++mn L1 L2

( 1)

X

1

2

n

n

( ) 1は

項式 Wi から成る。逆変換 T 0(T 0) 1 () =    e L e L eL3(T 0 ) 1 = 1 +(T 0)n 1 =

n=1

2

2

(111)

1

n (T 0 )n 1

(112)

Lmn    Lm2 Lm1
m +2m ++nm =n m1 !m2 !    mn !
X

1

3

1

X

3

n

2

2

1

(113)

n

で表される。例えば、 次まで求めると= L1 ;
1
T20 = L2 + L21 ;2
1 L3 ;
T30 = L3 + L1 L2
1
6(T 0)1 1 = L1 ;(T 0)2 1 = L2 + 21 L21;
(T 0)3 1 = L3 + L2 L1 + 61 L31
 (= n nH n) = (T 0) 1H はとなる。新しいハミルトニアン H
T10

(114)
(115)
(116)
(117)
(118)
(119)

P

H n =
で表される。（

nX
m=0(Tn0 m) 1Hm(120)58）式と Tn0 の n における項の数を比較すると、表 1 のように、5 次以降の項の数
表

1: 各 n 次における Tn と Tn0 の項数n

Tn Tn0

1 1 1
2 2 2
3 4 3
4 8 5
5 16 7
10 512 4218

が顕著に異なってくることがわかる。これは、たとえば、正準変換

T 1 ()f (z()) = eL f (z() = e(L +L + L ) f (z()) = h(z(); )
1

2

2

(121)

3

を一度に行うのではなく、(T 0) 1 ()f (z()) =    e L e L eL f (z())=    e L e L h(1) (z(); ) =    e L e(L ) h(1) (z(); )=    e L h(2) (z(); ) =    e( L ) h(2) (z(); )=    h(3) (z(); ) ! h(1) (z(); )
3
3

3

3

3

2

2

2

h

i

1

3

2

h

3

2

3

h

2

3

i

i

(122)
(123)
(124)
(125)

-Deprit の方法で指数部をべき級数展開する際、LiLj 6= Lj Li の非可換性 (i 6= j ) から生じた項数の増大を抑えられたためである（ここで h(i) は \ハミルトニアン
(i 1) Wi " によるリー正準変換により変換されたことを表している）。しかしながら、分子科学の興味の対象である、例えば、アルゴン 6 原子クラスター（12 内部自由度）、の少数多体系では、2と逐次正準変換するため、堀〜3 のオーダーを越えると、評価すべき項の数は膨大に膨れ上がってしまうため、（少なくとも5

既存の手法では） 次以上への摂動計算は困難のようである24 。

Birko -Gustavson の標準形 (Normal Form) Birko と Gustavson は、2 次を最低次とする


斉次多項式の級数展開で表されるハミルトン系は「H のみを求めるならば」新旧の両正準変数か(p; q) を用いても容易に高次まで導びくことができることを証明した。変換され は、調和振動子の作用変数の巾級数展開の形た新しいハミルトニアン Hら成る母関数 F

H (p; q) !

X

i

!i Ji +Xi;jcij Ji Jj +

X

i;j;k

dijk Ji Jj Jk + : : :(126)

Birko -Gustavson の標準形、単に Birko の標準形と呼ばれている。原著論文
では新旧の両正準変数から成る母関数を用いて証明しているが、ここでは Dragt-Finn の論文25に習って逐次リ─変換を用いて説明し、2 次を最低次とする斉次多項式の級数展開で表されるハをしており、今日、ミルトン系の性質を考えよう。
まず、任意の停留点近傍でハミルトニアンが巾級数で形式的に展開できると仮定する。

H (p; q) =

1

X

k=2

(k 2) Hk (p; q)(127)

Hk (p; q) は p および q に関する k 次の斉次多項式である。勾配は 0 であるから最低次は 2 次となる（ただし H0 (0) = 0）ので、それを  の 0 次とした。例えば、
H (p; q) = 0 H2 (p; q) + 1 H3 (p; q) + 2 H4 (p; q)   (128)
1
Cjklmqj qk ql qm + : : : : (129)= 0 2 (p2j + !j2qj2) + 1 Cjklqj qk ql + 2
jj;k;l;m
j;k;l
X

X

X

24 近年、Caltech から Bristol 大学へ移った応用数学者の Stephen Wiggins は 160CPU の並列計算機を構築して、ハードウェアの観点から大自由度系のリー変換を用いた

CPT を敢行する試みをしている。

25 A.J. Dragt and J.M. Finn, J. Math. Phys. 20, 2649-2660(1979)19

p(! ) ! pi; (!i)qi ! qi と置きなおすと、H2 はq

ここで、 pii

H2 (p; q) =更に、zi

= pi + iqi ; zi = pi

1 N ! p2 + q2 ;2 i=1 i i i


X

(130)



iqi のように複素変数を導入すると

1 N ! z z2 i=1 i i i) 表示でと表すことができる。また、Hk (k > 2) は一般に (z; zm nHk (z; z) =(m)(n) zi zi H2 (z; z) =

X

P

i

Y

(m +n )=k
i

(131)

X

i

X

i

(m)(n)ii(m)(n) z z

m n

(132)

( + ni ) = k（mi および ni は非負の任意の整数）を満たす項に対して取る。 (m)(n) は複素定数である。複素変数表示の場合、Poisson 括弧式はと表される。ここで、和は

P

i mi

@u @v @v @u
fu; vg = 2i ( @z)@ zi @zi
i @ zi
i
X

(133)

となる。

 = (T 0) 1H を使って新しいハミルトニアン H を求めてみよう。ここでは、Dragt-Finn の原著論文の表記法に従って書くことにする。
（122）-（125）と同様の操作をハ
さて、逐次リー変換 H

ミルトニアン H に対して行うとH (k) = eF H (k 1)

(134)

k

 = H とする）
(T 0) 1H =    eF eF eF H =    eF eF H (3) =    eF H (4) =    H (5) = H (1) (= H ): (135)とおくと（ただし H (2)
5

4

h

3

i

5

h

4

i

h

5

i

ここで、Fk は k 2 Lk 2 に相当し、k 2 Wk 2 を表す k 次の斉次多項式 fk を用いるとFk  ffk ; gと書くことができる26 。Wk 2 が k 次の斉次多項式である理由は（単一）リー変換の（
び（

71）式

(136)
68）およ

1 : H 1= H1 + fW1 ; H0g
112 : H 2 = H2 + fW1 ; H 1 + H1 g + fW2 ; H0 g
2
2
からも理解することができよう。一般に、Poisson 括弧式 ffk ; gj g（fk および gj は任意の k およびj 次の斉次多項式）は正準座標および運動量に関する 1 階偏微分をそれぞれ行い掛け合わせるので (j + k 2)( 1) 次の斉次多項式を産み出す。
26 原著論文では Fk  f ; f 0 g を用いているが、Poisson 括弧式の性質から fk = f 0 であり議論の内容になんら影
k
k

響しない。

20

H1 が 3 次の斉次多項式であるから、fW1 ; H0 g は（H1 の振動部分を取り除くためには）3 次の 1）も 3 次の斉次多項式でなければならない。
斉次多項式でなければならないので、W1（および H
 1 + H1 g は 4 次の斉次多項式となるから H2 + 12 fW1 ; H 1 + H1g の振動部分を同様に、H2 、fW1 ; H 2）も 4 次の斉次多項式でなければならない。
取り除くためには W2 （および H
 も H 同様、斉次多項式で級数展開できるとしよう。（134）式は
さて、新しいハミルトニアン H
1

X

j =2

H j(k) = eF

1

k

X

j =2

H j(k 1) = 1 +

1 1

X

i=1

(F )
i! k

!

i

2
H 2(k) = H 2(k 1) =    = H 2(0) = H2

1

X

j =2

H j(k 1)(137)と書くことができる。ここで、例えば 次の係数を等しくおくと、

(138)

( )
( ( 2) + j ) 次の斉次多項式を産み出し、k > 2 では i = 0（すなわち恒等変換）を除いて j よりも大きい次数の斉次多項となる。なぜならば Fk i gj （gj は任意の j 次の斉次多項式）は一般に i k式を与えるためである。一般に j 次の係数を等しくおく（すなわち、j 2 の各係数を等しくおく
ことに対応している）と

j < k : H j(k)j = k : H k(k)= H j(k 1)(139)= H k(k 1) + Fk H2
(140)
2
i
(F )
(F )
j > k : H j(k) = H j(k 1) + Fk H j(k (k1) 2) + k H j(k 2(1)k 2) +    + k H j(k i(1)k 2) +    : (141)2!
i!
i
ここで、恒等変換以外の変換 (Fk ) (i = 1; : : : ; 1) からの最低次数は k であるので、k 未満の次数は恒等変換の項しか残らない。また、k 以上の次数における Fk のべき数 i の上限は i  kj 22 を満 に含まれたす最小整数である。（135）、（139）および（140）式から、新しいハミルトニアン H
 k (= H k(1)) は eF H (k 1) のリー変換の時点で決定され、以後の k より高次
る k 次の斉次多項式 Hkのリー変換では変化しないことが分かる。すなわち、

H k(1) = H k(k) (= H k(k 1) + Fk H2 ):

(142)

 (k)

それゆえ、Fk は eF H (k 1) のリー変換時点で「Hk の関数が見通しがよくなるような」正準変換kを与えるように決めてやればよい。

k 次斉次多項式である H k(k 1) をH k(k 1) =X(m )(n ) z
0

0

m zn
0

(143)

0

(m )(n )Pと書きあらわそう。ただし、和は (m0 + n0 ) = k を満たす任意の非負整数 m0 ; n0 についてとる。
i

一方、

i

0

0

i

i

1 N ff ; ! z z g = N i! z @fk z @fk
Fk H2 =
2 i=1 k i i i i=1 i i @zi i @ zi
X

X

!

i

(144)

であるから、k 次の任意の斉次多項式 fk

fk (z; z) =X(m)(n)

21(m)(n) z z

m n

(145)

P

( + ni ) = k を満たす任意の非負整数 mi ; ni について和をとる）に対して（ただし、 i mi

Fk H2 =P"

X

(m)(n)(m)(n)

NX
i=1#ni ) zmzn

i!i (mi(146)( + ni ) = k のもと、任意の非負整数の組 (m)(n) に対してとなる。 i mi

NX
i=1!i (mini ) 6= 0

(147)

1

ならば、振動数 ! は次数 k まで有理数上 次独立であるという。Fk は任意であったので、その係数(k 1)(m)(n) を任意に選ぶことができる。振動数 ! は次数 k まで有理数上 次独立であるとき、Hkのうち m0i16= n0i であるすべての項を取り除くように (m)(n) を選ぶことができる。すなわち、
X(m )(n ) z0m zn + Fk H2 = 0

0

0

(148)

0

(m )6=(n ) (k 1)i と記すと、逆に消去できなかった部分を hH
0k

H k(k)0= hH k(k 1)iNm=(m )(n ) (zi zi )XX0(m )(n )0== (k )X(m )(n )00X(m )(n )00(m )(n )
0

i

i
N
X

(p2i + qi2)m(151)(2Ji)m

(152)

0

i

0

(m )(n )
0

0

0

0

(149)
(150)

i
N
X

0

i

0

i

147）式がすべての次数に対して成り立つならば、振動数 ! は有理数上 1 次独立であるといい、H は作用だけの関数としてべき級数が収束することが期待される。ここで重要なことは（151）式から明らかなように、となり、Hk は作用変数 Ji の巾級数として表されることが分かる。もし（標準形はすべて pi および qi の偶数次の項から構成され、奇数次の項を持たないことが分かる。

-Deprit の（単一）リー変換による正準変換、例えば、（68）および（71）式、でも同様に解くことができ、Dragt-Finn の手法のところでも述べたように、少なくとも 3 次以上の高次展開でなければ単一および逐次リー変換に違いはない。また、「（Birko と Gustavson の）標準形」と「（堀、Deprit および Dragt-Finn らの）リー正準変換摂動理論」は同一線上に比較すべ
きものでは全くなく、前者は最低次を 2 次とする斉次多項式展開で書き下せるハミルトン系に対する正準変換後の \形"、すなわち、「作用変数の巾級数展開」に対する呼称であって、後者は非摂動系を可積分とし、摂動パラメータ  の巾展開で近似可能なハミルトン系に対する \方法論" でちなみに、堀ある27 。また、標準形およびリー変換の手法・概念自体は保存力学系に限らず散逸力学系についても用いられている28 ことを付記しておく。

27 歴史的に見ると、新旧の正準変数を用いた母関数を用いる Poincare と Zeipel の正準変換摂動理論では高次まで

Birko

Gustavson

2とは最低次を 次とする斉次多項式で表されるハミルトン系
摂動展開が困難であったのに対し、に対しては新旧の正準変数を用いた母関数を用いても高次までの展開が容易であることを見出した。
28 例えば、
J. Sound and Vibration, 217, 795-806 (1998)

A.Y.T. Leung and Q.C. Zhang,222.2.3 小さい分母の問題
Birko -Gustavson の現実に存在するほとんどのハミルトン系に対し作用だけの関数、例えば、標準形、に無限次まで展開可能な系は存在しない。それは摂動計算を破綻させる（準）共鳴条件
が相空間全域を通じて密に分布していることによる。これはカオスの発生に密接に関係しており、かつ後章でも述べる化学反応のサドル領域の反応方向の「頑健な」規則性の起源とも深く関係している。
一般に、非摂動ハミルトニアン H0 が可積分である任意の N 次元ハミルトニアンは H0 の作用
変数

J と角変数  を用いて以下のように表すことができる。
H (J; ) = H0 (J) + n Hn(J; )n= H0(J)
1
1+  H1m(J)eim + 2 H2m(J)eim + : : : ;
X

例えば、

X

X

m

m

H1m(J)eim 

ここで、mk は任意の整数で、H1mで非摂動系 H0

(J) の振動数はによって与えられる。
さて、一般に !k は作用変数

k

N
Y
k=1

H1m (Jk )eim  ;
k

(153)

k

k

(Jk ) は k 番目の作用 Jk に依存するフーリエ係数である。ここ!(J) = @H@0J(J)(154)

J を通じて相空間上の場所（作用が保存する場合には初期点）に依存して変化し得るとしよう。例えば、N 次元モース振動子系
p2k
H0 (p; q) =+Dk e 2a (q q )k=1 2mk
!2
N
X
a
k Jk
! H0(J) =Dk 1 p
2mk Dk
k=1

N
X

"

k

k

eq
k

2e

ak (qk qkeq )

#


(155)
(156)

では、k 番目の振動数は@H (J)
2Dk!k = 0 = ak
@Jk
mk
s

1

p2amk JkD
k k

!

(157)

となり、Jk の値に依存して変化する。

0

リ─変換の処方箋に基づくと、 次のオーダーでは1

また、 次では、H 0 (J) = H0 (J);) +H 1 = H1 (J; 23)dW1 (J; 
d

(158)
(159)

 に依存していない  平均部分を H 1 とすると、
 )i ;
H 1 (J) = hH1 (J; (160)となり、平均化の原理に従って H1 のなかの

および

W1 =

Z 

 ( 0 ) g
d 0 fH1 J; 




(161)

となる。それゆえ、H1 の振動部分についてのフーリエ成分を積分することで W1 を解くことができる。

) = iW1 (J; X

H1m(J) ime :m6=0 m  !(J) に依存しない部分なので W1 に含まれない。
ここで m = 0 は (162)162）式から答えは一般には
)(=「否」であることが想像がつくであろう。基準振動数の任意のある整数倍の組み合わせ m  ! (JN m ! (J)) が W の分母に現れる。各振動数 ! は作用 J
 の連続関数であるから、（準）共鳴1k
k=1 k k
さて、この摂動展開は収束することが約束されているだろうか？（

P

条件

N
X
k=1

mk !k (J) ' 0(163)

が相空間上に存在する可能性がある。このように、相空間上の場所に依存して小さな分母の問題（すなわち、共鳴）が誘発される現象を一般に非線形共鳴(Nonlinear Resonance) と呼ぶ。W1 の分母が零に漸近すると W1 は発散し、摂動展開はその時点で破綻する。小さな分母の問題は W1 に

CPT 計算過程で現れる展開項のそれぞれに対し発生し得るため、作用だけ
 " は、一般にはその摂動展開は収に依存する形で求めようとした新しい多次元ハミルトニアン \H

限らず、一般に n の

束せず、その存在を否定されることになる。
仮りに ! が作用変数に依存しない可積分系を零次ハミルトニアンとして最初定義したとしても、

 (J) を零次ハミルトニアンとし、無視した項 H H (J) を摂動項として再 (J) が求まらない限り」この問題は普遍的に存定義することは常に可能であり、「無限次まで H在する。すなわち、Birko -Gustavson の標準形においても、真の意味で積分可能であるか否かはある有限次で求めた H

級数の収束性に問題をすり替えている。すなわち、実際の計算では有限次で級数展開を打ち切り
それ以上を無視するが、無視した（無限次に至る）高次項において発散する可能性は否めない。
すなわち、振動数の組が有理数比でなくても、あらゆる非負整数の組について和をとる過程で、
（

147）式

N
X
i=1

!i (mini ) 6= 0

が非常に小さな値をとるとき（準共鳴とよぶ）、標準形が収束する保証はなくなってくる。一方= 1）では、小さな分母の問題は存在せず、リ─摂動展開は収束することが期待される（すなわち、1 次元連続力学系は常に作用変数だけからで、すべての一次元（連続）力学系（すなわち、N記述されるハミルトニアンに変換可能（可積分）である）。
従来、（相空間全域に渡る大域的な運動不変量を問う視点から）発散源になっているいくつかの角変数を H に含めたり、あるいはある有限次で摂動展開を打ち切り、与えられた系を「積極

24

\

"

的に」積分可能系に近似して、小さな分母の問題からくる発散を回避しつつ、 見通しの良い ハ

ミルトニアンを抽出する試み29 が行われてきた。統計性と力学的決定性の中間が問われる系にお(p; q) の運動方程式による軌跡 (p(t); q(t)) に沿って、この新しい作用変数 Jk (p; q) および振動数 !
 k (p; q) を「運動の \局所" 不変量の検出器」
いて、最も有効な用途は変換前のハミルトニアン Hとして用いるものである。なぜならば、一般に自由度の数が多くなればなるほど、大域的な運動不変量の占める割合は急速に減少すると予想されるが、運動不変量はある「局所」においては有意に生き残っている可能性があるからである。反応過程において（特に自由度の数が増えれば増えるほど）作用が相空間のどの領域においてどのくらいの時間スケールで、保存し崩壊するのかという見方こそが「統計性と力学決定性の中間を問うべき」化学反応のダイナミックスを理解す
るうえで最も重要となる。発散源になっている特定の角変数を新しいハミルトニアン
る共鳴ハミルトニアンの方法は\H " に含め

\H " の発散を回避することができるが、同時に「その部分空間の局所的な運動不変量を検出できる可能性」を暗黙のうちに放棄していることに留意されたい。

2.2.4 サドル領域における反応方向の作用不変性
1（ヘシアン行列の純虚数固有値がひとつである）ランク のサドル領域における相空間はどの1ような構造をしているのであろうか？ランク の鞍部点はヘシアン行列の負の固有値がひとつある不安定固定点なので、定義より純虚数の振動数 !F

残りの実数の振動数 ! B

1(J)(2 <) をもつ安定な非反応自由度（まとめて B と呼ぶ）から構成される。この領域では 次のリ─母関数 W1（（(J)(2 =) をもつ不安定な反応自由度（F ）と) = iW1 (J; X

H1m(J) im
e

(164)

m6=0 m  !(J)162）式）の !F (J) を含むすべての項は小さな分母の問題は生じないことが分かる。一方、反応性モードを除く、それ以外のすべての項は一般に小さな分母の問題に常に晒らされている。すなわち、 !F(J) を含む項においては、分母はいかなる組み合わせに対しても有限の下限 j!F (J)jが存在するため発散は生じ得ない。

jm  !(J)j  j!F (J)j > O(n):

(165)



このような「部分」非線形共鳴が生じる場合、新しいハミルトニアン H はどのような構造を取り1得るのであろうか？ランク の鞍部点近傍でポテンシャルをテーラー展開したハミルトン系を取り扱えばよいので、

Uzer ら30 に習って Birko の標準形を使って考察してみよう。

この局所領域のハミルトニアンを以下のように表記しよう。

1 q p + N 1 ! z z + k 0
m n m n2 1 1 i=2 2 i i i k (m)(n) (m)(n) z z q1 p1= H2 + k Hk+2

H (p; q) =

X

X

X

1

1

X

k

(166)
(167)

29 換言すると、与えられた系の真の意味での積分可能か否かの問いは角変数を包むこと／摂動展開の収束性の問題に責任を転嫁する。

30 T. Uzer, C. Ja e, J. Palacian, P. Yanguas and S. Wiggins, Nonlinearity, 15, 957-992 (2002)

25

+ ni = k のもとの任意の非負整数 m、n に対する和をとることを意味する。q1 および p1 は反応性自由度の座標系31 を、zi および zi (2  i  N ) は非反応性自由度の座標系を表し、振動数 (= ijj) は純虚数、!i は実数である。この場合（146）式は
PP
ここで、 0 は N
i=1 miFk H2 ="

X(m)(n)(m)(n)となり、標準形の議論から

jj(m1 11) +

N
X
i=2

i!i (mi#ni ) zmzn q1m pn1
1

1

(168)

0
m n m1 n1(m)(n) (m)(n) z z q1 p1 のうち、Pjj(m1 n1) +

N
X
i=2

i!i (mini ) 6= 0(169)

である斉次多項式は Fk H2 によって取り除くことができるので、それらはすべて新しいハミルトニアン Hk から取り除かれることになる。ゆえに、あらゆる任意の振動数 ! の組に対して、常に反1( = n1 ) 項だけ（すなわち、（152）

応性自由度 に関しては q1 および p1 に掛るべき数が等しい m1
m1

式より作用変数 J1だけ）が Hk に含まれる。それゆえ、一般にサドルの近傍領域における非可積分ハミルトン系をH (p; q) = H 0 (J) +X

n=1

n H n (J1 ; B );

(170)



に変換する正準変換は任意の次数において存在する。ここで  B は非反応性自由度の作用および角変数(JB ;  B ) をまとめて表している。

級数の収束性について 確かにサドルの近傍領域で任意の非可積分ハミルトン系を（170）式に変換する正準変換は一般に存在するが、次の意味において一般に収束しない。不安定固定点周り
では  の高次摂動展開を実行すればするほど、一般にサドルの局所領域から離れ、より大きな非線形高次摂動項を摂動計算に取り込まなくてはならない。反応モード F に沿った変曲点に漸近するほど、非摂動ハミルトニアンが零に近づくので、摂動計算の前提自体が成り立たなくなる。しかし、サドル近傍の局所領域32 に関する限り、相空間上の反応自由度の抽出および非反応性自由度のあいだの相互作用を無視できる程度に小さくすることは可能である。

1993 年に Hernandez と Miller が Van Vleck 摂動理論に基づく半古典論の研究に関する速報誌

33 のなかで「サドル領域における反応自由度と非反応自由度のあいだに小さな分母の問題は現れ1999 年に小松崎と Berry が初めてリー変換をサドルを交差する化学反応ない」と指摘している。

ダイナミックスの解析に適用し、カオスの海のなかにあっても反応方向の自由度は作用が保存する規則的な弾道運動であることを数値的に立証した。これは、上述のようにサドルのある近傍領域では任意の非可積分系を次の形式にH (p; q) = H 0 (J) +X

n=1

n H n (J1 ; B );(171)

31 既に (p; q) ! (z; z) のような変換を行ったものと考えればよい。
32 サドルエネルギーよりも遥か上空の高エネルギー領域では振動数 !k (JF ; B ) を介して反応および非反応自由度の間の独立性が壊れる。例えば、
J. Chem. Phys.

T. Komatsuzaki and R.S. Berry, Proc. Nat. Acad. Sci. USA 78, 7666 (2001);

115, 4105 (2001).

33 R.Hernandez and W.H. Miller, Chem. Phys. Lett. 214129(1993)26

変換する正準変換が一般に存在することを意味している。（

 B )) に従う自由度 F の運動方程式
H (JF ; B ) (B  (JB ; 170）式で表されるハミルトニアン

qF (p; q) + !_ F q_ F (p; q) + ! F2 qF (p; q) = 0;!FpF (p; q) =!F _q (p; q);! F Fは容易に導くことができる。ただし! F = ! F (JF ; B ) =

_





(172)
(173)

@ H (JF ; B )
:
@ JF

(174)


であり、x および x は（H に従う）時間 t についての x の一次および二次微分を表す。JF は時間 t (  )さて、この場合、pF (p; q) および qF (p; q) はどのように求めることができるのであろうか？

に対して不変であるので振動数 !F JF ; B は非反応性自由度  B を介してのみ時間 t に依存する。
サドル領域においては、任意のハミルトニアンを不安定固定点近傍で形式的に以下のように巾

級数展開できる。

1 (p2 + !2q2) +  C q q q + 2(175)
2 j j j j j;k;l jkl j k l j;k;l;m Cjklmqj qk ql qm + : : : は一般に
この場合、新しいハミルトニアン H
1(176)
H = !  J + 2 hH2 + fW1 ; H1 gi + O(4 )2と表すことができる（Birko -Gustavson の標準形の議論から分かるように  の奇数次が現れない。
 1 も存在しない）。通常、平均化操作において、それゆえに、2 の項のなかに一般にふくまれる H
H (p; q) =

X

X

X

J のみから構成されるとして定める。ここでは非反応性自由度からなる部分空間の B も加内部の（準）共鳴が生じ、特定の摂動展開項の発散を誘発する可能性を除去するために 

非振動部分は (  )えた。ゆえに、振動数 !k JF ; B は

! k (JF ; B ) =と書くことができる。それゆえ、@ H
= !k + 2 ! k2nd(JF ; B ) + O(4 )
@ Jk
!_ F
' O(2)
! F

(177)
(178)

(172) 式は真っ直ぐのびた糸の長さが「ゆっくり」変化する一次元振り子の運動方程式に対_ F = 0 とおいた \補助方程式" の解の定数 、 および !F を時間の関数としてゆっく
応する34 。!
り変化する（ _ ; _ ; !
_ F  O(2 ) for !F  O(1)）とみなすと、qF (t) = (t)ei! (t)t + (t)e i! (t)t ;
(179)
から、

F

F

34 作用変数の断熱定理を説明するのによく用いられる例である。作用の断熱定理とは「作用が保存する系に対して、\"時間とともに ゆっくり 変化する外部パラメータ a で、その系に摂動を加えた場合、振動の様子は変わり、エネル
ギーも変わり得るが、作用変数 J は一定で変化しない」ことをさすものである。本例ではハミルトニアンが共役な角変数に依存しない形式で表せるために作用が保存する。

27

とおき、（178）式より付加条件として_ ei! t + _ e i! t + i!_ F t( ei! t
F

を用い（

F

F

e i! t ) ' 0

(180)

F

172）式を解いて得られた (t) および (t) を (179) 式に再び代入すると、
qF (t) =

および

Rt

Rt

(0)ei ! (t )dt + (0)e i ! (t )dt ;
0

F

Rt

0

pF (t) = i (0)!F ei ! (t )dt

0

F

(181)

0

Rt

(182)
を得ることができる（上式を（271）式に代入して成り立つことを確認せよ）。また（181）-（182）0F0式から作用が保存することも確認できる35 。

 ()i (0)!F e i ! (t )dt :F00 () =

ij! F (t)j）、反応分割面を横断するのに要する時間はおおよそ反応性双曲型軌道の周期の半分  =j!
F j にすぎない。更には、通過時間に比べて、振動数 !
F が曲率に対して実効的に変調する時間スケールは長い。すなわち、
=j! F j(= O(0 )) << j! F =!_ F j(= O( 2 )). それゆえ、ハミルトニアンが  F に依らず JF が保存する場合は、振動数 !
F (p; q) が運動不変量であるか否かに依らず、!F の時間変化は qF ダイナミッ
サドル近傍領域では、!F t は純虚数であり（!F t

クスの規則性にほとんど影響を与えないことが分かる。すなわち、qF (t) 'および(0)ej! (0)jt + (0)e j! (0)jt ;FpF (t) '

(183)

F

(0)j!F jej! (0)jt (0)j!F je j! (0)jt :(184)
例えば、非常に平坦なサドルで !
F (t) の変調が無視できないほど大きく (183)-(184) 式が近似的になりたたない状況にあっても i t !
 F (t0 )dt0 の「符号」がサドルを通過する過程で変わらない限り、
F

F

R

アプリオリに反応の終状態を予測することは可能である。すなわち、例えば、系が配位空間上の( = 0) を時刻 t0 で、 (0) > 0 の初期条件で交差した場合は、その終状態（ qF > 0 の

分割面 S qF方向）は「その時刻」t0 で予測可能である36 。

Keshavamurthy と Miller の半古典論 Keshavamurthy と Miller は「相空間」反応自由度に1

WKB 理論を

cq1 2
;
m! 2

(185)

沿った虚の作用変数を多次元相空間上から切り出し、トンネル現象に対する 次元

2 自由度ハミルトン系

H (p; q) =1 p2 + p2 + 1 aq2 1 bq3 + 1 m!2 q
2
2m 1 2 2 1 3 1 2に対し展開した37 （彼らは a、b、m、! および c の値は典型的な水素原子移動反応モデルになるように決めているが、c については幾つかの値を用いている）。彼らは結合定数 c の値に依らず、「サドル領域において」すべての新しい作用J は保存し、2 まで求めた新しい H は元の H と近似

35 戸田盛和著 「振動論」培風館 (1968)
36 T. Komatsuzaki and R.S. Berry, J. Phys. Chem. A 106 10945(2002)
37 S. Keshavamurthy and W.H. Miller, Chem. Phys. Lett. 205, 96 (1993); His thesis at University of California,

Berkeley, (1994), Chapter 2, pp 9-36.

28

 (p; q) がサドル領域で局所的に保存する時間を的に等しい値を取ることを見いだした（彼らは J1

トンネル時間と定義し、反応速度を見積もった結果、モード特異性が残らなくなる強結合カオス
領域（＝半古典遷移状態理論が破綻する領域）にあっても、量子力学的に厳密に見積もった反応
速度とよく合致することも発見した）。ここまでの章を読まれた読者は、彼らが見いだしたことは2「自由度の数が であるならば、系の種類および結合定数の強弱に依らず常に成り立つ」ことが想像できよう。なぜならば、サドル領域において、どの（部分）自由度空間をとってみても小さな3( 2) 空間では一般に分母の問題は現れ得ないためである。自由度数 以上では非反応性自由度 >共鳴が存在し得るため、もはやサドル領域では新しい H は元の H と等しくないが、小章で論じ

2

たように、相空間上の反応方向の作用は頑健に保存するのである。彼らの研究は残念ながら 自1

由度系以上見ていないが、彼らの多次元トンネル効果に対する相空間 次元

WKB 理論は彼らが想像した（であろう）以上に、（半古典的に量子効果が論じ得る範囲で）半古典遷移状態理論が使えない強結合カオス領域にあっても一般の多次元化学反応系に適用可能であると思われる38 。

2.3 ポアンカレ断面
我々は、通常、分子動力学計算の結果得られるダイナミックスの様相を理解するとき、実際に系を構成する原子群が反応過程のあいだどのように動くかを映画(movie) にする。これは、特に複雑な生体分子内や溶液内で生起する反応ダイナミックスの「分子論」を展開する際に必要不可
欠であるばかりか本質的に重要である。一方、力学屋さんの立場は「（ポテンシャルエネルギー関数がそうだからという以上に）なぜそのように動くのか」の仕組み・メカニズムを理解しない
ことには分かった気持ちにはなれない。特に局所平衡近似を越えた非平衡かつ非定常な化学反応
ダイナミックスや生体分子系の構造転移における長時間分子記憶効果などを理解するためには両方の視点に立って（精密な分子シミュレーションを得意としてきた）分子科学がこれから挑戦す\べき問題なのである。しかしながら、 どのようにして相空間上のダイナミックスを観たらよいの"
ここでは Henon と Heiles が提案した、銀河系内の円軌道に近い星の運動を記述する自由度数 n
が 2 であるハミルトン系
11
1 q3
H (p; q) = (p21 + p22 ) + (q12 + q22 ) + q12 q2
(186)
2
2
32
を取り上げて相空間構造の振舞いを考えることにする。図 4 に ポテンシャルエネルギー面ならびに q1 = 0 でのポテンシャルエネルギーの断面図を示す。このポテンシャルは正 3 角形の対称性を有し、全エネルギー E = 1=6 を越えたら、系は峠を越えて反応することが予想される。
さて、2 自由度系の相空間次元は 4 であるが、自明な運動の恒量（すなわち、全エネルギー）がひとつあるので、実際には運動は 4 次元相空間上の 3 次元等エネルギー曲面に束縛されているこであろうか ？

38 実は「サドル領域における反応自由度と非反応自由度のあいだに小さな分母の問題は現れない」ことを最初に指Hernandez Miller の論文と同時期に Keshavamurthy と Miller の論文は発表されている。1999 年 5 月米国・
Miller と議論した際、彼は Keshavamurthy の名を \絞り出すように" 思い出して Keshavamurthy
Hernandez の指摘を生かせなかったことを大変悔やんで、筆者の研究を高く評価してくれた。ちな
Keshavamurthy はその後 Ezra のポツドクになった。

と摘したアトランタで筆者がの研究に対して
みに、

29

0.3

V(q1,q2)

0.25
0.2
E=1/6

0.15

1
-1

-0.5

0
q1 0

0.5

-0.5
1

0.5
q2

-1

0.1
0.05
0 -0.5

E=1/8
E=1/12
E=1/24

0

0.5

1

1.5

q2 (at q1=0)

4: Henon-Heiles= 0 でのポテンシャ

図ハミルトニアンのポテンシャルエネルギー面（左図）。q1
ルエネルギーの断面：点線は調和振動子で近似した場合を表している（右図）

5(b) の (1)!(2) 参照）。いま、図 5(a) のように、ある適当な (2n 2) 次元面 R を相空間上に定義し、相空間内の系の時間発展を表す軌道がその面をある方向から（図では \裏" から \表" へ）交差する点だけを取り出して、 R 上の点の離散的な交点の時間発展 xi ! xi+1 を考える。ハミルトニアンとエネルギーの値を決めれば運動方程式から写像 xi ! xi+1 は一意的に決Pとになる（図Pまる。この写像をポアンカレ写像と呼び、このポアンカレ写像の不動点は周期解に対応する。また、「相空間の有効次元数（＝相空間次元2n から自明な運動の恒量の数を差し引いた数、ここで2n 1 = 3 である）から１を引いた次元をもつ曲面をポアンカレ断面と呼び、この場合、図5(b)(3) のようになる。また、n = 3 では図 5(c)(2) に示されるように 4 次元となる。
\このポアンカレ断面は全エネルギーの増加に伴ってどのように変わるだろうか。" 図 6 に
E = 241 ; 121 ; 18 および 16 において Henon-Heiles 系が q1 = 0 の断面を p1 > 0 で (p2 ; q2 ) 面（図中、(py ; y) 面）を交差した点の集まりを図中の右に示した（図中の左は後述する）。エネルギーが低
1 ; 1 ）では写像点列はすべてある滑らかな曲線上にあることが分かる39 。エネルい領域（E = 2412ギーが増大すると（E = 18 ）、曲線の交点近傍の曲線領域は消失し、ランダムな点列が現れるが、は元あった楕円の中心近くでは滑らかな楕円曲線は生き残り、その周りに（滑らかな曲線から成る）\ " \ "
エネルギー閾値（E = 16 ）では、ごく小さな領域を除いて、曲線群は消失し大部分がランダムな点列となる。更にエネルギーが増大し、系が一旦峠を越えて反応した後はポアンカレ断面 q1 = 0

新しい小さな 島 の 鎖 が生まれている。更にエネルギーが増大し、系のエネルギーが反応のを二度と交差することはない。
実は滑らかな曲線が存在することは（その存在する領域において）エネルギー以外に非自明な2運動の恒量が存在し、系は「運動の恒量が系を構成する自由度の数に等しい」 次元トーラス40 に束縛されていることを意味している。これは次のようにして理解することができる。

39 （Poincare-Birko の不動点定理に示されるように）、曲線の交点は不安定な周期解に対応し、楕円の中心には安定な周期解に対応する不動点が存在する
40 変形したドーナツ を連想していただきたい。

\

"

30

5: (a)

(b)2

(1) 4

図
ポアンカレ断面を一方向から軌道が交差する様相、自由度系のケース次元相空間上の 次元等エネルギー面に束縛されているある軌道
その軌道の p1 ; q1 ; q2 空間への射影
q2 一定 の 次元ポアンカレ断面にその軌道がつくる つの写像点列、自由度系のケース
次元相空間上の 次元等エネルギー面に束縛されているある軌道
q3 一定 の 次元
ポアンカレ断面にその軌道がつくる つの写像点列のポアンカレ断面上の点列を p1 ; q1
および p2 ; q2 面へ射影したもの

(3) =
(1) 6

(

3

2

)

5

(2)
3

3

(3) (2)

31

(

)

(c)3
(2) =

4
(

)

(p1 ; p2; q1 ; q2) = c(= 一定) があるとする。
H (p1 ; p2 ; q1 ; q2 ) = E を p1 について解いた p1(p2 ; q1 ; q2 ; E ) を  のなかの p1 に代入すると、q1 = 0
いま、エネルギー以外に非自明な運動の恒量の断面上では

(p1 (p2; 0; q2; E ); p2; 0; q2) =  (p2 ; q2; E ) ＝ c 2 次元面内の曲線の一般式

(187)



となり、 は E をパラメータとする p2 と q2 の関数として表され、ポアンカレ写像によって得られ (p2 ; q2; E ) = c で定義される曲線（すなわち 2 次元トーラス）に束縛されることを要請している。すなわち、2 次元トーラス上の点 xi の写像点列 xi+1 ; xi+2 ; : : : ; x1 もそのトー (p2 ; q2; E ) = c で定義される）独立なトーラスのあい
ラス上に恒久的に存在しなければならず、（る「無限」点列はだを一本の軌道が遍歴することは許されない。逆に、そういった恒量が存在しなければ写像点列は束縛される必要はなく、相空間の等エネルギー面を一様に埋め尽くすことが期待される。例え= 18 では、ポアンカレ断面はトーラスとランダム点列のカオス領域から構成されるが、カオス領域のなかを動く写像点列は決してトーラスに入り込むことはできない。(q1 q2 ) 平面に投影して軌道の振舞いをみると、ポアンカレ断面上にランダムな点列を与える軌道はみな正 3 角形
ば、Eの等エネルギー面をランダムに縦横埋めつくす一方、滑らかな曲線を与える軌道は規則的な（準）3周期軌道を与え、正 角形のある部分領域に束縛されることが確認される。

6

9

Birko -Gustavson の標準形に基づく
1 ; 1 では、標準形より得られたポポアンカレ写像を示している41 。トーラスが観測される E = 2412
図 の左側の図は実は（ 次の斉次多項式まで算出した）

アンカレ断面は実際のそれと良く一致し、このエネルギー領域ではハミルトニアンは実際に標準形に変換できると考えて良い。しかしながら、エネルギーが更に増大すると（E= 81 ; 16 ）、直接数値計算で得られたポアンカレ写像のランダムな点列は決して（有限次の）標準形では表されなくなり、標準形は無限次で発散することが予想される。
化学反応の文脈でこの相空間の振舞いを解釈すると、反応「前」の高振動エネルギー励起状態

Birko -Gustavson の標準形はなんの役にも立たず、従来の化学反応理論が立脚する局所平では、衡を仮定し得る状況にあることを示唆しているかもしれない。しかしながら、反応「前」の高エネ
ルギー振動励起状態にあっても実際の軌道は決してランダムな拡散的過程ではなく、軌道はトー(vague tori、remnant tori、cantori などと呼ばれている) を遍歴しながら等エネルギー
面を \ゆっくりと" 経巡り渡っている42 。このトーラス残骸を遍歴する特徴的な時間スケールが活ラスの残骸性化障壁を越える反応時間に比べて十分短い場合に限り、化学反応理論の局所平衡仮定が成立するのである。また、化学反応における真の反応経路・遷移状態を考察するうえで標準形の概念・
リー正準変換摂動理論は非常に大きな役割を担うことが後章で明らかになる。

n 次元トーラスは、n 個の独立な運動の恒量を伴うので、2n 次元相空間内の (2n 1) 次元の等エネルギー曲面上の n 次元面を構成していることになる。さて、2 次元面を 2 つの領域に分ける
41 Henon-Heiles 系では振動数が等しい（!1 = !2 = 1）ため、通常の標準形は発散項を有する。Gustavson(Astron.
 に含め発散を回避した。6.3 の（準）共鳴のJ. 71670(1966)) は発散項を産む項を積極的に新しいハミルトニアン Hリー正準変換摂動計算の手順がそれに相当する。このとき、H は作用 Ji のみならず (1 2 ) の関数として共役な
角変数 i にも依存することになるが、(1 2 ) = (!1 !2 ) + ( 12 ) = 0 + 1 2 = 一定より  には依存し において個々の作用 Ji は保存しなくなる（∵ @@Hi 6= 0）が、H0 は時間に対して不変である（∵ない。そのため H H 0 g = fH 0; H g = ddtH 0 ）。
0 = ddH = fH;
42 R.B. Shirts and W.P. Reinhardt, J. Chem. Phys. 77, 5204(1982);79, 3173(1983)32

6: Henon-Heiles 系の実際の数値計算によるポアンカレ断面（右）と標準形によるポアンカレ

図
断面（左）。

33

3
2
2
よい。2 次元トーラスは 3 次元等エネルギー面上の 2 次元面を与え、一般にエネルギー曲面を 2つの領域に完全に分離できる。それゆえ、2 自由度系ではカオス領域とトーラスが混在しているには線を引けばよいし、 次元の体積を つの領域に分けるには、 次元面をその境界に使えばとき、一般にカオス領域はトーラスの存在によって等エネルギー面上のある局所的な領域に（棲3

み分け的に）制限されることになる。しかしながら、 自由度のハミルトン系ではトーラスの次

3

5

元は であり、 次元の等エネルギー面を分離することができないため、（トーラス上の写像点列はトーラス上から抜け出せないが）相空間内の異なるカオス領域は少なくともトポロジカルには

3

2

繋がっている。つまり、 自由度以上の高次元系では「 自由度系に比べて」局所平衡が達成されやすいといえる43 。ちなみに、等エネルギー曲面のもつ次元から k を引いた次元のことを余次元k(codimension k) と呼び、2 次元トーラスは余次元 1 であり、3 次元トーラスは余次元 2 である。

1系を反応の「前」と「後」に分割する遷移状態を構築するうえで必要なのは余次元 の多次元曲
面となる。

2

さて、ポアンカレ写像にはもうひとつ重要な性質がある。それは、 自由度系のポアンカレ写1ことである。Lichtenberg-Lieberman の教科書44 に習って略証を示そう。
2 次元ポアンカレ断面（q1 = 0 および p1 > 0）を時刻 t に交差した点 (p2 (t); q2(t)) は次の交差時
刻 (t + 1) における点 (p2 (t + 1); q2 (t + 1)) に（2 自由度ハミルトン系の解の一意性から）一意的に写像される。いま、図 7 に示されるように、2 次元ポアンカレ断面上、時刻 t および (t + 1) での交点をそれぞれ (pt ; qt ) および (pt+1 ; qt+1 ) とする。いま、点 (pt ; qt ) を原点とする微小変位ベクト
ルを d~
pt および d~qt とするとき、時刻 (t + 1) ではその微小変位ベクトルは点 (pt+1 ; qt+1 ) を原点と像は変換のヤコビ行列式の値が となり、「面積」と「向き」がポアンカレ写像の前後で保存する

する d~pt+1 および d~qt+1 にマップされ、写像の前後は次の微分形式で関係づけられるとする。

d~qt+1d~pt+1= @q@qt+1 d~qt + @q@pt+1 d~pt
t

t

t

t

= @p@qt+1 d~qt + @p@pt+1 d~pt

(188)
(189)

pt ; qt ; pt+1 ; qt+1 の 4 変数のなかで独立な変数は 2 個だけであるが、(pt ; qt ) および (pt+1 ; qt+1 ) はハ
ミルトンの正準方程式に従う相軌道上の時刻 t および (t + 1) での点 (p2 ; q2 ) に相当するので、ともに正準共役な変数の組である。すなわち、正準変数の組 (pt ; qt ) から正準変数の組 (pt+1 ; qt+1 ) へ

Lichtenberg-Lieberman では新旧の正準座標系が混在した母関数 F (qt ; pt+1 ) を用いているが、ここではリー母関数 W (pt ; qt ) を用いて考察しよう45 。6.2 章の（293）-（294）式より、時間発展演算子を T とするとの正準変換を可能にするなんらかの母関数が必ず存在する。

qt+1
pt+1

= T qt = exp= T pt = exp

LW (p ;q ) qt = qt+1 (pt ; qt )hiLW (p ;q ) pt = pt+1 (pt ; qt )

h

i

t

t

t

t

(190)
(191)

43 が、実際のところ高次元相空間上の運動の振舞いはあまり良く分かっていない。例えば、C.C. Martens, M.J. Davis,

and G.S. Ezra, Chem. Phys. Lett. 142, 519(1987)。最近の本条氏と金子氏のアーノルド拡散と重なり共鳴の研究。
44 A. J. Lichtenberg and M. A. Lieberman, Regular and Chaotic Dynamics 2nd Ed. (SpringerVerlag, New York, 1992).
45 Lichtenberg-Lieberman の p19-20 の F (qt ; pt+1 ) を用いた方法と対比するとリー変換の簡便さが実感できる。

34

p
dpt+1
pt+1

dqt+1
dpt

dSt

pt

dqt
qt+1

qt
図

dSt+1

q

7: 2 次元ポアンカレ断面上の微小面積 dSt(= jd~qt ^ d~ptj) の写像 dSt+1(= jd~qt+1 ^ d~pt+1 j)
(

)

( +1)

より、 pt+1 ; qt+1 は pt および qt を独立変数とする関数として表すことができる。いま、時刻 t
における外積 d~qt+1^ d~pt+1 を計算すると、
d~qt+1 ^ d~pt+1) d~q ^ d~p= @ (@qt(+1q ;; ppt+1
) t t
t

t

(192)
!

= @q@qt+1 @p@pt+1 @q@pt+1 @p@qt+1 d~qt ^ d~pttttt
= fqt+1 ; pt+1gd~qt ^ d~pt
= fT qt; T pt gd~qt ^ d~pt
= T fqt; pt gd~qt ^ d~pt
= (T  1)d~qt ^ d~pt = d~qt ^ d~pt
;p +1 )となる。ここで @ (q@ (+1（
q ;p ) は相空間ヤコビ行列式であり、
t

t

t

t

(193)
(194)
(195)
(196)
(197)

195）式から（196）式への変換は（291）

式の時間発展演算子の性質を用いた。外積は、d~qt ^ d~pt
d~qt ^ d~qt

= d~pt ^ d~qt
= d~pt ^ d~pt = 0

2

(198)
(199)

などを満たし、異なる つのベクトル（例えば、 d~
qt 、d~pt ）で張られる平行四辺形の面積の大き197）式は 2 次元ポアンカレ写像において任意の向きづけられた面積は保存することを意味し、2 自由度ハミルトン系は 2 次元ポアンカレ断面上では（\Liouville の定理" が成り立つ）1 自由度離散力学系に還元されることになる。
証明は他46 に委ねるが、一般に、n 自由度ハミルトン系では 2n 次元相空間上の任意の点近傍の微小変位ベクトル (d~p; d~q) が張る超曲面の時間発展を考えたとき、各時刻においてこの超曲面をさをもつベクトル量である。それゆえ、（

46 アーノルド 「古典力学の数学的方法」 岩波書店 p226 (1980)35

(pi; qi) 平面に射影した図形の（向きづけられた）面積の和が保存する。すなわち、nXi=1d~pi (t) ^ d~qi (t) = 一定:(200)

ちなみに、この（向きづけられた）面積の総和をシンプレクティック形式、それが保存する写像をシ22ンプレクティック写像と呼ぶ。
（ 自由度系でみたように）エネルギー一定と断面の つの条件によ(2n 2) 次元のポアンカレ断面上では (n 1) 自由度の保存力学系（＝ i=11 d~pi(t) ^ d~qi (t)が保存）が存在する。すなわち、2 次元ポアンカレ写像は面積を保存するが、3 自由度系の 4 次元
ポアンカレ断面では (q1 ; p1 ) および (q2 ; p2 ) に射影した面積は一般に保存しない。
Pn

り決まる

2.4 ハミルトン写像
この章では

H0 (J) ! H0 (J) + H1 (J; )(201)(J) に摂動が加わった系のポアンカレ断面の様子を調べる上で便利なハミルトン写像を紹介する。以下では簡単のため 2 自由度系を例に話を進める。
捻れ写像（twist map） まず、非摂動ハミルトン系 H0 (J1 ; J2 ) において、図 8 のように 2 =const: のポアンカレ断面 (J1; 1) を通過する交差点列を考えよう。このとき交差の前後で時間はのように n 次元の可積分ハミルトン系 H0

図8: (a)2 次元トーラス上の運動 (b)2 = const: のポアンカレ断面 (J1 ; 1) 上の写像点列36t(= !2(J) ) だけ経過するので、1 は !1(J)t(= 2 !! ((JJ))  2 (J)) だけ進む。
2 自由度系のときはエネルギーを E とすると J2 = J2(J1; E ) なので は J1 にのみ依存すると考えてよい。それゆえ、ポアンカレ断面 (J1 ; 1 ) 上の、n 回目の交点 (Jn ; n ) から (n + 1) 回目の交点(Jn+1; n+1) への写像はJn+1 = Jn(202)
n+1 = n + 2 (Jn+1)(203)と書くことができる。ここで便宜上、添字 1 は省略し は Jn ではなく Jn+1 の関数とした。ちな
2

自由度 の一周期

1

2

2

みに、相空間ヤコビ行列に対して@ (Jn+1 ; n+1 )=1@ (Jn ; n )(204)
を満たすことからシンプレクティック形式が保存されていることを確認することができる。図 8(b)中、「外」と「内」の輪は (Jn+1 ) が無理数であるような初期条件から出発させた写像点列を、また Æ; 4 は (Jn+1 ) = 1=6 であるような異なる 2 つの初期条件から出発させた写像点列を n ! 1
までプロットしたものである。後者のように (Jn+1 ) が有理数のとき、写像をある回数施すと元の点に戻る「写像の固定点」が必ず存在することが分かる。（202）-（203）式の写像を捻れ写像（twist map）と呼ぶ。
これらの結果は N (> 2) 自由度系へそのまま拡張することができ、例えば、N = const: を交差する (N 1) 自由度の作用-角変数 (J; ) のポアンカレ写像はJn+1 = Jn(205)
n+1 = n + 2 (Jn+1); = f 1; 2; : : : ; N 1g
(206)と書くことができる。ここで i は !i と !N の振動数比 !i (J)=!N (J) である。
重要なことは可積分系においては「連続時間の軌道計算によるポアンカレ断面の交差点列の系
譜と捻れ写像は等価である」こと、従って「計算時間を大幅に短縮してポアンカレ断面の性質を解析することが可能である」ことである。
摂動を受けた捻れ写像（"

perturbed twist map）（

201）式のように、H0(J) に \僅かな摂( )
Jn+1 = Jn + f (Jn+1 ; n )
n+1 = n + 2 (Jn+1) + g(Jn+1; n)

動 が加わることによって、ポアンカレ断面 J1 ; 1 上に定義された捻れ写像は

(207)
(208)

となると期待される。ここで 、f および g は Jn ではなく、面積保存を満たす条件を容易に見出せ(J; ) が加わることによって、2 の一周期のあいだに
ポアンカレ断面上の 1 の飛びが  の \1 次" で近似できるか否かは自明ではなく、むしろ「（207）（208）式でポアンカレ断面の点列を近似できるような」摂動を受けた系のポアンカレ断面を考え、る Jn+1 の関数とした47 （一般に、摂動 H1

可積分系に摂動を加えた保側写像が有するであろう一般的な性質を論じる立場をとる）。

47 （207）-（208）式のポアンカレ写像は f および g を Jn+1 の関数としたことによって  の 1 次以上の項を部分的にくり込んでいることに当たる。

37

( )

(

 )

さて、 Jn ; n から Jn+1 ; n+1 への写像が向きづけられた面積を保存する写像（すなわち、正準変換）であるならば母関数が存在する。いま、母関数をF2 = Jn+1 n + 2 A(Jn+1 ) + G (Jn+1 ; n )(209)

n+1 = @J@F2 ; Jn = @@F2(210)= dJdA ;
n+1
@G
f =
;
@ n
@G
g =
@J

(211)

とおくと

n+1

であるから、

n

(212)
(213)

n+1

より、向きづけられた面積を保存する写像であるための f と g の条件は

@f
@g
+
=0
@Jn+1 @ n(214)と書き下すことができる。

(J; ) = H0(J) + H1(J; ) の運動方程式dJ1@H (J; )
=
 1
(215)
dt
@ 1
@H (J; )d1=!1 (J) +  1
(216)
dt
@J1
からも関数 f を \O () の次数で" 評価することができる。すなわち、ポアンカレ断面 (J1 ; 1 ) 上の n 番目の交点 (Jn ; n ) から (n + 1) 番目の交点 (Jn+1 ; n+1 ) への遷移のあいだの作用 J1 の変化
量 J1 (= Jn+1 Jn ) は、1 の次数ではT
1 (J ; J ;  + ! (J)t;  + ! (J)t)
J1 =  0 dt @H
(217)
20
2
@ 1 n+1 2 n 1と書くことができる。n 番目の交差から (n + 1) 番目の交差に掛かった時間（周期）T2 は一般に nに依存するが、 の 1 次の範囲では、トーラス面に沿って非摂動軌道を 2 の 1 周期（= ! 2(J) ）に渡って積分すればよい（すなわち、2 番目の自由度については平均操作を施したことに対応する）
ことが分かる。ここで 20 は n 番目の交差時の 2 である。このとき、作用の飛び (jump) は
f (Jn+1 ; n ) = J1 (Jn+1 ; n )(218)と表すことができる48 。また、対応する位相の飛び（jump）は面積保存の条件である（214）式
ハミルトニアン H

Z

2

2

から

g (J; ) =

Z

 @f

@J

d0

(219)

48 （217）式は非摂動軌道に対して積分を取るので被積分関数のなかの J1 は一定であり、 の 1 次の範囲では厳密（207）式の f (Jn+1 ; n ) に対応させるために、被積分関数の引数の J1 を Jn+1 と表した。
には Jn である。ここでは、

38

( 2) 自由度系の捻れ写像は、N = const: を交差する (N 1) 自由度のとなる。摂動を受けた N >-作用 角変数(J; ) のF2 = Jn+1  n + 2 A(Jn+1 ) + G (Jn+1 ; n )(220)

Jn+1 = Jn + f (Jn+1; n)n+1 = n + 2 (Jn+1) + g(Jn+1; n)

(221)
(222)

に対し、

と書くことができる。
素朴な疑問： 捻れ写像のときは自明であるが、非可積分系である
次元ポアンカレ断面の交差点列の時間発展を

2N 次元相空間上の (2N 2）(2N 2) 個の変数から成る写像として記述できるのはなぜか？
答：

2N 次元のハミルトン系のポアンカレ断面の次元は (2N 2）次元である。それゆえ、ポアンカレ断面上の任意の点は (2N 2) 個の変数で「一義的に定義する」ことができる。しかし「定義できること」とその交差点列の時間発展を「写像として書き下せること」は同じではない。交差点列の「時間発展」を記述するのに要する独立な変数の数（次元）は本来(2N 1) 個であり、(2N 2) 個の変数から成る写像の写像点列として交差点列を表現できるということは、ある 1 次元の情報が時間発展を記述するうえで必要ないことを意味する。
（217）式に示されるように、「少なくとも  の 1 次の範囲では」ポアンカレ断面を構成する
自由度の \残りの自由度 " の作用は保存し（ポアンカレ断面の場所に依らない）周期が定
義でき次元をひとつ下げることができる。そのため、(2N 2）次元のポアンカレ断面の交差点列の時間発展は (2N 2) 個の変数から成る写像として表すことができる。
このように、連続時間のハミルトン系に基づくポアンカレ断面の交差の時間発展を離散時間の1写像として定性的に（ の 次の範囲で）捉えることができるため、現実の化学反応ダイナミック

1

スの相空間構造を解析するうえで多用されている。しかしながら、 の 次の範囲でポアンカレ

断面の振る舞いを定性的に評価できるためには、少なくとも問題とする系は近可積分的でなくてはならない。このため、実際の連続時間のハミルトン系に基づくポアンカレ断面の数値計算と相補的に取り扱うべきであろう。
動径捻れ写像 さて「f が Jn+1 に依らず

n の関数として表される」摂動を受けた捻れ写像を取り上げよう。

= Jn + f (n)n+1 = n + 2 (Jn+1)
Jn+1

(223)
(224)

radial twist map）と呼ばれるこの捻れ写像はもっとも良く調べられているハミルトン写像で、この写像は（その写像点列を産み出す）\kicked" 系と呼ばれる時間依存のハミルト

動径捻れ写像（39()ニアンを容易に導出することができる。いま、周期的なデルタ関数 Æ1 t

Æ1 (t) =

1

X

n= 1

Æ (t n)(225)

n は時刻 t = n の直前の J (t) と (t)、Jn+1 と n+1 は時刻 t = n の直後の J (t)と (t) とみなすと（223）-（224）式はdJ@H=f ()Æ1 (t) !
(226)
dt
@
@H
d
=
2
 (J ) !
(227)
dt
@J
を導入し、Jn とと表すことができる。それゆえ、この写像点列は時間 t に依存するハミルトニアン H

H (J; ; t) = 2Z J

(J 0)dJ 0Æ1 (t)Zf (0 )d0(228)に従う時間発展とみなすことができる。

(J0) が整数であるような周期 1 の固定点（J0 = Jn+1 = Jn）近傍で（224）式
を線形化 (Jn = J0 + Jn ) すると、

例１ 標準写像n+1 = n + 2 (Jn+1) ! n+1 = n + 2（!

(J0) + ddt Jn+1 ;
J

(229)

0

223）-（224）式は= In + Kf (n)n+1 = n + In+1; mod 2
In+1

(230)
(231)

となる。ここで In および K は= 2 ddt Jn = 2 0Jn
(232)
= 2 0fmax(233)
であり、K は摂動の強さ  に比例し、f  (= f=fmax ) はポアンカレ断面上の作用変数の \飛び(jump)" をその最大値 fmax によって規格化したものである。特に f  = sinn のときを標準写像（standard map）と呼び、ハミルトニアンは
In
K

H=

I2

2 +K

1

X

n= 1

Æ (t n) cos 

\

"

(234)

となることから、可積分な回転子に周期的な外力（ キック ）が加わった系に対応すること= 0:5  4:0 の範囲の 4 つの K 値に対し、数多くの異なる初期点 (I0; 0) から出発させたとき得られる写像点列 Ik ( Ik modulo 2 ) と k ( k modulo 2 ) を図 9 に示した（図中では In ! p; n !  ）。K = 0 のときは、
In+1 = In(235)
n+1 = n + In+1; mod 2(236)が分かる。K40

図9: 標準写像：(a)K = 0:5; (b) K = 1:0; (c) K = 2:5; (d) K = 4:0

41

= I0 = const: の水平線上に束縛される。このとき、n = 0 + nI0 とな
1 つの初期点から出発した写像点列は In = I0 の線上を密に埋めつくすのに対し、有理数 r=s（r; s は整数）のときは s 回毎に同じ点に再帰し、1 つの初期点から出発した写像点列は In = I0 の線上を密に埋めつくすことはないことが分かる。この有理数のときの不変（曲）線を共鳴トーラスと呼ぶ。図 9(a) に見られるように、(1) このから写像点列は Inることから、 2I0 が無理数のとき系に僅かな摂動が加わった場合は「 2I0 ＝無理数」のトーラスの多くは頑健に生き残っていること（In' I0 ）、(2)「I0 modulo 2 = 0 の共鳴トーラスは破壊され楕円型不変曲線（トー\"ラス）が生まれ、その近傍には小さな島模様の 子供 のトーラスが生まれる。摂動を強く

9(b)）、もはや 0 から 2 を縦断するようなトーラスは消失し、\子供" のトーラ
スの周りに \孫"、\曽孫"、\玄孫" のトーラスが出現するとともにランダムな写像点列（カすると（図オス領域）が現れ始める。しかしながら個々のカオス領域はトーラスに遮られ、お互いの領域を系は行き来することができない。更に摂動を強くすると（図"9(c)(d)）、島模様の \子

供、孫、曽孫、: : : のトーラスは消失し、トーラスで仕切られ棲み分けていたカオス領域が大域的なカオス（確率論的な記述が許されるであろう状況）に発展する。このように、キッ2

ク回転子の相構造を眺めると、可積分系に摂動が加わった（ 自由度）近可積分系が保有する一般的な現象が観測される。

217）式から分かるように、可積分系
H0 (J) に正準座標 q にのみ依存する弱い摂動 V (q)（ O()）が加わったときのポアンカレ
断面の性質を \ の一次の範囲で" 捉えることができる（と期待される）。このため、化学反
応ダイナミックスとカオスの研究において (N 1) 自由度の自由粒子に化学反応を表すポテ
ンシャルエネルギー V (q) を起源に持つ周期 T で離散化された外力が加わったモデルハミル

例２ 化学反応シンプレクテック写像

動径捻れ写像は（

トニアン

（ここで、

H

= H0 + T V (q)

H0

=

N
X1

p2i
i=1 2miX

n=0;1;2;:::Æ (t nT )

(237)
(238)

q = fq1; q2; : : : ; qN 1g）がよく用いられる49 。なぜならば、このモデル系は、実際の N 自由度ハミルトニアン

H=p2i+ ()V 0(q1 ; q2; : : : ; qN 1; qN )
2
m
i
i=1

N
X

に基づくポアンカレ断面の性質を定性的に（(239)

 O()）捉えることができるハミルトン写像

(q )
(240)
! pn+1 = pn T @V
@q n
Tp
qn ! qn+1 = qn + n+1(241)m
49 例 え ば 、P.Gaspard and S.A.Rice, J. Phys. Chem. 93,6947(1989); R.E. Gillilan and G.S. Ezra,
J. Chem. Phys. 94, 2648 (1991)
pn

42

(

1) 個の各自由度に対し）存在するためである（N 番目の自由度の運動が周期的であるとする）。この写像がシンプレクティックであることは相空間ヤコビ行列が 1 であることから確認することができる。例えば、2 次元ハミルトン系の相空間ヤコビ行列は@q@qT
T @@qV mT
1+@ (J ; )
@q
@p
m
=
(242)
@p
@ (J ; ) = @p
T @@qV
1 =1
@q
@p
が（ N

n+1

n

n+1
n



n+1

n+1
n

n

n+1

n+1

n

n



2

2

2

2

となり、シンプレクティック形式が保存することが分かる。

2.5 セパラトリックス、安定・不安定不変多様体、ヘテロ（ホモ）クリニック
交差
この章では図

10 のような double well 型の 1 自由度化学反応系に対し、熱浴自由度が系に摂動を与えるときの相空間構造の変化を考察し、相空間上のポテンシャル井戸の「内」から「外」への移送経路、すなわち、化学反応のチャネル「遷移状態」を統計性を予め仮定しない力学の立場から論じよう。

2

まず、熱浴が存在しないときの 次元相空間を考えてみる。全エネルギー E が活性化障壁 Ea
よりも低いとき、系はポテンシャル井戸内に束縛され秤動(libration) 運動する。このとき、図 10（下）に示されるように相空間では楕円軌道を描くのが分かる（図中の矢印は時間の向きを表している。逆向きの矢印は存在し得ないことを確認すること）。一方、E > Ea では、系はポテンシャ
ル井戸に束縛されることなく、相空間全域を回転(rotation) 運動する。全エネルギー E が活性化

障壁 Ea にちょうど等しいとき、系は無限大の時間を掛けてサドル点に到達し、無限大の時間を掛けてサドル点から離れていく。この相空間上のサドル点を不安定固定点（固定点のほかに不動点、停留点ともいう）と呼ぶ。不安定固定点の「不安定」は相空間上の初期値がほんの僅かに狂うだけでその点から系が離れていくためである。
一般に固定点近傍の相空間の構造はその固定点 p(Æpt ; Æqt) ! (Æpt+1 ; Æqt+1) を考えればよい。
Æpt+1
Æqt+1

!

=

0

@pt+1
@pt
@
@qt+1
@pt

1

= p0 ; q = q0 のまわりで線形化された写像

@pt+1
@qt A
@qt+1
@qt
p0 ;q0

Æpt
Æqt

!

=M

Æpt
Æqt

!

(243)

1

ここでは面積保存の写像を考えているので M のヤコビ行列式は であるという制限がつく。そのため、特性方程式は

22trM + 1 = 0(244)

2 <もしくは=) と trM (= @p@p + @q@q =  + 1 ) の値に応じて固定点近傍の相構造が特徴づけられることが分かる。重根（ = 1）の特別なケースを除となるので、 つの固有値  および 1 （

t+1

t+1

t

t

くと、

1) jtrM j > 2 : 固有値は  1 以外の実数
2) jtrM j < 2 : 固有値は絶対値 1 の複素数、 = exp(i)
43

(245)
(246)

図

10: 1 次元 2 重井戸型ポテンシャルの 2 次元相空間図

2
1)

(

の つの場合が一般に存在する。すなわち、適当な座標変換 Æp; ÆqÆp0t+1
Æqt0+1

2)

!

 0

= 0 1


!

Æp0t
Æqt0

!

) ! (Æp0; Æq0) をすると、各々( > 1)

 sin Æp0t
= cos
sin  cos Æqt0となり、11 図に示すように、1）では双曲線族 Æp0 Æq 0 = 一定、2) では一様な回転を表す同心円族
Æp02 + Æq 02 = 一定 の不変曲線50 の群を構成し、1) に対応する固定点を双曲固定点 (hyperbolic xed
point)、2) に対応する固定点を楕円固定点 (elliptic xed point) という。双曲固定点は運動量零の鞍部点に、楕円固定点は運動量零のポテンシャルの極小点に対応する（6.4 章の固定点近傍の相空Æp0t+1
Æqt0+1!

!

!

間ヤコビ行列を参照のこと）。
このように双曲固定点「近傍」の様子は固定点に近づく不変曲線（安定不変多様体（

stable in-variant manifold）51 と呼ぶ）と固定点から離れる不変曲線（不安定不変多様体（unstable invariant
50 ある写像（この場合、M ）によってその全体の形が変わらない曲線のこと。最初にその曲線上に位置すると永久に同じ曲線上に写像される。
51 通常、n 次元多様体とは「任意の」各点近傍が n 次元ユークリッド空間で幾何学的に表現される空間を指す。例
えば、 次元多様体は曲線であり、その曲線上の任意の点近傍は幾何学的に線で表される。アルファベット文字 O やD は 次元多様体とみなすことができる。O のように任意の点で微分可能であるときを特に微分可能多様体と呼ぶ。
L また X は終点（境界）や線と線の交点の近傍系は幾何学的に線ではないから多様体とはみなせない。L のような境界をもつとき「境界をもつ多様体」と呼ぶことがある。また、ボールやドーナツの表面は 次元多様体である。

1

1

2

44

δq'

δq'

δp'

δp'

図

11: 双曲固定点（左）楕円固定点（右）manifold）と呼ぶ）によって特徴づけられる。一方、固定点から離れた遠方領域では、この 2 つの不変多様体は滑らかに接続し閉曲線をつくっている。これは解が有界であれば 2 自由度可積分系のポアンカレ写像に対しても同様に起きる。この閉曲線は 2 次元相空間をちょうど 2 つに分割するので、この閉曲線を分離閉曲線（セパラトリックス（separatrix））と呼ぶ。すなわち、1 自由度系ではポテンシャル井戸内の（秤動）運動から活性化障壁を越える（回転）運動への遷移は起こり得ない（つまり、恒久的にポテンシャル井戸内に留まるか、井戸内に留まることなく相空間を回転し続ける）ため、化学反応という概念そのものが存在しない。
さて、この系に対して、例えば

H=

1 p2
2 1


!12 q12 + bq14 ! H =1 p2 + p2
2 1 2
!12 q12 + !22 q22 + bq14 + cq12 q2



(247)

のような熱浴自由度 q2 が加わるとどのように相空間の構造は変化するのだろうか？双曲固定点

2

は熱浴自由度による摂動が加わっても、少なくとも 次元ポテンシャルエネルギー面上に鞍部点
@V
（ @q
1

=

@V
@q2

= 0 で、かつヘシアン行列の 2 つの固有値が正と負である点）が存在するならば、周期解に対応する不安定固定点がポアンカレ断面上に必ず存在するのが分かるであろう。また、その双曲固定点近傍において、その点にちかづいていく安定（不変）多様体 W s 52 およびその点から離れていく不安定（不変）多様体 W u もポアンカレ断面上に同様に必ず存在する。しかしながら、双曲固定点から離れた領域においては、熱浴自由度の影響により（無限大の時間を掛けて）双曲固定点からでてくる不安定多様体 W u と（無限大の時間を掛けて）双曲固定点に入ってくる

安定多様体 W s はもはや滑らかに接続することは一般にあり得ないであろう。そして、有界のと2

き一般に有限の角度をもって つの多様体は交差することになる。同じ双曲固定点から発する安(homoclinic point) とよばれ、異なった双曲固定点から発する安定・不安定多様体の交点はヘテロクリニック点 (heteroclinic point) とよばれる。
図 12（左）のように双曲固定点にはいっていく W s を逆発展させたものと双曲固定点から離れ定・不安定多様体の交点はホモクリニック点ていく W u が最初に交差した点を P0 と記そう53 。さて、双曲固定点から発した不安定多様体 W u
は P0 で安定多様体 W s と交差した後、どのようにのびてゆくだろうか？

52 この時点で多様体に境界（終点）が存在するか否か分からないので、厳密な意味で多様体か否かは分からない

\

"

\

が、この後の議論で、この不安定 多様体 は双曲固定点へ（無限大の時間掛けて）近づき接続し、この安定 多様
体 も双曲固定点から（無限大の時間掛けて）辿ってきたことが分かる。
53 第一（ホモクリニック）交差点（と呼ばれる。

"

primary intersection point,pip)

45

1)

12

図 （右）のように、P0 は安定多様体 W s 上に乗っているため P0 の写像点列 P1 ; P2 ; : : : ; P1

も安定多様体 W s 上になければならない。この写像点列は無限大の長さの時間を掛けて双曲固定点に近付いていくため、Pi と Pi+1 の安定多様体 W s 上の間隔も i! 1 につれて無限に小さくなっていくことになる。

P4

図

2)

P3

P2

P1

12: ホモクリニック点 P0; P1; P2; : : :また、P0 は不安定多様体 W u 上にもあるから、写像点列 P1 ; P2 ; : : : ; P1 は不安定多様体

W u 上にもなければならない。すなわち、不安定多様体 W u は振動しながら安定多様体 W sと交差しなければならない。まず図13（左）のような交差が考えられよう。しかしながら、

2 次元ポアンカレ写像は向きづけられた面積が保存しなければならないので、この写像はd~1 ^ d~1 = d~0 ^ d~0 に示されるように、面積が保存したとしても向きが逆になる写像となる。それゆえに、向きづけられた面積が保存するためには図 13（右）のように、ホモク

リニック点列 P1 ; P2 ; : : : ; P1 の間に別のホモクリニック点列 Q1 ; Q2 ; : : : ; Q1 が存在する必要がある。また、

lobe で囲まれた面積（影の部分）を保存しながら54 双曲固定点に近づいていくので影の領域は指数的に細長くなっていく必要がある。そして双曲固定点近傍では W u

はほぼ平行になって無限大に引き延ばされるが、同時にそれは等エネルギー面上に束縛されているので無限に折り畳まれていくことになる55 。つまり、相空間上の初期点がほんの僅かズレるだけでも双曲固定点近傍では指数関数的にズレが増大する。

3)

ホモクリニック点 P0 の逆写像による点列 P 1 ; P 2 ; : : : ; P 1 に対しても同様適用され、安定多様体 W s も不安定多様体 W u と無限回交差しながら P 1 の双曲固定点近傍では無限に引き延ばされ折れ畳まれている。図

14 のように、双曲固定点に近づくにつれて、安定・不安定多様体の無限個の交点ができ、より「稠密に」分布するが、写像の一意性から安定多様体および不安定多様体とはそれ自身と決して交差しないことに注意しよう。また、定性的にも手で描くのが極めて困難であることは、例えば図

深く観察すると良く分かるであろう56 。

14 中のホモクリニック点列 X を注意

54 影をつけていない「内」側の lobe と「外」側の lobe の面積は同じであるが、それはどう考えれば分かりますか（宿題）。
55 是非、一度コンピューターで実際に計算をして確認してほしい。
56 入門書として有名なののは対応するホモクリニック点 を探せば、あのが間違っていることが分かる。そでは図 （左）に基づく議論をしているのもよろしくないけど）。
れくらい大変という意味（

Michael Tabor \Chaos and Integrability in Nonlinear Dynamics" John Wiley & Sons,
NY. (1988) Figure 4.21
X
FigureFigre 4.21
13

46

δζ1

P4

P3

P2

δξ1

13:

P1 δξ0

δζ0

δξ1

lobe
^ =
^
^ = ^

δξ0 δζ0δζ1

図W s と W u のホモクリニック交差の様子：に囲まれた面積は等しいが向きが違う、間~~~~違った 次元ポアンカレ写像（d1 d1d0 d0）（左）、に囲まれた面積も向きも等しい、正しい 次元ポアンカレ写像（d~1 d~1 d~0 d~0 ）（右）

2

2

lobe

X

図

14: 無限個のホモクリニック点とその集積の様子

47

図15 に Gaspard と Rice57 が解析した分子解離反応の安定・不安定不変多様体のホモクリニッ

ク交差の様子を示した。彼らは分子の解離反応のモデルであるモース型ポテンシャル

図15: 分子解離反応の相構造（d = 1:8）：安定・不安定不変多様体のホモクリニック交差
2

V (q ) = D [1

exp(

q )]2

(248)

T
p ; qn = qn ;
m n

(249)

を用いた 次元動径捻れ写像を実際には用いている。

pn =と変換すると、（

240）-（241）式はpn+1 = pn + d [exp ( 2qn ) exp ( qn )] ;57 P.Gaspard and S.A.Rice, J. Phys. Chem. 93,6947(1989)

48

(250)

qn+1 = qn + pn+1 ;(251)となり、ハミルトン写像はひとつのパラメータ dd=

2 2T 2D (> 0)(252)のみで規定される。不安定不変多様体 W u および安定不変多様体 W s は、不安定固定点 (p; q ) =(0; 1) 周りに相空間ヤコビ行列を線形化して得られる行列の 2 つの固有ベクトル上の「微小」領m

域に密に分布させた初期点群を正および負の方向に時間発展することで近似的に描くことができ

15(b)）。図 15(a) に示されるように、(p; q) 平面上に t = 0 で分布させた初期分布（図中、0とラベルしている）が 0 ! 1 ! 2 ! 3 ! 4 ! 5 ! 6 と時間発展するに従い、引き延ばされ、折る（図れ畳まれていく様子が分かる。このように（トポロジー的には）「引き延ばしと折れ畳みの繰り返し」がカオスを作り出している。

2.5.1 ホモクリニック点の存在と 1 価の解析的積分の非存在
双曲固定点近傍はホモクリニック点が集積していることが分かった。このとき、（少なくとも）1

この近傍領域には 価の解析的積分（つまり、非自明な運動の恒量）は存在しないことが次のようにして理解することができる。

1)

1

まず、非自明な 価の解析的積分があると仮定し、それから導かれるポアンカレ写像の() = 一定 とする。W s と W u は不変 \曲線" なので（不変量 I (p1; q1) があれば）W s と W u は I (p1 ; q1 ) の 1 つの等高線を成すはずである。
2) つまり、W s または W u 上の任意の点において I (p1; q1) の W s または W u 方向の微分は 0不変曲線族を I p1 ; q1

となる。

3)

(

ホモクリニック点では W s と W u は有限の角度で交差する（平行ではない）から、I p1 ; q1

1

2

0

の 次独立な つの方向の微分は となる。つまり、勾配ベクトル

4))rI (p1; q1) は 0 となる。

rI (p1 ; q1) = 0 の点が双曲固定点近傍では集積しているので、双曲固定点近傍領域では恒等的に rI (p1 ; q1 ) = 0 でなければならない。つまり、その領域では I (p1 ; q1 ) は勾配ベクトル

変数 p1 と q1 に依らない定数関数でなければならず、最初の仮定と矛盾する。したがって

1価の解析的積分、すなわち、運動不変量、は存在しない。

6 Henon-Heiles 系のポアンカレ写像において E < 121 で滑らかな曲線点列のみ見つ

例えば、図 のかったことはこのエネルギー領域においては、双曲固定点から離れた領域で双曲固定点へ近づく

W s と双曲固定点から離れていく W u が滑らかに接続しているか、または、ほぼ角度 0 で交差していることを意味している（ただし、このときの双曲固定点および楕円固定点はポテンシャルの鞍部点および極小点ではなく、井戸の「内」領域の相空間上の固定点であることに注意されたい）。
また、E= 121 ! 81 と系のエネルギーを大きくした際、E = 121 において存在していた曲線写像点列の交点の近傍領域の曲線点列がランダムな点列となった理由は、エネルギー増大に伴って非線形性が増大することによって W s と W u が有意な大きさの角度をもって交差したためである。

49

2.5.2 相空間上の移送-大域的な遷移状態さて、安定・不安定多様体のホモクリニック交差は化学反応を考える上でどのような意味があるのであろうか？実はホモクリニック交差が存在することでポテンシャル井戸の「内」から「外」への移送（およびその逆移送）、すなわち、「化学反応」現象が創出されるのである。
図16 にポテンシャル井戸の「内」から「外」へ抜けていくポアンカレ写像点列 X1; X2; : : :（左図）とポテンシャル井戸の「外」から「内」へ入ってくる写像点列 Y1 ; Y2 ; : : :（右図）の一例を示す。双曲固定点から出てくる不安定多様体と入っていく安定多様体が P0 で交差をもつまでの両切片から作られる「セパラトリックス」に注目すると、写像点列は一旦ポテンシャル井戸の「内」から「外」へ、もしくは「外」から「内」へ、このセパラトリックスを越えれば再びそのセパラト

1

リックスを横断することがないことがわかる。このセパラトリックスの余次元は なので等エネ

2

ルギー面を つの領域に分割する。ゆえに、このセパラトリックスはまさに反応の前後を反応系
および生成系に分割する「再交差を与えない」遷移状態と位置づけることができる。ポテンシャ

Y1

Y2

Y3
Y4

X1

Y8
Y7

X5

Y9
Y6

Y5

X2
X4

X3

16:

図
ポテンシャル井戸の「内」から「外」への写像点列 X1 ; X2 ; : : :（左）とポテンシャル井戸の「外」から「内」への写像点列 Y1 ; Y2 ; : : :（右）ル井戸の「内」からセパラトリックス遷移状態を越えて「外」へ向かう写像点群がつくる、ポア

17（左）に示した（当然、各々の領域（＝影の部分）の面積は等しい）。
このとき、\一回の写像操作で" セパラトリックス遷移状態を（「内」から「外」へ）越えた（＝反応が進んだ）直後の写像点群がつくる領域は turnstile58 (＝回転扉) と呼ばれる領域の「耳たぶ」
領域 A（図 17（右））であることが理解できるであろう。もし図 18 に示す、反応「前」の領域
がエルゴード的であるならば、一回の写像につき反応が生起する確率 (E )（領域 A と 領域 のンカレ断面上の領域を図

面積はエネルギー E に依存する）は、

(E ) =

領域 A の面積
領域 の面積( )k(E ) = (E )  単位時間当たりの領域 A の生成率(253)となり、実際の反応速度定数 k E は

58 M.J. Davis and S.K. Gray, J. Chem. Phys. 84, 5389 (1986) で最初に名づけられた。

50

(254)

と表現される。

turnstile

B
A
17:

図
セパラトリックス遷移状態を越えて「外」からポテンシャル井戸の「内」へ向かう写像点群がつくる領域（左）と＝回転扉 と呼ばれる A B 領域（右）turnstile()

図
反応「前」の領域( + )

18: 反応「前」の領域

がエルゴード的であるか否かは一般に系およびエネルギー値に依存する。

2 自由度化学反応系 HeI2 に基づく具体的な解析は Davis と Gray に関する後の章（レジュメをお願いする）で議論する。

2.6 法双曲的不変多様体

2.5.2 章で、局所平衡を予め仮定しない力学の立場から
1) ホモクリニック交差、すなわち、カオス、が存在しなければ化学反応という概念そのものが存在しない2)

化学反応過程を相空間の双曲固定点から発する安定多様体および不安定多様体の不変構造
によって理解することができる

3)

反応系と生成系を分割する「再交差軌道を与えない」遷移状態は鞍部点近傍に存在する必要はなく、むしろ相空間の全域に広がっている51

ことが示された。もしこの考えを多自由度化学反応系に適用することができれば、IVR（Intra-

molecular Vibrational energy Redistribution）に始まる、すべての非統計性に由来する困難を解決し、「局所平衡」仮定および「非再交差」仮定を前提としない新しい化学反応理論を構築することができる。化学反応理論をカオスの立場から再構成する意義がここにある。

3

しかしながら、系の自由度の数 N が 以上の多自由度系において、N 次元トーラスが2(2N 1)

次元の等エネルギー面を つの領域に分けることはトポロジカルに不可能であることが知られる

2

ように、N 次元不変多様体を遷移状態とする理論は多自由度系には適用できない。すなわち、 N

2次元相空間上の等エネルギー面を つの領域に分ける分離曲面を構成するためには、なんらかの(2N 2) 次元の不変多様体を考えなければならない。
N > 2 の多自由度系において、2 次元ポアンカレ断面上の双曲固定点の不安定周期軌道に対応するのが、Wiggins59 が提唱する法双曲的不変多様体（Normally Hyperbolic Invariant Manifold,
NHIM）（ここでは M と表すことにする）である。法双曲的 (Normally Hyperbolic) とは法線方向の双曲性（＝軌道不安定性）が接線方向の双曲性よりもはるかに強いことを意味する。法双曲的不変多様体の重要な性質は、もしそれが存在するならば、摂動に対して、その不変多様体は頑健に存在する60 ことである（＝構造安定(structurally stable) であるという）。多次元ポテンシャ

NHIM を本稿の表記に従って書くと、
M = f(q1 ; p1; q2; p2; : : : ; qN ; pN )jqF = pF = 0; H 0(JB) + nH n(JB ;  B ) = E g

ル曲面上の鞍部点領域に存在する

X

n=1

 = pF = 0 と H 0 (JB) +となり、エネルギー E 一定のもと、qF

P

(255)

 (JB ;  B ) = E を満足

nn=1  Hn(p; q) の集合を意味する。すなわち、相空間上に切り出した反応自由度 F の座標qF (p; q) と運動量 pF (p; q) をともに零とするので、仮りに NHIM が共鳴条件 m  !B = 0 を内部で満たしているとしても、NHIM は不変集合をなすことが理解できるであろう。また式からも示されるように、NHIM の次元は (2N 3) であり、NHIM は余次元 2 の不変多様体である。
「法双曲的である」ことから帰結されるもうひとつの大事な性質は「余次元が 1 である」NHIM へ近づいていく安定多様体 W s (M) と NHIM から離れていく不安定多様体 W u (M) が存在し、それらも構造するすべての安定であることである。これも多次元相空間上に切り出した反応自由度 F に沿った作用保存の頑(NHIM) 内でカオスが発生していても、鞍部領域の相空間上では反応方向の次元は 1 次元的である。すなわち、NHIM

健性の議論からも分かるように、たとえ非反応性自由度のつくる部分空間から発する安定・不安定多様体はトポロジー的には等エネルギー面を分離することが可能である。

NHIM から発する W s(M) と W u(M) がホモクリニック交差をするならば、第一ホモクリニック交差に至るまでの両切片を繋ぐことで、多次元相空間を反応の「前」と「後」の 2 つの領域に分離する (2n 2) 次元の \セパラトリックス" を構成することができる。すなわち、2 自由度
仮りに系のときと同様の反応理論を展開することができることになる。ただし、このホモクリニック交差はもはや「点」ではなく、(2N 4) 次元の多様体61 を構成する。

59 S. Wiggins, Normally Hyperbolic Invariant Manifolds in Dynamical Systems, (Springer-Verlag, New York,

1991).

60 なめらかに変化し、急に消失することがない。
61 簡単な幾何学的考察から (2N 2) 次元のユークリッド空間上で (2N

不変多様体が交差を一般

3) 次元の不変多様体と (2N 3) 次元の(generic) に有するとき、交差の次元は (2N 4)[= (2N 3) + (2N 3) (2N 2)] とな

52

2

後章で議論するように、この法双曲的不変多様体の概念を用いて、 自由度系において成功を納Gillilan とEzra62 である。彼らは 3 自由度化学反応系 HeI2 に対して、あるパラメータに対しては、W s(M)と W u (M) のあいだにホモクリニック交差が存在し、2 自由度系のセパラトリックス遷移状態理

めたセパラトリックス遷移状態論を多次元化学反応系に拡張する試みを最初にしたのが論が適用可能であること、しかしながら、別のパラメータに対しては、ホモクリニック交差が存在しないことを見いだした。戸田63 は多自由度系においてはホモクリニック交差の存在は一般には保証されず、ホモクリニック接触を境に「カオスのトポロジー的な性質が変わる」分岐現象がUzer、Ja e と Wiggins ら64 が小松崎と Berry のリー正準変換摂動理論による Ar6 の構造異性化の研究にヒントを得て任意の鞍部点領域で NHIM を計算する手法を確立し、NHIM を実空間上に可視化することに成功した。
生じることを示した。近年、

3 鞍部領域の相空間構造
3.1 遷移状態分光および原子クラスターの \相転移" に見た鞍部領域の可積分性

1980 年代後半からの実験研究の進展は目覚ましいものがあった。例えば、遷移状態の直接観測により、反応速度はエネルギーの増加に対して階段的に増大し、近似的に作用が保存する規則的\' 高いエネルギー領域に至るな遷移が「反応の活性化障壁に対応するエネルギー閾値から あるまで」存在することが実験的に立証された65 。すなわち、遷移は時として単なる確率過程ではな
く近似的な運動不変量を伴う協同的な動力学的過程なのである。また、分子内振動エネルギー再分配（

Intra-molecular Vibrational energy Redistribution, IVR）のダイナミックスの階層性の発見66 も、相空間は決して一様でなく、フラクタル的な階層性を有していることを意味している。
いずれも、化学反応ダイナミックスは完全に統計的でもなければ決定論的でもなくその中間に位置することを強く示唆している。これら実験に関する、力学系カオスの専門家から見た詳細な紹介がカオスとの関わりに力点を置きながら文献67 （和文）になされているので参照されたい。

\"

有限少数多体系における 相転移 現象：鞍部点領域に発現する規則性

Berry らは希ガス原子ク

ラスターの少数多体有限系に観測される \擬相転移" 68 を非線形力学的手法に基づいてその非一様

31
1 + 1 3 < 0）2
1(= 2 + 2 3)
62 R.E. Gillilan and G.S. Ezra, J. Chem. Phys. 94, 2648 (1991)
63 M. Toda, Phys. Rev. Lett. 74, 2670 (1995);Phys. Lett. A 227, 232 (1997)
64 S. Wiggins, L. Wiesenfeld, C. Ja e, and T. Uzer, Phys. Rev. Lett. 86, 5478(2001); T. Uzer, C. Ja e, J.
Palacian, P. Yanguas and S. Wiggins, Nonlinearity, 15, 957-992 (2002)
65 E. R. Lovejoy, S. K. Kim and C. B. Moore, Science 256(1992)1541.
66 K. Yamanouchi, N. Ikeda, S. Tsuchiya, D.M. Jonas, J.K. Lundberg, G.W. Adamson and R.W. Field, J. Chem.
Phys. 95(1991)6330.
67 戸田幹人、「化学反応の動力学とカオス」物性研究 74(2000)597.
68 擬相転移とは 2N hEkini （hEkin i は軌跡に沿った運動エネルギー Ekin の平均値）として定義されるマイクロカ(3N 6)kB（一旦、動的な固液共存相を越ノニカルな \温度"T の増大とともに、系が固体類似相から液体類似相へ状態を変え、えると）各相に長時間滞在する「反応」現象を指す。無限自由度の多体系における相転移とは異なり、相転移 \温度"る。例えば、 次元空間上で（ 次元の）線と線が交差をもつことは一般にはあり得ない（つまり、である。
が、（ 次元の）面と面は一般的に交差をもち、その次元は近傍、ある有限の温度範囲で系が固体類似相と液体類似相の間を動的に移り変わる固液共存相が存在し、その温度範
囲は系を構成する原子の数が増大するに従って小さくなる。

53

Amitrano と Berry69 は Ar3 および Ar7 の擬相転移ダイナミッ

な動力学的特性を詳細に解析した。

クスにおいて、「局所」最大リヤプノフ指数（最大リヤプノフ指数を幾つかの有限時間毎に評価す\"る）のスペクトルを解析し、 相転移温度 付近において（最大リヤプノフ指数の）短時間スペク
トルに２つのピークを見いだした（長時間スペクトルでは単一ピークのみ）。彼らは短時間領域では、鞍点（

saddle 70 ）領域のダイナミックスとポテンシャルのベイスン領域のダイナミックスが分離し、前者のほうがより小さい最大局所リヤプノフ指数を伴う擬規則的ダイナミックスである

\

"

\

"

こと、他方、 相転移温度 以下71 および以上の 温度 領域では、どの時間領域でも単一ピーク\"見いだした。Hinde と Berry72 は希ガス原子クラスター (原子数 N =3-7) の構造転移のダイナミッだけが観測され、かつ 相転移温度 の鞍点領域のダイナミックスよりも強いカオスを示すことをクスを系統的に調べ、ポテンシャルエネルギーの幾何構造のどのような性質が鞍点領域の規則的ダイナミックスを決めているのかを、局所\"

Kolmogorov エントロピー (局所 K エントロピー)、（マイクロカノニカルな） 温度 、モード間の結合強度、瞬間の K エントロピー（鞍点における虚の振動数の絶対値の総和を取ったものに相当する）を解析し詳細に調べあげた。その結果、彼らは平坦な鞍点ほど、鞍点近傍におけるモード間の結合が部分的もしくは全体的にdecouple し、遷移

ダイナミックスを規則化する傾向にあること、また、活性化障壁よりも全エネルギーが大きくなるにつれて遷移状態の交差運動の規則性が減少することなど、を見いだした。彼らはトラジェク
トリーに沿った瞬間の基準座標に基づく局所作用解析（

Local Action Analysis）も開発し、系が鞍点に近づくと次第に作用が保存する傾向があることを解き明かした。

Berry らの手法は各瞬間の動力学のトラジェクトリーを局所安定構造に Quench することによってどの瞬間に系がポテンシャルエネルギー曲面上の異なるベイスンを乗り移っているかが帰属できる場合に有効であるが、原子数 N= 3; 4 に比べて N = 6 以上では、明確な動力学描像を描くことができなかった。原子の数が増大するに従って、一般に、より複雑に絡み合った非線型結合がモード空間全域を通して存在するようになるため、ほとんどすべての自由度に渡って作用が保存する擬規則的ダイナミックスから高次元カオスへと移行しやすくなる。このため、反応過程における作用保存性を解析する上でどのような「座標」を選択するかがより大きな問題となる。

3.2 遷移状態理論と Grote-Hynes 理論の等価性：反応における「系」とはなにか

1930 年代に提案されて以来、遷移状態理論のひとつの大きな ambiguity であ

非再交差仮説は、り、対処策として再交差を極力減らすように反応分割面を反応座標 q1 に沿って変分的に最適化( = 0) ! S (q1 = c),c:任意の定数）する変分型遷移状態理論や再交差挙動を反応の進行を阻害する \分子摩擦" として捉える Kramers-Grote-Hynes 理論が既に存在していた。van der Zwan

（ S q1

69 C. Amitrano and R. S. Berry, Phys. Rev. Lett. 68, 729 (1992); Phys. Rev. E 47, 3158 (1993)
70 「saddle（もしくは鞍点）」とは、断らない限り、ここでは任意の qi に対する勾配 @V (q) が零で、ヘシアン行列@qi2@ V (q)の固有値のうち負の成分がひとつある配位空間上の点もしくはその点の近傍を意味する。
@qi @qj
71 勿論、調和振動子近似が成り立つような温度に比べたら遥かに高温である温度領域。
72 R. J. Hinde and R. S. Berry, J. Chem. Phys. 99, 2942 (1993)

54

とHynes73 は、1983 年、非反応性の楕円型軌道と反応性の双曲型軌道が双一次結合を介して結合している多次元調和振動子ハミルトン系（すなわち、完全な可積分系）においては、反応性座標と非反応性座標の一次結合として構成される（ヘシアン行列を対角化して得られる）不安定固有ベクトルを反応座標としてみなす遷移状態理論とオリジナルの双曲型の反応性（基準）座標を反応座標とみなす

Kramers-Grote-Hynes 理論が、局所平衡が成立する条件下、厳密に等価であることを証明した。

3.3 高次元相空間上の遷移状態：カオスの海に埋もれた規則的経路
近年、小松崎らはリ─正準変換摂動理論を鞍部点領域に対し展開し、「系に依らず普遍的に」遷(p; q) の座標分離の程度）に関して階層構造が存在すること、殆んどの非反応性座標は互いに結合し強いカオスを呈するのに対し、相空間上の反応座標 q1 (p; q)移状態領域には運動の規則性（〜は有意な距離に渡って他のモードから独立であるエネルギー領域が比較的高いエネルギー領域に\"渡って存在し得ること、従って、 ランダム に揺らぐ運動に埋もれた規則的な遷移を与える「反 (p; q)」に基づく相空間遷移状態理論が導出できること、従来のいかなる方法でも不可

応座標 q1

能であった多次元相空間上の反応のボトルネックを可視化し、複雑な再交差パターンの起源を読み取ることができること、などをに示した74 。

Ar6（Lennard − Jones potential）の構造転移を用いて具体的

Ar6 は時定数が大きく異なる特定のモードがなく、容易に（準）共鳴条件を満たす

強いカオスが遷移状態近傍に出現し、局所リヤプノフ解析では遷移の規則性が見い出せなかった

2

系である75 。彼らの研究成果の一部を概説しよう。このクラスターは 種類の極小エネルギー構= 12:712"）は正八面体（octahedral）配位 (OCT) で、次安定構造（E = 12:303"）は 5 原子が三角の隣り合った２つのピラミッドを形成し、6 番目の原子がそのうちの一つの面を覆う capped trigonal bypyramid 配位 (CTBP) を形成している。更に、異なる２つ（負曲率モードをひとつ有する第一ランクの）鞍点を有し、最も低い saddle I（E = 12:079"）
は OCT と CTBP を、より高エネルギーの saddle II（E = 11:630"）は CTBP に関する２つの置換 （permutational）異性体を繋いでいる。
まず不安定停留点（saddle）周りにポテンシャルエネルギーをべき級数展開し、造を有し、最安定構造（EH = H0 +H0 (p; q) =
1

X

n=1

1

n Hn ;

(256)

1 (p2 + !2q2) = ! J = H (J);
2 j j j j j=1 j j 0(257)

n Hn (p; q) = 

X

n=1

X

X

j;k;l

+ 2

X

(258)

Cjklqj qk ql

Xj;k;l;mCjklmqj qk ql qm + : : : =73 G. van der Zwan, J.T. Hynes, J. Chem. Phys. 78(1983)4174
74 T. Komatsuzaki and R. S. Berry, Adv. Chem. Phys. 123,79(2002)
75 R. J. Hinde and R. S. Berry, J. Chem. Phys. 99, 2942 (1993)

55

X

n=1

n Hn (J; ):(259)

Potential energy [ε ]

-10.5

+1.0ε
+1.0ε

-11.0

+0.5ε

saddle II

saddle I
+0.5ε

-11.5

+0.1ε

+0.1ε

-12.0

saddle I

CTBP/CTBP’minima-12.5

OCT minimum"reaction path"

図19: 6 原子アルゴンクラスターのポテンシャルエネルギー面(3N -6) 自由度のハミルトニアンを求めた。ここで、安定な非反応性モード B (!B 2 <) および不安定な反応性モード F (!F 2 =) に対系全体の並進および回転の自由度を（近似的に）取り除いた

する作用─角変数は

1
i p2F
JF = ImpF dqF =j!F jqF2 ;2 barrier
2 j!F j = i tanh 1 j!F jqF ; !  j! ji:
!

Z

(260)

!

F

pF

F

(261)

F

である。反応性モードに対する作用変数は半古典遷移状態理論に対してMiller が最初導入したもの76 で量子力学における半古典論の透過積分に相当する。反応性および非反応性モードに対するJ

作用─角変数は各々純虚数および実数であり、不安定反応性モードも含めて および

 は正準であることが容易に証明される。このハミルトニアンに対し、我々は各鞍点エネルギーを基準とす= 0:1; 0:5; および 1:0" のそれぞれにおいて、10,000 本の `well-saddle-well' トラジェクトリー
を発生させて異性化反応の遷移ダイナミックスの規則性を詳細に調べた77 。(p; q) ／ (p; q) 空間
における反応座標を q1 ／ q1 、および非反応性座標を q2 ; q3 ; : : : ; q12 ／ q2 ; q3 ; : : : ; q12 と表記する。
ここで、番号付けは !2  !3  : : :  !12 ／ !
2  !3  : : :  ! 12 とした。
E =0.05 および 0.5" において saddle I を交差した軌跡を 0〜2 次の qj (p; q) と qk (p; q) から成る
２次元面に映してみよう。ここで、0 次の座標系 (p0th; q0th) は元の (p; q) 系に対応し、以下の図
20,22 は各軌跡をそれぞれの次数 i の (qjith (p; q), qkith (p; q)) 平面に投影したものを重ね描きした図であることに注意してほしい。まず最初に、非反応性自由度のなかから q3ith (p; q) と q8ith (p; q) を選んで、その 2 次元平面に投影してみよう。図 20（右）は E = 0:05" における「ひとつの」軌道をるE

76 W.H. Miller, Faraday Discussions Chem. Soc. 62,40 (1977)
77 両鞍点領域における基準モード間の 3 体および 4 体結合の総数は saddle I で 106 および 365、saddle II で 189 および

674 である

56

0.10

0.3

0th order

q8 (p,q) [σamu A]

0th order
0.05

1/2

1st
order

0.1
0.0

ith

−0.1

ith

1/2

q8 (p,q) [σamu A]

0.2

−0.2

0.00

−0.05

−0.3
−0.8

−0.4

0

ith

0.4

1/2

0.8

−0.10
−0.2

−0.1

ith

0

1/2

0.1

0.2

q3 (p,q) [σamu A]

q3 (p,q) [σamu A]

図

2nd
order

1st order

2nd order

20: (q3ith (p; q),q8ith(p; q)) から眺める（saddle I）。右：E = 0:05"、左：E = 0:5"

( (p; q), q8ith (p; q)) 平面に投影したものである。活性化障壁よりも僅かに高い 0.05" では、すべての摂動次数において近似的な Lissajous 図形が観測される。このことは、この２つの非反応性の自由度は近似的に分離しており「少なくともこの 2 自由度に関する」運動は単純な調和振動子
各 q3ithの重ね合わせとして近似できることを示唆している。「なぜ、化学反応系は行きつ戻りつしつつ、活性化障壁を越えるのであろうか？また、なぜ越えられずに元の状態に戻ってしまうのであろう

21 は E = 0:05 における「峠 S (q1 = 0)」を一旦越えたのちに、元来た方向へ戻ってしまった再交差軌道を、同様に、配位空間と相空間平面に平面に投影したものである。CTBP から登ってきたこの軌道はいずれの相空間上の反応分割面 S (q11st (p; q) = 0) もしくは S (q12nd (p; q) = 0) もか？」図

横切っていない。このことからひとつの重要な結論が導かれる。すなわち、「もし相空間上の高次の反応座標（例えば、q12nd ）に沿った作用が遷移過程において保存するならば」、配位空間上の( = 0) に対して観測された、あらゆるすべての \非反応性の再交差" は、相空間q12nd (p; q) = 0)）を横切らない軌道として解釈される:(p; q) 空間に
上の反応分割面（例えば、S (おいて作用が保存するということは qj (p; q) と qk (p; q) のあいだに結合がないことに対応し、一反応分割面 S q1

旦、系がその反応分割面を越えた場合、系を元来た方向へ戻す如何なる「力」も存在し得ない。
0.20

dividing surface

1st

0.10

ith

1/2

q4 (p,q) [σamu A]

0th

0.00

2nd

−0.10

OCT side
−0.20
−0.02

CTBP side
0

0.02

1/2

ith

0.04

q1 (p,q) [σamu A]

図21: (q1ith (p; q),q4ith(p; q)) から眺める。E = 0:05",saddle I

全エネルギーが活性化障壁から更に10 倍増加した場合はどうであろうか？このエネルギー領域
57

では近似的な作用不変量は消失し、非反応性モード空間上の運動は規則的なものから完全にカオ

20（左））。活性化障壁よりもエネルギーが高く、遷移が殆んどすべての自由度に渡ってカオスである場合は再交差挙動が極めて顕著に現れる。0:5" にス的なものに移行しているのが分かる（図0.6

dividing surface

1st

0.2

0th

0.0
−0.2

ith

1/2

q4 (p,q) [σamu A]

0.4

2nd

−0.4

CTBP side

OCT side
−0.6
−0.04

−0.02

0

ith

1/2

0.02

0.04

q1 (p,q) [σamu A]

図22: (q1ith(p; q),q4ith (p; q)) から眺める。E = 0:5",saddle I( (p; q), q4ith (p; q)) 平面に投影した図 22 に示されるように、座標空間上に定義された「峠（＝ S (q1 = 0)）」を何度か縦断したのち \やっと" 反応した軌道は、高次の「峠（＝ S (q12nd (p; q) = 0)）」に対しては唯一、一度だけ横切って生成系へと至っている。更には、反応座標 q12nd に沿った軌跡は、鞍点領域 0:04 < q12nd < 0:04 に渡って、その反応分割面 S (q12nd = 0)
へ戻ろうとする「力」が一切働いていない。他方、0 次および 1 次の反応座標に対しては、その反応分割面のごく近傍および面から離れた領域のいずれにおいても、他のモードと結合している.
このような広範囲なエネルギー領域（E = 0:05 0:5"）に渡り高次の相空間上の座標描像に立脚おけるカオス軌道を q1ith

することによって、配位座標系

q ないしは低次の相空間上の座標系 q において観測された反応性の「見掛けの」再交差は一回通過の決定論的な遷移として焼き直すことができる。つまり、殆どの再交差挙動は観測者の視点に依り、再交差仮説は観測者の視点を相空間に移すことによって 仮説

ではなくなるのである。
これは、非反応性モードに対する基準振動数は、複素振動数空間上における虚軸に直交な実軸に帰属されるため、反応性モードがひとつ含まれる場合は如何なるモードの組合わせを考えても（準）共鳴条件は満たされないことに起因する。このように、高エネルギーな高次元カオス領域であっても、遷移状態に至る運動と遷移状態を離れる運動のあいだの動的相関は（少なくとも遷移
状態領域において）顕著に存在する。

Kramers-Grote-Hynes 理論の等価性再考 従来の配位空間上の反応分割面
S (q1 = 0) の代わりに、相空間上の反応分割面 S (q1 (p; q) = 0) を用いれば、KGH 理論を包含す
遷移状態理論とる一般化遷移状態理論を再構成することができる。すなわち、「遷移の過程において遷移の終状\"態を決定するのに充分長い 短い時間 の間、相空間上の反応自由度方向に対して近似的な作用( = 0) を再交差する場合にあっても、一度横切った相空間上の S (q1 (p; q) = 0) へ再び系を戻す手段（すなわち、力）が存在し得ないため、q1
は \真" の反応座標と見做すことができる。再構成されたマイクロカノニカルな遷移状態理論によ不変量が保存する場合」は、系が配位空間上の S q158

( (p; q) = 0) を横切る一方向の流れ
j+ (= q_ 1 (p; q)h(q_ 1 (p; q))) の小正準集団平均として与えられる：る反応速度定数 k T ST は、与えられたエネルギー E の下、S q1

kT ST (E ) = hZj+ iE = hq_ 1 (Zp; q)Æ [q1 (p; q)] h [q_ 1 (p; q)]iE ;= 1 dq1 dp1 : : : N dqN dpN
 Æ [E H (p; q)] q_ 1 (p; q)Æ [q1(p; q)] h [q_ 1(p; q)] ;(262)
T ST (E ) の反
ここで、h(x) および Æ (x) はそれぞれ Heaviside 関数、Dirac のデルタ関数である。 k
応速度定数の実験値 kexp(E ) に対する比は、新しい透過係数 cc (E )  kexp(E )=kT ST (E )(263)として定義され、振動エネルギー緩和が速く反応系における局所平衡を仮定でき、かつ核による\ "量子効果が無視できる場合には、観測者の視点に殆んど依らない 真 の再交差がどの程度発現したか、すなわち、完全に発達したカオスへの漸近の程度、を示す指標と考えることができる78 。

( )( ; ( = 0)) をオリジナルの H (p; q) に従う軌跡計算に基づいて調べてみよう
79 。
「系のエネルギーが増大するにつれて MD
q1ith = 0)) はどのように振る舞うのだろうか？」c (t; S (q10th = 0)) は S (q1 = 0) に基づく従来の MD (t) に対応すること
ここで 0 次の反応分割面 MD
c (t; S (に注意しよう。saddle I および saddle II 領域、E = 0:1; 0:5; 1:0" における MD
q1ith = 0)) をc (t; S (図 23 に示す。配位空間上の反応分割面に基づく MD (t) は（非常に短い時間領域を除くと）１か
ここでは、c E のうち遷移状態領域の「局所」のダイナミックスに焦点を絞って「再交差」のみに由来する MDt S q1ithc

ら顕著にズレており、そのズレの大きさは系のエネルギーが増大するに従って大きくなることが()分かる。図中の MD t は時間が長くなると平坦な領域に漸近するが、これは再交差トラジェクト
リーが最終的に終状態へ到達しその遷移状態領域に滞在した時間よりも遥かに長い時間80 生成系に留まり、（再）交差が終結したことを反映している。この平坦値が（実験によって観測される）

saddle I における透過係数  はどのエネルギーにおいても、saddle II に比べて１からのズレが大きいことが分かる。このことは saddle I の領域では従来の配位空間上の反応
座標 q1 が他の非反応性座標 qk (k 6= 1) により強くカップルしていることを意味している。一方、
相空間描像に立脚すると、低いエネルギー (E = 0:1") および比較的高いエネルギー (E = 0:5")の両領域において、saddle に依らず、摂動計算の次数を高くすれば、MD
q1 (p; q) = 0)) は
c (t; S (1 に漸近することが分かる。活性化障壁よりもかなり高いエネルギー  0.5" にあっても、収斂値
c (S (q12nd (p; q) = 0)) は殆んど１である：0.9999 (saddle I) および 0.9996 (saddle II)。系の全エ
ネルギーが更に大きくなる（E = 1:0"）と、収斂値 c (S (q12nd (p; q) = 0)) は１よりも顕著にズレてくる。この「1 からのズレの大きさ」は従来の  に当たる。

78 後述するように、\完全に発達したカオス領域" は相空間の反応座標方向にあっても作用は終状態を決定するのに充分長い時間保存しない領域として定義される。

79 T. Komatsuzaki and R.S. Berry, J. Chem. Phys. 105(1999)10838.
80 生成系における「部分」エルゴード性が完全に成立する時間に比べればまだ短い時間である。

59

1.000E = 0.1ε
MD

κ (t)

0.5ε

0.994

1.0ε

0.992
0.990

0

2

4

6

8

10

MD

MD

ith

0.996

κc (t){S(q1 (p,q)=0)}

0.998

ith

κc (t){S(q1 (p,q)=0)}

1.000

12

E = 0.1ε

0.998

0.5ε
0.996

1.0ε

0.994

0.992

14

0

t [ps]

23:

MD

κ (t)

2

4

6

8

10

12

14

t [ps]( ; ( (p q) = 0)) ( = 0 1 2)

=0.1, 0.5, and 1.0

saddle図
新しい透過係数 MDt S q1ith ;i ; ; （E"）左：c
、右：実線、点線、太線はそれぞれ 0，1，2 次の反応分割面に対するものである。saddle I にI

saddle IIおける収斂値は 0.9988(0), 0.99996(1), 1.00000(2) (E = 0:1"); 0.9940(0), 0.9987(1), 0.9999(2) (E = 0:5");

0.9912(0), 0.9949(1), 0.9982(2) (E = 1:0")、saddle II における収斂値は 0.9991(0), 0.99995(1), 1.00000(2)(E = 0:1"); 0.9958(0), 0.9986(1), 0.9996(2) (E = 0:5"); 0.9931(0), 0.9948(1), 0.9955(2) (E = 1:0") である。括弧内の数字は CPT の次数を意味する。

\2 次" の CPT がどの程度不十分か、
２．仮に CPT が無限次まで展開でき q1 に沿った運動不変量が抽出できたとしても、その収束半
１．反応座標 q1 をその他の非反応性座標から分離する上で径は非常に小さくなる高エネルギー領域は必ず存在する（であろう）から、非再交差仮定が\"成り立つような 有意な長さの 相経路が消失する高エネルギー領域に系がどれくらい接近しているかを表している。

saddle I のほうが saddle II に比べてこの「ズレ」が小さいということは、saddle

I のほうが q1 に沿った近似的な運動不変量がより良く抽出できていることを物語っている。つまり、より弱いカオスを有する遷移ダイナミックス（saddle II）が必ずしも q1 に沿った良い運動不変量を有することに対応していないことに留意されたい。我々の結果は不安定停留点近傍の任意の非線型多次元ハミルトン系において、遷移状態理論とKramers-Grote-Hynes 理論が、系が殆ん

どカオス的で運動不変量が殆んどない、より一般的な条件下においても等価であることを示唆している。

\

相空間上の反応のボトルネック（本当の峠）を 見る" この相空間上の反応分割面は反応座

 (p; q) がある一定値（＝零）を持つという意味においては従来の反応分割面と同じであるが、

標 q1

各モードの運動量および系全体のエネルギーにも依存する非常に複雑かつ抽象的な相空間上の超\多次元曲面である。この抽象的な超曲面の構造を捉えることはできないだろうか？更には 見る"ことを通じて多次元配位空間上の複雑な再交差挙動を統一的に理解できないだろうか？

( = 0) を越えているのであろうか？ここ

その前に、実際に系はどのように配位空間上の峠 S q1\"

では、配位空間上の峠を 越えきれなかった 、すなわち、最終的に反応しなかった再交差軌跡の60

1 OCT と CTBP を繋ぐ saddle I に対し E =0.1,0.5および
1.0" における 10,000 本の well-saddle-well トラジェクトリーのうち、峠 S (q1 = 0) を \越えきれな
かった" 再交差軌道を交差のパターンによって分類した。
交差パターンを考察してみよう。表 に、

表2: 配位空間上の峠 (saddle I) を \越えきれなかった" 再交差軌道のパターンによる分類
E=0:1" (= 11:979") 0:5" (= 11:579") 1:0" (= 11:079")(1,1,{)
161
329
369(2,2,{)
7
23
22(3,3,{)
0
1
3(4,4,{)
00
0(1,1,+)
26
67
86(2,2,+)
1
6
12(3,3,+)
00
0(4,4,+)0
1
0()( = 0) を特定の方向から何回横切ったかを表す

ここで N+ ; N + ;  は与えられた分割面 S q1

インデックスである：もしある軌道が最初の交差の際に  符号の交差速度 p1 を有し、その分割

面を「正」から「負」の方向へ N+ 回交差し、「負」から「正」の方向へ N + 回交差した場合()
saddle I において、与えられた分割面を２度交差し OCT から CTBP 安定構造から最初交差している場合のインデックスは (1,1,+) となる。表から、CTBP 次安定構造から登ってきたトラジェクトリーは、OCT 最安定構造から登っ
は N+ ; N + ;  タイプの交差であった、とする。例えば、てきた場合に比べて、その峠を一旦越えたにも関わらず、最初来た安定構造に戻りやすいことが= 0:1" における 161(1,1, ) >> 26(1,1,+);E = 0:5" における 329(1,1, )>> 67(1,1,+)。なぜ、OCT から登る場合に比べて CTBP から登る \登山者" は峠 S (q1 = 0) を分かる81 。例えば、E越えづらいのだろうか？ 従来、如何なる方法を用いても「遷移状態を通過する際の遷移の容易\"さに対する、このような 登ってくる方向 に依存した違い」を説明することはできなかったが、\"我々は相空間の反応分割面をある少数自由度空間に投影して 見る ことによって、この問いに答

2

( )
( )
へ各エネルギーにおける S (q1 (p; q) = 0) を投影してみよう。式で書けば、(qi ; qj ) の場合、\超曲
面"S (
q1 (p; q) = 0) を投影した \分布"S(qj ; qk ; E ) はえることができる。ここでは、元の座標および運動量から成る 次元平面 q1 ; q2 および q1 ; p1

S(qj ; qk ; E ) = hÆ [q1 (p0 ; q0 )] Æ (qj0ZZqj )Æ (qk0

qk )iE= 1 dq10 dp01 : : : N dqN0 dp0N
 Æ [E H (p0; q0)] Æ [q1 (p0; q0 )] Æ(qj0qj )Æ (qk0

qk );81 構造が同じ置換異性体を繋ぐ saddle II の場合は統計誤差の範囲内で、当然、同じ頻度である。また、\wellsaddle-well" トラジェクトリーが相空間上の反応分割面 S (q1 (p; q) = 0) をどのように交差しているかについては T.

Komatsuzaki and R.S. Berry,Phys. Chem. Chem. Phys. 1, 1387 (1999) を参照のこと。

61

(264)( (p; q) = 0) が従来の配位空間上に定義され\ "と書くことができる82 。q1 への投影は 真 の峠 S q1\"

( = 0) からどのくらいずれているかを、また、p1 への投影は真の峠の形・分る 見掛け の峠 S q1

布が如何に鞍点を通過する速さに依存しているかを、把握する上で重要である。ここで、基準座\"( = 0)標描像が成立するくらいの活性化障壁 すれすれ のエネルギー領域では、相空間上の S q1
は従来の配位空間上の面 q1= 0 と一致するはずであることを想起してほしい。

図 24-25 に S (q12nd (p; q) = 0) を saddle I および saddle II に対して E =0.1 および 0.5" の条件下、(q1 ; q2) 平面に投影した S2nd (q1; q2 ; E ) を示す。 両 saddle に対して、全エネルギー E が増大する

(b)

(a)

1.0

1.0
0.5

0.5
-0.15 -0.10

0
-0.05

0

q1

図

0.05 0.10
-1.0
0.15

-0.5

-0.15 -0.10

q2

0
-0.05

0

q1

0.05 0.10
-1.0
0.15

-0.5

q2

24: saddle I における S2nd(q1 ; q2; E )、(a)E =0.1"; (b)E =0.5"

(a)

(b)
1.0

1.0

0.5
-0.15 -0.10

0
-0.05

0

q1

図

0.05 0.10
-1.0
0.15

-0.5

q2

0.5
-0.15 -0.10

-0.05

0

q1

0.05 0.10
-1.0
0.15

-0.5

0 q
2

25: saddle II における S2nd(q1 ; q2; E )、(a)E =0.1"; (b)E =0.5"

 (

; )

( = 0) から離れた領域へ大きくずれることが分かる。特筆すべき点は、saddle I では相空間上の真の峠 S (q12nd (p; q) = 0)、つまり、反応のボトルネックの重心は CTBP 次安定構造側 (q1 > 0) よりも OCT 最安定構造側 (q1 < 0) に存在している点である。q1 に対する S (q1 (p; q) = 0) の非対称性が峠 S (q1 = 0) に対する「\登山"のルートの違いによる \踏破" の難易度の差」の力学的起源である。つまり図 26 に示すように q1
につれて、分布 S 2nd q1 ; q2 E は拡がり配位空間上の分割面 S q1

に対して、

1. !: 系が「負」から「正」へ \見掛け" の峠 S (q1 = 0) を通過する場合 系は S (q1 = 0) へほ
とんど戻り得ない。なぜならば、S (q1 = 0) を通過した時点で真の峠 S (q1 (p; q) = 0) が多く
82 実際の計算は 10,000 本の \well-saddle-well" トラジェクトリーを用いて、それらが反応分割面 S (q1 (p; q) = 0) を交差した瞬間の座標・運動量情報から頻度分布を作成して算出した。
。

Chem. Phys.

1, 1387 (1999)

62

T. Komatsuzaki and R.S. Berry,Phys. Chem.

0( = 0) へ戻す \力" が系に働く確率が小分布している領域（q1 < ）を通過している為、S q1
さい。

2. : 系が「正」から「負」へ \見掛け" の峠 S (q1 = 0) を通過する場合 系は S (q1 = 0) を通過しても、まだ q1 < 0 の領域では真の峠 S (q1 (p; q) = 0) が多く分布している。その為、系が S (q1 = 0) に戻る確率が高い。この場合、真の峠を通過しきるのに充分な反応自由度方向の入射運動量 p1 (p(0)q(0)) を持っていない限り、系は S (q1 = 0) に必ず戻る。
S[q1(p,q)=0](0,0,+)(1,1,+)(2,2,+)(0,0,-)(1,1,-)(2,2,-)-+q1

S(q1=0)26: saddle I の真の峠と見掛けの峠に対する概念図

図( (p; q) = 0) が配位空間上の \見掛け" の峠 S (q1 = 0)から次第に離れた領域へ拡大してゆく様子が観測されたが、特に多自由度系において \真" の峠と \見掛け" の峠の違いを増幅させる全エネルギー以外の重要な物理量はないだろうか？。 図 27-28 に、\ "エネルギーが増大するに従って、 真 の峠 S q1

(b)

(a)

2.0

2.0
1.0
-0.15 -0.10

0
-0.05

0

q1

図

0.05 0.10
-2.0
0.15

-1.0

p1

1.0
-0.15 -0.10

0
-0.05

0

q1

0.05 0.10
-2.0
0.15

-1.0

p1

27: saddle I における S2nd(q1 ; p1; E )、(a)E =0.1"; (b)E =0.5"

saddle I および saddle II 近傍、E = 0:1; 0:5" 条件下、S (q12nd(p; q) = 0) を (q1; p1) 平面に投影したS2nd (q1 ; p1 ; E ) を示す83 。いずれの鞍点いずれのエネルギー下においても、\分布" S2nd (q1 ; p1 ; E )に対し q1 と p1 の間に正の相関が見られ、遷移状態を通過する速度 p1 が大きいほど \真" の反応83 ここでは、10,000 well-saddle-well (i.e., well1 /well2 ! saddle ! well2 /well1) トラジェクトリーの代わりに、
S (q1 = 0) から出発した 20,000 saddle-well (i.e., saddle ! well1 /well2 ) トラジェクトリーを用いて、交差時のq1 ; p1 を記録してその頻度分布から S2nd (q1 ; p1 ; E ) を評価した。well-saddle-well トラジェクトリーを用いた場合、S (q1 (p; q) = 0) の p1 軸上への投影は、 q1 (p; q) が時間反転に対して対称であるため、至るところ p1 および p1 の両側に本質的に同じである。

63

(a)

(b)
2.0

2.0

1.0
-0.15 -0.10

0
-0.05

0

q1

図

0.05 0.10
-2.0
0.15

-1.0

p1

1.0
-0.15 -0.10

0
-0.05

0

q1

0.05 0.10
-2.0
0.15

-1.0

p1

28: saddle II における S2nd(q1 ; p1; E )、(a)E =0.1"; (b)E =0.5"( (p; q) = 0) は q1 = 0 から p1 の符号によって示される \終状態" の側に多くのボトルネック S q12nd\分布" している。
図 29 は E = 0:5" において、系が S (q12nd (p; q) = 0) を再交差した際の q1 と p1 の数値をプロットしたものである。S (q12nd (p; q) = 0) に対する再交差は q1 の殆んど至るところで生起している
0.10

0.05

p1

saddle I
0.00

saddle II−0.05

naive dividing surface
−0.10
−0.10

−0.05

0.00

0.05

0.10

q1

図

29: saddle I における S2nd(q1 ; p1; E ) 再交差時の q1 と p1 (E =0.5")

0:1 < q1 < 0:05 (saddle I) および 0:1 < q1 < 0:1 (saddle II)）のに対し p1 では jp1j の非常に小さい範囲（ jp1 j  0:1）でのみ生起していることが読み取れる (図 27-28 の x 軸および y 軸の目盛りと比べよ)。このことは、鞍点領域を通過する系の速度（の絶対値）が大きいほど、系が S (
q12nd (p; q) = 0) を \再" 交差する確率は小さいことを意味している。また、saddle I に比べてsaddle II を \登る" ほうが、系は相空間上の峠 S (q12nd (p; q) = 0) をより頻繁に、かつ q1 の広範囲な領域に渡って再交差していることが分かる。このことは saddle II のほうが小さい負の曲率を有すること、ならびに相空間上の反応座標 q12nd (p; q) に沿った運動の規則性が弱いことに起因して（いる。

( = 0) から離

我々は反応のボトルネックが遷移状態を通過する速度の大きさと向きに従って S q1

れた領域に拡大してゆく様相を確認することができたが、全エネルギーに対しての類似の描像が共線型84

H + H2 , H + Cl2, および F + H2 反応に対する２自由度ハミルトン系に対する Pechukas

84 共線型とは 3 原子すべてが一つの軸上に束縛されている反応を指す。

64

らの

1970 年代後半の周期軌道分割面（periodic orbit dividing surface,PODS）の理論研究85 に見ることができる。彼らの手法は遷移状態領域において反応座標に直交する周期軌道が存在しかつ抽出可能であることを前提にしているが、多自由度系の場合、複雑に絡み合った非反応性モード
間の非線形結合による非線型共鳴が存在するため、多くの場合遷移状態を交差するダイナミック
スは多かれ少なかれカオスである。我々が知る限り、この研究は多自由度系（\"の反応のボトルネックを 見た 最初の例である。

\

"

\

Ar6 は 12 自由度）

" (

; )
るルートに対する峠を越える難易度の違い」を解き明かしたが、同時に真の峠 S (q1 (p; q) = 0) を\見る" ことの限界も示している。文献86 に同条件下の S1st (q1; q2 ; E ) の図を載せたが、(q1; q2 ) 平
面への S (q11st (p; q) = 0) および S (q12nd (p; q) = 0) の投影 \分布 は、\反応座標 q1 に沿った反応ダ
イナミックスが明らかに異なっている" のに、1 次と 2 次で殆んど区別することができない。第一
最後に、反応のボトルネックを 見る 限界を指摘したい。投影 分布 S qj ; qk E は確かに「登に、「エネルギー、運動量そして座標に複雑に絡み合っている相空間上の分割超曲面を完全に描
く、もしくは想像することが可能であろうか？」この問いは「多次元の鞍点交差のダイナミック
スの本質を記述するのに充分な少数次元へ縮約可能であろうか？」という普遍的な問いに直結している。
遷移状態領域の階層的規則性 前章で紹介してきたように、遷移状態近傍、活性化障壁に相当す
るエネルギー閾値より高いエネルギー領域において「系の種類に依らず普遍的に」少なくとも次の３種類の階層構造が存在する。すなわち、
擬規則（

Quasi-Regular）領域 ほとんどすべての作用が遷移状態近傍において局所的に保存し、遷移ダイナミックスは規則的である。反応系から生成系へ至る「遷移の過程」は 運動方程式の解析的な解が存在する力学的決定論に完全に従い、遷移状態に至る運動と遷移状態を離れる運動のあいだの動的相関は非常に強く、遷移の次元性は１である（(p; q) 相空間上の反応性モードq1 に対応）。配位空間上に定義された反応分割面を横切る見掛けの再交差運動は相空間上の反応
分割面 S (q1 (p; q) = 0) を横切る一回通過の非再交差運動に常に焼き直すことができ、その相空間上の反応座標に共役な運動量の絶対値が、" しきい値（jp1 (0)j > j!1 q1 (0)j ）" 以上をもつ軌跡はすべて反応系から遷移状態を越えて生成系へと辿る87 。ケテンの高振動励起状態の解離反応ダイ

ナミックスにおいて観測されるように、「量子化された階段的な増加」が存在するエネルギー領域
はこの擬規則領域に対応する。
中間的なカオス（Intermediate, semi-chaotic）領域 エネルギーを更に高くすると、ポテン

シャルエネルギー面の非線形性を系が強く感じるようになり、強い（準）共鳴条件に遭遇するためほとんどすべての作用はもはや保存しなくなり、遷移は擬規則的なものからカオス的なものに

 (p; q) に沿った運動に対する虚の作移行する。しかしながら、少なくとも相空間上の反応座標 q1

用は保存する。これは非反応性モードに対する基準振動数は、複素振動数空間上における虚軸に85 P. Pechukas, E. Pollak, J. Chem. Phys. 67(1977)5976: E. Pollak, in: The Theory of Chemical Reaction

Dynamics, D. C. Clary (Ed.) NATO ASI Series C 170, Reidel, Dordrecht, (1985)p135: E. Pollak, in: Theory ofChemical Reaction Dynamics, M. Baer (Ed.) CRS Press, Boca Raton, Florida, (1986)p123.
86 T. Komatsuzaki and R.S. Berry, J. Mol. Struct.(Theochem) 506, 55 (2000)
87 T. Komatsuzaki and R.S. Berry, J. Phys. Chem. A 106, 10945 (2002)

65

V(q)
Saddle"図"30: 遷移状態領域における階層的規則構造

直交な実軸に帰属されるため、反応性モードがひとつ含まれる場合は如何なるモードの組合わせを考えても（準）共鳴条件は満たされないためである：
M
X
k=1



yn !

k k

 j!1j > O(n )ここで、 y は双曲型モードを必ずひとつ含めた（ n1(265)

6= 0）、任意の整数 nk に対する任意の組み合わせを意味する。この領域では、遷移状態に至る運動と遷移状態を離れる運動のあいだの動的' M 1 となる。虚の振動数 !
 1(p; q) が、J1 (p; q) 同様、遷移過程を通じて近似的に保存する場合、相空間上に定義される反応座標 q1 は強いカオスを呈する非反応性座標から分離され、強いカオス系であるにも関
わらず q1 に沿ったダイナミックスは決定論に従う。すなわち、再交差運動を与えない反応分割面
S (q1 (p; q) = 0) が非常に広いエネルギー領域に渡って抽出可能であることを意味する。このカオ相関は弱くなるが（でも零ではない！）遷移の次元性は q1 の一次元を除く

スに埋もれた局所的な運動不変量は最大局所リヤプノフ指数の解析では「明らかに」検出できな
い。ひとつの可能性としては局所リヤプノフスペクトルを見積もる方法であるが、一次元双曲型
運動に沿ったリヤプノフ指数は可積分系であるにも関わらず正値を取る。これは「局所」リヤプ
ノフ解析の、ある意味、限界と思われる。
完全に発達したカオス（Fully-developed Chaotic）領域（\確率論的" 領域） 更にエネル\

ギーが高くなると、q1 に沿った作用が保存する領域（すなわち、摂動計算の収束半径）が 遷移"
は本質的に \確率論的" 過程に移行する。この領域では遷移状態に至る運動と遷移状態を離れる運に要する距離（範囲） に比べてはるかに小さくなるため、反応系から生成系へ至る遷移の過程

動のあいだに動的相関はなく、遷移の次元性は系を構成する全自由度 M に等しくなる。この領域

66

では運動量空間の情報は反応経路を記述するうえで重要ではなく、再交差運動を与えない反応経路および反応分割面を厳密に定義することは困難になる。
鞍部領域における長寿命滞在：新しいボトルネック.| 多くの化学反応理論88 では、反応の律速を決めるボトルネックは鞍部領域を生成系へ抜けてゆく交差に対応し、そのダイナミックスは、
弾道的（再交差が一切生じない遷移状態理論）もしくは拡散的（ランジュバン的な再交差が生じる

KGH 理論）なもののいずれかであると想定されてきた。これまで見てきたように、エネルギーが増大するにつれ、! 中間カオス領域 ! 完全に発達したカオス領域

擬規則領域の順で、鞍部領域における交差ダイナミックスが「質的に」変化すると予想できるが、交差は弾道的なもの（擬規則、中間カオス領域）から拡散的なもの（完全に発達したカオス領域）へと単調に転移するわけではない。最近、中間カオス領域から完全に発達したカオス領域へ至る途上の、あるエネルギー域において普遍的に、弾道的でも完全に拡散的でもない、長寿命の捕獲を伴う遷
移メカニズムが存在することも明らかになっている89 。

4 大域的な相空間構造
4.1 Davis-Gray のセパラトリックス遷移状態理論

M. J. Davis and S.K. Gray,

J. Chem. Phys.

84, 5389 (1986)

4.2 ホモクリニック交差の構造不安定性

R.E. Gillilan and G.S. Ezra,

J. Chem. Phys.

94, 2648 (1991)

4.3 「カオスのトポロジーが変わる」分岐現象

M. Toda,Phys. Rev. Lett.

74, 2670 (1995);

Phys. Lett. A

227, 232 (1997)

4.4 法双曲的不変多様体とその安定・不安定不変多様体の可視化

T. Uzer, C. Ja e, J. Palacian, P. Yanguas and S. Wiggins,, 15, 957-992 (2002)

Nonlinearity

88 唯一、Kramers 理論の低摩擦極限の表式がエネルギーを供給する媒質と系の相互作用が非常に小さい場合を取り扱っていて、系内部の作用空間（相互作用のない零次では作用は保存する）上での拡散を反応の律速段階と考えている。

89 T. Komatsuzaki and R.S. Berry, Proc. Nat. Acad. Sci. USA 78, 7666 (2001); J. Chem. Phys. 115, 4105

(2001). また、近年の清水寧氏の研究参照。

67

5 これからの化学反応動力学理論の課題
小章ではこれからの反応動力学理論における未解決問題をまとめてみよう。通常、化学反応をおこす系は量子多体系であり、本来、量子カオスの問題に踏み込んだダイナミックスの議論を展開するべきである。しかしながら、

Born-Oppenheimer 近似に基づく古典系に限って化学反応を見渡しても、多くの未解決な問題が山積みの状態である。
「再交差問題」、「部分エルゴード性90 」、「多谷遷移過程における動的相関」、「高次ランクサド
ルを含んだ系における遷移状態概念」、「系と環境、部分と全体」、「階層の異なる相の非隷属性」、「粗視化スケールの協同性」、「タンパク質フォールディング」、「細胞内の化学反応」、。。。思いつくままに、化学反応理論が今後取り組んで行かなければならないキーワードを列挙しただけでもこれだけのものが思い浮かぶ。いくつかについて述べよう。
「再交差問題」：狭義の再交差問題に関していえば、原理的にはほぼ解決した（と思われる）。
サドルの局所領域の相空間構造の解析を通じて、小松崎らは「比較的広範囲なエネルギー領域に\"渡って、化学反応理論における 非再交差仮説 は、観測者の視点を相空間に移すことによって"\

仮説 ではなくなり、殆どの再交差挙動は観測者の視点に依存するのであって、「再交差運動」に\"Kramers-Grote-Hynes 理論は相空間上に展開される遷移状態理論として再解釈することができる。また \摩擦" として振る舞う非反応性自由度（本よって誘起される分子 摩擦 描像に基づく

解説の正の曲率を持つ基準モード）は、「見掛け上」反応の進行を阻害するように振る舞い得るが、実際の遷移はカオスのなかに潜む規則的経路を辿っているのであり、「ものの状態が変わる」\"

ことの 決定性 を問い直す必要がある」を明らかにした。よく量子化学計算で云われる「遷移状態＝鞍部点の構造」という捉え方は、元々のEyring、Wigner らの定義からみても間違っている。

しかし、この捉え方は化学反応を逆に良く理解しているともいえる。すなわち、反応系と生成系をつなぐ再交差を与えない遷移状態はランク１のサドル近傍に局在化し、複雑な化学反応ダイナ
ミックスの過程は「いかに反応系から遷移状態へ至り」、ついで「いかに遷移状態を生成系へ離れ2るか」という つの独立した事象に分割できること、ならびに実験結果とその捉え方からのズレは（古典的には）局所平衡の破れ、もしくは非再交差仮定の破れから解釈できるということを暗黙理にきちんと理解していることを表している。実際、「典型的な化学結合の生成・解離を含む」
化学反応の活性化障壁はおよそ数10kcal/mol である。仮に部分エルゴード性が仮定できる場合、

1 自由度当りの平均エネルギー f は全エネルギー E の一定条件の下、E を総ての自由度の数 n で割った E=n 程度（温度一定条件下では、室温で〜0.6kcal/mol 程度）である：換言すれば、化学反応系の多くは反応「前」と「後」の状態は一様に強いカオスを呈し91 、反応に関わる自由度に活性化障壁以上のエネルギーが集中する確率は極めて小さい。このため、典型的な化学反応系の相空間は反応の「前」と「後」が細いチューブ（＝ボトルネック）によって繋がれ、系は反応「前」90 これは戸田氏（物性研究 74(2000)597）に習って、反応系に対応する相空間の部分においてエルゴード性を仮定

し、生成系については直接は問われていないことを強調するために「部分」という言葉を用いた。
91 実際に反応の「前後」のカオスの強弱を解析した例は少数自由度系以外に殆んどない。ちなみに、（化学結合の生成・解離を含まないが） 3 の異性化反応では、最安定構造付近の反応始状態のダイナミックスは全エネルギーが活性化障壁よりも低い零点振動エネルギーにおいて既に非常に強いカオスを示すことがリヤプノフスペクトル解析か
J. Chem. Phys.
など。
ら知られている。

Ar

T. L. Beck, D. M. Leitner, and R. S. Berry,

68

89, 1681 (1988)

の相空間の部分領域から、そのボトルネックを探し当てるまでに、「部分」エルゴード性が仮定できる程度に隈無く経巡る、と近似できると考えられる。非再交差仮説の根拠は、粗く云えば、遷

(b)

(a)

f

図

31: (a) 化学反応系に典型的な相空間描像, (b) 再交差軌跡の概念図

移状態を越えるために反応方向の自由度に一旦蓄積されたエネルギーが散逸して活性化障壁よりf 程度下がると、その遷移状態を再び交差する確率が小さい（と考えられる）ためである。つま

\

"

)

り、鞍点近傍に非再交差仮説を許すような 分割面 が存在し得る（であろう と考えられる。このような考えがEyring および Wigner 流の遷移状態概念の源流に存在する。上記の「狭義の再交差問題」というのはこのような遷移状態に対して、である。
しかしながら、活性化障壁の高さが１自由度当たりの平均エネルギーに\comparable" である

化学反応系においては、暗黙理に仮定されてきた前提そのものを再び吟味しなければならない。
このような状況においても、再交差を与えない反応分割面（すなわち遷移状態）が存在するか否

2かは決して自明なことではない。 自由度系に限ればDavis と Gray が明らかにしたように答えはYes」であるが、3 自由度系以上の高次元系では、存在する場合と存在しない両方のケースがあり得る（Gillilan と Ezra および Toda らの研究）。しかしながら、どのような分子系および初期条件「の場合に存在／消失するのかは必ずしも明らかではない。このような場合において、もし「真のDavis と Gray のセパラトリックス遷移状態のような相空間

意味の」遷移状態が存在するならば、に渡って非局在化したものであろう。また、こうした状況は狭義の化学反応系では生じない（であろう）が、クラスターの構造転移ならびにタンパク質フォールディングなどの、種々の高さの活性化障壁を越える、多谷遷移の過程においては、普遍的に存在することが容易に予想される。
「部分エルゴード性」： これまで遷移状態領域におけるカオスに埋もれた規則的遷移を議論し\"てきたが、どのように、系は反応系から規則的遷移が埋もれている 領域 に侵入してくるのであろうか？反応系の安定停留点から遷移状態領域の不安定停留点を繋ぐ反応経路の途上で必ずポテ
ンシャルエネルギーの２階微分が零になる変曲点を経由することになり、この領域では調和振動子系を零次に仮定した

CPT の前提そのものが壊れる。遷移の決定論的特性は一旦 \その領域" に侵入した場合は確かに保持されるが、侵入「前」および通過「後」のダイナミックスはまだ手付

Kramers-Grote-Hynes 理論といった統計的反応論の最後に残った未解決問題、反応系における「部分」エルゴード性、すなわち、\遷移

かずの問題として残っている。これは遷移状態理論および状態を通過する前に、系は反応系に対応する相空間の部分をエルゴード的に徘徊してそのボトル

69

"

ネックを探し当てるとする局所平衡 、の問題に深く関連している。高次元カオスと化学反応の1980 年台後半から米国を中心に精力的に行われていたが（例えば、Gray-Davisのセパラトリックス理論）、2 次元系に限定されており、3 次元系以上への拡張は本質的な困難を伴うため、Gillilan-Ezra および戸田らの一部の研究者を除いて頓挫していた。リー変換を高次元
動力学の研究は系鞍部領域に適用した小松崎らの研究は鞍部領域の局所から大域へ繋ぐセパラトリックス理論を多次元に拡張する上で大きな手掛かりを与えている。すなわち、彼らの研究はサドル領域（法双曲的不変多様体）へ近づく（から離れる）安定（不安定）不変多様体を任意の多自由度化学反応( (p; q) = 0) 面上、反応方向に大きい運動量系に対して計算可能であることを示している92 。S q1p1 (p; q) を持っている軌道の時間反転ダイナミックスは、反応を速く進行させるために反応開始

時にどのような初期状態を設定すれば良いかを探る絶好の手段であり、新しい反応制御を創出する可能性を秘めている。
（反応「前」の）高次元相空間構造の一般的な性質については、力学の立場からは、ほとんど理解されていないといってもよいであろう。相空間がエルゴード的でも規則的でもないことを示した特筆すべき実験研究として山内ら93 による高励起振動状態のアセチレンのスペクトルのフラクタル的な階層構造の研究が挙げられる。理論研究としては

3が挙げられるであろう。彼らは 自由度系

Martens、Davis と Ezra94 の解析

OCS 分子のトラジェクトリに沿って局所的な振動数の時系列を解析し、いくつかの振動数の対が有理数比を保ったまま長時間固定されること、すなわ
ち、相空間上の共鳴線上を辿っていること、ならびにある共鳴線と異なる共鳴線の交差領域で長12 自由度系の Ar6 の鞍部交差においても、相空間上の作用変数は保存しない状況下で、振動数 !
k2nd(p; q) が局所的に固定される現象が見出されている
時間相関が生起することなどを見出した。

95 。一般に次元数が増大すればするほど、共鳴を起こすような有理数比の組合せは膨大になり、

3

解析が殆んど不可能になる。本稿で述べたように、 次元以上の高次元トーラスは等エネルギー面を分割することができないため、カオス領域はトポロジー的に繋がっている。それゆえ、粗く云えば、「化学反応にとって重要なこと」は反応の時定数に比べて井戸内のカオス領域を経巡る時
間スケールが短いか否かであって、サドル領域から延びている不安定不変多様体を反応始原系に延ばしていく手法は「化学反応にとって重要な」ポテンシャル井戸内の相空間構造を探る上でも極めて有効な手段といえる。
「多谷遷移過程における動的相関」： 従来の反応理論は単一の（ランク１の）サドルを交差する過程を対象にしているが、これは、反応の進行を記述する自由度が感じるエネルギー障壁が、
例えば kB T に比べて有意に大きく、かつ反応に直接携わらないそれ以外の自由度はエネルギー障
壁を感じないか感じるとしても kB T 程度である系、ないしは、隷属化原理が成立する（反応の自由度とそれ以外の自由度の間に時間スケールの分離が存在する）系では、妥当な描像であると予想される96 が、（局所平衡が各ベイスン間遷移において常に成り立つ特殊な場合を除いて）同程\

度の複数のエネルギー障壁を越えるダイナミックスにおいては一般に破綻する。すなわち、 多谷

92 応用数学者の Stephen Wiggins らは 160CPU のクラスターを用いてこの問題に取り組み始めている。
93 K. Yamanouchi, N. Ikeda, S. Tsuchiya, D.M. Jonas, J.K. Lundberg, G.W. Adamson and R.W. Field, J. Chem.

Phys.

95(1991)6330.

94 C.C. Martens, M.J. Davis, and G.S. Ezra, Chem. Phys. Lett. 142, 519(1987)
95 T. Komatsuzaki and R.S. Berry, J. Chem. Phys. 105(1999)10838
96 統計性の問題や反応座標の定義の問題は別として。

70

"

\

"

1

性 ならびに 動的相関 を問う必要が生じる。本稿でみてきたように、ランク のサドルの「前」と「後」を繋ぐダイナミックスの間には、高エネルギー条件下にあっても、動的相関が顕著に存在することが原理上明らかであるが、それらがどのように保持し得るかは全くの未解決である。

2

このような動的相関が生じ得るダイナミックスに関する研究は恐らく 通りの戦略があり得る。


図32: 多遷移化学反応ダイナミックス

1. 動的相関が発生するための相空間構造の基礎論 リー変換正準変換摂動論をある鞍部点近傍に展開し、近似的に法双曲的不変多様体
隣接する別の鞍部点近傍の(NHIM) を抽出し NHIM から離れる不安定不変多様体と、

NHIM に入ってゆく安定不変多様体との間の多次元相空間上の交差現象を解析することである。これは非常に魅力的な制御問題を含んでいる。すなわち、安定／不安定多様体の次元はともに(2n 2) であるから、交差の次元は一般に (2n 3) 以下である。交差に関与する自由度ならびに交差に関与しない自由度が仮に存在するならば、法双曲的多様体の初期条件（交差に関与する自由度だけを励起）に応じて、交差の次元、すなわち、動的相関の強さを制御できる可能性がある。

2. ある大域的な集団座標を仮定し、その集団座標で張られる状態空間構造を解析する 凸凹なエネルギーランドスケープを伴うタンパク質のように、時定数が拮抗する自由度が多数カップル
して存在し、かつ、広い配位空間を調べる必要がある場合、なんらかのコヒーレンスが存在しているとしても「その集団座標をどのように切り出すか」という問題に対する一般的な解答はない。
ここでは、主成分解析と呼ばれる統計学の手法を用いたタンパク質ダイナミックスの研究例をいくつか紹介しよう、主成分解析とは、観測値の分散が最大になるように、観測値を線形変換した
方向を第一主成分座標とする（第二、第三主成分座標になるほど分散が順次小さくなる）もので、北尾らは小さな無数の起伏を有するベイスンの「内」の運動に対する第一主成分座標は単純な拡散運動ではなく、むしろコヒーレントな弾道運動であることを示した。これはまさに動的相関のGarcia と Hummer97 は大小のベイスンの「内」から「外」への遷移を伴うダイナミックスの解析を行った。彼らは水中の酸化シトクロム c 全原子 MD シュ
存在を端的に示した最初の例である。

ミレーションを敢行し蛋白質部位の主成分座標 Xi を抽出し、その平均自乗変位の自己相関関数

97 A.E. Garcia and G. Hummer, Proteins 36, 175 (1999)71i (t)(= h(Xi (t) hXi i)2 i) を各温度で見積もった。天然状態が安定に存在する 300K での第一、第二主成分の 1 (t) と 2 (t) を図 33 に示す。i (t)  t ( ' 1) であれば、観測量 Xi (t) がマルコフ過程を意味するブラウン運動に従うことを意味するが、これらの集団座標は明らかにブラウン運動

100 ピコ秒（ps）以下では ' 0:5 であり、（小さな）ベイスンの「内」の大振幅集団運動は通常のブラウン運動よりも抑制された \拡散"（sub di usion）であること；一方、100ps 以上の時間領域では ' 1:5 であり、大振幅集団運動は通常のブラウン運動よりも促進された協奏的 \拡散"（super di usion）であることを示している。これは（大きな）ベイスン
ではない。すなわち、

図33: 300Kにおける X1

X2 平面上の軌跡（左図）：1 (t)-t と 2 (t)-t の両対数プロット（右図）の「内」から「外」へ越える遷移の過程が非ブラウン的な弾道運動であることを反映しているも100ps 以上において観測された非ブラウン的な弾道運動が、より大きなベイスン間遷移を伴う高い温度 360K および 430K においても同様に観測されること；ns の時間領域で unfolding が観測されはじめる 500K 以上になると、これらの集団座標においてもブラウン運動（ ' 1）に漸近することを明らかにした。松永ら98 は最近カオス時系列解析の埋め込みのと考えられる。彼らは、論と呼ばれる手法をタンパク質フォールディングのダイナミックス解析に応用し、転移温度では他の温度領域に比べて、揺らぎの大きい数十成分の埋め込み次元が顕著に小さく、ある自由度に対しては動的相関の存在を支持する結果を見出している99 。

高次ランクサドルを含んだ系における遷移状態概念 負の曲率のモードが１つある、第一ランクの鞍点のポテンシャルエネルギーよりも系の全エネルギーが遥かに大きくなると、系はより広い
ポテンシャルエネルギー曲面の地形・形態を見渡せるようになる。この場合、（負の曲率を持つモードの数が１以上の）より高いランクの鞍点も経巡る可能性が生じ、ひとつの平衡状態から別の平衡状態へ転移する過程は必ずしも両者を繋ぐ通常の峠を通らずに、複数の方向への分岐を伴\ "う（多次元ポテンシャルエネルギー面固有の） 峠 を経由することも充分考えうる。系を構成す\"

るすべての自由度の数だけ分岐の方向がある場合がこの多次元面の 山頂 と見做すことができ

98 Y. Matsunaga, K.S. Kostov and T. Komatsuzaki, J. Phys. Chem. A 106,10898 (2002)
99 埋め込み論の最大の利点は任意の時系列データに関して有効である点であるが、単に次元を解析するだけでなく、メカニズムを抽出するための理論的道具の開発も必要となる。例えば、保存力学系に限って言えば、物理的な意味が不明瞭な遅延座標系は、リー変換により元の座標系の「なに」を見ていることに対応するかを抜出すことが原理上可
能である。

72

る。また、通常の峠を交差する場合に限っても高エネルギーの条件下では、遷移状態領域の局所的な相空間体積が膨張するため、「部分」エルゴード性が成り立つ前に系は遷移してしまい、更には、次のベイスン領域に充分留まることなく、また別の峠を次々と越えてゆく多遷移現象が誘起されるであろう。このように、高エネルギー下の化学反応は多遷移・多チャンネルな「動力学的」な過程であって、従来の「統計的」反応理論の範疇を遥かに越えている。例えば、鞍部領域にお\"反応は（有限多体系の）\相転移" や蛋白質の折れ畳みダイナミックスに普遍的に内在しており、
ける ! の虚空間における ある短時間 の共鳴と分岐に関する未開拓問題も生じる。このような

そこでは始状態と終状態を二分する反応分割面、すなわち、（単純な意味での）遷移状態概念そのものが問われることになる。また、多遷移現象は「反応の自由度とそれ以外」という描像が時間に対してもはや不変ではなくなるため、反応論の文脈における「部分と全体」の境界の生成・崩壊のダイナミックスの問題も包含している。このような未開拓問題に関する近年の幾つかの試み（

Nayak ら、志田、世古・高塚）に関しては、総説100 に紹介したのでそちらを参照して頂きたい。

6 付録
6.1 (39)-(40) 式の証明
無摂動ハミルトニアン H0 が調和振動子系の場合、角変数

 に依存しない H に従う (p; q) の運動方程式は以下のようにして解析的に解くことができる。すなわち、

dpk
dt

dqk
dt

=
=

@ H (J)@ H (J) @ Jk=@ qk@ Jk @ qk! k (J)!k qk ;

(266)
(267)


 
= @ H@ p(J) = @ H@ J(J) @@Jpkkkk= !k (J) p ;

(268)
(269)

k

!k

(269) 式を時間 t で微分した式に (267) 式を代入すると
d2 qk ! k (J) dpk
= ! dt = !k2(J)qk :
(270)
dt2
k
ここで (267),(269) 式の 3 番目の等号は (34)-(39) 式および
1 p dq = p2k + !k2qk2
Jk =
(271)
2 k k2!kの関係から導かれる。正準変換摂動理論では新しい作用 Jk の (p; q) に関する \関数の形" は非摂動ハミルトニアンにおける作用 Jk の (p; q) に関する関数の形を保存する。
 ) = H 0(p; q) = H0(p; q) = H0(J;  )H 0 (J; (272)= 12 (p2k + !k2qk2) = !k Jk :(273)
I100 小松崎民樹、物性研究 76(2001)1.

X

X

k

k

73

（実際に

(41)-(42) 式を (271) 式に代入すると (271) 式の 2 番目の等号が成り立っていることを確認できる）

6.2 リー変換
ほとんどすべての解析力学の教本では、新旧の正準座標が混ざった母関数を用いて正準変換を説明しているが、リー母関数 W を使った正準変換のほうが直感的に理解しやすく、また、正準変
換摂動理論を高次に展開してゆくうえで圧倒的に強力かつ必要不可欠である。

 および W を \時間" および \（ある正準座標 q とそれに共役な運動量 p の関数として表され"る）ハミルトニアン とする次式の正準方程式を考える。

dpi
@W
=
;
d
@ qi

dqi @W= ;d @ pi(pi; qi) をまとめて zi として表記（zi = (pi; qi)）すると、
dzi= fzi; W (z)g  LW zi :dとなる。ここで

(274)
(275)
(276)

f g は Poisson 括弧を意味し、任意の微分可能な関数 u; v および w ならびに任意の正準座標・運動量の組に対する偏微分操作に対して

fu; vg  ( @@uq @@vp @@vq @@up );
i
i
i
i
i
X

(277)

により定義され、次の関係式が成り立つ。

fu; vg = fv; ug;(278)fu; v + wg = fu; vg + fu; wg;
fu; vwg = fu; vgw + vfu; wg;
fu; fv; wgg + fv; fw; ugg + fw; fu; vgg = 0:

(279)
(280)
(281)

任意の正準座標系（例えば、(p; q)）に対して

fqi; pj g = Æij ; fqi ; qj g = fpi; pj g = 0（Æij は(282)

Kronecker デルタ）が成り立ち、逆に、(282) 式が成り立つ座標系は正準であることが証明される（詳細は力学の教科書を参考のこと）。

(276) 式で導入された演算子 LW ( fW; g) は次の性質をもつことが (278)-(281) 式がら容易に確認することができる：

LW ( u + v ) = LW u + LW v;

74

(283)

LW uv = uLW v + vLW u;

(284)
(285)
(286)

LW fu; v g = fu; LW v g + fLW u; v g;（,

LV LW = LfV;W g + LW LV ;および V は任意の数および任意の微分可能な関数である）。

子 LW を n 回操作すると(283)-(285) 式に対して、演算

LnW ( u + v ) = LnW u + LnW v;
LnW uv =

n
X

m=0
n
X

LnW fu; v g =が成り立つので、(

mn Cm LW v

m=0)(LnW m u);

(287)
(288)

g

(289)

f

m
n m
n Cm LW v; LW u ;e L uv = e L u e L v ;




W

W

e L fu; v g = fe L u; e L v g
W

\

(290)
(291)



W

W

が得られる。さて、以下に ハミルトニアン

W" W が \時間" に陽に依存しない \自律系" を取り上げて、リー変換による正準変換を解説してゆこう101 。

\自律系" に対する (276) 式の形式解は、簡単のために、発展演算子 T ( e L ) を導入すると、Wz() = e L z(0)  T z(0)(292)
と表される。この形式で表現されるあらゆる変換は、もし z(0) が正準であれば、z() も正準であり、逆もしかりであることが以下のようにして証明することができる：\時刻"0 における z を(p; q) として、\時刻" における z を (p; q) と便宜上表記する。
pi = T pi = e L p q pi = pi (p; q; );(293)qj = T qi = e L p q qj = qi (p; q; );
(294)
W

W(

および

3

2

;

)

W( ;

)

fqi; pj g = fT qi; T pj g = T fqi; pj g;
fqi; qj g = fT qi; T qj g = T fqi; qj g;
fpi; pj g = fT pi; T pj g = T fpi; pj g;

ここで、これら つの式の 番目の等号は(291) 式の性質を使った。

(295)
(296)
(297)

fqi; pj g = T Æij = e L Æij = Æij ;(298)

fqi; qj g = fpi; pj g = T  0 = 0

(299)

W

101 非自律系についても最終的に得られる式は同じ形をしている10275

(p; q) が正準であれば (p; q) も正準であることが証明された。一般に pi および qi への LW (p;q) の畳み込みは (p; q) の複雑な非線形関数を生成する。(293),(294) 式における pi およびqi はハミルトニアン W に従う時間発展に沿った \（経過した）時間"  と時刻 0 での (p; q) の初期値の関数として表されることを意味している。すなわち、pi (p; q; 0) = pi , qi (p; q; 0) = qi :
 から 0 への時間反転を考えると、
より、もしz(0) = T 1z();もしくはpi = T 1 pi = eLqi = T 1 qi = eL

p
;q)

W(

p
;q)

W(

(300)

pi = pi (p; q; );

(301)
(302)

qi = qi (p; q; );= T 1T = 1 より、(295)-(296) 式に時間反転演算子 T 1 を前から演算すると「もし (p; q) が正準であれば (p; q) も正準である」ことが容易に証明される。
図 34 に図解すると、発展演算子 T を z(0) に演算するという T z(0) 操作により、(p; q) 座標系の; q) 座標系の点 B へ移される。これは \ハミルトニアン" W (p; q) によって \時
任意の点 A は (p刻" 0 から  まで (pA ; qA ) を \時間発展" させたことに相当する。一方、T 1 z() により、 (p; q) の（任意の）点 B は (p; q) の点 A へ移される。これは \ハミルトニアン" W (p; q) によって \時刻" から 0 に (pB ; qB ) を時間反転したことに相当する。ここで、W (p; q) は W (p; q) の \引き数"pi お
よび qi を pi および qi に名称変更しただけで、関数の形は両者で保存している。
と表される。T T 1

p

p

"B"
x

x

t=ε

A t=0

B
x

z(ε)=Tz(0)qqpq

Hamiltonian dynamics by W( , )z(ε)=(p,q)z(0)=(p,q)
p

p

B
x

t=ε

z(0)=T z(ε)
-1

x

x

A

q

"A" t=0q
pq

Hamiltonian dynamics by W( , )図34: リー変換による (p; q) の正準変換

では、リー変換は関数に対してどのように働くのか？

76微分可能な任意の関数 f(p; q) に対するリー変換 T f (p; q) は p、q および  で記述される、新しい非線形関数（g と以後表記する）を生成する。

T f (p; q) = e L f (p; q) =W

n
( LW )nf (p; q)  g(p; q; ):
n
!
n=0
X

(303)

g (p; q; ) はどのような意味をもつのであろうか？LW f は
LW f = fW; f g =

=
=
更に、

L2 f = L(

W LW fW!

@W @f @W @f@ qi @ pi @ pi @ qii!
Xdpi @f dqi @f+d @ pi d @ qi
i
df
;
d

X

df
d

) = LW!

2= ddf2 ;

(304)
(305)
(306)
(307)

となることから、n 回 LW を f に演算すると( LW )nf = ddfn :(308)

X n dn fn(
LW (z(0)) )n f (z(0)) =nn=0 n! d =0n=0 n!

(309)nとなる。したがって、新しい関数 g(z(0); ) は

g (z(0); ) = T f (z(0)) =Xのように表すことができる。最後の項は関数 f を原点近傍  でテーラー展開したものであり、f は

 に対して z を介して陰に依存しているので、

g (z(0); ) = f (z()) = f (T z(0))(310)(z(0); ) は（ \ハミルトニアン" W に従って時間発展する時刻  での）「点」z() における f の関数「値」を「初期点」z(0) と「経過時間」 の関数として表している。関数 f へのリー変換およびその逆変換を図解すると（図 35）、
T f (p; q) は (p; q) 座標系の任意の点 A, (pA ; qA ), で評価される関数 f を「同じ」点 A で評価される新しい関数 g (p; q; ) に変換し、その関数「値」は (p; q) 座標系の点 B で評価される f (pB; qB)と等しい。点 B は (T pA ; T qA ) に等しく、\ハミルトニアン"W (p; q) によって A 点を \" だけ時間発展した点である。逆に、T 1 f (p; q) は (p; q) の点 B で評価される関数 f を「同じ」点 B で評価される新しい関数 h(p; q; ) に変換し、その関数「値」は (p; q) 座標系の点 A の f (p; q) のとなる。つまり、新しい関数 g値と等しい。
関数 f へのリー変換をまとめると、

T f (z(0)) = f (T z(0)) = f (z()) = g (z(0); )

77

(311)

f

f

A

f( p ,q )Tf(p,q)f( p ,q )
A

B

B

p

p

x

A

t=0

q

"B"
x

t=ε

f

x

q
f

f( p ,q )B

B

T f( p ,q )-1f( p ,q )
A

B

A

p

p

x

x

A

"A"

t=0

q

図
ないしは

t=ε

q

x

B

35: リー変換による関数 f の変換

T f (p; q) = f (T p; T q) = f (p; q) = g (p; q; ) ;(312)

および逆リー変換をまとめると、T 1 f (z()) = f T 1 z() = f (z(0)) = h (z(); ) ;(313)

T 1 f (p; q) = f T 1 p; T 1 q = f (p; q) = h (p; q; ) :(314)

W (p; q) = W (e L p; e L q)= e L W (p; q)= W (p; q)

(315)
(316)
(317)



ないしは







(p; q) で評価される、ある任意の微分可能（原理的には C 1 級）な非線形関数 f (p; q) に対して、リー変換は、すべての成分が pi = T pi および qi = T qi に従う (p; q) で評価される f (p; q) の関数「値」を与える新しい関数 g (p; q; ) を生成する。例えば、f が \ハミルトニアン"W の場合、 n = 0
除いて LnW W = 0 より
W

W

\

"

W

となる。すなわち 時間  に陽に依存しない W によるリー変換は時間発展に沿って W の関数形は保存される。

78

あとで必要になるので、発展演算子 T と LW の関係式を導出しておく。 について T ff (p; q) を微分すると、(p; q) =X
dT ()@f (p ; q) @ qi @f (p ; q) @ pif (p; q) =+ @ p @d@ qi @ii!
X@f (p ; q) @W @f (p; q) @W=@ qi @ pi@ pi @ qii= ff (p; q); W (p; q)g
= fT ()f (p; q); T ()W (p; q)g= T ()ff (p; q); W (p; q)g= T ()LW f (p; q):!

これはあらゆる任意の f同様に、 T T 1(p; q) に対して成り立つのでdT= T LW :
d

(318)

dT 1
=
LW T 1 ;d(319)= 1 を  について微分するととなる。なぜならば、

dT
T 1+T
0 = dT ()dT () = dT
d
d
1

=

T LW T 1 + Tとなり、最後に得られた式に対して、前から T 1 を掛けると

1

dT 1
;
d

(319) 式が得られる。

\ハミルトニアン"W が \時間"  に陽に依存する非自律系 W (p; q; ) におけるリー変換は自律
; q) で得られた性質を保存する。すなわち、
系 W (p

および任意の関数 fpi = T pi = expqj = T qi = expZ f (p; q) = exp

Z 

i

(320)

LW (p;q; )d0 q = q (p; q; );

(321)

0

(p; q) に対して


LW (p;q; )0Z LW (p;q; )0d0 p = p (p; q; );
i



i

i



d0 f (p; q) = T f (p; q) = g (p; q; );が成立し、非自律系のリー変換をまとめると、

T f (z(0)) = f (T z(0)) = f (z()) = g (z(0); )T f (p; q) = f (T p; T q) = f (p; q) = g (p; q; ) ;

79

(322)

および逆リー変換をまとめると、T 1 f (z()) = f T 1 z() = f (z(0)) = h (z(); ) ;T 1 f (p; q) = f T 1 p; T 1 q = f (p; q) = h (p; q; ) :となる。これはハミルトンの正準方程式が自律／非自律系を問わずに成立することから容易に理解することができる。図

"34 および 35 においても、非自律系の場合は \ハミルトニアン" W が \時

間  に陽に依存するため、仮想的時間発展に沿ってその形が一般に変化し得ることを除いて、同じ解釈が成り立つ。

6.3 リ─正準変換摂動計算 (LCPT) の応用例

2 自由度ハミルトン系
1
H (p; q) = (p21 + p22 + !12 q12 + !22 q22 ) q1 q22 ;2= (! J + ! J )  2J1 2J2 cos  cos2  ;
1



1 1

に対して
は

2 2

2

!1

!2

1

2

(323)
(324)

LCPT を実際に計算してみよう。LCPT の処方箋に従って、新しいハミルトニアン H 0(p; q)O(0 ) において、H 0 (p; q) = H0 (p; q)

(325)
(326)

= 21 (p21 + p22 + !12q12 + !22q22) = !1J1 + !2J2
 0 (p; q) に従う運動方程式の一般解はによって与えられる。 H
dpk@ H 0 dqk @ H 0=;
= @ p ;d@ qk dkqk ( ) =pk ( ) =である。さて、まず H1

v
u
u
t

2Jk

(327)

 k = 2Jk cos(!k  + k );
cos

!k
!k
2!k Jk sin  k = 2!k Jk sin(!k  + k ) (k = 1; 2)
v
u
u
t

q

q

(J;  ) は

 ) = 1 2J1 2J2 
H1 (J; 
4 !1 !2
2 cos  1 + cos( 1 + 2 2) + cos( 1 2 2)

(328)
(329)

!1

2





(330)
となるが、平均化の原理に従って、これを時間に顕に依存しない非振動部分 hH1 i と振動部分 fH1 gに分けて、新しいハミルトニアンをできる限り、「見通しのよい」形式に変換することにしよう。

80

1. 非共鳴の場合 : n1!1 + n2!2 6= 0

ここでは、!1 と !2 が有理数上一次独立ではない、すなわち、n1 !1+ n2!2 6= 0 (n1; n2 : 任意の

0 でない整数）、場合を扱う。(330) 式からH 1 = 0;W1 =fH1gd = 14 2!J1 2!J2 [ !2 sin  1
1
2 1
+ ! +1 2! sin( 1 + 2 2) + ! 1 2! sin( 1 2 2)];
1
2
1
2
!1

Z

もしくは、

2

(p; q) を用いるとW1 =ゆえに、H1

(2!22

(331)

!12 )p1 q22 + 2!12q1 q2 p2 + 2p1 p22:
!12 (!1 + 2!2 )(!1 2!2 )(332)(p; q) = q1 q22 を用いると
2
2 q 4 + 2
q22 p22 + 4!12q12 q22 + 8q1 p1 q2 p22
fW1; H1g = (2!2 !1 )
;!12 (!1 + 2!2 )(!1 2!2 )(333)(J;  ) 表記では
fW1; H1g = !2(! + 2!4)(! 2! ) [
2 1
2
1 1
1 !1 2 J2 (cos 4 + 4 cos 2 + 3) + J2 (1 + cos 2 )2
2
2
2
8 !2 2+ !!1 J1J2 1 + cos 2 1 + cos 2 2
2+ 12 cos 2( 1 +  2) + cos 2( 1  2)
 1 +  2) cos 2( 1  2))J1 J2 (cos 2(はとなる。O (2 ) までの次数で新しいハミルトニアン H
H = 0 (!1 J1 + !2 J2 )(334)
2+ 2 !2(! + 2!2)(! 2! ) 1 38 !!1 J22 + !!1 J1J2 :
(335)
2 1
2
2
2
1 1
 (p; q) に従う運動方程式ならびにその解はによって与えることができる。この 2 次までの H
dpk@ H dqk @ H
;
(336)
=
= ;
dt
@ qk dt @ pk










i

"

qk (t) =pk (t) =

v
u
u
t



 !

2Jk cos(! (J)t + );kk!k
2!k Jk sin(!k (J)t + k );

q

81

#

(337)
(338)

で与えられる。ここで、 k はモード k の初期位相であり、振動数は! 1 (J) = !1 + 2

2J2;
!1 !2 (!1 + 2!2 )(!1 2!2 )3 !1 2 J + !1 J :
2
2 ! 1
2!2)
4 !2
2

2
!2 + 2 2!1 (!1 + 2!2 )(!1! 2 (J) =

(339)

"



 !

#

(340)

で与えられる。リ─変換による関数変換はf (p; q) = T f (p; q)= 0f (p; q)= f(p; q)

1 fW1 ; f g +2

2 (fW1; fW1; f gg fW2; f g) + : : :(p; q) で表すと(!2 2!22)q22 2p22 + O(2 );1 2 1!1 (!1 + 2!2 )(!1 2!2 )2q2 p22
1(! + 2! )(! 2! ) + O( );2であるから、 次までの範囲で新しい座標、運動量、作用をq1 (p; q) = q1
p1 (p; q) = p1
J1 (p; q) =1

p21 + !12 q12



1



2!1
となり、!
1(p; q) を O(3 ) まで求めると

(!2
1 1

2

1

(341)
(342)

2

2!22)q1q22 2q1 p22 + 2p1q2 p2 + O(2);!1 (!1 + 2!2 )(!1 2!2 )2! 1 (p; q) = !1 + 2



!1 !2 (!1 + 2!2 )(!1 2!2 )"
1 p2 + !2q2 + 1 2!22q1 q22 2q1p22 + 2p1q2 p2 # :2!2 2 2 2
!2 (!1 + 2!2 )(!1 2!2 )f(p; q) を逆に (p; q) の関数として変換する逆リ─変換f(p; q) = T 1 f(p ; q)2= 0f(p; q) + 1 fW1; fg + 2 fW1; fW1; fgg + fW2 ; fg + : : := f (p; q)を用いると、例えば q1

(343)

(344)



(p; q) は(!12 2!22)q22 2p22
1
+
O
(
2 );
q1 (p; q) = q1 +  2!1 (!1 + 2!2 )(!1 2!2 )と直ちに書き下すこともできる。

82

(345)

2. （準）共鳴の場合: n1!1 + n2!2 ' 0 もし !1 2!2 が非常に小さい、すなわち、(346)
;  ) の式は発散することが分かる。このような場合は、(330) 式の H1(J;  ) の右の場合は W1 (J 1 に組み入れることで発散を回避することができ
辺の第 3 項を  に依存しない非振動項として H!1

2!2  O(1 )る。すなわち、

1 2J1 2J2 cos( 2 )H 1 =
1
2
4 !1 ! 2W1 =(H1 H 1)d
= 14 2!J1 2!J2 !2 sin  1 + ! +1 2! sin( 1 + 2 2)
1
2
1
1
2
!1

2

Z

!1
2

（詳細は堀による原著103

) 小さな分母の問題から発散する項を新しいハミルトニアンの一部として組み込んだが、厳密には共鳴ではないが、無限の発散ではないが有限の大きい値を産み出す準共鳴の場合、原理上は収束し得るが摂動級数を無限次まで展開する必要が生じる。

6.4 相空間ヤコビ行列の行列要素(p; q) のべきで表現可能な場合はリー変換を使えば n 次元相空間ヤコビ行列の行列要素も極めて容易に任意の次数まで算出することができる。例えば、1 自由度系の相空間
ヤコビ行列 M のトレース trM は、一般に、
ハミルトニアンがH (pt ;qt ) p

trM = @p@pt+ + @q@qt+ = @e @p
t

t

t

t

H (pt ;qt ) q+ @e @q

t

t

(347)

と書くことができる。ポテンシャルエネルギーの任意の停留点近傍では、ハミルトニアン H はH=

1 p2 + !2q2
2




(348)

と近似することができるので、例えば、 2 までリー変換を施すと、

trM = 2 !2 2 + O( 3)(349)となる。すなわちサドル点近傍（! 2 < 0）ならびにポテンシャル極小点近傍（! 2 > 0）では、確
かに jtrM j > 2 ならびに jtrM j < 2 となることが分かる。
7 おわりに
現在、化学反応とカオスに関する書籍は勿論、化学反応を理解する上で必要なカオスの基礎知
識に関して平易に説明されているものも見当たらない。本稿は、近年進展著しい「化学反応動力103 G. Hori, Pub. Astro. Soc. Japan 19, 229 (1967) を参照のこと。

83

学とカオス」の研究の一端を単に説明するだけではなく、将来の分子科学を背負って立つべき皆さんが進展著しい

ab initio MD 法などの分子動力学シミュレーションを通して、「カオス力学屋

からだけでは生まれてこない」分子科学・分子論的な視点からの新しい動的反応描像を描いてくPoincare-Birko の不動点定理、大れる手掛かりになることも願っている。時間と紙面の制約で、域的相空間構造に基づく化学反応動力学など書くことができなかったが、夏の学校当日の皆さんにお願いする文献紹介を通じて皆さんと一緒に考えていきたい。
本稿の作成に当って、神戸大学理学部物理学科の下野昌宣君
究科の堀田浩司君(B4)、東京大学大学院総合文化研(D2) に図 1 の作成（下野君）、原稿の査読などで大変お世話になりました。ここに感謝します。

8 参考文献
化学反応理論に関しては、「化学反応論」笛野 高之 著（朝倉書店,1975）：日本の研究者に依って書かれた化学反応理論に関する数少ない貴重な良書であるが、残念ながら絶版。この他、
「化学
動力学」

Francisco and Hase 著、佐藤 伸訳（東京化学同人,1995）

大学院講義 物理化学 近藤保編、小谷正博・幸田清一郎・染田清彦著（東京化学同人,1997）が、Grote-Hynes 理論に関してはJ.T. Hynes, (1985) in Theory of Chemical Reaction Dynamics ed. Baer, M. (CRC, Boca Raton,
FL), 171-234. を挙げておく。

まずまず纏まっている。また、

カオス力学系と化学反応理論に関しては、
戸田幹人、「化学反応の動力学とカオス」物性研究

74(2000)597.

小松崎民樹、「化学反応の「揺らぐ」世界における力学的決定性| 70 年来の「非再交差仮説」|」物性研究 76(2001)1.
T. Komatsuzaki and R. S. Berry, \Chemical Reaction Dynamics:Many-Body Chaos and Regularity"
123,79(2002).
M. Toda, \Dynamics of Chemcial Reactions and Chaos"
123,153(2002).
T. Uzer, C. Ja e, J. Palacian, P. Yanguas and S. Wiggins, \The Geometry of Reaction Dynamics", 15, 957-992 (2002).

の解決

Adv. Chem. Phys.

Adv. Chem. Phys.

Nonlinearity

カオス力学系の入門書としては

Michael Tabor \Chaos and Integrability in Nonlinear Dynamics" John Wiley & Sons, NY.
(1988)
じっくり読むには

A. J. Lichtenberg and M. A. Lieberman,
Verlag, New York, 1992).

Regular and Chaotic Dynamics 2nd Ed.

がよいであろう。和書では現代の物理学１「力学」、大貫義郎・吉田春夫著 第理論を解説している。

84(Springer-

II 部 (岩波講座,1994) がリー正準変換摂動

9 問題集
文献紹介リスト[1] M. J. Davis and S.K. Gray,84, 5389-5411 (1986)[2] R.E. Gillilan and G.S. Ezra,
94, 2648-2668 (1991)[3] M. Toda,
74, 2670-2673 (1995) (参考資料 M. Toda
123,153-198(2002))[4] T. Uzer, C. Ja e, J. Palacian, P. Yanguas and S. Wiggins,, 15, 957-992 (2002)[5] T. Komatsuzaki and R.S. Berry,
106 10945-10950 (2002)[6] C. Amitrano and R. S. Berry,47, 3158-3173 (1993)[7] D. Sugny, M. Joyeux and E.L. Sibert III
113, 7165-7177 (2000)
1) 6.3 章 「（準）共鳴の場合」を参考に、Henon-Heiles の 2 自由度ハミルトン系
11 q3
1H (p; q) = (p21 + p22 ) + (q12 + q22 ) +  q12 q22232
 を (J;  ) およびのリ─正準変換摂動計算を  の 2 次まで計算し、新しいハミルトニアン H(p; q) の関数として求めなさい。新しい作用 J1; J2 は H に従う時間 t に対して変化するが、J. Chem. Phys.

J. Chem. Phys.

Phys. Rev. Lett.

Adv.

Chem.

Phys.

Nonlinearity

J. Phys. Chem. A

Phys. Rev. E

J. Chem. Phys.



H0 は t に対して不変量になっていることを確認しなさい。

2)323）式の 2 自由度ハミルトン系における p1; q1 ; p2; q2 を新しい座標 pk とそれに共役な運動量 pk の関数として  の 2 次まで求めなさい。p2 = q2 = 0 とおき、q1 ,
p1 に（337）、（338）（

式を代入したときの p1 ; q1 ; p2 ; q2 の振る舞いを示しなさい。

3)17（右）の turnstile の A と B のローブ領域の面積は等しいことを示しなさい。同様に、ヘテロクリニック交差がつくる turnstile の 2 つのローブの面積も等しいことを示しなさい。
4) 正準変換 (p; q) ! (p; q) とは、
fqi; pj g(p;q) = Æij ; fqi; qj g(p;q) = fpi; pj g(p;q) = 0(350)
を満足する座標変換である。ここで、下付きの添字 (p; q) は Poisson 括弧は「(282) 式が成り立つ」正準座標系 (p; q) で計算することを強調するために用いた。このとき、正準変数の図

任意の（微分可能な）関数 A; B に対して

fA; B g(p;q) = fA; B g(p;q)(351)

Poisson 括弧は正準変換のもとで不変であり、添字は落せること）、であること（すなわち、
ならびに(p; q) も (p; q) 同様にハミルトンの正準方程式に従うことを示しなさい。
85

5)

6.2 章の（283）-（291）式が成り立つことを確認しなさい。

6)（25）式の Poisson 括弧の関係式が（260）および（261）式で定義される反応座標の作用および角変数を含めて、成り立っていることを確認しなさい。

7)

2 重井戸型ポテンシャル
H=
および

H=

1 p2
2 1


1 p2 + p2
2 1 2

q12 + q14 （1 自由度系）


q12 + q22 + q14 + q12 q2 （2 自由度系）



(
不安定多様体の描き方については M.J.Davis,

)

(352)
(353)

における安定多様体および不安定多様体を p1 ; q1 断面上に図示しなさい。（安定多様体と
J. Chem. Phys.

, 83,1016 (1985) の Appendix

を参照すること）可能であれば、片方の井戸の内部からもう片方の井戸へ遷移する交差点\"

列、どちらの井戸内部にも到達しないで 回転しつづける 交差点列、ずっと井戸の内部に留まる交差点列などを探しなさい。

8)

軌道不安定性を評価する指標として最大リヤプノフ数・リヤプノフスペクトルと呼ばれるものがある。その計算方法を概説するとともに、化学反応ダイナミックスを解析するうえでの(利点・欠点を論じなさい。 参考資料

96, 1376 (1992))

R.J. Hinde, R. S. Berry and D.J. Wales,

表

日

8/7

時間
午前

8/8

全体講演
午後
夜
午前
夜
午前
午後

8/9

3h

3h
2h
3h
2h
3h
3h

,

J.Chem.Phys.

3: 分科会予定

内容
講義（予稿で分からないところのおさらい）
正準変換摂動理論 ホモクリニック交差
化学反応論の理論、、蛋白質フォールディング
問題集の解答
論文紹介、
講義（化学反応とカオスの解説）
論文紹介
、
論文紹介
、
論文紹介
講義 化学反応動力学理論の向かうべき方向
高次ランクサドルの役割、 ボトルネックの存在、動的相関

Gaspard-Rice
Overview, Grote-Hynes
Zwan-Hynes

Marc Joyeux Amitrano-Berry(

Komatsuzaki-Berry Davis-Gray
Gillilan-Ezra Uzer&Wiggins
Toda86)