In Silico

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

スクリーニングの漏斗——天文学的な候補を数個まで削る

2026/6/21 (更新: 2026/8/2) シリーズ「コンピュータで材料を探す」 第2回 / 全4回

材料探しには、ひとつの壁がある。 ありうる材料は天文学的に多く、全部は作って試せない。でも、宝はどこかにある。

ではどうするか。やることは、宝探しというより「ふるい分け」だ。 膨大な候補を、何段もの“ふるい”にかけて、最後に実際に作る価値のある数個まで削っていく。コツは、ふるいをかける順番にある。

※ 概念図(フロー) 天文学的 → 数千 → 数百 → 数十 → 数個 ① 安いルール・制約でざっくり ② 機械学習で高速に性質を予測 ③ 精密な物理計算(DFT)で確認 ④ 実際に合成して実験 上:安い 速い 下:高い 遅い・正確
上ほど安く速くざっくり、下ほど高く遅く正確。安いふるいで大半を落とし、高いふるいは生き残りだけにかける。

なぜ「安いふるいを先に」なのか

漏斗の各段は、1個あたりのコストが違う。

ポイントは順番だ。安いふるいで大半を落としてから、高いふるいを少数にかける。 もし順番が逆で、いきなり全部を実験したら——何百年あっても終わらない。この作法が実際にどれだけ絞り込むかは、材料探索の実例が示している。ある大規模計算(GNoME)では、機械学習で絞り込んだ末に220万個の構造が「それまでの研究を基準にすれば安定」=当時の凸包より下と判定された。そのうち38万1000個が、更新後の凸包の上に残る新規材料として報告されている。ただしこれはDFT計算による判定だ。同じ論文は、この安定構造のうち736件が独立に実験で実現済みだった(=著者らが合成したのではなく、既存の実験報告と照合して一致した)とも報告している1。38万1000に対して0.2%——漏斗の下段ほど数が落ちるという本稿の構図が、そのまま数字に出ている。同じ「安い絞り込み→高い確認」の順番は、新薬探し(バーチャルスクリーニング)でも採られている。ある大規模スクリーニングの報告では、クラスタリングと簡易ドッキングで候補を1/10に減らしてから個別の精密ドッキングにかけ、最後に機械学習で真の当たりを選り分けている2

ふるいは「正確さ」と「安さ」を取引している

ここに、この分野のいちばん正直な難しさがある。

安くて速いふるい(②の機械学習)は、ざっくりだから間違える。 間違いには2種類ある。

だから機械学習は、最終結論ではなく“絞り込み役”として使うのが筋だ。安いふるいで候補を減らし、最後は精密計算と実験が確かめる。しかも二つの誤りは効き方が違う——空振り(偽陽性)は高い③④が後で捨ててくれるが、早い安いふるいで見逃した(偽陰性)宝は、二度と下流に現れない。だから早い段では“取りこぼさなさ”(再現率)を優先したくなる。

取りこぼし(偽陰性)だけを恐れていればいいわけでもない。Matbench Discovery のリーダーボード(2026年7月時点)では、機械学習モデルの一つ CHGNet(v0.3.0)を単独のふるいとして使った場合の適合率は0.51と報告されている3。適合率とは「安定」と判定した候補のうち実際にDFTでも安定だった割合のことだ。ただしこの0.51は、モデルが「安定」と言った候補を全部そのまま下流へ流す場合の値である。同じ論文は、実務ではそうせず予測の安定度で並べて予算の分だけ上から検証する、と書いている——偽陽性は候補列の「安定度が低い側の端」に溜まるので、そこで打ち切れば全体の発見率を損なわずに避けられる。実際、上位1万件だけを検証する想定で測り直すと、同じCHGNetの適合率は0.85に上がる3。∴ 効くのは「どのモデルか」だけではない。どこで打ち切るかが、同じモデルのまま歩留まりを変える。値がモデルによって大きく違うのも確かで、同じリーダーボードの上位モデルは全件評価でも0.9台に達する3。再現率と適合率のどちらを優先すべきかは、段ごとの1個あたりのコスト比と下流で処理できる件数で決まり、一律には決まらない。

もう一つ、正直に言っておくべきことがある。実験の一つ手前にある③精密計算(DFT)も、絶対的な正解ではない。DFTが弾き出す生成エネルギーは、実験値と比べると系統的なずれを持つ。あるプレプリントは、機械学習で補正した後でも生成エネルギーの平均絶対誤差は50 meV/原子を下回る水準——実験の不確かさと同程度——だと報告している4。補正前の標準的な計算(GGA)のずれはそれより大きい。凸包からの距離が数十meV/原子しか離れていない候補では、このずれが安定・不安定の判定を左右しうる。つまり「精密な計算で確認した」は「実験と完全に一致する」を意味しない——漏斗の各段は、下に行くほど正確になるが、計算でいちばん精密な段ですら近似だという留保がつく。


多段スクリーニング(安い近似で絞り、高精度計算・実験で確認する漏斗)は、計算材料科学や創薬のバーチャルスクリーニングで広く使われる考え方。DFT=密度汎関数理論、機械学習原子間ポテンシャルなどの位置づけは、本シリーズおよび姉妹シリーズ『凸包』で扱う一般的な知見に基づく。なお姉妹シリーズ『凸包』は GNoME の 220万→38万 という漏斗の実例を扱う。本稿はその一般形——なぜ安いふるいを先にかけるのか——を見ている。

出典

  1. A. Merchant et al., “Scaling deep learning for materials discovery,” Nature 624, 80–85 (2023)。先行研究(Materials Project)を基準にすれば安定な構造を220万個見つけ、うち38万1000件が更新後の凸包の上に残る新規材料だと報告した——220万と38万1000は基準の異なる別々の数である。本稿が一般形として引く実例(詳細は姉妹シリーズ『凸包』)。 https://www.nature.com/articles/s41586-023-06735-9

  2. “A Machine Learning-Enabled Pipeline for Large-Scale Virtual Drug Screening” (2021)。クラスタリングと限定的ドッキングで化合物ライブラリを1/10に絞り、残りを個別ドッキングにかけ、最後にニューラルネットで真陽性/偽陽性を分類する多段パイプラインを報告——材料探索と同型の「安い→高い」順序を創薬で採った一例(慣行の調査ではなく、あるチームのパイプラインの報告)。 https://pmc.ncbi.nlm.nih.gov/articles/PMC8478848/

  3. J. Riebesell ら, “Matbench Discovery — A framework to evaluate machine learning crystal stability predictions,” arXiv:2308.14920。評価はWBMテストセットに対して行われ、正解は Materials Project の参照構造から作ったDFT計算の凸包(凸包より上、すなわち0 eV/原子超は安定でないと見なす)。同フレームワークが継続的に更新・公開しているリーダーボードでは、2026年7月時点でCHGNet(v0.3.0)の適合率(precision)は0.514、上位モデルは0.928。適合率は「安定」と予測した候補のうちDFTでも安定だった割合で、単独運用時の誤検出リスクがモデルによって大きく異なることを示す(リーダーボードは更新されるため値は変わる)。 https://matbench-discovery.materialsproject.org/ 2 3

  4. “Correcting DFT formation energies towards experimental accuracy using foundational MLIPs and latent-feature delta-learning,” arXiv:2607.18092 (2026)。MC3D(Materials Cloud三次元結晶データベース)の生成エネルギーを、OQMD・Materials Project および実験の生成エンタルピーと比較した報告。r²SCAN級の基盤MLIP(PET-OMATPES)による補正でGGA比の平均絶対誤差を40%以上削減し、さらに潜在特徴のdelta-learningで50 meV/原子未満(実験の不確かさと同程度)まで下げたとしている。漏斗の「精密な」確認段階も、実験との比較では近似にとどまることを示す。 https://arxiv.org/abs/2607.18092

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