🔬 materials science

Comparative Assessment of Thermal Transport Theories: Dual-Channel Mechanism Dictates Heat Transport in Ultralow- $\kappa$ Materials

This study demonstrates that the Wigner transport equation (WTE) framework, which accounts for both particle-like and wave-like coherence channels alongside higher-order phonon scattering, successfully predicts the ultralow thermal conductivity of strongly anharmonic materials like TlAgSe and Cs $_2$ PbI $_2$ C $_2$ , resolving the limitations of conventional Green-Kubo and Boltzmann transport theories.

Original authors: Soham Mandal, Ashutosh Srivastava, Tanmoy Das, Manish Jain, Abhishek Kumar Singh, Prabal K. Maiti

Published 2026-08-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Soham Mandal, Ashutosh Srivastava, Tanmoy Das, Manish Jain, Abhishek Kumar Singh, Prabal K. Maiti

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Heat is a traveler that moves through solids in ways we are only just beginning to fully understand. In most solid materials, such as the silicon in a computer chip or the ceramic in a coffee mug, heat travels as a stream of tiny, vibrating waves called phonons. These waves behave somewhat like a gas of particles, bouncing off one another and carrying energy from hot spots to cold ones. Scientists have long used a standard set of rules to predict how fast this heat moves, a model that works very well for ordinary materials. However, there is a special class of crystals where the atoms do not vibrate in neat, predictable patterns. Instead, they jiggle chaotically, creating a situation where the standard rules break down. In these materials, heat moves incredibly slowly, making them perfect candidates for turning waste heat into electricity or for protecting engines from extreme temperatures. The challenge has been that our best mathematical tools for describing heat either fail to capture the chaotic nature of these atoms or require so much computing power that they become impractical to use.

A team of researchers at the Indian Institute of Science has now taken a fresh look at this problem by studying two specific materials: a metal compound called TlAgSe and a layered crystal known as Cs2PbI2Cl2. Both of these substances are famous for having exceptionally low ability to conduct heat, a property that makes them valuable for energy applications. The researchers wanted to understand exactly how heat moves through these chaotic structures. They found that the old way of thinking, which treats heat carriers as simple particles, was missing a crucial piece of the puzzle. While that traditional model could predict some of the heat flow, it consistently underestimated how much heat actually moved. The team discovered that heat in these materials does not just travel as a stream of particles; it also moves through a second, hidden channel where the vibrations of different atoms stay in sync with one another, acting more like a wave than a particle.

To uncover this hidden mechanism, the scientists used a powerful combination of modern computing techniques. They first simulated the movement of atoms at different temperatures using a type of artificial intelligence that mimics the behavior of real atoms with high accuracy. This allowed them to see how the crystal structures distorted as they got hotter, a detail that older, simpler models often missed. They then applied two different theoretical frameworks to interpret these simulations. The first was the traditional approach, which counts the number of heat-carrying waves and how often they collide. The second was a more advanced framework that also accounts for the wave-like connections between different vibrations. When they compared these calculations to real-world measurements, the traditional approach fell short, predicting a thermal conductivity of about 0.22 watts per meter-kelvin for TlAgSe, while experiments showed the actual value was around 0.30. The new, more complete model, which included the wave-like connections, predicted a value of 0.38, which was much closer to reality, though still slightly high.

The researchers realized that the slight overestimation came from ignoring a specific type of interaction where four vibrations collide at once, rather than just the usual three. When they added these four-way collisions into their advanced model, the predictions became remarkably accurate. For TlAgSe, the final calculated value dropped to 0.31, and for Cs2PbI2Cl2, it settled at 0.38, both matching experimental data almost perfectly. This result confirmed that heat transport in these materials is a dual-channel process. One channel involves the standard particle-like scattering, which becomes less effective as the material gets hotter. The other channel involves a wave-like coherence, where vibrations in different parts of the crystal move together in a coordinated fashion. This second channel actually becomes stronger as the temperature rises, a behavior that resembles how heat moves through glass rather than a crystal.

The study also revealed that in these materials, the vibrations are so chaotic that many of them lose their identity as distinct particles entirely. In the standard model, a vibration is treated as a clear, individual wave that travels until it hits something. In these ultralow-conductivity crystals, the vibrations are so heavily damped by the chaotic atomic motion that they blur together. The researchers found that for a large number of these vibrations, the chaos is so intense that the standard particle picture no longer applies at all. Instead, the heat moves through a diffusive process driven by the wave-like connections between these blurred vibrations. This finding explains why the old models failed: they tried to count individual particles in a system where the particles had effectively merged into a collective wave.

By combining artificial intelligence simulations with these advanced theoretical tools, the team provided a unified description of how heat moves in some of the most complex crystals known. Their work shows that to truly understand and design materials for thermoelectric devices or thermal barriers, scientists must look beyond the simple idea of bouncing particles. They must account for the subtle, wave-like cooperation between atomic vibrations and the intense, multi-way collisions that occur at high temperatures. This deeper understanding offers a solid foundation for engineers who want to create better materials for capturing energy or insulating against heat, turning a fundamental mystery of physics into a practical guide for future technology.

Technical Summary: Comparative Assessment of Thermal Transport Theories in Ultralow-κ Materials

Problem Statement
Theoretical understanding of heat transport in strongly anharmonic crystalline solids with ultralow lattice thermal conductivity ( $\kappa$ ) faces significant challenges. Conventional frameworks often fail to capture experimental observations in these materials. Specifically, the standard Peierls-Boltzmann transport equation (BTE), based on the phonon quasiparticle and phonon-gas models, proves inadequate for materials exhibiting strong anharmonicity, such as metal chalcogenides and halide perovskites. The BTE typically underestimates $\kappa$ because it neglects wave-like interbranch coherence and fails when phonons enter an overdamped regime. Conversely, while the Green-Kubo (GK) framework accurately reproduces experimental $\kappa$ at high temperatures by capturing nonperturbative anharmonic interactions, it offers limited microscopic insight, neglects Bose-Einstein (BE) statistics of lattice vibrations, and relies on classical molecular dynamics (MD) that misses nuclear quantum effects at lower temperatures. There is a critical need for a unified theoretical framework that accounts for arbitrary anharmonic strength, full temperature-dependent lattice dynamics, and both particle-like and wave-like transport mechanisms.

Methodology
This study investigates the thermal transport properties of two representative ultralow- $\kappa$ materials: TlAgSe (a metal chalcogenide) and Cs $_2$ PbI $_2$ Cl $_2$ (an all-inorganic layered Ruddlesden-Popper perovskite). The authors employ a multi-faceted computational approach:

Machine Learning Interatomic Potentials (MLIP): To bypass the prohibitive computational cost of ab initio molecular dynamics (AIMD) for these complex systems, MLIPs trained on density functional theory (DFT) data were used. These potentials enabled the simulation of finite-temperature lattice dynamics and the extraction of interatomic force constants (IFCs) up to the fourth order, capturing temperature-induced structural distortions and local symmetry breaking.
Perturbative Transport Frameworks: The study utilizes the Wigner transport equation (WTE) framework, which extends the conventional BTE. The WTE explicitly separates thermal conductivity into two channels:
- Populations ( $\kappa_P$ ): The diagonal term corresponding to particle-like phonon quasiparticle propagation (equivalent to BTE).
- Coherences ( $\kappa_C$ ): The off-diagonal term representing wave-like energy transfer between quasi-degenerate eigenstates (interbranch coherence).
  Calculations included contributions from three-phonon and four-phonon scattering processes, as well as phonon-isotope impurity scattering.
Green-Kubo (GK) Benchmarking: Results from the perturbative WTE and BTE were compared against nonperturbative GK predictions derived from equilibrium MD simulations, which serve as a benchmark for capturing full anharmonic interactions.

Key Results

Dual-Channel Mechanism: The study reveals that heat conduction in these materials is governed by a dual-channel mechanism. While the BTE-based $\kappa_P$ $κ_{P}$ underestimates experimental values, the inclusion of the coherence channel ( $\kappa_C$ $κ_{C}$ ) within the WTE framework significantly increases the predicted $\kappa$ $κ$ .
- For TlAgSe, $\kappa_C$ contributes approximately 42% to the total $\kappa$ at 300 K, rising to 54% at 500 K. The total WTE prediction ( $\sim$ 0.38 Wm $^{-1}$ K $^{-1}$ ) aligns well with experiments ( $\sim$ 0.30 Wm $^{-1}$ K $^{-1}$ ) only when four-phonon scattering is included.
- For Cs $_2$ PbI $_2$ Cl $_2$ , the coherence contribution is smaller (27% at 300 K) but still essential. The WTE prediction ( $\sim$ 0.38 Wm $^{-1}$ K $^{-1}$ ) matches experimental ranges ( $\sim$ 0.37–0.41 Wm $^{-1}$ K $^{-1}$ ) when combined with four-phonon scattering.
Necessity of Four-Phonon Scattering: Restricting perturbative calculations to three-phonon interactions leads to a systematic overestimation of $\kappa$ when coherence is included, or underestimation when only populations are considered. The inclusion of quartic anharmonicity (four-phonon scattering) is essential to correct the phonon scattering rates and achieve quantitative agreement with experimental data.
Temperature Dependence: The $\kappa_P$ component decreases with increasing temperature, whereas $\kappa_C$ exhibits a "glass-like" rise. This behavior results in a total $\kappa$ with very weak temperature dependence, consistent with experimental observations.
Breakdown of Quasiparticle Picture: Analysis of phonon scattering rates ( $\tau^{-1}$ ) shows that a substantial fraction of phonon modes in both materials exceed the Ioffe-Regel limit ( $\tau^{-1} > \nu$ ), indicating that phonons become overdamped. Furthermore, many modes exceed the Wigner limit ( $\tau^{-1} > \Delta\omega_{ave.}$ ), confirming the dominance of wave-like interbranch coherence over particle-like propagation in these regimes.

Significance and Claims
The paper claims to provide a unified description of heat transport in ultralow- $\kappa$ materials by reconciling the strengths of different theoretical frameworks. The authors demonstrate that:

Coherence is Critical: The wave-like interbranch coherence ( $\kappa_C$ ) is not a minor correction but a dominant transport channel in strongly anharmonic crystals, challenging the validity of the conventional phonon-gas model based solely on the BTE.
Higher-Order Anharmonicity is Essential: Accurate modeling of these systems requires the inclusion of four-phonon scattering processes to correctly capture scattering rates and avoid discrepancies with experimental data.
Methodological Validation: By combining MLIP-based finite-temperature dynamics with the WTE framework, the study offers a computationally feasible yet rigorous approach to extract microscopic insights (such as scattering rates and coherence contributions) that are inaccessible via nonperturbative GK methods alone, while retaining the accuracy of GK predictions.

This computational approach establishes a basis for the rational design of phononic and thermoelectric devices by accurately predicting thermal transport in materials where conventional theories fail.

More like this