{"id":"a6c4092e-8ebc-444f-a396-707f48924212","arxiv_id":"2511.13217","paper_version":2,"verdict":"REJECT","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"A claimed coercive variational principle for the Helmholtz equation fails: the exact impedance solution is not a stationary point of the proposed energies.","lead":"This paper proposes regularized energy functionals intended to make the Helmholtz equation with absorbing boundaries into a convex minimization problem. The central equivalence between the energy minimum and the true wave solution is invalid because boundary terms in the first variation are omitted.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The exact impedance Helmholtz solution is not stationary for Eγ or Fγ: the first variation retains a nonzero boundary term Re∫∂Ω iku v̄, so Theorems 2–3 characterize a different variational problem.","rationale":"The reader's weakest assumption is exactly the load-bearing gap: the stationarity conditions (35) and (41) are asserted to characterize the impedance Helmholtz solution without checking the boundary Euler–Lagrange term. My independent check confirms the concern. At an exact solution u, the first variation of the physical/regularized energy leaves a boundary term Re∫∂Ω iku v̄. This term vanishes only if one imposes an additional condition (effectively ∂nu=0 or a modified Robin condition) that is not part of (1) and is not enforced by V_BC. The counterexample v=iu is particularly clean because iu is always in V_BC if u is, and it isolates the boundary contribution with a strictly positive real part. Thus the minimizer of the proposed coercive energies is not the Helmholtz solution. The coercivity estimates in Section 4 may be correct for the modified variational problem, but they do not rescue the central claim. The paper's own Proposition 1 already hints at the gap by deriving only the interior equation from compactly supported tests, and the boundary condition is then assumed rather than derived. I therefore agree with the reader's reject verdict; no adjustment is needed.","tokens_in":22209,"tokens_out":8216,"duration_ms":76952,"concrete_test":"Take Ω=(0,1), k=1. Choose a smooth u (e.g. a quadratic) satisfying the Robin conditions −u′(0)=i u(0) and u′(1)=i u(1); set f=−u″−u. Then u∈V_BC and Lu=f. Compute the first variation of Eγ at u in direction v=iu. The difference D=A(u,iu)−∫f(iū+2γ L(iu)̄) has real part Re∫∂Ω ∂nu iū = |u(1)|²+|u(0)|² > 0, so (35) fails for the true solution. Repeating the check for Fγ gives the same positive boundary contribution. This is a direct, parameter-free counterexample to Theorem 2(ii) and Theorem 3(ii).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing gap is that the paper never verifies that the true solution of (1) satisfies the stationarity equations (35)/(41). For an exact solution u (Lu=f, ∂nu=iku), the first variation of Eγ at u in direction v∈V_BC is Re{A(u,v)−∫f(v̄+2γLv̄)}. By Green's identity this equals Re{∫∂Ω ∂nu v̄} = Re{∫∂Ω iku v̄}. This is not zero for all v∈V_BC: choosing v=iu, which lies in V_BC whenever u does, gives Re{∫∂Ω iku (−iū)} = k∥u∥²_{L²(∂Ω)} > 0. Hence u is not a stationary point of Eγ, and the same argument applies to Fγ because the boundary-penalty term vanishes on the exact solution. The derivation in Section 3.2 uses only compactly supported test functions to obtain the interior equation (28), then asserts that the impedance condition is 'encoded' in V_BC; but imposing ∂nv−ikv=0 on admissible variations does not make ∫∂Ω ∂nu v̄ vanish. Consequently, the unique minimizer of (30)/(38) solves a different boundary-value problem, and Lax–Milgram in Theorems 2–3 proves well-posedness of the wrong equation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes variational principles for the Helmholtz equation with impedance boundary conditions. Starting from a time-averaged Lagrangian, it derives an indefinite physical energy EP, then augments it with least-squares residual terms, obtaining Eγ over the impedance-constrained space VBC and Fγ with a weak boundary penalty over V. The main theoretical claim is that, for sufficiently large penalty parameters, these augmented energies are strongly coercive and their unique minimisers coincide with the weak solution of the Helmholtz problem (1). Coercivity is proved via Rellich/Morawetz identities, and the framework is illustrated with H²-conforming finite element and plane-wave neural network discretisations.","tokens_in":22511,"tokens_out":6323,"duration_ms":65326,"significance":"If the central equivalence were correct, the paper would offer a physically motivated, coercive variational formulation for the impedance Helmholtz problem with potential value for finite element and neural network methods. The coercivity estimates in Section 4 are nontrivial and point to a useful set of energy identities. However, the central claim fails: the true solution is not a stationary point of Eγ or Fγ because a boundary term survives integration by parts. This is a load-bearing mathematical error, not a presentation issue. The unique minimisers of the proposed energies solve a different boundary-value problem, so the main theorems and the numerical methods built on them do not deliver what is claimed. The paper is therefore not suitable for publication in its current form.","major_comments":[{"comment":"The stationarity condition (35) is not satisfied by the true solution u of (1). For u with Lu=f and ∂nu=iku, the first variation of Eγ at u is the same as that of the physical part, since the residual term vanishes. Using Green's identity, ⟨DEγ(u),v⟩ = Re{∫_Ω(∇u·∇v̄ - k²uv̄ - f v̄)} + boundary terms = Re{∫_∂Ω ∂nu v̄} = Re{∫_∂Ω iku v̄}. This does not vanish for all v∈VBC; choosing v=iu∈VBC gives Re{∫_∂Ω iku (-iū)} = k∫_∂Ω|u|² > 0. Thus u is not stationary, contradicting Theorem 2(ii)-(iii). The statement that the impedance condition is 'encoded in VBC' confuses the admissible test space with the boundary condition satisfied by the solution; the Euler-Lagrange derivation in Section 3.2 only tests with compactly supported functions and never checks the boundary contribution.","section":"Section 3.2, Eq. (33)-(35) and Proposition 1"},{"comment":"The same defect invalidates the weakly penalised energy Fγ. At the exact solution u, the boundary residual ∂nu-iku=0, so the γ2 term contributes nothing to the first variation; the physical part again leaves the boundary term Re{∫_∂Ω iku v̄}, which is nonzero for general v∈V. Hence the stationarity equation (41) is not equivalent to (1). In fact, for smooth solutions satisfying the interior equation, the boundary term would force a homogeneous Neumann-type condition rather than the impedance condition on the unconstrained space V. Consequently Theorem 3(ii)-(iii) is false, and the finite element and neural network minimisations in Section 5 target a different boundary-value problem.","section":"Section 3.2, Eq. (40)-(41) and Theorem 3"},{"comment":"The coercivity estimates (73) and (77) are statements about the quadratic parts Eγ(u)|_{f=0} and Fγ(u)|_{f=0}; at best they prove well-posedness of the variational problems (35) and (41). They do not establish that the unique solutions of those variational problems satisfy the impedance boundary condition. The paper's transition from coercivity to 'unique minimiser equals the solution of (1)' relies entirely on the false equivalence between the stationarity equations and problem (1). Since the true solution is not stationary for the proposed functionals, the coercivity results cannot rescue the central claim. The proofs of Theorem 2 and Theorem 3 therefore have a load-bearing gap at the point where Lax-Milgram is invoked.","section":"Section 4, Theorems 14 and 15"}],"minor_comments":[{"comment":"The real inner product is defined as ⟨φ,ψ⟩_R = Re∫_Ω φψ dx, but the first variation of |v|² in a complex Hilbert space should involve Re∫_Ω φ ψ̄ dx. Equations (25), (33) and (34) use unconjugated products; this inconsistency should be fixed, although the missing-boundary-term criticism above is independent of the conjugation convention.","section":"Section 3.1-3.2, Eq. (25) and Eq. (33)"},{"comment":"The proof of the low-order Morawetz identity appears to identify M0u with +ikβu at one step while the definition (56) sets M0u=-ikβu. Please check the sign consistency in the derivation; the final identity may still be correct, but the proof as written is confusing.","section":"Proposition 7, Eq. (57)"},{"comment":"The neural network experiments show snapshots of the field but no comparison against a reference solution or an exact solution. Since the variational problem solved is not the impedance Helmholtz problem, the numerical evidence does not support the claim that the method approximates (1).","section":"Section 5.4 and Figure 2"},{"comment":"The abstract states that the unique minimisers coincide with solutions of the Helmholtz equation, but Proposition 1 itself only proves the interior equation from compactly supported variations. The boundary condition is asserted without a stationarity calculation on the boundary. This gap should be acknowledged at the statement of Proposition 1.","section":"Abstract and Section 3.2"}],"recommendation":"reject","confidential_remarks":"The central mathematical claim is contradicted by the paper's own stationarity equations: a boundary term from integration by parts is omitted, and the true impedance solution is not a stationary point of Eγ or Fγ. This is not a matter of presentation or missing numerical validation; it invalidates Theorems 2 and 3 and the numerical methods derived from them. The coercivity machinery in Section 4 may be of independent interest, but it cannot justify the paper's main conclusions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The central claim of this paper is false. For the exact solution u of (1), the first variation of Eγ (and of Fγ, since the boundary penalty vanishes on u) reduces to Re{ik∫∂Ω u v̄}, which does not vanish for all admissible v. Taking v=iu gives k∥u∥²_{L²(∂Ω)}>0. So u is not a stationary point, and the unique minimizer of (30) or (38) solves a different boundary-value problem — one with an extra boundary source in the weak formulation. This invalidates Theorems 2 and 3 and the abstract's central statement.\n\nI want to give credit where it's due. The derivation of the physical energy from Hamilton's principle is clear, the Rellich/Morawetz identities are carefully developed, and the coercivity bounds for the quadratic parts of Eγ and Fγ — with explicit constants in Remark 17 — appear to be correct as statements about the functionals themselves. That part is a solid technical contribution, though the technique is not radically different from Moiola–Spence.\n\nThe soft spots are concentrated in the Euler–Lagrange analysis. The paper tests (33) only against compactly supported functions to obtain the interior equation (28), then asserts the impedance condition is 'encoded' in V_BC. But imposing ∂nv−ikv=0 on the variation v does not make ∫∂Ω ∂nu v̄ vanish. This is not a minor gap; it is the load-bearing step connecting the minimization problem to the Helmholtz equation. The numerical experiments cannot detect this because they are purely illustrative: no convergence tables, no reference solutions, only snapshots of the computed field. The novelty of the neural network formulation is also overstated — it is a least-squares residual formulation with a physics energy term.\n\nWho should read this? Researchers working on coercive formulations of Helmholtz would find the coercivity estimates useful, but only as a lemma, not as a variational principle. The paper deserves a serious referee because the authors are credible and the flaw is subtle enough to need expert scrutiny. My own verdict would be reject with a clear explanation of the missing boundary term; if the authors can modify the energy so that the exact solution is stationary (for instance by adding a suitable boundary term), the coercivity machinery might be salvageable.","headline":"The coercivity estimates look careful, but the central variational principle is false: a missing boundary term in the first variation means the exact Helmholtz solution is not stationary for Eγ or Fγ.","tokens_in":23046,"tokens_out":4980,"would_cite":false,"duration_ms":46770,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35J05","65N30","65N35"],"pacs":[],"model":"deepseek-v4-flash","headline":"Minimizing a regularized physical energy recovers the unique Helmholtz solution.","keywords":["Helmholtz equation","impedance boundary condition","variational principle","strong coercivity","Rellich identity","Morawetz identity","finite element method","physics-informed neural networks"],"falsifier":"Take the one-dimensional Helmholtz problem on an interval with a known solution, e.g. u = e^{ikx}, and a smooth test function w in V_BC that is not proportional to u. Evaluate the first variation of Eγ at u against w; if the surviving boundary contribution Re(∂n u ̄w)|_0^1 is nonzero, the exact solution is not a stationary point over V_BC, contradicting the theorem. A discrete version: solve the variational problem (35) with conforming H² elements for a manufactured solution and measure whether the boundary impedance residual ∥∂n u_h − iku_h∥_{L²(∂Ω)} tends to zero under mesh refinement.","tokens_in":22055,"feed_emoji":"🌊","tokens_out":10621,"duration_ms":89649,"temperature":0.7,"pith_summary":"The paper asks whether the Helmholtz equation with absorbing (impedance) boundary conditions can be derived from an energy principle. Starting from the least-action principle for the wave equation, the authors obtain the usual time-harmonic energy; it is sign-indefinite, so it cannot be minimized directly. The central claim is that adding least-squares penalties for the interior residual and for the impedance boundary condition, with sufficiently large weights, produces a strongly coercive energy whose unique minimizer is exactly the weak solution of the Helmholtz problem. Coercivity is proved with Rellich–Morawetz identities that control the indefinite bulk terms and boundary contributions, and the resulting energy supports both conforming finite-element discretizations with wavenumber-independent error bounds and a plane-wave neural network formulation.","feed_headline":"Penalized energy recovers Helmholtz solution uniquely","feed_subtitle":"Large residual and boundary penalties turn the sign-indefinite wave energy into a coercive one, enabling stable FEM and neural solvers.","key_machinery":"The key identity is the Rellich–Morawetz identity for the Helmholtz operator, which relates an integral of (x·∇u)Lu plus bulk gradient and boundary terms to zero; a second low-order identity with multiplier M0u = −ikβu controls the boundary terms. The paper uses these identities to show that a weighted combination of the physical energy and the least-squares residual has positive coefficients on ∥∇u∥², k²∥u∥², ∥Lu∥², and boundary terms whenever the penalty parameters γ1, γ2 are large enough. The regularized energies themselves carry the argument: Eγ restricts the trial space to functions satisfying the impedance condition, while Fγ adds a boundary residual term so the condition is enforced w","core_discovery":"The paper constructs two regularized energies, Eγ on the space of functions that already satisfy the impedance condition and Fγ on the larger space V where that condition is imposed weakly by a boundary penalty. Theorems 2 and 3 assert that for γ above computable thresholds the quadratic parts of these energies are strongly coercive in the respective norms, and therefore the minimization problems are well posed and the unique minimizers coincide with the weak solution of the Helmholtz equation. The proof uses Rellich and Morawetz identities, which convert the indefinite kinetic-minus-potential term into positive bulk and boundary contributions once the residual penalty is present. The same v","pith_inferences":["Editorial inference: the 'physical energy + residual penalty + boundary penalty' template should transfer to other sign-indefinite time-harmonic boundary problems whenever a Rellich-type multiplier identity is available; the paper does not state this extension.","Editorial inference: because the coercivity and continuity constants are explicit and independent of k, the energy could serve as a certified loss for adaptive sampling or a posteriori error estimation in neural solvers; the paper only reports proof-of-concept experiments.","Editorial inference: the streamed least-squares initialization of the plane-wave network suggests a general warm-start strategy for variational neural solvers, separating the linear high-frequency part from the nonlinear correction; this is not developed in the paper."],"forward_implications":["For γ above explicit thresholds, Eγ and Fγ are strictly convex and each has a unique minimizer, making the Helmholtz problem a well-posed minimization problem.","The coercivity constant is independent of k, removing the k-dependent inf-sup degradation that plagues standard Helmholtz finite element discretizations.","A conforming H² finite element discretization inherits coercivity and satisfies a quasi-optimal error bound with a constant independent of mesh size and wavenumber.","The same energy provides a neural-network loss that is an alternative to residual-based physics-informed formulations for high-wavenumber Helmholtz problems.","The regularization terms vanish on the solution manifold Lu = f, so at the solution the physical energy is recovered and the variational principle keeps its physical meaning."],"fun_headline_variants":["Coercive penalties make Helmholtz energy solvable","Indefinite wave energy fixed by boundary penalties","Regularized energy yields unique Helmholtz solution","Penalty-tuned energies stabilize FEM and neural solvers","From indefinite to coercive: Helmholtz variational principles"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that stationarity of the regularized energy over the impedance-constrained space forces both the interior Helmholtz equation and the impedance boundary condition; the proof tests only compactly supported variations and never evaluates the boundary Euler–Lagrange term, so the boundary part of that premise is assumed rather than established.","fun_headline_variants_meta":{"raw":{"variants":["Coercive penalties make Helmholtz energy solvable","Indefinite wave energy fixed by boundary penalties","Regularized energy yields unique Helmholtz solution","Penalty-tuned energies stabilize FEM and neural solvers","From indefinite to coercive: Helmholtz variational principles"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000147,"raw_usage":{"total_tokens":984,"prompt_tokens":666,"completion_tokens":318,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":410,"completion_tokens_details":{"reasoning_tokens":258}},"tokens_in":410,"tokens_out":318,"duration_ms":25823,"temperature":1.0,"reasoning_tokens":258,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T21:53:53.330947+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the one-dimensional Helmholtz problem on an interval with a known solution, e.g. u = e^{ikx}, and a smooth test function w in V_BC that is not proportional to u. Evaluate the first variation of Eγ at u against w; if the surviving boundary contribution Re(∂n u ̄w)|_0^1 is nonzero, the exact solution is not a stationary point over V_BC, contradicting the theorem. A discrete version: solve the variational problem (35) with conforming H² elements for a manufactured solution and measure whether the boundary impedance residual ∥∂n u_h − iku_h∥_{L²(∂Ω)} tends to zero under mesh refinement.","supporting_citations":[],"review_version":1}