pith. machine review for the scientific record. sign in

arxiv: 2604.27685 · v1 · submitted 2026-04-30 · ❄️ cond-mat.mtrl-sci · cs.AI· cs.LG· physics.comp-ph

Recognition: unknown

VibroML: an automated toolkit for high-throughput vibrational analysis and dynamic instability remediation of crystalline materials using machine-learned potentials

Gian-Marco Rignanese, Rog\'erio Almeida Gouv\^ea

Authors on Pith no claims yet

Pith reviewed 2026-05-07 06:25 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci cs.AIcs.LGphysics.comp-ph
keywords machine-learned potentialsgenetic algorithmdynamical stabilitycrystal polymorphshigh-throughput screeningstructural remediationperovskitesvibrational analysis
0
0 comments X

The pith

VibroML uses an energy-guided genetic algorithm on machine-learned potentials to automatically remediate dynamical instabilities and uncover stable crystal polymorphs.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces VibroML as a toolkit that shifts focus from simply checking for vibrational instabilities in computationally predicted crystals to actively fixing them. It drives an energy-guided genetic algorithm with foundational machine-learned interatomic potentials to explore the potential energy surface and locate diverse dynamically stable structures, outperforming traditional soft-mode methods. The workflow adds automated molecular dynamics runs to test whether candidate structures hold up at finite temperatures rather than only at zero kelvin. It also links to combinatorial structure prediction for alloying to stabilize otherwise frustrated crystal networks, such as certain perovskites. When applied to large databases with many missing or unstable high-symmetry entries, the approach yields thousands of lower-symmetry candidates that standard high-throughput pipelines would discard.

Core claim

VibroML couples machine-learned interatomic potentials with an energy-guided genetic algorithm that navigates the potential energy surface to identify diverse dynamically stable polymorphs for high-symmetry structures showing instabilities, includes molecular dynamics to confirm finite-temperature retention, and integrates with combinatorial prediction to stabilize frustrated topologies through targeted alloying, as shown by rescuing perovskite networks and mining databases for abandoned quaternary and quinary stoichiometries.

What carries the argument

Energy-guided genetic algorithm driven by machine-learned interatomic potentials that explores the potential energy surface to locate diverse dynamically stable polymorphs.

If this is right

  • Enables stabilization of functional perovskite networks such as Cs2KInI6 and KTaSe3 through targeted alloying.
  • Recovers dynamically stable low-symmetry candidates from databases where roughly half of quaternary and nearly all quinary elemental combinations lack any structural entries.
  • Moves high-throughput workflows from instability detection alone to full structural remediation plus thermal validation.
  • Supports systematic exploration of complex compositions that standard pipelines leave abandoned.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The method could lower the fraction of computationally generated materials that are discarded solely because of initial instabilities, opening routes to more candidates for applications.
  • Pairing the remediation step with additional prediction engines might help close coverage gaps in materials databases for multi-element systems.
  • The emphasis on finite-temperature checks suggests that purely harmonic zero-kelvin analyses may systematically miss viable structures in future screening campaigns.

Load-bearing premise

The machine-learned potentials must represent the true potential energy surface accurately enough for the genetic algorithm to reach genuine stable polymorphs instead of metastable or spurious minima.

What would settle it

Recomputing the phonon spectra and running extended molecular dynamics on the remediated polymorphs with density functional theory, then checking against any available experimental stability data, would confirm or refute whether the structures are truly stable.

Figures

Figures reproduced from arXiv: 2604.27685 by Gian-Marco Rignanese, Rog\'erio Almeida Gouv\^ea.

Figure 1
Figure 1. Figure 1: Evolution of dynamical stability across the view at source ↗
Figure 2
Figure 2. Figure 2: Thermodynamic stability analysis and phase evolution of the view at source ↗
Figure 3
Figure 3. Figure 3: Dynamical stability map for the targeted cubic view at source ↗
Figure 4
Figure 4. Figure 4: Thermodynamic stability analysis and structural evolution of the view at source ↗
read the original abstract

While machine-learned interatomic potentials (MLIPs) accelerate phonon dispersion calculations, merely identifying dynamical instabilities in computationally predicted materials is insufficient; automated pathways to resolve them are required. We introduce VibroML, an open-source Python toolkit driven by foundational MLIPs that shifts the paradigm from stability verification to automated structural remediation. VibroML employs an energy-guided genetic algorithm that vastly outperforms traditional soft-mode following, efficiently navigating the potential energy surface to uncover diverse, dynamically stable polymorphs. As 0 K harmonic stability does not guarantee macroscopic viability, an automated molecular dynamics workflow evaluates finite-temperature structural retention. VibroML also couples with ProtoCSP, our combinatorial structure prediction engine, to stabilize frustrated crystal topologies via targeted alloying, successfully rescuing functional perovskite networks like Cs$_2$KInI$_6$ and KTaSe$_3$. Demonstrating broader applicability, we mined the Alexandria database -- where ~50% of quaternary and 99.5% of quinary elemental combinations lack any structural entries -- to identify thousands of abandoned, high-symmetry stoichiometries. Deploying ProtoCSP's "cold start" retrieval and VibroML's evolutionary search on a sample, we successfully identified dynamically stable low-symmetry candidates. Through integrated structural remediation, thermal validation, and systematic compositional exploration, VibroML enables a comprehensive deep-screening approach, yielding physically sound structural propositions that far surpass standard high-throughput workflows.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript introduces VibroML, an open-source Python toolkit for high-throughput vibrational analysis and dynamic instability remediation in crystalline materials using foundational machine-learned interatomic potentials (MLIPs). It employs an energy-guided genetic algorithm claimed to outperform traditional soft-mode following for discovering diverse, dynamically stable polymorphs, incorporates automated molecular dynamics workflows for finite-temperature structural retention, and integrates with ProtoCSP for targeted alloying to stabilize frustrated topologies. Applications include remediation of specific perovskites (Cs₂KInI₆ and KTaSe₃) and mining the Alexandria database to identify stable low-symmetry candidates from high-symmetry stoichiometries lacking prior entries.

Significance. If the central claims hold with proper validation, VibroML could represent a practical advance in automated materials screening by shifting focus from instability detection to remediation, potentially improving efficiency in high-throughput workflows. Strengths include the open-source toolkit, integration of MLIPs with evolutionary search and MD validation, and demonstration on real databases. However, the overall significance hinges on demonstrating that MLIP-guided structures are not artifacts and remain stable under higher-accuracy methods.

major comments (3)
  1. [Results on specific material remediation] The claims of successful remediation for Cs₂KInI₆ and KTaSe₃ (abstract and results) provide no quantitative metrics such as MLIP vs. DFT energy differences, phonon frequency agreement, or false-positive rates for the final structures. This is load-bearing because MLIP inaccuracies in ionic or van der Waals interactions could produce spurious minima that appear stable only under the surrogate potential.
  2. [Methods and results on genetic algorithm performance] The assertion that the energy-guided genetic algorithm 'vastly outperforms' soft-mode following lacks any benchmarks, success rates, computational cost comparisons, or validation protocols on a test set (methods and results sections). Without these, the central performance claim cannot be assessed.
  3. [Application to Alexandria database] In the Alexandria database application, the identification of 'thousands of abandoned, high-symmetry stoichiometries' and 'dynamically stable low-symmetry candidates' is stated without specific examples, energy-above-hull values, or DFT cross-checks on the output structures. This undermines the broader applicability demonstration.
minor comments (2)
  1. [Abstract] The abstract would benefit from naming the specific foundational MLIPs employed and their training data sources to allow readers to assess transferability.
  2. [Methods on MD validation] Clarify in the methods whether the MD workflow includes anharmonic effects or merely checks retention of the 0 K structure, and ensure consistent terminology for 'dynamically stable' across sections.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their thorough review and valuable feedback on our manuscript. We address each of the major comments below and outline the revisions we will make to strengthen the paper.

read point-by-point responses
  1. Referee: [Results on specific material remediation] The claims of successful remediation for Cs₂KInI₆ and KTaSe₃ (abstract and results) provide no quantitative metrics such as MLIP vs. DFT energy differences, phonon frequency agreement, or false-positive rates for the final structures. This is load-bearing because MLIP inaccuracies in ionic or van der Waals interactions could produce spurious minima that appear stable only under the surrogate potential.

    Authors: We agree that additional quantitative validation is essential to substantiate the remediation claims. In the revised manuscript, we will include a new table or section detailing the MLIP-predicted energies compared to DFT calculations for the final structures of Cs₂KInI₆ and KTaSe₃. We will also report phonon frequency comparisons for these structures. Regarding false-positive rates, we will discuss our multi-stage validation protocol, including finite-temperature MD simulations, and acknowledge the limitations of MLIPs in capturing certain interactions. We will perform DFT cross-validations on the key remediated structures to address concerns about spurious minima. revision: yes

  2. Referee: [Methods and results on genetic algorithm performance] The assertion that the energy-guided genetic algorithm 'vastly outperforms' soft-mode following lacks any benchmarks, success rates, computational cost comparisons, or validation protocols on a test set (methods and results sections). Without these, the central performance claim cannot be assessed.

    Authors: We acknowledge that the performance claim requires more rigorous benchmarking. We will add a dedicated subsection in the Methods and Results sections that includes benchmarks on a curated test set of materials with known instabilities. This will report success rates in finding stable polymorphs, computational costs (e.g., number of energy evaluations), and direct comparisons to soft-mode following methods. The test set will consist of small-to-medium sized systems where DFT is feasible for validation. This will allow readers to quantitatively assess the advantages of the energy-guided GA. revision: yes

  3. Referee: [Application to Alexandria database] In the Alexandria database application, the identification of 'thousands of abandoned, high-symmetry stoichiometries' and 'dynamically stable low-symmetry candidates' is stated without specific examples, energy-above-hull values, or DFT cross-checks on the output structures. This undermines the broader applicability demonstration.

    Authors: We agree that providing concrete examples and validation metrics would enhance the demonstration of broader applicability. In the revised manuscript, we will include specific examples of the identified low-symmetry candidates from the Alexandria database sample, along with their MLIP-computed energy-above-hull values. Additionally, we will perform and report DFT calculations for a subset of these candidates to provide cross-checks. This will be presented in a new figure or table in the Results section, illustrating the remediation process for these cases. revision: yes

Circularity Check

0 steps flagged

No circularity in VibroML toolkit description or claims

full rationale

The manuscript describes an engineering toolkit (VibroML) that combines existing MLIPs with a genetic algorithm for structural remediation, MD validation, and coupling to the authors' prior ProtoCSP engine. No equations, fitted parameters, or predictions are presented that reduce to inputs by construction. Claims of outperforming soft-mode following rest on reported applications to specific compounds and the Alexandria database rather than any self-referential derivation. Self-reference to ProtoCSP is acknowledged but is not load-bearing for any mathematical result; the toolkit's outputs are framed as externally validated via MD and phonon checks. This is a standard methods contribution with no detectable circular steps in its derivation chain.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on the domain assumption that pre-trained machine-learned interatomic potentials are accurate enough for reliable phonon calculations and energy-guided structural search; no free parameters or invented entities are explicitly introduced in the abstract.

axioms (1)
  • domain assumption Machine-learned interatomic potentials provide a sufficiently accurate representation of the potential energy surface for the crystalline materials under study.
    All vibrational analysis, genetic algorithm guidance, and stability assessments rely on these potentials as the computational engine.

pith-pipeline@v0.9.0 · 5578 in / 1338 out tokens · 78614 ms · 2026-05-07T06:25:31.095138+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

66 extracted references · 16 canonical work pages

  1. [1]

    Shahzad, K., Mardare, A. I. & Hassel, A. W. Accelerating materials discovery: Combinatorial synthesis, high-throughput characterization, and computational advances.Science and Technology of Advanced Materials: Methods4, 2292486 (2024). 26

  2. [2]

    L., Takeuchi, I

    Joress, H., Green, M. L., Takeuchi, I. & Hattrick-Simpers, J. R. Applications of High Throughput (Combinatorial) Methodologies to Electronic, Magnetic, Structural, and Energy- Related Materials. InEncyclopedia of Materials: Metals and Alloys, 353–371 (Elsevier, 2022)

  3. [3]

    Chen, C.et al.Accelerating Computational Materials Discovery with Machine Learning and Cloud High-Performance Computing: From Large-Scale Screening to Experimental Validation.Journal of the American Chemical Society146, 20009–20018 (2024)

  4. [4]

    Bartel, C. J. Review of computational approaches to predict the thermodynamic stability of inorganic solids.Journal of Materials Science57, 10475–10498 (2022)

  5. [5]

    Rudin, S. P. Generalization of soft phonon modes.Physical Review B97, 134114 (2018)

  6. [6]

    & Tadano, T

    Xiao, E. & Tadano, T. High-throughput computational screening of Heusler compounds with phonon considerations for enhanced material discovery.arXiv preprint(2025). 2502.17946

  7. [7]

    D., Xia, Y

    Griesemer, S. D., Xia, Y. & Wolverton, C. Accelerating the prediction of stable materials with machine learning.Nature Computational Science3, 934–945 (2023)

  8. [8]

    & Mizukami, W

    Shiota, T., Ishihara, K. & Mizukami, W. Universal neural network potentials as descriptors: Towards scalable chemical property prediction using quantum and classical computers. Digital Discovery3, 1714–1728 (2024)

  9. [9]

    Wang, J., Huang, Y., Gao, Z. & Du, J. Lead-Free Antimony-Based Perovskite Nanocrystals as a Fluorescent Nanoprobe for the Turn-On Detection of Trace Water in Organic Solvents. ACS Applied Nano Materials7, 7958–7965 (2024)

  10. [10]

    Batzner, S.et al.E(3)-equivariant graph neural networks for data-efficient and accurate interatomic potentials.Nature Communications13, 2453 (2022)

  11. [11]

    P., Simm, G

    Batatia, I., Kov´ acs, D. P., Simm, G. N. C., Ortner, C. & Cs´ anyi, G. MACE: Higher Order Equivariant Message Passing Neural Networks for Fast and Accurate Force Fields (2022)

  12. [12]

    & Jung, Y

    Choi, J., Nam, G., Choi, J. & Jung, Y. A Perspective on Foundation Models in Chemistry. JACS Au5, 1499–1518 (2025)

  13. [13]

    Musaelian, A.et al.Learning local equivariant representations for large-scale atomistic dynamics.Nature Communications14, 579 (2023)

  14. [14]

    Fu, X.et al.Learning Smooth and Expressive Interatomic Potentials for Physical Property Prediction.arXiv preprint(2025).2502.12147

  15. [15]

    Wood, Misko Dzamba, Xiang Fu, Meng Gao, Muhammed Shuaibi, Luis Barroso-Luque, Kareem Abdelmaqsoud, Vahe Gharakhanyan, John R

    Wood, B. M.et al.UMA: A Family of Universal Models for Atoms.arXiv preprint(2025). 2506.23971

  16. [16]

    S., Ghorbani, K., Hatam-Lee, S

    Arabha, S., Aghbolagh, Z. S., Ghorbani, K., Hatam-Lee, S. M. & Rajabpour, A. Recent advances in lattice thermal conductivity calculation using machine-learning interatomic potentials.Journal of Applied Physics130, 210903 (2021)

  17. [17]

    & Marques, M

    Loew, A., Sun, D., Wang, H.-C., Botti, S. & Marques, M. A. L. Universal machine learning interatomic potentials are ready for phonons.npj Computational Materials11, 178 (2025)

  18. [18]

    & Cao, B

    Guo, L., Liu, Y., Yang, L. & Cao, B. Lattice dynamics modeling of thermal transport in solids using machine-learned atomic cluster expansion potentials: A tutorial.Journal of Applied Physics137, 081101 (2025). 27

  19. [19]

    Batatia, I.et al.A foundation model for atomistic materials chemistry.arXiv preprint (2025).2401.00096

  20. [20]

    & Cheng, Y

    Han, B. & Cheng, Y. Benchmarking universal machine learning interatomic potentials for rapid analysis of inelastic neutron scattering data.Machine Learning: Science and Technology6, 030504 (2025)

  21. [21]

    G., Allen, C

    Radova, M., Stark, W. G., Allen, C. S., Maurer, R. J. & Bart´ ok, A. P. Fine-tuning foundation models of materials interatomic potentials with frozen transfer learning.arXiv preprint(2025).2502.15582

  22. [22]

    Liu, X.et al.Fine-Tuning Universal Machine-Learned Interatomic Potentials: A Tutorial on Methods and Applications.arXiv preprint(2025).2506.21935

  23. [23]

    S., Trybel, F., Faber, F

    Parackal, A. S., Trybel, F., Faber, F. A. & Armiento, R. Screening 39 billion protostructures for materials discovery.arXiv preprint(2026).2601.21393

  24. [24]

    Hoellmer, P.et al.Open Materials Generation with Stochastic Interpolants.arXiv preprint (2025).2502.02582

  25. [25]

    & Lu, Z.-Y

    Han, X.-Q., Guo, P.-J., Gao, Z.-F. & Lu, Z.-Y. PhononBench:A Large-Scale Phonon- Based Benchmark for Dynamical Stability in Crystal Generation.arXiv preprint(2025). 2512.21227

  26. [26]

    & Tanaka, I

    Togo, A. & Tanaka, I. Evolution of crystal structures in metallic elements.Physical Review B87, 184104 (2013)

  27. [27]

    P., Simm, G., Ortner, C

    Batatia, I., Kov´ acs, D. P., Simm, G., Ortner, C. & Cs´ anyi, G. MACE foundations: Pre-trained foundation models for atomistic materials chemistry. https://github.com/ ACEsuit/mace-foundations(2024). Accessed: 2025-02-10

  28. [28]

    & Rabe, K

    Qi, Y. & Rabe, K. M. Competing Phases of HfO 2 from Unstable Flat Phonon Bands of an Unconventional High-Symmetry Structure.Physical Review Letters135, 046101 (2025)

  29. [29]

    Rahim, W.et al.Polymorph exploration of bismuth stannate using first-principles phonon mode mapping.Chemical Science11, 7904–7909 (2020)

  30. [30]

    A., Bentria, B., Dahame, T

    Fadla, M. A., Bentria, B., Dahame, T. & Benghia, A. First-principles investigation on the stability and material properties of all-inorganic cesium lead iodide perovskites CsPbI3 polymorphs.Physica B: Condensed Matter585, 412118 (2020)

  31. [31]

    & Navrotsky, A

    Wang, B., Novendra, N. & Navrotsky, A. Energetics, Structures, and Phase Transitions of Cubic and Orthorhombic Cesium Lead Iodide (CsPbI 3 ) Polymorphs.Journal of the American Chemical Society141, 14501–14504 (2019)

  32. [32]

    Lv, S.et al.Physical origin of hafnium-based ferroelectricity.Computational Materials Today4, 100010 (2024)

  33. [33]

    Pallikara, I., Kayastha, P., Skelton, J. M. & Whalley, L. D. The physical significance of imaginary phonon modes in crystals.Electronic Structure4, 033002 (2022)

  34. [34]

    Kamath, P. D. & Persson, K. A. The impact of spurious imaginary phonon modes on thermal properties of Metal-organic Frameworks.arXiv preprint(2026).2602.07295

  35. [35]

    A., Rignanese, G.-M

    Bakhsh, M., Trinquet, V., Gouvˆ ea, R. A., Rignanese, G.-M. & Ponc´ e, S. From Symmetry to Stability: Structural and Electronic Transformation in Cs$ 2$KInI$ 6$.arXiv preprint (2026).2602.06893. 28

  36. [36]

    & Yin, J

    Bian, T., Zhu, W., Lei, Q. & Yin, J. Machine Learning Potentials for Inorganic and Hybrid Lead Halide Perovskites: From Phase Stability to Defects and Interfaces.ACS Applied Materials & Interfacesacsami.5c24460 (2026)

  37. [37]

    J.et al.Dynamic Local Structure in Caesium Lead Iodide: Spatial Correlation and Transient Domains.Small20, 2303565 (2024)

    Baldwin, W. J.et al.Dynamic Local Structure in Caesium Lead Iodide: Spatial Correlation and Transient Domains.Small20, 2303565 (2024)

  38. [38]

    Braeckevelt, T.et al.Accurately Determining the Phase Transition Temperature of CsPbI 3 via Random-Phase Approximation Calculations and Phase-Transferable Machine Learning Potentials.Chemistry of Materials34, 8561–8576 (2022)

  39. [39]

    & Sleight, A

    Shannon, R., Bierlein, J., Gillson, J., Jones, G. & Sleight, A. Polymorphism in Bi2Sn2O7. Journal of Physics and Chemistry of Solids41, 117–122 (1980)

  40. [40]

    W.et al.An Exhaustive Symmetry Approach to Structure Determination: Phase Transitions in Bi 2 Sn2 O7.Journal of the American Chemical Society138, 8031–8042 (2016)

    Lewis, J. W.et al.An Exhaustive Symmetry Approach to Structure Determination: Phase Transitions in Bi 2 Sn2 O7.Journal of the American Chemical Society138, 8031–8042 (2016)

  41. [41]

    Z.et al.A Comprehensive Assessment and Benchmark Study of Large Atomistic Foundation Models for Phonons.Advanced Intelligent Discoverye202500075 (2025)

    Anam, M. Z.et al.A Comprehensive Assessment and Benchmark Study of Large Atomistic Foundation Models for Phonons.Advanced Intelligent Discoverye202500075 (2025)

  42. [42]

    A., Vandergheynst, P

    Babu, A., Gouvˆ ea, R. A., Vandergheynst, P. & Rignanese, G.-M. MEIDNet: Multimodal generative AI framework for inverse materials design.arXiv preprint(2026).2601.22009

  43. [43]

    Shannon, R. D. Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides.Acta Crystallographica Section A32, 751–767 (1976)

  44. [44]

    npj Computational Materials11, 9 (2025)

    Deng, B.et al.Systematic softening in universal machine learning interatomic potentials. npj Computational Materials11, 9 (2025)

  45. [45]

    https://github.com/ACEsuit/mace-foundations

    GitHub - ACEsuit/mace-foundations: MACE foundation models (MP, OMAT, Matpes). https://github.com/ACEsuit/mace-foundations

  46. [46]

    Barroso-Luque, L.et al.Open Materials 2024 (OMat24) Inorganic Materials Dataset and Models.arXiv preprint(2024).2410.12771

  47. [47]

    Hjorth Larsen, A.et al.The atomic simulation environment - A Python library for working with atoms.Journal of Physics Condensed Matter29, 273002 (2017)

  48. [48]

    PHON: A program to calculate phonons using the small displacement method

    Alf` e, D. PHON: A program to calculate phonons using the small displacement method. Computer Physics Communications180, 2622–2633 (2009)

  49. [49]

    First-principles Phonon Calculations with Phonopy and Phono3py.Journal of the Physical Society of Japan92, 012001 (2023)

    Togo, A. First-principles Phonon Calculations with Phonopy and Phono3py.Journal of the Physical Society of Japan92, 012001 (2023)

  50. [50]

    Oganov, A. R. & Glass, C. W. Crystal structure prediction using ab initio evolutionary techniques: Principles and applications.arXiv preprint(2009).0911.3186

  51. [51]

    Abraham, N. L. & Probert, M. I. J. Improved real-space genetic algorithm for crystal structure and polymorph prediction.Physical Review B77, 134117 (2008)

  52. [52]

    Curtis, F.et al.GAtor: A First-Principles Genetic Algorithm for Molecular Crystal Structure Prediction.Journal of Chemical Theory and Computation14, 2246–2264 (2018)

  53. [53]

    S., Wei, L., Hu, M

    Omee, S. S., Wei, L., Hu, M. & Hu, J. Crystal structure prediction using neural network potential and age-fitness Pareto genetic algorithm.Journal of materials Informatics4 (2024). 29

  54. [54]

    Piechota, J.et al.Ab Initio Molecular Dynamics Insight to Structural Phase Transition and Thermal Decomposition of InN.International Journal of Molecular Sciences25, 8281 (2024)

  55. [55]

    & Ozgen, S

    Kazanc, S., Ahmet Celik, F. & Ozgen, S. The investigation of solid–solid phase trans- formation at CuAlNi alloy using molecular dynamics simulation.Journal of Physics and Chemistry of Solids74, 1836–1841 (2013)

  56. [56]

    P.et al.Python Materials Genomics (pymatgen): A robust, open-source python library for materials analysis.Computational Materials Science68, 314–319 (2013)

    Ong, S. P.et al.Python Materials Genomics (pymatgen): A robust, open-source python library for materials analysis.Computational Materials Science68, 314–319 (2013)

  57. [57]

    Hart, G. L. W. & Forcade, R. W. Algorithm for generating derivative structures.Physical Review B77, 224115 (2008)

  58. [58]

    S., Hart, G

    Morgan, W. S., Hart, G. L. & Forcade, R. W. Generating derivative superstructures for systems with high configurational freedom.Computational Materials Science136, 144–149 (2017)

  59. [59]

    Ramlaoui, A.et al.LeMat-Traj: A Scalable and Unified Dataset of Materials Trajectories for Atomistic Modeling.arXiv preprint(2025).2508.20875

  60. [60]

    LeMat-BulkUnique dataset (2023)

    LeMaterial. LeMat-BulkUnique dataset (2023)

  61. [61]

    & Furthm¨ uller, J

    Kresse, G. & Furthm¨ uller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set.Physical Review B - Condensed Matter and Materials Physics 54, 11169–11186 (1996)

  62. [62]

    P., Burke, K

    Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys Rev Lett77, 3865 (1996)

  63. [63]

    & Joubert, D

    Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method.Physical Review B59, 1758–1775 (1999)

  64. [64]

    [object Object] (2024)

    Ganose, A.et al.Materialsproject/atomate2: V0.0.13. [object Object] (2024)

  65. [65]

    K.et al.Accelerated data-driven materials science with the Materials Project

    Horton, M. K.et al.Accelerated data-driven materials science with the Materials Project. Nature Materials24, 1522–1532 (2025)

  66. [66]

    population

    Kuhn, H. W. The Hungarian method for the assignment problem.Naval Research Logistics Quarterly2, 83–97 (1955). Acknowledgements Computational resources have been provided by the Consortium des ´Equipements de Calcul Intensif (C´ECI), funded by the Fonds de la Recherche Scientifique de Belgique (F.R.S.-FNRS) under Grant No. 2.5020.11 and by the Walloon Reg...