{"id":"2c347525-0028-4e81-8ff6-61c5c4e3fba9","arxiv_id":"2411.15572","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"Two energy-conserving HDG discretizations for the nonlinear Klein-Gordon equation are analyzed, with h^{k+1} error estimates, h^{k+2} superconvergent postprocessing, and second-order time discretization.","lead":"This paper proposes and analyzes two hybridizable discontinuous Galerkin (HDG) schemes for the nonlinear Klein-Gordon equation, with a priori error estimates, discrete energy conservation, and a postprocessing step that raises the convergence order. Generalist readers: the paper offers a recipe for high-order, structure-preserving numerical methods for nonlinear wave equations, though the fully discrete analysis contains proof gaps.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 5.3's O(Δt²) claim requires unstated u_tttt ∈ L²(0,T;L²); without it the bound (5.28) is uncontrolled and the fully discrete error estimate is unproved.","rationale":"The paper's main advertised contribution is the fully discrete analysis: the abstract promises optimal error estimates for the completely discrete method, and Theorem 5.3 states O(h^{k+1}+(Δt)²) with no regularity hypotheses. The proof of precisely this estimate invokes u_tttt in Eq. (5.28). The reader's weakest-assumption analysis identifies the same gap. I checked the surrounding argument: the semidiscrete estimates (Theorems 3.1–3.2) and the energy identities appear consistent with known HDG theory, and the postprocessing theorem follows the standard duality pattern, though it is sketched. The fully discrete error analysis is therefore the least secure part. The false Lemma 5.1 is a real defect, but a corrected quartic bound would still allow the kickback/Grönwall argument because the solution is uniformly bounded in L⁶ by Theorem 5.1 and Lemma 3.1; thus it is a proof gap, not a counterexample to the final estimate. In contrast, the missing u_tttt regularity is a condition on the data: without it, the temporal truncation error is not O((Δt)²) in general, and the theorem as stated is unsupported. The proposed numerical test with u=t^{5/2}sin(πx)sin(πy) directly probes whether the claimed temporal order can hold without this regularity. Since the reader already rejected the paper on essentially these grounds, my stress test does not change the verdict.","tokens_in":36857,"tokens_out":10738,"duration_ms":91779,"concrete_test":"Construct the exact solution u(x,t)=t^{5/2} sin(πx)sin(πy) on Ω=(0,1)², T=1, so that u_tttt ∉ L²(0,T;L²) because ∫_0^1 t^{-3} dt diverges. Run the fully discrete conservative scheme (5.1) with k=2, fixed small h (e.g., h=1/16), and Δt = 1/2, 1/4, 1/8, 1/16; compute the temporal order of ‖u(T)-U^N‖. If the measured order is significantly below 2, the claim in Theorem 5.3 is false as stated. Alternatively, analytically re-derive (5.28) for this u: since the integral of ‖u_tttt‖² is infinite, the presented proof cannot establish (5.25), so the theorem needs an explicit regularity hypothesis or a different argument.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of optimal O(h^{k+1}+(Δt)²) convergence for the fully discrete scheme rests on Theorem 5.3 (Eqs. 5.25–5.26), whose statement imposes no regularity hypotheses on the exact solution. In the proof, the temporal consistency term is expanded in Eq. (5.28) and bounded by C(h^{2(k+1)}‖u_tt‖²_{L∞(H^{k+1})} + (Δt)³∫_{t_{n-1}}^{t_{n+1}}‖u_tttt(s)‖² ds). The last integral requires u_tttt ∈ L²(0,T;L²), an assumption absent from Theorem 5.3 and also from the non-conservative Theorem 5.4, which uses the same expansion. For a solution such as u(x,t)=t^{5/2} sin(πx)sin(πy), the integral ∫_0^1 t^{-3} dt diverges, so the right-hand side of (5.28) is infinite and the proof cannot deliver (5.25). This is not a cosmetic omission: second-order time accuracy for this implicit-midpoint-type scheme is not a generic property of weak solutions; it depends on this unstated smoothness. A secondary gap is Lemma 5.1, whose claimed algebraic inequality is false (taking Φ1=Ψ1=M, Φ2=Ψ2=0 gives left-hand side ~ M⁶ and right-hand side ~ C M⁴), and it is used in the uniqueness and error analysis. That flaw looks repairable by replacing the prefactor with a quartic power and using the uniform L⁶ bounds from Theorem 5.1 and Lemma 3.1, whereas the missing u_tttt regularity is a genuine restriction on the problem data that the theorem must state explicitly.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops and analyzes two hybridizable discontinuous Galerkin (HDG) schemes for the nonlinear Klein-Gordon equation u_tt - Δu + f(u) = 0 with f(u) = u^3 - u. The main claims are: optimal O(h^{k+1}) semidiscrete error for both displacement and flux, O(h^{k+2}) superconvergence for a local postprocessed displacement when k ≥ 1, discrete energy conservation for a fully discrete conservative scheme, fully discrete errors of order O(h^{k+1} + (Δt)^2) for both conservative and non-conservative time discretizations, and a variant HDG method with O(h^{k+2}) displacement error. The analysis uses HDG projection estimates, L^p bounds for the discrete solution to handle the locally Lipschitz nonlinearity, a Baker-type argument for the displacement, and a Brouwer fixed-point argument for discrete well-posedness. Numerical experiments for k = 1, 2, 3 and for several exact solutions are reported and confirm the predicted orders and energy conservation.","tokens_in":37246,"tokens_out":7711,"duration_ms":69341,"significance":"If the results are correct, the paper would be a useful contribution to HDG methods for nonlinear hyperbolic problems: it addresses a locally Lipschitz polynomial nonlinearity without imposing a global Lipschitz condition, proves discrete energy conservation, and provides postprocessing superconvergence together with fully discrete second-order time error estimates. The semidiscrete part is structurally substantial, and the L^p technique used to control the nonlinearity is a genuine strength. However, the fully discrete claims currently rest on an unstated fourth-order time regularity assumption and on a false algebraic inequality in Lemma 5.1, while the superconvergence proof for the variant in Section 6 is only sketched. No fitting of parameters or circularity was detected; the concerns are about technical correctness of key proof steps.","major_comments":[{"comment":"The inequality stated in Lemma 5.1 and displayed in Eq. (5.8) is false as written. Taking Φ1 = Ψ1 = M, Φ2 = Ψ2 = 0 gives a left-hand side behaving like ‖M^3 - M‖^2, which is of order M^6 for large M, whereas the right-hand side is of order M^4 since each difference factor contributes M^2 and the bracketed factor contributes M^2. This lemma is used in the uniqueness proof at Eq. (5.16) and in the error estimate at Eq. (5.30) to reduce the nonlinear term to a difference squared in L^6 times a bounded factor. The lemma is repairable by putting a quartic power in the bracketed factor and then using the uniform L^6 bounds from Theorem 5.1, but as stated the estimate is not valid and the subsequent proofs relying on it lack support.","section":"Section 5.1, Lemma 5.1"},{"comment":"The fully discrete error analysis in Theorem 5.3 needs an unstated regularity assumption on the exact solution. In Eq. (5.28), the Taylor remainder terms R*_1 and R*_2 involve u_tttt, and the displayed bound contains ∫_{t_{n-1}}^{t_{n+1}} ‖u_tttt(s)‖^2 ds. The theorem statement, however, imposes no regularity hypotheses on u beyond the semidiscrete context. The same expansion is used in Theorem 5.4. Without an explicit assumption such as u_tttt ∈ L^2(0,T;L^2(Ω)), the right-hand side of (5.28) is not controlled, so the claimed O((Δt)^2) rates in (5.25), (5.26), (5.39), and (5.40) are not proved. This is not a cosmetic omission: second-order accuracy for this implicit midpoint-type scheme genuinely depends on this fourth-order time smoothness, and the theorem must state it.","section":"Theorem 5.3, Eq. (5.28)"},{"comment":"The proof of the O(h^{k+2}) superconvergence result (6.5) for the variant HDG method is only sketched. After the decomposition into Q1,...,Q8, several terms are bounded, but the text then says that the remaining terms are handled 'systematically' or 'similarly as done in subsection 4.2' without giving the actual estimates. In particular, the nonlinear terms Q4 and Q8 and the time-integrated φ terms require the local Lipschitz structure of f and are not straightforward copies of the earlier argument because the projection and error equations for the variant differ. Since the abstract explicitly advertises the error estimates for this variant, the proof needs to be completed rather than summarized.","section":"Section 6.2, proof of Eq. (6.5)"}],"minor_comments":[{"comment":"In the non-conservative scheme, the flux equation (5.37c) writes the test function as v instead of v_h; this should be corrected.","section":"Eq. (5.37c)"},{"comment":"In the proof of Theorem 5.1, the sentence 'By adding ... on both sides of (5.37c)' appears to refer to Eq. (5.5), not to (5.37c).","section":"Proof of Theorem 5.1"},{"comment":"In the proof of Lemma 3.1, the sentence 'an application of Theorem 2.1 implies' is misleading; the argument used is a Cauchy-Schwarz/Hölder bound, not the energy conservation theorem.","section":"Lemma 3.1"},{"comment":"The statement of Theorem 6.1 says 'u0, u1 ∈ L∞(0,T;H^{k+2})', but u0 and u1 are initial data, not time-dependent functions; the statement should read u0, u1 ∈ H^{k+2}.","section":"Theorem 6.1"},{"comment":"The time-difference notation is overloaded: both ∂_t and overbar-∂_t are introduced in Section 5.1, and the distinction between ∂_tϕ^n and ∂_tϕ^n is easy to miss; a summary table or consistent use of different symbols would improve readability.","section":"Notation, Section 5"}],"recommendation":"major_revision","confidential_remarks":"The paper is not circular and the semidiscrete analysis has substantial value, but the fully discrete part is not currently reliable: the missing u_tttt assumption and the false Lemma 5.1 affect central claims in Theorems 5.3 and 5.4. I nevertheless recommend major revision rather than rejection because both gaps appear repairable within the manuscript's scope: state the needed time regularity explicitly and repair Lemma 5.1 with the quartic factor. The Section 6 variant also needs a genuinely complete proof of (6.5) if it remains an advertised contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The part of this paper worth keeping is the semidiscrete analysis. Extending the Cockburn-Fu-Hungria-Ji-Sanchez-Sayas HDG framework to the locally Lipschitz Klein-Gordon nonlinearity, with the Lp bounds that handle f(u)=u^3-u and a k+2 postprocessing result, is a genuine contribution. The energy-conservation identities in Theorems 2.1 and 5.1 check out, and the a priori estimates for the semidiscrete scheme are coherent. I would cite that part.\n\nThe fully discrete part is where the paper breaks. The stress-test note is right. Lemma 5.1 as stated is false: take Phi1=Psi1=M, Phi2=Psi2=0. The left side grows like M^6 while the right side only grows like M^4. That lemma is used in the uniqueness proof and in the error analysis, so the proof of Theorem 5.2 and Theorem 5.3 is not valid as written. The fix is probably repairable - use a quartic prefactor and the uniform L6 bounds from the energy conservation - but it has to be done.\n\nMore serious is the missing u_tttt regularity in Theorem 5.3. The bound (5.28) contains (Delta t)^3 integral ||u_tttt||^2 ds, and the theorem statement imposes no regularity beyond what was assumed for the semidiscrete case. The counterexample u(x,t)=t^(5/2) sin(pi x) sin(pi y) makes that integral infinite, so the claimed O((Delta t)^2) error bound is not established for the class of solutions the theorem claims to cover. This is not cosmetic. The same expansion is used in Theorem 5.4, so the non-conservative scheme has the same gap. The authors need to state the additional smoothness explicitly and check whether the problem data admit it.\n\nSection 6 is also thin. The superconvergence proof for the second variant is a sketch, with multiple steps deferred to 'similarly' and 'combining all the terms.' The numerical tables are mostly believable but contain some impossible experimental order entries (e.g., jumps in error that cannot produce the reported EOC).\n\nBottom line: this is not a desk reject. The semidiscrete analysis is worth refereeing, and the fully discrete problems are likely repairable. But the current version does not support the abstract's fully discrete claims. Send it to a knowledgeable referee with instructions to focus on Sections 5 and 6, and ask for a revision that fixes Lemma 5.1 and states the missing regularity.","headline":"The semidiscrete HDG analysis and energy conservation are solid and new, but the fully discrete error claims rest on a false lemma and an unstated fourth-order time regularity assumption.","tokens_in":701,"tokens_out":737,"would_cite":false,"duration_ms":18048,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65N12","65N22"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves optimal O(h^{k+1}) HDG convergence for the nonlinear Klein-Gordon equation, O(h^{k+2}) after post-processing, and a fully discrete scheme that conserves energy while keeping error O(h^{k+1} + (Δt)^2).","keywords":["HDG","hybridizable discontinuous Galerkin","Klein-Gordon equation","energy conservation","a priori error estimates","superconvergence","post-processing","nonlinear wave equation"],"falsifier":"Take the fully discrete conservative scheme on the unit square with the manufactured solution u(x,y,t) = $t^{{3.2}}$ sin(πx) sin(πy), whose fourth time derivative is not square-integrable near t = 0; if the observed temporal order of convergence is still 2 in Δt, the u_tttt hypothesis is not necessary, whereas a drop below 2 would show that the missing hypothesis is real. A complementary test is to solve the same problem on a domain with a re-entrant corner, where the spatial solution is only $H^{{k+1}}$; the postprocessed O($h^{{k+2}}$) rate should degrade if it depends on elliptic regularity.","tokens_in":36633,"feed_emoji":"🌊","tokens_out":9658,"duration_ms":77652,"temperature":0.7,"pith_summary":"The paper analyzes two hybridizable discontinuous Galerkin (HDG) schemes for the nonlinear Klein-Gordon equation u_tt − Δu + f(u) = 0 with f(u) = u³ − u, a locally Lipschitz nonlinearity. It proves that the semidiscrete approximation converges at the optimal spatial rate O($h^{{k+1}}$) for both the displacement and the flux when piecewise polynomials of degree k are used, and that a cheap element-by-element post-processing of the displacement improves this to O($h^{{k+2}}$) for k ≥ 1. For the time discretization, a second-order conservative finite-difference scheme is shown to conserve a discrete energy exactly while keeping the overall error at O($h^{{k+1}}$ + (Δt)²); a non-conservative variant that is linear at each time step achieves the same rates. A second HDG variant, with displacement polynomials of degree k+1, gives O($h^{{k+2}}$) for the displacement. Numerical experiments with polynomial degrees 1 to 3 confirm the predicted convergence orders and show energy errors at machine precision.","feed_headline":"Nonlinear wave solver hits optimal rates and conserves energy","feed_subtitle":"A hybridizable Galerkin discretization proves O(h^{k+1}) convergence and exact discrete energy for the Klein-Gordon equation.","key_machinery":"The load-bearing machinery is the HDG spatial discretization together with its projection-based error analysis. The method seeks a triple (u_h, q_h, û_h) of piecewise polynomials of degree k, where q_h approximates the flux ∇u, u_h the displacement, and û_h a single-valued trace on the mesh faces; the numerical flux is defined by q̂_h·n = q_h·n − τ(u_h − û_h). The HDG projection (Π_W, Π_V) is defined locally by L² orthogonality conditions and a trace condition on element faces, and it transfers the regularity of the exact solution into sharp O($h^{{k+1}}$) projection errors. The discrete energy identity supplies the boundedness that tames the locally Lipschitz nonlinear term, and for the fully discrete scheme the symmetric quotient F($U^{{n+1}}$,$U^{{n−1}}$) = (F($U^{{n+1}}$) − F($U^{{n−1}}$))/($U^{{n+1}}$ − $U^{{n−1}}$) with F(s) = (1−s²)²/4 makes the discrete energy exactly conserved.","core_discovery":"The paper establishes that the HDG framework, which rewrites the second-order wave equation as a first-order system in a displacement variable u_h and a flux variable q_h together with a single-valued trace û_h on element faces, can handle the locally Lipschitz nonlinearity f(u) = u³ − u without any global Lipschitz condition. The key step is the discrete energy identity E(t) = 1/2(‖u_ht‖² + ‖q_h‖² + ‖u_h − û_h‖²₊ + 2(F(u_h),1)), with F(s) = (1−s²)²/4, which is conserved exactly and yields the L^p bounds needed to control the nonlinear term. On the strength of this identity and the HDG projection, the paper proves optimal O($h^{{k+1}}$) error estimates for displacement and flux, superconvergence O($h^{{k+2}}$) after local post-processing for k ≥ 1, and a fully discrete error bound O($h^{{k+1}}$ + (Δt)²) for a conservative time-stepping scheme that preserves the discrete energy exactly. The analysis also covers a non-conservative linear-per-step scheme and a mixed-degree variant, and the claims are backed by numerical experiments.","pith_inferences":["A natural extension is to apply the same symmetric-quotient energy-conserving time discretization to other Hamiltonian wave equations with a polynomial potential, such as the sine-Gordon equation; the proof structure would likely carry over as long as the potential gives a priori L^p bounds.","The gap between the stated regularity in the fully discrete theorem and the regularity actually used in the Taylor remainders suggests that the (Δt)² estimate could potentially be proved under weaker time regularity by adapting the paper's own nonstandard energy technique, and a numerical test with rough-in-time data would reveal whether the missing hypothesis is essential.","Because the discrete energy is an exact invariant, the conservative scheme is a candidate for long-time integration of the Klein-Gordon equation; one could test this by measuring phase and amplitude errors over very long intervals."],"forward_implications":["The semidiscrete HDG approximation converges at rate O(h^{k+1}) in L² for both displacement and flux, and the locally post-processed displacement converges at rate O(h^{k+2}) for polynomial degree k ≥ 1.","The fully discrete conservative scheme conserves the discrete energy exactly, E^{n+1/2} = E^{3/2}, and has error O(h^{k+1} + (Δt)²).","The non-conservative variant, which solves only linear systems at each time step, attains the same optimal error rates.","The mixed-degree HDG variant, with displacement polynomials of degree k+1 and flux polynomials of degree k, yields O(h^{k+2}) for the displacement and O(h^{k+1}) for the flux.","Numerical experiments with polynomial degrees 1 through 3 match the predicted convergence orders and give energy errors at machine precision."],"supporting_citations":[{"why":"Supplies the HDG projection and its approximation properties, which are the backbone of the semidiscrete error analysis.","marker":"[11]"},{"why":"Provides the energy-conservative HDG framework for acoustic waves and the Stormer-Numerov time discretization that the fully discrete scheme extends to the nonlinear Klein-Gordon equation.","marker":"[8]"},{"why":"Establishes conditions for HDG superconvergence that underpin the post-processing result.","marker":"[12]"},{"why":"Source of the mixed-degree HDG variant analyzed in Section 6.","marker":"[25]"},{"why":"Gives the earlier conforming finite element analysis with product approximation for the nonlinear Klein-Gordon equation against which the local-Lipschitz treatment is positioned.","marker":"[30]"},{"why":"Supplies the nonstandard energy technique used in Lemma 3.5 to handle the displacement error under reduced regularity.","marker":"[3]"}],"fun_headline_variants":["HDG for nonlinear Klein-Gordon: optimal error, exact energy","Energy-conserving HDG achieves optimal convergence rates","Nonlinear wave HDG: O(h^{k+1}) error, conserved energy","Optimal convergence and exact energy for HDG wave solver","Conservative HDG for Klein-Gordon proves optimal error bounds"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"For the fully discrete error bound to hold, the exact solution must have a square-integrable fourth time derivative, a regularity condition that is used in the Taylor remainder estimates but not stated in the theorem; if the solution is less smooth in time, the (Δt)² part of the claimed convergence is not proven.","fun_headline_variants_meta":{"raw":{"variants":["HDG for nonlinear Klein-Gordon: optimal error, exact energy","Energy-conserving HDG achieves optimal convergence rates","Nonlinear wave HDG: O(h^{k+1}) error, conserved energy","Optimal convergence and exact energy for HDG wave solver","Conservative HDG for Klein-Gordon proves optimal error bounds"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000257,"raw_usage":{"total_tokens":1629,"prompt_tokens":1043,"completion_tokens":586,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":659,"completion_tokens_details":{"reasoning_tokens":498}},"tokens_in":659,"tokens_out":586,"duration_ms":5813,"temperature":1.0,"reasoning_tokens":498,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:10:16.509166+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the fully discrete conservative scheme on the unit square with the manufactured solution u(x,y,t) = $t^{{3.2}}$ sin(πx) sin(πy), whose fourth time derivative is not square-integrable near t = 0; if the observed temporal order of convergence is still 2 in Δt, the u_tttt hypothesis is not necessary, whereas a drop below 2 would show that the missing hypothesis is real. A complementary test is to solve the same problem on a domain with a re-entrant corner, where the spatial solution is only $H^{{k+1}}$; the postprocessed O($h^{{k+2}}$) rate should degrade if it depends on elliptic regularity.","supporting_citations":[{"cited_title":"Cockburn, J","cited_arxiv_id":null,"evidence_quote":"Supplies the HDG projection and its approximation properties, which are the backbone of the semidiscrete error analysis."},{"cited_title":"Cockburn, Z","cited_arxiv_id":null,"evidence_quote":"Provides the energy-conservative HDG framework for acoustic waves and the Stormer-Numerov time discretization that the fully discrete scheme extends to the nonlinear Klein-Gordon equation."},{"cited_title":"Cockburn, W","cited_arxiv_id":null,"evidence_quote":"Establishes conditions for HDG superconvergence that underpin the post-processing result."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Source of the mixed-degree HDG variant analyzed in Section 6."},{"cited_title":"Tourigny","cited_arxiv_id":null,"evidence_quote":"Gives the earlier conforming finite element analysis with product approximation for the nonlinear Klein-Gordon equation against which the local-Lipschitz treatment is positioned."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the nonstandard energy technique used in Lemma 3.5 to handle the displacement error under reduced regularity."}],"review_version":1}