Sampling-accelerated prediction of phonon scattering rates for converged thermal conductivity and radiative properties

Ziqi Guo,Zherui Han,Dudong Feng,Guang Lin,Xiulin Ruan
DOI: https://doi.org/10.1038/s41524-024-01215-8
IF: 12.256
2024-02-08
npj Computational Materials
Abstract:The prediction of thermal conductivity and radiative properties is crucial. However, computing phonon scattering, especially for four-phonon scattering, could be prohibitively expensive, and the thermal conductivity for silicon after considering four-phonon scattering is significantly under-predicted and not converged in the literature. Here we propose a method to estimate scattering rates from a small sample of scattering processes using maximum likelihood estimation. The calculation of scattering rates and associated thermal conductivity and radiative properties are dramatically accelerated by three to four orders of magnitude. This allows us to use an unprecedented q -mesh (discretized grid in the reciprocal space) of 32 × 32 × 32 for calculating four-phonon scattering of silicon and achieve a converged thermal conductivity value that agrees much better with experiments. The accuracy and efficiency of our approach make it ideal for the high-throughput screening of materials for thermal and optical applications.
materials science, multidisciplinary,chemistry, physical
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the computational challenges faced in predicting thermal conductivity and radiative properties in materials science, especially for materials such as silicon. Specifically: 1. **High computational cost of four - phonon scattering**: The calculation of four - phonon scattering is very expensive, especially when a large number of scattering processes need to be considered. This has led to the significantly underestimated and non - convergent predicted values of silicon's thermal conductivity after considering four - phonon scattering in the existing literature. 2. **Accurate prediction of thermal conductivity and radiative properties**: Thermal conductivity (\(\kappa\)) and radiative properties (such as the dielectric function \(\varepsilon\)) are important properties of materials and are crucial for thermal management and thermoelectric applications. However, accurately predicting these properties requires considering three - phonon and four - phonon scattering mechanisms, which are very time - consuming computationally. To solve these problems, the authors propose a method based on sampling and maximum likelihood estimation (MLE) to significantly reduce the computational cost and accelerate the prediction of thermal conductivity and radiative properties. The specific methods are as follows: - **Sampling method**: By randomly selecting a small part from all possible scattering processes to estimate the scattering rates (\(\tau_\lambda^{3ph}\) and \(\tau_\lambda^{4ph}\)), the amount of calculation is greatly reduced. - **Maximum likelihood estimation**: Use maximum likelihood estimation to estimate the scattering rates from the sampled data and evaluate the uncertainty of the estimate through the confidence interval. - **Efficient calculation**: This method enables the q - grid for calculating four - phonon scattering to reach an unprecedented 32 × 32 × 32, thereby obtaining convergent thermal conductivity values that are closer to the experimental results. Through this method, the authors achieve a computational acceleration of three to four orders of magnitude, making it possible to perform high - throughput screening of materials for thermal and optical applications.