{"id":"b3759941-108b-47c1-93c1-47e4362433bb","arxiv_id":"1909.02498","paper_version":3,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A new family of exact analytic viscous hydrodynamic solutions is derived for Hubble-flow fireballs, in which shear viscosity cancels and only bulk viscosity affects the temperature evolution.","lead":"This paper presents exact analytic solutions of relativistic viscous hydrodynamics for expanding, ellipsoidal fireballs with a Hubble-type flow. The solutions isolate the effect of bulk viscosity, because shear viscosity and heat conduction drop out, which could help constrain bulk viscosity from heavy-ion collision data.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The physical bulk-viscosity solutions are not actually ellipsoidal: in Cases A, C, D, F the temperature and energy density depend only on τ, and the only S-dependent case (B) is unphysical; the title and abstract overstate the ellipsoidal-symmetry claim.","rationale":"The algebraic derivation from Eq. (26) to Eqs. (30)-(31) and the ODE solutions are internally consistent; I verified the parallel projection and the Case A/B particular solution. The Case F printed formula (54) contains a likely typo (Γ(B,τ/τΠ) should be Γ(B,τ0/τΠ)), but this is minor and easily repaired. The load-bearing issue is the ellipsoidal-symmetry claim. The definition of ellipsoidal symmetry enters through the scaling variable S and the conserved-density function V(S) in Eq. (25). In the physical scenarios D and E the authors explicitly set V(S) = 1 (Section III); scenarios A, C and F have no conserved density and therefore no V(S). Only scenario B supports an arbitrary V(S), and the authors themselves state its temperature diverges as τ^{d-1}. Hence the exact solutions used to illustrate bulk-viscosity effects are spherically symmetric functions of the Hubble proper time, not ellipsoidal temperature and density profiles. The orthogonality condition Eq. (29) shows why: any S-dependence of ζ must be canceled by S-dependence of p, which the ansatz p = p(τ) forbids. This does not invalidate the mathematical solutions, but it materially narrows the central claim: the paper delivers exact bulk-viscous Hubble-flow solutions, but not the advertised ellipsoidal fireballs for the physical cases. I therefore recommend conditional acceptance with a revision that either supplies a physical ellipsoidal example or qualifies the title and abstract accordingly, and corrects Eq. (54).","tokens_in":13745,"tokens_out":34587,"duration_ms":359819,"concrete_test":"Re-derive the orthogonality condition Eq. (29) for Case D without assuming V(S) = 1, taking n = n0(τ0/τ)^d V(S) and ζ = ζ0 n/n0. If the transverse part of (d/τ)∂νζ − ∂νp does not vanish for nonconstant V(S), then no ellipsoidal density or temperature profile is compatible with the presented physical bulk-viscosity cases; this settles whether the title's ellipsoidal-symmetry claim is actually realized.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central claim, as stated in the title and abstract, is exact viscous solutions for fireballs with ellipsoidal symmetry. But the explicit physical solutions (Cases A, C, D, F) have p = p(τ) and T = T(τ); the only possible S-dependence is the arbitrary V(S) in Eq. (25) for the conserved density. In Cases D and E the authors impose V(S) = 1 (Section III after Eq. (25)); Cases A and C have no conserved n, so no V(S) exists; Case F also has no conserved density. Only Case B retains an arbitrary V(S), and it is the case the authors themselves reject as unphysical because T ∝ τ^{d-1} diverges. Consequently, none of the physical bulk-viscosity solutions exhibits a nontrivial ellipsoidal temperature or energy-density profile. The orthogonality condition Eq. (29), which requires p − dζ/τ to be a function of τ alone, is what eliminates S-dependence: whenever ζ depends on S through n or T, V(S) is forced to a constant. Thus the ellipsoidal-symmetry aspect of the central claim is not delivered by the scenarios used for the bulk-viscosity study; the solutions are effectively spherically symmetric Hubble-flow solutions with bulk viscosity.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper constructs exact analytic solutions of relativistic viscous hydrodynamics for the Hubble-type velocity field u^mu = x^mu/tau. For this profile the shear tensor vanishes identically and, when thermal conductivity is neglected, energy-momentum conservation reduces to the ordinary differential equations (36)-(38). The authors solve these equations for six scenarios for the bulk viscosity (constant, proportional to s, n, T^kappa, or Pi) and for two equations of state (p = nT or epsilon = kappa p), in both the Navier-Stokes and one Israel-Stewart case, and they present explicit p(tau) and T(tau) formulas together with illustrative plots. They find that bulk viscosity slows the cooling of the fireball relative to the perfect-fluid Hubble flow, and they identify Case B and, conditionally, Case E as unphysical.","tokens_in":13998,"tokens_out":18153,"duration_ms":169154,"significance":"The central ODE reduction is sound, and the resulting explicit solutions are useful as benchmarks for numerical viscous hydrodynamics, particularly because the shear viscosity cancels identically for the Hubble profile and the thermal conductivity also cancels under the condition T = T(tau). Exact analytic viscous solutions with bulk viscosity are rare, so this is a genuine contribution. The paper is also honest in flagging the unphysical cases. However, the advertised ellipsoidal symmetry is not realized in the physical solutions, which reduces the significance relative to the title's promise.","major_comments":[{"comment":"The central claim that the solutions describe fireballs with ellipsoidal symmetry is not supported by the explicit solutions. In the physical Cases A, C, D, and F, all thermodynamic quantities are functions of tau alone: p(tau) and T(tau) are given by Eqs. (41), (43), (46), and (53), with no dependence on the scaling variable S of Eq. (21). The only solution that can retain an arbitrary V(S) in Eq. (25) is Case B, which the authors themselves reject as unphysical because T diverges as tau^(d-1) (Section IV, after Eq. (41), and Section V). In Cases D and E the paper sets V(S)=1 (Section III, after Eq. (25)), and Cases A, C, and F have no conserved density n, so no V(S) exists. The orthogonality condition Eq. (29) forces p - d zeta/tau to be a function of tau alone, and this is what eliminates S-dependence whenever zeta depends on T or n. Consequently, the title and abstract overstate the ellipsoidal-symmetry content of the new solutions; the physical viscous solutions presented are spherically symmetric Hubble-flow solutions with bulk viscosity. Please either construct a genuine ellipsoidal family (for example by allowing S-dependent p in the reduction leading to Eq. (36)) or revise the title, abstract, and Section VI to describe the solutions as Hubble-flow, spherically symmetric solutions.","section":"Title, abstract, and Section VI (Summary); see also Eqs. (21), (25), (29), (36) and Table I"}],"minor_comments":[{"comment":"Equation (54a) contains a likely typo: the right-hand side contains Gamma(B, tau/tau_Pi), which makes p_A tau-dependent and inconsistent with Eq. (53). From Eq. (53) and the boundary condition p(tau_0)=p_0, the argument should be Gamma(B, tau_0/tau_Pi).","section":"Eq. (54)"},{"comment":"The exponent of (tau_0/tau) in Eq. (41) is typeset ambiguously as \"(d kappa+1)/kappa\"; it should be d(kappa+1)/kappa to match the homogeneous solution of Eq. (36). Please clarify.","section":"Eq. (41)"},{"comment":"The sentence stating that the solutions are valid for any type of thermal conductivity is too broad without the qualification given two paragraphs earlier: the cancellation requires T = T(tau), which in Cases B, D, and E forces V(S)=1. This condition should be stated in the summary sentence as well.","section":"End of Section IV"},{"comment":"The captions of Figures 1-4 do not state that the plotted solutions have no nontrivial S-dependence (V(S)=1 or no conserved density). Adding this information would make the ellipsoidal-symmetry limitation visible to the reader.","section":"Figures 1-4 and Table I"},{"comment":"There is a spelling error: \"reseach directions\" should be \"research directions\".","section":"Section VI"}],"recommendation":"major_revision","confidential_remarks":"The mathematical derivations appear sound, and the paper can be made publishable by correcting the overstatement about ellipsoidal symmetry. The typo in Eq. (54) should also be fixed. If the authors reframe the paper as exact viscous Hubble-flow solutions rather than ellipsoidal fireball solutions, the contribution is appropriate for the journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. The central reduction is correct: with u^μ = x^μ/τ, the shear tensor vanishes identically, and the projection of energy-momentum conservation gives Eqs. (28)-(31). The ODEs (36)-(38) follow, and I spot-checked the ζ = const solution (41): it satisfies the ODE. The explicit solutions for Cases A-F are new, and the observation that shear viscosity cancels for Hubble flow is clean and useful. As closed-form checks for numerical viscous hydro codes, this is real value.\n\nNow the soft spot. The title and abstract advertise ellipsoidal symmetry, but the physical solutions do not have it. The orthogonality condition (29) requires p − dζ/τ to be a function of τ alone. Cases A, C, F have no conserved density, and their T and ε come out as functions of τ only. In Cases D and E, the previously arbitrary V(S) in the density solution (25) is forced to V(S) = 1. Only Case B retains an arbitrary V(S), and Case B is unphysical (T diverges). So no physical solution has non-trivial S-dependence in the thermodynamics. The ellipsoids are there in the variable S, but not in the solutions actually presented. The perfect-fluid limit is the spherically symmetric Hubble flow, not the ellipsoidal CCHK solutions cited as the starting point. This is an overstatement in the abstract rather than a mathematical error, but it is significant: the paper promises something the solutions do not deliver.\n\nMinor items: Eq. (54) almost certainly has a typo — Γ(B, τ/τ_Π) should read Γ(B, τ₀/τ_Π) for p_A to be a constant; and the abstract's suggestion that one can infer bulk viscosity from measurements goes beyond what a simple, symmetric, acceleration-free Hubble solution can support.\n\nBottom line: solid technical work, worth refereeing. In revision, the authors should either soften the ellipsoidal language or exhibit an ellipsoidal viscous case. Good for the exact-solutions community and for code benchmarks; the headline claim just needs to be brought in line with the content.","headline":"Solid exact viscous Hubble-flow solutions, but the advertised ellipsoidal symmetry is not realized in any physical solution — the thermodynamics ends up depending on τ alone.","tokens_in":14532,"tokens_out":8191,"would_cite":true,"duration_ms":74730,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["20.24","20.25"],"model":"deepseek-v4-flash","headline":"New exact analytic solutions of relativistic viscous hydrodynamics show that for Hubble-type ellipsoidal fireballs the shear viscosity cancels completely, leaving bulk viscosity as the only dissipative influence on the temperature…","keywords":["relativistic viscous hydrodynamics","Hubble flow","bulk viscosity","exact analytic solutions","ellipsoidal symmetry","heavy-ion collisions","shear viscosity cancellation","Israel-Stewart hydrodynamics"],"falsifier":"Measure or simulate the velocity field of an expanding heavy-ion fireball and compute the shear tensor $\\sigma^{\\mu\\nu}$: the paper's cancellation of shear viscosity holds only where $\\sigma^{\\mu\\nu}=0$ exactly, so a single spacetime region with non-vanishing shear would break the claimed applicability there.","tokens_in":13544,"feed_emoji":"🔥","tokens_out":8278,"duration_ms":75148,"temperature":0.7,"pith_summary":"The paper constructs new exact, analytic solutions of relativistic viscous hydrodynamics for expanding fireballs whose velocity is the Hubble-type profile $u^\\mu = x^\\mu/\\tau$ and whose density and temperature profiles have ellipsoidal symmetry. In these solutions the shear-viscosity terms vanish identically, and when the temperature depends only on the proper time $\\tau$, the heat-conduction terms cancel as well, so bulk viscosity is the only dissipative mechanism left. The paper reduces the full partial differential equations to ordinary differential equations for the pressure and solves them for six concrete assumptions about how the bulk viscosity coefficient depends on temperature, entropy density, particle density, or bulk pressure. The physical solutions all show that bulk viscosity slows the cooling relative to a perfect fluid, and different scenarios for $\\zeta$ produce qualitatively different temperature histories, raising the possibility that measured energy-density and temperature evolution in heavy-ion collisions could constrain the bulk viscosity.","feed_headline":"Shear vanishes: exact viscous solutions isolate bulk viscosity","feed_subtitle":"For Hubble-type fireballs, temperature history depends only on bulk viscosity, making it measurable.","key_machinery":"The central object is the Hubble-type four-velocity profile $u^\\mu = x^\\mu/\\tau$, with proper time $\\tau=\\sqrt{t^2-r_x^2-r_y^2-r_z^2}$. Its defining identities are zero acceleration, $u^\\nu\\partial_\\nu u^\\mu=0$, zero shear, $\\sigma^{\\mu\\nu}=0$, and $u^\\mu\\partial_\\mu S=0$ for the ellipsoidal scaling variable $S=r_x^2/X^2+r_y^2/Y^2+r_z^2/Z^2$. These identities collapse the divergences of the viscous stress tensor into equations that depend only on $\\tau$; the machinery of the argument is the projection of $\\partial_\\mu T^{\\mu\\nu}=0$ into parts parallel and pseudo-orthogonal to $x^\\mu$, reducing a multi-dimensional PDE problem to the ordinary differential equations (36)-(38).","core_discovery":"For a relativistic fireball with Hubble-like velocity field $u^\\mu=x^\\mu/\\tau$ and ellipsoidal symmetry, the energy-momentum conservation equations of first-order viscous hydrodynamics reduce to a single ordinary differential equation for the pressure, Eq. (36), in the Navier-Stokes case and to the pair (37)-(38) in the Israel-Stewart case. Because the shear tensor $\\sigma^{\\mu\\nu}$ vanishes identically for this velocity field, all shear-viscosity effects cancel; if the temperature depends only on the proper time $\\tau$, all heat-conduction terms cancel as well. The paper solves the reduced equations for six scenarios for the bulk viscosity coefficient, obtaining closed-form pressure and temperature histories. Four of these cases are thermodynamically physical, while two produce unphysical reheating or divergence and are rejected. In the physical cases bulk viscosity always slows the decrease of energy density and temperature relative to the perfect-fluid limit, and the shape of the cooling curve depends on the assumed scaling of $\\zeta$.","pith_inferences":["If a real heavy-ion system ever approaches Hubble flow, shear viscosity would be hidden from the temperature profile by geometry rather than by smallness of $\\eta$, so extracting bulk viscosity would require independent evidence that the flow is truly boost-invariant.","The technique of choosing a zero-shear velocity field to cancel shear viscosity could be carried over to other flow geometries, such as boosted or accelerating profiles, to construct analogues that isolate bulk viscosity in less symmetric settings.","The unphysical cases suggest a model-discrimination rule: a bulk-viscosity ansatz that grows too rapidly as the conserved density dilutes will produce runaway heating, so the observed absence of reheating in heavy-ion data can set a lower bound on how fast $\\zeta$ may grow with $T$ or $n$."],"forward_implications":["Bulk viscosity changes the cooling law of the fireball: constant $\\zeta$, $\\zeta\\propto s$, and $\\zeta\\propto n$ produce noticeably different temperature histories, so measured cooling curves can in principle distinguish these scenarios.","Shear viscosity drops out of these solutions for any $\\eta(T)$ or $\\eta/s$ function, making the solutions exact benchmarks for numerical viscous-hydro codes that must reproduce the same $\\tau$-dependent evolution.","Two of the six scenarios, Case B (constant $\\zeta$ with conserved particle number) and Case E ($\\zeta\\propto T^\\kappa$ with $p=nT$), lead to unphysical reheating or divergence and are therefore excluded within this flow family.","Because thermal conductivity also cancels whenever $T=T(\\tau)$, the solutions remain valid for arbitrary $\\lambda$ in cases with no conserved charge (or with $V(S)=1$), so they isolate bulk viscosity as the only active dissipative transport coefficient."],"supporting_citations":[{"why":"Supply the perfect-fluid ellipsoidal Hubble solutions that are the starting point for the viscous generalization.","marker":"[18, 19]"},{"why":"Introduce the Hubble-type Hwa-Bjorken flow profile whose zero-shear property carries the argument.","marker":"[12, 13]"},{"why":"Provide the Israel-Stewart relaxation equations used for the bulk-pressure dynamics in Case F.","marker":"[26]"},{"why":"Show that Hubble flow develops in high-energy heavy-ion collisions, motivating the velocity ansatz.","marker":"[43]"},{"why":"Show compatibility of Hubble-type solutions with data, supporting the physical relevance of the profile.","marker":"[44, 45]"},{"why":"Identifies the class of asymptotically perfect fluid solutions that Case F belongs to and supplies the entropy-production condition $\\Pi_0<0$.","marker":"[46]"}],"fun_headline_variants":["Shear drops out: exact Hubble-flow solutions expose bulk viscosity","Hubble flow cancels shear: new viscous solutions isolate bulk","Exact solutions: shear irrelevant, bulk viscosity determines cooling","Bulk viscosity measurable from fireball cooling: exact solutions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire construction rests on assuming that the expanding matter moves with a Hubble-type velocity profile whose shear and acceleration vanish; if a real fireball's flow deviates from this boost-invariant, acceleration-free form, the cancellation of shear viscosity and heat conduction fails and the solutions do not apply.","fun_headline_variants_meta":{"raw":{"variants":["Shear drops out: exact Hubble-flow solutions expose bulk viscosity","Hubble flow cancels shear: new viscous solutions isolate bulk","Exact solutions: shear irrelevant, bulk viscosity determines cooling","Bulk viscosity measurable from fireball cooling: exact solutions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001021,"raw_usage":{"total_tokens":4267,"prompt_tokens":868,"completion_tokens":3399,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":484,"completion_tokens_details":{"reasoning_tokens":3331}},"tokens_in":484,"tokens_out":3399,"duration_ms":22020,"temperature":1.0,"reasoning_tokens":3331,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:48:03.124267+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure or simulate the velocity field of an expanding heavy-ion fireball and compute the shear tensor $\\sigma^{\\mu\\nu}$: the paper's cancellation of shear viscosity holds only where $\\sigma^{\\mu\\nu}=0$ exactly, so a single spacetime region with non-vanishing shear would break the claimed applicability there.","supporting_citations":[{"cited_title":"On the formation of Hubble flow in Little Bangs","cited_arxiv_id":"nucl-th/0410036","evidence_quote":"Show that Hubble flow develops in high-energy heavy-ion collisions, motivating the velocity ansatz."},{"cited_title":"Israel, and J","cited_arxiv_id":null,"evidence_quote":"Provide the Israel-Stewart relaxation equations used for the bulk-pressure dynamics in Case F."}],"review_version":1}