Accounting for the vibrational contribution to the configurational entropy in disordered solids with machine learned forcefields: a case study of garnet electrolyte Li <sub>7</sub> La <sub>3</sub> Zr <sub>2</sub> O <sub>12</sub>
Accounting for lattice vibrations to accurately determine the phase stabilities of site-disordered solids is a long-standing challenge in computational material designs, due to the high computational cost associated with sampling the vast configurational space to obtain the converged thermodynamic quantities. One example is the garnet electrolyte Li<sub>7</sub>La<sub>3</sub>Zr<sub>2</sub>O<sub>12</sub>, the high-temperature and high-ion-mobility cubic phase of which is disordered in its Li<sup>+</sup> site occupations, such that both the vibrational and configurational entropic contributions to its phase stability cannot be ignored. Understanding the subtle interplay between vibrational and configurational entropies in this material will therefore play a critical role in the rational manipulation of dopants and defects to stabilise cubic Li<sub>7</sub>La<sub>3</sub>Zr<sub>2</sub>O<sub>12</sub> at room temperature for practical applications. Here, by developing machine learned forcefields based on an equivariant message-passing neural network SO3KRATES, we follow a strict statistical thermodynamic protocol to quantify the phase stability of cubic Li<sub>7</sub>La<sub>3</sub>Zr<sub>2</sub>O<sub>12</sub> through structural optimisations, as well as molecular dynamics simulations at 300 and 1500 K, for a total of 70 120 configurations of cubic Li<sub>7</sub>La<sub>3</sub>Zr<sub>2</sub>O<sub>12</sub>. Although this only covers a tiny fraction of the configurational space (∼7 × 10<sup>34</sup> configurations in total), we are able to deterministically show that the vibrational contributions to the total configurational free energy at 1500 K are significant (on the order of 1 eV per atom) in correctly ordering the stability of the cubic Li<sub>7</sub>La<sub>3</sub>Zr<sub>2</sub>O<sub>12</sub> over its tetragonal counterpart, thanks to the high data efficiency, accuracy, stability and good transferability of the transformer-based equivariant network architecture behind SO3KRATES. Therefore, our work opens up new avenues to accelerate the accurate computational designs of disordered solids, such as solid electrolytes, for technologically important applications.
Authors
- Jack Y. Yang (ORCID: https://orcid.org/0000-0003-3589-6807)
- Sean Li (ORCID: https://orcid.org/0000-0003-4437-8817)
- Ziqi Yin (ORCID: https://orcid.org/0000-0001-8708-358X)
Institutions
- UNSW Sydney (AU)
Publication Details
- Journal
- Physical Chemistry Chemical Physics
- Published
- 2025-01-01
- DOI
- https://doi.org/10.1039/d5cp00138b
- Citations
- 1
- Primary Topic
- Solid-state spectroscopy and crystallography
- Type
- article
- Field-Weighted Citation Impact
- 0.52
Funders
- National Computational Infrastructure
- University of New South Wales
- National Cancer Institute