[論文レビュー] AutoMeKin2021: An open-source program for automated reaction discovery

Emilio Martínez‐Núñez, George L. Barnes+9 more|arXiv (Cornell University)|May 7, 2021

Mass Spectrometry Techniques and Applications被引用数 1

ひとこと要約

AutoMeKin2021 は、希少事象の分子動力学、グラフ理論に基づく遷移状態検出、化学的知識ルール、および結合次数時間系列解析を統合した、反応機構の自動発見を目的としたオープンソースのソフトウェアパッケージである。高スルーレットで自動化された複雑な反応ネットワークの探索を可能にし、強化サンプリングと統計解析を組み合わせることで、多様な化学系における予期しないまたは非共有結合経路の検出が顕著に向上する。

ABSTRACT

AutoMeKin2021 is an updated version of tsscds2018, a program for the automated discovery of reaction mechanisms (J. Comput. Chem. 2018, 39, 1922-1930). This release features a number of new capabilities: rare-event molecular dynamics simulations to enhance reaction discovery, extension of the original search algorithm to study van der Waals complexes, use of chemical knowledge, a new search algorithm based on bond-order time series analysis, statistics of the chemical reaction networks, a web application to submit jobs, and other features. The source code, manual, installation instructions and the website link are available at: https://rxnkin.usc.es/index.php/AutoMeKin

研究の動機と目的

• 従来の手動的手法をはるかに超える複雑な反応機構を自動的かつスケーラブルに発見するフレームワークの開発。
• 従来の遷移状態探索の限界を克服し、分子動力学軌道からの反応的イベント検出を統合すること。
• ファンデルワールス錯体や非共有結合相互作用を含む自動機構発見への拡張。
• 化学的知識の統合、結合次数時間系列解析、統計的ネットワーク特徴化を通じた効率性および正確性の向上。
• ジョブ送信と再現可能な計算化学ワークフローを可能にするユーザーフレンドリーでWebアクセス可能なプラットフォームの提供。

提案手法

• 反応的イベントを誘発するため、適応的振動エネルギー配置を用いた短時間の反応的分子動力学（MD）シミュレーションを用いる。
• 原子間距離を共有結合半径に基づいて表現する隣接行列を用いて、分子構造をグラフ理論で表現する。
• 結合の切断/形成を検出する基準（式 2）を採用し、非近接原子が最近接に移行した場合に反応的イベントを特定する。
• MD軌道からの遷移状態検出を目的とした、結合次数時間系列解析に基づく新規探索アルゴリズムを実装する。
• BXDE法を用いて希少事象の加速を実現し、まれな反応的イベントのサンプリングを向上させる。
• NetworkXを用いて化学反応ネットワークを構築・解析し、経路およびエネルギー障壁の統計的特徴付けを行う。

実験結果

リサーチクエスチョン

• RQ1 MDスナップショットからの遷移状態の自動検出は、複雑で直感的でない反応機構を信頼性高く明らかにできるか？
• RQ2 複数の同時進行する結合の形成・切断を伴う系において、結合次数時間系列解析は遷移状態をどの程度効果的に特定できるか？
• RQ3 化学的知識の統合が、発見された機構の正確性および関連性をどの程度向上させるか？
• RQ4 BXDEのような希少事象サンプリング技術は、高エネルギー障壁または低確率経路の発見を顕著に向上させられるか？
• RQ5 反応ネットワークの統計的性質（例：接続性、分岐性）は、機構の複雑さおよび選択性とどの程度相関するか？

主な発見

• AutoMeKin2021 は、プロトン化ウラシルのような複雑な系においても遷移状態を効果的に同定でき、高レベルの経路から電荷を帯びた断片化経路が明らかになった。
• 結合次数時間系列解析により、複数の結合形成/切断を伴う反応においても遷移状態の検出が可能となり、検出のロバスト性が向上した。
• BXDEを用いた希少事象手法の統合により、高エネルギー障壁反応の観測に必要なシミュレーション時間が著しく短縮され、まれなイベントの効率的サンプリングが可能になった。
• NetworkXを用いた反応ネットワークの統計的解析により、分岐点や動力学的ボトルネックといった構造的特徴が明らかになり、機構解釈の支援が得られた。
• WebアプリケーションおよびSingularityを用いたコンテナ化デプロイメントにより、多様な計算環境で再現可能でスケーラブルかつアクセス可能なジョブ送信が実現した。
• Entos Qcore84 のサポートにより、スレート的および高水準の計算が効率的に行えるようになり、パリレン + NO2 やプロトン化ウラシルといった系への応用が成功裏に実施された。