{"id":"ac5789c2-4619-4b09-810f-0fb6f10838b3","arxiv_id":"1908.07403","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"New point-weighting 17-point and 25-point finite difference schemes for the 2D Helmholtz equation with PML are fourth-order consistent and are tuned to reduce numerical dispersion.","lead":"The paper develops two fourth-order finite difference schemes for the 2D Helmholtz equation with perfectly matched layer absorbing boundaries, using a point-weighting strategy. It proves the schemes are consistent with the PML and tunes their free weights to reduce numerical dispersion.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The fourth-order consistency proof assumes PML coefficients are smooth up to sixth order, but the actual damping profile (2)-(3) is only C^1, so the central PML claim is unproven exactly where it matters.","rationale":"The paper makes a constructive and plausible contribution: it proposes a point-weighting 17-point stencil aimed at fixing the PML inconsistency of the earlier derivative-weighting 17-point scheme, and the algebraic construction is internally consistent for smooth coefficients. In the constant-coefficient, no-PML setting, Tables 1-4 give reasonable fourth-order evidence, and the dispersion analysis is a useful addition. The reader's CONDITIONAL verdict is appropriate because the consistency proof is only sketched, with typographical slips, unreported fitted parameters, and no code. The load-bearing technical weakness is the unstated regularity assumption: Proposition 2.1 uses Taylor expansions that require A, B, and C to be smooth up to sixth order, while the paper's own PML profiles (2)-(3) are only C^1 at the interface. This is not an objection to the method's usefulness; it is a precise gap between the theorem as stated and the numerical setting in which the method is promoted. The proposed fix is straightforward: either add a smoothness hypothesis and prove a localized interface-error estimate, or verify fourth-order convergence separately inside the PML. Because such a verification is feasible and the authors' numerical examples do not currently provide it, the correct verdict remains CONDITIONAL rather than a clean accept or a rejection.","tokens_in":30432,"tokens_out":17711,"duration_ms":172128,"concrete_test":"Implement the refined PW 17p scheme in a PML-truncated 2D homogeneous model with a known exact interior solution (e.g., the free-space Green's function or a manufactured plane-wave solution) and the C^1 damping profile (2)-(3) with LPML and a0 as in Example 4.3. Refine the grid from roughly N=101 to N=801, keep the PML thickness fixed in physical units, and compute the C-norm error in the interior excluding a buffer of a few wavelengths next to the PML. Repeat the identical experiment with a C^infty damping profile that has zero derivatives of all orders at the interface, matched to the same absorption strength.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Proposition 2.1 claims unconditional pointwise consistency and fourth-order accuracy for the 25p and 17p schemes with the Helmholtz-PML equation (1). The proof expands p, A, B, C by Taylor's theorem, and the stated remainder coefficients in Eqs. (37) and (39) contain derivatives of A and B up to order five and six, e.g. terms like partial^5_x A and partial^6_x A. The PML damping profiles in Eqs. (2)-(3) are only C^1 at the interface l_x=0 and l_z=0: the second derivative of sigma_x jumps from a nonzero value inside the layer to zero outside, so A = sz/sx and B = sx/sz have no fifth or sixth derivatives there. Consequently, the h^4 truncation-error terms are unbounded near the interface, and the Taylor-based proof does not cover grid points whose 17- or 25-point stencil crosses a PML boundary. The paper does not state a smoothness hypothesis in Proposition 2.1, does not acknowledge that the example profile (2)-(3) lacks the required regularity, and the numerical experiments do not isolate the PML layer: Example 1 uses Dirichlet boundary conditions with no PML, Example 2 uses a single grid size without a convergence study, and Example 3 provides only qualitative images. Thus the central claim that the 17p point-weighting scheme is fourth-order for Helmholtz with this PML is unproven in exactly the setting motivating the paper. If the unsmooth profile degrades the local truncation error to O(h^2) or O(h) at the interface, the practical fourth-order behavior and the claimed replacement of the inconsistent 17p scheme are not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a point-weighting strategy for constructing fourth-order finite difference schemes for the two-dimensional Helmholtz equation with a perfectly matched layer (PML). Two schemes are proposed: an optimal point-weighting 25-point scheme and an optimal point-weighting 17-point scheme. Proposition 2.1 claims that both schemes are pointwise consistent with the Helmholtz-PML equation and are fourth-order accurate under only the sum constraints on the weighting parameters. A dispersion analysis is then used to select the free parameters by minimizing the numerical dispersion, giving the 'refined PW 25p' and 'refined PW 17p' schemes. Numerical experiments compare these schemes with existing rotated 9p, 17p, and 25p methods, and a layered-model example is used to argue that the new 17-point scheme removes the PML inconsistency of the earlier 17-point scheme of Dastour and Liao.","tokens_in":30814,"tokens_out":5953,"duration_ms":60005,"significance":"If the central claim is established, the refined 17-point scheme is a practically valuable contribution: it offers a narrower matrix bandwidth than the 25-point scheme while retaining fourth-order accuracy and PML consistency, which is important for large frequency-domain seismic modeling. The paper provides explicit stencils, dispersion relations, parameter-selection algorithms, and numerical comparisons, and the non-PML convergence experiments in Tables 1-4 show the expected roughly 16-fold error reduction on halving the mesh size. The main gap is that the consistency proof for the PML setting assumes more smoothness than the PML profile specified in the paper actually has, and the numerical examples do not provide a PML-specific convergence study. With that gap addressed, the contribution would justify publication.","major_comments":[{"comment":"The Taylor-expansion proof of Proposition 2.1 assumes that p, A, B, and C are sufficiently smooth for the stated h^4 remainders to be controlled; Eqs. (37) and (39) contain derivatives of A and B up to order five and six (e.g., partial^5_x A and partial^6_x A). However, the PML damping profiles in Eqs. (2)-(3) are only C^1 at the interfaces l_x=0 and l_z=0, since the second derivative of sigma_x and sigma_z jumps from a nonzero value inside the layer to zero outside. Therefore A=sz/sx and B=sx/sz are not C^5 across the interface, and for any grid point whose 17-point or 25-point stencil includes points on both sides of a PML interface, the Taylor remainders in the proof are not bounded. Proposition 2.1 thus does not establish pointwise consistency or fourth-order accuracy in the PML setting that motivates the paper, unless the proof adds a smoothness hypothesis, treats interface-crossing stencils separately, or the PML profile is replaced by a smooth one.","section":"Section 2, Proposition 2.1 and Eqs. (33)-(39)"},{"comment":"The numerical evidence does not yet fill the smoothness gap in the PML case. Example 4.1, which produces the convergence rates in Tables 1-4, uses Dirichlet boundary conditions with no PML, so the observed fourth-order error reduction does not test the PML claim. Example 4.2 uses a single grid size (nx=nz=51) with no refinement study, and Example 4.3 provides only qualitative wavefield snapshots. To support the central claim, the authors should report a convergence study for the Helmholtz-PML problem with the profile (2)-(3), for example by comparing against a known solution on a truncated domain with PML and measuring the C-norm error as h is halved.","section":"Section 4, Tables 1-4 and Examples 4.2-4.3"}],"minor_comments":[{"comment":"The fourth term on the right-hand side of Eq. (6) is evaluated at x_{m-3\\Delta x/2}, but it should be x_{m+3\\Delta x/2}; as written the formula is incorrect.","section":"Section 2, Eq. (6)"},{"comment":"The Taylor expansion for p^*_{m-2,n} in Eq. (32) uses b1 instead of a1; since this expansion is for the 25-point scheme, the coefficient should be a1.","section":"Section 2, Eq. (32) and Proposition 2.1 proof"},{"comment":"The symbol 'g' appears in place of 'gamma' in several places, for example '(g4 + 1)' in Eq. (61) and '(g − 1)(g + 1)' in Eq. (66).","section":"Section 3, Eqs. (61) and (66)"},{"comment":"The 'Result' line of Algorithm 2 says 'a1, . . . ,c4', but the algorithm solves for b1, d2, d3; this should be corrected.","section":"Section 3, Algorithm 2"},{"comment":"The sentence 'From [? ], the analytical solution of this homogeneous model is available' contains a missing citation placeholder; the reference should be supplied.","section":"Section 4.2"},{"comment":"The discussion states that 'zeta and eta depend on k, A and B' where eta was previously defined as 1+1/gamma^2; the intended quantities are likely zeta and xi, the h^4 coefficients in Eqs. (41)-(42).","section":"Section 2, after Eq. (42)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a development of the authors' earlier work [25] and would benefit from a thorough proofreading pass; the number of typographical errors is high. The main substantive issue is the smoothness gap in the PML consistency proof, which the numerical experiments do not currently address. I recommend major revision rather than rejection because the proposed schemes may well be valid and useful, but the central claim needs either a more careful analysis for the non-smooth PML profile or a convincing PML-specific convergence study."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Briefly: the new 17-point point-weighting scheme is a real contribution. The prior 17p scheme from [25] was known to be inconsistent with PML; this paper shows that the point-weighting variant restores consistency, at least formally, and the bandwidth saving vs 25p is real. The interior error tables (Example 1) show clean fourth-order behavior, and the dispersion analysis is standard but careful. I'd give credit for the basic construction and for demonstrating the practical difference between the old and new 17p in the layered example.\n\nThe soft spots are not minor, though. Proposition 2.1 is proved by Taylor expanding A and B and keeping derivatives up to fifth/sixth order in the remainders (Eqs. 37, 39). But the PML profile in (2)-(3) is only C^1 at the interface—the second derivative of sigma jumps. So A and B do not have the derivatives the proof needs, and the Taylor remainder is not under control for stencils crossing the PML boundary. The proposition as stated is therefore unproven in exactly the regime the paper is about. The paper does not flag this, and none of the numerical experiments isolate the PML layer: Example 1 has no PML, Example 2 uses one grid size, Example 3 is qualitative. I would not say the scheme is wrong; but the central fourth-order claim for the Helmholtz-PML equation is not established.\n\nThere are also smaller issues: Eq. (6) has a sign/index typo, Eq. (32) uses b1 where a1 is meant, and Proposition 3.1 states (k_N)^2 = k^2(1 + O(k^6 h^6)) while its own expansion (61) gives a leading k^6 h^4 term—so the stated dispersion order is off by a power of h. The fitted weights used in the experiments are not reported, and no code is provided, so the numerical section is not independently reproducible.\n\nWho is this for? Someone working on frequency-domain finite differences for wave propagation who wants a low-bandwidth fourth-order stencil and is willing to verify the PML behavior themselves. The construction is worth knowing; the analysis is not the last word. I would send it to a serious referee, but with a request to fix the smoothness issue, add a convergence study with PML active, and report the parameters.\n\nI would not cite it as it stands for the PML consistency claim, though I might mention the stencil as a lead to follow.","headline":"A genuinely new 17-point point-weighting stencil for Helmholtz-PML with solid interior convergence, but the PML consistency proof assumes smoothness the paper's own damping profile does not have.","tokens_in":31393,"tokens_out":5136,"would_cite":false,"duration_ms":53402,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65N06","65M06","35J05"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves two point-weighting finite-difference stencils are fourth-order consistent with the Helmholtz equation with PML, fixing the earlier 17-point scheme's inconsistency.","keywords":["Helmholtz equation","perfectly matched layer","point-weighting scheme","fourth-order accuracy","numerical dispersion","17-point scheme","25-point scheme","wavenumber error"],"falsifier":"Run a convergence study for the refined point-weighting 17-point scheme on a manufactured solution with the wavefield supported inside the PML, where the damping profiles (2)-(3) are active; if the C-norm error over the layer does not shrink by roughly a factor of 16 when the grid is halved, the fourth-order claim fails exactly where the PML matters.","tokens_in":30199,"feed_emoji":"🌊","tokens_out":9061,"duration_ms":84300,"temperature":0.7,"pith_summary":"The paper seeks a general way to build fourth-order finite-difference schemes for the two-dimensional Helmholtz equation equipped with perfectly matched layer (PML) absorbing boundaries. Its main constructions are two point-weighting stencils, an optimal 25-point and an optimal 17-point scheme, in which off-center grid values are replaced by weighted averages before the difference operator is applied. The central result is that both stencils are pointwise consistent with the Helmholtz-PML equation and fourth-order accurate, and that the 17-point version repairs an inconsistency in the authors' earlier 17-point scheme while also working when grid spacings differ in the two directions. The authors also derive the numerical-wavenumber error, choose weights by minimizing dispersion over a user-specified range of grid points per wavelength, and demonstrate the benefit in layered-model examples.","feed_headline":"New 17-point stencil fixes PML inconsistency at 4th order","feed_subtitle":"A point-weighting trick keeps Helmholtz solvers fourth-order accurate even where absorbing layers modify the equation.","key_machinery":"The central mechanism is the point-weighting substitution: before applying the fourth-order difference operator, each off-center grid value is replaced by a weighted arithmetic average of neighboring values, with the weights required to sum to one. The 25-point scheme uses the substitutions in (12)-(13); the 17-point scheme uses (23)-(24). These substitutions preserve the $O(h^4)$ Taylor expansion of the operator while coupling the stencil in the transverse direction, which is what restores consistency for the PML equation and permits unequal spacings $\\Delta x \\neq \\Delta z$. The mass term $k^2 C p$ is approximated by the operators $I^{(j)}$ with weights $c_j$ and $d_j$ that also sum to one. Least-squares minimization of the numerical-dispersion functional over propagation angle and grid-points-per-wavelength interval then fixes the free weights.","core_discovery":"On its own terms, the paper's discovery is that a fourth-order finite-difference discretization of the PML-modified Helmholtz equation can be made pointwise consistent by replacing individual grid values in the difference operator with weighted averages of neighboring values, and that this repairs the inconsistency of the authors' earlier 17-point scheme. The two resulting stencils, optimal point-weighting 25p and optimal point-weighting 17p, satisfy the Helmholtz-PML equation up to $O(h^4)$ error terms for essentially arbitrary weights, provided the relevant weights sum to one; the numerical wavenumber obeys $(k_N)^2 = k^2(1+O(k^6h^6))$. The authors present the 17-point version as a practical replacement for the inconsistent 17-point scheme, with a smaller matrix bandwidth and support for unequal spacings $\\Delta x \\neq \\Delta z$.","pith_inferences":["Because the Taylor proof assumes smoothness of the PML coefficients, the observed order may drop to less than four inside the absorbing layer, where the damping profiles in (2)-(3) are only piecewise smooth; the paper does not report a separate convergence test restricted to the layer.","The construction's parameter freedom suggests a natural stress test: fix the weights with the refined strategy on one velocity and frequency range, then evaluate the schemes outside that range, where the dispersion curves imply accuracy may degrade more sharply.","The same point-weighting template should extend to three dimensions, as the authors list as future work, but the wider stencil and corner regions where both damping profiles vary will make the smoothness issue more pronounced."],"forward_implications":["The refined point-weighting 17-point scheme can be used in place of the earlier 17-point scheme for PML problems, keeping a narrower matrix bandwidth than the 25-point scheme while restoring consistency.","Both stencils keep fourth-order accuracy when the step sizes in the two coordinate directions differ, removing the earlier 17-point scheme's restriction to equal spacing.","With weights tuned on a given grid-point-per-wavelength interval, the schemes' normalized phase and group velocity curves stay close to one, so fewer grid points per wavelength are needed than with conventional second-order stencils.","The wavenumber error relation $(k_N)^2 = k^2(1+O(k^6h^6))$ means the pollution error at fixed $kh$ is reduced by two orders relative to second-order schemes."],"supporting_citations":[{"why":"Supplies the prior derivative-weighting 17-point and 25-point schemes whose inconsistency and limitations the new point-weighting schemes are designed to fix.","marker":"[25]"},{"why":"Introduces the point-weighting strategy and flexible weight-selection rule that the new stencils generalize.","marker":"[18]"},{"why":"Provides the PML-consistent 9-point scheme and the refined choice strategy used for selecting the new schemes' parameters.","marker":"[6]"},{"why":"Supplies the rotated 9-point scheme whose dispersion-minimization framework is extended to the new stencils.","marker":"[5]"},{"why":"Gives the extended 25-point operator and the grid-point-per-wavelength interval conventions used in the dispersion analysis.","marker":"[7]"},{"why":"Supplies the PML-modified Helmholtz equation (1) that the new stencils discretize.","marker":"[16,17]"}],"fun_headline_variants":["Fourth-order Helmholtz solver fixes PML inconsistency","Point-weighting stencils achieve 4th order with PML","Optimal 17-point scheme cuts dispersion in PML Helmholtz","New stencils make Helmholtz-PML fourth-order accurate","PML no longer degrades Helmholtz solver to 2nd order"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The fourth-order proof expands $p$ and the PML coefficients $A$, $B$, $C$ by Taylor's theorem at every grid point, so it assumes these fields are smooth everywhere, but the damping profiles used in the numerical examples are only piecewise smooth inside the PML.","fun_headline_variants_meta":{"raw":{"variants":["Fourth-order Helmholtz solver fixes PML inconsistency","Point-weighting stencils achieve 4th order with PML","Optimal 17-point scheme cuts dispersion in PML Helmholtz","New stencils make Helmholtz-PML fourth-order accurate","PML no longer degrades Helmholtz solver to 2nd order"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00024,"raw_usage":{"total_tokens":1507,"prompt_tokens":923,"completion_tokens":584,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":539,"completion_tokens_details":{"reasoning_tokens":496}},"tokens_in":539,"tokens_out":584,"duration_ms":5756,"temperature":1.0,"reasoning_tokens":496,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:19:54.324232+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a convergence study for the refined point-weighting 17-point scheme on a manufactured solution with the wavefield supported inside the PML, where the damping profiles (2)-(3) are active; if the C-norm error over the layer does not shrink by roughly a factor of 16 when the grid is halved, the fourth-order claim fails exactly where the PML matters.","supporting_citations":[{"cited_title":"Dastour, W","cited_arxiv_id":null,"evidence_quote":"Supplies the prior derivative-weighting 17-point and 25-point schemes whose inconsistency and limitations the new point-weighting schemes are designed to fix."},{"cited_title":"Cheng, X","cited_arxiv_id":null,"evidence_quote":"Introduces the point-weighting strategy and flexible weight-selection rule that the new stencils generalize."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the PML-consistent 9-point scheme and the refined choice strategy used for selecting the new schemes' parameters."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the rotated 9-point scheme whose dispersion-minimization framework is extended to the new stencils."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the extended 25-point operator and the grid-point-per-wavelength interval conventions used in the dispersion analysis."}],"review_version":1}