REVIEW 3 major objections 5 minor 43 references
General relativistic heat flow from first order hydrodynamics
T0 review · 3 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read A non-viscous fluid in normal flow on static or stationary spacetimes obeys a heat equation for the redshifted temperature, with analytic solutions for four black hole backgrounds.
desk verdict A clean, correct extension of heat-flow hydrodynamics to stationary black holes, with a narrow scope that the authors mostly acknowledge. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the first-order heat-flux constitutive relation together with the normal-flow condition $u_a=-N\nabla_a t$. Under static or stationary symmetry and radial heat flow, the $\tau_q$ corrections in that relation drop out, leaving $q^a=-\kappa T\nabla^a\ln(T\sqrt{N^2})$. Substituting this into the energy-momentum conservation law $\nabla_a q^a+q^a\dot{u}_a=0$ turns the conservation statement into $\nabla_a(\sqrt{N^2}\,q^a)=0$ and then into the scalar heat equation for the redshifted temperature, $\nabla_a[\kappa\nabla^a(T\sqrt{N^2})]=0$. The factor $N^2$ is exactly the Tolman--Ehrenfest redshift factor.
What would settle it
Reduce the paper's framework to a concrete test on Schwarzschild: solve the full first-order hydrodynamic equations for a non-viscous fluid whose four-velocity has a nonzero radial component. If $\nabla_a(\sqrt{N^2}\,q^a)=0$ and the profile (31) fail under that tilted flow, the normal-flow assumption is load-bearing; if they survive, the result extends beyond the stated restriction.
Extended reading notes
Core claim
The paper's central claim is that gravity modifies heat conduction through the same redshift factor that governs thermal equilibrium: in the static and stationary cases considered, it is not $q^a$ but $\sqrt{N^2}\,q^a$ that is divergence-free, $\nabla_a(\sqrt{N^2}\,q^a)=0$. Combined with the reduced heat-flux law $q^a=-\kappa\,T\,\nabla^a\ln(T\sqrt{N^2})$, this yields the heat equation $\nabla_a[\kappa\,\nabla^a(T\sqrt{N^2})]=0$ for the redshifted temperature. The paper solves this equation for radial heat flow on Schwarzschild, Schwarzschild--de Sitter, Kerr, and Kerr--de Sitter backgrounds and obtains explicit temperature and chemical-potential profiles. In the limit of vanishing heat flux the profiles reproduce the Tolman--Ehrenfest relation and Klein's law.
Load-bearing premise
The derivation assumes the fluid four-velocity is normal to $t=\text{constant}$ slices and that heat flows only radially; these restrictions are what make the $\tau_q$ terms in the heat-flux law vanish, and the paper explicitly notes that relaxing the normal-flow condition would require a separate analysis.
Editorial extensions
If this is right
- On static and stationary backgrounds the conserved heat current is $\sqrt{N^2}\,q^a$, so the same redshift factor that fixes equilibrium temperature also controls heat transport away from equilibrium.
- Vanishing heat flux reduces the heat equation to $T\sqrt{N^2}=\text{constant}$, recovering the Tolman--Ehrenfest relation, and the chemical-potential equation recovers Klein's $\mu/T=\text{constant}$.
- Radial temperature profiles are obtained analytically for Schwarzschild, Schwarzschild--de Sitter, Kerr, and Kerr--de Sitter spacetimes, with the Kerr and de Sitter results reducing to the Schwarzschild and Schwarzschild--de Sitter expressions in the $a\to0$ and $\Lambda\to0$ limits.
- The entropy production rate between two surfaces at different redshifted temperatures is proportional to $q_0(1/(T_2\sqrt{N^2_2})-1/(T_1\sqrt{N^2_1}))$, so second-law consistency is tied to the redshift factors.
- The analytic profiles can be used as benchmarks for numerical codes in general-relativistic dissipative hydrodynamics.
Reading between the lines
- The same normal-flow reduction should apply to any static or stationary spacetime with a timelike Killing vector, so the profile method likely extends to neutron stars and accretion tori, a direction the paper hints at but does not work out.
- Testing whether $\nabla_a(\sqrt{N^2}\,q^a)=0$ survives angular heat flux on an axisymmetric background would decide whether the conservation law is a general geometric identity or an artifact of radial slicing.
- The paper's first-order-hydrodynamics profiles can be compared directly with the earlier static derivation in [24]; different predictions near the horizon would quantify how much causality and stability conditions alter heat transport.
- Coupling the same heat equation to a charged fluid would add a source term from the electromagnetic field, which could produce temperature profiles relevant to magnetized accretion flows.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper derives a relativistic heat-flow equation for a non-viscous fluid using the BDN first-order hydrodynamics, restricted to static and stationary backgrounds, to a hypersurface-orthogonal 'normal flow' four-velocity, and to purely radial heat flux. The central results are the conservation law \nabla_a(\sqrt{N^2}q^a)=0 (Eq. 14) and the heat equation \nabla_a[\kappa\nabla^a(T\sqrt{N^2})]=0 (Eq. 19), together with analytic temperature profiles for Schwarzschild, Schwarzschild-de Sitter, Kerr, and Kerr-de Sitter spacetimes. The Tolman-Ehrenfest law is recovered as the q_0=0 limit, and chemical-potential profiles are discussed.
Significance. If the central derivation is correct, the paper provides explicit analytic benchmarks for heat transport in BDN first-order hydrodynamics on black-hole spacetimes, including stationary backgrounds, which goes beyond earlier Eckart-based treatments. A strength of the paper is that the main derivation is transparent and does not require an equation of state; the limitations of normal flow and radial heat flux are acknowledged in the conclusions. However, the manuscript contains a number of algebraic and notational errors that currently prevent the results from being taken at face value, especially in the chemical-potential section and in the labeling of inverse metric components for the stationary applications.
major comments (3)
- [Sec. III.A, Eq. (29)] The second equality in Eq. (29) does not follow from Eqs. (26)-(28). For Schwarzschild, expanding the first equality of Eq. (29) to first order in q0 gives terms proportional to A/T0 and to T0, namely d(mu sqrt(-g00))/dr = -q0 A/(kappa T0 r^2 sqrt(f)) + q0 T0/(alpha r^2 f^{3/2}) + O(q0^2), whereas the displayed second equality is O(q0^2). Thus the simplified form is false unless the integration constants A and T0 are set to zero, which they are not. The chemical-potential results in Appendix C are based on Eq. (30), so they may survive, but the claimed simplification in Eq. (29) must be corrected or removed.
- [Sec. III.B, Eqs. (41) and (49)] The quantity labeled g00 in Eq. (41) is actually the inverse time-time component g^00. Substituting the displayed expression into N^2 = -1/g^00 gives the correct Kerr lapse and reproduces the Schwarzschild limit N^2 = 1 - 2m/r only if the expression is read as g^00. As printed with g00, the definition of N^2 and all subsequent Kerr temperature formulas are inconsistent. The same care is needed for Eq. (49), whose a=0 limit appears to give the metric component g00 rather than g^00; the KdS formulas should be re-verified with the inverse metric component.
- [Appendix A.2, Eq. (A11)] The argument that the last two terms of Eq. (A11) vanish is written with incorrect index placement: the text states qb xi^b_(t) = qr xi^r_(t) = 0 because xi^r_(t)=0, but q_b xi^b_(t) = q_t, not q_r xi^r_(t). The intended vanishing terms are q^b xi_(t)b = g_{0b} q^b = g_{0r} q^r and q^b xi_(phi)b = g_{phi r} q^r, which are zero only because the stationary metric (37) has no g_{0r} or g_{phi r} cross terms and q^a is purely radial. The conclusion \pounds_u q_c=0 is correct under these stated restrictions, but the proof should be rewritten and the no-r-cross-term condition made explicit.
minor comments (5)
- [Sec. IV, Eqs. (62) and (63)] The proper-time factors mix g00 and g^00 notation: for a static metric d\tau = dt/u^0 = dt/\sqrt{-g^{00}} = dt\sqrt{-g_{00}}, but the text writes this in a way that appears to equate dt/\sqrt{-g_{00}} with dt\sqrt{-g_{00}}. Please use explicit superscripts to distinguish g^{00} from g_{00} throughout this section.
- [Sec. III.A, Eq. (27)] The displayed formula is ambiguous in the typeset text: the correct factor is \kappa r^2\sqrt{-g_{00}}/\sqrt{g_{rr}} times d(T\sqrt{N^2})/dr, not \kappa r^2\sqrt{-g_{00}}\sqrt{g_{rr}}. Please check the typesetting so the denominator is unambiguous.
- [Sec. III.B, text before Eq. (41)] After Eq. (39), the statement that N^2 = -1/g00 should read N^2 = -1/g^{00}; in static spacetimes the two coincide because g^{00}=1/g_{00}, but in stationary spacetimes they do not. This is directly related to the major comment on Eq. (41).
- [Sec. I and Sec. V] The comparison with Ref. [24] is fair, but the abstract and title should state more explicitly that the stationary results are limited to normal flow and radial heat flux; the conclusions already acknowledge this, and the abstract should not imply general stationary heat flow.
- [Appendix C, Eq. (C2)] The displayed expression for (mu sqrt(1-2m/r))_{sch} has a parenthesis structure that is hard to parse; please check the bracket matching and the sign of the q0/(m alpha) term against Eq. (C1).
Circularity Check
No significant circularity: the heat flow equation is derived from the BDN conservation equations, and the Tolman-Ehrenfest relation emerges as the q0=0 limit of the derived solutions rather than being assumed.
full rationale
The paper's derivation chain is self-contained relative to its stated inputs. It begins with the BDN constitutive relations (Eqs. 1-6), imposes a non-viscous fluid and normal flow u_a = -N grad_a t, and then uses Killing symmetries and the Frobenius hypersurface-orthogonality condition to show that L_u q^c = 0, giving Eq. (A10) and hence the simplified heat flux q^a = alpha grad^a(mu/T) in Eq. (16). Appendix B converts this to q^a = -kappa T grad^a ln(T sqrt(N^2)) in Eq. (18) using only thermodynamic identities and the already-derived Eq. (A10). Substitution into Eq. (12) yields the conservation law grad_a(sqrt(N^2) q^a) = 0 (Eq. 14) and the heat equation (Eq. 19). None of these steps assumes the target result: the Tolman-Ehrenfest relation is not an input but is recovered only in the q0 = 0 limit of the solved temperature profiles, e.g. Eq. (31) and Eq. (45), and Klein's law is similarly a q0 -> 0 consequence of Eq. (30). No parameter is fitted to a subset of data and then renamed a prediction; the integration constants T0 and q0 are genuine constants of integration. The only self-citation (Ref. 22) appears in the introduction as background on earlier work obtaining the TE and Klein relations from first-order hydrodynamics, and it is not used as a load-bearing premise in the derivation of Eqs. (14), (18), or (19). The paper explicitly acknowledges in the Conclusions that the results are restricted to normal flow and, in the stationary case, to radial heat flux; Appendix A.2 shows that the latter assumption is needed to make the last two terms of Eq. (A11) vanish. This is an honest scope restriction, not a circular step. The central derivation is therefore independent of any fitted output or self-citation chain, and no circularity is present.
Assumptions & free parameters
free parameters (4)
- q0 (heat flux integration constant) =
arbitrary real
- T0 (temperature integration constant) =
positive constant
- kappa (heat conductivity) assumed constant =
constant
- alpha = sigma T (epsilon+p)/n assumed constant =
constant
assumptions (7)
- domain assumption BDN first-order constitutive relations (Eqs. 1-6) describe causal and stable relativistic hydrodynamics.
- domain assumption The fluid is non-viscous, with dissipation only through heat flux.
- ad hoc to paper The fluid four-velocity is hypersurface-orthogonal normal to t=constant slices.
- domain assumption Fluid scalars respect the background spacetime symmetries, giving u^a grad_a X=0 and Lie derivatives of q along Killing vectors equal to zero.
- ad hoc to paper Heat flow is purely radial.
- domain assumption Equilibrium thermodynamic relations (Euler relation and first law) apply to the out-of-equilibrium fluid.
- domain assumption The background metric is fixed and the fluid's backreaction is neglected.
Cite this review
Pith. "Pith review of General relativistic heat flow from first order hydrodynamics." pith.science (2026). https://pith.science/paper/6LT62OYK
@misc{pith2026241202364,
author = {Pith},
title = {Pith review of: General relativistic heat flow from first order hydrodynamics},
year = {2026},
howpublished = {\url{https://pith.science/paper/6LT62OYK}},
note = {Machine review of arXiv:2412.02364}
}
read the original abstract
Following the recently proposed stable and causal first-order relativistic hydrodynamics by Bemfica, Disconzi, and Noronha, we find the heat flow equation in the presence of gravity for a non-viscous fluid, which suffers heat dissipation. The derivation is confined to static and stationary backgrounds. We find that in the presence of gravity, the heat flux times a redshift factor is conserved. Then for radial heat flow, the temperature profiles are obtained from the heat equation when the gravity is sourced by Schwarzschild, Schwarzschild-dS, Kerr and Kerr-dS black holes, respectively. Consequently, the chemical potential profile is also discussed.
Reference graph
Works this paper leans on
-
[24]
C. Rovelli and M. Smerlak, Class. Quant. Grav. 28, 075007 (2011), 1005.2985
arXiv 2011
-
[22]
Gravitational instability caused by the weight of heat
Z. Roupas, Symmetry 11, 1435 (2019), 1809.04408
work page Pith review arXiv 2019
-
[1]
Application to static backgrounds To get a feel of the above results we apply them on a few specific static metrics. Schwarzschild Geometry: – In this section, we take the Schwarzschild geometry as an example of the static spherically symmetric case. On solving Eq. (26) for the Schwarzschild case with g00 = −1 + 2m r and grr = −1 g00 , we get T (r) = 1√ 1 ...
-
[2]
Application to stationary metrics Kerr Geometry: – Let us now turn to an example of a stationary Kerr metric. The metric elements in this case in Boyer-Lindquist coordinates are given by g00 = − ( 1 − 2mr ρ2 ) ; g0φ = − 2mar sin2 θ ρ2 ; grr = ρ2 ∆ ; gθθ = ρ2; gφφ = ( r2 + a2 + 2ma2r sin2 θ ρ2 ) sin2 θ,(39) with ρ2 ≡ r2 + a2 cos2 θ, ∆ ≡ r2 − 2mr + a2, and ...
work page 2020
-
[3]
Static background As discussed in Section III, the equations governing the conservation of the energy-momentum tensor, along with the conditions for a non-viscous fluid, are given by Eq. (9) and Eq. (10). For a static background given by Eq. (11), we have ua∇aX = 0 and ∇aua = 0, where X represents any scalar fluid variable. Consequently, under these conditi...
-
[4]
(10) for the stationary background (37), simplifies to Eq
Stationary background Again the non-viscous fluid equation, given by Eq. (10) for the stationary background (37), simplifies to Eq. (A1) due to ua∇aε = 0 and ∇aua = 0. Since in this case our choice of ua also satisfies the Frobenius theorem, we will have Eq. (A7) for the stationary background as well. The background has two Killing vectors, ξb (t) = (1, 0, 0...
-
[5]
Schwarzschild black hole Solving Eq. (30) for the Schwarzschild case with g00 = −1/g rr = − ( 1 − 2m r ) and assuming α to be constant, we have µ T = A + q0 α ∫ 1 r2( 1 − 2m r ) 3/ 2 dr . (C1) Further, on using Eq. (31), we have ( µ (r) √ 1 − 2m r ) (sch) = ( T0 − q0 2mκ ln ( 1 − 2m r ) ) × ( A − q0 m α 1√ 1 − 2m r ) , (C2) where A is the integration constant
-
[6]
Schwarzschild-de Sitter black hole Again for the Schwarzschild-de Sitter case, Eq. (30) assuming α to be constant reduces to, µ T = A + q0 α ∫ 1 r2( 1 − 2m r − Λ r2 3 ) 3/ 2 dr , (C3) which doesn’t have a simplified or compact form. How- ever, it can be solved near the horizon, using r → r′ + rb, where rb is the black hole horizon. The above equation near ...
Show all 43 references
-
[7]
Tolman and P
R. Tolman and P. Ehrenfest, Phys. Rev. 36, 1791 (1930)
1930
-
[8]
R. C. Tolman, Phys. Rev. 35, 904 (1930)
1930
-
[9]
H. A. Buchdahl, Phys. Rev. 76, 427 (1949), URL https://link.aps.org/doi/10.1103/PhysRev.76.427.2
1949 doi
- [10]
-
[11]
Klein, Rev
O. Klein, Rev. Mod. Phys. 21, 531 (1949), URL https://link.aps.org/doi/10.1103/RevModPhys.21.531
1949 doi
-
[12]
W. J. Cocke, Annales de l’I.H.P. Physique th´ eorique 2, 283 (1965), URL http://eudml.org/doc/75503
1965
-
[13]
Katz and Y
J. Katz and Y. Manor, Phys. Rev. D 12, 956 (1975), URL https://link.aps.org/doi/10.1103/PhysRevD.12.956
1975 doi
-
[14]
R. D. Sorkin, R. M. Wald, and Z. J. Zhang, Gen. Rel. Grav. 13, 1127 (1981)
1981
-
[15]
Gao, Phys
S. Gao, Phys. Rev. D 84, 104023 (2011), [Addendum: Phys.Rev.D 85, 027503 (2012)], 1109.2804
2011 arXiv
-
[16]
Gao, Phys
S. Gao, Phys. Rev. D 85, 027503 (2012), URL https://link.aps.org/doi/10.1103/PhysRevD.85.027503
2012 doi
-
[17]
S. R. Green, J. S. Schiffrin, and R. M. Wald, Class. Quant. Grav. 31, 035023 (2014), 1309.0177
2014 arXiv
-
[18]
K. Shi, Y. Tian, X. Wu, H. Zhang, and C. Zhu, Class. Quant. Grav. 39, 085004 (2022), 2108.08729
2022 arXiv
-
[19]
K. Shi, Y. Tian, X. Wu, H. Zhang, and J. Zhang (2022), 2211.11574
2022 arXiv
-
[20]
Roupas, Class
Z. Roupas, Class. Quant. Grav. 30, 115018 (2013), [Erra- tum: Class.Quant.Grav. 32, 119501 (2015)], 1411.0325
2013 arXiv
- [21]
- [23]
- [25]
-
[26]
J. A. S. Lima, A. Del Popolo, and A. R. Plastino, Phys. Rev. D 100, 104042 (2019), 1911.09060
2019 arXiv
-
[27]
J. A. S. Lima and J. Santos, Phys. Rev. D 104, 124089 (2021), 2112.12282
2021 arXiv
-
[28]
Ram and B
B. Ram and B. R. Majhi, Laws of thermodynamic equilib- rium within first order relativistic hydrodynamics (2024), 2304.11843, URL https://arxiv.org/abs/2304.11843
2024 arXiv
- [29]
-
[30]
Miranda, M
M. Miranda, M. Rinaldi, and V. Faraoni, Phys. Rev. D 110, 104065 (2024), URL https://link.aps.org/doi/10.1103/PhysRevD.110.104065
2024 doi
-
[31]
L. Cui, X. Hao, and L. Zhao (2024), 2411.03094
2024 arXiv
-
[32]
Landau and E
L. Landau and E. Lifshitz, 2nd ed. (Butterworth- Heinemann, London) 6, 552 (1987)
1987
-
[33]
Eckart, Phys
C. Eckart, Phys. Rev. 58, 919 (1940)
1940
-
[34]
W. A. Hiscock and L. Lindblom, Annals Phys. 151, 466 (1983)
1983
-
[35]
Israel, Annals Phys
W. Israel, Annals Phys. 100, 310 (1976)
1976
-
[36]
Israel and J
W. Israel and J. M. Stewart, Annals Phys. 118, 341 (1979)
1979
-
[37]
F. S. Bemfica, M. M. Disconzi, and J. Noronha, Phys. Rev. X 12, 021044 (2022), URL https://link.aps.org/doi/10.1103/PhysRevX.12.021044
2022 doi
-
[38]
F. S. Bemfica, M. M. Disconzi, and J. Noronha, Phys. Rev. D 98, 104064 (2018), URL https://link.aps.org/doi/10.1103/PhysRevD.98.104064
2018 doi
-
[39]
F. S. Bemfica, M. M. Disconzi, and J. Noronha, Phys. Rev. D 100, 104020 (2019), URL https://link.aps.org/doi/10.1103/PhysRevD.100.104020
2019 doi
- [40]
-
[41]
R. E. Hoult and P. Kovtun, JHEP 06, 067 (2020), 2004.04102
2020 arXiv
-
[42]
S. W. Hawking, Commun. Math. Phys. 43, 199 (1975), [Erratum: Commun.Math.Phys. 46, 206 (1976)]
1975
- [43]
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.