pith. sign in

arxiv: 2507.05682 · v2 · submitted 2025-07-08 · ❄️ cond-mat.mtrl-sci · cond-mat.mes-hall

Unified Statistical Theory of Heat Conduction in Nonuniform Media

Yi Zeng , Jianjun Dong This is my paper

Pith reviewed 2026-05-19 06:30 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci cond-mat.mes-hall
keywords heat conductionnonlocal transportspatiotemporal kernelprojection operatorGreen-Kubophonon transportinterfacial resistancenonuniform media
0
0 comments X

The pith

A single causal kernel derived from heat-flux correlations unifies conduction across diffusive, ballistic, and hydrodynamic regimes in nonuniform media.

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

The paper derives a causal two-point spatiotemporal kernel for heat conduction from the Zwanzig projection-operator formalism applied to equilibrium correlations. This kernel is defined microscopically as the space-resolved heat-flux time-correlation function and treats temporal memory, spatial nonlocality, and material heterogeneity on equal footing. Classical diffusion, nonlocal transport, and hydrodynamic descriptions appear as controlled limits of the same kernel. The framework incorporates interfaces by making Kapitza resistance a coarse-grained limit and connects directly to atomistic simulations without extra fitting parameters.

Core claim

Using the Zwanzig projection-operator formalism, we derive a causal two-point spatiotemporal kernel for heat conduction, defined microscopically as a space-resolved equilibrium heat-flux time-correlation function, that encodes temporal memory, spatial nonlocality, and material heterogeneity on equal footing. Classical diffusion, nonlocal transport, and hydrodynamic models emerge as controlled asymptotic limits of this kernel, providing a unified constitutive description across diffusive, quasi-ballistic, and hydrodynamic regimes. Interfacial heat transfer is incorporated through a spatially resolved kernel formulation, in which the conventional Kapitza resistance arises as a coarse-grained l

What carries the argument

The causal two-point spatiotemporal kernel, defined microscopically as the space-resolved equilibrium heat-flux time-correlation function from the Zwanzig projection-operator formalism, which encodes memory, nonlocality, and heterogeneity together.

Load-bearing premise

The Zwanzig projection-operator formalism applied to heat-flux time correlations in nonuniform media yields a causal kernel whose controlled asymptotic limits recover classical diffusion, nonlocal transport, and hydrodynamic models without additional ad-hoc closures.

What would settle it

If molecular dynamics calculations of the space-resolved heat-flux correlation function for silicon in a transient thermal grating setup produce a kernel that, when inserted into the transport equation, fails to match measured temperature decay profiles across varying grating periods, the claim of unification through this kernel would be refuted.

read the original abstract

Using the Zwanzig projection-operator formalism, we derive a causal two-point spatiotemporal kernel for heat conduction, defined microscopically as a space-resolved equilibrium heat-flux time-correlation function, that encodes temporal memory, spatial nonlocality, and material heterogeneity on equal footing. Classical diffusion, nonlocal transport, and hydrodynamic models emerge as controlled asymptotic limits of this kernel, providing a unified constitutive description across diffusive, quasi-ballistic, and hydrodynamic regimes. Interfacial heat transfer is incorporated through a spatially resolved kernel formulation, in which the conventional Kapitza resistance arises as a coarse-grained limit. The kernel admits a spatiotemporal Green--Kubo representation and can, in principle, be evaluated from atomistic simulations for bulk media, providing a direct connection between microscopic dynamics and continuum transport without empirical closure. For crystalline solids, we derive explicit kernel forms in the hydrodynamic and attenuated-streaming limits and introduce a hybrid reduction that captures the coexistence of collective and quasi-ballistic transport. For disordered harmonic solids, the framework recovers a spatial diffusion kernel consistent with the Allen--Feldman limit. To illustrate the theory, we construct the kernel for silicon at room temperature within the relaxation-time approximation and apply it to transient thermal grating configurations. Spatial nonlocality associated with the phonon mean-free-path distribution is the primary source of deviation from Fourier transport under these conditions, while temporal memory mainly influences short-time dynamics. These findings identify the spatiotemporal kernel as a unifying constitutive descriptor whose coarse-grained limits recover conventional transport coefficients.

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 / 2 minor

Summary. The manuscript develops a unified statistical theory of heat conduction in nonuniform media by applying the Zwanzig projection-operator formalism to derive a causal two-point spatiotemporal kernel, defined microscopically as a space-resolved equilibrium heat-flux time-correlation function. This kernel is asserted to encode temporal memory, spatial nonlocality, and material heterogeneity on equal footing, with classical diffusion, nonlocal transport, and hydrodynamic models emerging as controlled asymptotic limits. Explicit kernel forms are derived for crystalline solids in hydrodynamic and attenuated-streaming limits, a hybrid reduction for collective and quasi-ballistic transport, and recovery of the Allen-Feldman spatial diffusion kernel for disordered harmonic solids. The framework is illustrated for silicon at room temperature within the relaxation-time approximation applied to transient thermal grating configurations, identifying spatial nonlocality from the phonon mean-free-path distribution as the main deviation from Fourier transport, while interfacial heat transfer is incorporated via a spatially resolved formulation in which Kapitza resistance appears as a coarse-grained limit. The kernel admits a spatiotemporal Green-Kubo representation evaluable from atomistic simulations.

Significance. If the central derivation is rigorous and the asymptotic limits are shown to be controlled without additional closures, the work would provide a valuable microscopic-to-continuum bridge for heat transport in heterogeneous and nanostructured materials, enabling direct computation of transport kernels from simulations and a parameter-free unification across diffusive, quasi-ballistic, and hydrodynamic regimes.

major comments (2)
  1. [Abstract and theory derivation] The skeptic concern regarding unstated choices in the slow-variable subspace for the Zwanzig projection in nonuniform media is not resolved in the provided derivation outline; the abstract and summary sections assert a unique causal kernel from space-resolved heat-flux correlations, but without an explicit definition of the position-dependent projection operator (e.g., in the methods or theory section), it is unclear whether implicit scale separations are avoided, which directly affects the claim of controlled asymptotic limits without ad-hoc closures.
  2. [Silicon illustration section] § on silicon example: the statement that spatial nonlocality is the primary source of deviation from Fourier transport in transient thermal grating relies on the relaxation-time approximation kernel; quantitative comparison to full molecular dynamics or error estimates on the mean-free-path distribution effects are needed to substantiate this as load-bearing for the unified claim.
minor comments (2)
  1. The notation distinguishing the microscopic kernel from its coarse-grained limits (e.g., Kapitza resistance) could be made more explicit to aid readability.
  2. A brief discussion of how the Green-Kubo representation connects to existing simulation protocols would strengthen the reproducibility claim.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and insightful comments. We provide point-by-point responses to the major comments and indicate the changes we will implement in the revised manuscript.

read point-by-point responses

  1. Referee: [Abstract and theory derivation] The skeptic concern regarding unstated choices in the slow-variable subspace for the Zwanzig projection in nonuniform media is not resolved in the provided derivation outline; the abstract and summary sections assert a unique causal kernel from space-resolved heat-flux correlations, but without an explicit definition of the position-dependent projection operator (e.g., in the methods or theory section), it is unclear whether implicit scale separations are avoided, which directly affects the claim of controlled asymptotic limits without ad-hoc closures.

    Authors: We thank the referee for raising this important issue concerning the explicit definition of the projection operator. In our derivation, the slow-variable subspace consists of the local energy density at each position, and the Zwanzig projector is defined accordingly to be spatially resolved. This construction ensures that the resulting kernel is causal and that the asymptotic limits are controlled by the separation of timescales and lengthscales inherent to the physical regimes, without introducing additional closures. To resolve any ambiguity, we will expand the theory section to include the explicit mathematical form of the position-dependent projection operator and a step-by-step outline of the derivation. This addition will directly address the concern about unstated choices and reinforce the claim of a unified framework without ad-hoc assumptions. revision: yes

  2. Referee: [Silicon illustration section] § on silicon example: the statement that spatial nonlocality is the primary source of deviation from Fourier transport in transient thermal grating relies on the relaxation-time approximation kernel; quantitative comparison to full molecular dynamics or error estimates on the mean-free-path distribution effects are needed to substantiate this as load-bearing for the unified claim.

    Authors: We agree with the referee that additional substantiation would strengthen the silicon illustration. The current analysis employs the relaxation-time approximation to construct the kernel from known phonon properties. While performing full molecular dynamics simulations for the transient thermal grating setup would be valuable, it lies beyond the scope of the present theoretical development. In the revised manuscript, we will include error estimates by examining the sensitivity of the nonlocality effects to uncertainties in the mean-free-path distribution, using literature values for phonon lifetimes in silicon. We will also add a paragraph discussing the validity of the relaxation-time approximation in this context and how it supports the identification of spatial nonlocality as the dominant deviation from Fourier's law at the considered length and time scales. revision: partial

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained from microscopic correlations

full rationale

The paper applies the standard Zwanzig projection-operator formalism to define a causal kernel directly as the space-resolved equilibrium heat-flux time-correlation function. This is a first-principles microscopic definition rather than a fit, renaming, or self-referential construction. Asymptotic limits (diffusive, nonlocal, hydrodynamic) are recovered as controlled limits of this kernel without additional closures or ad-hoc parameters. No load-bearing self-citations, uniqueness theorems from prior author work, or smuggled ansatzes are indicated in the abstract or description. The connection to atomistic simulations and recovery of known limits (e.g., Allen-Feldman) provides independent content. The derivation chain does not reduce the claimed result to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

Ledger constructed from abstract alone; full text may list additional assumptions or parameters.

axioms (1)
  • domain assumption Zwanzig projection-operator formalism is applicable to deriving a causal heat-flux kernel in nonuniform media
    Invoked at the start of the derivation to obtain the two-point spatiotemporal kernel.
invented entities (1)
  • spatiotemporal heat-conduction kernel no independent evidence
    purpose: Encodes temporal memory, spatial nonlocality, and material heterogeneity on equal footing
    Introduced as the central derived object; abstract provides no independent falsifiable test outside the derivation itself.

pith-pipeline@v0.9.0 · 5795 in / 1331 out tokens · 63227 ms · 2026-05-19T06:30:43.874145+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

62 extracted references · 62 canonical work pages

  1. [1]

    As a result, all phonon properties, includ- ing mode frequency ω, group velocity ⃗ v, and the phonon scattering matrix L, are spatially homogeneous and independent of position ⃗ r

    Homogeneous Phonon Gases We begin by considering an infinitely large and defect-free crystal in which phonon–phonon scattering arises solely from lattice anharmonicity. As a result, all phonon properties, includ- ing mode frequency ω, group velocity ⃗ v, and the phonon scattering matrix L, are spatially homogeneous and independent of position ⃗ r. Under t...

    work page

  2. [2]

    Inhomogeneous Phonon Gases Real materials consist of finite-sized lattices with surfaces and may exhibit spatial inho- mogeneities in atomic structure and composition. When these variations occur gradually in space, without sharp discontinuities, the conventional phonon gas model can be extended to an inhomogeneous phonon framework in which material prope...

    work page

  3. [3]

    A detailed mathematical derivation of this expansion is provided in Appendix B 4

    Hierarchy of Nonlocality Equation (28) reveals an intrinsic hierarchical structure in the mode-resolved dynamics of ηλ(⃗ r, t), which admits an iterative expansion. A detailed mathematical derivation of this expansion is provided in Appendix B 4. The zeroth-order solution, η(0) λ (⃗ r), is given by η(0) λ (⃗ r) = − γ−1 λ (⃗ r) kBT 2 0 ⃗Λ(⃗ r; λ) · ⃗∇rT (⃗...

    work page

  4. [4]

    Kernel Function Representation While Eq. (33) provides a general and complete mathematical description of heat con- duction in the diffusive limit by using macroscopically measurable generalized transport coefficients, its infinite series form is often cumbersome to implement numerically. To ad- dress this, we reformulate the theory using a more compact b...

    work page

  5. [5]

    Interfacial Heat Conduction An important feature of the spatial kernel function formalism presented in Eq. (35) is its universality for nonuniform media including both inhomogeneous materials with varying compositions and finite dimensions heterogeneous materials featuring sharply defined inter- faces, where temperature gradients may become discontinuous ...

    work page

  6. [6]

    The time evolution off(Γ, t) follows the Liouville equation: ∂f (Γ, t) ∂t = −iˆL(Γ)f(Γ, t), (A1) where ˆL is the Liouville operator

    Zwanzig formulism Let f(Γ, t) be the full probability distribution in the complete phase space Γ, which includes all microscopic degrees of freedom. The time evolution off(Γ, t) follows the Liouville equation: ∂f (Γ, t) ∂t = −iˆL(Γ)f(Γ, t), (A1) where ˆL is the Liouville operator. The formal solution is f(Γ, t) = e−itˆLf(Γ, 0) for a given initial conditio...

    work page

  7. [7]

    In heat transport, the local energy density e(⃗ r, t) and the local heat flux ⃗j(⃗ r, t) are fundamental and satisfy: ∂ ∂t e(⃗ r, t) + ⃗∇r · ⃗j(⃗ r, t) = 0

    Relevant Variables at Local Thermal Equilibrium A key step is the identification of relevant variables that describe irreversible processes. In heat transport, the local energy density e(⃗ r, t) and the local heat flux ⃗j(⃗ r, t) are fundamental and satisfy: ∂ ∂t e(⃗ r, t) + ⃗∇r · ⃗j(⃗ r, t) = 0. (A6) However, these fields are not directly accessible in e...

    work page

  8. [8]

    Local Homogeneous Approximation Many real materials exhibit spatial variations in structures and compositions. When the spatial resolution of measured fields T (⃗ r, t) and ⃗j(⃗ r, t) is coarser than the relevant inhomo- geneity length scales, a mean-field or coarse-grained approximation can be adopted. In this case, the medium is treated as piecewise wit...

    work page

  9. [9]

    Thermal Equilibrium At thermal equilibrium at temperature T0, the average occupation number of each phonon mode α follows the Bose–Einstein distribution: neqα(T0) = 1 eℏωα/kBT0 − 1 , (B1) where ℏ is the reduced Planck constant, kB is the Boltzmann constant, and ωα is the angular frequency of mode α. The equilibrium thermal energy density at temperature T0...

    work page

  10. [10]

    (9), can be recast can be recast in the Zwanzig formalism by expanding the phonon distribution in a complete basis of eigenmodes, as introduced in Eqs

    Solution of Eigenmode-Resolved phBTE The phBTE, given in Eq. (9), can be recast can be recast in the Zwanzig formalism by expanding the phonon distribution in a complete basis of eigenmodes, as introduced in Eqs. (5) and (6). This approach describes the nonequilibrium phonon dynamics in terms of the local T (⃗ r, t) and two sets of auxiliary variables: th...

    work page

  11. [11]

    Position-Dependent Formalism For spatially inhomogeneous materials, the Zwanzig variables acquire spatial dependence. The coupled equations read C · ∂T (⃗ r, t) ∂t = − MOX λ=1 ⃗Λ(⃗ r; λ) · ⃗∇rηλ(⃗ r, t), (B12a) ∂ηλ(⃗ r, t) ∂t + γλ(⃗ r)ηλ(⃗ r, t) = − 1 kBT 2 0 ⃗Λ(⃗ r; λ) · ⃗∇rT (⃗ r, t) − MEX ι=1 ⃗Π(⃗ r; λ, ι) · ⃗∇rθι(⃗ r, t),(B12b) ∂θι(⃗ r, t) ∂t + ει(⃗ r...

    work page

  12. [12]

    (15) allows systematic inclusion of memory and spatial nonlocality

    Iterative Expansion of Modal Solutions The iterative expansion of Eq. (15) allows systematic inclusion of memory and spatial nonlocality. For example, the flux-mode variable ηλ(⃗ r, t) can be written as: ηλ = η(0) λ + ∆η(1) λ + ∆η(2) λ + . . . , (B15) 34 where the zeroth-order term η(0) λ is η(0) λ (⃗ r, t) = − 1 kBT 2 0 Z ∞ 0 dt1e−γλ(⃗ r)t1 ⃗Λ(⃗ r; λ) · ...

    work page

  13. [13]

    Conduction Kernel Formalism To capture all orders, an infinite set of time-domain kernels is required. For practical computation, a more compact representation is achieved by introducing a spatiotemporal kernel ← →Z (⃗ r,⃗R, t): ⃗j(⃗ r, t) = − Z V d3 ⃗R Z ∞ 0 dt′ ← →Z (⃗ r,⃗R, t′) · ⃗∇rT ( ⃗R, t − t′). (B27) Here, we adopt the Taylor expansion series: ⃗∇r...

    work page

  14. [14]

    J. B. J. Fourier, Th´ eorie analytique de la chaleur, Vol. 1 (Gauthier-Villars, 1888)

    work page

  15. [15]

    Carslaw and J

    H. Carslaw and J. Jaeger, Conduction of Heat in Solids , 2nd ed. (Oxford University Press, Oxford, 1959)

    work page 1959

  16. [16]

    M. S. Green, Markoff random processes and the statistical mechanics of time-dependent phe- nomena. ii. irreversible processes in fluids, J. Chem. Phys. 22, 398 (1954)

    work page 1954

  17. [17]

    J. Dong, O. F. Sankey, and C. W. Myles, Theoretical study of the lattice thermal conductivity in ge framework semiconductors, Phys. Rev. Lett. 86, 2361 (2001)

    work page 2001

  18. [18]

    Kapitza, The study of heat transfer in helium ii, J

    P. Kapitza, The study of heat transfer in helium ii, J. Phys.(Moscow) 4, 181 (1941)

    work page 1941

  19. [19]

    C. C. Ackerman, B. Bertman, H. A. Fairbank, and R. Guyer, Second sound in solid helium, Physical Review Letters 16, 789 (1966)

    work page 1966

  20. [20]

    McNelly, S

    T. McNelly, S. Rogers, D. Channin, R. Rollefson, W. Goubau, G. Schmidt, J. Krumhansl, and R. Pohl, Heat pulses in naf: onset of second sound, Physical Review Letters 24, 100 (1970)

    work page 1970

  21. [21]

    H. E. Jackson, C. T. Walker, and T. F. McNelly, Second sound in naf, Physical Review Letters 25, 26 (1970)

    work page 1970

  22. [22]

    D. W. Pohl and V. Irniger, Observation of second sound in naf by means of light scattering, Physical Review Letters 36, 480 (1976)

    work page 1976

  23. [23]

    Narayanamurti and R

    V. Narayanamurti and R. Dynes, Observation of second sound in bismuth, Physical Review Letters 28, 1461 (1972)

    work page 1972

  24. [24]

    Danil’Chenko, V

    B. Danil’Chenko, V. Poroshin, and O. Sarbei, An observation of second sound in sapphire, JETP Lett 30 (1979)

    work page 1979

  25. [25]

    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 temperatures above 100 k, Science 364, 375 (2019)

    work page 2019

  26. [26]

    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 communications 13, 285 (2022)

    work page 2022

  27. [27]

    Jeong, X

    J. Jeong, X. Li, S. Lee, L. Shi, and Y. Wang, Transient hydrodynamic lattice cooling by 38 picosecond laser irradiation of graphite, Physical Review Letters 127, 085901 (2021)

    work page 2021

  28. [28]

    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, et al., Observation of second sound in a rapidly varying temperature field in ge, Science advances 7, eabg4677 (2021)

    work page 2021

  29. [29]

    A. A. Maznev, J. A. Johnson, and K. A. Nelson, Onset of nondiffusive phonon transport in transient thermal grating decay, Physical Review B—Condensed Matter and Materials Physics 84, 195206 (2011)

    work page 2011

  30. [30]

    A. J. Minnich, J. A. Johnson, A. J. Schmidt, K. Esfarjani, M. S. Dresselhaus, K. A. Nelson, and G. Chen, Thermal conductivity spectroscopy technique to measure phonon mean free paths, Physical review letters 107, 095901 (2011)

    work page 2011

  31. [31]

    J. A. Johnson, A. Maznev, J. Cuffe, J. K. Eliason, A. J. Minnich, T. Kehoe, C. M. S. Torres, G. Chen, and K. A. Nelson, Direct measurement of room-temperature nondiffusive thermal transport over micron distances in a silicon membrane, Physical Review Letters 110, 025901 (2013)

    work page 2013

  32. [32]

    J. A. Johnson, J. K. Eliason, A. A. Maznev, T. Luo, and K. A. Nelson, Non-diffusive thermal transport in gaas at micron length scales, Journal of Applied Physics 118 (2015)

    work page 2015

  33. [33]

    C. Hua, L. Lindsay, X. Chen, and A. J. Minnich, Generalized fourier’s law for nondiffusive thermal transport: Theory and experiment, Physical Review B 100, 085203 (2019)

    work page 2019

  34. [34]

    Beardo, J

    A. Beardo, J. L. Knobloch, L. Sendra, J. Bafaluy, T. D. Frazer, W. Chao, J. N. Hernandez- Charpak, H. C. Kapteyn, B. Abad, M. M. Murnane, et al., A general and predictive under- standing of thermal transport from 1d-and 2d-confined nanostructures: theory and experi- ment, ACS nano 15, 13019 (2021)

    work page 2021

  35. [35]

    Machida, A

    Y. Machida, A. Subedi, K. Akiba, A. Miyake, M. Tokunaga, Y. Akahama, K. Izawa, and K. Behnia, Observation of poiseuille flow of phonons in black phosphorus, Science advances 4, eaat3374 (2018)

    work page 2018

  36. [36]

    Markov, J

    M. Markov, J. Sjakste, G. Barbarino, G. Fugallo, L. Paulatto, M. Lazzeri, F. Mauri, and N. Vast, Hydrodynamic heat transport regime in bismuth: a theoretical viewpoint, Physical Review Letters 120, 075901 (2018)

    work page 2018

  37. [37]

    Martelli, J

    V. Martelli, J. L. Jim´ enez, M. Continentino, E. Baggio-Saitovitch, and K. Behnia, Ther- mal transport and phonon hydrodynamics in strontium titanate, Physical review letters 120, 125901 (2018). 39

    work page 2018

  38. [38]

    Machida, N

    Y. Machida, N. Matsumoto, T. Isono, and K. Behnia, Phonon hydrodynamics and ultrahigh– room-temperature thermal conductivity in thin graphite, Science 367, 309 (2020)

    work page 2020

  39. [39]

    Huang, Y

    X. Huang, Y. Guo, Y. Wu, S. Masubuchi, K. Watanabe, T. Taniguchi, Z. Zhang, S. Volz, T. Machida, and M. Nomura, Observation of phonon poiseuille flow in isotopically purified graphite ribbons, Nature Communications 14, 2044 (2023)

    work page 2044

  40. [40]

    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, Nature , 1 (2024)

    work page 2024

  41. [41]

    J. C. Maxwell, On the dynamical theory of gases, in The kinetic theory of gases: an anthology of classic papers with historical commentary (World Scientific, 2003) pp. 197–261

    work page 2003

  42. [42]

    Cattaneo, Atti semin

    C. Cattaneo, Atti semin. mat. fis. univ. modena 3, 3 1948, CR Acad. Sci. Paris 247, 431 (1948)

    work page 1948

  43. [43]

    Vernotte, Les paradoxes de la theorie continue de l’equation de la chaleur, Comptes rendus 246, 3154 (1958)

    P. Vernotte, Les paradoxes de la theorie continue de l’equation de la chaleur, Comptes rendus 246, 3154 (1958)

    work page 1958

  44. [44]

    Chester, Second sound in solids, Physical Review 131, 2013 (1963)

    M. Chester, Second sound in solids, Physical Review 131, 2013 (1963)

    work page 2013

  45. [45]

    R. A. Guyer and J. Krumhansl, Solution of the linearized phonon boltzmann equation, Physical Review 148, 766 (1966)

    work page 1966

  46. [46]

    Guyer and J

    R. Guyer and J. Krumhansl, Thermal conductivity, second sound, and phonon hydrodynamic phenomena in nonmetallic crystals, Physical Review 148, 778 (1966)

    work page 1966

  47. [47]

    D. D. Joseph and L. Preziosi, Heat waves, Reviews of modern physics 61, 41 (1989)

    work page 1989

  48. [48]

    Simoncelli, N

    M. Simoncelli, N. Marzari, and A. Cepellotti, Generalization of fourier’s law into viscous heat equations, Physical Review X 10, 011019 (2020)

    work page 2020

  49. [49]

    Hua and L

    C. Hua and L. Lindsay, Space-time dependent thermal conductivity in nonlocal thermal trans- port, Physical Review B 102, 104310 (2020)

    work page 2020

  50. [50]

    Chiloyan, S

    V. Chiloyan, S. Huberman, Z. Ding, J. Mendoza, A. A. Maznev, K. A. Nelson, and G. Chen, Green’s functions of the boltzmann transport equation with the full scattering matrix for phonon nanoscale transport beyond the relaxation-time approximation, Physical Review B 104, 245424 (2021)

    work page 2021

  51. [51]

    Malviya and N

    N. Malviya and N. K. Ravichandran, Efficient calculation of phonon dynamics through a low-rank solution of the boltzmann equation, arXiv preprint arXiv:2502.00337 (2025)

    work page arXiv 2025

  52. [52]

    Zwanzig, Ensemble method in the theory of irreversibility, The Journal of Chemical Physics 40 33, 1338 (1960)

    R. Zwanzig, Ensemble method in the theory of irreversibility, The Journal of Chemical Physics 40 33, 1338 (1960)

    work page 1960

  53. [53]

    Zwanzig, Memory effects in irreversible thermodynamics, Physical Review 124, 983 (1961)

    R. Zwanzig, Memory effects in irreversible thermodynamics, Physical Review 124, 983 (1961)

    work page 1961

  54. [54]

    Zwanzig, Nonequilibrium statistical mechanics (Oxford university press, 2001)

    R. Zwanzig, Nonequilibrium statistical mechanics (Oxford university press, 2001)

    work page 2001

  55. [55]

    D. E. Crawford, Y. Zeng, J. Vidal, and J. Dong, Theory of transient heat conduction, arXiv preprint arXiv:2501.00679 (2024)

    work page arXiv 2024

  56. [56]

    Sussmann and A

    J. Sussmann and A. Thellung, Thermal conductivity of perfect dielectric crystals in the absence of umklapp processes, Proceedings of the Physical Society 81, 1122 (1963)

    work page 1963

  57. [57]

    Mahan and F

    G. Mahan and F. Claro, Nonlocal theory of thermal conductivity, Physical Review B 38, 1963 (1988)

    work page 1963

  58. [58]

    Sobolev, Equations of transfer in non-local media, International Journal of Heat and Mass Transfer 37, 2175 (1994)

    S. Sobolev, Equations of transfer in non-local media, International Journal of Heat and Mass Transfer 37, 2175 (1994)

    work page 1994

  59. [59]

    P. T. Alvarez, Thermal transport in semiconductors: first principles and phonon hydrodynam- ics (Springer, 2018)

    work page 2018

  60. [60]

    A. P. Bart´ ok, J. Kermode, N. Bernstein, and G. Cs´ anyi, Machine learning a general-purpose interatomic potential for silicon, Physical Review X 8, 041048 (2018)

    work page 2018

  61. [61]

    Onsager, Reciprocal relations in irreversible processes

    L. Onsager, Reciprocal relations in irreversible processes. i., Physical review 37, 405 (1931)

    work page 1931

  62. [62]

    Onsager, Reciprocal relations in irreversible processes

    L. Onsager, Reciprocal relations in irreversible processes. ii., Phys. Rev. 38, 2265 (1931). 41

    work page 1931