Validity of relativistic hydrodynamics beyond local equilibrium
Pith reviewed 2026-05-18 20:53 UTC · model grok-4.3
The pith
Relativistic hydrodynamics remains valid far from equilibrium because it interpolates smoothly between free streaming and collective behavior.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The author claims that formal solutions of the Boltzmann moment equations in the relaxation time approximation decompose into a divergent gradient series and exponentially decaying non-perturbative modes that encode initial conditions. These non-perturbative contributions control causality and the divergence of the gradient series. In the 0+1D Bjorken scenario the exact evolution of non-equilibrium terms has the same structural form as the gradient expansion, differing only through modified transport coefficients that reflect both initial data and free-streaming dynamics. Extending the analysis to 3+1D shows that hydrodynamics remains effective not because the system is close to equilibrium,
What carries the argument
Decomposition of formal solutions of the Boltzmann moment equations in the relaxation time approximation into a divergent gradient series plus exponentially decaying non-perturbative modes that encode initial conditions.
If this is right
- In the 0+1D Bjorken scenario the exact evolution of non-equilibrium terms shares the structural form of the gradient expansion but with transport coefficients modified by initial data and free-streaming dynamics.
- The non-perturbative modes are required to understand causality and the divergence of the gradient series.
- In 3+1D hydrodynamics remains effective because it interpolates smoothly between free streaming and collective behavior.
- This interpolation supplies a natural explanation for the success of hydrodynamic modeling of quark-gluon plasma evolution far from equilibrium.
Where Pith is reading between the lines
- Initial-condition information carried by the non-perturbative modes could be folded into hydrodynamic simulations through adjusted transport coefficients rather than full kinetic evolution.
- The same decomposition might be applied to other kinetic-theory approximations to test whether smooth interpolation between streaming and flow is generic.
- Comparisons with full 3+1D kinetic simulations of heavy-ion collisions could quantify how far the modified coefficients improve agreement with data.
Load-bearing premise
The formal solutions of the Boltzmann moment equations in the relaxation time approximation can be cleanly decomposed into a divergent gradient series plus exponentially decaying non-perturbative modes that encode initial conditions.
What would settle it
A direct numerical solution of the Boltzmann equation in the relaxation time approximation for Bjorken flow that fails to reproduce the claimed structural match between exact non-equilibrium evolution and the modified gradient expansion, or that shows non-exponential decay of the additional modes, would falsify the decomposition and its implications.
read the original abstract
We examine the applicability of relativistic hydrodynamics far from equilibrium by constructing formal solutions of the Boltzmann moment equations in the relaxation time approximation. These solutions naturally decompose into a divergent gradient series and exponentially decaying non-perturbative modes that encode initial conditions. The non-perturbative contributions are essential for understanding causality, the divergence of the gradient series, and the unexpected effectiveness of relativistic hydrodynamics far from equilibrium. In the 0+1D Bjorken scenario, we demonstrate that the exact evolution of non-equilibrium terms shares the same structural form as the gradient expansion, differing only through modified transport coefficients that reflect both initial data and free-streaming dynamics. Extending to 3+1D, we find that hydrodynamics remains effective not because the system is close to equilibrium, but because it interpolates smoothly between free streaming and collective behavior. This perspective offers a natural explanation for the remarkable success of hydrodynamics in modeling quark-gluon plasma evolution.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript constructs formal solutions of the Boltzmann moment equations in the relaxation time approximation (RTA). These solutions decompose into a divergent gradient series plus exponentially decaying non-perturbative modes that encode initial conditions. In the 0+1D Bjorken scenario the exact evolution of non-equilibrium terms is shown to share the same structural form as the gradient expansion, differing only through modified transport coefficients that incorporate initial data and free-streaming dynamics. The work extends the analysis to 3+1D and concludes that relativistic hydrodynamics remains effective far from local equilibrium because it interpolates smoothly between free-streaming and collective regimes, thereby explaining its success in modeling quark-gluon plasma evolution.
Significance. If the decomposition and its implications hold, the paper supplies a concrete mathematical framework that reframes the applicability of hydrodynamics beyond local equilibrium. The explicit 0+1D construction and the identification of non-perturbative modes that control causality and the divergence of the gradient series constitute a genuine advance. The perspective that hydrodynamics succeeds by interpolating regimes rather than by proximity to equilibrium has clear relevance for heavy-ion phenomenology.
major comments (2)
- [3+1D extension and concluding discussion] The central claim that hydrodynamics interpolates smoothly between free streaming and collective behavior rests on the clean separation of the formal solutions into a divergent gradient series and exponentially decaying modes. This separation is demonstrated within the linear, memoryless RTA collision term; the manuscript does not examine whether the same structure persists for nonlinear or nonlocal collision kernels relevant to the QGP. Because the broad validity statement beyond local equilibrium depends on this separation, an explicit discussion or counter-example in a more general kinetic setting is required.
- [0+1D Bjorken scenario] In the 0+1D Bjorken analysis the exact solution is stated to share the structural form of the gradient expansion with only modified transport coefficients. Without the explicit derivation of those coefficients, error estimates, or direct comparison to known hydrodynamic limits, it is not possible to confirm that the non-perturbative modes are correctly isolated and that the claimed structural equivalence is not an artifact of the RTA.
minor comments (1)
- [Abstract] The abstract and introduction would benefit from a brief statement of the precise moment equations being solved and the definition of the relaxation time used throughout.
Simulated Author's Rebuttal
We thank the referee for the positive assessment of our work and for the constructive major comments. We address each point below and outline the revisions we will make.
read point-by-point responses
-
Referee: The central claim that hydrodynamics interpolates smoothly between free streaming and collective behavior rests on the clean separation of the formal solutions into a divergent gradient series and exponentially decaying modes. This separation is demonstrated within the linear, memoryless RTA collision term; the manuscript does not examine whether the same structure persists for nonlinear or nonlocal collision kernels relevant to the QGP. Because the broad validity statement beyond local equilibrium depends on this separation, an explicit discussion or counter-example in a more general kinetic setting is required.
Authors: We agree that the decomposition is derived specifically within the linear RTA. The separation into gradient series plus non-perturbative modes follows directly from the formal solution of the moment equations under this collision term. In the revised manuscript we will add a dedicated paragraph in the concluding section that explicitly discusses the assumptions of the RTA, notes that nonlinear or nonlocal kernels would require a different treatment, and outlines how the non-perturbative modes could be expected to persist or be modified. A full counter-example in a general kinetic theory lies beyond the scope of the present study. revision: partial
-
Referee: In the 0+1D Bjorken analysis the exact solution is stated to share the structural form of the gradient expansion with only modified transport coefficients. Without the explicit derivation of those coefficients, error estimates, or direct comparison to known hydrodynamic limits, it is not possible to confirm that the non-perturbative modes are correctly isolated and that the claimed structural equivalence is not an artifact of the RTA.
Authors: The explicit construction of the modified transport coefficients is given in Section III, where the exact solution of the RTA moment equations for Bjorken flow is matched term-by-term to the gradient expansion. To make this clearer we will add error estimates for the truncation of the non-perturbative sector and direct numerical comparisons against the ideal hydrodynamic limit in the revised version. revision: yes
Circularity Check
No circularity: derivation proceeds from direct formal solutions of Boltzmann equations
full rationale
The paper constructs formal solutions of the Boltzmann moment equations in the relaxation time approximation and shows that these solutions decompose into a divergent gradient series plus exponentially decaying non-perturbative modes. This decomposition is presented as arising directly from the structure of the RTA equations themselves rather than from any fitted parameter, self-referential definition, or load-bearing self-citation. No equation is shown to equal another by construction, no prediction is obtained by renaming a fit, and the central claims about hydrodynamics interpolating free streaming and collective behavior follow from the explicit solutions in 0+1D and 3+1D without reducing to prior inputs. The derivation chain is therefore self-contained.
Axiom & Free-Parameter Ledger
free parameters (1)
- relaxation time
axioms (1)
- domain assumption The Boltzmann equation in the relaxation time approximation provides an adequate description of the microscopic dynamics.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
These solutions naturally decompose into a divergent gradient series and exponentially decaying non-perturbative modes... hydrodynamics remains effective... because it interpolates smoothly between free streaming and collective behavior.
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leancostAlphaLog_fourth_deriv_at_zero unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
δ df/dt = −(f(t)−g(t))... exact solution f(t)=e^{−t/δ}f(0)+∫... γ(n+1,w) expansion
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
-
[1]
[f ′f ′ 1 − f f1] (11) where W (p, p1 | p′, p′
-
[2]
The term in brackets, [ f ′f ′ 1 − f f1] rep- resents gain and loss due to collisions
is the Lorentz invariant transi- tion rate for the process p + p1 → p′ + p′ 1 and dP = d3p/(2π)32p0. The term in brackets, [ f ′f ′ 1 − f f1] rep- resents gain and loss due to collisions. The vanishing of the collision kernel defines local thermal equilibrium distribution C[feq, feq] = 0, feq = 1 e(u·p−µ)/T + r , (12) where r takes values −1, 0, 1 corresp...
-
[3]
f(p1) . (28) and L[f] = Z dP ′dP ′ 1dP1 W (p′, p′ 1 | p, p1) f ′f ′ 1 (29) A direct integration of (27) using the method of char- acteristics yields a solution [30] (subject to appropriate boundary conditions), f(pµ, xµ) =e− 1 δ R t t0 v(t,t′)f(t, t0) + Z t t0 v(t, t′) δ e− 1 δ R t t′ v(t,t′′)dt′′ g(t, t′)dt′ (30) 5 where g = L/v. Here we have introduced ...
-
[4]
Solutions to the moment equation Consider the free streaming equation, h ∂τ + ˆF i ⃗ ρ= 0, (78) for which we can write down a propagator which satisfies, ∂τ ˆFτ,τ0 = −ˆF ˆFτ,τ0 . (79) The free streaming solution can be written as, ⃗ ρF (τ) = ˆFτ,τ0 ⃗ ρ(τ0) (80) Now consider the damped equation, h ∂τ + ˆF + ˆνR i ⃗ ρ= 0 . (81) To solve (81), we treat the d...
-
[5]
Exact solutions to the moment equation We rewrite the moment equations as " ∂τ + 2l + 1 τ ˆI − 2l − n τ ˆS + ˆI τR # ρn(τ, l) = ˆI τR ρeq n (τ, l) (B1) where we have defined the shift operator, ˆS [h(l)] = h(l + 1) . (B2) This essentially encodes the linear matrix algebra of the set of coupled differential equations. But, this allows us to use formal alge...
-
[6]
0 + 1D Gradient expansion Consider the hydrodynamic generator, ρG n,l(τ) = Z τ τ0 dτ ′ 1 τR e−ξ′ e− ˆKτ,τ ′ ρeq n,l(τ ′) (B23) where ξ′ = R τ τ ′ dτ ′ 1 τR . The free streaming propagator Kτ,τ ′ = e− ˆKτ,τ ′ (B24) satisfies ∂τ ′Kτ,τ ′ρn,l(τ) = Kτ,τ ′ ∂τ ′ + ˆF ρn,l(τ ′) (B25) Consider the change of variables, τ ′ → ξ′. We now expand e− ˆF(0,ξ′)ρeq n,l(ξ′)...
-
[7]
Hydrodynamics Hydrodynamics in 0 + 1D requires us to find the dynamical equations for the energy density ϵ and the two anisotropies, bulk pressure Π and shear pressure π. From energy conservation we get the dynamical equation for energy density (ϵ = ρ1,1), ˆDρ1,0 = ˆDn = 0 (B30) ˆDρ1,1 = ˆDϵ = 0 (B31) which gives us the familiar energy density and number ...
-
[8]
Collision picture We will choose to decompose the propagator into free streaming and collision part and treat the collision as the interaction. Consider the free streaming equation, h ∂τ + ˆL i = 0 (C2) with the propagator ˆF satisfying, ∂τ ˆFτ,τ0 = −ˆL ˆFτ,τ0 (C3) We can now define the collision picture operator and variables, ⃗ ρc(τ) = ˆF −1 τ,τ0 ⃗ ρ(τ)...
-
[9]
Consider the damping equation, [∂τ + ˆνR] ⃗ ρ(τ) = 0
F ree Streaming picture To derive the gradient expansion in the 3 + 1D case, we adopt an alternative approach in which free streaming is treated as a perturbation. Consider the damping equation, [∂τ + ˆνR] ⃗ ρ(τ) = 0 . (C15) 20 with the associated propagator obeying, ∂τ Wτ,τ0 = −ˆνR Wτ,τ0 . (C16) The propagator has the simple form, Wτ,τ0 = e− R τ τ0 dτ ′′...
-
[10]
R. B. Clare and D. Strottman, Physics Reports 141, 177 (1986)
work page 1986
- [11]
-
[12]
P. F. Kolb and U. Heinz, in Quark Gluon Plasma 3 , edited by R. C. Hwa and X. N. Wang (World Scientific, 2003) pp. 634–714, also available as arXiv preprint nucl-th/0305084
work page internal anchor Pith review Pith/arXiv arXiv 2003
-
[13]
D. A. Teaney, Viscous Hydrodynamics and the Quark Gluon Plasma, in Quark-gluon plasma 4 , edited by R. C. Hwa and X.-N. Wang (2010) pp. 207–266, arXiv:0905.2433 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2010
- [14]
-
[15]
V. Khachatryan et al. (CMS), Phys. Rev. Lett. 116, 172302 (2016), arXiv:1510.03068 [nucl-ex]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[16]
G. Aad et al. (ATLAS), Phys. Rev. Lett. 116, 172301 (2016), arXiv:1509.04776 [hep-ex]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[17]
Fluid Dynamics far from Local Equilibrium
P. Romatschke, Phys. Rev. Lett. 120, 012301 (2018), arXiv:1704.08699 [hep-th]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[18]
M. P. Heller, R. A. Janik, and P. Witaszczyk, Phys. Rev. Lett. 108, 201602 (2012), arXiv:1103.3452 [hep-th]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[19]
M. P. Heller and M. Spalinski, Phys. Rev. Lett. 115, 072501 (2015), arXiv:1503.07514 [hep-th]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[20]
A. Buchel, M. P. Heller, and J. Noronha, Phys. Rev. D 94, 106011 (2016), arXiv:1603.05344 [hep-th]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[21]
M. P. Heller, R. A. Janik, and P. Witaszczyk, Phys. Rev. Lett. 110, 211602 (2013)
work page 2013
-
[22]
A. Dash and V. Roy, Phys. Lett. B 806, 135481 (2020), arXiv:2001.10756 [nucl-th]
-
[23]
Z. Chen and L. Yan, Phys. Rev. C 105, 024910 (2022), arXiv:2109.06658 [nucl-th]
- [24]
-
[25]
S. Jaiswal, C. Chattopadhyay, L. Du, U. Heinz, and S. Pal, Phys. Rev. C 105, 024911 (2022), arXiv:2107.10248 [hep-ph]
-
[26]
R. Gangadharan, A. K. Panda, and V. Roy, Phys. Rev. D 109, 074020 (2024), arXiv:2311.10547 [nucl-th]
-
[27]
Fluid dynamics of out of equilibrium boost invariant plasmas
J.-P. Blaizot and L. Yan, Phys. Lett. B 780, 283 (2018), arXiv:1712.03856 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[28]
J.-P. Blaizot and L. Yan, Annals Phys. 412, 167993 (2020), arXiv:1904.08677 [nucl-th]
-
[29]
J.-P. Blaizot and L. Yan, Phys. Lett. B 820, 136478 (2021), arXiv:2006.08815 [nucl-th]
-
[30]
J.-P. Blaizot and L. Yan, Phys. Rev. C 104, 055201 (2021), arXiv:2106.10508 [nucl-th]
-
[31]
S. Jaiswal, J.-P. Blaizot, R. S. Bhalerao, Z. Chen, A. Jaiswal, and L. Yan, Phys. Rev. C 106, 044912 (2022), arXiv:2208.02750 [nucl-th]
-
[32]
Transport coefficients for bulk viscous evolution in the relaxation time approximation
A. Jaiswal, R. Ryblewski, and M. Strickland, Phys. Rev. C 90, 044908 (2014), arXiv:1407.7231 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[33]
W. A. Hiscock and L. Lindblom, Phys. Rev. D 35, 3723 (1987)
work page 1987
-
[34]
W. A. Hiscock, Phys. Rev. D 33, 1527 (1986)
work page 1986
-
[35]
G. S. Denicol, T. Kodama, T. Koide, and P. Mota, J. Phys. G 35, 115102 (2008), arXiv:0807.3120 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[36]
Israel, Annals of Physics 100, 310 (1976)
W. Israel, Annals of Physics 100, 310 (1976)
work page 1976
-
[37]
J. M. Stewart, Proceedings of the Royal Society of London A 10.1098/rspa.1977.0155
-
[38]
C. Cercignani and G. M. Kremer, The Relativistic Boltzmann Equation: Theory and Applications (Birkh¨ auser, 2002)
work page 2002
-
[39]
H. Grad, in Handbuch der Physik / Encyclopedia of Physics, Thermodynamik der Gase / Thermodynamics of Gases , Vol. 3/12, edited by S. Fl¨ ugge (Springer, Berlin, Heidelberg, 1958) pp. 205–294
work page 1958
-
[40]
S. Chapman and T. G. Cowling, The Mathematical Theory of Non-uniform Gases: An Account of the Kinetic Theory of Viscosity, Thermal Conduction and Diffusion in Gases , 3rd ed. (Cambridge University Press, Cambridge, 1990)
work page 1990
-
[41]
J. L. Anderson and H. R. Witting, Physica 74, 466 (1974)
work page 1974
-
[42]
G. S. Denicol, T. Koide, and D. H. Rischke, Phys. Rev. Lett. 105, 162501 (2010), arXiv:1004.5013 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[43]
G. S. Denicol, H. Niemi, E. Molnar, and D. H. Rischke, Phys. Rev. D 85, 114047 (2012), [Erratum: Phys.Rev.D 91, 039902 (2015)], arXiv:1202.4551 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[44]
Relativistic dissipative hydrodynamics from kinetic theory with relaxation time approximation
A. Jaiswal, Phys. Rev. C 87, 051901 (2013), arXiv:1302.6311 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[45]
Relativistic third-order dissipative fluid dynamics from kinetic theory
A. Jaiswal, Phys. Rev. C 88, 021903 (2013), arXiv:1305.3480 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2013
- [46]
- [47]
- [48]
-
[49]
G. S. Denicol and J. Noronha, (2016), arXiv:1608.07869 [nucl-th]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[50]
R. Gangadharan and V. Roy, Phys. Rev. D 111, L071901 (2025), arXiv:2405.10846 [nucl-th]
-
[51]
L. Gavassino, M. M. Disconzi, and J. Noronha, Phys. Rev. Lett. 132, 162301 (2024), arXiv:2307.05987 [hep-th]
-
[52]
J. D. Bjorken, Phys. Rev. D 27, 140 (1983)
work page 1983
- [53]
-
[54]
R. Gangadharan, S. Mitra, and V. Roy, (2025), arXiv:2504.00126 [nucl-th]
-
[55]
M. McNelis and U. Heinz, Phys. Rev. C 101, 054901 (2020), arXiv:2001.09125 [nucl-th]
-
[56]
M. Strickland and U. Tantary, JHEP 10, 069, arXiv:1903.03145 [hep-ph]
- [57]
- [58]
- [59]
- [60]
- [61]
- [62]
- [63]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.