MOSP: A User-interface Package for Simulating Metal Nanoparticle Structure and Reactivity under Operando Conditions

MOSPの技術要約：オペランド条件下における金属ナノ粒子構造と反応性のシミュレーションのためのユーザーインターフェース・パッケージ

問題提起
金属ナノ粒子（NP）は不均一系触媒において中心的な役割を果たしており、その触媒反応性は原子構造および活性サイトの集合状態と密接に関連している。平均場マイクロキネティックモデリング（MKM）やキネティックモンテカルロ法（KMC）といった従来の速度論的シミュレーションは、固定されたモデル構造に基づいた触媒挙動の予測には成功してきたが、近年のin situ実験による観察では、金属ナノ粒子の構造が変化する反応条件に応じて動的に再構成されることが明らかになっている。静的な構造に依存する従来のアプローチは、このような反応条件に起因する形状変化を考慮できておらず、合理的な触媒設計に不可欠な「動的な構造と特性の関係」の理解における空白を生じさせている。

手法
この制限に対処するため、著者らは、2つの異なる計算モデルを統合したオープンソースのPythonベースのグラフィカルユーザーインターフェース（GUI）であるMOSP（Multi-scale Operando Simulation Package）を提示する。

1. マルチスケール構造再構成（MSR）モデル:

   • Wulff理論に基づき、特定の反応雰囲気下における金属ナノ粒子の平衡形状を予測する。
   • 各結晶面の界面張力（$\gamma_{hkl}^{int}$）を、ガス被覆率、吸着エネルギー、および吸着種間の横方向相互作用を考慮した項を用いて、真空表面張力を修正することで算出する。
   • 吸着種被覆率は、分圧、温度、およびエントロピー変化を組み込んだFowler-Guggenheim吸着等温式を解くことで決定される。
   • このモデルは、有限サイズ効果を考慮しつつ、3 nmより大きい粒子の平衡構造を定量的に予測する。

2. キネティックモンテカルロ（KMC）モデル:

   • MSRモジュールによって生成された特定のナノ粒子構造（またはインポートされた既存の構造）上での反応ダイナミクスをシミュレートする。
   • システムの状態は、元素、位置、および占有値によって特徴付けられる離散サイトの表面構成として定義される。
   • シミュレーションには、時間発展に関するマスター方程式を解くために、Gillespieアルゴリズムの直接法（Direct Method）を用いる。
   • イベントには、吸着、脱離、拡散、および表面反応が含まれる。速度定数は、衝突理論（吸着・脱離用）および遷移状態理論（反応・拡散用）から導出され、BEP関係および一般化配位数（GCN）を利用して、サイト固有のエネルギープロファイルをマッピングする。

主な貢献

• 統合されたワークフロー: MOSPは、熱力学的な構造予測（MSR）と速度論的な反応シミュレーション（KMC）を結合する統一されたフレームワークを提供し、動的な構造・活性関係の探索を可能にする。
• ユーザーのアクセシビリティ: パッケージはwxpythonで構築された直感的なGUIを備えており、ユーザーはコードを書くことなく、複雑なパラメータ（横方向相互作用行列、スケーリング関係、反応ネットワークなど）を入力できる。また、pyOpenGLを用いたナノ粒子構造のインタラクティブな可視化パネルや、pandasおよびmatplotlibを用いたシミュレーション結果の表示機能も備えている。
• モジュール性: 結合運用を目的として設計されているが、MSRとKMCの各モジュールは独立して機能させることができる。MSRは平衡構造をエクスポートでき、KMCは既存のXYZ形式の構造をロードできる。
• オープンソース: ソースコードは公開されており、コミュニティによる採用と拡張を促進している。

結果とデモンストレーション
論文では、以下の2つの例証を通じてMOSPの妥当性を検証している。

1. Ptナノ粒子上でのCO酸化:

   • MSRモジュールは、変動する圧力（0–2000 Pa）および温度（800 K）の下でのPtナノ粒子（半径2.5 nm）の再構成に成功した。
   • 結果は、圧力依存的な構造進化を示した。圧力が上昇するにつれて、（111）面の割合が減少し、（110）面の割合が増加し、ナノ粒子の形状が切頂八面体から準菱形十二面体へと変化した。
   • その後の500 PaにおけるKMCシミュレーション（1000万ステップ）により、エッジおよびコーナーが活性の主要な寄与因子であることが明らかになった。計算された被覆率はCOが0.01、Oが0.11であり、ターンオーバー頻度（TOF）は$8.33 \times 10^5$ site$^{-1}$s$^{-1}$であった。
   • パフォーマンスのベンチマークでは、MSRの実行時間は原子数（$N_{atoms}$）に対して$N_{atoms}^{2.4}$でスケールし、KMCはほぼ線形（$N_{atoms}^{0.93}$）でスケールすることが示された。

2. Cuナノ粒子上での水ガスシフト（WGS）反応:

   • 650 K、4000 Paにおける2190原子の切頂八面体Cuナノ粒子に対してKMCシミュレーションを実施した。
   • モデルには、カルボキシル機構に従う5種類の吸着種と14種類のイベントが含まれている。
   • 結果は、COが支配的な吸着種であることを示した（被覆率 0.67）。
   • 決定的なことに、シミュレーションは原子レベルでの速度論的結合を浮き彫りにした。低配位サイトはCO$_2$生成に有利であり、高配位サイトはH$_2$生成に有利であることが特定されたが、両生成物の全体的なTOFは同程度であった。

意義
本論文は、MOSPを、オペランド条件下における金属ナノ粒子の構造と触媒特性を予測するための効率的なツールとして位置付けている。静的な構造の仮定と動的な実験の実態との間の溝を埋めることで、MOSPは最適な反応条件のスクリーニングと、触媒の合理的設計を支援する。著者らは、本パッケージが標準的なノートPCで動作するように設計されており、初期スクリーニングにおいてハイパフォーマンス・コンピューティング・クラスターに依存することなく、高度なオペランド・シミュレーションへのアクセスを可能にしていることを強調している。