High-throughput Search for Transition Metal Chalcogenides

High-throughput Computational + Experimental Search for Transition Metal Chalcogenides

Phys. Rev. B 94, 045105 (2016) — Narayan, Bhutani, Rubeck, Eckstein, Shoemaker, Wagner | DOI: 10.1103/PhysRevB.94.045105

⚠️ 無機材料科学論文（凝縮系物理）。高スループット計算スクリーニング + 実験検証の方法論として参照。創薬への直接応用は限定的。

① 背景と目的

遷移金属カルコゲナイド（Se/S 系）は Iron-based 超伝導体・遷移金属ダイカルコゲナイドなど機能的電子・磁気特性を示すが、酸化物（ICSD 2300+）に比べてセレン化物・硫化物（ICSD ~350）は著しく未探索。

ターナリー相空間 27 系が ICSD に登録なし → "空"の相図が実際に存在するのか不明

DFT 単独の予測は誤差が大きく（既知ターナリーの 40% のみ凸包上）実験検証が不可欠

→ 確率的イオン置換モデル + DFT + 実験合成の3段構えで 27 相図を包括的に探索

② 計算手法

ICSD データベース → イオン置換確率モデル（MLE）
↓ 上位100候補を選択
↓ DFT 構造緩和（QE + PBE）
↓ PYMATGEN で凸包 ΔH 計算
↓ 24 候補に絞り込み

27 → 24

相図 27 系から 24 候補材料に絞り込み

② 確率モデルの要点

イオン対 (α,β) の置換重み λ_α,β を ICSD から MLE で学習。新規組成の確率スコアを計算し、候補をランク付け。

③ 実験結果

全 24 候補: ターナリー相の合成に失敗

結論: これら 27 相図はバルク合成では「空」。ただし薄膜・高圧合成では可能性残存。

④ 探索対象系

④ DFT 計算パラメータ

④ 限界点

DFT (PBE) の精度限界 — 既知ターナリーの 60% が hull 上にない

スピン配置: 強磁性初期化のみ（反強磁性未考慮）

バルク合成のみ — 薄膜/MBE/水熱合成は未試行

無機材料向け → 有機分子への直接移植不可

④ 方法論的教訓

否定的結果を誠実に報告することの価値

24 候補全て失敗という否定的結果を系統的に記録・公開。計算予測と実験の乖離を定量化した点は、仮想スクリーニングの偽陽性問題と同質の教訓。

0/24

合成成功率（全候補が実験で失敗）

⑤ ケムインフォマティクスへの応用可能性

🔄 バイオイソスター置換スコアラー
イオン置換確率 → 有機官能基置換確率モデル。ChEMBL から学習した置換重みを lib/molgen の MolgenYaml スコアラーに統合。

📊 多段フィルタリング戦略
確率モデル → DFT の二段式は docking score → MM-GBSA → FEP の階層ワークフローに対応。計算コスト最適化の手法として参考に。

⚖️ False Negative 最小化
Bayes 的な閾値設定で「見逃し」を減らす考え方は、バーチャルスクリーニングの recall 最大化（hit を見逃さない設計）に直接応用可能。

⚠️ 計算予測の限界認識
DFT でも大量の false positive が生じた事実は、DL ドッキング手法の予測精度過信への警鐘。実験検証の重要性をリマインド。

応用先モジュール

適用先	提案（類推レベル）
lib/molgen	DB駆動バイオイソスター置換確率スコアラー
lib/fep	多段フィルタリング（確率→粗→精密評価）
lib/docking	false negative 最小化の閾値設定ロジック

※ 無機材料論文のため直接移植不可。設計思想・方法論の参考として活用

本研究のインパクト