🔬 materials science

Microscopic Origin of the Ultralow Lattice Thermal Conductivity in Vacancy-Ordered Halide Double Perovskites Cs $_2BX_6$ ( $B$ = Zr, Pd, Sn, Te, Hf, and Pt; $X$ = Cl, Br, and I)

This study employs first-principles calculations and machine learning to reveal that the ultralow lattice thermal conductivity in vacancy-ordered Cs $_2BX_6$ double perovskites stems primarily from intrinsically weak chemical bonding leading to low sound velocities, rather than the rattling phonon modes typically associated with their structural voids.

Original authors: Lingzhi Cao, Yateng Wang, Zhonghao Xia, Jiangang He

Published 2026-02-09

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Lingzhi Cao, Yateng Wang, Zhonghao Xia, Jiangang He

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

Imagine a bustling city where the "citizens" are tiny particles called atoms, and the "traffic" is heat trying to move from one side of the city to the other. In most solid materials, this heat traffic flows smoothly and quickly, like cars on a wide, straight highway. This is why metals feel hot to the touch so fast.

However, the researchers in this paper studied a special family of materials called Cs2BX6. Think of these materials as a city built with a very specific, slightly broken blueprint. They are "vacancy-ordered," which means the city planners intentionally left some empty lots (vacancies) in a perfect pattern.

Here is the simple breakdown of what they found:

1. The Mystery of the "Super-Insulator"

The scientists wanted to know why these materials are so good at stopping heat. In the world of physics, a material that stops heat well is like a thick winter coat. These specific crystals are so effective that they are "ultra-low" conductors. They are so good at blocking heat that they could be used to keep things cold or hot without wasting energy (like in thermoelectric devices or thermal insulation).

2. The Old Theory vs. The New Discovery

For a long time, scientists had a favorite theory for why these materials were so good at blocking heat. They thought it was because of "rattling."

The Old Idea (The Rattling Cage): Imagine a large, empty cage with a small ball inside. If you shake the cage, the ball rattles around wildly, bumping into the walls and stopping any smooth movement. Scientists thought the large empty spaces in these crystals acted like cages, and the atoms inside were "rattling" so much that they blocked the heat traffic.
The New Discovery (The Weak Roads): The researchers used super-powerful computer simulations (like a high-tech traffic camera) to look closely at what was actually happening. They found that the "rattling" wasn't the main culprit. Instead, the problem was the roads themselves.

The chemical bonds holding these atoms together are naturally weak. Imagine trying to drive a car on a road made of jelly instead of asphalt. The road is so soft and wobbly that the car (the heat) can't pick up speed. Because the "roads" (chemical bonds) are weak, the heat moves incredibly slowly. This is the main reason these materials are such good insulators.

3. The "Machine Learning" Detective

To prove this, the researchers used a "machine learning" detective. They fed the computer data about the atoms, the size of the empty lots, and how fast heat moved. The computer learned a simple rule: The weaker the bond, the slower the heat moves.

It turned out that the speed of heat in these materials is directly tied to how "soft" the chemical connections are. The heavier the atoms and the weaker the bonds, the slower the heat travels.

4. The One Exception: The "Traffic Jam"

There was one special case in their study: Cs2SnI6 (a compound with Tin and Iodine).
In this specific material, the "weak road" theory wasn't the only thing happening. Here, the Iodine atoms were vibrating in a way that created a massive traffic jam. It wasn't just that the roads were soft; the atoms were bumping into each other so chaotically (a phenomenon called "strong scattering") that it created the ultimate traffic stop. This made Cs2SnI6 the absolute best at blocking heat of all the materials they studied.

5. Why It Matters (According to the Paper)

The paper concludes that if we want to design new materials that are excellent at blocking heat (for things like thermal insulation or energy conversion), we shouldn't just look for big empty spaces to make atoms "rattle." Instead, we should look for materials with intrinsically weak chemical bonds.

In a nutshell:
These materials are like a city where the streets are made of jelly. Heat tries to drive through, but the streets are so soft and wobbly that the heat cars can't go fast. In one specific neighborhood (the Iodine one), there's also a massive, chaotic traffic jam that stops everything completely. This discovery helps scientists understand how to build better "thermal coats" for our technology.

Technical Summary: Microscopic Origin of Ultralow Lattice Thermal Conductivity in Vacancy-Ordered Halide Double Perovskites

Problem Statement
Lead-free halide vacancy-ordered double perovskites with the formula $Cs_2BX_6$ (where B = Zr, Pd, Sn, Te, Hf, Pt and X = Cl, Br, I) have emerged as promising candidates for optoelectronic and thermoelectric applications due to their intrinsic ultralow lattice thermal conductivity ( $\kappa_L$ ). While previous studies have identified specific compounds like $Cs_2SnI_6$ and $Cs_2SnBr_6$ as having exceptionally low $\kappa_L$ , a systematic understanding of how $\kappa_L$ evolves across the entire $Cs_2BX_6$ family remains lacking. Specifically, the microscopic mechanisms governing these ultralow values—whether they arise from strong phonon scattering due to "rattling" modes in large structural voids or from intrinsically weak chemical bonding—have not been comprehensively clarified. Furthermore, the relative contributions of three-phonon, four-phonon, and coherent phonon transport processes to the total thermal resistance in these materials require rigorous quantification.

Methodology
The authors employed state-of-the-art first-principles calculations based on Density Functional Theory (DFT) using the Vienna Ab initio Simulation Package (VASP). The study utilized the PBEsol exchange–correlation functional and the projector augmented-wave (PAW) method. Key computational steps included:

Force Constants: Harmonic force constants were calculated using the finite-displacement method (Phonopy) with large supercells (4×4×4). Higher-order anharmonic force constants were extracted using the compressive sensing lattice dynamics (CSLD) method.
Phonon Renormalization: To account for temperature effects, second-order force constants were renormalized using the self-consistent phonon (SCPH) method (ALAMODE package).
Thermal Transport: Lattice thermal conductivity was calculated by iteratively solving the Peierls–Boltzmann transport equation. The study incorporated both three-phonon and four-phonon scattering processes (FourPhonon package) and evaluated the coherent contribution ( $\kappa^C_L$ ) using a unified theory of thermal transport for crystals and glasses.
Machine Learning: A predictive model for $\kappa_L$ was constructed using the Sure Independence Screening and Sparsifying Operator (SISSO) algorithm to identify the most critical physical descriptors governing thermal transport.
Chemical Bonding Analysis: Crystal Orbital Hamilton Population (COHP) analysis was performed to quantify chemical bonding strengths and their relationship to force constants.

Key Results

Ultralow Thermal Conductivity: All investigated $Cs_2BX_6$ compounds exhibit intrinsically ultralow $\kappa_L$ values below 1.0 W m $^{-1}$ K $^{-1}$ at 300 K. The values deviate significantly from the conventional $T^{-1}$ temperature dependence, showing a much weaker temperature dependence ( $\kappa_L \propto T^{-\alpha}$ with $\alpha \ll 1$ ).
Dominant Mechanism: Contrary to the expectation that large structural voids lead to strong "rattling" phonon scattering, the analysis reveals that the primary factor suppressing heat transport in most compounds is the low phonon group velocity ( $\nu_g$ ), ranging from 1100 to 1600 m s $^{-1}$ . This low velocity stems from intrinsically weak chemical bonding within the halide framework.
Scattering Contributions:
- Three-phonon scattering is the dominant resistance mechanism at room temperature for most compounds.
- Four-phonon scattering reduces $\kappa_L$ by approximately 20–25% in most cases, though this varies with composition (e.g., up to 35% in $Cs_2ZrCl_6$ ).
- Coherent transport contributes roughly 15% to the total $\kappa_L$ in most systems, rising to 21% in $Cs_2SnI_6$ .
Specific Anomalies:
- $Cs_2SnI_6$ : Exhibits the lowest $\kappa_L$ (0.28 W m $^{-1}$ K $^{-1}$ ) among the studied compounds. Unlike the general trend, its ultralow conductivity is driven by exceptionally strong three-phonon scattering in the low-frequency range (0.5–1.0 THz) due to a large scattering phase space associated with flat optical branches dominated by Iodine vibrations.
- $Cs_2TeX_6$ : Shows anomalous behavior where $\kappa_L$ does not strictly decrease with increasing halide mass (I > Br > Cl), attributed to stronger three-phonon scattering in the chloride and bromide variants compared to the iodide.
Machine Learning Insights: The SISSO model, which outperformed traditional Slack and modified Slack models (achieving an $R^2$ of 0.96), identified that $\kappa_L$ scales linearly with sound velocity ( $\nu_g$ ) rather than $\nu_g^3$ . The model highlighted that the combination of sound velocity, bulk modulus, and the mean-square displacement of the B-site atom is crucial for prediction, reinforcing the role of weak chemical bonding.
Chemical Bonding Trends: COHP analysis confirmed that B–X bonding strength decreases from Cl to I. The electronic configuration of the B-site cation significantly influences bonding; for instance, Pd-based compounds show weaker bonding (lower integrated COHP) than Zr-based ones due to the occupation of antibonding $t_{2g}$ states.

Significance and Claims
The paper claims to provide a unified microscopic understanding of thermal transport in vacancy-ordered double perovskites. It challenges the prevailing assumption that "rattling" modes in large voids are the primary cause of ultralow $\kappa_L$ in this family. Instead, the authors assert that the intrinsically weak chemical bonding leading to low sound velocities is the dominant mechanism for the majority of $Cs_2BX_6$ compounds.

The study clarifies the distinct roles of B-site and X-site chemistry:

The halide species (X) plays a dominant role in determining $\kappa_L$ in chloride-based compounds.
The influence of the B-site cation is more pronounced in chlorides but less significant in bromides and iodides, with the notable exception of Tellurium-based compounds.

By establishing that weak bonding and low sound velocities are the root causes of ultralow thermal conductivity, the work offers guiding principles for the rational design of halide-based materials with tailored thermal transport properties for thermoelectric and thermal management applications. The findings also highlight $Cs_2SnI_6$ as a unique case where strong anharmonic scattering, rather than just low velocity, drives the record-low thermal conductivity.

1. The Mystery of the "Super-Insulator"

2. The Old Theory vs. The New Discovery

3. The "Machine Learning" Detective

4. The One Exception: The "Traffic Jam"

5. Why It Matters (According to the Paper)

More like this