Hierarchy of Exchange-Correlation Functionals in Computing Lattice Thermal Conductivities of Rocksalt and Zincblende Semiconductors

2024/04/09 by Jiacheng Wei, Wei, Jiacheng, Zhonghao Xia +5 · 5 citations
Materials Science · #Advanced Thermoelectric Materials and Devices #FOS: Physical sciences #Materials Science (cond-mat.mtrl-sci) #Thermal Expansion and Ionic Conductivity #Thermal properties of materials

paper · pdf · doi:10.48550/arxiv.2404.06346

openalex publication_date 2024/04/09 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28

Abstract

Lattice thermal conductivity (κ_{\rm L}) is a crucial characteristic of crystalline solids with significant implications for practical applications. While the higher order of anharmonicity of phonon gas model is commonly used for explaining extraordinary heat transfer behaviors in crystals, the impact of exchange-correlation (XC) functionals in DFT on describing anharmonicity has been largely overlooked. Most XC functionals in solids focus on ground state properties that mainly involve the harmonic approximation, neglecting temperature effects, and their reliability in studying anharmonic properties remains insufficiently explored. In this study, we systematically investigate the room-temperature κ_{\rm L} of 16 binary compounds with rocksalt and zincblende structures using 8 XC functionals such as LDA, PBE, PBEsol, optB86b, revTPSS, SCAN, rSCAN, r²SCAN in combination with three perturbation orders, including phonon within harmonic approximation (HA) plus three-phonon scattering (HA+3ph), phonon calculated using self-consistent phonon theory (SCPH) plus three-phonon scattering (SCPH+3ph), and SCPH phonon plus three- and four-phonon scattering (SCPH+3,4ph). Our results show that the XC functional exhibits strong entanglement with perturbation order and the mean relative absolute error (MRAE) of the computed κ_{\rm L} is strongly influenced by both the XC functional and perturbation order, leading to error cancellation or amplification. The minimal (maximal) MRAE is achieved with revTPSS (rSCAN) at the HA+3ph level, SCAN (r²SCAN) at the SCPH+3ph level, and PBEsol (rSCAN) at the SCPH+3,4ph level. Among these functionals, PBEsol exhibits the highest accuracy at the highest perturbation order. The SCAN-related functionals demonstrate moderate accuracy but are suffer from numerical instability and high computational costs.

Cited by

• A Review of AI-Driven Approaches for Nanoscale Heat Conduction and Radiation
• Optimizing Cross-Domain Transfer for Universal Machine Learning Interatomic Potentials
• FourPhonon_GPU: A GPU-accelerated framework for calculating phonon scattering rates and thermal conductivity
• Revisiting thermal transport in CuCl: First-principles calculations and machine learning force fields
• Thermoelectric Properties of Copper-based Chalcopyrite Semiconductors CuMX₂ (M = Al, Ga, and In; X = S, Se, and Te) from First-Principles Calculations

Related

• High-throughput computational framework for lattice dynamics and thermal transport including high-order anharmonicity: an application to cubic and tetragonal inorganic compounds
• Phonon properties and thermal conductivity from first principles, lattice dynamics, and the Boltzmann transport equation
• Ab initio Investigation of Thermal Transport in Insulators: Unveiling the Roles of Phonon Renormalization and Higher-Order Anharmonicity
• Lattice Thermal Transport Beyond the Quasiparticle Approximation: Nontrivial Spectral Competition between Three- and Four-Phonon Interactions
• Quantitative assessment of perturbation theory-based lattice thermal conductivity models using quasi-continuum approximation