pith. sign in

arxiv: 2510.07934 · v6 · submitted 2025-10-09 · ❄️ cond-mat.mtrl-sci · physics.comp-ph

Higher-Order-Phonon Scattering Governs Targeted Control of Heat Conduction in Bulk Boron Arsenide

Tianhao Li , Yangjun Qin , Dongkai Pan , Shixian Liu , Han Meng , Nuo Yang This is my paper

Pith reviewed 2026-05-18 09:16 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci physics.comp-ph
keywords boron arsenidethermal conductivityphonon scatteringfour-phonon processestargeted excitationheat conduction modulationboltzmann transport
0
0 comments X

The pith

Four-phonon scattering turns targeted phonon excitation into predominantly suppressive control of thermal conductivity in bulk boron arsenide.

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

The paper establishes that targeted phonon excitation can dynamically modulate thermal conductivity in a three-dimensional bulk material such as boron arsenide. Within a three-phonon scattering model the effect remains weak and bidirectional at room temperature. Adding four-phonon scattering changes the outcome to strong suppression that varies sharply with excitation frequency. The change occurs because four-phonon processes raise the overall scattering background and amplify the additional scattering imposed on low-frequency heat-carrying modes. A reader would care because this route offers reversible, in-situ tuning of heat flow without any permanent alteration to the material structure.

Core claim

Based on first-principles calculations and phonon Boltzmann transport analysis, targeted phonon excitation modulates the thermal conductivity of bulk BAs in a strongly frequency-dependent manner. Within the three-phonon-only framework the modulation at 300 K is weak but clearly bidirectional. However, once four-phonon scattering is included the modulation changes qualitatively to a predominantly suppressive behavior. In the combined three-phonon plus four-phonon framework the strongest suppression occurs at 20.5 THz, where the relative thermal conductivity decreases to 0.828 and 0.415 for excitation intensities of 5 and 25, respectively. Four-phonon scattering raises the intrinsic scattering

What carries the argument

The combined three-phonon plus four-phonon scattering framework that raises background scattering rates and amplifies excitation-induced scattering of low-frequency phonons under targeted drive.

If this is right

  • Modulation becomes predominantly suppressive rather than bidirectional once four-phonon scattering is included.
  • Suppression reaches its maximum at 20.5 THz with relative thermal conductivity falling to 0.415 at the higher excitation intensity.
  • Four-phonon processes increase scattering of low-frequency heat-carrying phonons beyond the rise in background rates.
  • The approach enables reversible in-situ regulation of thermal transport in bulk three-dimensional crystals without structural modification.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same frequency-selective suppression could appear in other bulk materials that exhibit strong four-phonon scattering.
  • Experimental mapping of thermal conductivity versus excitation frequency around 20.5 THz would directly test the predicted peak.
  • If the mechanism holds, analogous targeted excitation might modulate other phonon-limited properties such as thermoelectric figure of merit.

Load-bearing premise

Targeted phonon excitation at the stated intensities can be realized experimentally while leaving the lattice temperature and other scattering channels essentially unchanged.

What would settle it

A measurement of thermal conductivity in boron arsenide under selective excitation at 20.5 THz that shows an increase or no change rather than the predicted drop to 0.415 relative conductivity at intensity 25 would falsify the central claim.

Figures

Figures reproduced from arXiv: 2510.07934 by Dongkai Pan, Han Meng, Nuo Yang, Shixian Liu, Tianhao Li, Yangjun Qin.

Figure 3
Figure 3. Figure 3: (a) The spectral thermal conductivity and (b) the relative thermal conductivity spectrum normalized by total thermal conductivity of BAs under the excitation of phonons of 9.2–9.3 THz with intensities of 5 and 25 at 300 K. To understand the different impacts on fundamental phonon scattering induced by the excitations of phonons of 9.2–9.3 THz at different intensities, the contributions of different three-p… view at source ↗
Figure 4
Figure 4. Figure 4: T [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
read the original abstract

Conventional approaches for modulating thermal conductivity usually rely on structural modifications and therefore cannot achieve reversible in situ regulation. Targeted phonon excitation has recently emerged as a promising strategy for dynamically tuning thermal transport, but its applicability has so far been demonstrated mainly in two-dimensional systems. Here, we extend this strategy to a three-dimensional bulk material by taking boron arsenide (BAs) as a representative example. Based on first-principles calculations and phonon Boltzmann transport analysis, we show that targeted phonon excitation modulates the thermal conductivity of bulk BAs in a strongly frequency-dependent manner. Within the three-phonon-only framework, the modulation at 300 K is weak but clearly bidirectional. However, once four-phonon scattering is included, the modulation changes qualitatively to a predominantly suppressive behavior. In the combined three-phonon plus four-phonon (3ph+4ph) framework, the strongest suppression occurs at 20.5 THz, where the relative thermal conductivity decreases to 0.828 and 0.415 for excitation intensities of 5 and 25, respectively. By comparing the 3ph-only and 3ph+4ph results, we show that four-phonon scattering plays a decisive role in determining the net modulation effect by raising the intrinsic scattering background and promoting a more systematic excitation-induced increase in the scattering of low-frequency heat-carrying phonons.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The paper claims that targeted phonon excitation provides a reversible way to modulate thermal conductivity in bulk boron arsenide (BAs), extending prior 2D work to 3D. First-principles calculations combined with the phonon Boltzmann transport equation show weak bidirectional modulation in the three-phonon-only framework, but predominantly suppressive behavior once four-phonon scattering is included. The strongest suppression occurs at 20.5 THz, with relative thermal conductivity dropping to 0.828 and 0.415 for excitation intensities of 5 and 25, respectively; the authors attribute the qualitative change to four-phonon processes raising the intrinsic scattering background and enhancing scattering of low-frequency heat carriers.

Significance. If the quantitative results survive scrutiny of the non-equilibrium assumptions, the work would be significant for demonstrating frequency-dependent, in-situ thermal control in a bulk 3D material without structural changes. It highlights the decisive role of higher-order (four-phonon) scattering in determining the net modulation sign and magnitude, providing concrete predictions that could inform experiments on dynamic heat management.

major comments (1)
  1. [phonon Boltzmann transport analysis and results] BTE implementation under targeted excitation (described in the phonon Boltzmann transport analysis and results): the reported suppression ratios (0.828 and 0.415 at 20.5 THz for intensities 5 and 25 in the 3ph+4ph case) are obtained by selectively increasing occupations near the target frequency while holding lattice temperature fixed at 300 K and leaving all other scattering rates unchanged. This assumption is load-bearing for the central claim that four-phonon scattering governs the net suppressive effect, because anharmonic decay of the pumped phonons would deposit energy, raising the effective temperature and increasing intrinsic 3ph and 4ph rates across the spectrum; the 3ph-only vs. 3ph+4ph comparison therefore risks misattributing part of the suppression to temperature shifts rather than four-phonon processes alone.
minor comments (2)
  1. [Abstract and numerical results] No error bars, convergence tests with respect to q-grid density, or details on how the non-equilibrium excitation is numerically implemented in the BTE solver are provided for the quoted relative conductivity values.
  2. [Discussion] The manuscript would benefit from a brief discussion or supplementary calculation estimating the temperature rise expected from the stated excitation intensities to quantify the validity of the fixed-T approximation.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their insightful comments on our manuscript. We address the major comment point by point below.

read point-by-point responses

  1. Referee: [phonon Boltzmann transport analysis and results] BTE implementation under targeted excitation (described in the phonon Boltzmann transport analysis and results): the reported suppression ratios (0.828 and 0.415 at 20.5 THz for intensities 5 and 25 in the 3ph+4ph case) are obtained by selectively increasing occupations near the target frequency while holding lattice temperature fixed at 300 K and leaving all other scattering rates unchanged. This assumption is load-bearing for the central claim that four-phonon scattering governs the net suppressive effect, because anharmonic decay of the pumped phonons would deposit energy, raising the effective temperature and increasing intrinsic 3ph and 4ph rates across the spectrum; the 3ph-only vs. 3ph+4ph comparison therefore risks misattributing part of the suppression to temperature shifts rather than four-phonon processes alone.

    Authors: We appreciate the referee pointing out the implications of our fixed lattice temperature assumption in the BTE calculations. Our approach of selectively increasing occupations at the target frequency while keeping the lattice temperature fixed at 300 K is designed to isolate the effects of the non-equilibrium phonon distribution on the thermal transport properties. This enables a clear demonstration of how the inclusion of four-phonon scattering qualitatively changes the modulation from bidirectional to predominantly suppressive by raising the background scattering and enhancing scattering of low-frequency phonons. We recognize that anharmonic decay of the excited phonons would in practice cause heating, thereby increasing scattering rates throughout the spectrum. Nevertheless, the fixed-temperature model provides a useful baseline for understanding the role of higher-order scattering in the modulation mechanism. In the revised manuscript, we will include additional discussion clarifying this modeling choice and its limitations, as well as suggestions for experimental verification that accounts for possible temperature changes. revision: partial

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained first-principles BTE solution

full rationale

The paper computes phonon lifetimes and thermal conductivity via first-principles methods, then solves the phonon Boltzmann transport equation after selectively increasing occupations of modes near a target frequency to model excitation. The reported relative conductivity values (e.g., 0.828 and 0.415 at 20.5 THz) are direct numerical outputs of this procedure under the stated assumptions, not quantities fitted to data within the paper or defined in terms of the result itself. No self-citations are invoked as load-bearing uniqueness theorems, no ansatz is smuggled via prior work, and no renaming of known results occurs. The 3ph vs. 3ph+4ph comparison follows from including additional scattering channels computed independently. The derivation chain therefore remains non-circular and externally falsifiable against independent phonon calculations or experiments.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claims rest on the validity of the phonon Boltzmann transport equation under targeted excitation and on the accuracy of first-principles force constants for both three- and four-phonon processes; no new particles or forces are postulated.

free parameters (1)
  • excitation intensity
    Values of 5 and 25 are chosen to illustrate the modulation strength; these are external drive parameters rather than fitted material constants.
axioms (2)
  • domain assumption Phonon Boltzmann transport equation remains applicable under selective mode excitation
    Invoked throughout the analysis to compute the modified thermal conductivity.
  • domain assumption Four-phonon scattering matrix elements can be computed accurately from first-principles force constants
    Required for the 3ph+4ph framework to produce the reported suppression.

pith-pipeline@v0.9.0 · 5791 in / 1372 out tokens · 28309 ms · 2026-05-18T09:16:39.722013+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

43 extracted references · 43 canonical work pages

  1. [1]

    Anitha, A

    A. Anitha, A. Hemalatha, P. Udhayakumar. Heat resistant coatings—an overview, in: Coatings for High -temperature Environments: Anti -corrosion and Anti-wear Applications, Cham: Springer Nature Switzerland, 2024: 403-430

    work page 2024

  2. [2]

    Zheng, M

    Q. Zheng, M. Hao, R. Miao, J. Schaadt, Chris Dames. Advances in thermal conductivity for energy applications: a review. Progress in Energy, 2021, 3(1): 12002

    work page 2021

  3. [3]

    Henry, R

    A. Henry, R. Prasher, Arun Majumdar. Five thermal energy grand challenges for decarbonization. Nature Energy, 2020, 5(9): 635-637

    work page 2020

  4. [4]

    S. Wu, T. Yan, Z. Kuai, Weiguo Pan. Thermal conductivity enhancement on phase change materials for thermal energy storage: a review. Energy Storage Materials, 2020, 25: 251-295

    work page 2020

  5. [5]

    Z. He, Y . Yan, Zhien Zhang. Thermal management and temperature uniformity enhancement of electronic devices by micro heat sinks: a review. Energy, 2021, 216: 119223

    work page 2021

  6. [6]

    X. L. Shi, J. Zou, Zhi -Gang Chen. Advanced th ermoelectric design: from materials and structures to devices. Chemical Reviews, 2020, 120(15): 7399-7515

    work page 2020

  7. [7]

    J. Mao, G. Chen, Z. Ren. Thermoelectric cooling materials. Nature Materials, 2021, 20(4): 454-461

    work page 2021

  8. [8]

    J. Pei, B. Cai, H. L. Zhuang, Jing -Feng Li. Bi2Te3-based applied thermoelectric materials: research advances and new challenges. National Science Review, 2020, 7(12): 1856-1858

    work page 2020

  9. [9]

    Z. Zong, S. Deng, Y . Qin, X. Wan, J. Zhan, D. Ma, et al. Enhancing the interfacial thermal conductance of Si/PVDF by strengthening atomic couplings. Nanoscale, 2023, 15(40): 16472-16479

    work page 2023

  10. [10]

    L. Dong, B. Liu, Y . Wang, X. Xu. Tunable thermal conductivity of ferroelectric P(VDF-TrFE) nanofibers via molecular bond modulation. Chinese Physics Letters, 2022, 39(12): 127201

    work page 2022

  11. [11]

    J. Chen, J. He, D. Pan, X. Wang, N. Yang, J. Zhu, et al. Emerging theory and phenomena in thermal conduction: a selective review. Science China Physics, Mechanics & Astronomy, 2022, 65(11): 117002

    work page 2022

  12. [12]

    S. Deng, J . Yuan, Y . Lin, X. Yu, D. Ma, Y . Huang, et al. Electric-field-induced modulation of thermal conductivity in poly(vinylidene fluoride). Nano Energy, 2021, 82: 105749

    work page 2021

  13. [13]

    S. Deng, D. Ma, G. Zhang, N. Yang. Modulating the thermal conductivity of crystalline nylon by tuning hydrogen bonds through structure poling. Journal of Materials Chemistry A, 2021, 9(43): 24472-24479

    work page 2021

  14. [14]

    Davis, M

    Bruce L. Davis, M. I. Hussein. Nanophononic metamaterial: thermal conductivity reduction by local resonance. Physical Review Letters, 2014, 112(5): 55505

    work page 2014

  15. [15]

    Jin-Wu Jiang, J. S. Wang, B. Li. Topological effect on thermal conductivity in graphene. Journal of Applied Physics, 2010, 108(6): 64307

    work page 2010

  16. [16]

    Mitrovic, D

    Jen-Kan Yu, S. Mitrovic, D. Tham, J. Varghese, J. R. Heath. Reduction of thermal conductivity in phononic nanomesh structures. Nature Nanotechnology, 2010, 5(10): 718-721

    work page 2010

  17. [17]

    N. Yang, X. Ni, J. W. Jiang, B. Li. How does folding modulate thermal conductivity of graphene. Applied Physics Letters, 2012, 100(9): 93107

    work page 2012

  18. [18]

    S. Hu, J. Chen, N. Yang, B. Li. Thermal transport in graphene with defect and doping: phonon modes analysis. Carbon, 2017, 116: 139-144

    work page 2017

  19. [19]

    J. Y . Cho, X. Shi, J. R. Salvador, G. P. Meisner, J. Yang, H. Wang, et al. Thermoelectric properties and investigations of low thermal cond uctivity in Ga- doped Cu 2 GeSe 3. Physical Review B, 2011, 84(8): 85207

    work page 2011

  20. [20]

    Shanshan Chen, Q. Wu, C. Mishra, J. Kang, H. Zhang, K. Cho, et al. Thermal conductivity of isotopically modified graphene. Nature Materials, 2012, 11(3): 203-207

    work page 2012

  21. [21]

    Y . Yang, D. Ma, L. Zhang. Introduction of asymmetry to enhance thermal transport in porous metamaterials at low temperature. Chinese Physics Letters, 2023, 40(12)

    work page 2023

  22. [22]

    Bera, Natalio Mingo, S

    C. Bera, Natalio Mingo, S. V olz. Marked effects of alloying on the thermal conductivity of nanoporou s materials. Physical Review Letters, 2010, 104(11): 115502

    work page 2010

  23. [23]

    Chunlei Wan, Z. Qu, A. Du, W. Pan. Order–disorder transition and unconventional thermal conductivities of the (Sm1−x ybx )2 Zr2 O7 series. Journal of the American Ceramic Society, 2011, 94(2): 592-596

    work page 2011

  24. [24]

    M. An, H. Wang, Y . Yuan, D. Chen, W. Ma, S. W. Sharshir, et al. Strong phonon coupling induces low thermal conductivity of one -dimensional carbon boron nanotube. Surfaces and Interfaces, 2022, 28: 101690

    work page 2022

  25. [25]

    Yangjun Qin, L. Mu, X. Wan, Z. Z ong, T. Li, H. Fang, et al. Deep potential for interaction between hydrated cs + and graphene. Langmuir, 2025, 41(18): 11506 - 11514

    work page 2025

  26. [26]

    X. Wan, D. Ma, D. Pan, L. Yang, N. Yang. Optimizing thermal transport in graphene nanoribbon based on phonon resonance hy bridization. Materials Today Physics, 2021, 20: 100445

    work page 2021

  27. [27]

    S. Shen, A. Henry, J. Tong, R. Zheng, Gang Chen. Polyethylene nanofibres with very high thermal conductivities. Nature Nanotechnology, 2010, 5(4): 251-255

    work page 2010

  28. [28]

    D. Pan, T. Li, X. Wan, Z. Zong, Y . Qin, N. Yang. Using targeted phonon excitation to modulate thermal conductivity of boron nitride. Chin. Phys. Lett., 2025

    work page 2025

  29. [29]

    Y . Yoon, Z. Lu, C. Uzundal, R. Qi, W. Zhao, S. Chen, et al. Terahertz phonon engineering with van der waals heterostructures. Natu re, 2024, 631(8022): 771 - 776

    work page 2024

  30. [30]

    X. Wan, Z. Zong, Y . Qin, J. T. Lü, S. V olz, L. Zhang, et al. Modulating thermal conductivity via targeted phonon excitation. Nano Letters, 2024, 24(23): 6889 - 6896

    work page 2024

  31. [31]

    Sekiguchi, H

    F. Sekiguchi, H. Hirori, G. Yumoto, A. Shimazaki, T. Nakamura, A. Wakamiya, et al. Enhancing the hot -phonon bottleneck effect in a metal halide perovskite by terahertz phonon excitation. Physical Review Letters, 2021, 126(7): 77401

    work page 2021

  32. [32]

    Lindsay, D

    L. Lindsay, D. A. Broido, T. L. Reinecke. First -principles determination of ultrahigh thermal conductivity of boron arsenide: a competitor for diamond? Physical Review Letters, 2013, 111(2): 25901

    work page 2013

  33. [33]

    H. Ma, C. Li, S. Tang, J. Yan, A. Alatas, L. Lindsay, et al. Boron arsenide phonon dispersion from inelastic x -ray scattering: potential for ultrahigh thermal conductivity. Physical Review B, 2016, 94(22): 220303

    work page 2016

  34. [34]

    Sang Kang Joon, M. Li, H. Wu, H. Nguyen, Y . Hu. Experimental observation of high thermal conductivity in boron arsenide. Science, 2018, 361(6402): 575-578

    work page 2018

  35. [35]

    Zheng, Y

    Sheng Li, Q. Zheng, Y . Lv, X. Liu, X. Wang, Y . Huang Pinshane, et al. High thermal conductivity in cubic boron arsenide crystals. Science, 2018, 361(6402): 579-581

    work page 2018

  36. [36]

    Kresse, J

    G. Kresse, J. Furthmüller. Efficiency of ab -initio total energy calculations for metals and semiconductors using a plane-wave basis set. Computational Materials Science, 1996, 6(1): 15-50

    work page 1996

  37. [37]

    Kresse, J

    G. Kresse, J. Furthmüller. Efficient iterative schemes for ab initio total -energy calculations using a plane-wave basis set. Physical Review B, 1996, 54(16): 11169

    work page 1996

  38. [38]

    J. A. Perri, S. La Placa, B. Post. New group III-group V compounds: BP and BAs. Acta Crystallographica, 1958, 11(4): 310-310

    work page 1958

  39. [39]

    First-principles phonon calculations with phonopy and Phono3py

    Atsushi Togo. First-principles phonon calculations with phonopy and Phono3py. Journal of the Physical Society of Japan, 2023, 92(1): 12001

    work page 2023

  40. [40]

    W. Li, J. Carrete, N. A. Katcho, N. Mingo. ShengBTE: a solver of the boltzmann transport equation for phonons. Computer Physics Communications, 2014, 185(6): 1747-1758

    work page 2014

  41. [41]

    Z. Han, X. Yang, W. Li, T. Feng, X. Ruan. FourPhonon: an extension module to ShengBTE for computing four-phonon scattering rates and thermal conductivity. Computer Physics Communications, 2022, 270: 108179

    work page 2022

  42. [42]

    T. Feng, X. Ruan. Four-phonon scattering reduces intrinsic thermal conductivity of gr aphene and the contributions from flexural phonons. Physical Review B, 2018, 97(4)

    work page 2018

  43. [43]

    R. G. Greene. Pressure induced metastable amorphization of BAs: evidence for a kinetically frustrated phase transformation. Physical Review Letters, 1994, 73(18): 2476-2479

    work page 1994