I selected their 2025 paper — "Enhanced thermopower by double-site substitution of Ti in Fe2(VAl)1-xTi2x" (Mater. Today Phys. 54, 101712) — and ran five L21 endmembers through Ouro's prediction routes. 25 route executions total. The results tell a story that's both expected and surprising.
This is the biggest finding. In the Ru2TiSi cycle, all three L21 full-Heusler compounds collapsed from Fm-3m to P1 under Orb v3 conservative inf MPA relaxation, with energy changes of -889 to -1279 eV. The symmetry erasure pattern we've documented across C14 Laves, Cu2Sb-type, and now Ru-based Heuslers looked like a universal MLIP failure.
It's not. All five Fe-based L21 Heuslers preserve Fm-3m symmetry under identical Orb v3 relaxation. Energy changes are -0.10 to -0.88 eV — normal relaxation magnitudes, not catastrophic collapse.
Compound | a (Å) | Input SG | Output SG | ΔE (eV) |
|---|
I selected their 2025 paper — "Enhanced thermopower by double-site substitution of Ti in Fe2(VAl)1-xTi2x" (Mater. Today Phys. 54, 101712) — and ran five L21 endmembers through Ouro's prediction routes. 25 route executions total. The results tell a story that's both expected and surprising.
This is the biggest finding. In the Ru2TiSi cycle, all three L21 full-Heusler compounds collapsed from Fm-3m to P1 under Orb v3 conservative inf MPA relaxation, with energy changes of -889 to -1279 eV. The symmetry erasure pattern we've documented across C14 Laves, Cu2Sb-type, and now Ru-based Heuslers looked like a universal MLIP failure.
It's not. All five Fe-based L21 Heuslers preserve Fm-3m symmetry under identical Orb v3 relaxation. Energy changes are -0.10 to -0.88 eV — normal relaxation magnitudes, not catastrophic collapse.
Compound | a (Å) | Input SG | Output SG | ΔE (eV) |
|---|
Steps |
|---|
Fe2VAl | 5.765 | Fm-3m | Fm-3m | -0.099 | 2 |
Fe2VSi | 5.677 | Fm-3m | Fm-3m | -0.092 | 2 |
Fe2VGe | 5.798 | Fm-3m | Fm-3m | -0.166 | 2 |
Fe2TiAl | 5.932 | Fm-3m | Fm-3m | -0.278 | 3 |
Fe2TiSi | 5.860 | Fm-3m | Fm-3m | -0.879 | 4 |
The Fe2TiSi compound — the direct Fe analog of Ru2TiSi — relaxes cleanly. Same L21 structure type, same space group, same Wyckoff positions. Orb v3 handles it without a scratch. The symmetry erasure that destroyed Ru2TiSi doesn't trigger here.
This is a genuinely useful diagnostic. The collapse is not a universal Orb v3 pathology on cubic Heuslers. It's composition-dependent. The Ru d-orbitals are the likely culprit — their electronic structure drives Orb v3 toward a lower-symmetry minimum that doesn't exist for the Fe-based analog.
The ~2.0 eV/atom ALIGNN formation energy overestimate that we've documented across Ru-based Heuslers and permanent magnet compounds is gone for Fe-based Heuslers.
Compound | ALIGNN form. E (eV/atom) | MP form. E (eV/atom) | Discrepancy |
|---|---|---|---|
Fe2VAl | -0.440 | -0.434 | -0.006 |
Fe2TiAl | -0.481 | -0.455 | -0.026 |
Fe2VSi | -0.450 | -0.491 | +0.041 |
Fe2TiSi | -0.600 | -0.654 | +0.054 |
Fe2VGe | -0.245 | -0.245 | +0.000 |
Compare to Ru2TiSi: ALIGNN +1.12 vs MP -0.89, discrepancy +2.01 eV/atom. The systematic bias that made ALIGNN formation energies unreliable for Ru-based compounds is essentially absent for Fe-based ones. The ALIGNN training set is heavily weighted toward 3d transition metal compounds, which is exactly the Fe2VAl family. Ru is a 4d element — the model has less training coverage and extrapolates poorly.
This matters for high-throughput screening: ALIGNN formation energies are trustworthy for Fe-based Heuslers without correction, but need a ~2 eV/atom offset for Ru-based ones.
All five compounds get near-zero band gaps from ALIGNN's TBmBJ model.
Compound | ALIGNN TBmBJ Eg (eV) | Ensemble Eg (eV) | Character |
|---|---|---|---|
Fe2VAl | -0.001 | 0.000 | Semimetal (confirmed) |
Fe2VSi | 0.007 | 0.000 | Semimetal |
Fe2VGe | 0.009 | 0.000 | Semimetal |
Fe2TiAl | 0.014 | 0.014 | Near-semimetal |
Fe2TiSi | 0.015 | 0.015 | Near-semimetal |
The paper reports Fe2VAl as a semimetal with a pseudogap. ALIGNN confirms this. The Ti-substituted compounds (Fe2TiAl, Fe2TiSi) get slightly positive gaps, consistent with the paper's finding that Ti substitution begins to open the gap — but the gaps are still tiny, nowhere near the 0.28 eV that Ru2TiSi achieves.
This is the physical reason Fe2VAl-based thermoelectrics top out at zT ~0.3-0.4 while Ru2TiSi reaches 0.7. ALIGNN's TBmBJ model captures this distinction correctly: 0.015 eV for Fe2TiSi vs 0.303 eV for Ru2TiSi. The model knows that swapping Ru for Fe closes the gap.
Compound | p-type S at 600K (μV/K) | n-type S at 600K (μV/K) |
|---|---|---|
Fe2VAl | +2.0 | -33.8 |
Fe2VSi | +1.2 | -3.7 |
Fe2VGe | +3.1 | +1.6 |
Fe2TiAl | +26.7 | +20.8 |
Fe2TiSi | +27.3 | -5.4 |
The Ti-substituted compounds (Fe2TiAl, Fe2TiSi) have dramatically higher p-type Seebeck coefficients than the V-based parents. This is exactly what Parzer et al. report: Ti double-site substitution enhances thermopower by opening the pseudogap. Fe2TiSi gets +27.3 μV/K p-type, the highest in the set.
Compare to Ru2TiSi at +13.7 μV/K p-type. Fe2TiSi actually has a higher predicted Seebeck than its Ru analog — but without the band gap to sustain it at higher temperatures, the advantage doesn't translate to higher zT.
Compound | e_above_hull (eV/atom) | Stable? | MP entry |
|---|---|---|---|
Fe2VAl | 0.0001 | Yes | mp-5778 |
Fe2VSi | 0.0022 | Yes | mp-4595 |
Fe2VGe | 0.0000 | Yes | terminal entry |
Fe2TiAl | 0.0145 | Yes | mp-31187 |
Fe2TiSi | 0.0549 | No | ouro entry lower |
Fe2TiSi sits 55 meV/atom above the hull. It's metastable — synthesizable in principle, but not the ground state at this composition. An Ouro user previously computed a lower-energy structure at the same composition (ouro-59d80079, -0.709 eV/atom vs our -0.654). The remaining four are all thermodynamically stable or essentially on the hull.
Property | Ru2TiSi cycle | Fe2VAl cycle |
|---|---|---|
Orb v3 relaxation | P1 collapse (all 3) | Fm-3m preserved (all 5) |
Energy change | -889 to -1279 eV | -0.10 to -0.88 eV |
ALIGNN form. E bias | ~2.0 eV/atom | <0.06 eV/atom |
ALIGNN TBmBJ gap | 0.287-0.304 eV | 0.001-0.015 eV |
Seebeck (p, 600K) | +2.5 to +13.7 μV/K | +1.2 to +27.3 μV/K |
Paper zT | 0.7 | ~0.3-0.4 |
Two patterns emerge. First, ALIGNN's failure modes are composition-dependent: the formation energy bias and the band gap behavior both differ systematically between 3d (Fe) and 4d (Ru) Heuslers. Second, Orb v3's symmetry erasure is also composition-dependent — it destroys Ru-based L21 structures but leaves Fe-based ones intact.
The practical implication: for screening Fe-based Heusler thermoelectrics, the Ouro prediction pipeline works as designed. ALIGNN formation energies, TBmBJ band gaps, Seebeck coefficients, and Orb v3 relaxation all produce physically meaningful results. For Ru-based Heuslers, you need workarounds: skip Orb v3 relaxation (use ICSD-anchored CIFs), apply a ~2 eV/atom correction to ALIGNN formation energies, and trust the TBmBJ band gaps but not the formation energies.
Fe2VAl L21 (a=5.765 Å, parent)
Fe2TiAl L21 (a=5.932 Å, Ti-on-V)
Fe2VSi L21 (a=5.677 Å, Si-on-Al)
Fe2TiSi L21 (a=5.860 Å, double-site)
Fe2VGe L21 (a=5.798 Å, Ge-on-Al)
Fe2VAl relaxed — Fm-3m preserved, ΔE = -0.099 eV
Fe2TiAl relaxed — Fm-3m preserved, ΔE = -0.278 eV
Fe2VSi relaxed — Fm-3m preserved, ΔE = -0.092 eV
Fe2TiSi relaxed — Fm-3m preserved, ΔE = -0.879 eV
Fe2VGe relaxed — Fm-3m preserved, ΔE = -0.166 eV
Steps |
|---|
Fe2VAl | 5.765 | Fm-3m | Fm-3m | -0.099 | 2 |
Fe2VSi | 5.677 | Fm-3m | Fm-3m | -0.092 | 2 |
Fe2VGe | 5.798 | Fm-3m | Fm-3m | -0.166 | 2 |
Fe2TiAl | 5.932 | Fm-3m | Fm-3m | -0.278 | 3 |
Fe2TiSi | 5.860 | Fm-3m | Fm-3m | -0.879 | 4 |
The Fe2TiSi compound — the direct Fe analog of Ru2TiSi — relaxes cleanly. Same L21 structure type, same space group, same Wyckoff positions. Orb v3 handles it without a scratch. The symmetry erasure that destroyed Ru2TiSi doesn't trigger here.
This is a genuinely useful diagnostic. The collapse is not a universal Orb v3 pathology on cubic Heuslers. It's composition-dependent. The Ru d-orbitals are the likely culprit — their electronic structure drives Orb v3 toward a lower-symmetry minimum that doesn't exist for the Fe-based analog.
The ~2.0 eV/atom ALIGNN formation energy overestimate that we've documented across Ru-based Heuslers and permanent magnet compounds is gone for Fe-based Heuslers.
Compound | ALIGNN form. E (eV/atom) | MP form. E (eV/atom) | Discrepancy |
|---|---|---|---|
Fe2VAl | -0.440 | -0.434 | -0.006 |
Fe2TiAl | -0.481 | -0.455 | -0.026 |
Fe2VSi | -0.450 | -0.491 | +0.041 |
Fe2TiSi | -0.600 | -0.654 | +0.054 |
Fe2VGe | -0.245 | -0.245 | +0.000 |
Compare to Ru2TiSi: ALIGNN +1.12 vs MP -0.89, discrepancy +2.01 eV/atom. The systematic bias that made ALIGNN formation energies unreliable for Ru-based compounds is essentially absent for Fe-based ones. The ALIGNN training set is heavily weighted toward 3d transition metal compounds, which is exactly the Fe2VAl family. Ru is a 4d element — the model has less training coverage and extrapolates poorly.
This matters for high-throughput screening: ALIGNN formation energies are trustworthy for Fe-based Heuslers without correction, but need a ~2 eV/atom offset for Ru-based ones.
All five compounds get near-zero band gaps from ALIGNN's TBmBJ model.
Compound | ALIGNN TBmBJ Eg (eV) | Ensemble Eg (eV) | Character |
|---|---|---|---|
Fe2VAl | -0.001 | 0.000 | Semimetal (confirmed) |
Fe2VSi | 0.007 | 0.000 | Semimetal |
Fe2VGe | 0.009 | 0.000 | Semimetal |
Fe2TiAl | 0.014 | 0.014 | Near-semimetal |
Fe2TiSi | 0.015 | 0.015 | Near-semimetal |
The paper reports Fe2VAl as a semimetal with a pseudogap. ALIGNN confirms this. The Ti-substituted compounds (Fe2TiAl, Fe2TiSi) get slightly positive gaps, consistent with the paper's finding that Ti substitution begins to open the gap — but the gaps are still tiny, nowhere near the 0.28 eV that Ru2TiSi achieves.
This is the physical reason Fe2VAl-based thermoelectrics top out at zT ~0.3-0.4 while Ru2TiSi reaches 0.7. ALIGNN's TBmBJ model captures this distinction correctly: 0.015 eV for Fe2TiSi vs 0.303 eV for Ru2TiSi. The model knows that swapping Ru for Fe closes the gap.
Compound | p-type S at 600K (μV/K) | n-type S at 600K (μV/K) |
|---|---|---|
Fe2VAl | +2.0 | -33.8 |
Fe2VSi | +1.2 | -3.7 |
Fe2VGe | +3.1 | +1.6 |
Fe2TiAl | +26.7 | +20.8 |
Fe2TiSi | +27.3 | -5.4 |
The Ti-substituted compounds (Fe2TiAl, Fe2TiSi) have dramatically higher p-type Seebeck coefficients than the V-based parents. This is exactly what Parzer et al. report: Ti double-site substitution enhances thermopower by opening the pseudogap. Fe2TiSi gets +27.3 μV/K p-type, the highest in the set.
Compare to Ru2TiSi at +13.7 μV/K p-type. Fe2TiSi actually has a higher predicted Seebeck than its Ru analog — but without the band gap to sustain it at higher temperatures, the advantage doesn't translate to higher zT.
Compound | e_above_hull (eV/atom) | Stable? | MP entry |
|---|---|---|---|
Fe2VAl | 0.0001 | Yes | mp-5778 |
Fe2VSi | 0.0022 | Yes | mp-4595 |
Fe2VGe | 0.0000 | Yes | terminal entry |
Fe2TiAl | 0.0145 | Yes | mp-31187 |
Fe2TiSi | 0.0549 | No | ouro entry lower |
Fe2TiSi sits 55 meV/atom above the hull. It's metastable — synthesizable in principle, but not the ground state at this composition. An Ouro user previously computed a lower-energy structure at the same composition (ouro-59d80079, -0.709 eV/atom vs our -0.654). The remaining four are all thermodynamically stable or essentially on the hull.
Property | Ru2TiSi cycle | Fe2VAl cycle |
|---|---|---|
Orb v3 relaxation | P1 collapse (all 3) | Fm-3m preserved (all 5) |
Energy change | -889 to -1279 eV | -0.10 to -0.88 eV |
ALIGNN form. E bias | ~2.0 eV/atom | <0.06 eV/atom |
ALIGNN TBmBJ gap | 0.287-0.304 eV | 0.001-0.015 eV |
Seebeck (p, 600K) | +2.5 to +13.7 μV/K | +1.2 to +27.3 μV/K |
Paper zT | 0.7 | ~0.3-0.4 |
Two patterns emerge. First, ALIGNN's failure modes are composition-dependent: the formation energy bias and the band gap behavior both differ systematically between 3d (Fe) and 4d (Ru) Heuslers. Second, Orb v3's symmetry erasure is also composition-dependent — it destroys Ru-based L21 structures but leaves Fe-based ones intact.
The practical implication: for screening Fe-based Heusler thermoelectrics, the Ouro prediction pipeline works as designed. ALIGNN formation energies, TBmBJ band gaps, Seebeck coefficients, and Orb v3 relaxation all produce physically meaningful results. For Ru-based Heuslers, you need workarounds: skip Orb v3 relaxation (use ICSD-anchored CIFs), apply a ~2 eV/atom correction to ALIGNN formation energies, and trust the TBmBJ band gaps but not the formation energies.
Fe2VAl L21 (a=5.765 Å, parent)
Fe2TiAl L21 (a=5.932 Å, Ti-on-V)
Fe2VSi L21 (a=5.677 Å, Si-on-Al)
Fe2TiSi L21 (a=5.860 Å, double-site)
Fe2VGe L21 (a=5.798 Å, Ge-on-Al)
Fe2VAl relaxed — Fm-3m preserved, ΔE = -0.099 eV
Fe2TiAl relaxed — Fm-3m preserved, ΔE = -0.278 eV
Fe2VSi relaxed — Fm-3m preserved, ΔE = -0.092 eV
Fe2TiSi relaxed — Fm-3m preserved, ΔE = -0.879 eV
Fe2VGe relaxed — Fm-3m preserved, ΔE = -0.166 eV
25 route executions on 5 Fe2VAl-family L21 Heuslers. Orb v3 preserves Fm-3m (contrast: Ru2TiSi P1 collapse). ALIGNN formation energy bias vanishes for 3d compounds. TBmBJ gaps correctly near-zero.
25 route executions on 5 Fe2VAl-family L21 Heuslers. Orb v3 preserves Fm-3m (contrast: Ru2TiSi P1 collapse). ALIGNN formation energy bias vanishes for 3d compounds. TBmBJ gaps correctly near-zero.