{"id":"a3035eac-9dba-445e-81ae-18e87821180d","arxiv_id":"2509.08864","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Second- and higher-order hydrodynamic attractors converge to the same universal curve before the Navier-Stokes limit, with the leading correction controlled only by second-order coefficients.","lead":"This paper compares hydrodynamic theories of different orders and finds that their attractor curves merge before the Navier-Stokes limit. This supports using simpler second-order equations, with correct coefficients, to model the expanding quark-gluon plasma.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Cross-theory attractor universality is only tested at c2,2=0; the paper's own Fig. 5 shows physical c2,2 values shift the attractor, so the w~20 merging claim is not established for actual DNMR/MIS coefficients.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing concern: the cross-theory comparison sets c2,2=0 via the relaxation-time approximation, while Figure 5 shows that nonvanishing c2,2 changes convergence. This is the most important issue because the paper's headline claim--that all theories with identical second-order coefficients converge to a universal attractor well before the Navier-Stokes limit--is numerically tested only on the c2,2=0 slice. The slow-roll analysis shows that second-order coefficients control the leading correction, so the claim is internally consistent; the gap is in the numerical generality. The paper is honest about the choice and includes Figure 5, which is creditworthy, but the central quantitative statement (w~20 merging) is not demonstrated for physical nonzero c2,2 values. A concrete rerun of the random-coefficient scan at fixed nonzero c2,2 would settle the issue. The reader's CONDITIONAL verdict is appropriate; no adjustment is needed, though the condition should explicitly require the nonzero-c2,2 test.","tokens_in":10625,"tokens_out":5956,"duration_ms":63929,"concrete_test":"Repeat the random-coefficient scan of Fig. 2 (10^4 samples) with c2,2 fixed to 27/70 (DNMR hard-sphere) and separately to -3 C_lambda/(8 C_eta) with C_lambda=1/pi, keeping all higher-order coefficient ranges as in the paper. Compute the band of attractor curves phi(w) for standalone Eq. (5) and coupled Eq. (6) for each fixed c2,2. If the relative spread phi/phi_ref - 1 remains within, say, 1% for w>=20 for both nonzero c2,2 values, the 'attractor of attractors' is robust to the RTA choice; if the spread is comparable to the gap between the c2,2=0 and c2,2!=0 curves in Fig. 5, the claimed universality is an artifact of the relaxation-time approximation and the conclusion should be restricted to matched c2,2=0 models.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The numerical demonstration of an 'attractor of attractors' is built on setting c2,2=0 for DNMR, PJP, BD, and YJG via the relaxation-time approximation (paragraph before Eq. 5), and the random-coefficient bands in Fig. 2 also fix c2,2=0. The slow-roll coefficient A2 (Eqs. 8 and 10) depends explicitly on c2,2 through (2/3 - c2,1)C_tau + (16/9)c2,2 C_eta, so changing c2,2 at fixed c2,1 changes the leading sub-Navier-Stokes correction. The paper's own Fig. 5 shows that with nonvanishing c2,2 (DNMR hard-sphere 27/70; MIS -3 C_lambda/(8 C_eta)), the attractors do not coincide with the c2,2=0 curves at w~20; convergence to the common curve is postponed. Therefore the central claim that higher-order theory converges to the same attractor as a second-order theory 'well before the Navier-Stokes limit' is only demonstrated for one specially chosen second-order coefficient. The more general statement--'all theories with identical second-order coefficients converge to a universal attractor'--is plausible but untested for nonzero common c2,2 values. Since the abstract's practical conclusion tells users to adopt 'correct transport coefficients,' the coefficient choice is not cosmetic: if physical c2,2 values differ across theories, there is no single universal attractor, only coefficient-dependent attractors. The paper does disclose the RTA choice and Figure 5, but disclosure does not remove the dependence of the central quantitative claim on this choice.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies attractor solutions of transient relativistic hydrodynamics in a conformal, boost-invariant, transversely homogeneous Bjorken flow. For several second-order (MIS, DNMR) and third-order (PJP, BD, YJG) theories, written as ordinary differential equations for the shear anisotropy φ and, for BD/YJG, a fourth-rank moment Φ, it compares the exact attractor with a late-time slow-roll expansion in powers of 1/w. It derives the expansion coefficients A_s^3 and A_c^3, shows that the A2/w^2 subleading term depends only on second-order transport coefficients, and demonstrates numerically that the listed theories with common second-order coefficients merge at w~20 (MIS at w~100) before the Navier-Stokes limit at w~400. The conclusion is that the energy-momentum tensor evolution of higher-order theories is reproduced by properly matched second-order hydrodynamics from an intermediate time onward.","tokens_in":11078,"tokens_out":15451,"duration_ms":171573,"significance":"If correct, the result is practically useful: it would make second-order hydrodynamic simulations adequate at intermediate and late times even when far-from-equilibrium corrections require higher-order moments. The slow-roll coefficients are obtained by direct substitution into the stated equations and are not fitted to the target attractor, and the numerical attractors are computed separately from the evolution equations, so the main line is not circular. The random-coefficient scan over 10^4 samples is a genuine stress test of sensitivity to higher-order coefficients. However, the claim as stated needs additional support for nonzero c2,2 values before the practical conclusion is fully established.","major_comments":[{"comment":"The quantitative cross-theory merging claim at w≳20 is only demonstrated for the relaxation-time value c2,2=0. A2 in Eq. (8) is linear in (16/9)c2,2Cη, so at fixed c2,1 the subleading correction changes with c2,2. The paper's own Supplemental Fig. 5 shows that the DNMR hard-sphere value c2,2=27/70 and the MIS value c2,2=-3Cλ/(8Cη) delay convergence to the c2,2=0 attractor. Since the abstract's practical conclusion invokes \"correct transport coefficients,\" this is not a cosmetic choice. The central claim has not been tested for common nonzero c2,2 values. I request either a direct numerical scan over c2,2 (with fixed common c2,1 and random higher-order coefficients) with reported merging times, or a restriction of the advertised conclusion to the c2,2=0/RTA class.","section":"Eqs. (8)-(10), Supplemental Fig. 5"},{"comment":"The statement that theories with the same combination c2,1 - (16Cη/(9Cτ))c2,2 exhibit identical late-time behavior is stronger than Eq. (8) supports. A_s^3 depends on c2,1^2, c2,2^2 and c2,1c2,2 separately, not only through this combination. Thus two theories can have the same A2 but different A3 and different convergence before φNS; this is also consistent with MIS's delayed convergence in Fig. 2. Please qualify the claim to \"identical at the A2/w^2 order\" or prove that the extra terms cancel at the relevant w.","section":"Section \"Attractors and asymptotic behaviors\", paragraph after Eq. (10)"}],"minor_comments":[{"comment":"Please define Cη, Cτ, CΦ, and Cλ explicitly at first use; Cτ appears in Eq. (5), CΦ in Eq. (6), and Cλ later in the text without formal definitions.","section":"Definitions around Eqs. (5)-(6)"},{"comment":"DMMZ is listed as a third-order theory in the introduction but has no coefficients in Table I and is absent from Fig. 2; state explicitly that it is not included in the quantitative comparison because its coefficients are unspecified.","section":"Introduction and Table I"},{"comment":"The numerical procedure for extracting the attractor (e.g., backward integration from large w or solving the boundary condition at w→0) is not described. A short paragraph in the Supplemental Material would make the figures reproducible.","section":"Numerical methods"},{"comment":"The \"converge by w≳20\" statement is assessed via visual percent-deviation insets. Give a quantitative convergence threshold and list the deviations for each curve, including the c2,2≠0 cases in Fig. 5.","section":"Figs. 2 and 5"}],"recommendation":"major_revision","confidential_remarks":"The slow-roll analysis appears algebraically sound and the paper is not circular; the main gap is the mismatch between the demonstrated c2,2=0 case and the abstract's broader practical message. This is fixable with additional numerical tests and wording changes. I would not reject for novelty or scope; the \"attractor of attractors\" idea is publishable once the nonzero-c2,2 case is addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the cross-theory comparison: prior work studied attractors inside individual theories, and showing that the leading slow-roll correction A2 depends only on second-order coefficients is a clean, useful observation. The algebra is consistent—I spot-checked A1 and A2 in the standalone equation—and the numerical attractor solutions are computed independently, so the core of the paper is trustworthy.\n\nWhat the paper does well: it states its assumptions openly, including the relaxation-time approximation that sets c2,2=0 for DNMR and MIS, and it includes Fig. 5 showing that nonzero c2,2 shifts the attractor. That disclosure is honest, and it frames the main caveat. The conclusion that second-order hydrodynamics with correctly chosen coefficients suffices from intermediate times onward is actually robust: the slow-roll expansion shows higher-order coefficients enter only at order 1/w^3, so the leading sub-Navier-Stokes behavior is governed by c2,1 and c2,2 regardless of what c3's or moment couplings are.\n\nThe soft spots are real but not fatal. The headline claim that all theories converge to a single universal attractor by w~20 is only demonstrated for c2,2=0. Fig. 5 shows that with the physical hard-sphere DNMR value or the MIS value, the curves separate and convergence is postponed. So the strong quantitative statement is conditional on a coefficient choice; the more general point that theories with identical second-order coefficients share a common attractor remains plausible, and the slow-roll analysis actually supports it for any common c2,2. The paper just doesn't test it there. I'd also like to see code or at least tabulated attractor curves—the numerical solutions are central to the w~20 claim and nothing is shipped. That's a fixable omission.\n\nMinor points: DMMZ is listed in Table I but no coefficients are given, so it drops out of the comparison; fine but worth a sentence. The random-coefficient bands also fix c2,2=0, so they corroborate insensitivity to third-order terms only under that same choice.\n\nWho is this for: anyone using viscous hydrodynamics in heavy-ion phenomenology, and people working on attractor theory. It deserves a serious referee. I would not desk-reject; I'd send it out and ask for a sensitivity scan over c2,2 (with physical values) and either code or data. The central insight is likely correct, and the caveats are addressable.\n\nRecommendation: peer review, with moderate revision expected.","headline":"Useful cross-theory attractor comparison; the structural result is sound, but the quantitative merge-by-w~20 claim depends on a zero-c2,2 choice that the paper's own Fig. 5 shows matters.","tokens_in":11468,"tokens_out":1657,"would_cite":true,"duration_ms":20817,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Higher-order hydrodynamic attractors converge to the second-order attractor before the Navier-Stokes limit.","keywords":["hydrodynamic attractor","heavy-ion collisions","quark-gluon plasma","second-order hydrodynamics","slow-roll expansion","transport coefficients","Bjorken flow","hydrodynamization"],"falsifier":"Solve Eq. (5) for DNMR and PJP with c2,2=27/70 instead of 0 and record where their attractors merge with the coupled-theory attractor; the paper's own Figure 5 indicates a later merging, so the claim that all theories with identical second-order coefficients merge by w~20 would be falsified if the merging point shifts beyond the intermediate window.","tokens_in":1557,"feed_emoji":"","tokens_out":6731,"duration_ms":130530,"temperature":0.7,"pith_summary":"The paper asks whether higher-order hydrodynamic theories, which capture far-from-equilibrium physics, eventually agree with the commonly used second-order equations. Studying boost-invariant conformal plasmas, it claims that all theories sharing the same second-order transport coefficients converge to the same attractor by scaled time w~20, well before any of them reaches the Navier-Stokes solution at w~400. This means the intermediate-time evolution of the energy-momentum tensor is controlled by second-order coefficients alone, regardless of higher-order structure. If true, second-order hydrodynamics with correctly chosen coefficients is adequate from an intermediate time onward, even when the far-from-equilibrium stage would require higher-order theory.","feed_headline":"Hydrodynamic attractors merge before viscosity sets in","feed_subtitle":"Theories with matching second-order coefficients converge by w~20, so lower-order codes stay valid from intermediate times.","key_machinery":"The key objects are the normalized shear anisotropy phi = Pi/epsilon and the scaled time w = tau T, which turn the conservation and relaxation equations into first-order ODEs: a standalone equation for phi (Eq. 5) or a coupled system with a four-index non-equilibrium moment Phi (Eq. 6). The load-bearing identity is the slow-roll expansion phi = sum A_n/w^n: A1 is the Navier-Stokes term, A2 depends only on the second-order coefficient combination (c2,1 - 16 C_eta/(9 C_tau) c2,2), and third-order coefficients appear only in A3 and beyond. Because the compared theories share c2,1=10/21 and c2,2=0, their attractors coincide before approaching the Navier-Stokes limit.","core_discovery":"The paper demonstrates that, for boost-invariant, transversely homogeneous conformal systems, the attractor of a higher-order hydrodynamic theory converges to the attractor of a second-order theory before either reaches the Navier-Stokes limit. Writing the shear anisotropy in a slow-roll expansion in 1/w, the leading non-trivial coefficient depends only on the second-order coefficients c2,1 and c2,2; higher-order transport coefficients enter only at the next order. Since the compared second- and third-order theories share the same second-order coefficients, their attractors merge by w~20. The paper calls this an 'attractor of attractors' and concludes that second-order hydrodynamics with pro","pith_inferences":["The quantitative merging point w~20 likely depends on the choice c2,2=0; with the hard-sphere value c2,2=27/70, Figure 5 suggests the convergence window shifts, so the universal claim is really a statement about the relaxation-time approximation.","The slow-roll argument probably extends to non-conformal and transversely inhomogeneous systems, which the paper lists as a next step; if so, second-order codes would remain sufficient for a wider class of observables.","The attractor-of-attractors hierarchy implies that the leading late-time stress-tensor behavior is fixed by the second-order transport sector alone, which could guide how kinetic-theory results are matched to hydrodynamic codes."],"forward_implications":["All theories with identical second-order coefficients share the same attractor by w~20, well before any reaches the Navier-Stokes limit at w~400.","Higher-order transport coefficients and additional four-index moments affect only the A3 and higher slow-roll terms, so they do not change intermediate-time stress-tensor evolution.","Second-order hydrodynamics with correctly chosen transport coefficients is adequate from an intermediate time onward, even when the far-from-equilibrium regime would require a higher-order description.","Slow-roll truncations at third order reproduce the full attractor to about 10^-3 relative accuracy for w>~5, while Navier-Stokes converges only for w>~100.","Practical heavy-ion simulations can retain second-order equations without losing accuracy for late-time collective flow."],"supporting_citations":[{"why":"Introduced the hydrodynamic attractor and the slow-roll expansion method on which the paper's asymptotic analysis is built.","marker":"[20]"},{"why":"Supplies the relaxation-time-approximation coefficients, including c2,2=0, adopted for cross-theory comparison.","marker":"[52]"},{"why":"Provides the DNMR second-order derivation and the second-order coefficient values used in the comparison.","marker":"[14]"},{"why":"Gives the Muller-Israel-Stewart second-order theory whose attractor is compared in Figure 2.","marker":"[12, 13]"},{"why":"Provides the third-order theory (PJP) and its transport coefficients used in the standalone comparison.","marker":"[54]"},{"why":"Supplies the coupled four-index-moment third-order theory whose attractor and coefficients enter the comparison.","marker":"[56]"},{"why":"Supplies the coupled two-moment theory used as the main numerical example and reference target for ratio comparisons.","marker":"[57]"},{"why":"Provides the analytical slow-roll expansion framework used to determine which transport coefficients enter at each order.","marker":"[51]"}],"fun_headline_variants":["Attractors of hydrodynamics merge by w~20","Second-order hydro survives far-from-equality","Higher-order attractors collapse to second-order","Why second-order hydro wins from mid-times","Attractor convergence saves lower-order codes"],"cache_read_input_tokens":13184,"weakest_assumption_plain":"The cross-theory comparison sets the nonlinear coefficient c2,2 to zero for the DNMR and PJP theories by adopting the relaxation-time approximation (stated just before Eq. 5), yet Figure 5 shows that choosing the hard-sphere value c2,2=27/70 changes how quickly each theory's attractor converges.","fun_headline_variants_meta":{"raw":{"variants":["Attractors of hydrodynamics merge by w~20","Second-order hydro survives far-from-equality","Higher-order attractors collapse to second-order","Why second-order hydro wins from mid-times","Attractor convergence saves lower-order codes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":7.9e-05,"raw_usage":{"total_tokens":729,"prompt_tokens":646,"completion_tokens":83,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":390,"completion_tokens_details":{"reasoning_tokens":25}},"tokens_in":390,"tokens_out":83,"duration_ms":2145,"temperature":1.0,"reasoning_tokens":25,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T20:52:51.596137+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve Eq. (5) for DNMR and PJP with c2,2=27/70 instead of 0 and record where their attractors merge with the coupled-theory attractor; the paper's own Figure 5 indicates a later merging, so the claim that all theories with identical second-order coefficients merge by w~20 would be falsified if the merging point shifts beyond the intermediate window.","supporting_citations":[],"review_version":1}