Cooperative Suppression Strategy for Dual Thermal Transport Channels in Crystalline Materials

Chenhan Liu; Hao Zhang; Liujiang Zhou; Shuming Zeng; Su-Huai Wei; Ying Chen; Yu Wu

REVIEW 4 major objections 7 minor 40 references

Cooperative Suppression Strategy for Dual Thermal Transport Channels in Crystalline Materials

T0 review · 4 major / 7 minor · reviewed 2026-08-15 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A 'heavy-light, soft-stiff' atomic arrangement can suppress both particle-like and wave-like heat conduction in crystals, yielding lattice thermal conductivities near 0.15 W/mK and bypassing the usual trade-off between the two channels.

desk verdict Plausible design rule and three new candidate materials, but the central κc values rest on a computation the paper never describes, so the main claim is not yet supportable. read the letter →

arxiv 2508.17318 v1 pith:EAY5V3IR submitted 2025-08-24 cond-mat.mtrl-sci physics.comp-ph

Yu Wu , Ying Chen , Shuming Zeng , Hao Zhang , Liujiang Zhou , Chenhan Liu , Su-Huai Wei This is my paper · ORCID

classification cond-mat.mtrl-sciphysics.comp-ph PACS 63.20.-e66.70.-f

keywords heavy-lightsoft-stiffmotifultralowlatticethermalconductivityparticle-likephonontransportwave-liketunnelinghierarchicalspectrumdensity-of-statessparsityhigh-throughputmaterialsscreeningTl4SiS4

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper proposes that ultralow lattice thermal conductivity can be designed by a 'heavy-light, soft-stiff' atomic motif: heavy atoms on soft bonds create dense low-frequency vibrations that scatter heat, while light atoms on stiff bonds create sparse high-frequency vibrations that avoid the near-degeneracies behind wave-like heat tunneling. First-principles calculations and two screening passes support the claim that this hierarchical phonon spectrum suppresses both the particle-like channel $\kappa_p$ and the wave-like channel $\kappa_c$ at once, with Tl4SiS4 and Tl4GeS4 predicted near $\kappa_L \approx 0.15$ W/mK. A minimal 1D triatomic chain reproduces the same decoupling, suggesting the mechanism is not tied to a specific chemistry. A sympathetic reader would care because it offers a concrete route past the usual trade-off that has limited how low crystalline heat conduction can go.

What carries the argument

The load-bearing object is the dimensionless phonon-density-of-states sparsity $$\xi_{\mathrm{phDOS}} = \frac{\sqrt{N} - \|\tilde{g}\|_1/\|\tilde{g}\|_2}{\sqrt{N}-1},$$ where $\tilde{g}_i \equiv g(\omega_i)/\max[g(\omega)]$ is the normalized phonon density of states, $\|\tilde{g}\|_1$ its $\ell^1$ norm, and $\|\tilde{g}\|_2$ its $\ell^2$ norm. It ranges from 0 for a uniform spectrum to 1 for a maximally sparse single-peak spectrum, and it is the descriptor that selects materials with concentrated low-frequency density but sparse high-frequency branches. Screening combines this sparsity condition with a maximum acoustic frequency below 0.8 THz and, in a second pass, with structural signatures of isolated atoms beside finite clusters and a mass contrast of at least 150 atomic mass units. The same hierarchy is rationalized with a 1D triatomic chain in which mass contrast and bond-stiffness contrast independently shift the low- and high-frequency branches.

What would settle it

Compute $\kappa_c$ for Tl4SiS4 with the Wigner transport formulation from the unified theory cited in the paper, using a dense q-mesh and including four-phonon scattering, and compare the resulting $\kappa_L$ with a single-crystal measurement; a value well above 0.15 W/mK would falsify the central claim.

Watch

Extended reading notes

Core claim

The central claim is that the antagonism between particle-like and wave-like phonon transport is not a fundamental lower bound but a spectral-design problem. In the proposed motif, the heavy, softly bonded species supply the scattering channels that cut $\kappa_p$, while the light, stiffly bonded species push optical branches to high frequency and keep them sparse, so few phonon pairs have the small frequency differences needed for coherent tunneling; hence $\kappa_c$ stays low. The paper identifies Rb6Re6S8I8 as a natural example, introduces a phonon-density sparsity metric to screen thousands of compounds, and reports Tl4SiS4 ($\kappa_p=0.10$, $\kappa_c=0.06$ W/mK) and Tl4GeS4 ($\kappa_p=0.09$, $\kappa_c=0.06$ W/mK) as candidate materials in which both channels are simultaneously suppressed. It also shows in a 1D triatomic chain that mass contrast and bond-stiffness contrast move the low- and high-frequency branches independently, which is the mechanism's core.

Load-bearing premise

The load-bearing premise is that the reported values of the wave-like (coherent) channel $\kappa_c$ are correct, yet the paper does not state the equation or implementation used to compute that channel; if those numbers are wrong, the claim of simultaneous suppression loses its foundation.

Editorial extensions

If this is right

If the motif is sufficient, crystals built from isolated heavy cations and stiff covalent clusters should reach lattice thermal conductivities below 0.2 W/mK without needing strong anharmonicity from soft frameworks alone.
The screening descriptors (acoustic cutoff below 0.8 THz, sparsity above 0.4, mass contrast above 150 amu) can be applied to any phonon database, making the search for ultralow-thermal-conductivity materials largely structural.
The 1D chain result implies that one-dimensional or layered materials with the same mass and bond hierarchy should show the same dual suppression, extending the claim beyond the three-dimensional compounds computed.
Because the suppression is built into an ordered crystal rather than into disorder, the candidate materials may preserve clean electronic bands, which matters for thermoelectric applications.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

One testable extension the authors do not spell out: the stiffness ratio, not just the mass ratio, could be tuned continuously by pressure, so hydrostatic pressure should change the balance between $\kappa_p$ and $\kappa_c$ in a predictable way.
The sparsity metric could be applied to existing databases of molecular and hybrid crystals, where heavy counterions and stiff covalent groups are common, to see whether this is already a hidden route to ultralow conductivity.
If the predicted values are right, a straightforward single-crystal measurement of Tl4SiS4 or Tl4GeS4 would be a sharp test, since no nanostructuring or alloying would be needed to see the effect.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Colleague,

The paper has a genuinely interesting idea: use heavy atoms on soft bonds to generate dense low-frequency modes that suppress the particle-like channel, and light atoms on stiff bonds to create sparse high-frequency branches that suppress the wave-like channel. That qualitative picture is worth taking seriously. The paper also adds something concrete: a sparsity metric ξphDOS, a screening pipeline over PhononDB and Materials Project, and three new Tl–S candidates (Tl4SiS4, Tl4GeS4, Tl2S5) with predicted κL around 0.15 W/mK. If those numbers hold, it would break the usual κp–κc trade-off, so the significance is real within phonon engineering.

The problem is that the central quantitative claim is not verifiable from the text. κc contributes 29–40% of total κL in the headline compounds, but the Methods section only says the thermal conductivity is computed with ShengBTE, which solves the particle-like Boltzmann equation. There is no equation, citation, or code for the Wigner/coherent channel—no mention of off-diagonal velocity matrix elements, no description of the implementation, nothing. The 1D triatomic chain model is illustrative but never computes κc for the chain. This is not a minor documentation gap; if κc is miscomputed, the "co-suppression" conclusion loses its foundation. The stress-test note is on target.

There are also smaller internal inconsistencies: the text says "top five" but Table 1 lists eight materials, and the abstract highlights two candidates while the screening section names three Tl–S compounds. The ξphDOS threshold was calibrated on Rb6Re6S8I8, the same compound used later for validation. That is not fatal—the new candidates are independent predictions—but it weakens the universality claim.

What the paper does well is articulate a mechanism and back it with a descriptor that captures something real: sparse high-frequency phDOS correlates with suppressed coherent transport. The structural analysis of Tl4SiS4 makes the heavy-light/soft-stiff picture concrete, and the physics is plausible and falsifiable.

Who should read this? People hunting for ultralow-κ or thermoelectric candidates and anyone interested in phonon design rules. It deserves a serious referee, but only with the firm requirement that the authors disclose the κc calculation in full—formula, implementation, convergence, and ideally code or data. I would not cite the specific κc values until that is out.

Send it to review, but conditional on major revision.

Referee Report

4 major / 7 minor

Summary. The manuscript proposes a 'heavy-light, soft-stiff' structural design principle intended to circumvent the usual trade-off between particle-like (κp) and wave-like (κc) lattice thermal conductivity in crystals. Using TDEP and ShengBTE calculations, the authors characterize phonon dispersions and two-channel conductivities of Rb6Re6S8I8 and MgAgSb, introduce a phDOS sparsity metric ξphDOS, screen PhononDB and Materials Project databases with frequency, mass, and structural descriptors, and identify candidate compounds—most prominently Tl4SiS4 and Tl4GeS4—with predicted κp ≈ 0.1 W/mK and κc ≈ 0.06 W/mK. A 1D triatomic chain model is used to illustrate the proposed spectral decoupling mechanism.

Significance. If the results hold, the design principle would be a valuable contribution: it offers a concrete, physically motivated route to materials whose total lattice conductivity falls below the usual κp–κc trade-off, and the use of simple structural descriptors makes the screening strategy transferable. The paper is also commendable for testing the idea on known reference compounds and for providing differential frequency-resolved κp and κc analyses. However, the quantitative case rests entirely on values of κc that are not reproducible from the manuscript as written, so the significance is conditional on a complete specification of the coherent-transport calculation.

major comments (4)

[Methods, p. 2] The manuscript nowhere gives the formula, implementation, or code used for the coherent wave-like channel κc, although this channel is the central quantity behind the claim of breaking the κp–κc trade-off. Sheng-BTE (refs. [26,27]) solves the particle-like Boltzmann equation and does not by itself produce a Wigner-type κc; the distinction matters because Table I lists κc = 0.05–0.08 W/mK, which is 29–40% of the predicted total for Tl4SiS4 and Tl4GeS4. Please add the Wigner transport equations used (e.g., the unified formalism of ref. [9]), the treatment of off-diagonal velocity matrix elements and phonon linewidths, the numerical implementation or code, and q-mesh/convergence benchmarks; without these, the central quantitative claim cannot be verified.
[Fig. 2(b) and Table I] The text says 'the top five candidates listed in Table 1' and Fig. 2(b) labels the output 'Materials with Ultralow κL (5)', yet Table I contains eight materials, and later the text says that the Materials-Project screen yields '226 qualified materials, including three Tl–S-based compounds ... as listed in Table 1'. The manuscript must state explicitly which candidates came from the PhononDB stage and which from the Materials Project stage, whether the five and eight are compatible as 5 + 3, and provide κL totals for all entries; otherwise the screening results are not auditable.
[Eq. (1) and Fig. 2(b)] The screening thresholds (ξphDOS > 0.4, ωa-max < 0.8 THz, Mmin ≤ 35, ΔM ≥ 150, bond-distance factor 1.15) are introduced with benchmark values from Rb6Re6S8I8 and MgAgSb, and Rb6Re6S8I8 (ref. [2]) is also the motivating example used to justify the metric; this creates a partial circularity. A sensitivity analysis with respect to these thresholds, plus a test on a material not used to calibrate the descriptors, is needed to show that the identification of ultralow-κc candidates is robust rather than tuned to a single known example.
[Fig. 2(d)] The 'normalized κc (κp/0.2×κc)' rescaling is not a proper statistical normalization; because it multiplies the quantity of interest by a factor that correlates with the opposite channel, the apparent decreasing trend with ξphDOS is not by itself evidence of a physical correlation. Report the raw data, the materials included, the correlation coefficient and significance, and preferably also a partial-correlation analysis that controls for κp; otherwise this figure does not support the central sparsity-coherence claim.

minor comments (7)

[Main text and Table I] The text states that Tl4SiS4 and Tl4GeS4 achieve κL = 0.15 W/mK, but Table I implies κp+κc = 0.16 W/mK for Tl4SiS4; clarify the total and add a κL column to Table I.
[Methods, p. 2] The computational methods do not state the temperature at which κp and κc are evaluated; the figures say 300 K, but this should be in the Methods text.
[Fig. 3(a)–(c)] The 1D triatomic-chain discussion is purely qualitative; it shows how masses and springs change the dispersion, but no κp or κc is computed for the chain, so the phrase 'demonstrates the generality' in the abstract is stronger than the evidence.
[Methods, p. 2] The Methods would benefit from reporting the TDEP expansion order (harmonic + third-order, or fourth-order), supercell sizes, and q-mesh convergence tests; the current description is too brief to reproduce κp either.
[Fig. 2(c)] The caption says circle sizes encode primitive-cell atom counts, but no legend or size scale is provided; add one.
[Eq. (1)] In Eq. (1), define N explicitly as the number of frequency bins and state the smearing used to build g(ω).
[General] There is no data-availability or code-availability statement; given the centrality of the κc values, please add one.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the headline predictions are independent first-principles transport calculations; the only self-referential element is the motivating Rb6Re6S8I8 example and descriptor threshold calibration, which does not force the reported κp/κc values.

full rationale

The design principle is inferred from Rb6Re6S8I8 (ref 2, by overlapping authors), and the ξphDOS threshold is benchmarked against that example, but the central predictions for Tl4SiS4 and Tl4GeS4 are computed from first-principles IFCs and phonon transport calculations, not derived from the descriptor. Equation (1) defines ξphDOS directly from the phonon density of states and does not contain κp or κc, so screening on ξphDOS does not by construction produce the reported transport numbers. The 1D triatomic chain is an illustrative model with parameters chosen to exhibit the desired dispersion features; it is not used to compute κc or to predict the candidate materials. No fitted parameter is renamed as a prediction, and no uniqueness theorem or ansatz is imported from the authors' prior work to forbid alternatives. The absence of a disclosed Wigner/coherent-channel formula in the methods is a reproducibility and correctness gap, but the missing equation is not shown to reduce the prediction to the input, so it is not circularity within the definition used here. The new compounds are independent of the calibration example, and the quantitative claim stands or falls on first-principles calculations that are self-contained rather than definitionally equivalent to the screening input.

Assumptions & free parameters 4 free parameters · 4 assumptions · 0 invented entities

The central claim rests on hand-picked screening thresholds and an unstated κc method. No new physical entity is introduced; the only new construct is the ξphDOS metric, which is a mathematical descriptor. The ledger shows that the paper's contribution is a proposed mechanism plus a screen, not a parameter-free derivation.

free parameters (4)

ξphDOS threshold > 0.4 = 0.4
Sparsity cutoff in the PhononDB screening workflow; chosen to retain materials similar to Rb6Re6S8I8 (ξ=0.79) and exclude MgAgSb (ξ=0.26), with no sensitivity analysis.
Maximum acoustic frequency cutoff ωa-max < 0.8 THz = 0.8 THz
Used to ensure low acoustic group velocities for κp suppression; threshold is ad hoc and no justification is given.
Mass contrast filters Mmin ≤ 35 and Mmax-Mmin ≥ 150 = 35 and 150 (atomic mass units)
Applied in the Materials Project screening to impose heavy-light mass contrast; values hand-picked to isolate Tl and similar compounds.
Bond-absence distance factor 1.15(Rr+Rn) = 1.15
Used to identify isolated atoms and finite clusters in the structural descriptor; no calibration or physical derivation is provided.

assumptions (4)

domain assumption DFT with the PBEsol functional gives accurate forces and phonons for the soft, heavy-atom compounds studied.
All results depend on PBEsol anharmonic force constants; no validation against experiment or higher-level theory is shown.
domain assumption TDEP force constants from 30 ps AIMD at 300 K capture the anharmonicity relevant to κp and κc.
The transport pipeline rests on converged AIMD-derived IFCs, but convergence with simulation length and supercell size is not demonstrated.
domain assumption A Wigner/coherent transport decomposition (unified theory, ref [9]) is the correct way to separate κp and κc, and it was applied correctly.
The paper reports κc values but does not give the formula, implementation, or code; this is the weakest unstated input.
ad hoc to paper The sparsity metric ξphDOS and structural descriptors are predictive of low coherent thermal conductivity.
The metric is defined in Eq. (1) and calibrated using only a few known compounds; no large-scale blind test is reported.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Cooperative Suppression Strategy for Dual Thermal Transport Channels in Crystalline Materials." pith.science (2026). https://pith.science/paper/EAY5V3IR

@misc{pith2026250817318,
  author       = {Pith},
  title        = {Pith review of: Cooperative Suppression Strategy for Dual Thermal Transport Channels in Crystalline Materials},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/EAY5V3IR}},
  note         = {Machine review of arXiv:2508.17318}
}

abstract

We propose a novel design principle for achieving ultralow thermal conductivity in crystalline materials via a "heavy-light and soft-stiff" structural motif. By combining heavy and light atomic species with soft and stiff bonding networks, both particle-like ($\kappa_p$) and wave-like ($\kappa_c$) phonon transport channels are concurrently suppressed. First-principles calculations show that this architecture induces a hierarchical phonon spectrum: soft-bonded heavy atoms generate dense low-frequency modes that enhance scattering and reduce $\kappa_p$, while stiff-bonded light atoms produce sparse high-frequency optical branches that disrupt coherence and lower $\kappa_c$. High-throughput screening identifies Tl$_4$SiS$_4$ ($\kappa_p$ = 0.10, $\kappa_c$ = 0.06 W/mK) and Tl$_4$GeS$_4$ ($\kappa_p$ = 0.09, $\kappa_c$ = 0.06 W/mK) as representative candidates with strongly suppressed transport in both channels. A minimal 1D triatomic chain model further demonstrates the generality of this mechanism, offering a new paradigm for phonon engineering beyond the conventional $\kappa_p$-$\kappa_c$ trade-off.

Figures

Figures reproduced from arXiv: 2508.17318 by the authors.

**Figure 2.** FIG. 2. (a) Characteristic phDOS profile for materials enabl [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

**Figure 3.** (b,c). It can also be seen from the differential of cumulative κc that the high-frequency optical branches exhibit relatively small contribution to κc. In summary, this work advances the fundamental understanding of phonon engineering by establishing three key insights: (1) The ”heavy-light and softstiff” design principle represents a universal strategy to decouple κp and κc suppression mechanisms, overcoming thei… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

40 extracted references · 35 canonical work pages

[9]

Y. Wu, A. Huang, L. Ji, J. Ji, Y. Ding, and L. Zhou, Origin of intrinsically low lattice thermal conductivity i n solids, J. Phys. Chem. Lett. , 11525 (2024)

work page 2024
[2]

Finite Cluster Criterion /g7/g3 atom r, /g5 nearest neighbors n D(r, n) ≥ 1.15 (Rr+Rn) A connected atom group C where: /g5 i, j /g143/g3 C ğpath-connected /g5 i /g143/g3 C, k/g144 CğD(i, k) ≥ 1.15 (Ri+Rk) 139031 11478 906 226 First Principles Validation: Tl 2S5, Tl 4SiS4, Tl 4GeS4 Г NP 9 NP 13 NP 14 XГ 0 4 8 16 12 STl Si κp (ω) κc (ω) κc (ω) [W/(mK)/THz] ...

work page
[1]

Isolated Atom Criterion

work page
[3]

Z. Zeng, X. Shen, R. Cheng, O. Perez, N. Ouyang, Z. Fan, P. Lemoine, B. Raveau, E. Guilmeau, and Y. Chen, Pushing thermal conductivity to its lower limit in crystals with simple structures, Nat. Commun. 15, 10.1038/s41467-024-46799-3 (2024)

work page doi:10.1038/s41467-024-46799-3 2024
[4]

Xiong, H

W. Xiong, H. Huang, Y. Wu, X. Xu, G. Li, Z. Gu, and S. Zeng, Decoupling ultralow coherent and particle-like phonon transport via bonding hierarchy in soft superionic crystals, Adv. Sci. 10.1002/advs.202506807 (2025)

work page doi:10.1002/advs.202506807 2025
[5]

J. Li, X. Li, Y. Zhang, J. Zhu, E. Zhao, M. Kofu, K. Nakajima, M. Avdeev, P.-F. Liu, J. Sui, H. Zhao, F. Wang, and J. Zhang, Crystal-liquid duality driven ul- tralow two-channel thermal conductivity in α-mgagsb, Appl. Phys. Rev. 11, 10.1063/5.0173680 (2024)

work page doi:10.1063/5.0173680 2024
[6]

Zheng, C

J. Zheng, C. Lin, C. Lin, G. Hautier, R. Guo, and B. Huang, Unravelling ultralow thermal conductivity in perovskite cs2agbibr6: dominant wave-like phonon tun- nelling and strong anharmonicity, npj Comput. Mater. 10, 10.1038/s41524-024-01211-y (2024)

work page doi:10.1038/s41524-024-01211-y 2024
[7]

Y. Wu, Y. Chen, Z. Fang, Y. Ding, Q. Li, K. Xue, H. Shao, H. Zhang, and L. Zhou, Ultralow lat- tice thermal transport and considerable wave-like phonon tunneling in chalcogenide perovskite bazrs3, J. Phys. Chem. Lett. 14, 11465 (2023)

work page 2023

Show all 40 references

[8]

L. Ji, A. Huang, Y. Huo, Y.-m. Ding, S. Zeng, Y. Wu, and L. Zhou, Inﬂuence of four-phonon scattering and wave- like phonon tunneling eﬀects on the thermal transport properties of tlbise2, Phys. Rev. B 109, 214307 (2024)

2024
[10]

Di Lucente, M

E. Di Lucente, M. Simoncelli, and N. Marzari, Crossover from boltzmann to wigner thermal trans- port in thermoelectric skutterudites, Phys. Rev. Res. 5, 10.1103/physrevresearch.5.033125 (2023)

2023 doi
[11]

Simoncelli, N

M. Simoncelli, N. Marzari, and F. Mauri, Uniﬁed theory of thermal transport in crystals and glasses, Nat. Phys. 15, 809 (2019)

2019
[12]

Zhang, T

Y. Zhang, T. Frauenheim, T. Dumitric˘ a, and Z. Tong, Lonesome ag atoms drive ultralow thermal con- ductivity in argyrodite ag8snse6, PRX Energy 4, 10.1103/prxenergy.4.013012 (2025)

2025 doi
[13]

X. Zeng, G. Niu, X. Wang, J. Jiang, L. Sui, Y. Zhang, A. Chen, M. Jin, and K. Yuan, Strain-engineered ther- moelectric performance in superatomic semiconductor re6se8cl2: The role of four-phonon scattering and coher- ent phonons, J. Phys. Chem. Lett. 16, 3122 (2025)

2025
[14]

T. M. Tritt, Thermal conductivity: theory, properties, and applicatio ns (Springer US, 2004)

2004
[15]

Jia-Jun, F

S. Jia-Jun, F. Teng, F. Tie-Zheng, X. Jia- Zhan, Z. Xin-Bing, and Z. Tie-Jun, Lattice ther- mal conductivity in thermoelectric materials, J. Inorg. Mater. 34, 260 (2019)

2019
[16]

C. Liu, C. Wu, T. Song, Y. Zhao, J. Yang, P. Lu, G. Zhang, and Y. Chen, An eﬃcient strategy for searching high lattice thermal conductivity materials, ACS Appl. Energy Mater. 5, 15356 (2022)

2022
[17]

Christensen, A

M. Christensen, A. B. Abrahamsen, N. B. Christensen, F. Juranyi, N. H. Andersen, K. Lefmann, J. An- dreasson, C. R. H. Bahl, and B. B. Iversen, Avoided crossing of rattler modes in thermoelectric materials, Nat. Mater. 7, 811 (2008)

2008
[18]

Tadano, Y

T. Tadano, Y. Gohda, and S. Tsuneyuki, Impact of rattlers on thermal conductivity of a thermoelectric clathrate: A ﬁrst-principles study, Phys. Rev. Lett. 114, 10.1103/physrevlett.114.095501 (2015)

2015 doi
[19]

M. K. Jana, K. Pal, A. Warankar, P. Mandal, U. V. Waghmare, and K. Biswas, Intrinsic rattler- 7 induced low thermal conductivity in zintl type tlinte2, J. Am. Chem. Soc. 139, 4350 (2017)

2017
[20]

M. Wu, S. Yan, and J. Ren, Hierarchy-boosted fun- nel learning for identifying semiconductors with ultralow lattice thermal conductivity, npj Comput. Mater. 11, 10.1038/s41524-025-01583-9 (2025)

2025 doi
[21]

M. D. Nielsen, V. Ozolins, and J. P. Heremans, Lone pair electrons minimize lattice thermal conductivity, Energy Environ. Sci. 6, 570 (2013)

2013
[22]

E. J. Skoug and D. T. Morelli, Role of lone-pair electron s in producing minimum thermal conductivity in nitrogen- group chalcogenide compounds, Phys. Rev. Lett. 107, 10.1103/physrevlett.107.235901 (2011)

2011 doi
[23]

S. Lee, K. Esfarjani, T. Luo, J. Zhou, Z. Tian, and G. Chen, Resonant bonding leads to low lattice thermal conductivity, Nat. Commun. 5, 3525 (2014)

2014
[24]

J. Ji, Q. Tang, M. Yao, H. Yang, Y. Jin, Y. Zhang, J. Xi, D. J. Singh, J. Yang, and W. Zhang, Functional-unit- based material design: Ultralow thermal conductivity in thermoelectrics with linear triatomic resonant bonds, J. Am. Chem. Soc. 144, 18552 (2022)

2022
[25]

Kresse and J

G. Kresse and J. Furthm¨ uller, Eﬃcient iterative schem es forab initiototal-energy calculations using a plane-wave basis set, Phys. Rev. B 54, 11169 (1996)

1996
[26]

J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vy- drov, G. E. Scuseria, L. A. Constantin, X. Zhou, and K. Burke, Restoring the density-gradient expansion for exchange in solids and surfaces, Phys. Rev. Lett. 100, 136406 (2008)

2008
[27]

Hellman and I

O. Hellman and I. A. Abrikosov, Temperature-dependent eﬀective third-order interatomic force constants from ﬁrs t principles, Phys. Rev. B 88, 144301 (2013)

2013
[28]

W. Li, J. Carrete, N. A. Katcho, and N. Mingo, Sheng- bte: A solver of the boltzmann transport equation for phonons, Comput. Phys. Commun. 185, 1747 (2014)

2014
[29]

Z. Han, X. Yang, W. Li, T. Feng, and X. Ruan, Four- phonon: An extension module to shengbte for computing four-phonon scattering rates and thermal conductivity, Comput. Phys. Commun. 270, 108179 (2022)

2022
[30]

Togo and I

A. Togo and I. Tanaka, First principles phonon calcula- tions in materials science, Scr. Mater. 108, 1 (2015)

2015
[31]

F. Li, X. Liu, J. Yang, X. Wang, Y. Yang, N. Ma, L. Chen, and L. Wu, Long-range anion correlations mediating dynamic anharmonicity and contributing to glassy thermal conductivity in well-ordered k2ag4se3, Small 21, 2409524 (2024)

2024
[32]

J. Li, L. Wei, Z. Ti, J. Zhu, W. Zhai, H. Li, and P.-F. Liu, Ultralow two-channel thermal conductivity in aikinite, Mat. Today Phys. 45, 101454 (2024)

2024
[33]

Xie, P.-F

Q.-Y. Xie, P.-F. Liu, J.-J. Ma, L.-M. Wu, K.-W. Zhang, and B.-T. Wang, Microscopic mechanisms of glasslike lattice thermal conductivity in tetragonal α-cscu5se3, Phys. Rev. B 108, 014302 (2023)

2023
[34]

Pandey, M.-H

T. Pandey, M.-H. Du, D. S. Parker, and L. Lind- say, Origin of ultralow phonon transport and strong anharmonicity in lead-free halide perovskites, Mater. Today Phys. 28, 100881 (2022)

2022
[35]

Z. Tong, T. Dumitric˘ a, and T. Frauenheim, In- trinsically low thermal conductivity in the most lithium-rich binary stannide crystalline li5sn, J. Phys. Chem. Lett. 14, 8139 (2023)

2023
[36]

Y. Xia, V. Ozoli¸ nˇ s, and C. Wolverton, Microscopic mechanisms of glasslike lattice thermal transport in cubic cu12sb4s13 tetrahedrites, Phys. Rev. Lett. 125, 10.1103/physrevlett.125.085901 (2020)

2020 doi
[37]

Z. Tong, A. Pecchia, C. Yam, T. Dumitric˘ a, and T. Frauenheim, Glass-like transport dominates ultralow lattice thermal conductivity in modular crystalline bi4o4secl2, Nano Lett. 23, 9468 (2023)

2023
[38]

Z. Tong, Y. Zhang, T. Frauenheim, and T. Du- mitric˘ a, Layered germanium–selenium compounds as phonon–glass electron–crystals: A path- way to enhance the thermoelectric performance, Nano Lett. 25, 7283 (2025)

2025
[39]

X. Shen, N. Ouyang, Y. Huang, Y. Tung, C. Yang, M. Faizan, N. Perez, R. He, A. Sotnikov, K. Willa, C. Wang, Y. Chen, and E. Guilmeau, Amorphous- like ultralow thermal transport in crystalline argyrodite cu7ps6, Adv. Sci. 11, 2400258 (2024)

2024
[40]

A. Jain, S. P. Ong, G. Hautier, W. Chen, W. D. Richards, S. Dacek, S. Cholia, D. Gunter, D. Skinner, G. Ceder, and K. A. Persson, Commentary: The materials project: A materials genome approach to accelerating materials innovation, APL Mater. 1, 011002 (2013)

2013

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.

[1] [9]

Y. Wu, A. Huang, L. Ji, J. Ji, Y. Ding, and L. Zhou, Origin of intrinsically low lattice thermal conductivity i n solids, J. Phys. Chem. Lett. , 11525 (2024)

work page 2024

[2] [2]

Finite Cluster Criterion /g7/g3 atom r, /g5 nearest neighbors n D(r, n) ≥ 1.15 (Rr+Rn) A connected atom group C where: /g5 i, j /g143/g3 C ğpath-connected /g5 i /g143/g3 C, k/g144 CğD(i, k) ≥ 1.15 (Ri+Rk) 139031 11478 906 226 First Principles Validation: Tl 2S5, Tl 4SiS4, Tl 4GeS4 Г NP 9 NP 13 NP 14 XГ 0 4 8 16 12 STl Si κp (ω) κc (ω) κc (ω) [W/(mK)/THz] ...

work page

[3] [1]

Isolated Atom Criterion

work page

[4] [3]

Z. Zeng, X. Shen, R. Cheng, O. Perez, N. Ouyang, Z. Fan, P. Lemoine, B. Raveau, E. Guilmeau, and Y. Chen, Pushing thermal conductivity to its lower limit in crystals with simple structures, Nat. Commun. 15, 10.1038/s41467-024-46799-3 (2024)

work page doi:10.1038/s41467-024-46799-3 2024

[5] [4]

Xiong, H

W. Xiong, H. Huang, Y. Wu, X. Xu, G. Li, Z. Gu, and S. Zeng, Decoupling ultralow coherent and particle-like phonon transport via bonding hierarchy in soft superionic crystals, Adv. Sci. 10.1002/advs.202506807 (2025)

work page doi:10.1002/advs.202506807 2025

[6] [5]

J. Li, X. Li, Y. Zhang, J. Zhu, E. Zhao, M. Kofu, K. Nakajima, M. Avdeev, P.-F. Liu, J. Sui, H. Zhao, F. Wang, and J. Zhang, Crystal-liquid duality driven ul- tralow two-channel thermal conductivity in α-mgagsb, Appl. Phys. Rev. 11, 10.1063/5.0173680 (2024)

work page doi:10.1063/5.0173680 2024

[7] [6]

Zheng, C

J. Zheng, C. Lin, C. Lin, G. Hautier, R. Guo, and B. Huang, Unravelling ultralow thermal conductivity in perovskite cs2agbibr6: dominant wave-like phonon tun- nelling and strong anharmonicity, npj Comput. Mater. 10, 10.1038/s41524-024-01211-y (2024)

work page doi:10.1038/s41524-024-01211-y 2024

[8] [7]

Y. Wu, Y. Chen, Z. Fang, Y. Ding, Q. Li, K. Xue, H. Shao, H. Zhang, and L. Zhou, Ultralow lat- tice thermal transport and considerable wave-like phonon tunneling in chalcogenide perovskite bazrs3, J. Phys. Chem. Lett. 14, 11465 (2023)

work page 2023