計算化学ノート / 04 GRRM

Reaction Route Mapping

GRRM

局所最適化だけでは見つけにくい反応経路、遷移状態、中間体を、ADDF・AFIR familyによって系統的に探索します。

位置づけ
反応経路・停留点探索
主な入力
初期構造、電荷・多重度、計算engine
主な出力
EQ、TS、DC、path、reaction network

1. 概要

Global Reaction Route Mapping（GRRM）はpotential energy surface上の平衡構造、遷移状態、解離経路を自動・半自動で探索する戦略とプログラムです。通常のOpt=TSが遷移状態の良い初期guessを必要とするのに対し、GRRMは局所minimumから未知経路を広く探すことを目的とします。

versionを明記する

GRRM17、20、23、26では機能、keyword、対応interface、並列化が異なります。本ページは概念を共通化し、実際の入力は使用版の公式manualに合わせます。

2. Gaussian・xTBとの連携

GRRMは経路探索を制御し、各geometryで必要なenergy、gradient、Hessianを電子状態計算engineへ要求します。Gaussianは代表的な連携先です。

GRRM minimum optimization concept
# MIN/B3LYP/6-31G(d)

0 1
O    0.000000    0.000000    0.000000
H    0.758602    0.000000    0.504284
H   -0.758602    0.000000    0.504284

xTBなどGaussian以外を使う場合は、専用interfaceまたはgeneral external interfaceを介して要求されたenergy・gradient・HessianをGRRM形式へ返します。利用可能性と設定はGRRM版、xTB版、計算環境に依存するため、wrapper名や非公式scriptを一般化せず、実際のversionとinterfaceを記録します。

3. Potential energy surfaceとGRRM出力

記号意味検証
EQequilibrium structure、局所minimum虚振動なし
TStransition state、一次鞍点目的modeの虚振動が一つ
DCdissociation channelfragment、電荷・spin、energy asymptote
PT/path経路上の点またはpath情報接続先とenergy profile

構造同一性の判定はenergy差と原子間距離差のthresholdに依存します。柔軟分子、対称分子、ほぼ縮退した配座では、別構造の統合や同一構造の重複が起こり得ます。

4. ADDF

Anharmonic Downward Distortion Following（ADDF）は、局所minimumの周囲で調和近似から下方へ外れる方向、すなわちanharmonic downward distortionを検出し、その方向を追跡してTSやDCへ至る経路候補を探します。

ADDF path自体はIRCと同一ではありません。path最高点付近をTS guessとして最適化し、最終的にはfrequencyとIRCで停留点と接続を確認します。小分子の網羅的探索には強力ですが、自由度増加に伴い探索方向とforce evaluation数が急増します。

5. AFIRの基本概念

Artificial Force Induced Reaction（AFIR）は、選択したフラグメント間に人工的な引力または斥力を加えた修正potentialを最適化し、反応経路候補を生じさせます。

FAFIR(Q)=E(Q)+ραR(Q)

Eは元のpotential energy、Rはフラグメント間距離に基づく人工力項、αは人工衝突energyに関係するparameter、ρは力の向きです。人工力で得たpathをそのまま物理的反応pathとせず、元のPES上で精密化します。

6. SC-AFIR

Single-Component AFIRは、単一分子または一つのaggregate内部でfragment pairを選び、分子内の結合生成・切断・転位経路を探索します。予期しない反応経路を広く調べる用途に向きます。

  • fragment selection反応中心を狭く指定すると効率は上がるが、未知経路を除外し得ます。
  • collision energy parameter大きいほど高い障壁候補へ到達し得る一方、非現実的pathも増えます。
  • path priority既知経路や重要fragment pairを先に計算できます。
  • kinetic treatment利用版によってRCMC等との連携機能が異なります。

7. MC-AFIR

Multi-Component AFIRは複数のreactant component間に人工力を加え、会合、結合形成、解離を伴う経路候補を探索します。componentの相対配置と向きに依存するため、初期配置、探索回数、fragment定義を記録します。

二分子反応では、接近方向、立体障害、非共有結合性pre-complexが選択性へ影響します。MC-AFIRで候補を生成し、重要なpathを高い理論レベルで再最適化する階層的運用が現実的です。

8. 経路の精密化と検証

  1. 低コスト探索
    xTBまたは小さいbasisのDFTでpath候補を生成します。
  2. 構造整理
    重複EQ/TS、fragment、spin stateを確認します。
  3. RePATH・TS最適化
    候補pathのmaximumから元PES上のTSを求めます。
  4. Frequency・IRC
    停留点の次数と両側の接続を確認します。
  5. 高精度energy
    重要構造を共通条件で再評価し、熱化学・溶媒・配座を加えます。
「見つからなかった」は「存在しない」ではない

探索範囲は初期構造、fragment定義、人工力parameter、energy cutoff、計算level、停止条件に依存します。negative resultには探索条件を併記します。

9. 参考文献・公式資料

  1. AFIR-web Support Site, GRRM manuals.
  2. S. Maeda, K. Ohno, K. Morokuma, “Systematic Exploration of the Mechanism of Chemical Reactions,” Phys. Chem. Chem. Phys. 2013. DOI
  3. S. Maeda et al., “Artificial Force Induced Reaction (AFIR) Method,” Chem. Rec. 2016. DOI
  4. GRRM manual, ADDF introduction.

次章: PySCF