High-throughput computational framework for lattice dynamics and thermal transport including high-order anharmonicity: an application to cubic and tetragonal inorganic compounds
arXiv:2507.11750 · doi:10.1038/s41524-025-01920-y
Abstract
Accurately predicting lattice thermal conductivity (kL) from first principles remains a challenge in identifying materials with extreme thermal behavior. While modern lattice dynamics methods enable routine predictions of kL within the harmonic approximation and three-phonon scattering framework (HA+3ph), reliable results, especially for low-kL compounds, require higher-order anharmonic effects, including self-consistent phonon renormalization, four-phonon scattering, and off-diagonal heat flux (SCPH+3,4ph+OD). We present a high-throughput workflow integrating these effects into a unified framework. Using this, we compute kL for 773 cubic and tetragonal inorganic compounds across diverse chemistries and structures. From 562 dynamically stable compounds, we assess the hierarchical effects of higher-order anharmonicity. For about 60% of materials, HA+3ph predictions closely match those from SCPH+3,4ph+OD. However, SCPH corrections often increase kL, sometimes by over 8 times, while four-phonon scattering universally reduces it, occasionally to 15% of the HA+3ph value. Off-diagonal contributions are minor in high-kL systems but can be comparable to the diagonal ones in highly anharmonic, low-kL compounds. We highlight four cases-Rb2TlAlH6, Cu3VS4, CuBr, and KTlCl4-exhibiting distinct anharmonic behaviors. This work delivers not only a robust workflow for high-fidelity kL dataset but also a quantitative framework to determine when higher-order effects are essential. The hierarchy of kL results, from the HA+3ph to SCPH+3,4ph+OD level, offers a scalable, interpretable route to discovering next-generation extreme thermal materials.
40 pages, 5 figures
References in corpus (17)
- Restoring the density-gradient expansion for exchange in solids and surfaces
- Generalized gradient approximation for solids and their surfaces
- AFLOW: An automatic framework for high-throughput materials discovery
- Implementation strategies in phonopy and phono3py
- A Universal Graph Deep Learning Interatomic Potential for the Periodic Table
- Unified theory of thermal transport in crystals and disordered solids
- Self-consistent phonon calculations of lattice dynamical properties in cubic SrTiO with first-principles anharmonic force constants
- Anharmonic free energies and phonon dispersions from the stochastic self-consistent harmonic approximation: application to platinum and palladium hydrides
- Quantum mechanical prediction of four-phonon scattering rates and reduced thermal conductivity of solids
- Lattice anharmonicity and thermal conductivity from compressive sensing of first-principles calculations
- The Stochastic Self-Consistent Harmonic Approximation: Calculating Vibrational Properties of Materials with Full Quantum and Anharmonic Effects
- Four-phonon scattering diminishes the optical phonon contribution and isotope effect to thermal conductivity of III-V semiconductors
- Experimental Observation of High Intrinsic Thermal Conductivity of AlN
- Compressive sensing lattice dynamics. II. Efficient phonon calculations and long-range interactions
- Compressive sensing lattice dynamics. I. General formalism
- Anharmonic stabilization and lattice heat transport in rocksalt -GeTe
- A Unified Understanding of Minimum Lattice Thermal Conductivity