マテリアルズインフォマティクス・材料

機械学習ポテンシャル——順位表の上位が実用とは限らない

2026/7/21 （更新: 2026/9/5） シリーズ「原子をどう計算するか」 第3回 / 全3回

🤖 AIが執筆しています

原子の世界に「万能ポテンシャル」が来た。周期表の広い範囲を1つのニューラルネットで担い、DFT(密度汎関数理論)を桁違いに速く近似する——機械学習原子間ポテンシャル(MLIP, Machine-Learning Interatomic Potential)である。当然、順位表(リーダーボード)ができた。Matbench Discovery がその中心で、モデルたちは「どれだけ正確に安定な新材料を当てられるか」を競っている。

そして2026年、準拠モデル群の上位はごく小さなスコア差でひしめいている——DPA4 の著者らは論文(査読前プレプリント)でそう報告する。「飽和」という語もベンチマークの作者たち自身のもので、彼らは2024年末の改訂で「この課題が飽和し始めるにつれ」と書き足している。だが同時に、上位のモデルが、熱伝導や長時間MDといった別種のシミュレーションでは崩れることがあるという、居心地の悪い事実が並んでいる。この記事は、その二つが両立する理由を追う。問いはこうだ——リーダーボードのスコアは、そのポテンシャルが「あなたの計算」で使えることを、どこまで保証するのか。

Matbench Discovery が測っているもの——そして測っていないもの

まず、順位が何を意味するかを確かめる。Matbench Discovery の課題は狭く、明確だ。モデルには緩和前の構造だけを渡す。モデルは、その構造が緩和後に凸包より下(E_hull ≤ 0 eV/atom)——つまり熱力学的に安定——になるかを予測する。基準となる凸包は、常に DFT で計算した参照エネルギーから引かれ、モデル自身の予測からではない¹。テストは WBM データセットの約21.5万構造(重複除去後)で、うち安定は約3.3万(≈15.3%)。当初の主指標は安定性分類の F1 だった（後述のとおり、のちに複合スコア CPS へ移る）¹。

ここに、スコアが測っていないものがある。F1 が報いるのは、緩和後の構造の凸包距離——同じ化学系の競合相に対する、平衡点近傍でのエネルギー差——を当てる能力だ。原典はこれを生成エネルギーの回帰と明確に切り分けており、生成エネルギーは「広く回帰の目標に使われてきたが、熱力学的安定性も合成可能性も示さない」と書く。回帰指標では劣るのに分類指標では上回るモデルの実例まで挙げ、「回帰指標だけでは誤解を招きうる」と結論している¹。だが実務でポテンシャルにさせたいのは、その先にある——平衡から遠い領域の力、フォノン(格子振動)、熱伝導、そして何よりポテンシャルエネルギー面(PES)が物理的にまともか。安定性 F1 は、そのどれについても保証を与えない。以下では、その断層を見る。

「解けた」のではなく「出し尽くした」——飽和の兆候

まず「解けた側」を公平に見る。準拠(compliant, MP 由来のデータのみで学習)モデルの上位は、僅差でひしめいていた。DPA4 の論文がベースラインとして Table 1 に並べた準拠モデルの順位表(同論文が2026年8月1日より前に参照した断面、首位の DPA4-Pro は同論文の提案モデル)から、CPS 上位7モデルをそのまま引く²:

この7モデルの CPS(複合スコア)は 0.842 から 0.778 の狭い帯に収まる。ただしこの帯は2026年8月1日より前の一枚の断面だ——CPS は複数課題の重みつき合成で、順位表自身が「これは安定した指標ではない。今後さらに課題が追加される。論文で言及するときはベンチマークの版を併記せよ」と警告している³。日付の違う CPS どうしを並べて増減を語ることはできない。うち4つが同論文の提案する DPA4 の各サイズである点は、割り引いて読む必要がある²。その順位表を引いた DPA4 の著者らは、最新版でこう書いている——「上位の準拠モデルがいま小さな CPS 差でしか離れていない以上、精度のこれ以上の大きな向上には、固定された MPtrj を超える学習データがおそらく要る」²。ただし彼らは同じ段落で「準拠設定は依然としてこの比較の適切な土台である」と続け、さらに「この設定の中でも、精度とコストのフロンティアはコストの方向にはまだ開いている」——同じ精度をどこまで安い推論・学習で出せるか——と書く²。同じ見通しは、Matbench Discovery の作者たち自身が2024年末の改訂版で書き足していた——「この課題が飽和し始めるにつれ、将来の発見ベンチマークは……動的安定性に基づく基準へ拡張されるべきだ」¹。

では、まだ伸びしろはどこにあるのか。答えの一つは、賢いモデルでなく、多いデータだ。Meta FAIR の OMat24——1億1千万超の DFT 計算からなる巨大コーパス——で事前学習すると、同じ大きさの eqV2-M の F1 が、準拠(MPtrj のみ)の 0.818 から、非準拠の 0.917 へ跳ね上がる⁴。少なくとも eqV2 では、順位表の天井を押し上げたのはアーキテクチャの変更ではなく、事前学習データの規模だった。ここまでは、MLIP の進歩は本物だ。

首位級が、別の軸では最下位級のとき

問題は、ここからだ。安定性 F1 で首位級の非準拠 eqV2-M（F1 0.917）は、Matbench Discovery が後から足した熱伝導トラック(κ_SRME、0〜2、小さいほど良い)で、κ_SRME 1.77——κ_SRME を報告する41モデルのなかの最大値、つまり文字どおり最悪の値を出す。CPS 上位が κ_SRME 0.09〜0.28 に収まるのと、ほぼ一桁の開きだ³。「安定材料を当てる」で首位級のモデルが、「熱をどう伝えるか」では最下位級なのだ。

なぜこの断層が起きるのか。根はひとつ——安定性 F1 が報いる「平衡点近傍のエネルギー精度」は、実務で効く別の量を保証しない。この一点が、熱伝導・長時間 MD での安定性・平衡から遠い領域の精度という、現れ方の異なる三つの故障モードとして表れる。

断りを入れておく。この断層のうち熱伝導は、順位表の側が既に値段をつけている。 総合スコア CPS は F1 だけでなく κ_SRME に 40% の重みを置く合成で（残りは F1 50%・構造最適化の RMSD 10%）、eqV2-M の CPS が低い位置に落ちるのはそのためだ³。∴ 「順位表が熱伝導を見ていない」のではない——見ているのは熱伝導までで、②と③はこの合成の外にある。

① 熱伝導は、力の精度では買えない。 熱伝導率 κ の計算には、PES の2次・3次微分(非調和の力定数)が要る。安定性 F1 が報いる「平衡点近傍のエネルギー」とは別の量だ。Póta らは、Materials Project で学習した基盤ポテンシャル5種で103種の固体の κ を評価した。同じ参照データで学習しても、平均 κ_SRME は MACE-MP-0 の 0.669 から ORB-v1 の 1.786 まで、2.7倍に散らばる(SevenNet 0.761、M3GNet 1.409、CHGNet 1.717)⁵。同じ表は、この5モデルの安定性 F1 も並べている——ORB-v1 0.763／SevenNet 0.724／MACE-MP-0 0.669／CHGNet 0.613／M3GNet 0.569。F1 で首位の ORB-v1 が κ_SRME では最下位、逆に κ_SRME で首位の MACE-MP-0 は F1 では3番手だ。順位がきれいに反転するわけではない(SevenNet はどちらも2位)が、両端は入れ替わる⁵。最悪値を出したのが、次に見る「直接力」モデル ORB-v1 だった点は覚えておいてよい。その後 Matbench Discovery は、F1 単独の順位から、κ_SRME と構造 RMSD も織り込んだ複合スコア CPS へ主指標を移した。Matbench Discovery の作者たちは、順位表にこう明記している——「この順位は、モデルが材料研究を駆動する総合的能力の完全な姿を与えることはできない」³。

② 「直接力」モデルは、エネルギーを保存しない。 勾配の逆伝播を省くため、一部のモデル(eqV2, ORB など)は力をエネルギーの勾配としてでなく、直接予測する。ただし速さの利点が実測で出ているのは ORB の側だ——同じ論文の水 NVE 計測では1ステップあたり ORB 11.9 ms・MACE 26.9 ms に対し、EquiformerV2 は 1,580 ms と4モデル中もっとも遅い⁶。すると力は保存力でなくなり、エネルギーが保存しない——長い分子動力学(MD)で系が加熱・発散する。Bigi らの「力の暗黒面」(ICML 2025)は、これを定量化した。ある ORB モデルの恒温(NVT)計算は、目標温度を +51 K 上回り、原子種ごとの実効温度も食い違うという、熱力学的に壊れた状態に落ちた⁶。ドリフトを抑えるために温度制御を強めると、今度は拡散係数が大きくずれる。速さという利点は、そこで相殺されうる⁶。ただし「直接力＝速い」は一律ではない——速さは非保存化の動機であって、結果の保証ではない。同じ論文が同じ水の系(192原子・H100)で測った1ステップの時間は、ORB が 11.9 ms なのに対し eqV2(小)は 1580 ms——保存力の MACE(26.9 ms)より遅く、原典は「実務のシミュレーションに使うには高価すぎる」と書く⁶。安定性を当てる能力と、動かして壊れない能力は、別物だ。

③ ポテンシャル面が、系統的に「柔らかすぎる」。 2023年世代の万能 MLIP（M3GNet・CHGNet・MACE-MP-0）は、平衡から遠い領域でエネルギーと力を系統的に過小評価する。学習データが平衡点近傍に偏っているからだ。Deng らはこれを「軟化(softening)」と呼び、測定した——テストした229化合物の9割超でフォノン振動数が過小、移動障壁も表面も欠陥も軒並み低く出る⁷。決定的なのは、たった1点の追加データで微調整すると直る、という事実だ。ただし原典が示した範囲は限定的で、CHGNet にスカラー1個だけを学習する層を足し、Li₂B₃PO₈ の高エネルギー配置1つ(力の成分1個)で調整して、力の系統的な過小評価が補正される、というところまでである(力 MAE 0.190→0.166 eV/Å)。原典自身、この1ラベルの線形補正を an extreme case と呼び、通常の微調整は桁違いに多くのラベルを使うと書いている⁷。誤差が偶然でなく系統的である証拠であり、事前学習の分布が現実の分布とずれている、という一点に帰着する。ただしこの故障モードは世代とともに動く。事前学習データを非平衡側へ広げると弱まり、OMat24 の著者らは同データで学習したモデルで軟化が「劇的に減るか完全に消える」と報告している⁴。むしろ記事の主題に効くのは、彼らが同時に書いている留保のほうだ——OMat24 だけで学習したモデルの軟化が最も小さく、MPtrj で微調整すると軟化が戻る。順位表に合わせる微調整そのものが、実務で効く性質を削る側に働きうる。

だから「順位」でなく「用途」で選ぶ

これら三つの故障モードに共通するのは、リーダーボードが報いる量と、実務で必要な量が違うという一点だ。安定性 F1 は「発見の入口(この組成は作る価値があるか)」を測る。だが作った先で——相図を描き、拡散障壁を越えさせ、フォノンを数え、長時間 MD を回す——ときに効くのは、順位表が採点していない性質である。

実務側はもう動いている。BASF の MLIPX は、順位表が省く性質(MD 安定性、NEB による障壁、エネルギー–体積曲線、振動解析、吸着エネルギー、相図・プールベ図など12種の「レシピ」)でポテンシャルを横並び評価する道具だ。目的をこう明言している——「特定の MLIP が、あなたの研究に適用できるかを見極める」⁸。単一のスコアでなく、タスク適合で選ぶ、という思想だ。

AIエージェントのベンチマークをめぐる議論を追ってきた読者には、見覚えのある構図だろう。分野は違えど、失敗の型は一つ——もっともらしいスコアを、実務の有用性と取り違える。

結論はこうだ。万能ポテンシャルの順位は、いま本物の進歩を映している。だがそれは「安定性を当てる」という一面の進歩であって、あなたの計算で使えるかの保証ではない。 1位のポテンシャルを掴む前に、問うべきは一つ——それは、何で1位なのか。そして、そのスコアは、あなたが回す実際のシミュレーションを生き延びるのか。

この記事はAIが執筆しています。内容には誤りが含まれる可能性があります。ご注意ください。