skip to main content


This content will become publicly available on December 1, 2024

Title: Unlocking phonon properties of a large and diverse set of cubic crystals by indirect bottom-up machine learning approach
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
Award ID(s):
2030128 2110033
NSF-PAR ID:
10465913
Author(s) / Creator(s):
; ; ; ; ; ; ;
Date Published:
Journal Name:
Communications Materials
Volume:
4
Issue:
1
ISSN:
2662-4443
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Existing machine learning potentials for predicting phonon properties of crystals are typically limited on a material-to-material basis, primarily due to the exponential scaling of model complexity with the number of atomic species. We address this bottleneck with the developed Elemental Spatial Density Neural Network Force Field, namely Elemental-SDNNFF. The effectiveness and precision of our Elemental-SDNNFF approach are demonstrated on 11,866 full, half, and quaternary Heusler structures spanning 55 elements in the periodic table by prediction of complete phonon properties. Self-improvement schemes including active learning and data augmentation techniques provide an abundant 9.4 million atomic data for training. Deep insight into predicted ultralow lattice thermal conductivity (<1 Wm −1  K −1 ) of 774 Heusler structures is gained by p–d orbital hybridization analysis. Additionally, a class of two-band charge-2 Weyl points, referred to as “double Weyl points”, are found in 68% and 87% of 1662 half and 1550 quaternary Heuslers, respectively. 
    more » « less
  2. Full Heusler compounds have long been discovered as exceptional n-type thermoelectric materials. However, no p-type compounds could match the high n-type figure of merit ( ZT ). In this work, based on first-principles transport theory, we predict the unprecedentedly high p-type ZT = 2.2 at 300 K and 5.3 at 800 K in full Heusler CsK 2 Bi and CsK 2 Sb, respectively. By incorporating the higher-order phonon scattering, we find that the high ZT value primarily stems from the ultralow lattice thermal conductivity ( κ L ) of less than 0.2 W mK −1 at room temperature, decreased by 40% compared to the calculation only considering three-phonon scattering. Such ultralow κ L is rooted in the enhanced phonon anharmonicity and scattering channels stemming from the coexistence of antibonding-induced anharmonic rattling of Cs atoms and low-lying optical branches. Moreover, the flat and heavy nature of valence band edges leads to a high Seebeck coefficient and moderate power factor at optimal hole concentration, while the dispersive and light conduction band edges yield much larger electrical conductivity and electronic thermal conductivity ( κ e ), and the predominant role of κ e suppresses the n-type ZT . This study offers a deeper insight into the thermal and electronic transport properties in full Heusler compounds with strong phonon anharmonicity and excellent thermoelectric performance. 
    more » « less
  3. null (Ed.)
    Thermal anisotropy/isotropy is one of the fundamental thermal transport properties of materials and plays a critical role in a wide range of practical applications. Manipulation of anisotropic to isotropic thermal transport or vice versa is in increasing demand. However, almost all the existing approaches for tuning anisotropy or isotropy focus on structure engineering or materials processing, which is time and cost consuming and irreversible, while little progress has been made with an intact, robust, and reversible method. Motivated by the inherent relationship between interatomic interaction mediated phonon transport and electronic charges, we comprehensively investigate the effect of external electric field on thermal transport in two-dimensional (2D) borophene by performing first-principles calculations along with the phonon Boltzmann transport equation. Under external electric field, the lattice thermal conductivity of borophene in both in-plane directions first increases significantly to peak values with the maximum augmentation factor of 2.82, and the intrinsic anisotropy (the ratio of thermal conductivity along two in-plane directions) is boosted to the highest value of 2.13. After that, thermal conductivities drop down steeply and anisotropy exhibits oscillating decay. With the electric field increasing to 0.4 V Å −1 , the thermal conductivity is dramatically suppressed to 1/40 of the original value at no electric field. More interestingly, the anisotropy of the thermal conductivity decreases to the minimum value of 1.25, showing almost isotropic thermal transport. Such abnormal anisotropic to isotropic thermal transport transition stems from the large enhancement and suppression of phonon lifetime at moderate and high strength of electric field, respectively, and acts as an amplifying or reducing factor to the thermal conductivity. We further explain the tunability of phonon lifetime of the dominant acoustic mode by an electron localization function. By comparing the electric field-modulated thermal conductivity of borophene with the dielectric constant, it is found that the screened potential resulting from the redistributed charge density leads to phonon renormalization and the modulation of phonon anharmonicity and anisotropy through electric field. Our study paves the way for robust tuning of anisotropy of phonon transport in materials by applying intact, robust, and reversible external electric field without altering their atomic structure and would have a significant impact on emerging applications, such as thermal management of nanoelectronics and thermoelectric energy conversion. 
    more » « less
  4. The Mg 3 Sb 2− x Bi x family has emerged as the potential candidates for thermoelectric applications due to their ultra-low lattice thermal conductivity ( κ L ) at room temperature (RT) and structural complexity. Here, using ab initio calculations of the electron-phonon averaged (EPA) approximation coupled with Boltzmann transport equation (BTE), we have studied electronic, phonon and thermoelectric properties of Mg 3 Sb 2− x Bi x (x = 0, 1, and 2) monolayers. In violation of common mass-trend expectations, increasing Bi element content with heavier Zintl phase compounds yields an abnormal change in κ L in two-dimensional Mg 3 Sb 2− x Bi x crystals at RT (∼0.51, 1.86, and 0.25 W/mK for Mg 3 Sb 2 , Mg 3 SbBi, and Mg 3 Bi 2 ). The κ L trend was detailedly analyzed via the phonon heat capacity, group velocity and lifetime parameters. Based on quantitative electronic band structures, the electronic bonding through the crystal orbital Hamilton population (COHP) and electron local function analysis we reveal the underlying mechanism for the semiconductor-semimetallic transition of Mg 3 Sb 2-− x Bi x compounds, and these electronic transport properties (Seebeck coefficient, electrical conductivity, and electronic thermal conductivity) were calculated. We demonstrate that the highest dimensionless figure of merit ZT of Mg 3 Sb 2− x Bi x compounds with increasing Bi content can reach ∼1.6, 0.2, and 0.6 at 700 K, respectively. Our results can indicate that replacing heavier anion element in Zintl phase Mg 3 Sb 2− x Bi x materials go beyond common expectations (a heavier atom always lead to a lower κ L from Slack’s theory), which provide a novel insight for regulating thermoelectric performance without restricting conventional heavy atomic mass approach. 
    more » « less
  5. Utilizing first-principles calculations combined with phonon Boltzmann transport theory up to fourth-order anharmonicity, we systematically investigate the thermal transport properties of the biphenylene network [BPN, recently synthesized experimentally by Fan et al. , Science , 2021, 372 , 852–856] and hydrogenated BPN (HBPN). The calculations show that four-phonon scattering significantly affects the lattice thermal conductivity ( κ ) of BPN. At room temperature, the κ of BPN is reduced from 582.32 (1257.07) W m −1 K −1 to 309.56 (539.88) W m −1 K −1 along the x ( y ) direction after considering the four-phonon scattering. Moreover, our results demonstrate that the thermal transport in BPN could also be greatly suppressed by hydrogenation, where the κ of HBPN along the x ( y ) direction is merely 16.62% (10.14%) of that of pristine BPN at 300 K. The mechanism causing such an obvious decrease of κ of HBPN is identified to be due to the enhanced phonon scattering rate and reduced group velocity, which is further revealed by the increased scattering phase space and weakened C–C bond. The results presented in this work shed light on the intrinsic thermal transport features of BPN and HBPN, which will help us to understand the phonon transport processes and pave the way for their future developments in the thermal field. 
    more » « less