１．はじめに
マテリアルズインフォマティクス（MI）は、材料科学における新たなアプローチとして注目されており、第一原理計算や分子動力学（MD）計算を活用しながら、機械学習を用いた予測モデルの構築が進められてきました。しかし、MIの実用化にはいくつかの大きな課題が存在します。第一原理計算（DFT）は非常に高精度ですが、計算コストが極めて高く、大規模な材料探索には適していません。また、分子動力学（MD）計算では適切な力場が必要であり、新しい材料系には適用が難しいことが多いです。その解決策として、機械学習を用いたニューラルネットワークポテンシャル（NNP）が注目されています。本記事では、M3GNET（Materials Graph Neural Network）がどの程度の精度で実用的に使えるのかを評価することを目的とします。特に、第一原理計算ソフトウェアであるQuantum ESPRESSO（QE）と比較しながら、M3GNETの利点と限界を明確にし、どのような場面で活用できるのかを探ります。
２．比較条件および概要
本研究では、白金（Pt）スラブ表面における水（H₂O）の吸着挙動を評価するために、吸着エネルギーおよび原子間距離の変化を解析し、QEとM3GNETの計算精度と計算時間を比較しました。
使用したモデル
Ptスラブ
H₂O分子（吸着種）
Pt-H₂O吸着モデル（Pt表面上の水分子の分離シミュレーション）
Figure 1: Structure of Pt slab, H₂O molecule, and Pt-H₂O adsorption model
使用した計算環境
- Quantum ESPRESSO（QE）: 32 vCPUs @ Intel Xeon Platinum 8473C
- M3GNET: NVIDIA A100 1GPU
シミュレーションに用いたソルバー
- Quantum ESPRESSO（QE）: バージョン7.3を使用し、k点はガンマ点、カットオフエネルギーは適切に設定。使用した擬ポテンシャルは標準的なPBE-GGAに基づくものを採用。
- M3GNET: NanoLaboが提供するLAMMPS連携のM3GNETを利用し
ポテンシャルにはM3GNet-MP-2021.2.8-PESを使用した。
３．比較結果
本研究では、Ptスラブ上における水の吸着エネルギーを評価することで、M3GNETとQEの精度を比較しました。具体的には、スラブPt、単体H₂O、そしてPtスラブ上にH₂Oを配置したPt-H₂Oモデルのエネルギー値を計算し、それらの差分を求めました。以下の表は、QEとM3GNETにおけるエネルギー値（eV）を比較したものです。ΔEは式（１）から算出しています。
ΔE = E（Pt-H₂O） - E（スラブPt） - E（H₂O)（1）
これらの結果から、M3GNETは定性的な傾向を再現できるものの、エネルギー値の精度には若干の差があることが分かりました。ΔEには小さな差異が見られたものの、その符号は一致していました。Ptスラブ上にH₂Oを配置したPt-OH-Hモデルの最適化結果を比較したところ、OHとHの分離後に付着するPtの位置が異なるため、視点によって構造の見え方が異なりますが、最適化された構造自体は類似しており、M3GNETの最適化構造はQEとよく一致していました（Figure 2）。
Figure 2: Optimized Structures of Pt-H₂O Model Using QE and M3GNET
計算時間では、顕著な差が見られました。以下の図は、QEとM3GNETにおける計算時間（Wall Time）の比較を示しています。
Figure 3: Speedup Comparison of QE and M3GNET in Wall Time
このグラフからも明らかなように、M3GNETはQEと比較して大幅に短い時間でエネルギー計算を実行できることが確認されました。
- Ptモデルでは約 356倍 の高速化
- H₂Oモデルでは約 23倍 の高速化
- Pt-OH-Hモデルでは約 610倍 の高速化
このように、M3GNETは第一原理計算と比較して数百倍の計算速度向上を実現しており、大規模な材料探索や分子動力学（MD）シミュレーションにおいて極めて有用であると考えられます。
４．まとめ
本記事では、M3GNETの精度と計算時間の観点から、DFTと比較しました。その結果、M3GNETは大規模な材料探索において非常に有用である。
ΔEには若干の差が見られたものの、M3GNETの結果はQEと定性的に一致しており、最適化構造も類似していることが確認されました。特に、計算時間に関してはM3GNETが圧倒的に高速であり、Ptモデルで約356倍、H₂Oモデルで約23倍、Pt-OH-Hモデルでは約610倍の高速であったことから、電子構造やバンド計算などの詳細な物性解析には依然としてDFTが必要ですが、大規模な分子動力学シミュレーションやスクリーニングにはM3GNETが大いに役立つことが確認されました。M3GNETを用いることで、これまで時間的・計算資源的に不可能だった大規模シミュレーションが可能となる。
さらに、M3GNETはMIにおいても有用であり、材料候補のスクリーニングや特性予測を迅速に行うための強力なツールとなり得ます。従来の第一原理計算では時間がかかる大規模な材料探索を効率的に進めることで、新しい材料の発見や設計を加速する可能性を持っています。
次回の記事では、M3GNETを用いたスクリーニングや材料予測の具体的な活用方法について紹介します。
キーワード：
#M3GNET #GNN #QuantumESPRESSO #LAMMPS #MaterialsInformatics #MI