{"id":"2e7f13dc-a60b-431c-850e-15ae40abac32","arxiv_id":"2506.06430","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The BDNK conformal viscous fluid's full constraint set satisfies a homogeneous strongly-hyperbolic system, guaranteeing constraint preservation, and 1D numerical evolutions confirm stability.","lead":"Researchers prove that the constraint equations of a modern first-order viscous fluid model (BDNK) propagate correctly, because the constraints obey a strongly hyperbolic system. They then show in 1D simulations that these constraints stay small over time, supporting the framework's use in numerical relativity.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The constraint-evolution system (21)-(27) is asserted via 'one can show' with only final source terms in Appendix A; if this hidden algebra is wrong, the strong-hyperbolicity proof targets the wrong system and constraint preservation does not follow.","rationale":"The reader's weakest_assumption identifies exactly the load-bearing premise: the constraint-evolution system (21)-(27) is presented as the result of an unshown algebraic reduction. Everything downstream—the characteristic polynomial (30), the eigenspace count of 29+1+1, and the conclusion that zero initial constraints propagate—is conditional on that system being the true one. The hyperbolicity calculation itself is plausible: the block structure of the principal part does yield the factorized polynomial U^29(U^2 - (η/λ) P(k,k)) if (21)-(27) is accepted, and the null-space dimension in Appendix B checks out under that assumption. But the bridge from the fluid equations to (21)-(27) is not checkable from the manuscript as written. The numerical evolutions and convergence tests in Section V are real supporting evidence for the behavior of two differential-constraint components in a plane-symmetric setting, but they cannot validate the full covariant derivation or the algebraic-constraint sector. Because the missing derivation is the same concern the reader flagged, my read does not change the verdict: conditional acceptance pending an independent verification of (21)-(27).","tokens_in":22860,"tokens_out":34497,"duration_ms":317799,"concrete_test":"Use a symbolic algebra system to rederive (21)-(27) from the definitions (14)-(20) and the evolution equations (7)-(12). For each of the 31 constraints C_A, form the residual u^a ∇_a C_A plus the principal terms displayed in (21)-(27), substitute the definitions of A, Q^a, S^a, S^a_b, θ, u^a, and eliminate all first derivatives of the dynamical fields using Eqs. (7)-(12). The simplified residual must equal -v_A as given in Appendix A; a mismatch in any of the 31 components means the hyperbolicity analysis applies to a different system. The check should also force a resolution of the index ambiguity in Eq. (A14). If the identity is verified for all components, the central analytic claim is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim—that constraints (14)-(20) propagate because (21)-(27) is strongly hyperbolic—rests on two parts: the hyperbolicity calculation in Section III.B, and the derivation of (21)-(27) itself. The derivation is the weaker link. Section III.A states 'one can show' but displays no intermediate steps; only the final source terms v_i are given in Appendix A. As written, those source terms are not fully self-consistent: in Eq. (A14), for instance, the first two terms u^a S_b C3 and S_b C5^a carry a free lower index b, whereas the left-hand side v7^{ab} requires an upper index b. If the true constraint system differs from (21)-(27) in any principal term, or if the zero-constraint limit of the sources is not exactly the stated v_i, then the characteristic-polynomial and eigenspace-count arguments prove strong hyperbolicity for a different system, and the conclusion that zero initial constraints remain zero does not follow. The numerical section does not close this gap: it only monitors two plane-symmetric components of C6 and C7 (Section V.B), not the full 31-component system.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the constraint-propagation problem for the first-order reduction of conformally invariant BDNK viscous relativistic fluids introduced in [65]. The authors identify a 31-component constraint system consisting of the algebraic constraints (14)-(18) and the differential constraints (19)-(20). They claim that these constraints obey a homogeneous first-order quasi-linear evolution system (21)-(27), that the principal part of this system is strongly hyperbolic with a 29-fold eigenvalue and two simple real eigenvalues, and that therefore vanishing initial constraints remain zero throughout the evolution. The paper then presents plane-symmetric numerical evolutions of the BDNK equations, monitoring two components of the differential constraints plus convergence tests in space and time.","tokens_in":23046,"tokens_out":15010,"duration_ms":146030,"significance":"If the analytic result is correct, it fills a real gap in the BDNK program: local well-posedness was proved in [65] without explicitly analyzing propagation of the constraints needed for the first-order reduction, and the present paper would supply exactly that missing ingredient. The hyperbolicity computation is compact and parameter-free, depending only on the BDNK inequalities (13), and the numerical experiments are a reasonable first exploration of constraint stability. The main reservation is that the derivation of the constraint evolution system is not actually displayed, and a concrete index inconsistency appears in the printed source terms; the central theorem is therefore currently supported by an unverified algebraic reduction.","major_comments":[{"comment":"The central claim of the paper rests on the assertion that the constraints (14)-(20) satisfy the evolution system (21)-(27), but the derivation is not shown. Section III.A states 'one can show' and refers to Appendix A for the source terms, without displaying any intermediate steps. Since the strong-hyperbolicity analysis in Section III.B is applied to the principal part of (21)-(27), any error in this reduction would mean that the proof targets a system different from the true constraint system. This is a load-bearing step, not a presentation detail. In addition, Eq. (A14) as printed is not index-consistent: the first two terms, u^a S_b C3 and S_b C5^a, carry a free lower index b, whereas the left-hand side v7^{ab} requires an upper b. The authors should provide a complete derivation or a computer-algebra verification of (21)-(27) and correct the source terms.","section":"Section III.A and Appendix A"},{"comment":"The proof of strong hyperbolicity depends on the claim that the eigenvalue w29 has a 29-dimensional eigenspace, but the verification in Appendix B is not usable as written. The displayed block matrix mixes blocks labeled '07×7' and '024×24' in a way that does not assemble into a 31×31 matrix, and the list of spanning vectors does not obviously contain 29 vectors. Since the completeness of the eigenspaces is exactly what distinguishes strong hyperbolicity from mere weak hyperbolicity, this computation needs to be rewritten cleanly so that the nullity of N^A_B is actually verified for every relevant covector.","section":"Section III.B and Appendix B"},{"comment":"The numerical section monitors only the two plane-symmetric components CQ1 and CS11 of the differential constraints C6 and C7. It does not follow the full 31-component constraint system, nor does it monitor the algebraic constraints C1-C5 directly. The text does explain in Section IV.A that the algebraic constraints are used to reduce the number of evolution variables, so the omission may be acceptable for an exploratory study, but the paper should state explicitly that the numerical evidence is partial and is not a substitute for the analytic constraint-propagation proof.","section":"Section V.B, Eqs. (36)-(37)"}],"minor_comments":[{"comment":"The sentence 'p(w) \\equiv 0 if and only if ...' should read 'p(w) = 0 if and only if ...'; the identical-equality symbol is not what is meant.","section":"Section III.B, Eq. (30)"},{"comment":"The numerical convergence tests in Figure 4 and the surrounding text are reported as supporting second-order convergence of the constraints, but the figure shows only the difference between resolutions, not the ratio itself; adding the measured convergence factor would make the claim easier to verify.","section":"Section IV.C and V.C"},{"comment":"The relation between the inequality (13) and the condition η < λ used in Section III.B is asserted but not derived; a short derivation would help the reader see that (13) indeed implies a2 > 1 and hence η < λ.","section":"Section II, Eq. (13)"},{"comment":"The uniform-velocity argument in Appendix C is a separate application and is not needed for the main theorem; if kept, it should be clearly marked as a physical consequence rather than as part of the proof of constraint preservation.","section":"Appendix C"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses an important question and the hyperbolicity calculation appears plausible, but the unshown derivation of (21)-(27) and the index inconsistency in (A14) are serious gaps in the current manuscript. If the authors can supply a verified derivation (ideally with a supplementary Mathematica notebook or an appendix with the intermediate algebra), the result would be a useful contribution. The editorial bar should be that the constraint system whose principal part is analyzed is actually the constraint system induced by the BDNK reduction."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nThe short version: this paper gives the first full constraint-propagation analysis for the BDNK conformal first-order reduction. That is a genuinely missing piece, and the hyperbolicity argument looks credible: the characteristic polynomial has a 29-fold eigenvalue plus two simple real eigenvalues, the eigenvector count in Appendix B supports diagonalizability, and the parameter condition (η < λ) sits inside the BDNK inequalities. If the constraint evolution system (21)-(27) is correct, the main conclusion follows.\n\nWhat is actually new and good: the paper states the complete set of algebraic and differential constraints, derives a homogeneous evolution system for them, proves strong hyperbolicity, and supports the result with plane-symmetric numerical evolutions including convergence tests for the two constraints they monitor. The appendix on uniform-velocity configurations is a nice bonus.\n\nWhere it is soft: the derivation of (21)-(27) is the load-bearing step, and it is not shown. Section III.A says 'one can show' and the reader is left with the final source terms in Appendix A. The stress-test note is right: if that reduction is off in any principal term, the strong-hyperbolicity proof targets the wrong system. I also confirm the index issue in Eq. (A14): the first two terms have a lower b, while v7^{ab} needs an upper b. That is likely a typo, but it is exactly the sort of thing that should have been caught, and it makes it hard to trust the appendix as-is. The numerics don't close the gap: they only monitor two components of the differential constraints, not the full 31-component system, and they are smooth 1D data. The convergence tests are good but not a substitute for the algebra.\n\nWho benefits: anyone working on BDNK theory or first-order viscous fluids in numerical relativity. This paper deserves a serious referee. I recommend sending it to review with the request that the authors provide the derivation steps (or a machine-checkable notebook) and fix the index errors. My own verdict would be conditional acceptance, not rejection.\n\nBest,\n[Name]","headline":"First full constraint-propagation analysis for the BDNK conformal reduction, with a plausible strong-hyperbolicity proof, but the central derivation is asserted and the source terms have an index inconsistency that needs fixing.","tokens_in":23584,"tokens_out":3373,"would_cite":true,"duration_ms":31736,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35L45","35L60","35Q75","76Y05"],"pacs":["47.75.+f"],"model":"deepseek-v4-flash","headline":"All constraints of the BDNK first-order viscous fluid reduction satisfy a strongly hyperbolic homogeneous system, so initially satisfied constraints persist for all times, and simulations stay stable.","keywords":["relativistic viscous hydrodynamics","BDNK first-order theory","constraint propagation","strong hyperbolicity","well-posedness","conformal fluids","numerical relativity","algebraic and differential constraints"],"falsifier":"Recompute the time derivatives of the 31 constraints directly from (7)-(12) and check whether the principal part and the linear-or-quadratic structure of the source terms match (21)-(27) with the stated $v_i$; any omitted term that is not homogeneous, or any change in $M_a{}^A{}_B$, would invalidate the conclusion. Numerically, seed initial data that satisfy the fluid equations but violate the constraints by a small perturbation and look for exponential growth of the constraint residuals rather than bounded, stationary behavior.","tokens_in":22629,"feed_emoji":"🌊","tokens_out":9865,"duration_ms":93519,"temperature":0.7,"pith_summary":"This paper fills a gap in the theory of first-order relativistic viscous fluids: whether the algebraic and differential constraints used to reduce the BDNK equations to first order remain satisfied once they hold at one instant. The authors show that these constraints obey their own homogeneous evolution system and prove that this system is strongly hyperbolic, which makes the identically zero constraint configuration the unique solution for zero initial data. They then evolve plane-symmetric configurations in flat spacetime and report that the constraint residuals stay small and stable, settling near zero after a short transient. If correct, the result turns the BDNK first-order reduction into a sound starting point for numerical simulations of dissipative relativistic fluids, not just a well-posedness proof.","feed_headline":"Viscous fluid constraints survive the whole evolution","feed_subtitle":"First-order BDNK fluid keeps all algebraic and differential constraints, and plane-symmetric runs stay stable.","key_machinery":"The central object is the $31\\times 31$ principal symbol $M_a{}^A{}_B$ of the constraint system and its characteristic polynomial $p(w)=\\det\\big((w_a-w t_a)M_a{}^A{}_B\\big)$. Its factorisation yields the repeated eigenvalue $w_{29}=u_a w^a/(u_c t^c)$ and the pair $w_{\\pm}$; the load-bearing checks are the reality condition $E^2_{wt}>E_{ww}E_{tt}$, equivalent to $\\eta<\\lambda$ and consistent with the stability conditions on the transport coefficients, and the 29-dimensional eigenspace computed in the appendix. Strong hyperbolicity of the homogeneous system is what upgrades 'zero is a solution' into 'zero is the unique solution with zero initial data'.","core_discovery":"The paper's central claim is that the complete set of constraints (14)-(20) of the BDNK first-order reduction --- the unit-norm and orthogonality conditions on the four-velocity, plus the differential definitions of heat flux and velocity-gradient variables --- forms a homogeneous first-order quasi-linear system (21)-(27). The paper computes the $31\\times 31$ principal part and shows that the characteristic polynomial factorises as $[u_a(w^a-w t^a)]^{29}\\{[u_a(w^a-w t^a)]^2-(\\eta/\\lambda)\\bar{\\Pi}_{ab}(w^a-w t^a)(w^b-w t^b)\\}$, giving eigenvalues $w_{29}$ and $w_{\\pm}$, all real when $\\eta<\\lambda$, with eigenspaces of dimensions 29, 1, and 1 that span the whole space. Strong hyperbolicity plus homogeneity means that the identically zero constraint configuration is the unique solution for zero initial data, so the fluid equations themselves preserve the constraints. Numerical evolutions of plane-symmetric smooth data confirm that the two differential constraints, measured by $L^2$ norms, settle to stationary values close to zero after a short transient.","pith_inferences":["The homogeneity of the constraint system leaves open a practical route the paper does not pursue: adding constraint-damping terms to (21)-(27) that are lower order in the constraints would not alter the principal part, so it could suppress numerical violations without giving up hyperbolicity.","The factorisation of the characteristic polynomial suggests that the ratio $\\eta/\\lambda$ controls the geometry of constraint propagation; testing the boundary case $\\eta=\\lambda$ could reveal whether the 29-dimensional eigenspace degenerates.","Since only plane-symmetric, flat, smooth configurations are simulated, a natural next test is to seed small constraint violations in 3D or curved spacetimes and check whether the $L^2$ norms still settle toward zero."],"forward_implications":["Constraint preservation becomes a consequence of the fluid equations themselves: initial data satisfying (14)-(20) continue to satisfy them for all times, so no extra constraint projection step is needed during evolution.","The local well-posedness result for the conformal BDNK system in Sobolev spaces is reinforced by showing that the full constraint set is itself locally well-posed.","Uniform-velocity configurations with non-constant temperature are ruled out for conformally invariant fluids.","Plane-symmetric numerical evolutions with smooth Gaussian data show constraint residuals whose $L^2$ norms remain bounded and settle near zero after a short transient.","Convergence tests show the expected second-order spatial scaling for the constraints, indicating that the observed preservation is consistent with an accurate discretization."],"supporting_citations":[{"why":"Supplies the first-order reduction, its evolution equations (7)-(12), and the local well-posedness result whose constraint propagation this paper completes.","marker":"[65]"},{"why":"Provides the classical existence-uniqueness technique used in the original well-posedness proof, the context that makes constraint preservation necessary.","marker":"[59]"},{"why":"Introduces the BDNK first-order theory of relativistic viscous fluids with causality and existence results that the conformal reduction builds on.","marker":"[40]"},{"why":"Establishes linear stability of first-order theories under transport-coefficient conditions, including the parameter restrictions used here ($\\eta<\\lambda$).","marker":"[43]"},{"why":"Documents instability and acausality of earlier first-order formulations, the baseline problem that the BDNK framework and this constraint analysis address.","marker":"[33]"},{"why":"Discusses strong hyperbolicity and constraint preservation for constrained first-order systems, the framework in which the paper's hyperbolicity check is posed.","marker":"[62]"}],"fun_headline_variants":["Constraints survive the whole viscous fluid evolution","Viscous fluid constraints prove stable and self-consistent","BDNK fluid constraints hold under numerical simulation","First-order viscous fluid keeps every constraint intact","Hyperbolic constraints self-preserve in viscous fluid runs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing step is the algebraic reduction that turns the fluid equations (7)-(12) into the constraint evolution system (21)-(27); the paper states that 'one can show' this and gives only the final source terms, so if that reduction is wrong, the strong-hyperbolicity proof would apply to a different system.","fun_headline_variants_meta":{"raw":{"variants":["Constraints survive the whole viscous fluid evolution","Viscous fluid constraints prove stable and self-consistent","BDNK fluid constraints hold under numerical simulation","First-order viscous fluid keeps every constraint intact","Hyperbolic constraints self-preserve in viscous fluid runs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00019,"raw_usage":{"total_tokens":1355,"prompt_tokens":980,"completion_tokens":375,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":596,"completion_tokens_details":{"reasoning_tokens":305}},"tokens_in":596,"tokens_out":375,"duration_ms":4444,"temperature":1.0,"reasoning_tokens":305,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:57:39.374036+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the time derivatives of the 31 constraints directly from (7)-(12) and check whether the principal part and the linear-or-quadratic structure of the source terms match (21)-(27) with the stated $v_i$; any omitted term that is not homogeneous, or any change in $M_a{}^A{}_B$, would invalidate the conclusion. Numerically, seed initial data that satisfy the fluid equations but violate the constraints by a small perturbation and look for exponential growth of the constraint residuals rather than bounded, stationary behavior.","supporting_citations":[{"cited_title":"On constraint preservation and strong hyperbolicity,","cited_arxiv_id":null,"evidence_quote":"Supplies the first-order reduction, its evolution equations (7)-(12), and the local well-posedness result whose constraint propagation this paper completes."},{"cited_title":"Initial boundary value problems for hy- perbolic systems,","cited_arxiv_id":null,"evidence_quote":"Provides the classical existence-uniqueness technique used in the original well-posedness proof, the context that makes constraint preservation necessary."},{"cited_title":"Fixing extensions to general relativity in the nonlinear regime,","cited_arxiv_id":null,"evidence_quote":"Introduces the BDNK first-order theory of relativistic viscous fluids with causality and existence results that the conformal reduction builds on."},{"cited_title":"Hyperbolic theory of relativistic conformal dissipative fluids,","cited_arxiv_id":null,"evidence_quote":"Establishes linear stability of first-order theories under transport-coefficient conditions, including the parameter restrictions used here ($\\eta<\\lambda$)."},{"cited_title":"Lectures on hydrodynamic fluctuations in relativistic theories,","cited_arxiv_id":null,"evidence_quote":"Documents instability and acausality of earlier first-order formulations, the baseline problem that the BDNK framework and this constraint analysis address."},{"cited_title":"Choquet-Bruhat, C","cited_arxiv_id":null,"evidence_quote":"Discusses strong hyperbolicity and constraint preservation for constrained first-order systems, the framework in which the paper's hyperbolicity check is posed."}],"review_version":1}