{"id":"2b0048b6-2ac6-478f-bdff-60867c0d49c1","arxiv_id":"1908.09145","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"The standard L1 scheme is shown to be O(tau^(3-alpha)) accurate in the L2 norm for fractional wave equations with nonsmooth data, and a modified version reaches O(tau^2).","lead":"This paper proves error bounds for the L1 scheme, a common method for solving fractional wave equations, when the initial data are nonsmooth. It also introduces a modified L1 scheme with second-order accuracy, which matters for anyone computing fractional wave simulations with singular data.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Modified L1 scheme's O(τ^2) estimate is asserted without proof: Section 3.3 delegates Theorem 3.3 to a 'simple modification', but the required sector and symbol estimates for the β coefficients are never established.","rationale":"The reader's identified weakest assumption, the uniform-bound condition τ^α/h_min^2, is a genuine limitation and is acknowledged by the authors in the introduction, in Remarks 3.6, 3.10 and 4.4, and by Experiment 3. I agree that the constants C_{α,μ0} degenerate as μ0→∞, so the advertised rates are conditional on a quantitatively bounded mesh ratio. However, this is a stated hypothesis rather than a hidden gap, and it does not by itself undermine the conditional theorems. The more load-bearing weakness I see is the unproved status of Theorem 3.3 for the modified L1 scheme. That theorem is the sole basis for the modified scheme's O(τ^2) claim in the abstract and for Theorem 4.3, yet its proof is delegated to a 'simple modification' without any verification that the modified symbol satisfies the sector condition, the lower bound, or the asymptotic cancellation used in Theorems 3.1 and 3.2. The numerical tables provide empirical support, and the claim may well be true, but the central second-order result is not fully established as written. Since the reader's verdict is already CONDITIONAL, this concern does not change the verdict; it strengthens the condition by identifying a specific proof step that should be supplied before acceptance. I therefore recommend UNCHANGED, with the requested verification of the modified-symbol estimates as an explicit condition.","tokens_in":25153,"tokens_out":25981,"duration_ms":277803,"concrete_test":"Independently derive the modified-symbol analogue of Lemma 3.7: for Ψ(z)=e^{-z}(e^z-1)^3\\hatβ(z), determine whether inf_{z∈Υ1\\{0}} |Ψ(z)+μ(1+e^z)|/(μ+|z|^α) > 0 for the same sector angle θ_{α,μ0} as in Lemma 3.3, uniformly over 0<μ≤μ0. Then expand the integrand difference from the proof of Theorem 3.1 using β, and verify the claimed bound |z|^{α+2}+μ|z|^2+μ^2|z|^{2-α} on Υ1. If either verification fails, Theorem 3.3 lacks a proof; if both pass, the gap is closed by inserting the calculation into Section 3.3.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's second headline result, the modified L1 scheme with temporal accuracy O(τ^2), rests on Theorem 3.3, and its PDE counterpart Theorem 4.3. Theorem 3.3 is introduced with the sentence: 'Finally, by a simple modification of the proofs of Theorems 3.1 and 3.2, we readily obtain the following error estimate.' No proof is given. This is not a cosmetic omission: the proofs of Theorems 3.1 and 3.2 depend on several delicate, symbol-specific facts that the modified coefficients must also satisfy. In particular, Lemma 3.3 establishes a sector in which ψ(z)+μ(1+e^z) has no zeros, Lemma 3.7 proves the lower bound |ψ(z)+μ(1+e^z)| > C_{α,μ0}(μ+|z|^α), and the proof of Theorem 3.1 relies on a precise cancellation of the integrand difference, with the leading error term |z|^{α+2}+μ|z|^{3-α}+μ^2|z|^{2-α} producing the O(τ^{3-α}) rate. For the modified symbol β, only the single coefficient β1 is changed, and Remarks 3.8 and 3.9 assert the desired alternate bound with μ|z|^2 in place of μ|z|^{3-α}. But the corresponding modified kernel Ψ(z)=e^{-z}(e^z-1)^3\\hatβ(z) must be shown separately to satisfy the same no-zero sector condition, the same lower bound, and the claimed cancellation. None of these checks appears in the manuscript. The numerical tables are consistent with O(τ^2), but they do not supply the missing proof. If the modified symbol fails any of these estimates, Theorems 3.3 and 4.3, and hence the abstract's modified-scheme claim, would be unsupported.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes the L1 scheme and a modified L1 scheme for a fractional wave equation D_{0+}^{\\alpha-1}(u'-u_1)-\\Delta u=f with 1<\\alpha<2 and nonsmooth initial data. For a full discretization consisting of the L1 rule in time and P1 finite elements in space, the authors claim a new stability estimate and temporal accuracy O(\\tau^{3-\\alpha}) in the L2 norm at positive times under a uniform temporal grid and a uniformly bounded ratio \\tau^\\alpha/h_{\\min}^2. A modified L1 scheme with one altered coefficient is proposed and claimed to achieve temporal accuracy O(\\tau^2). Convergence in the inhomogeneous case is also stated. The proofs for the standard L1 scheme use Laplace transform representations, sector arguments, and contour estimates following the framework of Jin-Lazarov-Zhou. Numerical experiments in Section 5 verify the predicted orders and also demonstrate the deterioration of accuracy for large \\tau^\\alpha/h^2.","tokens_in":25575,"tokens_out":11764,"duration_ms":115676,"significance":"If the main results are valid, the paper fills a genuine gap: it establishes L2-norm temporal accuracy O(\\tau^{3-\\alpha}) for the L1 scheme for fractional wave equations with nonsmooth data, a result not previously available in the literature. The paper also identifies the condition \\tau^\\alpha/h_{\\min}^2\\lesssim 1 as important for the full discretization, and it documents numerically that large values of this ratio degrade accuracy, which is a useful practical warning. The analysis of the standard L1 scheme is detailed and follows the established Laplace-transform route; the stability proofs for both schemes are given, and the numerical experiments are extensive and match the predicted rates, including the predicted suboptimal behavior in the large-ratio regime. The modified-scheme result, however, is a headline claim that is currently not actually proved in the manuscript: Theorem 3.3 is asserted by reference to a 'simple modification' of earlier proofs, with no verification that the modified symbol satisfies the necessary sector, lower-bound, and cancellation estimates.","major_comments":[{"comment":"The modified L1 scheme's O(\\tau^2) error estimate is not proved. The text says only that 'by a simple modification of the proofs of Theorems 3.1 and 3.2, we readily obtain the following error estimate.' The proof of Theorem 3.1 relies on Lemma 3.3 (a zero-free sector for \\psi(z)+\\mu(1+e^z)), Lemma 3.7 (the lower bound |\\psi(z)+\\mu(1+e^z)| \\gtrsim \\mu+|z|^\\alpha), and a delicate cancellation estimate bounding the integrand difference by |z|^{\\alpha+2}+\\mu|z|^{3-\\alpha}+\\mu^2|z|^{2-\\alpha}. For the modified symbol \\hat\\beta, the corresponding kernel \\Psi(z)=e^{-z}(e^z-1)^3\\hat\\beta(z) must be shown to satisfy the same sector condition, the same lower bound, and the improved cancellation error \\mu|z|^2+|z|^{2+\\alpha}. Remarks 3.8 and 3.9 only state the desired bounds; no proof is supplied for any of these three requirements. The numerical tables in Section 5 are consistent with O(\\tau^2), but they do not replace the missing proof. Because Theorem 4.3 and the abstract's modified-scheme claim rest on Theorem 3.3, the full proof must be provided, or the claim must be explicitly downgraded to a conjecture supported by numerics.","section":"Section 3.3, Theorem 3.3"},{"comment":"The passage from the scalar ODE theorems to the full discretizations is compressed into the sentence 'By the above procedure, we have the following two theorems.' For the standard L1 scheme this is probably routine but should still be written out at least as a lemma: one needs to apply Theorems 3.1 and 3.2 mode by mode with \\lambda=\\lambda_i, to verify that the constants are uniform under the standing assumption \\tau^\\alpha/h_{\\min}^2\\le \\mu_0, and to justify summing the f-dependent terms in Theorem 3.2 over the eigenbasis of -\\Delta_h. For the modified scheme the procedure is not routine at all, because Theorem 3.3 itself is unproved. The paper should either include the detailed eigenvalue decomposition argument or state clearly which part of the proof is being deferred.","section":"Section 4, Theorems 4.2 and 4.3"},{"comment":"The proof of the semidiscrete spatial estimate (56) uses the resolvent bound \\|(z^\\alpha-\\Delta)^{-1}-(z^\\alpha-\\Delta_h)^{-1}P_h\\|_{\\mathcal{L}(L^2)} \\lesssim h^2 for all z\\in\\Upsilon\\setminus\\{0\\}, citing [21, Theorem 2.1]. The uniformity in z on the whole contour \\Upsilon is essential for the subsequent integration (\\int e^{sz} h^2 dz \\sim s^{-1}h^2). Please state the exact form of the resolvent bound that is being quoted and explain why it holds uniformly along the contour, including near the origin. If the known bound has an additional factor such as |z|^{-\\alpha}, the integrals in the proof of Lemma 4.1 would need to be re-examined.","section":"Section 4.2, Lemma 4.1"}],"minor_comments":[{"comment":"The abstract and introduction contain grammatical slips, e.g. 'The convergence of these schemes in inhomogeneous case are also established' and 'As a extension of integer order equation'; these should be corrected in a final pass.","section":"Abstract and Introduction"},{"comment":"After defining \\mathcal{E}(t):=\\int_0^t(E-\\tilde E)(s)ds, the proof of Lemma 3.9 immediately writes E(t) for the same quantity, which conflicts with the kernel E(t) defined in (30). Please use a distinct symbol for the integrated error consistently throughout the proof.","section":"Section 3.2.3"},{"comment":"In the proof of Lemma 3.7 the expansion of (1+re^{i\\theta})^{-1}\\psi(re^{i\\theta}) is written with a factor 1/2, namely r^\\alpha e^{i\\alpha\\theta}/2 + r^{\\alpha+1}g(r), whereas (19) states \\psi(re^{i\\theta})/(r^\\alpha(\\cos\\alpha\\theta+i\\sin\\alpha\\theta))\\to 1, which suggests the leading term should be r^\\alpha e^{i\\alpha\\theta} without the factor 1/2. Please check and correct this, or explain the convention used for the expansion.","section":"Section 3.2.2, Lemma 3.7"},{"comment":"In Table 6, for problem (e) with \\alpha=1.2 and h=2^{-8}, the reported Error2 value 3.33e-5 appears inconsistent with the previous row (1.29e-5 at h=2^{-7}) and the stated order 1.96; it should probably be 3.33e-6. Please verify the entry.","section":"Section 5.2, Table 6"},{"comment":"The theorems involving f(0) assume f'\\in L^1; in that case f is absolutely continuous and f(0) is well defined as a trace, but this should be stated explicitly for the reader, especially in the PDE setting where f takes values in L^2(\\Omega).","section":"Section 3 and Section 4"}],"recommendation":"major_revision","confidential_remarks":"The paper contains a genuinely valuable result for the standard L1 scheme, and the main proofs for that part appear sound and are accompanied by thorough numerics. The issue is the modified L1 scheme: Theorem 3.3 is asserted without proof, and the proof is not a routine one-line consequence of the previous arguments because the modified symbol must satisfy several nontrivial sector and symbol estimates. If the authors can supply a complete proof of Theorem 3.3 (or of a slightly weaker but precisely stated version), the paper could be suitable for publication. Otherwise, the authors should consider restructuring the paper around the O(\\tau^{3-\\alpha}) result and presenting the modified scheme's O(\\tau^2) behavior as a numerical observation, with the claim clearly marked as unproved."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague—\n\nHere’s the short version: the paper has one solid new result and one unsupported headline. The solid part is the L2-norm analysis of the standard L1 scheme for fractional wave equations with nonsmooth data. That was genuinely open, and the authors prove O(τ^{3−α}) at positive times through a careful Laplace-transform argument, in addition to a new stability estimate and a clear statement of the τ^α/h^2 condition. The numerical experiments follow the predicted rates and—more tellingly—reproduce the predicted deterioration when τ^α/h^2 is large, which gives the analysis real credibility.\n\nThe soft spot is exactly where the stress-test note points. The modified L1 scheme, advertised as achieving O(τ^2), is not actually proved. Theorem 3.3 is introduced as “a simple modification of the proofs of Theorems 3.1 and 3.2,” and then nothing. The modified symbol β only changes β1 and gets the derivative condition (β̂−z^{α−3})(0)=0, but the hard work in the standard proof—Lemma 3.3’s no-zero sector, Lemma 3.7’s lower bound on |ψ+μ(1+e^z)|, and the exact cancellation that yields the rate—is symbol-specific. None of those checks is shown for Ψ(z)=e^{-z}(e^z−1)^3β̂(z). Maybe the checks are routine; maybe they aren’t. As written, Theorem 3.3 and its PDE analogue Theorem 4.3 are assertions, not proven statements. The numerical tables are consistent with O(τ^2), but they don’t fill the gap.\n\nOther concerns are minor by comparison. No code or data is shipped, so the experiments are hard to reproduce exactly, though the tables are detailed. The constants C_{α,μ0} blow up as μ0→∞; the authors are upfront about it and Experiment 3 confirms the practical effect. The self-reference to their earlier Petrov–Galerkin paper [22] is appropriate, since they are extending exactly that result.\n\nWho gets value: anyone working on L1-type discretizations for fractional wave problems. The standard-scheme contribution alone makes the paper worth refereeing. But I would not accept as is. The referee should demand a complete proof of Theorem 3.3, or a clearly labeled derivation of the β symbol estimates; if those checks don’t go through, the modified-scheme claim should be downgraded to numerical observation. Recommend: send to peer review, with the missing proof as a blocking issue.","headline":"Genuine new result for the standard L1 scheme; the modified-scheme O(τ^2) claim is asserted, not proved.","tokens_in":26053,"tokens_out":2903,"would_cite":true,"duration_ms":31000,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M60","65M12","35R11","65M15"],"pacs":[],"model":"deepseek-v4-flash","headline":"For nonsmooth data, the L1 scheme reaches O(tau^{3-alpha}) in L2 at positive times, if tau^alpha/h^2 stays bounded.","keywords":["L1 scheme","fractional wave equation","nonsmooth data","finite element method","discrete Laplace transform","convergence analysis","modified L1 scheme"],"falsifier":"Run the standard L1 scheme on the nonsmooth-data test u0(x)=$x^{{-0.49}}$ in (0,1), keep tau fixed, and halve h repeatedly so tau^$\\alpha$/$h^{2}$ leaves the bounded regime; the observed temporal error at t=1 should stop decaying like $tau^{{3-alpha}}$ and should worsen as h shrinks. Alternatively, repeat with a nonuniform time grid: a clean O($tau^{{3-alpha}}$) rate there would contradict the paper's stated restriction, since the analysis does not cover that case.","tokens_in":24959,"feed_emoji":"🧮","tokens_out":6226,"duration_ms":60935,"temperature":0.7,"pith_summary":"This paper pins down what the well-known L1 time-stepping scheme actually delivers for fractional wave equations when the initial data are only square-integrable, not smooth. The central result is that a full discretization combining the L1 scheme in time with piecewise-linear finite elements in space has L2 error at positive times of order O($tau^{{3-alpha}}$), and a slightly modified L1 scheme attains the rate O($tau^{2}$). Both rates hold with nonsmooth data, and the inhomogeneous case with nonzero forcing is covered as well. The practical price is a grid condition: the time step tau and the smallest spatial mesh diameter h_min must satisfy tau^$\\alpha$/$h_min^{2}$ uniformly bounded. If that condition fails, the paper's analysis says temporal accuracy deteriorates, and its Experiment 3 shows exactly that.","feed_headline":"L1 scheme reaches O(tau^{3-alpha}) if tau^alpha/h^2 stays bounded","feed_subtitle":"Standard and modified L1 time-stepping now have proven convergence rates with nonsmooth initial data.","key_machinery":"The load-bearing object is the discrete Laplace transform of the time-stepping kernel. For the standard L1 scheme the kernel is b_j = $j^{{2-alpha}}$/Gamma(3-$\\alpha$), whose transform satisfies \\hat b(z) = \\sum_{k=-\\infty}^\\infty (z + 2k\\pi i)^{\\$\\alpha$-3} in a strip, and the companion symbol \\psi(z) = $e^{{-z}}$(e^z-1)^3 \\hat b(z) approximates z^\\$\\alpha$. The analysis compares \\psi(z)+\\mu(e^z+1) with the continuum symbol z^\\$\\alpha$ + 2\\mu, where \\mu = \\$\\lambda$ \\tau^\\$\\alpha$/2, on a pair of rays avoiding the eigenvalues. The standard kernel leaves a nonzero constant in \\hat b(z) - $z^{{\\alpha-3}}$ at z=0, which is what costs the 3-$\\alpha$ rate. The modified scheme changes only beta_1, adding a constant that forces (\\hat\\$\\beta$(z)-$z^{{\\alpha-3}}$)(0)=0 and removes that leading error, yielding second order.","core_discovery":"The paper's target is the L2($\\Omega$)-norm convergence of two fully discrete time-stepping methods for the fractional wave equation. For the standard L1 scheme, Theorem 4.2 gives, at each positive time t_k, an error bounded by ($t_k^{{alpha-3}}$ $tau^{{3-alpha}}$ + $t_k^{{-alpha}}$ $h^{2}$) times the L2 norm of u0, plus analogous contributions from u1 and the forcing term; in particular, the temporal order is 3-$\\alpha$. For the modified L1 scheme, Theorem 4.3 gives temporal order 2. The proofs write both the exact and numerical solutions as contour integrals over a pair of rays, using the discrete Laplace transform of the time-stepping kernel, and then compare the integrands. The modified kernel is constructed so that the leading symbol error at z=0 vanishes, which upgrades the temporal rate from 3-$\\alpha$ to 2.","pith_inferences":["Beyond the paper: the symbol-correction trick that fixes beta_1 could be tried on nonuniform time grids by replacing the global correction with a local one, though the paper explicitly says its techniques do not apply there.","Beyond the paper: the tau^alpha/h^2 condition behaves like a CFL-type restriction, so one could test whether local time stepping with tau roughly proportional to h^{2/alpha} preserves the proved rates at lower computational cost.","Beyond the paper: an analogous modified L1 scheme might yield O(tau^2) for variable-coefficient or nonlinear fractional wave equations, but the present analysis gives no evidence beyond the linear constant-coefficient setting.","Beyond the paper: because the limiting step is a one-dimensional symbol comparison, the same Laplace-transform machinery could be reused to derive sharp constants or to predict when the 3-alpha rate can be improved for smoother data."],"forward_implications":["Practitioners can use the standard L1 scheme for fractional wave equations with nonsmooth initial values and still expect temporal convergence in L2 of order O(tau^{3-alpha}), rather than only first order.","The modified L1 scheme gives a genuine O(tau^2) temporal rate with nonsmooth data, which is useful when long-time accuracy is limited by the initial-data singularity.","For inhomogeneous problems, the error is controlled by f(0) and the L1 norm of f', with the same temporal exponents and a logarithmic factor at alpha = 3/2.","Spatial discretization with piecewise-linear elements contributes the expected O(h^2) term up to a logarithmic factor, so temporal and spatial errors balance naturally when the two contributions are comparable.","A user must monitor tau^alpha/h_min^2; refining space alone while keeping tau fixed can destroy the temporal accuracy rather than improve the total error."],"supporting_citations":[{"why":"Supplies the contour-integral representation of the exact solution and the semidiscrete spatial error estimate used in Lemma 4.1.","marker":"[12]"},{"why":"Introduces the L1 discretization scheme for the diffusion-wave system on which Discretization 1 is built.","marker":"[33]"},{"why":"Identifies the low-order Petrov-Galerkin method with the L1 scheme and provides the H^1-norm convergence result this paper extends to the L2 norm.","marker":"[22]"},{"why":"Gives the Hankel integral representation and series form of the discrete Laplace transform \\hat b(z) used throughout the analysis.","marker":"[35]"},{"why":"Provides the modified-L1 idea and its O(tau^{2-alpha}) analysis for subdiffusion, which the second scheme adapts to the fractional wave case.","marker":"[36]"},{"why":"Supplies the nonsmooth-data convolution-quadrature framework and resolvent bounds used to control the spatial semidiscretization.","marker":"[21]"},{"why":"Establishes the L1 scheme's O(tau) rate for subdiffusion, serving as the baseline that the new 3-alpha rate improves upon.","marker":"[11]"}],"fun_headline_variants":["Fractional wave L1 scheme: proven 3-alpha order with nonsmooth data","Fractional wave L1: order 3-alpha proven, modified reaches order 2","L1 scheme for fractional waves: 3-alpha accuracy, modified gives 2 for rough data","Fractional wave schemes: L1 gets 3-alpha, modified gets 2 for nonsmooth data"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claimed rates collapse if the time grid is not uniform or if tau^$\\alpha$/$h_min^{2}$ is not uniformly bounded, since the proof constants blow up as that ratio grows.","fun_headline_variants_meta":{"raw":{"variants":["Fractional wave L1 scheme: proven 3-alpha order with nonsmooth data","Fractional wave L1: order 3-alpha proven, modified reaches order 2","L1 scheme for fractional waves: 3-alpha accuracy, modified gives 2 for rough data","Fractional wave schemes: L1 gets 3-alpha, modified gets 2 for nonsmooth data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000828,"raw_usage":{"total_tokens":3558,"prompt_tokens":822,"completion_tokens":2736,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":438,"completion_tokens_details":{"reasoning_tokens":2653}},"tokens_in":438,"tokens_out":2736,"duration_ms":18730,"temperature":1.0,"reasoning_tokens":2653,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:20:15.372012+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the standard L1 scheme on the nonsmooth-data test u0(x)=$x^{{-0.49}}$ in (0,1), keep tau fixed, and halve h repeatedly so tau^$\\alpha$/$h^{2}$ leaves the bounded regime; the observed temporal error at t=1 should stop decaying like $tau^{{3-alpha}}$ and should worsen as h shrinks. Alternatively, repeat with a nonuniform time grid: a clean O($tau^{{3-alpha}}$) rate there would contradict the paper's stated restriction, since the analysis does not cover that case.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the contour-integral representation of the exact solution and the semidiscrete spatial error estimate used in Lemma 4.1."},{"cited_title":"Sun and X","cited_arxiv_id":null,"evidence_quote":"Introduces the L1 discretization scheme for the diffusion-wave system on which Discretization 1 is built."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identifies the low-order Petrov-Galerkin method with the L1 scheme and provides the H^1-norm convergence result this paper extends to the L2 norm."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Hankel integral representation and series form of the discrete Laplace transform \\hat b(z) used throughout the analysis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the modified-L1 idea and its O(tau^{2-alpha}) analysis for subdiffusion, which the second scheme adapts to the fractional wave case."},{"cited_title":"Lubich, I","cited_arxiv_id":null,"evidence_quote":"Supplies the nonsmooth-data convolution-quadrature framework and resolvent bounds used to control the spatial semidiscretization."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the L1 scheme's O(tau) rate for subdiffusion, serving as the baseline that the new 3-alpha rate improves upon."}],"review_version":1}