{"id":"d0c4582d-6332-4c7e-8f43-a0e102411e75","arxiv_id":"2507.21552","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A Petrov-Galerkin discretization framework preserves discrete dissipation inequalities for a general class of energy-based models, including circuits, Cahn-Hilliard, and doubly nonlinear diffusion.","lead":"This paper builds a single numerical recipe for simulating a broad class of physics models that conserve or dissipate energy. The recipe preserves a discrete version of the energy balance in time, space, and model reduction, and comes with code for circuits, phase separation, and degenerate diffusion.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Proposition 3.6's dissipation identity is correct, but the scheme's central claim is conditional on unproved solvability of the per-interval nonlinear systems; a concrete existence check should be added.","rationale":"I re-derived the proof of Proposition 3.6 step by step. Lemma 3.1 is valid because ∂tz2τh lies in Pτk−1(Xh2), the X2-Z2-X2* Gelfand triple makes the dual pairing agree with the Z2 inner product, and the orthogonal projections Πh2 and Πτ2 are idempotent on the relevant subspaces. The localization to [tj−1,tj] is justified because Wτh2 consists of piecewise polynomials without continuity constraints, so the global projection decomposes over intervals. The test function (∂tz1τh, Πτ2Πh2∂z2H, z3τh) is in Wτh, and the chain rule applies because H is continuously Fréchet differentiable and z1τh,z2τh are continuous piecewise polynomials. The same argument works for the quadrature variant because positive weights preserve the sign of −⟨r(ζ),ζ⟩. Thus the central dissipation claim is internally sound. The reader's weakest assumption, existence and uniqueness of discrete solutions, is indeed the essential precondition: the paper assumes unique solvability of the continuous model (Section 2) and merely cites references for the discrete case, while listing well-posedness as future work (Section 5). Without solvability, the dissipation inequality is conditional to the point of vacuousness. I also noted the known constant error in Proposition 3.14 and the unproven convergence orders, but these are secondary to the main structure-preservation claim and do not change the verdict. The conditional acceptance remains appropriate, and the proposed existence check would either close the gap or show exactly where the scheme requires additional assumptions.","tokens_in":31097,"tokens_out":12382,"duration_ms":144215,"concrete_test":"For the semi-explicit index-1 case of Proposition 3.14 (system (3.15)) with polynomial degree k=1, write out the nonlinear system that Scheme 3.2 produces on the first time interval. Prove, via the implicit function theorem or a contraction argument, that for all sufficiently small τ there is a unique solution in a ball around the extrapolated previous-step state, using det(J33−R33)≠0 and Lipschitz continuity of ∇H. If the Jacobian at the previous state is singular, or the contraction constant fails to vanish as τ→0, the existence assumption is violated even in this benign subclass, and Proposition 3.6 must be restricted to cases where solvability is established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central identity in Proposition 3.6 is verified correctly: Lemma 3.1's projection identity holds locally on each time interval; choosing φτh = (∂tz1τh, Πτ2Πh2∂z2H, z3τh) is legitimate because all components lie in Wτh; j(ζ)·ζ = 0 and r(ζ)·ζ ≥ 0 yield the dissipation inequality. The quadrature version (3.21)-(3.23) also works because the weights are positive and the left-hand side remains an exact integral. The weak point is not the algebra but the precondition: Scheme 3.2 is stated for 'the solution zτh', and Section 3.3 defers existence/uniqueness of the resulting per-interval nonlinear algebraic systems to references [AF25, EG21b, HW18] rather than proving it. If for some admissible combination of spaces, time step, and consistent initial data the interval system has no solution, then Proposition 3.6 is vacuous for that case and the framework does not deliver a structure-preserving method in practice. The paper's own Section 5 lists well-posedness of the scheme as future research, confirming this gap. Since the abstract and conclusions advertise the scheme as structure-preserving without this caveat, the unproved solvability is the most load-bearing concern about the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a unified Petrov–Galerkin discretization framework for a recently introduced class of energy-based models, covering temporal discretization, spatial discretization, and model order reduction. The central claim is that all schemes in the framework preserve the continuous dissipation inequality exactly, without requiring symplecticity. The key technical tool is the projection identity of Lemma 3.1, used to define a projected gradient that is then inserted into the discrete equations. The authors prove the dissipation identity in Proposition 3.6, discuss special cases (space discretization, model reduction, quadrature variants), and support the theory with numerical experiments for a nonlinear circuit, the Cahn–Hilliard equation, and a doubly nonlinear parabolic equation. Code is publicly available.","tokens_in":31376,"tokens_out":6054,"duration_ms":68243,"significance":"If the central claims hold, the paper provides a useful unified viewpoint on structure-preserving discretization for a broad class of dissipative systems, including differential-algebraic and infinite-dimensional settings. The framework is attractive because it avoids mixed finite element tailoring and symplectic integrators, and it naturally includes model order reduction. The numerical experiments and the availability of code are strengths. The core dissipation identity in Proposition 3.6 is a short and essentially correct argument from the projection identity, and the quadrature extension is also sound. However, the paper leaves two load-bearing points open: the existence of solutions to the per-interval nonlinear algebraic systems is not proved, and the error bound in Proposition 3.14 contains a power-of-C2 algebra error. These issues do not invalidate the central dissipation argument but do affect the strength of the advertised claims.","major_comments":[{"comment":"Scheme 3.2 is stated as 'find zτh' and Proposition 3.6 is conditional on the existence of such a solution, but the paper does not prove that the per-interval nonlinear systems are solvable. Section 3.3 refers only to related approaches [AF25, EG21b, HW18], and Section 5 explicitly lists well-posedness of the scheme as future research. This is a load-bearing gap because the abstract and conclusions claim that the scheme is dissipation-preserving without this caveat. I recommend stating the existence assumption explicitly in the main theorems and either proving a local existence result under natural monotonicity hypotheses or clearly labeling the well-posedness as an open condition that the user must verify.","section":"Section 3.3 and Section 5"},{"comment":"The proof derives the pointwise bound ∥z3(t)−ez3(t)∥ ≤ C2(∥(I−Πh2)∂z2H(z12(t))∥+K∥ed(t)∥) and then concludes a squared L2 bound involving C2 to the first power. Squaring the pointwise inequality yields a factor C2², not C2. As written, the stated error bound for the algebraic variable is not established by the proof. This is a localized algebra error in a secondary but advertised contribution; it can likely be repaired by replacing C2 with C2² and adjusting constants, but the current statement is incorrect.","section":"Proposition 3.14 and its proof"}],"minor_comments":[{"comment":"In the displayed error bound, the term ∥(I−Πh2)∂z2H(z1,z2)∥ is written with norm L2(0,T;R^{n1+n2}), but the projection error (I−Πh2)∂z2H takes values in R^{n2}; the norm should be L2(0,T;R^{n2}).","section":"Proposition 3.14 statement"},{"comment":"The caption says 'left, middle' for the convergence histories, but the figure has panels (a), (b), and (c); please refer to the panels explicitly to avoid ambiguity.","section":"Figure 4.2 caption"},{"comment":"The text notes that sQ = 2k is used for the Cahn–Hilliard example because it improves Newton convergence, whereas Section 4.2 uses sQ = k+1. A brief general remark on how sQ is selected in practice would improve reproducibility.","section":"Section 4.3"},{"comment":"Reference [GM26] is listed with the arXiv identifier 2603.08656, which appears to contain an impossible date; please verify the identifier.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The central dissipation-preservation argument is sound and the paper makes a useful contribution. The main reasons for major revision are the missing discrete solvability proof/reference and the C2-vs-C2² error in Proposition 3.14, both of which are localized and likely fixable. I would encourage the authors to state the well-posedness condition as an explicit hypothesis and to correct the error bound; after those changes the paper could be acceptable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Colleague,\n\nThe short version: this paper is a credible and useful contribution. The central claim—that the proposed Petrov–Galerkin scheme preserves a dissipation inequality for time, space, and model reduction in one framework—is correct. The caveats are a solvability gap and a small constant error in one secondary bound.\n\nWhat's new: the infinite-dimensional formulation of the AS25 model class, the unified discrete dissipation inequality covering temporal and spatial discretization plus model reduction, the Dirac–Frenkel/nonlinear-manifold preservation result (Prop 3.9), and the semi-explicit index-one DAE reduction error bound (Prop 3.14). The paper is upfront that it combines existing ideas—projections from GKT25 and EHS21, reduction from CBG16—and it does so cleanly.\n\nProposition 3.6 is the heart. I checked the argument: the projection identity (3.3), the inclusion of the projected gradient terms in the test space, and the sign properties of j and r give the local dissipation inequality exactly. The quadrature variant (3.21)–(3.23) also works, since positive weights and an exact left side preserve the balance. That is a genuinely useful result for numerical analysts.\n\nThe main gap is solvability of the per-interval nonlinear systems. Scheme 3.2 is stated for 'the solution zτh', but Section 3.3 defers existence to references, and Section 5 lists well-posedness as future work. The dissipation inequality is therefore conditional on a solution existing. That is not unusual for nonlinear time-stepping schemes, but it should be stated honestly in the abstract and conclusions. As written, the scheme is advertised unconditionally.\n\nProposition 3.14 contains a small algebra slip: squaring the pointwise bound should produce C2², not C2, in the algebraic-variable error term. The proof strategy and qualitative form of the bound are fine; the constant is just wrong as printed.\n\nRelatedly, the observed convergence orders are empirical. The concluding 'arbitrary order of convergence' overstates what is proved; the authors do flag proofs as future research, so this is a wording issue rather than a hidden flaw.\n\nCredit where due: the code and data are available, the three numerical examples are run with care, and the implementation details (quadrature node counts, Newton settings, POD dimensions) are reported.\n\nWho this is for: numerical analysts working on structure-preserving time integration, port-Hamiltonian systems, gradient flows, or DAE model reduction. They will get a ready-made dissipation-preserving recipe and a clear comparison with prior approaches. It deserves a serious referee and likely publication after revision.\n\nMy recommendation: send it to peer review. The central result is sound, the framework is useful, and the fixes are routine: add a solvability caveat, correct the constant, and tone down the convergence claim to 'observed order'.","headline":"A sound and useful unification of dissipation-preserving discretization for energy-based models; core identity correct, but solvability gap and a constant error need attention.","tokens_in":31951,"tokens_out":3364,"would_cite":true,"duration_ms":37628,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["37J06","65P10","65M60"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that a single Petrov–Galerkin discretization framework for energy-based models preserves the continuous dissipation inequality exactly, simultaneously for temporal discretization, spatial discretization, and model…","keywords":["energy-based modeling","dissipation inequality","structure preservation","Petrov-Galerkin discretization","model order reduction","port-Hamiltonian systems","gradient flows","differential-algebraic equations"],"falsifier":"Take a model in the class (2.4), choose finite-dimensional subspaces and a time step for which Scheme 3.2 has a computable solution, evaluate the energy difference and the right-hand side of inequality (3.6) with quadrature exact to machine precision, and check whether the inequality holds; a single violation would refute the central claim. A second decisive experiment is to exhibit any model, ansatz spaces, and time step in the stated generality for which the nonlinear system of Scheme 3.2 has no solution, because the dissipation statement is vacuous there.","tokens_in":30849,"feed_emoji":"⚡","tokens_out":9998,"duration_ms":104676,"temperature":0.7,"pith_summary":"This paper claims that a single discretization recipe, a Petrov–Galerkin ansatz with piecewise polynomials in time and suitable projections, preserves the energy balance of a broad class of energy-based models exactly, in the sense of a discrete dissipation inequality. The class includes port-Hamiltonian systems, generalized gradient flows, and differential-algebraic systems with algebraic constraints. The recipe covers temporal discretization, spatial discretization, and linear model order reduction, and the paper proves the dissipation inequality in each setting. Numerical tests on a nonlinear circuit, the Cahn–Hilliard equation, and a doubly nonlinear parabolic equation show the discrete energy balance holding to machine precision. If correct, this gives one generic way to keep simulated energy dissipation physical without relying on symplectic integrators.","feed_headline":"Petrov–Galerkin scheme preserves dissipation in every discretization","feed_subtitle":"Covers port-Hamiltonian systems, gradient flows, and differential-algebraic models without symplectic integrators","key_machinery":"The load-bearing object is the projected Petrov–Galerkin discretization of the energy-based model (2.4), with piecewise-polynomial ansatz spaces in time and test spaces whose first two blocks are time derivatives of the ansatz spaces ($W_i^{\\tau h}=\\partial_t V_i^{\\tau h}$ for $i=1,2$). The second gradient $\\partial_{z_2}H$ is replaced by its combined orthogonal projection $\\Pi_2^{\\tau}\\Pi_2^{h}$ onto the corresponding test space, while the $z_1$-gradient needs no explicit projection because it is tested against $\\partial_t V_1^{\\tau h}$. Lemma 3.1, the identity $\\int\\langle\\partial_t v,w\\rangle=\\int\\langle\\partial_t v,\\Pi_2^{\\tau}\\Pi_2^{h} w\\rangle$ for $v$ in the ansatz space, is what turns the chain-rule energy difference into the right-hand side of the scheme and yields the dissipation inequality. The operators $j$ and $r$, which satisfy $\\langle j(v),v\\rangle=0$ and $\\langle r(v),v\\rangle\\ge0$, carry the conservative and dissipative parts of the balance.","core_discovery":"The central assertion, Proposition 3.6, is that for any solution $z^{\\tau h}=(z_1^{\\tau h},z_2^{\\tau h},z_3^{\\tau h})$ of Scheme 3.2, the energy difference over one time interval equals the integral of $-\\langle r(\\zeta),\\zeta\\rangle + \\langle y^{\\tau h},u\\rangle$, where $\\zeta$ is built from $\\partial_t z_1^{\\tau h}$, the projected gradient of $H$, and $z_3^{\\tau h}$. Since $r$ is dissipative, with $\\langle r(v),v\\rangle\\ge 0$, this gives the discrete dissipation inequality on every interval: the energy at the end minus the energy at the start is at most the supplied power $\\int \\langle y^{\\tau h},u\\rangle$. This holds without any symplectic structure; it follows from choosing the test space in the second block to be the time derivative of the ansatz space and inserting the orthogonal projection of $\\partial_{z_2}H$. The same mechanism is shown to preserve dissipation for the space-discrete scheme, the reduced-order model, and the quadrature variant, and the paper adds an error bound for a semi-explicit index-one reduced model.","pith_inferences":["A natural reading of Proposition 3.6 is that the temporal projection is the only genuinely new ingredient: any spatial discretization that supplies a well-defined $\\partial_{z_2}H$ can be made dissipation-preserving by adding this one projection, so the method should transfer readily to existing finite-element or neural-network codes.","The paper's separation of the state into $z_1,z_2,z_3$ suggests that the same projection trick applies componentwise in more general energy-based couplings, such as multi-physics or network interconnections, once the conservative and dissipative splitting (2.3) is identified.","Because only the dissipative part $r$ enters the correction term in the energy balance, one could test the method's robustness by deliberately mis-specifying $r$ slightly and checking whether the inequality still bounds energy growth; the paper does not perform such a perturbation study.","Proposition 3.9's nonlinear-manifold version gives state-dependent $j,r$; extending Proposition 3.14's error bound to that nonlinear setting would be a direct next step, but is not proved here."],"forward_implications":["Every solution of Scheme 3.2 reproduces the continuous dissipation inequality on each time interval, so simulations of the covered models cannot create energy from nothing.","The space-discrete coefficient system can be rewritten in the original model form (2.2), so finite-element implementations inherit the energy balance without tailored mixed methods.","Galerkin model reduction with orthogonal basis matrices produces a reduced system of the same form, hence a reduced model that dissipates exactly as the full model does.","For semi-explicit index-one systems, the reduced-order state error is bounded by projection errors of the state and of $\\partial_{z_2}H$, which makes the basis choice and the gradient approximation the two sources of error.","The quadrature variant keeps the inequality for any positive quadrature weights, while the equality holds only when the quadrature is exact."],"supporting_citations":[{"why":"Defines the finite-dimensional energy-based model class that Section 2 generalizes to infinite-dimensional systems (2.4).","marker":"[AS25]"},{"why":"Supplies the energy-consistent Petrov–Galerkin time discretization and the temporal projection underlying Scheme 3.2.","marker":"[GKT25]"},{"why":"Provides the energy-stable discretization idea for Hamiltonian and gradient systems and the observation that the $z_1$-component needs no explicit projection.","marker":"[EHS21]"},{"why":"Introduces the gradient-projection technique for structure-preserving model reduction and the error-bound strategy adapted for Proposition 3.14.","marker":"[CBG16]"},{"why":"Presents the related structure-preserving time discretization via auxiliary variables that Scheme 3.2 recovers as a special case.","marker":"[AF25]"},{"why":"Underpins existence and uniqueness of the nonlinear time-stepping systems and the collocation interpretation in Remark 3.4.","marker":"[EG21b]"},{"why":"Contains the technical evaluation theory that makes the projected time-stepping scheme well-defined and localizes it to individual time steps.","marker":"[Kar26]"}],"fun_headline_variants":["No symplectic needed: dissipation preserved across all discretizations","Petrov–Galerkin choice preserves dissipation in time, space, and reduction","Discrete dissipation inequality: one framework for energy-based models","Unified Petrov–Galerkin scheme guarantees dissipation preservation","Energy-based models keep dissipation under all discretizations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim rests on the assumptions that the continuous model (2.4) has a unique solution and that the nonlinear algebraic systems produced by Scheme 3.2 have a solution for the chosen discrete spaces and time steps; without those solutions the discrete dissipation inequality cannot even be evaluated.","fun_headline_variants_meta":{"raw":{"variants":["No symplectic needed: dissipation preserved across all discretizations","Petrov–Galerkin choice preserves dissipation in time, space, and reduction","Discrete dissipation inequality: one framework for energy-based models","Unified Petrov–Galerkin scheme guarantees dissipation preservation","Energy-based models keep dissipation under all discretizations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000681,"raw_usage":{"total_tokens":3072,"prompt_tokens":904,"completion_tokens":2168,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":520,"completion_tokens_details":{"reasoning_tokens":2082}},"tokens_in":520,"tokens_out":2168,"duration_ms":20641,"temperature":1.0,"reasoning_tokens":2082,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T12:39:04.312833+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a model in the class (2.4), choose finite-dimensional subspaces and a time step for which Scheme 3.2 has a computable solution, evaluate the energy difference and the right-hand side of inequality (3.6) with quadrature exact to machine precision, and check whether the inequality holds; a single violation would refute the central claim. A second decisive experiment is to exhibit any model, ansatz spaces, and time step in the stated generality for which the nonlinear system of Scheme 3.2 has no solution, because the dissipation statement is vacuous there.","supporting_citations":[],"review_version":1}