{"id":"318ee127-a2d9-4af8-8bf7-95e469e83be7","arxiv_id":"2608.13339","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A unified theoretical review shows that viscous heat equations, derived from the phonon Boltzmann equation, reduce to standard models and predict exotic hydrodynamic thermal phenomena such as backflow and vortices.","lead":"This paper reviews and connects the quantum, semiclassical, and mesoscopic theories of heat flow in crystals, from the Dyson equation to the viscous heat equations used to simulate phonon hydrodynamics in devices. It explains how a unified framework predicts unusual thermal phenomena such as heat flowing backward and forming vortices, which could help engineers design better thermal management in nanoscale electronics.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The VHE derivation rests on a finite spectral gap for the normal-scattering operator and on lambda/g << 1; this is asserted rather than established, and the gap can fail in the thermodynamic limit for acoustic phonons.","rationale":"The reader's weakest-assumption analysis identifies the same broad step: the spectral-gap separation and the smallness of Umklapp relative to that gap. I agree that this is the central condition for the VHE derivation. My stress test sharpens the concern in two ways not explicit in the reader's verdict: (i) the existence of a finite gap is not guaranteed by the absence of additional conservation laws, because gapless acoustic phonons can produce eigenvalue accumulation at zero; and (ii) the high-temperature diffusive limit, where the VHE should reduce to Fourier's law, is precisely a regime where lambda/g << 1 is unlikely to hold, so the decomposition of the conductivity into momentum and diffusion-damped parts may be uncontrolled even though the final steady-state equation is formally correct. These points do not overturn the review's positive evidence: the derivation from the Kadanoff-Baym equation is internally consistent, the VHE reduce to the DPLE and GKE in transparent limits, and the graphite results are benchmarked against both first-principles LBTE solutions and experimental second-sound and lattice-cooling measurements. The concern is therefore a genuine caveat about scope and convergence, not a demonstrated failure of the central argument. Since the reader already assigned CONDITIONAL with medium correctness risk, my analysis does not move the verdict; it reinforces the conditionality. The proposed numerical test would settle whether the spectral-gap assumption holds in practice beyond the cited benchmarks.","tokens_in":53541,"tokens_out":6833,"duration_ms":84391,"concrete_test":"Compute the spectrum of the symmetrized normal-scattering operator from first-principles scattering matrices for graphite and one additional high-conductivity material (e.g., diamond) at representative temperatures (e.g., 70 K and 300 K), using increasing q-grid densities (24^3 to 64^3). Evaluate g_min as the smallest nonzero eigenvalue, lambda = ||Omega~_U||, and the off-diagonal M-D contribution to kappa that is omitted in Eq. (137). If lambda/g_min is not small for a converged grid, or if the mixed-block contribution shifts kappa by more than about 10%, the Schrieffer-Wolff decomposition underlying the VHE coefficients is not justified in that regime. Additionally, check whether g_min decreases systematically with increasing grid density; if it tends to zero in the thermodynamic limit, the finite-gap assumption itself fails, independent of the Umklapp strength.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is in Section 4.3.1: after Eq. (136), the review assumes that the projected normal-scattering operator has a finite spectral gap g separating the three momentum relaxons from all remaining modes, and that Umklapp scattering is a small perturbation, lambda/g << 1, with lambda = ||Omega~_U||. This justifies the Schrieffer-Wolff block diagonalization and the conductivity decomposition kappa = kappa_MM + kappa_DD in Eq. (137), with corrections O(lambda/g) in Eq. (138). The claim that g must be finite, 'otherwise additional conserved quantities would exist,' is not a proof: for acoustic phonons, normal-scattering eigenvalues can accumulate at zero as q approaches 0 because long-wavelength lifetimes diverge, without any extra exactly conserved mode. Thus g can vanish in the thermodynamic limit even though the only conserved quantities are energy and momentum. In the high-temperature diffusive regime, where the VHE are supposed to reduce to Fourier's law, Umklapp need not be a small perturbation, so the mixed M-D blocks neglected in Eq. (137) can contribute to the conductivity; defining kappa_D as the exact remainder in Eq. (142) may restore the steady-state conductivity while still altering the transient and backflow predictions that are central to the review's claims. The graphite benchmarks support the assumption in one material and temperature window, but they do not establish the general validity needed for the unified mesoscopic framework.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This review article aims to present a unified theoretical description of thermal transport in insulating crystals from the quantum regime to mesoscopic device modeling. It derives the phonon Boltzmann transport equation from the Kadanoff-Baym equation using the Wigner representation and the generalized Kadanoff-Baym ansatz, introduces the relaxon formalism for the linearized Boltzmann equation, and then constructs the viscous heat equations (VHE) by projecting the LBTE onto the subspace of energy and momentum relaxons. The central claim is that the VHE constitute a general mesoscopic framework that reduces to Fourier's law, the dual-phase-lag equation, and the Guyer-Krumhansl equations in appropriate limits, and that they predict Poiseuille heat flow, second sound, steady-state heat backflow, vortices, and transient lattice cooling. The parameters entering the VHE are computed from first principles, and predictions are benchmarked against full LBTE solutions and against transient-thermal-grating and pump-probe experiments in graphite.","tokens_in":53772,"tokens_out":5943,"duration_ms":61742,"significance":"If the central claims hold, this review would provide a valuable bridge between nonequilibrium Green's function theory and practical mesoscopic modeling of phonon hydrodynamics. The manuscript is strongest where it makes the approximation chain explicit: from the Kadanoff-Baym equation to the BTE, from the BTE to the relaxon picture, and from the relaxon picture to the VHE. The introduction of thermal viscosity from even relaxons and the analytical reduction of steady-state VHE to modified Helmholtz and biharmonic equations are genuinely useful contributions. The benchmarks against spatially and temporally resolved LBTE solutions and against experimental data in graphite are also a clear strength, as is the explicit discussion of boundary conditions and finite-size corrections. The main caveat is that the 'unified framework' claim rests on a spectral-gap separation that is asserted rather than established, which limits the demonstrated generality of the VHE beyond the specific benchmarks presented.","major_comments":[{"comment":"The assertion that the projected normal-scattering operator has a finite spectral gap g because 'otherwise additional conserved quantities would exist' is not a valid proof in the infinite-dimensional phonon space. Eigenvalues of the normal-scattering operator can accumulate at zero as q approaches zero for acoustic phonons with divergent long-wavelength lifetimes, without any additional exactly conserved mode. Since the Schrieffer-Wolff block diagonalization and the decomposition kappa = kappa_MM + kappa_DD in Eq. (137) rely on a positive gap, the authors should either prove the existence of the gap from the structure of the normal-scattering operator, provide numerical spectra for representative materials, or explicitly restrict the VHE derivation to regimes or geometries where a minigap is known to exist.","section":"Section 4.3.1, Eqs. (136)-(138)"},{"comment":"The smallness condition lambda/g << 1, with lambda = ||Omega~_U||, is stated but not quantified. In the high-temperature diffusive regime, where the VHE are claimed to reduce exactly to Fourier's law following Eq. (143), Umklapp scattering is typically not a small perturbation to the normal-scattering spectrum. In that regime the mixed M-D blocks omitted in Eq. (137) can contribute at order one to the conductivity and, more importantly, can modify the transient and backflow predictions that are central to the review's claims. The authors should report numerical estimates of lambda/g across the temperature range for the materials considered, or qualify the generality of the VHE framework accordingly.","section":"Section 4.3.1, after Eq. (138)"},{"comment":"The supporting benchmarks are limited to graphite in a relatively narrow temperature window around 70-120 K. The claim that the VHE constitute a general mesoscopic framework for phonon hydrodynamics should be supported by evidence that the spectral-gap condition and the smallness of lambda/g also hold for other materials discussed in the text, such as diamond, silicon, and hBN, or by providing explicit criteria for when the projection fails. Without such evidence, the 'unified' status of the VHE remains a conjecture whose demonstrated range is narrower than the review's title and abstract suggest.","section":"Section 4, Figs. 13, 18, 19"}],"minor_comments":[{"comment":"The text states 'sigma_x = 2 micron and sigma_x = 2.8 micron'; the second occurrence should presumably be sigma_y, since the Gaussian heating profile is defined with widths in both the x and y directions.","section":"Section 4.10, text after Eq. (185)"},{"comment":"The claim that 'all transport coefficients ... without relying on experimental inputs or fitting procedures' is slightly overstated, because the finite slip length b = 0.4 micron used in Section 4.10.1 is obtained by matching VHE results to LBTE solutions with diffusive boundaries. The text should distinguish between first-principles transport coefficients and fitted boundary parameters.","section":"Section 4.3 and Section 4.10.1"},{"comment":"Equation (138) reports off-diagonal corrections of order O(lambda/g^2) and diagonal corrections O(lambda^2/g^4) in the inverse of the transformed operator, while Eq. (140) states kappa_MM = kappa_M + O(lambda/g). These order-of-magnitude notations should be harmonized so that the size of the neglected terms is unambiguous.","section":"Eqs. (138) and (140)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a carefully written review, but it relies very heavily on the authors' own previous work, Refs. [18] and [89], for the central VHE construction and its benchmarks. This is natural in a review, but it means that the load-bearing spectral-gap assumption currently lacks independent support. I would encourage the editor to require an explicit treatment of the spectral-gap condition, including numerical evidence for its range of validity, before publication. The paper otherwise fits the scope of a review journal and contains several valuable derivations and benchmarks."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe short version: this is a genuine review, not a disguised research paper, and it does its job well. The authors trace a coherent line from phonon Dyson and Kadanoff-Baym equations through the linearized Boltzmann equation to the viscous heat equations (VHE) and their limiting cases (DPLE, Guyer-Krumhansl, Fourier). The derivations are internally consistent, and the paper includes a new analytical recasting of the steady-state VHE as decoupled Helmholtz and biharmonic equations, which is a useful addition. The benchmarks against first-principles LBTE solutions and experimental second-sound data in graphite give the framework credible empirical grounding.\n\nCredit where due: the review is clearly written, the hierarchy of models is illuminating, and the discussion of when Fourier fails is concrete. The explicit limitations section is honest about linear-response and gradient assumptions, though it does not address the deeper spectral-gap question.\n\nThe soft spot, as you would expect, is the projection that produces the VHE. The review assumes a finite spectral gap g separating the three momentum relaxons from all other modes of the normal-scattering operator, and that Umklapp strength lambda satisfies lambda/g << 1. This is asserted rather than established. The stress-test raises a fair mathematical point: for acoustic phonons, normal-scattering eigenvalues can accumulate at zero as q approaches zero because long-wavelength lifetimes diverge, so the gap can vanish in the thermodynamic limit without any extra conserved quantity. That means the Schrieffer-Wolff block diagonalization and the clean split into momentum and diffusion-damped conductivities are not guaranteed beyond the specific benchmarked windows. The paper's material-specific benchmarks (graphite, around 70 K) are suggestive but not a general proof.\n\nAnother caveat: the abstract says 'we demonstrate' for effects that mostly come from earlier papers by the same group. The review is not a new result; it is a synthesis. That is fine for a review, but readers should not mistake the exotic predictions (backflow, vortices, negative thermal resistance) for experimentally confirmed phenomena. The steady-state backflow is a VHE prediction, backed by LBTE comparison in a specific device geometry, but no experiment has seen it yet.\n\nWho benefits: newcomers and practitioners who want a map from the many-body formalism to continuum hydrodynamics, and a clear statement of where each model lives. The analytical section is handy for device modeling. I would assign it to a referee with a background in kinetic theory; that referee should push on the gap condition or at least ask for a supplementary analysis of its range of validity.\n\nVerdict: accept for peer review, with revisions expected. The paper is a useful reference and deserves publication after the gap issue is addressed, even if only by sharpening the caveats and supplying more evidence of where the condition holds.","headline":"A competent and useful review of phonon hydrodynamics, but the load-bearing spectral-gap assumption in the VHE derivation is asserted rather than proven, and the reviewer should push on it.","tokens_in":54356,"tokens_out":2569,"would_cite":true,"duration_ms":27403,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["66.70.-f","63.20.kg"],"model":"deepseek-v4-flash","headline":"Projecting the Boltzmann equation onto conserved phonon modes yields viscous heat equations that unify Fourier diffusion with second sound, backflow, and vortices.","keywords":["phonon hydrodynamics","viscous heat equations","thermal transport","relaxons","Boltzmann transport equation","second sound","thermal backflow","Kadanoff-Baym equations"],"falsifier":"Take a material where hydrodynamic transport is claimed (for example graphite between 70 K and 225 K), compute the full normal-scattering and Umklapp scattering matrices from first principles, and evaluate the spectral gap $g$ of the normal operator and the operator norm $\\lambda$ of the Umklapp part; if $\\lambda/g$ is not much smaller than one in exactly the temperature and length-scale window where the review predicts backflow, second sound, or lattice cooling, the Schrieffer-Wolff step is invalid and the VHE predictions in that window lack their stated foundation. A complementary experimental check: in the tunnel-chamber graphite device, time-resolved thermal imaging should reveal the predicted steady-state temperature inversion inside the chamber, and its absence at the predicted temperature would directly contradict the framework.","tokens_in":53250,"feed_emoji":"🔥","tokens_out":12173,"duration_ms":107698,"temperature":0.7,"pith_summary":"This review tries to establish that the practical levels of thermal transport theory in crystals — from quantum nonequilibrium Green's functions down to Fourier's law — are connected by a chain of controlled approximations, whose mesoscopic endpoint is a single pair of partial differential equations called the viscous heat equations. The authors derive these equations by projecting the linearized phonon Boltzmann transport equation onto the subspace spanned by the collective modes that conserve energy and crystal momentum, which introduces a phonon drift velocity alongside the temperature field. If the derivation holds, one framework reproduces all the familiar limits (Fourier, Cattaneo, dual-phase-lag, Guyer-Krumhansl) and additionally predicts hydrodynamic phenomena Fourier's law cannot see, including steady-state thermal backflow and heat vortices in suitably shaped devices, with those predictions benchmarked against full Boltzmann solutions and recent transient-grating experiments in graphite. The payoff is a continuum model whose coefficients are computable from first principles, so device engineers could predict where heat flows backward or forms vortices without solving the full Boltzmann equation.","feed_headline":"A single theory predicts heat flowing backward in crystals","feed_subtitle":"Derived from conserved phonon modes, viscous heat equations unify diffusion with backflow and vortices.","key_machinery":"The load-bearing object is the relaxon formalism: the symmetrized scattering matrix of the linearized Boltzmann equation is diagonalized exactly, and its eigenvectors, the relaxons, are collective phonon excitations with definite parity and well-defined relaxation times. Four of them sit at zero eigenvalue: the energy relaxon, an eigenvector of the full scattering operator, and three momentum relaxons, zero-eigenvalue eigenvectors of the normal-scattering operator alone. Projecting the LBTE onto this four-dimensional subspace, then applying a Schrieffer-Wolff transformation to block-diagonalize the collision operator — valid when the strength $\\lambda$ of Umklapp scattering is small relative to the spectral gap $g$ separating the momentum modes from all other normal-scattering modes — decomposes the thermal conductivity into a momentum contribution $\\kappa^M$ and a diffusion-damped contribution $\\kappa^D$, and yields the closed-form coefficients of the viscous heat equations. The new coefficient this machinery produces is the thermal viscosity tensor $\\eta$, built from even-parity relaxons, which is what allows the heat flux to carry nonzero curl and thereby exhibit vortices and backflow.","core_discovery":"The paper argues that projecting the linearized phonon Boltzmann transport equation onto the subspace spanned by the four zero-eigenvalue collective modes of the scattering operator — the energy mode and the three crystal-momentum modes of normal scattering, which in the relaxon formalism are the conserved special relaxons — yields a closed pair of partial differential equations, the viscous heat equations, coupling the temperature field $T(\\mathbf{r},t)$ to a phonon drift velocity $\\mathbf{u}(\\mathbf{r},t)$. The resulting heat flux splits into a drifting component carried by momentum-conserving collisions and a diffusion-damped component proportional to the temperature gradient, and thermal viscosity, defined microscopically from the even-parity relaxons, governs momentum diffusion. On the authors' own terms, this single framework contains Fourier's law, Cattaneo's equation, the dual-phase-lag equation, and the Guyer-Krumhansl equations as special limits, and it predicts Poiseuille heat flow, second sound, negative thermal resistance, steady-state thermal backflow, and heat vortices; the review backs these predictions with solutions of the full LBTE in device geometries, analytical Helmholtz and biharmonic solutions, and agreement with transient-grating experiments in graphite.","pith_inferences":["The review leaves implicit that the spectral-gap condition ($\\lambda/g \\ll 1$) is itself a screening descriptor: materials with a wide separation between slow momentum-conserving and fast diffusive collision modes should display hydrodynamic heat transport, so the ratio could be computed and tabulated for candidate crystals before any device simulation.","Because the analytical solutions decompose the temperature field into compressibility and vorticity contributions, phonon devices could plausibly be designed as thermal analogs of fluidic circuits, with chambers and openings acting as vortex generators; testing this would require extending the analytical biharmonic solutions beyond the graphite strip studied here, which the paper does not do.","The review states that current phonon hydrodynamics is confined to the laminar, zero-Reynolds-number regime and that phonon turbulence has not been modeled; a natural open problem this implies is a nonlinear extension of the VHE that restores the advective term, which would predict whether and when collective heat transport becomes turbulent."],"forward_implications":["For device design, the Fourier deviation number, a dimensionless combination of the VHE coefficients, device size, and applied gradients, tells engineers when Fourier's law fails, so that hot spots and temperature inversions can be predicted from the VHE alone.","For theory, the dual-phase-lag and Guyer-Krumhansl equations are demoted from independent models to limiting cases of the VHE (inviscid, and linear-isotropic drift-dominated steady state, respectively), so disagreements between those models are resolved inside one framework.","For experiments, the VHE predict resonant temperature-wave amplification, steady-state temperature inversion in tunnel-chamber geometries, and lattice cooling after pulsed heating, at temperatures and length scales where the full Boltzmann equation agrees, giving observable signatures for thermal-imaging and pump-probe setups.","For computation, all VHE coefficients come from ab initio solutions of the LBTE with no fitting, so the same framework scales to screening materials and device geometries for hydrodynamic heat transport."],"supporting_citations":[{"why":"Original derivation of the viscous heat equations and of the relaxon-based first-principles formulas for thermal conductivity and thermal viscosity that parameterize them.","marker":"[18]"},{"why":"Benchmarks the VHE against full space-dependent LBTE solutions in device geometries and supplies the tunnel-chamber heat-backflow and lattice-cooling predictions.","marker":"[89]"},{"why":"Introduces the relaxon formalism, the exact diagonalization of the scattering operator whose zero-eigenvalue modes define the projection subspace.","marker":"[117]"},{"why":"Provide the Schrieffer-Wolff transformation used to block-diagonalize the collision operator and split the conductivity into momentum and diffusion-damped contributions.","marker":"[211, 212]"},{"why":"Experimental transient-thermal-grating observation of second sound in graphite above 100 K; the quantitative target the VHE must reproduce.","marker":"[24]"},{"why":"Time-resolved measurements extending temperature-wave signatures in graphite to about 200-225 K, benchmarks for the VHE hydrodynamic window.","marker":"[23]"},{"why":"Source of the dual-phase-lag equation, which the review derives as the inviscid limit of the VHE.","marker":"[220]"},{"why":"Reference for the Guyer-Krumhansl equations, shown to follow from the VHE in the linear-isotropic, drift-dominated steady-state limit.","marker":"[228]"},{"why":"Variational exact solution of the LBTE used to validate the relaxon thermal-conductivity expression underlying the VHE coefficients.","marker":"[109]"}],"fun_headline_variants":["Heat backflow and vortices from one phonon theory","Unified phonon theory predicts backward heat flow","Viscous heat equations: from Fourier to backflow","Phonon hydrodynamics: heat flows backward and swirls"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on the premise that momentum-destroying Umklapp scattering is weak compared with the gap separating the three momentum-conserving modes from all other collision modes ($\\lambda/g \\ll 1$); if that separation is not large, the block diagonalization, and with it the closed-form coefficients of the viscous heat equations, is not mathematically guaranteed.","fun_headline_variants_meta":{"raw":{"variants":["Heat backflow and vortices from one phonon theory","Unified phonon theory predicts backward heat flow","Viscous heat equations: from Fourier to backflow","Phonon hydrodynamics: heat flows backward and swirls"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00058,"raw_usage":{"total_tokens":2807,"prompt_tokens":1094,"completion_tokens":1713,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":710,"completion_tokens_details":{"reasoning_tokens":1647}},"tokens_in":710,"tokens_out":1713,"duration_ms":12608,"temperature":1.0,"reasoning_tokens":1647,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:21:56.276933+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a material where hydrodynamic transport is claimed (for example graphite between 70 K and 225 K), compute the full normal-scattering and Umklapp scattering matrices from first principles, and evaluate the spectral gap $g$ of the normal operator and the operator norm $\\lambda$ of the Umklapp part; if $\\lambda/g$ is not much smaller than one in exactly the temperature and length-scale window where the review predicts backflow, second sound, or lattice cooling, the Schrieffer-Wolff step is invalid and the VHE predictions in that window lack their stated foundation. A complementary experimental check: in the tunnel-chamber graphite device, time-resolved thermal imaging should reveal the predicted steady-state temperature inversion inside the chamber, and its absence at the predicted temperature would directly contradict the framework.","supporting_citations":[{"cited_title":"Derivation of a hydrodynamic heat equation from the phonon boltzmann equation for general semiconductors.Physical Review B, 103(14):L140301, 2021","cited_arxiv_id":null,"evidence_quote":"Reference for the Guyer-Krumhansl equations, shown to follow from the VHE in the linear-isotropic, drift-dominated steady-state limit."}],"review_version":1}