皆さんは、分子が踊っているところを見たことがありますか？
有機反応を扱う研究者として、
X線結晶構造解析で分子構造を「静止画」として眺めるだけでは、
どうしても物足りなさを感じてしまうときがあります。
反応とは、分子が動き、形を変え、
エネルギー障壁を越えて別の構造へと変換される“過程”そのものです。
その過程の中で、反応の行方を決定づけているのが、
一瞬だけ存在する 遷移状態（Transition State, TS） です。
このTSは寿命が極めて短く、
実験的に直接観測することはほぼ不可能です。
それにもかかわらず、
反応速度、収率、立体選択性、さらには副反応の有無に至るまで、
そのすべてがTS構造によって支配されていると言っても過言ではありません。
では、私たち研究者はどのようにして、
この「見えない構造」を覗き込むことができるのでしょうか。
その答えの一つが、計算化学による遷移状態探索です。
量子化学計算を用いれば、反応座標上を分子がどのように移動し、
どのような構造を経由して生成物へ至るのかを、
原子レベルで追跡することが可能になります。
とはいえ、実務的には、TS探索は決して簡単ではありません。
初期構造の与え方を少し誤るだけで、
計算は容易に別の極小構造へと落ち込み、
本来求めたいTSには一向に収束しません。
「TS初期構造をどう作るか」という問題は、
計算化学を扱う多くの研究者が一度は直面する壁でしょう。
今回は、その壁を越えるための実践的手法として、
xTBに実装されたメタダイナミクス（MTD）法を用い、
遷移状態に近い構造を高速にスクリーニングする方法を紹介します。
本手法の狙いは、
厳密なTSをいきなり最適化することではありません。
まず「TSっぽい構造」を大量に、短時間で集め、
それを出発点としてTS探索を行う――
このアプローチにより、試行錯誤と計算時間を大幅に削減しつつ、
TS探索の成功率を飛躍的に高めることができます。
有機合成研究者の方が、
自らの反応機構を自分の手で可視化するための第一歩として、
本稿が役立てば幸いです。
それでは、長い前振りもここまでで、初回は「単分子」における
xTBを用いたMTDシミュレーションによる遷移状態候補構造の高速探索法について、
具体的な操作手順とパラメータ設計の考え方を交えながら解説していきます。
1. 計算戦略の全体像
本手法の基本方針は以下である：
高温・強バイアス・短時間のMTDで構造遷移を強制的に起こし、
遷移状態付近の構造を大量に回収 → それをTS最適化の初期構造に使う
このアプローチにより：
初期構造の試行錯誤を激減
TS探索成功率の劇的向上
計算時間の大幅短縮
が実現できる。
2. 初期構造の準備
まず、初期構造を xyz形式 で準備する。
Avogadro上で作成できるのであれば、可視化もAvogadroで行うので、そちらがベターです。
3. 計算環境へのログインと作業ディレクトリの作成
作業ディレクトリの作成
mkdir mtd_test
cd mtd_test
ここに：
初期構造：xx.xyz
条件ファイル：metadyn.inp
を配置する。
もし計算サーバを利用する場合はファイル転送は：
sftpFileZilla
等を利用しましょう。
4. metadyn.inp の基本構成と意味
公式マニュアル：
https://xtb-docs.readthedocs.io/en/latest/mtd.html
基本構造は以下：
$metadyn
save=10
kpush=0.1
alp=0.2
$end
この3つのパラメータが MTD挙動のほぼ全てを決定する。
5. MTDパラメータの物理的意味と設計指針
5.1 save：ヒル追加頻度
save = 10
→ 10ステップごと（dt=1 fsなら10 fsごと）にヒルを追加
| save | 意味 |
|---|---|
| 5 | 超高速・荒い |
| 10 | 高速探索用（推奨） |
| 50 | 精密探索 |
| 100以上 | 低速・高精度 |
TSスクリーニング目的では：
save = 10 が最適。
5.2 kpush：バイアス強度（最重要）
kpush = 0.1
→ 構造をどれだけ強く押し出すか
| kpush | 挙動 |
|---|---|
| 0.001 | 非常に弱い |
| 0.01 | 穏やか |
| 0.05–0.1 | 高速スクリーニング |
| 0.2以上 | 地形破壊 |
TSスクリーニング目的では：
kpush = 0.05 – 0.1
が最も効率が良いが、分子サイズなどにもよりけり。
最初は少し走らせて止めて、適宜チューニングする。
5.3 alp：ヒルの幅（探索解像度）
alp = 0.2
| alp | 特性 |
|---|---|
| 0.1 | 超精密 |
| 0.2–0.3 | 実用的 |
| 0.4以上 | 雑・荒い |
高速探索では：
alp = 0.2 – 0.3
が最適。
6. TSスクリーニング向け・推奨初期MTD設定
遷移状態っぽい構造を高速収集するための実用設定：
metadyne.inp$metadyn
save = 10
kpush = 0.05
alp = 0.3
$end
追加計算条件の目安
| 項目 | 推奨値 |
|---|---|
| 温度 | 330–350 K |
| タイムステップ | 1 fs |
| 計算時間 | 20–50 ps |
metadyne.inp$md
time=50
step=1
temp=333
$end
$metadyn
save=10
kpush=0.005
alp=0.3
$end
7. MTDシミュレーションの実行
xtb --md --input metadyn.inp xx.xyz
正常終了すると：
xtb.trj
が生成される。
8. エネルギー解析と構造の可視化
8.1 エネルギー変化の抽出
grep energy xtb.trj > mtd.txt
→ エネルギーのみを抽出
→ Excel等でエネルギープロファイル作成可
これにより：
反転タイミング
エネルギー極大点
を特定できる。
8.2 トラジェクトリの可視化
mv xtb.trj xtb.xyz
として：
Avogadro
で読み込む。
File → Open → xtb.xyz
Extensions → Animation → 再生
fpsを下げるとゆっくり再生でき、
構造遷移の瞬間が直感的に把握できる。
9. xtb.xyz ファイル形式の理解
構造は以下の繰り返しで記述される：
原子数
エネルギー行
座標
例：
186
energy: -265.141180200059 ...
B 0.000000 0.000000 -2.326118
N -0.585299 -2.333618 -0.817113
...
これが1ステップ分の構造を表す。
10. TS候補構造の抽出戦略（最重要）
よくある間違い：
❌ 最大エネルギー構造をそのままTSと見なす
→ これは高確率で失敗する。
正しい構造抽出法：
構造遷移が起きる瞬間（研究者としての直感で）を特定
反応座標（距離・角度・二面角）が：
極値を通過
速度が反転
する瞬間の構造を複数抽出
5構造程度をTS候補として保存
11. 抽出構造からTS最適化への接続
TS候補構造を保存後：
原子間距離
二面角
結合角
の 1〜3個を固定 してGaussianなどで厳密な構造最適化を行い、
そこから通常のTS探索手順に移行する。
この方法により：
TS収束率：劇的に向上
試行回数：1/5〜1/10
となる。
12. パラメータ調整指針（実務者向け）
| 現象 | 対処 |
|---|---|
| 動かない | kpush ↑ |
| 変形しすぎ | kpush ↓ |
| 変化が粗すぎ | alp ↓ |
| 遷移が遅い | save ↓ + kpush ↑ |
まとめ
xTBのMTDは、
「単分子のTS探索を人力から半自動化するための最強ツール」
であると個人的には思っています。
特に：
構造変化が大きい反応（ヘリセンの反転障壁など）
反転・環反転・異性化
配座変換
において、劇的な効率改善をもたらすでしょう。
以上、【遷移状態の探し方】として、xTBを用いたメタダイナミクスによるTS候補構造の高速探索法を紹介しました。
本手法は、厳密なTS最適化そのものではなく、
TSに近い初期構造を効率よく準備することに主眼を置いています。
初期構造の質がTS探索の成否を大きく左右することを考えると、
本手法は、計算化学に不慣れな有機合成研究者にとっても、
非常に強力な武器になるはずです。
その２では、2分子間反応におけるTS探索の高速化法について紹介する予定です。
近年では、北海道大学の前田先生らによって開発された AFIR法 をはじめ、
反応経路探索を自動化・高速化する優れた手法が次々と提案されています。
これらの手法をどのように実務に落とし込むか、
具体的な計算戦略とともに解説できればなと。
それでは、また。