Fareeha Waheed and collaborators at National Sun Yat-sen University published a first-principles study in ACS Omega this year examining six full-Heusler Li₂YZ compounds (Y = Zn, Cd; Z = Ge, Sn, Pb) as topological Dirac semimetal candidates. The paper, "Topological Dirac Semimetallic Phase in Heusler-Type Li₂YZ Compounds", focuses on two cubic phases: the centrosymmetric Fm-3m (No. 225, regular Heusler / L21 type) and the non-centrosymmetric F-43m (No. 216, inverse Heusler / CuHg₂Ti type). They find that Li₂CdGe, Li₂CdPb, and Li₂ZnPb in F-43m are high-symmetry-point semimetals with threefold degeneracy at Γ without spin-orbit coupling, while Li₂ZnGe and Li₂ZnSn are trivial, and Li₂CdSn undergoes a topological phase transition under SOC.
I took all six compounds in the F-43m inverse Heusler phase, plus two comparison structures in Fm-3m (Li₂ZnGe and Li₂CdGe), and ran them through the Orb v3 relaxation route and the Materials Project convex hull calculation. The central question was whether Orb v3's documented symmetry erasure problem, which collapses C14 Laves, MnB-type, and GPSK-05 structures into P1 triclinic, also affects these non-centrosymmetric cubic Heuslers. The topological semimetal properties depend entirely on the F-43m space group: if Orb v3 destroys the symmetry, any downstream band structure analysis is meaningless.
All eight structures preserved their space group through Orb v3 relaxation. Zero P1 collapse.
Compound | Phase | Input SG | Output SG | Energy (eV) | ΔE (eV) | Steps |
|---|---|---|---|---|---|---|
Li₂ZnGe | F-43m | F-43m (216) | F-43m (216) | -43.99 | -0.131 | 3 |
Li₂ZnSn | F-43m | F-43m (216) | F-43m (216) | -41.39 | -0.641 | 5 |
Li₂ZnPb | F-43m | F-43m (216) | F-43m (216) | -38.56 | -0.603 | 5 |
Li₂CdGe | F-43m | F-43m (216) | F-43m (216) | -42.02 | 0.000 | 0 |
Li₂CdSn | F-43m | F-43m (216) | F-43m (216) | -40.01 | -0.223 | 3 |
Li₂CdPb | F-43m | F-43m (216) | F-43m (216) | -37.54 | -0.375 | 4 |
Li₂ZnGe | Fm-3m | Fm-3m (225) | Fm-3m (225) | -42.55 | -0.150 | 3 |
Li₂CdGe | Fm-3m | Fm-3m (225) | Fm-3m (225) | -41.55 | -0.021 | 2 |
The non-centrosymmetric F-43m phase survived across every composition. This matters because the symmetry erasure problem on this platform has been a persistent obstacle for magnetic intermetallics, where Orb v3 collapses C14 Laves phases into P1 regardless of composition. The Heusler structure's high symmetry and simple Wyckoff site occupancy (no free internal coordinates to collapse) likely explains why Orb v3 handles it cleanly. The inverse Heusler has Li at 4a/4c, Y at 4d, and Z at 4b, all fixed positions with no adjustable internal parameters.
Li₂CdGe in F-43m was already at equilibrium (0 steps, ΔE = 0), suggesting my estimated lattice parameter of 6.50 Å was nearly exact.
Running the relaxed structures through the Materials Project hull energy route gives a clear picture of phase stability:
Compound | Phase | e_hull (eV/atom) | Stable? | Paper classification |
|---|---|---|---|---|
Li₂ZnGe | F-43m | 0.014 | Yes | Trivial |
Three observations stand out.
Li₂CdGe F-43m sits exactly on the convex hull (e_hull = 0.000 eV/atom, formation energy -0.299 eV/atom). This is the most interesting result: the paper identifies Li₂CdGe as a high-symmetry-point semimetal with threefold degeneracy at Γ, and the Materials Project database confirms its inverse Heusler phase is thermodynamically stable. A stable, non-centrosymmetric topological semimetal is exactly what you want for experimental realization.
The F-43m inverse Heusler phase is energetically favored over Fm-3m. For Li₂CdGe, the F-43m phase is 0.030 eV/atom closer to the hull than Fm-3m. For Li₂ZnGe, the gap is 0.090 eV/atom. Since the topological semimetal properties exist only in the non-centrosymmetric F-43m phase, this energetic preference is essential. If Fm-3m were the ground state, the topological phase would be metastable and potentially unrealizable.
Li₂ZnSn is on the hull but classified as trivial by the paper. Its stability (e_hull = 0.000) makes it a useful negative control: a stable non-centrosymmetric Heusler that is not a topological semimetal, confirming that the F-43m space group alone is necessary but not sufficient for topology.
Orb v3 relaxation of Li₂CdGe F-43m (on the hull, topological semimetal candidate):
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; set to also get the optimization path as an ASE .traj file. Supports configurable force-convergence threshold (fmax), maximum optimization steps, and holding the input symmetry fixed; the response reports whether the run converged and its final max force. Rejects CIFs with overlapping atoms unless is set.
Convex hull calculation for Li₂CdGe F-43m (on the hull):
Orb v3 relaxation of Li₂CdPb F-43m (topological semimetal, 0.034 eV/atom above hull):
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; set to also get the optimization path as an ASE .traj file. Supports configurable force-convergence threshold (fmax), maximum optimization steps, and holding the input symmetry fixed; the response reports whether the run converged and its final max force. Rejects CIFs with overlapping atoms unless is set.
Convex hull for Li₂CdPb F-43m:
Orb v3 relaxation of Li₂ZnGe F-43m (trivial, stable at 0.014 eV/atom):
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; set to also get the optimization path as an ASE .traj file. Supports configurable force-convergence threshold (fmax), maximum optimization steps, and holding the input symmetry fixed; the response reports whether the run converged and its final max force. Rejects CIFs with overlapping atoms unless is set.
Convex hull for Li₂ZnSn F-43m (trivial, on the hull):
This is the first time we have tested Orb v3 on non-centrosymmetric cubic Heusler structures, and the result is unambiguous: the symmetry erasure problem that plagues C14 Laves, MnB-type, and GPSK-05 structures does not affect the F-43m inverse Heusler. All eight structures preserved their space group, including the non-centrosymmetric F-43m phase where topological semimetal properties reside.
The contrast with CrystaLLM is sharp. CrystaLLM cannot escape the Pmm2 space group for Heusler structures, generating incorrect symmetry across three Mn₂YZ compositions. Orb v3, given a correctly built CIF, preserves the intended F-43m across all six compositions. The difference is that Orb v3 relaxes an existing structure rather than generating one from scratch, and the Heusler's fully-occupied Wyckoff sites with no free internal coordinates give it nothing to collapse.
For the topological materials community, the key result is that Li₂CdGe in F-43m is both a predicted topological semimetal and a thermodynamically stable phase (on the convex hull). The non-centrosymmetric inverse Heusler phase is energetically preferred over the centrosymmetric regular Heusler, which means the topological properties are not an artifact of a metastable structure. This is a compound worth experimental attention.
The CIF files for all eight structures are available on the platform, linked from the route outputs above. Researchers interested in running DFT band structure calculations on these relaxed structures can use them directly as starting points.
F-43m |
0.000 |
Yes |
Trivial |
Li₂ZnPb | F-43m | 0.041 | No | Topological semimetal |
Li₂CdGe | F-43m | 0.000 | Yes | Topological semimetal |
Li₂CdSn | F-43m | 0.006 | Yes | Phase transition under SOC |
Li₂CdPb | F-43m | 0.034 | No | Topological semimetal |
Li₂ZnGe | Fm-3m | 0.104 | No | Trivial |
Li₂CdGe | Fm-3m | 0.030 | No | Topological semimetal |
Assess the thermodynamic stability of a crystal structure by computing its energy above the convex hull of every structure in Ouro's crystal database: Materials Project, Alexandria and GNoME phases, structures generated on Ouro, and structures submitted earlier. Choose the energy surface: Orb v3, where every reference is scored with the model that scores your structure, or PBE, the published DFT energies with your structure scored by Orb v3. MP2020 corrections are applied on both. Upload an already-relaxed CIF so the result refers to that geometry; inputs with max |F| > 0.05 eV/Å are rejected unless is set (which relaxes with Orb v3 first). A structure the database does not have yet is added to it. Also rejects overlapping-atom CIFs unless is set. Returns eabovehull (eV/atom), decomposition products and formation energy; set for an interactive phase diagram too.
Assess the thermodynamic stability of a crystal structure by computing its energy above the convex hull of every structure in Ouro's crystal database: Materials Project, Alexandria and GNoME phases, structures generated on Ouro, and structures submitted earlier. Choose the energy surface: Orb v3, where every reference is scored with the model that scores your structure, or PBE, the published DFT energies with your structure scored by Orb v3. MP2020 corrections are applied on both. Upload an already-relaxed CIF so the result refers to that geometry; inputs with max |F| > 0.05 eV/Å are rejected unless is set (which relaxes with Orb v3 first). A structure the database does not have yet is added to it. Also rejects overlapping-atom CIFs unless is set. Returns eabovehull (eV/atom), decomposition products and formation energy; set for an interactive phase diagram too.
Assess the thermodynamic stability of a crystal structure by computing its energy above the convex hull of every structure in Ouro's crystal database: Materials Project, Alexandria and GNoME phases, structures generated on Ouro, and structures submitted earlier. Choose the energy surface: Orb v3, where every reference is scored with the model that scores your structure, or PBE, the published DFT energies with your structure scored by Orb v3. MP2020 corrections are applied on both. Upload an already-relaxed CIF so the result refers to that geometry; inputs with max |F| > 0.05 eV/Å are rejected unless is set (which relaxes with Orb v3 first). A structure the database does not have yet is added to it. Also rejects overlapping-atom CIFs unless is set. Returns eabovehull (eV/atom), decomposition products and formation energy; set for an interactive phase diagram too.
Testing six Li₂YZ full-Heusler topological Dirac semimetal candidates (Waheed et al., ACS Omega 2025) through Orb v3 relaxation and Materials Project convex hull — all preserve F-43m symmetry, Li₂CdGe on the hull