pith. sign in

arxiv: 2605.17947 · v1 · pith:IK52GQJXnew · submitted 2026-05-18 · ❄️ cond-mat.mtrl-sci

Indicators for phonon hydrodynamics from first principles predictions of thermal conductivity

Nikhil Malviya , Navaneetha K. Ravichandran This is my paper

Pith reviewed 2026-05-20 09:56 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords phonon hydrodynamicsthermal conductivityPeierls-Boltzmann equationrelaxation time approximationfirst-principles calculationsphonon transportBrillouin zone samplingnon-Fourier heat flow
0
0 comments X

The pith

The ratio of thermal conductivity from the full phonon transport solution to the relaxation time approximation acts as a low-cost indicator for phonon hydrodynamics.

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

The paper establishes that the ratio of thermal conductivity calculated from the complete solution of the linearized Peierls-Boltzmann equation to the value from the relaxation time approximation serves as a practical indicator for when phonon heat flow becomes hydrodynamic. A sympathetic reader would care because this offers an inexpensive computational screen to identify materials that might exhibit collective phonon drift at accessible temperatures, expanding beyond graphite. The work shows that collective drifting of non-equilibrium phonons boosts the ratio while small values align with ordinary diffusive scattering. It further demonstrates that simpler approximations fail to capture the hydrodynamic regime and that denser sampling of the Brillouin zone can lower the ratio in high-conductivity materials at low temperatures.

Core claim

We show that the ratio of thermal conductivity obtained from the complete solution of the linearized Peierls-Boltzmann equation for phonon transport to that from the relaxation time approximation is a low cost indicator for phonon hydrodynamics. Collectively drifting non-equilibrium phonons amplify the ratio of full solution to approximation values, while a small ratio correlates with predominantly diffusive phonon transport. Conventional approaches relying only on normal and Umklapp scattering rates, such as the relaxation time approximation and Callaway methods, are inadequate to predict phonon hydrodynamics. The indicator ratio and therefore the strength of hydrodynamic signatures also be

What carries the argument

The ratio of thermal conductivity from the complete linearized Peierls-Boltzmann equation solution to the relaxation time approximation value, which increases when non-equilibrium phonons drift collectively instead of scattering independently.

If this is right

  • Materials showing a large ratio exhibit hydrodynamic phonon flow driven by collective drift of non-equilibrium phonons.
  • The relaxation time approximation and Callaway approximations cannot reliably predict the hydrodynamic regime.
  • Hydrodynamic signatures weaken with increasing Brillouin zone sampling density in ultrahigh-conductivity materials at low temperatures.
  • Careful convergence studies with respect to Brillouin zone sampling are required for robust predictions of phonon hydrodynamics.

Where Pith is reading between the lines

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

  • The ratio could serve as a quick pre-screen in computational searches for additional materials beyond graphite that support hydrodynamic flow at room temperature.
  • Prior predictions of hydrodynamics in certain ultrahigh-conductivity compounds may require rechecking after finer Brillouin zone sampling.
  • Similar ratio diagnostics might help identify conditions or material classes where other non-Fourier transport effects dominate.

Load-bearing premise

An amplified ratio is caused specifically by collectively drifting non-equilibrium phonons rather than other numerical or physical effects in the transport calculation.

What would settle it

A material computed to have a large ratio but experimentally observed to follow standard diffusive heat flow without collective drift signatures such as unusual boundary scattering or temperature dependence.

Figures

Figures reproduced from arXiv: 2605.17947 by Navaneetha K. Ravichandran, Nikhil Malviya.

Figure 1
Figure 1. Figure 1: (b) that exactly three nearly-degenerate eigen￾modes contribute to κLP BE entirely for a BZ discretiza￾tion of 213 while the spectral contribution of the eigen￾modes to κLP BE is more distributed across three differ￾ent sets of eigenmodes with different eigenvalues for a 353 discretization, reflecting a weaker hydrodynamic sig￾nature in the latter. This conclusion is also strengthened by the smaller eigenv… view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
read the original abstract

Hydrodynamic heat flow, where out-of-equilibrium phonons collectively drift in response to an applied temperature differential, has attracted renewed interest following its experimental observation in graphite at temperatures as high as 300 K. To accelerate discovery of material alternatives to graphite and suitable experimental conditions for realizing this non-Fourier heat flow regime, computationally efficient indicators derived from predictive first principles approaches are necessary. Here we show that the ratio of thermal conductivity ($\kappa$) obtained from the complete solution of the linearized Peierls-Boltzmann equation (LPBE) for phonon transport ($\kappa_{{LPBE}}$), to that from the relaxation time approximation (RTA) for phonon decay ($\kappa_{{RTA}}$), is a low cost indicator for phonon hydrodynamics. We show that collectively drifting non-equilibrium phonons amplify the ratio of $\kappa_{{LPBE}}$ to $\kappa_{{RTA}}$, while a small $\kappa_{{LPBE}}/\kappa_{{RTA}}$ correlates with predominantly diffusive phonon transport. On the other hand, we find that conventional approaches that rely only on momentum-conserving Normal and momentum-dissipating Umklapp scattering rates, such as the RTA and the Callaway approximations to the LPBE, are inadequate to predict phonon hydrodynamics. Furthermore, our study reveals that the indicator ratio - $\kappa_{{LPBE}}/\kappa_{{RTA}}$, and therefore the strength of hydrodynamic signatures, decrease with increasing Brillouin zone (BZ) sampling density for several ultrahigh-$\kappa$ materials at low temperatures, thus underscoring the need for careful BZ sampling convergence studies to ensure robust predictions of phonon hydrodynamics. This computationally inexpensive indicator of phonon hydrodynamics will accelerate the search for new materials that exhibit such unconventional heat flow regimes.

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

2 major / 1 minor

Summary. The manuscript proposes that the ratio of thermal conductivity obtained from the full solution of the linearized Peierls-Boltzmann equation (κ_LPBE) to that from the relaxation-time approximation (κ_RTA) serves as a computationally inexpensive indicator for phonon hydrodynamics. It argues that collective drifting of non-equilibrium phonons amplifies this ratio, while small values correlate with diffusive transport; RTA and Callaway approximations are shown to be inadequate; and the ratio (hence hydrodynamic signatures) decreases with denser Brillouin-zone sampling in ultrahigh-κ materials at low temperature, underscoring the need for convergence checks.

Significance. If the central claim holds after addressing numerical robustness, the ratio would provide a practical first-principles screen for materials that may exhibit hydrodynamic phonon transport, accelerating discovery beyond the graphite examples already observed experimentally. The explicit call for BZ-sampling convergence studies is a constructive contribution to reliable predictions in this field.

major comments (2)
  1. [Abstract] Abstract: the central claim that an elevated κ_LPBE/κ_RTA ratio specifically signals collective phonon drift (as opposed to other effects) is load-bearing, yet the abstract itself reports that the ratio decreases with increasing BZ sampling density for ultrahigh-κ materials at low T. This directly raises the possibility that the reported hydrodynamic indicator is sensitive to k-point discretization of normal-process matrix elements rather than intrinsic momentum-conserving physics; an explicit test isolating the drift contribution (e.g., controlled perturbation of normal scattering while holding Umklapp fixed) is required to substantiate the interpretation.
  2. [Abstract] Abstract: the statement that RTA and Callaway approximations are inadequate to predict phonon hydrodynamics is asserted without quantitative benchmarks (e.g., for which materials and at what temperatures the ratio deviates significantly from unity or from independent hydrodynamic metrics such as Poiseuille flow signatures). Without such data, it is unclear how much the full LPBE solution is strictly necessary versus merely sufficient for the proposed indicator.
minor comments (1)
  1. [Abstract] The abstract would benefit from a brief definition or reference to the precise form of the RTA used (e.g., whether it includes only Umklapp or also normal processes) to avoid ambiguity in the ratio definition.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. We address each major comment below and have revised the manuscript to improve clarity and provide additional supporting details where feasible.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the central claim that an elevated κ_LPBE/κ_RTA ratio specifically signals collective phonon drift (as opposed to other effects) is load-bearing, yet the abstract itself reports that the ratio decreases with increasing BZ sampling density for ultrahigh-κ materials at low T. This directly raises the possibility that the reported hydrodynamic indicator is sensitive to k-point discretization of normal-process matrix elements rather than intrinsic momentum-conserving physics; an explicit test isolating the drift contribution (e.g., controlled perturbation of normal scattering while holding Umklapp fixed) is required to substantiate the interpretation.

    Authors: We agree that convergence with respect to Brillouin-zone sampling is essential and that the observed decrease in the ratio with denser sampling warrants careful discussion. This trend, which we already report, reflects the improved resolution of momentum-conserving normal processes on finer grids, allowing a more accurate representation of collective phonon drift in the full LPBE solution. In the converged limit the ratio remains distinctly larger than unity for materials known to support hydrodynamic transport. While a controlled perturbation that artificially modifies normal scattering rates while holding Umklapp rates fixed would be a valuable numerical experiment, it lies outside the scope of the present first-principles framework and would require substantial additional methodological development. We have therefore revised the abstract and added a concise paragraph in the main text that explicitly links the sampling dependence to the drift term in the LPBE, thereby strengthening the physical interpretation without claiming an exhaustive isolation test. revision: partial

  2. Referee: [Abstract] Abstract: the statement that RTA and Callaway approximations are inadequate to predict phonon hydrodynamics is asserted without quantitative benchmarks (e.g., for which materials and at what temperatures the ratio deviates significantly from unity or from independent hydrodynamic metrics such as Poiseuille flow signatures). Without such data, it is unclear how much the full LPBE solution is strictly necessary versus merely sufficient for the proposed indicator.

    Authors: We acknowledge that the original abstract was concise and did not include explicit numerical benchmarks. The full manuscript already demonstrates, across multiple ultrahigh-κ materials at low temperatures, that the κ_LPBE/κ_RTA ratio exceeds 2 in several cases while the RTA and Callaway approximations remain close to unity. To make this quantitative evidence immediately visible, we have added a new table summarizing the ratios for representative materials and temperatures and have expanded the discussion to note the absence of Poiseuille-flow signatures in the approximate solutions. These revisions clarify that the full LPBE solution is required to capture the hydrodynamic enhancement captured by the proposed indicator. revision: yes

Circularity Check

0 steps flagged

No significant circularity; ratio derived from independent LPBE and RTA solutions

full rationale

The paper defines the indicator as the ratio of thermal conductivity from the complete linearized Peierls-Boltzmann equation solution (κ_LPBE) to the relaxation-time approximation (κ_RTA). These are two standard, independently formulated methods for solving the phonon Boltzmann transport equation; the ratio is obtained by direct numerical evaluation on first-principles phonon dispersions and scattering rates rather than by fitting a parameter or re-expressing one quantity in terms of the other. The physical interpretation linking large ratios to collective drift is presented as an observed outcome of the calculations, not as a definitional equivalence or self-referential construction. No load-bearing step reduces to a self-citation chain, an ansatz smuggled via prior work, or a fitted input renamed as a prediction. The result remains self-contained against external benchmarks of the two transport approximations.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard phonon transport theory without introducing new free parameters or postulated entities; the indicator is simply a ratio of two existing computational approaches.

axioms (2)
  • domain assumption The linearized Peierls-Boltzmann equation provides an accurate reference solution for phonon transport in the materials considered.
    Invoked when the full LPBE solution is taken as the benchmark against which the RTA ratio is evaluated.
  • domain assumption The relaxation time approximation supplies a reliable baseline for predominantly diffusive phonon transport.
    Used to define the denominator of the proposed indicator ratio.

pith-pipeline@v0.9.0 · 5862 in / 1513 out tokens · 57477 ms · 2026-05-20T09:56:47.976767+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

37 extracted references · 37 canonical work pages · 1 internal anchor

  1. [1]

    1 is solved forf λ using iterative solvers [20, 21], which have become computa- tionally inexpensive in recent years due to the improve- ments in computing hardware and algorithms

    OnceΩis computed, Eq. 1 is solved forf λ using iterative solvers [20, 21], which have become computa- tionally inexpensive in recent years due to the improve- ments in computing hardware and algorithms. Fromf λ, we obtain the heat fluxJ= 1 V P λ ℏωλvλfλ for a crys- tal volumeV, and cast it into a linear response relation: J=−κ LP BE∇Tto obtainκ LP BE. In ...

    work page 2022

  2. [2]

    R. A. Guyer and J. A. Krumhansl, Thermal Conductiv- ity, Second Sound, and Phonon Hydrodynamic Phenom- ena in Nonmetallic Crystals, Physical Review148, 778 (1966)

    work page 1966

  3. [3]

    R. J. Hardy, Phonon Boltzmann Equation and Second Sound in Solids, Physical Review B2, 1193 (1970)

    work page 1970

  4. [4]

    S. Lee, D. Broido, K. Esfarjani, and G. Chen, Hydrody- namic phonon transport in suspended graphene, Nature Communications6, 6290 (2015)

    work page 2015

  5. [5]

    Cepellotti, G

    A. Cepellotti, G. Fugallo, L. Paulatto, M. Lazzeri, F. Mauri, and N. Marzari, Phonon hydrodynamics in two-dimensional materials, Nature Communications6, 6400 (2015)

    work page 2015

  6. [6]

    Simoncelli, N

    M. Simoncelli, N. Marzari, and A. Cepellotti, General- ization of Fourier’s Law into Viscous Heat Equations, Physical Review X10, 011019 (2020)

    work page 2020

  7. [7]

    Huang, R

    X. Huang, R. Anufriev, L. Jalabert, K. Watanabe, T. Taniguchi, Y. Guo, Y. Ni, S. Volz, and M. Nomura, A graphite thermal Tesla valve driven by hydrodynamic phonon transport, Nature634, 1086 (2024)

    work page 2024

  8. [8]

    Shi, Nonresistive heat transport by collective phonon flow, Science364, 332 (2019)

    L. Shi, Nonresistive heat transport by collective phonon flow, Science364, 332 (2019)

    work page 2019

  9. [9]

    Chen, Non-Fourier phonon heat conduction at the mi- croscale and nanoscale, Nature Reviews Physics3, 555 (2021)

    G. Chen, Non-Fourier phonon heat conduction at the mi- croscale and nanoscale, Nature Reviews Physics3, 555 (2021)

    work page 2021

  10. [10]

    H. E. Jackson, C. T. Walker, and T. F. McNelly, Second Sound in NaF, Physical Review Letters25, 26 (1970)

    work page 1970

  11. [11]

    H. E. Jackson and C. T. Walker, Thermal Conductivity, Second Sound, and Phonon-Phonon Interactions in NaF, Physical Review B3, 1428 (1971)

    work page 1971

  12. [12]

    Narayanamurti and R

    V. Narayanamurti and R. C. Dynes, Observation of sec- ond sound in bismuth, Physical Review Letters28, 1461 (1972)

    work page 1972

  13. [13]

    Kawabata, K

    T. Kawabata, K. Shimura, Y. Ishii, M. Koike, K. Yoshida, S. Yonehara, K. Yokoi, A. Subedi, K. Behnia, and Y. Machida, Phonon hydrodynamic regimes in sap- phire, Physical Review Research7, 033017 (2025)

    work page 2025

  14. [14]

    Beardo, M

    A. Beardo, M. L´ opez-Su´ arez, L. A. P´ erez, L. Sendra, M. I. Alonso, C. Melis, J. Bafaluy, J. Camacho, L. Colombo, R. Rurali, F. X. Alvarez, and J. S. Reparaz, Observation of second sound in a rapidly varying temper- ature field in Ge, Science Advances7, eabg4677 (2021)

    work page 2021

  15. [15]

    Huberman, R

    S. Huberman, R. A. Duncan, K. Chen, B. Song, V. Chiloyan, Z. Ding, A. A. Maznev, G. Chen, and K. A. Nelson, Observation of second sound in graphite at tem- peratures above 100 K, Science364, 375 (2019)

    work page 2019

  16. [16]

    Z. Ding, K. Chen, B. Song, J. Shin, A. A. Maznev, K. A. Nelson, and G. Chen, Observation of second sound in graphite over 200 K, Nature Communications13, 285 (2022)

    work page 2022

  17. [17]

    Z. Xie, Y. Zhang, X. Huang, Z. Ding, J. Wei, D. Dong, K. Cao, T. Lai, K. Watanabe, T. Taniguchi, X. Qian, M. Nomura, and K. Chen, Room-temperature second sound in isotopically pure graphite, Nature Communi- 6 cations 10.1038/s41467-026-70807-3 (2026)

    work page doi:10.1038/s41467-026-70807-3 2026

  18. [18]

    N. K. Ravichandran and D. Broido, Unified first- principles theory of thermal properties of insulators, Physical Review B98, 085205 (2018)

    work page 2018

  19. [19]

    N. K. Ravichandran and D. Broido, Phonon-phonon in- teractions in strongly bonded solids: selection rules and higher-order processes, Physical Review X10, 021063 (2020)

    work page 2020

  20. [20]

    N. K. Ravichandran, Elasticity reshapes heat flow in graphene (2026), arXiv:2604.03910 [cond-mat]

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  21. [21]

    D. A. Broido, M. Malorny, G. Birner, N. Mingo, and D. A. Stewart, Intrinsic lattice thermal conductivity of semiconductors from first principles, Applied Physics Letters91, 231922 (2007)

    work page 2007

  22. [22]

    Fugallo, M

    G. Fugallo, M. Lazzeri, L. Paulatto, and F. Mauri, Ab initio variational approach for evaluating lattice thermal conductivity, Physical Review B88, 045430 (2013)

    work page 2013

  23. [23]

    Malviya and N

    N. Malviya and N. K. Ravichandran, Failure of the Call- away description of heat flow in boron arsenide and boron antimonide driven by phonon scattering selection rules, Physical Review B108, 155201 (2023)

    work page 2023

  24. [24]

    Callaway, Model for Lattice Thermal Conductivity at Low Temperatures, Physical Review113, 1046 (1959)

    J. Callaway, Model for Lattice Thermal Conductivity at Low Temperatures, Physical Review113, 1046 (1959)

    work page 1959

  25. [25]

    P. B. Allen, Improved Callaway model for lattice thermal conductivity, Physical Review B88, 144302 (2013)

    work page 2013

  26. [26]

    Guo and M

    Y. Guo and M. Wang, Heat transport in two-dimensional materials by directly solving the phonon Boltzmann equation under Callaway’s dual relaxation model, Phys- ical Review B96, 134312 (2017)

    work page 2017

  27. [27]

    Y. Zhou, X. Zhang, and M. Hu, Nonmonotonic Diam- eter Dependence of Thermal Conductivity of Extremely Thin Si Nanowires: Competition between Hydrodynamic Phonon Flow and Boundary Scattering, Nano Letters17, 1269 (2017)

    work page 2017

  28. [28]

    Ezzahri, Thermal decay in the one-dimensional tran- sient thermal grating experiment using modified Debye- Callaway model, Physical Review B106, 125203 (2022)

    Y. Ezzahri, Thermal decay in the one-dimensional tran- sient thermal grating experiment using modified Debye- Callaway model, Physical Review B106, 125203 (2022)

    work page 2022

  29. [29]

    X. Qian, C. Zhang, T.-H. Liu, and R. Yang, Analyti- cal Green’s function of the multidimensional Boltzmann transport equation for modeling hydrodynamic second sound, Physical Review B111, 035406 (2025)

    work page 2025

  30. [30]

    Malviya and N

    N. Malviya and N. K. Ravichandran, Callaway approxi- mation to the Boltzmann equation fails to predict phonon hydrodynamics, Journal of Physics: Condensed Matter 37, 265701 (2025)

    work page 2025

  31. [31]

    Cepellotti and N

    A. Cepellotti and N. Marzari, Thermal Transport in Crystals as a Kinetic Theory of Relaxons, Physical Re- view X6, 041013 (2016)

    work page 2016

  32. [32]

    Malviya and N

    N. Malviya and N. K. Ravichandran, Efficient calcula- tion of phonon dynamics through a low-rank solution of the Boltzmann equation (2025), arXiv:2502.00337 [cond- mat]

    work page arXiv 2025

  33. [33]

    Hua and L

    C. Hua and L. Lindsay, Space-time dependent thermal conductivity in nonlocal thermal transport, Physical Re- view B102, 104310 (2020)

    work page 2020

  34. [34]

    L. P. Pitaevskii and E. M. Lifshitz,Physical kinetics: vol- ume 10(Butterworth-Heinemann, 2012)

    work page 2012

  35. [35]

    Y. Zhu, A. L. N. Kondusamy, K. Chen, P. Koirala, H. Wu, M. S. Patel, E. R. Glaser, S. White, J. C. Culbertson, J. L. Lyons, N. K. Ravichandran, M. Ze- tati, X. Wang, R. Addou, R. M. Wallace, S. Hou, S. Vaziri, X. Bao, D. Broido, G. Chen, and B. Lv, Vapor- Flux Growth of c-BP Single Crystals With Concurrently High Electrical Resistivity and Isotope-Enhance...

    work page 2026

  36. [36]

    R. B. Lehoucq, D. C. Sorensen, and C. Yang,ARPACK Users’ Guide, Software, Environments, and Tools (Soci- ety for Industrial and Applied Mathematics, 1998)

    work page 1998

  37. [37]

    N. K. Ravichandran, Low bending rigidity and large Young’s modulus drive strong flexural phonon renormalization in two-dimensional monolayers (2026), arXiv:2603.16717 [cond-mat]

    work page arXiv 2026