{"id":"68d2267b-3e55-4f32-b3df-f87b8346b619","arxiv_id":"2608.01900","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"The leading-order VAH equations stay causal exactly when the three dimensionless ratios in Eq. (21) satisfy the stated chain inequalities.","lead":"This paper derives a set of simple inequalities that guarantee the equations of viscous anisotropic hydrodynamics (VAH) do not let information travel faster than light. The result gives heavy-ion physicists a precise way to tell when the early-stage hydrodynamic model is physically valid.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified; Eq. (21) is correctly derived for the truncated leading-order VAH system and the NLO caveat is explicitly scoped.","rationale":"The paper's central claim is Eq. (21). My independent review of the algebra: the characteristic polynomial for the local-rest-frame 6x6 system in Eq. (A6) has a double b=0 root arising from the constraints; after removing it, Eq. (14) is a quadratic in v^2. Conditions (a)-(d) in Sec. III are the standard necessary and sufficient conditions for a monic quadratic to have both roots in [0,1], and Appendices B and C reduce them to Eq. (21). I found no algebraic error. The derivation depends only on the state-dependent coefficients Gamma, and the inequality chain follows because condition (a) forces the cross product to be no larger than the diagonal product, while condition (d) forces it to be nonnegative. The only limitation is that pi_perp and W_perp_z, generated by the NLO correction delta f, are omitted; the authors explicitly state this in Sec. V and list the complete-VAH extension as future work. Therefore, the theorem is correct for the truncated system, and the caveat is explicit. This does not change the ACCEPT verdict.","tokens_in":18279,"tokens_out":21575,"duration_ms":205840,"concrete_test":"Implement a symbolic check of det(A^alpha xi_alpha) for the 6x6 matrix in Eq. (A6) to confirm it factorizes as a nonzero constant times b^2 times the left-hand side of Eq. (14), leaving no additional nondegenerate characteristics beyond those in Eq. (16).","verdict_should_be":"UNCHANGED","load_bearing_attack":"I found no internal flaw in the derivation of Eq. (21). For the local-rest-frame characteristic problem in Eq. (A6), the constraint-related b=0 modes factor out, leaving the quartic (14); requiring both roots of Eq. (16) to lie in [0,1] for every direction is exactly the condition that the monic quadratic has nonnegative roots not exceeding 1, and the reductions in Appendices B and C are algebraically consistent. The only substantive limitation is the deliberate omission of the NLO correction delta f and the associated pi_perp and W_perp_z, which means Eq. (21) certifies causality of the truncated system (5)-(9), not of the complete VAH theory. This limitation is stated explicitly in Sec. V, so the paper does not overclaim its scoped result.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives necessary and sufficient conditions for nonlinear causality in viscous anisotropic hydrodynamics (VAH), restricted to the leading-order truncation that neglects the next-to-leading-order correction δf. The equations of motion (5)-(9) are treated as a first-order quasilinear system, and characteristic analysis in the local rest frame reduces the characteristic determinant to the quartic (14), then to a quadratic in v_c^2. Requiring that for every propagation direction both roots lie in [0,1] yields the chain of inequalities (21) on the combinations ζ̄L_z/(E+P_L), ζ̄⊥⊥/(E+P_⊥), and the product (ζ̄L_⊥/(E+P_⊥))(ζ̄⊥_z/(E+P_L)). Section IV gives a physical interpretation in terms of two coupled modes, and Section V notes that incorporating NLO corrections is future work.","tokens_in":18419,"tokens_out":20976,"duration_ms":179291,"significance":"If correct, the result is practically valuable: it converts the causality requirement for the VAH equations into three simple inequalities that can be checked locally in simulations, analogous to the nonlinear causality constraints already used for conventional viscous hydrodynamics. The work goes beyond the earlier conformal first-order analysis of Ref. [74] by treating the second-order nonconformal anisotropic hydrodynamics. The derivation is verifiable, and I independently confirmed the reduction from Eq. (A6) to Eq. (14). The paper is honest about its scope, explicitly limiting the claim to the truncated system (5)-(9), and the final conditions are falsifiable: any simulation violating (21) will exhibit superluminal propagation. The main limitation, that the neglected δf and its associated π⊥, W⊥ modes could in principle introduce additional superluminal characteristic modes, is stated clearly in Sec. V, so the paper does not overclaim its scoped result.","major_comments":[],"minor_comments":[{"comment":"The reduction from the 6x6 determinant (A6) to the characteristic equation (14) is stated rather than shown. I verified the expansion independently, but please include the intermediate calculation or provide it as supplementary material so that the central equation is fully reproducible.","section":"Appendix A, Eq. (14)"},{"comment":"The derivation of condition (d) divides by (ΓL_z)^2 when writing the symmetry axis of k(t), so the case ΓL_z=0, which is allowed by the final inequality (21a), is not covered by the argument as written; the same concern applies to Γ⊥⊥=0. The final inequalities remain correct (they can be obtained by a separate direct check or by a continuity argument), but the proof should treat these boundary cases explicitly.","section":"Appendix B, step 5, Eq. (B24)"},{"comment":"The denominator in the displayed formula for v_c^2 is typeset as a_3^3; it should read a_3^2.","section":"Eq. (13)"},{"comment":"The sentence beginning 'we derived the necessary and sufficient conditions' is grammatically incomplete, and the paragraph should explicitly restate that the conditions apply to the truncated system (5)-(9), in line with the qualification in the abstract and in the final paragraph of Section V.","section":"Sec. V, first paragraph"},{"comment":"The ratios in Eq. (21) are meaningful only when E+P_L and E+P_⊥ are positive; the paper should state this standard assumption explicitly.","section":"Sec. II, after Eq. (4)"}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is a solid, verifiable contribution. I independently checked the determinant expansion and found no substantive technical error. The main limitations, namely the neglect of NLO corrections and the degenerate Γ=0 case in the proof, are explicit or easily fixed. I recommend minor revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe genuinely new result is Eq. (21), a compact chain of inequalities involving the four anisotropic transport coefficients that is necessary and sufficient for all characteristic speeds of the leading-order VAH equations to lie between zero and the speed of light. If you work with VAH simulations, this is immediately useful: it turns a causality check into a simple evaluation of four numbers. I checked the 6x6 determinant in Appendix A and reproduced Eq. (14); the two derivations in Appendices B and C are consistent, and the product-mixing inequality in (21c) is not in the earlier literature.\n\nThe paper earns its keep through transparency. The assumptions are stated up front: the NLO correction δf and the π⊥, W⊥z currents it generates are dropped, so the conditions certify causality of the reduced system (5)-(9), not of the complete VAH theory. This limitation is explicit in Sec. V, which I appreciate. The physical discussion in Sec. IV is heuristic but helpful, especially the two-fluid interpretation when the mixing terms vanish.\n\nSoft spots, in proportion. First and main: the scope is the truncated system. The authors acknowledge this, but a reader applying (21) to full VAH should be careful; the omitted NLO modes could in principle be superluminal. That is a real gap between the advertised 'regime of validity of VAH' in the abstract and what is proven. Second, the stability sentence at the end of Sec. V (\"both thermodynamic stability and causality are maintained\") goes beyond the derivation: real non-negative characteristic speeds give hyperbolicity and causality, not thermodynamic stability. That should be softened. Third, Appendix A states the determinant result rather than showing the expansion; a referee will want the intermediate algebra, though the final expression checks out.\n\nThe method itself is the established characteristic-velocity program from Bemfica et al.; the novelty is the application to this second-order anisotropic system and the clean inequality chain that emerges. That is a modest but solid step, not a paradigm shift. The citation pattern looks appropriate and includes the prior conformal first-order analysis that this paper improves on.\n\nWho is this for? People doing VAH phenomenology or causality constraints in heavy-ion hydrodynamics. For them it is worth a serious referee; I would send it to review and expect it to be accepted after the stability sentence is fixed and the NLO scope is made equally prominent in the abstract. If you sit on that kind of reading group, it is a reasonable 'maybe'—worth the hour, but not urgent. I would cite it if I were running VAH simulations.","headline":"Clean, correct derivation of the first necessary-and-sufficient nonlinear causality conditions for leading-order VAH; genuinely useful, but scoped to the truncated system and slightly overclaimed on stability.","tokens_in":18915,"tokens_out":4577,"would_cite":true,"duration_ms":38923,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that the leading-order viscous anisotropic hydrodynamics (VAH) equations are causally valid exactly when three simple inequalities on the anisotropic transport coefficients hold, and it proves these inequalities are both…","keywords":["viscous anisotropic hydrodynamics","nonlinear causality","characteristic velocities","relativistic heavy-ion collisions","quark-gluon plasma","pressure anisotropy","transport coefficients"],"falsifier":"Include the next-to-leading-order correction $\\delta\\tilde f$ in the equations of motion, recompute the characteristic equation, and search for parameter values satisfying the three inequalities that yield any characteristic speed $|v_c|>1$; a single such root would show the conditions do not certify causality of the complete theory. Alternatively, solve the truncated VAH equations as a Riemann problem with parameters just outside the allowed ranges and observe a signal propagating outside the light cone.","tokens_in":18108,"feed_emoji":"⚡","tokens_out":7432,"duration_ms":60690,"temperature":0.7,"pith_summary":"The paper claims that the leading-order viscous anisotropic hydrodynamics (VAH) equations are causally valid exactly when three simple inequalities on the anisotropic transport coefficients hold, and it proves these inequalities are both necessary and sufficient. These inequalities bound the longitudinal and transverse pressure-relaxation coefficients by the local energy and pressure scales, and impose a chain condition on their cross-coupling. The result matters because VAH is the framework proposed for describing the strongly anisotropic, far-from-equilibrium quark-gluon plasma created in the earliest moments of heavy-ion collisions, a regime where causality is not automatically guaranteed. The conditions give a practical test: if these ratios stay in the allowed ranges, the nonlinear evolution equations preserve relativistic causality in every propagation direction.","feed_headline":"Three inequalities certify causality in anisotropic hydrodynamics","feed_subtitle":"In the leading-order equations used for heavy-ion collisions, the bounds mark the causal regime.","key_machinery":"The load-bearing object is the characteristic equation of the quasilinear system, $\\det(A^\\alpha\\xi_\\alpha)=0$, reduced in the local rest frame to $$$v_c^{4}$+\\big[\\Gamma^\\perp_\\perp(\\hat $a_1^{2}$+\\hat $a_2^{2}$)+\\Gamma^L_z\\hat $a_3^{2}$\\big]$v_c^{2}$+(\\Gamma^L_z\\Gamma^\\perp_\\perp-\\Gamma^L_\\perp\\Gamma^\\perp_z)(\\hat $a_1^{2}$+\\hat $a_2^{2}$)\\hat $a_3^{2}$=0.$$ Because the equation is quadratic in $v_c^2$, causality for every direction is equivalent to four graphical conditions on the parabola $f(v_c^2)$: nonnegative at $0$ and $1$, axis within $[0,1]$, and nonnegative discriminant. Solving those four conditions successively yields the three inequalities of Eq. (21). The angle parametrization $\\chi=\\hat a_3^2$ turns the quartic into $v_c^4+[\\Gamma^\\perp_\\perp(1-\\chi)+\\Gamma^L_z\\chi]v_c^2+(\\Gamma^L_z\\Gamma^\\perp_\\perp-\\Gamma^L_\\perp\\Gamma^\\perp_z)\\chi(1-\\chi)=0$, exposing two independent sound modes in the decoupled limit and a smooth interpolation between fixed endpoint velocities when the coupling is active.","core_discovery":"By writing the VAH equations of motion as a first-order quasilinear system $A^\\alpha(\\Psi)\\partial_\\alpha\\Psi=F(\\Psi)$ and computing the characteristic equation $\\det(A^\\alpha\\xi_\\alpha)=0$, the paper shows that all characteristic velocities lie between zero and the speed of light in every direction exactly when $$-1\\le \\frac{\\bar\\zeta^L_z}{E+P_L}\\le0,\\qquad -1\\le \\frac{\\bar\\zeta^\\perp_\\perp}{E+P_\\perp}\\le0,\\qquad 0\\le \\frac{\\bar\\zeta^L_\\perp}{E+P_\\perp}\\frac{\\bar\\zeta^\\perp_z}{E+P_L}\\le \\frac{\\bar\\zeta^L_z}{E+P_L}\\frac{\\bar\\zeta^\\perp_\\perp}{E+P_\\perp}\\le 1.$$ In the decoupled limit the system separates into two sound-like modes with squared speeds $v_\\perp^2=-\\Gamma^\\perp_\\perp(1-\\chi)$ and $v_L^2=-\\Gamma^L_z\\chi$; the first two inequalities enforce causal boundary speeds, and the third extends the constraint to all mixing angles. The authors interpret the coupled system as two fluids exchanging pressure fluctuations, with $\\Gamma^L_\\perp\\Gamma^\\perp_z$ acting as a diffusion coupling.","pith_inferences":["If the next-to-leading-order correction $\\delta\\tilde f$ is not negligible, the three inequalities are only a necessary condition for causality of the complete VAH theory; adding transverse shear and longitudinal diffusion modes could introduce additional characteristic speeds that the current conditions do not constrain.","The inequality chain has a positivity-and-determinant structure, suggesting it may also imply linear stability of the anisotropic rest state, though that connection is not proven in the paper.","One could test the conditions numerically by running shock-type problems with parameters straddling the boundary of the inequalities and looking for superluminal precursors or unstable growth.","Because the conditions depend only on local ratios $\\bar\\zeta/(E+P)$, they can be mapped into spacetime for realistic heavy-ion collision simulations, potentially identifying when and where an anisotropic hydrodynamic description first becomes causally admissible."],"forward_implications":["For any parameter set satisfying the three inequalities, numerical solutions of the VAH equations remain within the relativistic causal region, so the conditions can be checked locally during a simulation to certify each cell.","Because the inequalities are necessary and sufficient, parameters outside the allowed range will produce at least one propagation direction with characteristic speed below zero or above the speed of light, marking the boundary of VAH's validity.","The simple algebraic form makes it feasible to impose causality as a constraint in parameter estimation or model calibration of anisotropic transport coefficients.","In the decoupled limit, the two eigenmodes provide a diagnostic: causality of the whole system collapses to causality of the longitudinal and transverse modes separately plus a coupling bound.","The results set the stage for extending the analysis to the complete VAH theory including the next-to-leading-order correction $\\delta\\tilde f$, which the authors identify as the immediate next step."],"supporting_citations":[{"why":"Defines the VAH decomposition and the equations of motion (5)-(9) whose causality is analyzed.","marker":"[60–63, 66, 68]"},{"why":"Supplies the nonlinear characteristic-velocity method and the notion of nonlinear causality constraints that this paper applies to VAH.","marker":"[86]"},{"why":"Earlier causality analysis for conformal first-order anisotropic hydrodynamics, contrasted with the present second-order VAH result.","marker":"[74]"},{"why":"Establishes the linear-regime causality framework for dissipative relativistic fluids that the nonlinear analysis extends.","marker":"[81, 82]"},{"why":"Provides the second-order relaxation structure with relaxation times $\\tau_\\pi$ and $\\tau_\\Pi$ used in the pressure relaxation equations.","marker":"[80]"},{"why":"Shows how nonlinear causality constraints are used to restrict hydrodynamic initial conditions, motivating the VAH parameter bounds.","marker":"[92]"}],"fun_headline_variants":["Causality in anisotropic hydro: three inequalities suffice","VAH causal only within these simple bounds","Exact causality conditions for viscous anisotropic hydro","Early-time hydro: causal limits in three inequalities"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof treats the truncated equations as the whole theory, assuming the omitted small correction term carries no faster-than-light signals.","fun_headline_variants_meta":{"raw":{"variants":["Causality in anisotropic hydro: three inequalities suffice","VAH causal only within these simple bounds","Exact causality conditions for viscous anisotropic hydro","Early-time hydro: causal limits in three inequalities"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000186,"raw_usage":{"total_tokens":1346,"prompt_tokens":988,"completion_tokens":358,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":604,"completion_tokens_details":{"reasoning_tokens":300}},"tokens_in":604,"tokens_out":358,"duration_ms":3691,"temperature":1.0,"reasoning_tokens":300,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T15:05:24.935272+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Include the next-to-leading-order correction $\\delta\\tilde f$ in the equations of motion, recompute the characteristic equation, and search for parameter values satisfying the three inequalities that yield any characteristic speed $|v_c|>1$; a single such root would show the conditions do not certify causality of the complete theory. Alternatively, solve the truncated VAH equations as a Riemann problem with parameters just outside the allowed ranges and observe a signal propagating outside the light cone.","supporting_citations":[{"cited_title":"Instability of Boost-invariant hydrodynamics with a QCD inspired bulk viscosity","cited_arxiv_id":"0805.0442","evidence_quote":"Supplies the nonlinear characteristic-velocity method and the notion of nonlinear causality constraints that this paper applies to VAH."}],"review_version":2}