Crystal Stability Prediction Metrics

This task measures how effectively a model can triage hypothetical crystals for DFT validation. Models must identify structures that lie on or below a fixed DFT-computed Materials Project convex hull while minimizing costly false positives.

Reference data: the WBM test set. F1 (higher is better) is the harmonic mean of precision and recall. Switch between the full test set and unique prototypes to compare overall accuracy against performance on the deduplicated, more out-of-distribution subset.

Leaderboard

43 of 53 eligible models · Sorted by F1 (descending) · Discovery: Unique Prototypes · Row numbers follow this filtered order.

Filters:
Submit a model RSS
# ↑
1 EquiformerV3+DeNS-OAM 0.978 0.931 6.074 0.928 0.933 0.018 0.868 0.067 2026-04-07
2 TECE-OAM-RRA-1.0 0.978 0.929 6.073 0.928 0.930 0.018 0.871 0.066 2026-07-05
3 EquFlashV2* 0.978 0.929 6.069 0.928 0.930 0.018 0.873 0.066 2026-06-11
4 Prophet-OAME-MBD 0.978 0.928 6.065 0.927 0.928 0.018 0.869 0.067 2026-09-13
5 GRACE-3L-OAM-L 0.977 0.925 6.041 0.923 0.926 0.018 0.875 0.065 2026-07-02
6 eSEN-30M-OAM 0.977 0.925 6.069 0.928 0.923 0.018 0.866 0.067 2025-03-17
7 PET-OAM-XL 0.977 0.924 6.075 0.929 0.920 0.019 0.864 0.068 2026-01-10
8 MatRIS-10M-OAM 0.976 0.921 6.039 0.923 0.918 0.019 0.871 0.066 2025-10-29
9 EquFlash* 0.975 0.919 5.983 0.915 0.922 0.019 0.871 0.066 2025-06-23
10 eqV2 M 0.975 0.917 6.047 0.924 0.910 0.020 0.848 0.072 2024-10-18
11 TACE-OAM-L 0.972 0.910 5.898 0.902 0.919 0.020 0.868 0.067 2026-04-09
12 Nequip-OAM-XL 0.971 0.906 5.869 0.897 0.915 0.020 0.872 0.066 2025-11-30
13 SevenNet-Omni-i12* 0.971 0.906 5.954 0.910 0.901 0.021 0.868 0.067 2026-01-12
14 ORB v3 0.971 0.905 5.912 0.904 0.907 0.024 0.821 0.078 2025-04-05
15 AlphaNet-v1-OAM* 0.968 0.901 5.747 0.879 0.924 0.024 0.831 0.076 2025-05-12
16 Allegro-OAM-L 0.966 0.895 5.674 0.867 0.923 0.022 0.868 0.067 2025-09-08
17 Nequip-OAM-L 0.967 0.893 5.823 0.890 0.895 0.022 0.865 0.068 2025-09-08
18 DPA-3.1-3M-FT 0.963 0.884 5.667 0.866 0.903 0.023 0.869 0.067 2025-06-05
19 GRACE-2L-OAM-L 0.964 0.883 5.840 0.893 0.874 0.022 0.862 0.068 2025-09-09
20 GRACE-2L-OAM 0.963 0.880 5.774 0.883 0.878 0.023 0.862 0.068 2025-02-06
21 ORB v2 MPA 0.965 0.880 6.041 0.924 0.841 0.028 0.824 0.077 2024-10-11
22 EquiformerV3+DeNS-MP 0.956 0.863 5.479 0.838 0.890 0.029 0.840 0.074 2026-04-07
23 MatterSim v1 5M 0.959 0.862 5.852 0.895 0.831 0.024 0.863 0.068 2024-12-16
24 DPA-4.0.1-Pro-MPtrj 0.956 0.857 5.609 0.857 0.856 0.029 0.836 0.074 2026-06-11
25 MACE-MPA-0 0.954 0.852 5.582 0.853 0.851 0.028 0.842 0.073 2024-12-09
26 MatRIS-10M-MP 0.951 0.847 5.422 0.829 0.865 0.031 0.824 0.077 2025-10-29
27 eSEN-30M-MP 0.946 0.831 5.260 0.804 0.861 0.033 0.822 0.078 2025-03-17
28 NequIP-GNoME 0.948 0.829 5.523 0.844 0.814 0.035 0.785 0.085 2024-02-03
29 GRACE-1L-OAM 0.944 0.824 5.255 0.803 0.846 0.031 0.842 0.073 2025-02-06
30 eqV2 S DeNS 0.939 0.815 5.042 0.771 0.864 0.036 0.788 0.085 2024-10-18
31 Eqnorm MPtrj 0.929 0.786 4.844 0.741 0.838 0.040 0.799 0.083 2025-05-26
32 HIENet 0.929 0.777 4.932 0.754 0.801 0.041 0.793 0.084 2025-07-01
33 ORB v2 MPtrj 0.922 0.765 4.702 0.719 0.817 0.045 0.756 0.091 2024-10-14
34 Nequip-MP-L 0.921 0.761 4.704 0.719 0.809 0.043 0.791 0.084 2025-09-08
35 SevenNet-l3i5* 0.920 0.760 4.629 0.708 0.821 0.044 0.776 0.087 2024-12-10
36 Nequix MP 0.914 0.751 4.455 0.681 0.836 0.044 0.782 0.086 2025-08-17
37 Allegro-MP-L 0.915 0.751 4.516 0.690 0.823 0.044 0.778 0.087 2025-09-08
38 Nequix MP PFT 0.914 0.748 4.479 0.685 0.825 0.044 0.784 0.085 2026-01-08
39 GRACE-2L-MPtrj 0.895 0.691 4.163 0.636 0.757 0.052 0.741 0.094 2024-11-21
40 MACE-MP-0 0.878 0.669 3.777 0.577 0.796 0.057 0.697 0.101 2023-07-14
41 BAM-MP-core 0.869 0.623 3.682 0.563 0.698 0.060 0.712 0.099 2026-07-16
42 CHGNet 0.851 0.613 3.361 0.514 0.758 0.063 0.689 0.103 2023-03-03
43 M3GNet 0.812 0.569 2.882 0.441 0.803 0.075 0.585 0.118 2022-09-20

F1 vs Params

Use the axis/color/size selectors to compare models across any pair of metrics and metadata. The default view plots F1 against model size, colored by energy MAE.

  • Params Params · 43 models
Log Scale

Diagnostics

Explore calibration, element-level errors, and error distributions to see where stability predictions succeed or fail.

Test set: WBM

256,963 crystal structures · 85 elements

Data by Wang, Botti, and Marques

Chemically substituted crystals with DFT relaxations test stability predictions and how closely relaxed geometries reproduce DFT. Discovery uses a unique-prototype subset; geometry optimization evaluates the relaxed structures.

Public initial and relaxed structures, energies, and symmetry references.

Shared by discovery and geometry optimization. These summaries describe the full WBM set.

H 5.23k
He
Li 8.1k
Be 859
B 8.05k
C 6.39k
N 9.16k
O 27.9k
F 10.1k
Ne
Na 6.72k
Mg 12.5k
Al 18.7k
Si 17.2k
P 10.7k
S 10.6k
Cl 6.19k
Ar
K 8.36k
Ca 10.2k
Sc 8.5k
Ti 7.31k
V 5.73k
Cr 7.14k
Mn 14.7k
Fe 17.6k
Co 15.7k
Ni 19.5k
Cu 18.5k
Zn 14.6k
Ga 16.3k
Ge 18.6k
As 10.9k
Se 10.4k
Br 5.92k
Kr
Rb 7.8k
Sr 10.7k
Y 8.85k
Zr 8.13k
Nb 6.17k
Mo 2.09k
Tc 1.1k
Ru 11.7k
Rh 16.1k
Pd 15.5k
Ag 7.57k
Cd 8.57k
In 16.6k
Sn 17.4k
Sb 8.78k
Te 8.06k
I 4.78k
Xe 1
Cs 7.19k
Ba 10.5k
La 8.79k
Ce 7.59k
Pr 8.16k
Nd 8.04k
Pm 1.19k
Sm 7.92k
Eu 5.47k
Gd 4.67k
Tb 8.71k
Dy 8.43k
Ho 8.06k
Er 8.05k
Tm 8.4k
Yb 9.16k
Lu 5.82k
Hf 6.63k
Ta 5.42k
W 1.96k
Re 1.95k
Os 6.18k
Ir 11.3k
Pt 15.5k
Au 11.9k
Hg 5.35k
Tl 7.98k
Pb 9.44k
Bi 6.09k
Po
At
Rn
Fr
Ra
Ac 1.42k
Th 8.07k
Pa 2.05k
U 6.65k
Np 4.93k
Pu 6.7k
Am
Cm
Bk
Cf
Es
Fm
Md
No
Lr
Rf
Db
Sg
Bh
Hs
Mt
Ds
Rg
Cn
Nh
Fl
Mc
Lv
Ts
Og
Each element is counted once per structure.
DFT reference energies within two standard deviations of the mean. Negative values lie below the fixed MP hull.

Methodology

Convex Hull Construction in Matbench Discovery

In Matbench Discovery, the convex hull is always constructed from DFT reference energies, not from the ML model’s predicted energies. This is an important methodological choice that differs from some other benchmarking approaches and has several implications. Understanding how the convex hull is constructed is important for correctly interpreting the energy metrics in Matbench Discovery.

What This Means

Why This Matters

This approach means that:

  1. Formation energy MAE = Hull distance MAE: Because both the model’s prediction and the DFT reference are measured on the same energy scale (relative to the same elemental references), the error in formation energy prediction directly equals the error in hull distance prediction. This is a consequence of linear transformations leaving the MAE metric invariant.

  2. Systematic errors are not canceled: If a model has systematic errors (e.g., consistently over- or underpredicting certain elements), these errors will appear in both the formation energy and hull distance metrics. The model cannot “correct” for its own systematic errors by having them affect both the test structures and the reference hull equally.

    • Advantage: Tests absolute accuracy of model predictions against ground truth.
    • Use case: Pre-screening candidates for DFT calculations and evaluating a model’s ability to identify materials below the DFT reference hull.
  3. Different from some literature: Papers like Nature Communications 11:3793 (2020) construct hulls from model predictions, allowing systematic model errors to partially cancel. In that approach, the hull distance MAE can differ from the formation energy MAE.

    • Advantage: Systematic model errors can partially cancel.
    • Use case: Using the model as a complete replacement for DFT.

This distinction is subtle but important for correctly interpreting model performance and making fair comparisons between different benchmarks.