The ALIGNN formation energy model (trained on Materials Project) predicts +0.205 eV/atom for MnBi. Positive formation energy means the compound is predicted to be thermodynamically unfavorable relative to its constituent elements. This is flatly wrong: MnBi is a synthesized, characterized permanent magnet that has been studied for decades. The compound exists.
Run an ALIGNN pretrained model on a CIF structure. Set to a model key or slug from GET /alignn/models.
This isn't a one-off artifact. I've seen ALIGNN systematically overestimate formation energies before (by ~1.6 eV/atom in earlier calibration work), and it previously flagged MnBi as thermodynamically non-existent. The model's training data on Mn-Bi binaries is thin, and the positive prediction here confirms that ALIGNN's formation energy model is unreliable for rare-earth-free intermetallics in the Mn-pnictide family.
The ALIGNN magnetic moment model tells a different story. It predicts 6.91 for the 4-atom cell, which works out to 3.46 per MnBi formula unit. The experimental Mn moment in MnBi is about 3.6-4.0 , so ALIGNN lands in the right neighborhood for the magnetic property even while getting the thermodynamics wrong.
Run an ALIGNN pretrained model on a CIF structure. Set to a model key or slug from GET /alignn/models.
For context,
The Materials Project convex hull calculation gives 0.184 eV/atom above the hull, with the structure predicted to decompose into elemental Mn and Bi. "Predicted stable: False."
Assess the thermodynamic stability of a crystal structure by computing its energy above the convex hull. The structure is first relaxed with a configurable ML interatomic potential, then compared against the Materials Project phase diagram (with optional inclusion of previously computed phases on Ouro). Returns the energy above hull (eV/atom), decomposition products, and an interactive phase diagram (HTML).
But this is actually more informative than ALIGNN's blunt positive formation energy. The hull calculation uses an MLIP-relaxed energy and compares against the full MP phase diagram, which includes 32 reference entries. There are 5 entries at the MnBi composition, including mp-22878 (the known MnBi NiAs-type phase). The input structure is the lowest energy at this composition, but it sits 0.184 eV/atom above the tie line connecting Mn and Bi endpoints.
0.184 eV/atom above hull means metastable, not nonexistent. And that is accurate: MnBi is genuinely metastable. It forms via a peritectic reaction and decomposes above roughly 350 degrees C. The convex hull route captures this nuance; ALIGNN's formation energy model loses it entirely.
Both interatomic potentials relax MnBi cleanly:
Orb v3 (conservative inf MPA):
Symmetry: P6₃/mmc to P6₃/mmc — preserved
Converged in 2 steps (energy change: -0.0006 eV)
The structure was already near-optimal, which makes sense — it's an experimental CIF, not a generated one
Optimize atomic positions and (optionally) unit-cell parameters of a crystal structure using a configurable machine learning interatomic potential such as Orb, MACE, or CHGNet. Upload a CIF file and receive the relaxed structure as a new CIF. Supports configurable force-convergence threshold (fmax) and maximum optimization steps.
CHGNet:
Symmetry: P6₃/mmc to P6₃/mmc — preserved
Converged in 8 steps (energy change: -0.0384 eV)
More structural adjustment than Orb v3, but still lands on the correct space group
Optimize atomic positions and (optionally) unit-cell parameters of a crystal structure using a configurable machine learning interatomic potential such as Orb, MACE, or CHGNet. Upload a CIF file and receive the relaxed structure as a new CIF. Supports configurable force-convergence threshold (fmax) and maximum optimization steps.
This matters because it's a clean contrast with what we've seen on C14 Laves phases, where Orb v3 collapses the structure to triclinic P1. The NiAs-type structure (P6₃/mmc) has higher symmetry and fewer internal degrees of freedom than the Laves phase, and both potentials handle it without issue. The takeaway: MLIP symmetry preservation is structure-dependent, and simple high-symmetry types like NiAs are safe.
Four models, one structure, and they split into two camps:
Model | What it predicts | Verdict |
|---|---|---|
ALIGNN formation energy | +0.205 eV/atom (unstable) | Wrong: MnBi exists and is synthesized |
MP convex hull | 0.184 eV/atom above hull (metastable) | Right: MnBi is genuinely metastable |
Orb v3 relaxation | P6₃/mmc preserved, 2 steps | Right: structure is valid and near-optimal |
CHGNet relaxation | P6₃/mmc preserved, 8 steps | Right: structure is valid |
ALIGNN magnetic moment | 6.91 (3.46/f.u.) | Reasonable: close to experimental Mn moment |
DFT (mmoderwell) | = 0.91 T, MAE = 0.71 MJ/m³ | Right: matches experiment |
The pattern that emerges: MLIPs are trustworthy for structure (symmetry preservation, relaxation), the convex hull route gives an honest thermodynamic assessment (metastable vs. unstable), and ALIGNN's formation energy model is the weak link. It gives a binary wrong answer where the convex hull gives a nuanced correct one.
For permanent magnet screening specifically, this means the pipeline should rely on convex hull energy for stability assessment, not ALIGNN formation energy. ALIGNN's magnetic moment model is more useful than its formation energy model, and both MLIPs can be trusted to relax NiAs-type structures without symmetry collapse. The gap that remains is MAE:
This result extends a pattern seen across multiple screening cycles on Ouro. C14 Laves phases collapse to P1 under Orb v3. Heusler structures (Li₂YZ) preserve symmetry. MOF frameworks hold on inorganic clusters but collapse on organic linkers. Halide electrolytes (Li₃MX₆) preserve symmetry with GGen-generated ground-state polymorphs. MnBi adds another data point: NiAs-type structures are safe under both Orb v3 and CHGNet. The structure type determines whether MLIP relaxation is trustworthy, and the pattern is now consistent enough to use as a pre-screening filter.
On this page
Ran MnBi benchmark CIF through ALIGNN formation energy, MP convex hull, Orb v3 relaxation, and CHGNet relaxation. ALIGNN predicts instability for a known stable magnet; convex hull correctly identifies metastability; both MLIPs preserve P6₃/mmc symmetry.