{"id":"6d4ba788-33c4-4838-88a3-c2cb0ef55ee2","arxiv_id":"2412.02405","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Nonlinear causality constrains one-dimensional Bjorken-expanding viscous fluids to small inverse Reynolds numbers, giving minimum initial times of about 0.5 to 1 fm and maximum initial energy densities of about 5 to 35 GeV/fm3 in central heavy-ion collisions.","lead":"Relativistic dissipative hydrodynamics for a simplified one-dimensional expanding fluid can violate nonlinear causality when the fluid starts far from local equilibrium. Using known causality conditions and RHIC/LHC transverse-energy data, the paper derives minimum starting times and maximum initial energy densities for a causal hydrodynamic description.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (37) does not follow from Eqs. (34)-(36): with C_tau_pi=(2-ln2)/(2pi) and Nf=3 ideal gas the prefactor is about 2.8, not 1.6, shifting the reported tau0,min/e0,max by roughly a factor 1.7.","rationale":"The reader's weakest assumption was the one-dimensional boost-invariant reduction and its applicability to RHIC/LHC conditions. That is a legitimate model-dependence concern. However, the most load-bearing issue I find is internal: the paper's own equations, when combined, do not yield Eq. (37). The abstract prominently advertises specific numbers (about 1 fm and 5 GeV/fm^3 at RHIC; 0.7 fm and 20 GeV/fm^3 at LHC for the conformal case). If the prefactor in Eq. (37) is actually about 2.8 rather than 1.6, these numbers shift substantially (e.g., RHIC tau0,min becomes about 1.7 fm and e0,max about 3 GeV/fm^3). This is a concrete, checkable arithmetic inconsistency, not a matter of fitting to a different theory or changing geometric assumptions. The qualitative central claim - that large inverse Reynolds number violates nonlinear causality and hence initial conditions are constrained - is supported by the inequalities and numerical solutions, so the paper should not be rejected outright. But because the headline quantitative results are currently based on an inconsistent derivation, the article needs revision; the CONDITIONAL verdict remains appropriate. My concrete test would settle whether the prefactor should be 2.8 or whether the readout in Eq. (36) or the degrees of freedom in Eq. (35) were miscalculated, and it would allow the authors to correct the reported constraints.","tokens_in":18581,"tokens_out":21751,"duration_ms":220005,"concrete_test":"Recompute Eq. (37) from Sec. III C: set tau0,min/tau_pi = 5.3, tau_pi = C_tau_pi/T with C_tau_pi=(2-ln2)/(2*pi), e0 = E0/tau0, and e0 = 47.5*pi^2/30 * T^4. Derive the prefactor (5.3*C_tau_pi)^{4/3} * (47.5*pi^2/30)^{1/3} and check whether it equals 1.6. If it yields about 2.8, re-plot the conformal-EoS results in Fig. 5 and recompute the RHIC/LHC values of tau0,min and e0,max accordingly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative headline constraints rest on Eq. (37), tau0,min = 1.6 * E0^{-1/3}, which is stated to follow from Eqs. (34)-(36). Combining Eq. (36) (tau0,min/tau_pi = tau0,min T0,max / C_tau_pi about 5.3) with tau_pi = C_tau_pi/T gives tau0,min T0,max = 5.3 * C_tau_pi = 5.3 * (2-ln2)/(2*pi) ≈ 1.10. Using the Nf=3 ideal-gas EoS (Eq. 35), e0 = 47.5*pi^2/30 * T^4, and the Bjorken relation e0 = E0/tau0, one obtains tau0,min = (5.3*C_tau_pi)^{4/3} * (47.5*pi^2/30)^{1/3} * E0^{-1/3} ≈ 2.8 * E0^{-1/3}. This is not 1.6. For RHIC E0=5 GeV/fm^2, the revised values are tau0,min≈1.7 fm and e0,max≈3.0 GeV/fm^3, versus the reported 1 fm and 5 GeV/fm^3; for LHC E0=14, tau0≈1.2 fm and e0≈12 GeV/fm^3, versus 0.7 fm and 20 GeV/fm^3. The qualitative conclusion that nonlinear causality constrains initial conditions is unaffected, but the quantitative claims in the abstract are not justified by the stated equations. The mismatch indicates an arithmetic slip, likely in the conversion from tau0/tau_pi to the final formula, and must be corrected before the numbers are used.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript analyzes the nonlinear causality conditions of Bemfica et al. (Ref. [20]) for one-dimensional boost-invariant expansion described by the BRSSS constitutive equation. After transcribing the necessary and sufficient conditions into inequalities in terms of the shear pressure φ and thermodynamic variables, the authors determine the allowed range of the inverse Reynolds number, map the causal/acausal regions in the initial-condition plane, and combine these restrictions with the measured transverse energy per rapidity at RHIC and LHC to obtain minimum initial proper times and maximum initial energy densities. For the conformal EoS with N=4 SYM transport coefficients, the reported values are τ0,min ≈ 1 fm and e0,max ≈ 5 GeV/fm^3 at RHIC and 0.7 fm and 20 GeV/fm^3 at LHC; for the lattice EoS, the corresponding numbers are about 0.5 fm and 11 GeV/fm^3 (RHIC) and 0.4 fm and 35 GeV/fm^3 (LHC). The central claim is that large inverse Reynolds numbers violate nonlinear causality, so the initial stage of heavy-ion collisions may require a non-equilibrium description.","tokens_in":19014,"tokens_out":11547,"duration_ms":109609,"significance":"If the central claim holds, the paper offers a new, parameter-free (in the sense of no fitting to data) constraint on hydrodynamic initial conditions that is complementary to Bayesian extraction and to attractor studies. The derivation is transparent: the causality inequalities are taken from an external group, the 1D equations are evolved explicitly, and the conformal-limit results are cross-checked against Ref. [20]. The use of measured dET/dy is a concrete falsifiable step. However, the quantitative headline numbers are affected by an arithmetic error in Sec. III C, and the heavy-ion application rests on the 1D Bjorken approximation; these issues require correction before the numbers can be used.","major_comments":[{"comment":"The prefactor in Eq. (37) does not follow from Eqs. (34)–(36). Combining Eq. (36) (τ0,minT0,max = 5.3 Cτπ) with Eq. (35) (e = 47.5 π^2/30 T^4) and the Bjorken relation e0 = E0/τ0 gives τ0,min = (5.3 Cτπ)^{4/3} (47.5 π^2/30)^{1/3} E0^{-1/3} ≈ 2.8 E0^{-1/3} fm, with Cτπ = (2−ln2)/(2π). The reported coefficient 1.6 would correspond to τ0,minT0,max ≈ 0.74, inconsistent with Eq. (36) and with Fig. 3. Consequently, the numbers in the abstract and Sec. IV for the conformal case change to τ0,min ≈ 1.6 fm and e0,max ≈ 3.0 GeV/fm^3 for RHIC (E0 = 5 GeV/fm^2) and τ0,min ≈ 1.2 fm and e0,max ≈ 12 GeV/fm^3 for LHC (E0 = 14 GeV/fm^2). The lattice-EoS values in Sec. III D are obtained by an analogous read-out and would also shift. The qualitative conclusion that nonlinear causality constrains initial conditions is unaffected, but the quantitative claims are not justified by the stated equations.","section":"§III C, Eq. (37)"},{"comment":"The conversion from the 1D Bjorken constraints to RHIC/LHC initial conditions assumes that the early-stage expansion is effectively one-dimensional and boost invariant at the relevant times (τ0 ≈ 0.5–1.6 fm). Transverse expansion and finite-size gradients are neglected. Since the derived τ0,min and e0,max are applied to heavy-ion collisions through the Bjorken formula, this is a load-bearing premise: in a more 3D expansion, the causality constraints could be weaker or stronger. The manuscript should discuss the expected magnitude of this effect or explicitly limit the heavy-ion conclusions to the 1D approximation.","section":"§III C; §IV"},{"comment":"The quantitative bounds for the lattice EoS (Sec. III D, Figs. 6–7) use the N=4 SYM transport coefficients (5) together with a nonconformal lattice EoS. This is acknowledged in Sec. II A, but the resulting systematic uncertainty is not estimated; the lattice-EoS numbers are presented in the abstract as if they were more realistic than the conformal ones. The authors should either provide an estimate of how the bounds vary with Cη, Cτπ, Cλ or clearly flag that the numbers are only illustrative examples.","section":"§II A; §III D"}],"minor_comments":[{"comment":"The notation 'E0 (GeV/fm^2) = E0/(ℏc) (1/fm^3)' is dimensionally confusing; please rewrite the unit conversion explicitly, e.g., state that in natural units with ℏc = 0.197 GeV fm, 1 GeV/fm^2 ≈ 5.07 fm^{-3}.","section":"Eq. (34)"},{"comment":"The extrapolation of the lattice EoS to temperatures outside the fitted range (0.13–0.4 GeV) is stated, but its impact on the constraints is not quantified; a sentence noting that the results for very high T are extrapolations would be helpful.","section":"Sec. II B"},{"comment":"The captions contain typos such as 'fmτ0=0.20' where a space is missing; please fix these typographical issues.","section":"Figs. 4, 5"},{"comment":"The conclusion that the initial stage needs a non-equilibrium description is stated in absolute terms; adding the qualifier 'within the 1D Bjorken approximation' would be more accurate.","section":"Abstract and Sec. IV"},{"comment":"The bounds in Eqs. (31) and (32) are specific to the conformal EoS and N=4 SYM transport coefficients; this should be stated in the main text immediately after the equations, not only in the caption of Fig. 2.","section":"Sec. III B"},{"comment":"The notation f = 1 ± 1/2 Re^{-1} is slightly confusing because the sign depends on the sign of φ; it may help to write f = 1 + 3φ/(8e) explicitly.","section":"Appendix A, Eq. (A2)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the central idea is interesting. The arithmetic error in Eq. (37) is easily corrected, but until it is fixed the quantitative claims should not be used. The experimental application is an extrapolation from 1D; the authors should be encouraged to either include a caveat or perform a more realistic 3D check. Overall, I see no issue of circularity or fabrication; the work is a legitimate application of known causality conditions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things you should know. The paper does something genuinely useful: it takes the nonlinear causality conditions of Bemfica et al. and turns them into inverse-Reynolds-number windows and initial-condition exclusion maps for one-dimensional Bjorken flow. The transcription of the inequalities looks careful, the ODE evolution is straightforward, and the conformal-limit checks against Ref. [20] are consistent. The second thing is less pleasant: the quantitative headline does not follow from the equations. Combining their Eq. (36), tau0,min/tau_pi = tau0,min T0,max / C_tau_pi ~ 5.3, with tau_pi = C_tau_pi/T and the Nf=3 ideal gas EoS gives tau0,min T0,max = 5.3 C_tau_pi ~ 1.10, and hence tau0,min ~ 2.8 E0^{-1/3}, not 1.6 E0^{-1/3}. I checked this twice; the stress-test note is right. The corrected numbers for RHIC are tau0,min ~ 1.6-1.7 fm and e0,max ~ 3 GeV/fm3, not 1 fm and 5; for LHC, ~1.2 fm and ~12 GeV/fm3, not 0.7 and 20. That changes the abstract substantially.\n\nOther soft spots are real but less severe. The one-dimensional boost-invariant geometry is assumed to describe the early expansion, and transverse expansion could shift constraints; the paper does not discuss this much. The N=4 SYM transport coefficients are used with the lattice EoS, which is inconsistent in principle, and the lattice EoS is extrapolated outside its fitted temperature range. There are no propagated uncertainties. I do not see circularity: the inequalities come from an external group, and the dET/dy input is measured.\n\nCredit where it is due: the paper is honest about the difference between necessary and sufficient conditions, gives the maps so a reader can see what drives the constraint, and cites the prior one-dimensional causality work. The qualitative conclusion — large inverse Reynolds number violates these causality conditions, so reasonable initial conditions are seriously constrained — survives the arithmetic slip. The slip is in the conversion to the phenomenological numbers, not in the causality inequalities themselves.\n\nWho this is for: people setting Bayesian priors or building pre-hydro models. They should not use the current numbers. The paper deserves a serious referee; I would send it out with a clear request to fix Eq. (37), rerun the figures and abstract, and preferably propagate uncertainties. After that it would be a useful contribution.","headline":"The causality-to-initial-condition pipeline is coherent and worth taking seriously, but the headline tau0,min and e0,max numbers are off by about a factor 1.7 because Eq. (37) does not follow from the paper's own equations.","tokens_in":19509,"tokens_out":5462,"would_cite":false,"duration_ms":58184,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["25.75.-q","47.75.+f"],"model":"deepseek-v4-flash","headline":"One-dimensional expanding viscous fluids violate causality when the inverse Reynolds number is large, so hydrodynamic initial conditions are sharply constrained.","keywords":["nonlinear causality","inverse Reynolds number","BRSSS equation","Bjorken expansion","relativistic heavy-ion collisions","initial conditions","quark-gluon plasma","second-order hydrodynamics"],"falsifier":"Run a (2+1)- or (3+1)-dimensional BRSSS simulation with the same transport coefficients, starting from an initial condition the one-dimensional criterion labels acausal (for example $\\tau_0 = 0.2$ fm with $\\varphi_0 = 0$ at RHIC energy); if no characteristic velocity ever exceeds the speed of light during the first few fm/c, then the one-dimensional reduction, rather than the nonlinear causality conditions, produces the claimed bound.","tokens_in":2239,"feed_emoji":"⚳️","tokens_out":3245,"duration_ms":110592,"temperature":0.7,"pith_summary":"The paper asks how soon after two heavy nuclei collide the hot, dense matter can be treated as a relativistic fluid, and answers that nonlinear causality sets a sharp limit. For a one-dimensional, boost-invariant (Bjorken) expansion of a viscous fluid obeying the BRSSS constitutive equation, the authors show that the necessary and sufficient conditions for causal propagation fail whenever the inverse Reynolds number—the size of the shear stress relative to enthalpy, $|\\varphi|/(e+p)$—becomes large. Requiring the solution to remain causal for its entire evolution therefore forbids many otherwise plausible initial conditions. Combined with measured transverse energy per rapidity at RHIC and LHC, this yields minimum initial proper times around 0.4–1 fm and maximum initial energy densities around 5–35 GeV/fm$^3$, depending on the equation of state. The conclusion is that the earliest stage of heavy-ion collisions needs a description other than this form of dissipative hydrodynamics.","feed_headline":"Expanding viscous fluids violate causality far from equilibrium","feed_subtitle":"Heavy-ion flow cannot start before about 0.4–1 fm, depending on the equation of state.","key_machinery":"The central object is the BRSSS constitutive equation for the shear stress tensor, which promotes dissipative currents to dynamical variables and is the basis of the causal second-order hydrodynamic framework. The argument proceeds by translating the general necessary and sufficient conditions of nonlinear causality, obtained from the characteristic equation $\\det(A^\\alpha \\xi_\\alpha) = 0$ of the quasi-linear system, into inequalities involving the energy density $e$, pressure $p$, shear pressure $\\varphi$, sound velocity $c_s^2$, shear viscosity $\\eta$, and relaxation time $\\tau_\\pi$ in the specific one-dimensional boost-invariant flow. The inverse Reynolds number $\\mathrm{Re}^{-1} = |\\varphi|/(e+p)$ is the organizing variable: the causality conditions collapse into a narrow allowed band for this dimensionless measure of departure from local equilibrium, and the time evolution of $\\varphi/(e+p)$ through the hydrodynamic equations determines whether a given initial condition stays causal.","core_discovery":"The paper's central claim is that for boost-invariant one-dimensional expansion, the BRSSS constitutive equation has a finite causal window in the inverse Reynolds number. With the conformal equation of state and the AdS/CFT transport coefficients, the necessary conditions of nonlinear causality require $-0.47 \\le \\varphi/(e+p) \\le 0.23$, while the sufficient conditions require $-0.07 \\le \\varphi/(e+p) \\le 0.07$; any solution that leaves this window at any time is acausal. Because a Bjorken-expanding system can be driven far from local equilibrium even from $\\varphi_0 = 0$ at very early proper times, the expansion itself forces a violation unless the initial proper time is large enough. Combining the causal-region boundary with the Bjorken energy-density formula $e_0 = (1/S)\\,dE_T/dy$ measured at RHIC and LHC gives, for the conformal case, $\\tau_{0,\\min} \\approx 1$ fm and $e_{0,\\max} \\approx 5$ GeV/fm$^3$ at RHIC and $\\tau_{0,\\min} \\approx 0.7$ fm and $e_{0,\\max} \\approx 20$ GeV/fm$^3$ at LHC; with the lattice equation of state the bounds become $\\tau_{0,\\min} \\approx 0.5$ fm and $e_{0,\\max} \\approx 11$ GeV/fm$^3$ at RHIC and $\\tau_{0,\\min} \\approx 0.4$ fm and $e_{0,\\max} \\approx 35$ GeV/fm$^3$ at LHC.","pith_inferences":["If the early expansion is genuinely three-dimensional rather than one-dimensional Bjorken flow, the expansion scalar is generally smaller at a given longitudinal proper time, so the derived bounds could weaken or strengthen depending on the velocity gradients; the 1D assumption is the main caveat in applying these numbers to real collisions.","The gap between the necessary and sufficient conditions leaves a band of initial conditions whose causal status is undetermined; directly computing the characteristic velocities of simulated 3D flows could shrink that band and sharpen the lower bound on $\\tau_0$.","Because the bounds depend on the chosen constitutive equation and transport coefficients, repeating the analysis with other second-order hydrodynamic schemes would quantify how much of $\\tau_{0,\\min}$ and $e_{0,\\max}$ is scheme-specific.","The inverse-Reynolds-number criterion could be applied locally to the output of pre-equilibrium models to decide, event by event and time by time, where a causal hydrodynamic description becomes legitimate rather than assuming a single global initial time."],"forward_implications":["Initial conditions for hydrodynamic simulations of heavy-ion collisions must be chosen so that the entire trajectory, not just the initial state, stays inside the causal region.","Starting from local equilibrium at sufficiently early times is ruled out, because the large expansion rate $\\partial_\\mu u^\\mu = 1/\\tau$ drives the system out of the causal window even when $\\varphi_0 = 0$.","For a conformal equation of state with $N_f = 3$, the minimum initial proper time satisfies $\\tau_{0,\\min} \\approx 1.6 E_0^{-1/3}$ fm with $E_0 = (1/S)\\,dE_T/dy$ in GeV/fm$^2$, giving the RHIC and LHC bounds quoted above.","The lattice equation of state broadens the acceptable region, shifting the bounds to smaller $\\tau_{0,\\min}$ and larger $e_{0,\\max}$ than the conformal case.","The pre-hydrodynamic stage of relativistic heavy-ion collisions must be described by a non-equilibrium framework other than the BRSSS/Müller–Israel–Stewart dissipative hydrodynamics."],"supporting_citations":[{"why":"Supplies the general necessary and sufficient conditions of nonlinear causality that the paper translates into the one-dimensional expanding-fluid inequalities.","marker":"[20]"},{"why":"Defines the BRSSS constitutive equation whose solutions are tested for causality.","marker":"[26]"},{"why":"Provides the boost-invariant Bjorken flow solution and the formula $e_0 = (1/S)\\,dE_T/dy$ connecting data to initial conditions.","marker":"[29]"},{"why":"Provides the lattice QCD equation-of-state parametrization used in the realistic case.","marker":"[28]"},{"why":"Supplies the conformal transport coefficients used in all numerical evaluations.","marker":"[27]"},{"why":"Measured transverse energy per rapidity in central Au+Au collisions at 200 GeV, fixing $e_0\\tau_0 = 5$ GeV/fm$^2$.","marker":"[37]"},{"why":"Measured transverse energy per rapidity in central Pb+Pb collisions at 2.76 TeV, fixing $e_0\\tau_0 = 14$ GeV/fm$^2$.","marker":"[38]"}],"fun_headline_variants":["Causality forces a minimum start time for heavy-ion hydrodynamics","Viscous fluid expansion goes acausal before 0.4 fm","Heavy-ion flow: earliest start set by nonlinear causality","No hydrodynamic flow before 0.4 fm in heavy-ion collisions","Nonlinear causality delays the onset of hydrodynamics in heavy-ion collisions"],"cache_read_input_tokens":21504,"weakest_assumption_plain":"The one-dimensional, boost-invariant Bjorken flow is assumed to capture the expansion rate of the earliest fluid stage in heavy-ion collisions, and the derived bounds are applied to RHIC and LHC through the Bjorken formula; if transverse expansion and finite-size gradients are important at those times, the causality constraints and the resulting $\\tau_{0,\\min}$ and $e_{0,\\max}$ could be weaker or stronger.","fun_headline_variants_meta":{"raw":{"variants":["Causality forces a minimum start time for heavy-ion hydrodynamics","Viscous fluid expansion goes acausal before 0.4 fm","Heavy-ion flow: earliest start set by nonlinear causality","No hydrodynamic flow before 0.4 fm in heavy-ion collisions","Nonlinear causality delays the onset of hydrodynamics in heavy-ion collisions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001379,"raw_usage":{"total_tokens":5657,"prompt_tokens":1086,"completion_tokens":4571,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":702,"completion_tokens_details":{"reasoning_tokens":4481}},"tokens_in":702,"tokens_out":4571,"duration_ms":32775,"temperature":1.0,"reasoning_tokens":4481,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T23:31:24.140574+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a (2+1)- or (3+1)-dimensional BRSSS simulation with the same transport coefficients, starting from an initial condition the one-dimensional criterion labels acausal (for example $\\tau_0 = 0.2$ fm with $\\varphi_0 = 0$ at RHIC energy); if no characteristic velocity ever exceeds the speed of light during the first few fm/c, then the one-dimensional reduction, rather than the nonlinear causality conditions, produces the claimed bound.","supporting_citations":[{"cited_title":"Measurements of transverse energy distributions in Au+Au collisions at $\\sqrt{s_{NN}}= 200$ GeV","cited_arxiv_id":"nucl-ex/0407003","evidence_quote":"Measured transverse energy per rapidity in central Au+Au collisions at 200 GeV, fixing $e_0\\tau_0 = 5$ GeV/fm$^2$."}],"review_version":1}