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
Lattice thermal conductivity is important for many applications, but experimental measurements or first principles calculations including three-phonon and four-phonon scattering are expensive or even unaffordable. Machine learning approaches that can achieve similar accuracy have been a long-standing open question. Despite recent progress, machine learning models using structural information as descriptors fall short of experimental or first principles accuracy. This study presents a machine learning approach that predicts phonon scattering rates and thermal conductivity with experimental and first principles accuracy. The success of our approach is enabled by mitigating computational challenges associated with the high skewness of phonon scattering rates and their complex contributions to the total thermal resistance. Transfer learning between different orders of phonon scattering can further improve the model performance. Our surrogates offer up to two orders of magnitude acceleration compared to first principles calculations and would enable large-scale thermal transport informatics.
more » « less- NSF-PAR ID:
- 10418484
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- npj Computational Materials
- Volume:
- 9
- Issue:
- 1
- ISSN:
- 2057-3960
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
Abstract 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. -
Abstract The thermal conductivity of boron arsenide (BAs) is believed to be influenced by phonon scattering selection rules due to its special phonon dispersion. Compression of BAs leads to significant changes in phonon dispersion, which allows for a test of first principles theories for how phonon dispersion affects three‐ and four‐phonon scattering rates. This study reports the thermal conductivity of BAs from 0 to 30 GPa. Thermal conductivity vs. pressure of BAs is measured by time‐domain thermoreflectance with a diamond anvil cell. In stark contrast to what is typical for nonmetallic crystals, BAs is observed to have a pressure independent thermal conductivity below 30 GPa. The thermal conductivity of nonmetallic crystals typically increases upon compression. The unusual pressure independence of BAs's thermal conductivity shows the important relationship between phonon dispersion properties and three‐ and four‐phonon scattering rates.
-
Abstract Although first principles based anharmonic lattice dynamics is one of the most common methods to obtain phonon properties, such method is impractical for high-throughput search of target thermal materials. We develop an elemental spatial density neural network force field as a bottom-up approach to accurately predict atomic forces of ~80,000 cubic crystals spanning 63 elements. The primary advantage of our indirect machine learning model is the accessibility of phonon transport physics at the same level as first principles, allowing simultaneous prediction of comprehensive phonon properties from a single model. Training on 3182 first principles data and screening 77,091 unexplored structures, we identify 13,461 dynamically stable cubic structures with ultralow lattice thermal conductivity below 1 Wm −1 K −1 , among which 36 structures are validated by first principles calculations. We propose mean square displacement and bonding-antibonding as two low-cost descriptors to ease the demand of expensive first principles calculations for fast screening ultralow thermal conductivity. Our model also quantitatively reveals the correlation between off-diagonal coherence and diagonal populations and identifies the distinct crossover from particle-like to wave-like heat conduction. Our algorithm is promising for accelerating discovery of novel phononic crystals for emerging applications, such as thermoelectrics, superconductivity, and topological phonons for quantum information technology.more » « less
-
Abstract Nanosized perovskite ferroelectrics are widely employed in several electromechanical, photonics, and thermoelectric applications. Scaling of ferroelectric materials entails a severe reduction in the lattice (phonon) thermal conductivity, particularly at sub‐100 nm length scales. Such thermal conductivity reduction can be accurately predicted using the information of phonon mean free path (MFP) distribution. The current understanding of phonon MFP distribution in perovskite ferroelectrics is still inconclusive despite the critical thermal management implications. Here, high‐quality single‐crystalline barium titanate (BTO) thin films, a representative perovskite ferroelectric material, are grown at several thicknesses. Using experimental thermal conductivity measurements and first‐principles based modeling (including four‐phonon scattering), the phonon MFP distribution is determined in BTO. The simulation results agree with the measured thickness‐dependent thermal conductivity. The results show that the phonons with sub‐100 nm MFP dominate the thermal transport in BTO, and phonons with MFP exceeding 10 nm contribute ≈35% to the total thermal conductivity, in significant contrast to previously published experimental results. The experimentally validated phonon MFP distribution is consistent with the theoretical predictions of other complex crystals with strong anharmonicity. This work paves the way for thermal management in nanostructured and ferroelectric‐domain‐engineered systems for oxide perovskite‐based functional materials.
-
Abstract Hexagonal boron nitride (h-BN) has been predicted to exhibit an in-plane thermal conductivity as high as ~ 550 W m−1K−1at room temperature, making it a promising thermal management material. However, current experimental results (220–420 W m−1K−1) have been well below the prediction. Here, we report on the modulation of h-BN thermal conductivity by controlling the B isotope concentration. For monoisotopic10B h-BN, an in-plane thermal conductivity as high as 585 W m−1K−1is measured at room temperature, ~ 80% higher than that of h-BN with a disordered isotope concentration (52%:48% mixture of10B and11B). The temperature-dependent thermal conductivities of monoisotopic h-BN agree well with first principles calculations including only intrinsic phonon-phonon scattering. Our results illustrate the potential to achieve high thermal conductivity in h-BN and control its thermal conductivity, opening avenues for the wide application of h-BN as a next-generation thin-film material for thermal management, metamaterials and metadevices.