REVIEW 4 major objections 4 minor 41 references
Second order causal hydrodynamics in Eckart frame: using gradient expansion scheme
T0 review · 4 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read The paper constructs the general causal second-order constitutive relations for a relativistic non-ideal fluid in the Eckart frame, including spacetime-curvature terms in shear viscosity, bulk viscosity, and heat flow, and shows that…
desk verdict A plausible Eckart-frame extension of the standard gradient expansion, but the 'general' constitutive relations are undercut by omitted terms and the causality section has a dimensional slip. 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 gradient expansion scheme: enumerate all independent second-order scalars, vectors, and symmetric traceless tensors built from transverse gradients of $u^\mu$, $T$, and $\mu$, plus curvature terms arising from the non-commutativity of covariant derivatives (Tables I-III). First-order ideal-fluid equations are used to replace time derivatives by spatial gradients, and first-order flux relations are substituted into second-order terms, converting the algebraic constitutive relations into relaxation-type dynamical equations that make the theory causal.
What would settle it
Check directly whether every structure in Tables I-III can be expressed as a combination of the terms retained in eqs. (17), (22), and (23) using only the ideal-fluid identities; if, for example, $M_1 = D^\alpha_\perp D_{\perp\alpha}\ln T$, $O_1=(D_\perp^{<\mu}\ln T)(D_\perp^{\nu>}\ln T)$, or $N_1^\nu=(D_{\perp\alpha}\ln T)(\nabla\cdot u)\Delta^{\alpha\nu}/3$ cannot be eliminated, then the formulas are not the most general second-order forms. Independent confirmation would also require solving the characteristic equations of the full nonlinear system in a curved background and verifying that all maximum propagation speeds remain real and subluminal.
Extended reading notes
Core claim
The central claim is that eqs. (17), (22), and (23) are, respectively, the general causal second-order constitutive relations for the shear viscosity tensor, bulk viscosity, and heat flow vector in the Eckart frame. Each flux explicitly carries curvature terms through the Ricci tensor, scalar curvature, and Riemann-type projections, and each is accompanied by a relaxation-time term in the spirit of the Israel-Stewart formalism. The paper argues that linearizing the resulting Navier-Stokes equations around Minkowski spacetime yields dispersion relations with finite large-wavenumber group velocities given by eqs. (37) and (42), provided the relaxation times satisfy conditions such as $\tau_\pi > \eta/(e_0+p_0)$. In the absence of conserved charges and heat flow, the sound-mode speed reduces to known results.
Load-bearing premise
The paper assumes the lists in Tables I-III exhaust all independent second-order gradient terms and that the terms omitted from the final formulas are redundant or absorbable, without showing a reduction; it also assumes that linearized causality results obtained around Minkowski spacetime carry over to the full curved-spacetime theory.
Editorial extensions
If this is right
- If eqs. (17), (22), and (23) are the general Eckart-frame forms, then second-order dissipative relativistic fluids with a conserved charge and heat flow can be modeled without giving up causality or stability in the linear regime.
- The explicit curvature terms mean bulk viscosity, shear viscosity, and heat flow respond directly to local spacetime geometry, not just to fluid gradients, which matters in strong-gravity environments.
- The finite maximum speeds (37) and (42) furnish concrete bounds on relaxation times, e.g., $\tau_\pi > \eta/(e_0+p_0)$, that any microscopic theory must respect if the macroscopic theory is to remain causal.
- The formalism gives a foundation for studying viscous thick accretion disks, modified hydrostatic equilibrium in neutron stars (with bulk-viscous curvature corrections to the effective pressure), and curvature-driven cosmological evolution in the Eckart frame.
- All first- and second-order transport coefficients remain in principle computable from an underlying microscopic theory, e.g., through Kubo formulas, so the general forms can be turned into predictive models once those coefficients are known.
Reading between the lines
- If the completeness of Tables I-III is later established rigorously, the same enumeration method should generalize to third-order gradient terms, and the flat-space causality calculation could be repeated on a curved background to test whether the same relaxation-time inequalities still guarantee subluminal propagation.
- The curvature terms $\kappa_1 R^{\langle\mu\nu\rangle}$ and $\kappa_2 u_\alpha u_\beta R^{\alpha\langle\mu\nu\rangle\beta}$ imply that a shear stress can be generated purely by spacetime curvature even when all fluid gradients vanish; this is a testable prediction for numerical relativity simulations of tidal encounters or gravitational-wave-driven fluid configurations.
- Because the Eckart frame freezes charge diffusion, the heat-flow sector here is tied to temperature and chemical-potential gradients; a natural extension would compare these results with Landau-frame calculations to isolate genuine frame dependence of the new second-order coefficients.
- The claimed curvature-driven bulk-viscous pressure at zero expansion ($\nabla\cdot u=0$, as in static stars) suggests that hydrostatic equilibrium configurations could shift slightly due to Ricci-tensor terms; computing the magnitude of this shift with realistic equations of state would provide a quantitative observational target.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a second-order gradient-expansion formulation of relativistic dissipative hydrodynamics in the Eckart frame. It claims to give the most general second-order constitutive relations for the shear viscosity tensor, bulk viscosity, and heat flow vector, eqs. (17), (22), and (23), including curvature couplings and relaxation-time terms. The paper then analyzes linearized perturbations around Minkowski spacetime to derive the shear and sound mode dispersion relations and concludes that the theory is causal because the propagation speeds are finite, eqs. (37) and (42). It also sketches astrophysical applications to thick accretion disks and static neutron stars. The central claims are the generality of the constitutive relations and the causality of the resulting theory in curved spacetime.
Significance. If the claims were established, the paper would provide a covariant second-order Eckart-frame hydrodynamics with explicit curvature terms, which would be useful for strong-gravity astrophysical applications such as accretion disks and neutron star interiors. However, the manuscript does not currently establish either of its two central claims. The derivation of the constitutive relations omits several terms from the paper's own tables without showing their dependence, and the shear-mode causality analysis rests on a dispersion relation that is algebraically inconsistent with the formula used to obtain the finite propagation speed. The transfer of a flat-space linearized analysis to the full curved-spacetime theory is also not justified. The paper contains no machine-checked proofs, reproducible code, or falsifiable predictions; its value at present is primarily programmatic.
major comments (4)
- [Sec. III, eqs. (33)-(35)] The central claim that eqs. (17), (22), and (23) are the general second-order constitutive relations in the Eckart frame is not supported, because several entries from the paper's own lists are omitted without any demonstration that they are dependent or reducible. For example, the scalar M1 = D^alpha_perp D^perp_alpha ln T from Table I does not appear in the bulk viscosity (22); the vector N1^nu = (D^alpha_perp ln T)(nabla.u)/3 Delta_alpha^nu from Table II is absent from the heat flux (23); and the tensor O1^mu nu = (D^<mu_perp ln T)(D^nu>_perp ln T) from Table III is absent from the shear tensor (17). For the shear tensor, eq. (19) eliminates O1 by inverting a basis relation, but the linear independence and invertibility of that basis are not shown. For bulk and heat, no reduction is attempted at all. The paper itself states in Sec. IV that determining how many of the second-order coefficients are independent is left for future work. If any omitted term is independent, the corresponding transport coefficient is missing and the claimed 'general form' is incomplete; this is a load-bearing gap in the paper's main result.
- [Sec. III, linearized analysis and curved spacetime] The shear-mode dispersion relation used for the causality claim is internally inconsistent. Eq. (34) is written as eta k^2 + i omega (1 - i omega tau_pi) = 0, which is equivalent to tau_pi omega^2 + i omega + eta k^2 = 0. For real wavenumber k, the solutions for omega are purely imaginary, so there is no propagating shear mode and no finite group velocity. The formula in eq. (35) contains the combination 4 k^2 tau_pi eta/(e0+p0), which would follow from an equation of the form tau_pi omega^2 + i omega - [eta/(e0+p0)] k^2 = 0, or its sign/e0+p0-normalized equivalent. The factor (e0+p0) is absent from eq. (34), and the sign of the k^2 term differs from the standard kinematic-viscosity term. Consequently, the expression for the maximum shear velocity in eq. (37), which is quoted as evidence of causality, does not follow from the dispersion relation actually derived. This is a concrete algebraic error in a load-bearing part of the paper's causality argument.
- [Sec. II, after eq. (21)] The causality proof is performed exclusively for linearized perturbations around Minkowski spacetime, with the specific coordinate choice delta g_mu nu (t,z) and a single Fourier mode. The constitutive relations, however, are general covariant expressions that explicitly include curvature terms, and the advertised applications target the strong-gravity regime. A causal theory in curved spacetime requires an analysis of the characteristic cone or hyperbolicity of the full nonlinear system, or at least an argument showing that the flat-space linear modes control the local characteristic speeds in an arbitrary background. No such justification is provided. The manuscript therefore does not establish that the proposed second-order equations are causal in the curved spacetimes used for the astrophysical applications in Sec. IV.
- [Sec. II, after eq. (21)] The derivation of the relaxation equation for the shear tensor relies on a change of basis whose coefficients c_i in eq. (18) are never defined or constrained. In particular, the combination <D sigma_mu nu> + T grad_lambda(u^lambda/(4T)) sigma_mu nu is expressed in the O_i basis with coefficients c_i, but the invertibility of this expansion and the identification tau_pi = lambda_1/(2 c_1 eta) are not derived. If c_1 = 0 or the combination is not linearly independent from the other O_i, the elimination step in eq. (19) is invalid. This is part of the same completeness issue as the omitted Table III terms and should be addressed explicitly in any revision.
minor comments (4)
- [Sec. II, text after eq. (23)] In the expression for delta q^mu, the term 'D u^mu t' contains an apparent typographical defect ('t' should not be present). Please correct this.
- [Sec. II, text after eq. (23)] The paper states that there are nine second-order transport coefficients for shear, but the list tau_pi, xi_2 through xi_8, kappa_1, and kappa_2 contains ten. The subsequent enumeration for heat flux also uses xi_i rather than chi_i for the coefficients of eq. (23).
- [Sec. III, eq. (36)] The symbol h appears in the group velocity formula without definition; it is presumably the enthalpy density e0+p0, but it should be defined explicitly.
- [Sec. III, eqs. (32)-(40)] Several correlator expressions, including eqs. (32) and (40), contain terms whose dimensions and signs are not checked in the manuscript. Given the error identified in eqs. (33)-(35), the entire linear-response computation should be re-examined for consistency.
Circularity Check
No significant circularity: the constitutive relations are assembled by explicit gradient expansion; omitted terms are completeness gaps, not circular reductions.
full rationale
Walking the derivation chain, I find no step in which a claimed output is identical by construction to an input, and no load-bearing self-citation. Equations (17), (22), and (23) are constructed by listing second-order scalars, vectors, and tensors in Tables I–III and applying the standard ideal-fluid order reduction (9)–(10) plus first-order substitutions such as σμν = −πμν/(2η) to cast the fluxes in relaxation form; this is an explicit algorithmic construction rather than a fit or a renamed input. The causality section computes linearized collective modes from those constitutive relations, so the finite speeds in (37) and (42) are consequences of the relaxation terms by design; that is a consistency check of the model, not a circular prediction. The self-citations ([6]/[37] and [32]) point to applications such as Ricci cosmology and an accretion-disk toy model and do not justify the central derivation. The genuine weaknesses are completeness/independence gaps: terms M1, N1, and O1 from Tables I–III are dropped or eliminated without showing reducibility, and Sec. IV concedes “It would also be worth investigating how many of these second order transport coefficients are really independent, following the formalism developed in [20, 36]. We leave these issues for future work.” Those are correctness/generality concerns, not circularity, so the paper receives score 0.
Assumptions & free parameters
free parameters (4)
- First-order transport coefficients η, ζ, κ
- Second-order shear coefficients τπ, ξ2-ξ8, κ1, κ2
- Second-order bulk coefficients τΠ, ζ2-ζ10
- Second-order heat coefficients τq, χ2-χ11
assumptions (4)
- domain assumption Gradient expansion: dissipative corrections are organized by powers of gradients and terms beyond second order are negligible.
- ad hoc to paper The second-order tensor basis in Tables I-III is complete and independent.
- domain assumption First-order solutions can be substituted into second-order terms (order reduction), e.g. σμν = -πμν/(2η), without changing the second-order accuracy.
- domain assumption Linearized perturbations around Minkowski spacetime are sufficient to establish causality of the full curved-spacetime theory.
Cite this review
Pith. "Pith review of Second order causal hydrodynamics in Eckart frame: using gradient expansion scheme." pith.science (2026). https://pith.science/paper/5FJHWGSF
@misc{pith2026190809462,
author = {Pith},
title = {Pith review of: Second order causal hydrodynamics in Eckart frame: using gradient expansion scheme},
year = {2026},
howpublished = {\url{https://pith.science/paper/5FJHWGSF}},
note = {Machine review of arXiv:1908.09462}
}
read the original abstract
In the present work, we develop a causal theory of relativistic non-ideal fluids up to the second order in the Eckart frame using gradient expansion scheme. Keeping the spirit of Mueller-Israel-Stewart formalism, the general forms of bulk viscosity, shear viscosity tensor and the heat flow vector are presented. Since each of the flux quantities explicitly carry curvature terms, we show that our formalism finds application in astrophysics in particular in the strong gravity regime. We elucidate two such applications namely in viscous thick accretion disks also known as Polish doughnuts and in addressing non-rotating equilibrium configuration, like neutron stars.
Reference graph
Works this paper leans on
-
[1]
Collective flow and viscosity in relativistic h eavy-ion collisions,
U. Heinz and R. Snellings, “Collective flow and viscosity in relativistic h eavy-ion collisions,” Ann. Rev. Nucl. Part. Sci. 63 (2013) 123 [arXiv:1301.2826 [nucl-th]]
arXiv 2013
-
[2]
Hydrodynamic Modeling of Heav y-Ion Collisions,
C. Gale, S. Jeon and B. Schenke, “Hydrodynamic Modeling of Heav y-Ion Collisions,” Int. J. Mod. Phys. A 28 (2013) 1340011 [arXiv:1301.5893 [nucl-th]]
arXiv 2013
-
[3]
Relativistic Fluid Dynamics In and Out of Equilibrium,
P. Romatschke and U. Romatschke, “Relativistic Fluid Dynamics In and Out of Equilibrium,” doi:10.1017/9781108651998 arXiv:1712.05815 [nucl-th]
-
[4]
W. Zimdahl, “Bulk viscous cosmology,” Phys. Rev. D 53, 5483 (1996) [astro-ph/9601189]
arXiv 1996
-
[5]
Causal thermodynamics in relativity,
R. Maartens, “Causal thermodynamics in relativity,” astro-ph/ 9609119
-
[6]
R. Baier, S. Lahiri and P. Romatschke, “Ricci cosmology,” arXiv:1 907.02974 [gr-qc]
-
[7]
The runaway instability of thick discs around black holes. I. The constant angular momentum case
J. A. Font and F. Daigne, “The Runaway instability of thick discs ar ound black holes. 1. The Constant angular momentum case,” Mon. Not. Roy. Astron. Soc. 334, 383 (2002) [astro-ph/0203403]
work page Pith review arXiv 2002
-
[8]
Euler, Principes generaux du mouvement des fluides, Mem
L. Euler, Principes generaux du mouvement des fluides, Mem. Aca d. Sci. Berlin 11 (1755) [printed in 1757]. Also in Opera omnia, ser. 2, 12, (1907) 54-91, E226
work page 1907
Show all 41 references
-
[9]
Navier, Memoire sur les lois du mouvement des fluides, Mem
C.L.M.H. Navier, Memoire sur les lois du mouvement des fluides, Mem. A cad. Sci. Inst. France, 6, (1822), 389-440
-
[10]
Stokes, On the theories of the internal friction of fluids in m otion, and of the equilibrium and motion of elastic solids, Trans
G.G. Stokes, On the theories of the internal friction of fluids in m otion, and of the equilibrium and motion of elastic solids, Trans. Camb. Philos. Soc. 8:287-319 (1845)
-
[11]
Eckart, Phys
C. Eckart, Phys. Rev. 58, 919 (1940)
1940
-
[12]
Landau, and E.M
Fluid Mechanics, Second Edition: Volume 6 (Course of Theoretica l Physics) L.D. Landau, and E.M. Lifshitz; Butterworth-Heinemann, 2 edition, (Jan 15, 1987)
1987
-
[13]
Generic instabilities in first order dis sipative relativistic fluids
W.A. Hiscock and L. Lindblom; “Generic instabilities in first order dis sipative relativistic fluids”, Phys. Rev. D 31 725 (1985)
1985
-
[14]
Muller, Z
I. Muller, Z. Phys. 198 (1967) 329
1967
-
[15]
Non-stationary Irreversible Thermodynamics: a C ausal Relativistic Theory
W. Israel, “Non-stationary Irreversible Thermodynamics: a C ausal Relativistic Theory”, Ann. Phys (N.Y.) 100 310 (1976); W. Israel, “Thermodynamics of relativistic sy stems”, Physica 106A 209 (1981)
1976
-
[16]
W. A. Hiscock and L. Lindblom, Ann. Phys. (N.Y.) 151, 466 (1983) ; Phys. Rev. D 31, 725 (1985); Phys. Rev. D 35, 3723 (1987); Phys. Lett. A 131, 509 (1988); Ph ys. Lett. A 131, 509 (1988)
1983
-
[17]
Relativistic viscous hydrodynamics, conformal invariance, and holography,
R. Baier, P. Romatschke, D. T. Son, A. O. Starinets and M. A. S tephanov, “Relativistic viscous hydrodynamics, conformal invariance, and holography,” JHEP 0804, 100 (2008)
2008
-
[18]
New Developments in Relativistic Viscous Hydro dynamics,
P. Romatschke, “New Developments in Relativistic Viscous Hydro dynamics,” Int. J. Mod. Phys. E 19 (2010) 1, [arXiv:0902.3663 [hep-ph]]. 19
2010 arXiv
-
[19]
Causal theories of dissipative relativistic fluid dyn amics for nuclear collisions,
A. Muronga, “Causal theories of dissipative relativistic fluid dyn amics for nuclear collisions,” Phys. Rev. C 69, 034903 (2004)
2004
-
[20]
Relativistic Viscous Fluid Dynamics and Non-Equilib rium Entropy,
P.Romatschke, “Relativistic Viscous Fluid Dynamics and Non-Equilib rium Entropy,” Class. Quant. Grav. 27, 025006 (2010)
2010
-
[21]
Second order dissipative fluid dynamics for ultra- relativistic nuclear collisions,
A. Muronga, “Second order dissipative fluid dynamics for ultra- relativistic nuclear collisions,” Phys. Rev. Lett. 88 (2002) 062302 Erratum: [Phys. Rev. Lett. 89 (2002) 159901], [nucl-th/0104064]
2002 arXiv
-
[22]
Fluid dyna mics of R-charged black holes,
J. Erdmenger, M. Haack, M. Kaminski and A. Yarom, “Fluid dyna mics of R-charged black holes,” JHEP 0901 (2009) 055, [arXiv:0809.2488 [hep-th]]
2009 arXiv
-
[23]
Kubo Formulae for Second-Ord er Hydrodynamic Coefficients,
G. D. Moore and K. A. Sohrabi, “Kubo Formulae for Second-Ord er Hydrodynamic Coefficients,” Phys. Rev. Lett. 106 (2011) 122302 [arXiv:1007.5333 [hep-ph]]
2011 arXiv
-
[24]
Transport coefficients of gluon plas ma,
A. Nakamura and S. Sakai, “Transport coefficients of gluon plas ma,” Phys. Rev. Lett. 94 (2005) 072305 [hep-lat/0406009]
2005 arXiv
-
[25]
Viscosity in strongly interacting quantum field theories from black hole physics,
P. Kovtun, D. T. Son and A. O. Starinets, “Viscosity in strongly interacting quantum field theories from black hole physics,” Phys. Rev. Lett. 94 (2005) 111601 [hep-th/0405231]
2005 arXiv
-
[26]
Determin ation of the Shear Viscosity Relax- ation Time at Weak and Strong Coupling,
G. S. Denicol, J. Noronha, H. Niemi and D. H. Rischke, “Determin ation of the Shear Viscosity Relax- ation Time at Weak and Strong Coupling,” J. Phys. G 38 (2011) 124177, [arXiv:1108.6230 [nucl-th]]
2011 arXiv
-
[27]
and the hyperbolic property of the MIS theory is extensi vely studied [28]. The gradient expansion scheme is based on the idea that hydro dynamics of non-ideal fluids can be constructed in a systematic way by considering gradients of hydrodynamical variables and the source gµν ...
-
[28]
Dissipative fluid dynamics in the 3 + 1 formalism
J. Peitz and S. Appl, “Dissipative fluid dynamics in the 3 + 1 formalism ” Classical Quantum Grav., 16, 979-989 (1999)
1999
-
[29]
In this image, the black hole can be seen to be surrounded by a d isk-like luminous structure known as the accretion disk
has published the first image of a rotating black hole loc ated at the centre of M87 galaxy. In this image, the black hole can be seen to be surrounded by a d isk-like luminous structure known as the accretion disk. A key mechanism involved within the disk is the inward mass tr...
-
[30]
3+1 formulation of non-ideal hydrodynam ics
J. Peitz and S. Appl, “3+1 formulation of non-ideal hydrodynam ics”, Mon. Not. R. Astron. Soc. 296, 231-244 (1998)
1998
-
[31]
First M87 Event Horizon Te lescope Re- sults. I. The Shadow of the Supermassive Black Hole,
K. Akiyama et al. [Event Horizon Telescope Collaboration], “First M87 Event Horizon Te lescope Re- sults. I. The Shadow of the Supermassive Black Hole,” Astrophys. J . 875 (2019) no.1, L1
2019
-
[32]
Shakura and R.A
N.I. Shakura and R.A. Sunyaev; Astronomy and Astrophysics, Vol. 24, p. 337 - 355
-
[33]
Advection dominated accretion flo ws in the Kerr metric: 1. Basic equations,
C. F. Gammie and R. Popham, “Advection dominated accretion flo ws in the Kerr metric: 1. Basic equations,” Astrophys. J. 498 (1998) 313, [astro-ph/9705117]
1998 arXiv
-
[34]
A toy model of viscous relativistic g eometrically thick disk in Schwarzschild geometry,
S. Lahiri and C. Lmmerzahl, “A toy model of viscous relativistic g eometrically thick disk in Schwarzschild geometry,” arXiv:1909.10381 [gr-qc].)
1909 arXiv
-
[35]
R. C. Tolman, Relativity, Thermodynamics, and Cosmology (Oxfo rd, 1934)
1934
-
[36]
J. R. Oppenheimer and G. M. Volkoff, On Massive neutron cores, Phys. Rev. 55, 374381 (1939)
1939
-
[37]
Landau and Eckart frames for relativistic fluids in nu clear collisions,
A. Monnai, “Landau and Eckart frames for relativistic fluids in nu clear collisions,” Phys. Rev. C 100 (2019) no.1, 014901 [arXiv:1904.11940 [nucl-th]]
2019 arXiv
-
[38]
Constraints on the second order trans port coefficients of an uncharged fluid,
S. Bhattacharyya, “Constraints on the second order trans port coefficients of an uncharged fluid,” JHEP 1207 (2012) 104 [arXiv:1201.4654 [hep-th]]
2012 arXiv
-
[39]
Ricci cosmology,
R. Baier, S. Lahiri and P. Romatschke, “Ricci cosmology,” arXiv :1907.02974 [gr-qc]
1907 arXiv
-
[40]
Quark-Hadron Phase Transitions in Vis cous Early Universe,
A. Tawfik and T. Harko, “Quark-Hadron Phase Transitions in Vis cous Early Universe,” Phys. Rev. D 85 (2012) 084032 [arXiv:1108.5697 [astro-ph.CO]]
2012 arXiv
-
[41]
Viscous Quark -Gluon Plasma in the Early Universe,
A. Tawfik, M. Wahba, H. Mansour and T. Harko, “Viscous Quark -Gluon Plasma in the Early Universe,” 20 Annalen Phys. 523 (2011) 194 [arXiv:1001.2814 [gr-qc]]
2011 arXiv
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.