化学反応系の微分幾何学：Tachibana-Fukuiによる『Meta-IRC』と拡張ヘシアンの定式化

4440 words

22 minutes

化学反応系の微分幾何学：Tachibana-Fukuiによる『Meta-IRC』と拡張ヘシアンの定式化

2026-01-03

Physical Chemistry

Theoretical Chemistry

Differential Geometry

Intrinsic Reaction Coordinate

Meta-IRC

Kenichi Fukui

Riemannian Geometry

last_modified: 2026-01-03

免責事項: 本記事は、Akitomo Tachibana, Kenichi Fukuiによる論文 Theoretica Chimica Acta 1978, 49, 321-347 を基に、生成AIによって作成された解説記事です。数式や論理の正確な理解のためには、必ず原著論文を参照してください。

1. 序論：反応経路理論の深化と座標不変性への挑戦#

1970年、福井謙一によって提唱された固有反応座標 (Intrinsic Reaction Coordinate; IRC) は、化学反応の経路を「ポテンシャルエネルギー曲面 (PES) 上の最急降下経路」として定義することで、遷移状態理論に明確な幾何学的基礎を与えた。初期のIRC理論は、主に直交座標系（質量重み付きカルテシアン座標）に基づいて定式化されていたが、実際の化学反応、特に巨大分子や複雑な有機反応系を扱う場合、結合距離や結合角といった内部座標（一般曲線座標）を用いる方が直観的かつ効率的である場合が多い。

しかし、曲線座標を用いた場合、従来の定義による「勾配」や「力の定数（ヘシアン）」の記述は、座標変換によってその物理的意味が変わってしまうという問題（共変性の欠如）を孕んでいた。StantonとMcIver (1975) らは、曲線座標におけるポテンシャル面の記述が座標系に依存する危険性を指摘していた。

1978年、立花明知と福井謙一は、論文 “Differential Geometry of Chemically Reacting Systems” において、化学反応系を $(n+1)$ 次元のリーマン空間 に埋め込まれた多様体として捉え直すことで、これらの問題を根本的に解決した。本稿では、彼らが構築した「化学反応系の微分幾何学」の枠組みと、そこから導かれる Meta-IRC、拡張ヘシアン、そして反応パターンのトポロジー解析について詳説する。

2. 化学反応系の幾何学的定式化#

2.1 $(n+1)$ 次元ユークリッド空間とリーマン空間の構築#

$N$ 個の原子核からなる化学反応系の配置は、通常 $3N$ 次元の質量重み付きカルテシアン座標 $x = (x^1, \dots, x^{3N})$ で記述される。著者らは、ここにエネルギー次元に対応する座標 $x^{3N+1}$ を加えた、 $(3N+1)$ 次元のユークリッド空間 $E_{3N+1}$ を導入した。

この空間内での系の運動は、位置ベクトル $r = (x, x^{3N+1})$ の軌跡として記述される。ここで、 $x^{3N+1}$ は全エネルギー $E$ （保存量）または断熱ポテンシャルエネルギー $U(x)$ に対応する。

さらに、系を記述するために $n (\le 3N-6)$ 個の一般化座標 $q^1, \dots, q^n$ を導入する。これにより、系は $E_{3N+1}$ 内に埋め込まれた $(n+1)$ 次元のリーマン空間 $R_{n+1}$ として定義される。この空間の線素 $ds^2$ は以下のように与えられる。

ds^2 = (dx^1)^2 + \dots + (dx^{3N})^2 + (dx^{3N+1})^2 = \sum_{i,j=1}^{n} a_{ij} dq^i dq^j + (dx^{3N+1})^2 \quad (2.4)

ここで、 $a_{ij}$ は計量テンソル（metric tensor）であり、座標変換のヤコビアンを用いて次のように定義される。

a_{ij} = \sum_{t=1}^{3N} \frac{\partial x^t}{\partial q^i} \frac{\partial x^t}{\partial q^j} \quad (2.5)

2.2 「飛行機と影」のアナロジー#

論文では、反応系のダイナミクスを直観的に理解するために、非常に示唆に富むアナロジーが提示されている。

ダイナミカル平面 (Dynamical Plane): $x^{3N+1} = E = \text{const.}$ で定義される超平面。
ポテンシャル曲面 (Potential Surface): $x^{3N+1} = U(x(q))$ で定義される超曲面。

化学反応の進行は、ダイナミカル平面という「一定の高度」を保って飛ぶ 飛行機 の運動に例えられる。一方、従来の反応座標の議論で扱われるポテンシャル面上の経路は、その飛行機が眼下の山岳地帯（ポテンシャル曲面）に落とす影に相当する。飛行機（系）は、山岳（ポテンシャル）の地形によって乱気流（力）を受け、揺れ動く。この幾何学的視点により、著者は以下の2つのリーマン面を定義した。

$R_n(a)$ : ダイナミカル平面上の幾何学。計量テンソルは $a_{ij}$ 。運動エネルギーの記述に関連する。
$R_n(b)$ : ポテンシャル曲面上の幾何学。計量テンソルは $b_{ij}$ 。

b_{ij} = a_{ij} + \frac{\partial U}{\partial q^i} \frac{\partial U}{\partial q^j} \quad (2.9)

平衡点（ $\nabla U = 0$ ）においては $b_{ij} = a_{ij}$ となり、2つの空間は接触する。この定式化の最大の利点は、エネルギー座標 $x^{3N+1}$ を明示的に幾何学に取り込むことで、座標変換に対して不変な議論が可能になる点にある。

3. ポテンシャル曲面の微分幾何学#

3.1 単位法線ベクトルと曲率#

ポテンシャル曲面の局所的な形状を議論するために、 $(n+1)$ 次元空間における単位法線ベクトル $n$ が導入される。

n = \frac{N}{|N|}, \quad N = \left( -\frac{\partial x^{3N+1}}{\partial x^1}, \dots, -\frac{\partial x^{3N+1}}{\partial x^{3N}}, 1 \right) \quad (3.15)

特に $R_n(a)$ （等エネルギー面）においては、法線ベクトルは単にエネルギー軸方向を向く（ $n = (0, \dots, 0, 1)$ ）。曲面上の曲線の曲がり具合を表す曲率ベクトル $d^2r/ds^2$ は、法線成分（法曲率 $\rho$ ）と接平面成分（測地線曲率 $\kappa_G$ ）に分解される。

\frac{d^2r}{ds^2} = \rho n + \kappa_G t \quad (3.19)

測地線 (Geodesic curves): $\kappa_G = 0$ となる曲線。曲面上の2点間を結ぶ最短距離（固有距離）を与える。
曲率線 (Lines of curvature): 法曲率 $\rho$ が極値をとる方向（主方向）に沿った曲線。

3.2 デュパンの標形 (Dupin Indicatrix)#

曲面上の点 $P_0$ 近傍のトポロジーを視覚化するツールとして、デュパンの標形 が導入された。これは接平面をわずかに平行移動させた際の、曲面との交線の形状である。

楕円点 (Elliptic point): 第二基本形式が定符号の場合。安定平衡点（極小点）や極大点はこれに該当する。標形は楕円となる。
双曲点 (Hyperbolic point): 第二基本形式が不定符号の場合。遷移状態（鞍点）はこれに該当する。標形は双曲線となり、その漸近線（Asymptotic cones）が重要な意味を持つ。

この幾何学的分類は、後述する反応経路の分岐や合流を理解する上で決定的な役割を果たす。

4. Meta-IRCと拡張ヘシアン#

本論文の核心部分は、従来のIRCを一般の非平衡点からの運動へと拡張し、かつ座標不変な形に定式化した Meta-IRC の導出、および 拡張ヘシアン の定義にある。

4.1 Meta-IRCの一般解#

従来のIRCは平衡点（遷移状態）を出発点として定義されていたが、著者はこれを任意の非平衡点から出発可能な経路へと一般化した。これを Meta-IRC と呼ぶ。 Meta-IRCは、「一般化運動量 $p$ の方向が常にポテンシャルの勾配場 $\nabla U$ の方向に接するような軌跡」 として定義される。これを共変成分（covariant components）を用いて記述すると、以下の微分方程式が得られる。

\frac{dq_1}{\partial U / \partial q^1} = \frac{dq_2}{\partial U / \partial q^2} = \dots = \frac{dq_n}{\partial U / \partial q^n} \quad (4.4)

ここで $dq_i$ は変位ベクトルの共変成分であり、 $dq_i = a_{ij} dq^j$ である。この方程式は、座標変換に対して不変である（テンソル方程式としての性質を持つ）。著者はこの解を具体的に構成し、Meta-IRCが 等ポテンシャル面（Equi-potential surface）に対して直交する曲線 （正確には、等ポテンシャル面上の計量が特定条件を満たす場合）として解釈できることを示した。これは、反応の進行を担う「促進モード (Promoting mode)」を幾何学的に抽出することに他ならない。

4.2 拡張ヘシアン (Extended Hessian) と基準振動#

通常の量子化学計算で用いられるヘシアン行列（力の定数行列） $H_{ij}^{usual} = \partial^2 U / \partial q^i \partial q^j$ は、曲線座標系においてはテンソルとして振る舞わない。すなわち、座標変換によってその固有値（振動数）が変化してしまうという致命的な欠陥がある（平衡点においては $\nabla U=0$ となるため、この問題は顕在化しないが、非平衡点では無視できない）。

TachibanaとFukuiは、共変微分（covariant derivative）を用いることで、任意の点において座標不変な 拡張ヘシアン (Extended Hessian) $H_{ij}$ を定義した。

H_{ij} = U_{ij}^{(2)} = \frac{\partial^2 U}{\partial q^i \partial q^j} - \Gamma_{ij}^k \frac{\partial U}{\partial q^k} \quad (4.30)

ここで、 $\Gamma_{ij}^k$ はクリストッフェル記号（第二種）であり、座標系の曲がり具合を表す項である。第二項 $-\Gamma_{ij}^k (\partial U / \partial q^k)$ が、座標の湾曲に由来する見かけの力を補正する役割を果たす。

この拡張ヘシアンを用いて定義される固有値方程式：

(H_{ij} - \mu a_{ij}) v^j = 0 \quad (4.31)

の解（固有値 $\mu$ および固有ベクトル $v$ ）は、座標変換に対して不変である。これにより、反応経路上の任意の点（非平衡点）において、物理的に意味のある「基準振動」を定義することが初めて可能になった。これは、反応進行中の振動モードの変化や、エネルギー移動（IVR）を議論する上で極めて重要な基盤となる。

4.3 ポテンシャルの不変なTaylor展開#

拡張ヘシアンを用いることで、ポテンシャル $U$ のTaylor展開もまた、任意の曲線座標系において不変な形式で記述できる。

U(x + \Delta x) = U(x) + U_i^{(1)} \Delta u^i + \frac{1}{2} H_{ij} \Delta u^i \Delta u^j + O((\Delta x)^3)

この展開式は、反応経路近傍のダイナミクスを摂動論的に扱う際の基礎式となる。特に、Meta-IRCに沿った局所座標系を用いることで、反応進行方向（促進モード）とそれに直交する振動モード間の結合（カップリング）を厳密に評価することが可能となる。

5. 遷移点におけるトポロジーとパターン認識#

論文の後半では、モデルポテンシャルを用いた具体的な解析が行われ、特に遷移状態（Transition Point, $P_{tr}$ ）近傍でのMeta-IRCの振る舞いが幾何学的に解明された。

5.1 漸近曲線とMeta-IRC#

遷移点 $P_{tr}$ は鞍点（双曲点）であり、そのDupin標形は双曲線となる。この双曲線の漸近線（Asymptotic lines）は、ポテンシャル曲面 $U$ と、その点の高さにある等エネルギー平面 $x^{3N+1} = U(P_{tr})$ との交差曲線に対応する。

著者は、「Meta-IRCは、遷移点においてこれら2本の漸近曲線のなす角を二等分する」 ことを示した（Fig. 8(c)）。幾何学的には、これはMeta-IRCがDupin標形の主軸（Principal axis）に沿って走ることを意味する。反応系がエネルギーの峠を越える際、その経路は「峡谷の底」を這うように進むが、遷移点そのものでは「尾根」を横切る形になる。このとき、Meta-IRCは尾根の最も低い点を、尾根の背骨に対して直交する方向に通過するのである。

5.2 反応パターンの分類#

このトポロジー的視点は、化学反応のパターン認識に応用できる。

異性化反応 (Isomerization): 2つの極小点（反応物と生成物）が、1つの鞍点を介して結ばれるパターン（Fig. 8(a)）。ポテンシャル面の等高線は「8の字」のような交差（separatrix）を持つ。
解離反応 (Fragmentation): 一方の極小点が無限遠へ伸びる谷へと変化するパターン。

著者は、これらの基本パターンに対する「折り畳み (Folding)」や「ねじれ (Twisting)」といった幾何学的操作によって、より複雑な多段階反応のメカニズムを分類・理解できる可能性を示唆した。

5.3 モデル系での実証#

論文では、以下の2次元モデルポテンシャルを用いた解析が行われている。

U(x_1, x_2) = \frac{1}{a^2} x_1^2 (x_1 - 1)^2 + \frac{1}{b^2} x_2^2 \quad (5.1)

この系において、Meta-IRCの厳密解が導出され、遷移点（ $x_1=1/2, x_2=0$ ）において経路が分岐する様子や、漸近線との関係が数値的にも実証された。特に、安定平衡点近傍では無数のMeta-IRCが一点に収束する（confluent）のに対し、遷移点では特定の経路のみが物理的に意味を持つことが示されている。

6. 結論と波及効果#

TachibanaとFukuiによる1978年の研究は、それまで直観的あるいは近似的に扱われていた反応経路の概念を、リーマン幾何学という堅牢な数学的基盤の上に再構築したマイルストーンである。

本研究の主要な成果は以下の通りである。

座標不変性の確立: Meta-IRCと拡張ヘシアンの導入により、どのような曲線座標を用いても物理的に同一の反応経路と振動解析が得られることが保証された。
非平衡ダイナミクスへの拡張: 平衡点のみならず、任意の核配置からの反応進行を記述できるMeta-IRCは、光化学反応やホットアトム化学のような非平衡過程の記述に道を開いた。
幾何学的直観の深化: 「飛行機と影」のアナロジーやDupin標形を用いた解析は、複雑な多次元ポテンシャル曲面のトポロジーを理解するための強力な概念的ツールを提供した。

この理論は、後の “Reaction Ergodography”（反応エルゴード法）の発展の基礎となり、ヘシアンを用いた信頼性の高い探索アルゴリズムの設計思想にも通じる、極めて先駆的な業績であると言える。

参考文献#

Tachibana, A.; Fukui, K. “Differential Geometry of Chemically Reacting Systems”, Theor. Chim. Acta 1978, 49, 321-347.
Fukui, K. J. Phys. Chem. 1970, 74, 4161.
Stanton, R. E.; McIver, J. W. J. Am. Chem. Soc. 1975, 97, 3632.
Glasstone, S.; Laidler, K.; Eyring, H. The theory of rate processes; McGraw-Hill: New York, 1941.
Laugwitz, D. Differential and Riemannian geometry; Academic Press: New York, 1965.

補足#

4.2 拡張ヘシアンと非平衡点における基準振動の定義#

1. 背景：なぜ通常のヘシアンでは駄目なのか？#

ポテンシャルエネルギー曲面 $U(q)$ 上の点における曲率（力の定数）は、通常、ヘシアン行列（二階偏微分行列）によって評価される。

H_{ij}^{\text{usual}} = \frac{\partial^2 U}{\partial q^i \partial q^j}

問題点：テンソル性の欠如#

直交座標系（デカルト座標）では問題ないが、結合距離や結合角といった**曲線座標系（一般化座標）**を用いる場合、この単純な二階微分は座標変換に対してテンソルとして振る舞わない。

座標変換 $q \to \bar{q}$ を行った際、二階微分の変換則は以下のようになる（連鎖律）：

\frac{\partial^2 U}{\partial \bar{q}^\alpha \partial \bar{q}^\beta} = \sum_{i,j} \frac{\partial q^i}{\partial \bar{q}^\alpha} \frac{\partial q^j}{\partial \bar{q}^\beta} \frac{\partial^2 U}{\partial q^i \partial q^j} + \sum_{k} \frac{\partial^2 q^k}{\partial \bar{q}^\alpha \partial \bar{q}^\beta} \frac{\partial U}{\partial q^k}

第一項: テンソル本来の変換（物理的な成分の変換）。
第二項: 座標変換の非線形性（ $\partial^2 q / \partial \bar{q}^2$ ）に由来する項。

平衡点と非平衡点の違い#

平衡点 ( $\nabla U = 0$ ): 勾配 $\partial U / \partial q^k$ がゼロになるため、第二項が消滅する。したがって、通常のヘシアンでも問題なく振動数が定義できる。
非平衡点 ( $\nabla U \neq 0$ ): 反応経路の途中などでは勾配がゼロでないため、第二項（余分な項）が残り、計算する座標系の選び方によって振動数（固有値）が変わってしまうという物理的に受け入れがたい結果を招く。