last_modified: 2026-01-04

免責事項: 本記事は、Mina Jafari and Paul M. Zimmermanによる論文 Phys. Chem. Chem. Phys. 2018, 20, 8860-8871 (doi: 10.1039/c8cp00044a)を基に、生成AIによって作成された解説記事です。数式や論理の正確な内容は、必ず原著論文を参照してください。

1. 序論：表面化学反応における探索空間の複雑性とその課題#

不均一触媒反応や薄膜成長プロセス（CVD, ALDなど）において、固体表面上での化学反応メカニズムを原子レベルで解明することは、材料設計およびプロセス最適化の中核をなす課題である 。量子化学計算、特に密度汎関数理論（DFT）の発展により、反応経路（Reaction Path; RP）や遷移状態（Transition State; TS）の第一原理的な解析が可能となった。しかし、表面反応系は、気相反応と比較して圧倒的に高い自由度と複雑性を有している。

1.1 既存手法の限界と歴史的背景#

従来の反応経路探索は、研究者の「化学的直感（Chemical Intuition）」に大きく依存していた 。すなわち、反応物と生成物をあらかじめ想定し、その間を結ぶ経路をNudged Elastic Band (NEB) 法などで探索するか、あるいは想定される反応座標に沿って構造最適化を行う手法が主流であった。

しかし、未知の反応や複雑なネットワークを持つ系においては、以下の課題が顕在化する。

1. 高次元ポテンシャルエネルギー局面（PES）の呪縛: 可能な吸着サイト、配向、分子内自由度の組み合わせが爆発的に増加するため、網羅的な探索が困難である 。
2. バイアスと予断: 研究者が想定していない反応経路（例えば、直感に反する協奏的プロセスや複雑な転位）は見過ごされるリスクが高い 。
3. 表面特有の困難: 表面には多様な吸着サイト（Top, Bridge, Hollow等）が存在し、吸着質の配向や表面原子との相互作用が反応障壁に劇的な影響を与える 。

これらの課題を克服するため、Metadynamics やArtificial Force Induced Reaction (AFIR) 法 など、PESの低次元部分空間を探索する自動探索手法が開発されてきた。しかし、これらの多くは反応座標の定義を必要とするか、あるいは表面反応特有の吸着サイトサンプリング機能が不足している場合があった 。

1.2 本稿の目的：S-ZStruct法の位置づけ#

JafariとZimmerman（2018）は、気相反応向けに開発されたグラフ理論ベースの探索手法「ZStruct」を拡張し、表面反応に特化した**S-ZStruct（Surface ZStruct）**を提案した。本手法は、Single End Growing String法（Single-Ended Growing String Method; SE-GSM）と組み合わせることで、初期状態のみを入力とし、中間体や生成物を事前知識なしに自動探索することを可能にする 。

本稿では、S-ZStructのアルゴリズム的構造、グラフ理論による反応性の定義、および実際の表面反応（プロパン酸分解およびTiN ALD）への適用事例を通じて、その学術的意義と実用性を詳細に論じる。

2. 理論的枠組みとアルゴリズムの構造#

S-ZStructは、化学反応を「原子間の結合（グラフの辺）の組み換え」として捉えるグラフ理論的アプローチを採用している 。PES全体を無作為に探索するのではなく、化学的に妥当な結合変化の組み合わせを組織的に生成し、それをSE-GSMの駆動座標（Driving Coordinates）として利用することで、効率的なTS探索を実現する。

2.1 グラフ理論による反応性の定義#

分子構造は、原子をノード、化学結合をエッジとするグラフ $G(V, E)$ として表現される。ここで、2つの原子 $i, j$ 間の結合性（Connectivity）は、原子間距離 $d_{ij}$ と各原子の共有結合半径 $r_i, r_j$ を用いて以下のように判定される 。

ここで $C$ は判定のための閾値係数である。S-ZStructでは、ユーザーが指定した「反応中心原子（Reactive Atoms）」に対して、この結合性が変化する操作（結合の切断および生成）を適用することで、反応候補を列挙する 。

また、非物理的な構造の生成を防ぐため、各原子種に対して配位数（Coordination Number）の制約条件（上限・下限）が課される 。例えば、水素原子の配位数は最大1、炭素原子は4といった化学的妥当性がグラフ操作時にチェックされる。

2.2 S-ZStructの5段階アルゴリズム#

S-ZStructは以下の5つのステージから構成され、これを反復することで反応ネットワークを拡張していく 。

Stage 1: 吸着サイトの同定 (Binding Site Finding)#

表面反応において、吸着質がどのサイトに結合するかは反応性を決定づける重要な因子である。S-ZStructは、Atomic Simulation Environment (ASE) を用いて、入力された表面スラブ構造（FCC, BCC, HCP等の結晶面）上の標準的な吸着サイト（Atop, Bridge, FCC/HCP Hollow）を自動的に同定する 。

Stage 2: 構造の整列とサンプリング (Aligning and Sampling)#

発見された吸着サイトに対し、吸着質を配置する。表面反応の特異性として、吸着質の回転配向（Orientation）が活性化エネルギーに大きく影響する 。S-ZStructは、吸着質の結合軸周りに特定の角度（例：$0^\circ, 30^\circ, 60^\circ$）で回転させた複数の初期構造を生成し、さらに指定された半径内の近傍サイトに対してもサンプリングを行うことで、初期構造の多様性を確保する 。

Stage 3: 幾何構造最適化 (Geometry Optimization)#

生成された多数の初期構造候補に対し、DFT計算による構造最適化を実行する。これは計算コストが最も高いステップの一つであるが、局所安定点からの探索を開始するために不可欠である 。

Stage 4: ユニークな構造の選別 (Unique Structure Identification)#

最適化の結果、同一の極小点に収束した構造が多数生成される可能性がある。計算資源の浪費を防ぐため、構造間の二乗平均平方根偏差（RMSD）を計算し、閾値（例えば $0.01 \text{ \AA}$）以下の構造を重複として削除する 。残ったユニークな構造はエネルギー順にソートされ、低エネルギー構造から順に次のステップへ進む 。

Stage 5: 異性体生成とSE-GSMの実行 (Isomer Generation & SE-GSM)#

最も重要なステップである。グラフ理論に基づき、反応中心原子における結合の切断・生成の組み合わせ（駆動座標）を生成する 。

• 単分子反応（解離など）: 結合している原子対の結合を切断し、解離フラグメントを近傍の空き吸着サイト（反応中心から $3.0 \text{ \AA}$ 以内）へと誘導するベクトルを生成する 。
• 二分子反応: 2つの分子間での結合生成と、それに伴う結合切断の組み合わせを生成する。一度に最大4つの結合変化（生成2、切断2）までを許容することで、協奏的な反応ステップも探索対象とする 。

生成された駆動座標を用いて Single-Ended GSM (SE-GSM) が起動される。SE-GSMは、反応物構造と接線ベクトル（駆動座標）のみから反応経路を成長させ、TSと生成物を同時に探索する手法である 。これにより、最終的な生成物が未知であっても探索が可能となる。

3. ケーススタディ1：Pd(111)上のプロパン酸分解反応#

S-ZStructの単分子反応探索能力と吸着配向サンプリングの有効性を検証するため、バイオマス変換において重要なプロパン酸（Propanoic Acid）のPd(111)表面での分解反応が解析された 。

3.1 実験的・理論的背景#

この系はHeydenらによって広範に研究されており、C-H結合やO-H結合の切断を経てエタンなどの脱酸素生成物に至る複雑なネットワークが知られている 。従来研究では化学的直感に基づいて経路が構築されていたが、S-ZStructが自動的に既知の経路を再現できるか、さらには新規経路を発見できるかが検証の焦点となった。

3.2 探索結果と考察#

S-ZStructは、中間体として重要な**プロピオニル（$\text{CH}_3\text{CH}_2\text{CO}$）およびカルボニルエチリデン（$\text{CH}_3\text{CH}\text{CO}$）**からの分解経路を探索した。

配向依存性の発見#

特筆すべき成果として、初期吸着構造の配向（Rotation Angle）が反応バリアに与える劇的な影響が明らかになった。例えば、カルボニルエチリデンの分解において、吸着角度を $0^\circ, 30^\circ, 60^\circ$ と変化させた際、それぞれの配向から異なるTSと生成物サイトが得られた 。

• エネルギー差: $60^\circ$ の配向から出発する経路は、$0^\circ$ の場合と比較して一様に低い活性化エネルギーを示した 。これは、吸着質と表面原子の立体的配置が遷移状態の安定化に寄与していることを示唆する。
• 生成物の吸着サイト: 分解により生じる $\text{CO}$ フラグメントは高配位サイトを好む一方、$\text{CH}_3$ フラグメントは低配位のブリッジサイト等に安定化される傾向が自動的に見出された 。

新規経路の同定#

S-ZStructは、Heydenらの先行研究で報告されていなかった**メチル基脱離経路（Methyl Fragmentation）**を自動的に検出した 。 従来研究では「重要ではない」と仮定され除外されていたこの経路（P1, P6経路）が、実際には計算可能なエネルギー障壁を持つことが示された 。これは、化学的直感によるバイアスを排除し、網羅的な探索を行う本手法の優位性を如実に示している。

4. ケーススタディ2：Si(100)上のTiN原子層堆積（ALD）#

次に、より複雑な二分子反応および多段階プロセスへの適用例として、TDMAT（テトラキスジメチルアミドチタン）とアンモニアを用いたTiNのALDプロセスが解析された 。

4.1 ALDプロセスの未解明問題#

TiN ALDにおいて、前駆体であるTi(IV)錯体がどのようにして還元され、Ti(III)を含む導電性膜を形成するかは、長年の未解明問題であった 。提案されているメカニズムには、ヒドラジン誘導体の脱離や、$\beta$-水素脱離などが存在するが、詳細な素過程やその競合関係は明らかではなかった 。

4.2 反応サイクルの自動構築#

S-ZStructを用いた多段階の探索により、以下の反応サイクルが構築された。

1. 核発生: 表面OH基とアンモニアのプロトン交換による $\text{NH}_x$ 末端の形成 。
2. Ti堆積: トランスアミノ化（Transamination）によるTDMATの吸着とジメチルアミンの脱離 。
3. 架橋形成: 隣接する $\text{NH}_2$ 基とのプロトン移動・配位子交換を経て、Ti-Nネットワークが形成される 。

ここまでは既知の知見と整合するが、S-ZStructはさらにTiの還元ステップに踏み込んだ探索を行った。

4.3 Ti(IV)からTi(III)への還元メカニズム#

S-ZStructは、Ti還元に関する複数の反応候補を生成し、その中から**隣接する配位子間での$\beta$-水素移動（$\beta$-hydride transfer）**が最も低い活性化エネルギーを持つことを突き止めた 。

立体障害と反応性#

特に興味深いのは、配位数と立体構造の関係である。

• 反応16: 両方のTi中心が四面体構造を持つ場合、H受容体となるNと移動するHの距離は $3.5 \text{ \AA}$ と遠く、反応障壁は高い 。
• 反応17（最優経路）: 一方のTiが歪んだ三角両錐構造（5配位）をとる場合、立体的な歪みによりH-N距離が $3.1 \text{ \AA}$ まで短縮し、反応障壁が大幅に低下する 。
• 角度解析: 反応物からTSに至る過程で、$\text{C}-\text{H}\cdots\text{N}$ の角度変化は、反応16では $104.3^\circ \to 170.2^\circ$ と劇的であるのに対し、反応17では $116.7^\circ \to 161.4^\circ$ と変化が小さく、初期構造がすでに反応に適した配置（Predisposed）にあることが判明した 。

この結果は、TiN膜中でのTi還元が、特定の中間体構造における立体障害によって促進されるという新たなメカニズムを提唱するものであり、S-ZStructが単なる経路探索を超えて、反応の微視的起源に迫る能力を持つことを示している 。

5. 議論：計算効率と「化学的直感」の役割#

5.1 計算コストと効率#

S-ZStructの特筆すべき点はその効率性である。分子動力学（MD）法が1つの反応経路を見つけるために約10,000回の勾配計算を必要とするのに対し、SE-GSMと統合されたS-ZStructは、1つのTSあたり1,000回以下の勾配計算で済むとされる 。さらに、グラフ理論によるプレスクリーニング（配位数制限など）により、DFT計算を行うべき構造を事前に絞り込むことで、無駄な計算資源の消費を抑えている。

5.2 化学的直感との協働#

本手法は「化学的直感」を排除することを目的としているが、それは研究者の洞察を否定するものではない。むしろ、S-ZStructは人間が思いつきにくい経路（例：プロパン酸のメチル脱離やTi還元の具体的な立体配置）を提示することで、研究者の直感を補完・拡張するツールとして機能する 。研究者は、自動生成された膨大なネットワークの中から、速度論的に支配的な経路を特定し、その物理的意味を解釈することに注力できる。

6. 結論#

JafariとZimmermanによるS-ZStructは、表面化学反応の探索におけるパラダイムシフトをもたらす手法である。

1. 体系的な探索: グラフ理論と幾何学的サンプリングを組み合わせることで、吸着サイト、配向、結合変化の多様な組み合わせを網羅的に探索可能とした 。
2. 実用的な適用性: 単分子分解反応（プロパン酸）から複雑な成膜プロセス（TiN ALD）に至るまで、多様な表面反応系に対して有効に機能し、実験事実を説明しうる詳細なメカニズムを提示した 。
3. 発見的価値: 既知の経路の再現にとどまらず、還元ステップにおける$\beta$-水素移動の立体的な起源など、新規かつ深い洞察を提供することに成功した 。

S-ZStructは、SE-GSMの強力なTS探索能力を表面反応へと拡張し、計算化学による材料設計や触媒開発の加速に大きく貢献する基盤技術であると結論付けられる。

参考文献#

1. Jafari, M.; Zimmerman, P. M. Phys. Chem. Chem. Phys. 2018, 20, 8860-8871.
2. Zimmerman, P. M. J. Comput. Chem. 2013, 34, 1385-1392.
3. Zimmerman, P. M. J. Chem. Phys. 2013, 138, 184102.
4. Maeda, S.; Morokuma, K. J. Chem. Theory Comput. 2011, 7, 2335-2345.
5. Henkelman, G.; Uberuaga, B. P.; Jónsson, H. J. Chem. Phys. 2000, 113, 9901.
6. Behn, A.; Zimmerman, P. M.; Bell, A. T.; Head-Gordon, M. J. Chem. Phys. 2011, 135, 224108.
7. Lu, J.; Behtash, S.; Heyden, A. J. Phys. Chem. C 2012, 116, 14328-14341.
8. Musschoot, J.; Xie, Q.; Deduytsche, D.; Van den Berghe, S.; Van Meirhaeghe, R.; Detavernier, C. Microelectron. Eng. 2009, 86, 72-77.
9. Elam, J.; Schuisky, M.; Ferguson, J.; George, S. Thin Solid Films 2003, 436, 145-156.
10. Bouman, M.; Zaera, F. J. Electrochem. Soc. 2011, 158, D524.