pith. sign in

arxiv: 2503.13126 · v2 · pith:MITCRL6Lnew · submitted 2025-03-17 · 🧮 math.NA · cs.NA· math.AP

Error analysis of the Strang splitting for the 3D semilinear wave equation with finite-energy data

Pith reviewed 2026-05-22 23:41 UTC · model grok-4.3

classification 🧮 math.NA cs.NAmath.AP
keywords semilinear wave equationStrang splittingerror analysisfinite-energy dataStrichartz estimatesconvergence rates3D torusnumerical methods for PDEs
0
0 comments X

The pith

Strang splitting achieves almost second-order L2 convergence for the 3D cubic semilinear wave equation with finite-energy data.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper studies error bounds for a Strang splitting time integrator applied to the semilinear wave equation on the three-dimensional torus when initial data have only finite energy. For cubic nonlinearities it proves an error that is almost of order two in the L2 norm and almost of order one in the H1 norm; for quartic nonlinearities the orders drop by one half. The proofs combine continuous and discrete Strichartz estimates with integration and summation by parts to obtain cancellations in the local error terms. Bounds are also derived for the fully discrete scheme that adds a Fourier pseudo-spectral spatial discretization, and a numerical test is presented to indicate sharpness of the rates.

Core claim

For the Strang splitting applied to the 3D semilinear wave equation with finite-energy data on the torus, the global error is almost O(τ²) in L² and almost O(τ) in H¹ when the nonlinearity is cubic, and almost O(τ^{3/2}) in L² and almost O(τ^{1/2}) in H¹ when the nonlinearity is quartic; these rates are obtained from continuous-time and discrete-time Strichartz estimates together with summation-by-parts cancellations and extend to the fully discrete Fourier pseudo-spectral version.

What carries the argument

The Strang splitting time-stepping operator, whose local truncation error is controlled by discrete Strichartz estimates and summation-by-parts identities that cancel leading error contributions.

If this is right

  • The same error bounds apply after replacing the continuous spatial operator by a Fourier pseudo-spectral discretization.
  • The analysis supplies explicit constants that depend only on the finite-energy norm of the data.
  • The numerical example confirms that the predicted rates cannot be improved in general under the finite-energy assumption.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If discrete Strichartz estimates can be verified for other time-splitting schemes, the same cancellation technique may yield convergence results for related dispersive equations.
  • The finite-energy setting makes the results directly relevant to physical models where higher Sobolev regularity is unavailable.
  • The approach could be tested on non-periodic domains once appropriate discrete Strichartz estimates are established there.

Load-bearing premise

Discrete-time Strichartz estimates are assumed to hold for the Strang splitting operator on the torus with finite-energy data.

What would settle it

A sequence of numerical runs on a fixed finite-energy solution in which the observed L2 error fails to improve by nearly a factor of four when the time step is halved would contradict the claimed almost second-order rate for the cubic case.

Figures

Figures reproduced from arXiv: 2503.13126 by Maximilian Ruff.

Figure 1
Figure 1. Figure 1: Errors for α = 3 [PITH_FULL_IMAGE:figures/full_fig_p036_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Errors for α = 4 We also did the experiment for the critical power α = 5, see [PITH_FULL_IMAGE:figures/full_fig_p037_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Errors for α = 5 References [1] L. Alexandrov and V. Popov. “Onesided Trigonometrical approximation of Periodic Multivariate Functions”. In: Math. Balkanica, New Series 2. Fasc. 2–3 (1988), pp. 230– 243. [2] H. Bahouri, J. Chemin, and R. Danchin. Fourier Analysis and Nonlinear Partial Differential Equations. Springer, Berlin Heidelberg, 2011 [PITH_FULL_IMAGE:figures/full_fig_p037_3.png] view at source ↗
read the original abstract

We study a variant of the Strang splitting for the time integration of the semilinear wave equation under the finite-energy condition on the torus $\mathbb{T}^3$. In the case of a cubic nonlinearity, we show almost second-order convergence in $L^2$ and almost first-order convergence in $H^1$. If the nonlinearity has a quartic form instead, we show an analogous convergence results, where the order is reduced by $1/2$ in both cases. To our knowledge these are the best convergence results available for the 3D cubic and quartic wave equations under the finite-energy condition. Our approach relies on continuous- and discrete-time Strichartz estimates. We also make use of the integration and summation by parts formulas to exploit cancellations in the error terms. Moreover, error bounds for a full discretization using the Fourier pseudo-spectral method in space are given. Finally, we discuss a numerical example indicating the sharpness of our theoretical results.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper analyzes error bounds for a variant of the Strang splitting scheme applied to the 3D semilinear wave equation on the torus T^3 with finite-energy (H^1 × L^2) initial data. For cubic nonlinearity it claims almost second-order convergence in L^2 and almost first-order convergence in H^1; for quartic nonlinearity the orders drop by 1/2. The proofs rely on continuous- and discrete-time Strichartz estimates together with integration/summation by parts to cancel error terms. Bounds are also given for the fully discrete Fourier pseudo-spectral version, and a numerical example is presented to indicate sharpness.

Significance. If the discrete-time Strichartz estimates for the Strang splitting operator close at the stated regularity, the results would constitute the best available convergence theory for these low-regularity 3D wave problems. The combination of Strichartz estimates with summation-by-parts cancellation is a technically interesting approach that could extend to other splitting methods.

major comments (2)
  1. [Abstract and the section establishing discrete Strichartz estimates] The central convergence claims rest on discrete-time Strichartz estimates for the Strang splitting flow itself at H^1 × L^2 regularity (abstract). The manuscript must supply a self-contained statement and proof of these estimates (including the precise dependence on the time step and any implicit smoothing), because standard continuous Strichartz estimates do not automatically transfer to the numerical flow.
  2. [Error analysis and summation-by-parts argument] After summation by parts, the error accumulation must be tracked explicitly to confirm that the “almost” orders are obtained without hidden losses that depend on the finite-energy norm; the current abstract description leaves the constant tracking and the precise cancellation mechanism unclear.
minor comments (2)
  1. [Introduction and preliminaries] Clarify the precise variant of Strang splitting used (e.g., the ordering of the linear and nonlinear steps) and state all assumptions on the nonlinearity (growth, smoothness) in a single preliminary section.
  2. [Numerical results] The numerical example should include a table or plot comparing observed rates against the theoretical almost-orders for both cubic and quartic cases.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major point below and will incorporate revisions to strengthen the presentation of the discrete Strichartz estimates and the error analysis.

read point-by-point responses
  1. Referee: [Abstract and the section establishing discrete Strichartz estimates] The central convergence claims rest on discrete-time Strichartz estimates for the Strang splitting flow itself at H^1 × L^2 regularity (abstract). The manuscript must supply a self-contained statement and proof of these estimates (including the precise dependence on the time step and any implicit smoothing), because standard continuous Strichartz estimates do not automatically transfer to the numerical flow.

    Authors: We agree that a fully self-contained statement and proof of the discrete-time Strichartz estimates for the Strang splitting flow at H^1 × L^2 regularity is essential. The manuscript currently derives these estimates from the structure of the splitting operator and references the continuous-time versions, but we will expand the relevant section to include an explicit statement of the discrete estimates, a complete proof, the precise dependence on the time step, and any smoothing effects. This will ensure the argument does not rely on automatic transfer from the continuous case. revision: yes

  2. Referee: [Error analysis and summation-by-parts argument] After summation by parts, the error accumulation must be tracked explicitly to confirm that the “almost” orders are obtained without hidden losses that depend on the finite-energy norm; the current abstract description leaves the constant tracking and the precise cancellation mechanism unclear.

    Authors: We will revise the error analysis to track the accumulation of terms after summation by parts in full detail. The revised version will explicitly bound all constants in terms of the finite-energy norm, demonstrate the absence of hidden losses, and clarify the precise cancellation mechanism that yields the almost second-order (cubic) and reduced (quartic) rates in both L^2 and H^1. The abstract will be updated to reflect these clarifications. revision: yes

Circularity Check

0 steps flagged

No circularity: convergence rates derived from external Strichartz estimates plus summation by parts

full rationale

The paper establishes almost second-order L2 and first-order H1 convergence (cubic case) or reduced orders (quartic) for Strang splitting on the 3D semilinear wave equation with finite-energy data. The derivation explicitly invokes continuous- and discrete-time Strichartz estimates together with integration/summation by parts to control error terms after summation. These estimates are drawn from the literature rather than fitted to the target result or justified solely by self-citation. No equation reduces the claimed rates to a parameter fit, no ansatz is smuggled via prior work by the same author, and the central claims retain independent content once the external estimates are granted. The analysis is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only review; the central claims rest on the validity of continuous and discrete Strichartz estimates for the wave propagator on T^3 and on the ability to exploit cancellations via integration/summation by parts.

axioms (1)
  • domain assumption Continuous- and discrete-time Strichartz estimates hold for the linear wave equation on the 3-torus with finite-energy data.
    Invoked to bound the error terms in the splitting analysis.

pith-pipeline@v0.9.0 · 5697 in / 1306 out tokens · 39943 ms · 2026-05-22T23:41:25.366678+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Improved error estimates for low-regularity integrators using space-time bounds

    math.NA 2025-03 unverdicted novelty 7.0

    Proves optimal order-1 and order-2 convergence for an exponential integrator on NLS and corrected Lie splitting on NLW by applying Strichartz and null-form estimates in continuous time.

Reference graph

Works this paper leans on

30 extracted references · 30 canonical work pages · cited by 1 Pith paper

  1. [1]

    Onesided Trigonometrical approximation of Periodic Multivariate Functions

    L. Alexandrov and V. Popov. “Onesided Trigonometrical approximation of Periodic Multivariate Functions”. In:Math. Balkanica, New Series2. Fasc. 2–3 (1988), pp. 230– 243

  2. [2]

    Bahouri, J

    H. Bahouri, J. Chemin, and R. Danchin.Fourier Analysis and Nonlinear Partial Differential Equations. Springer, Berlin Heidelberg, 2011. 38 REFERENCES

  3. [3]

    On averaged exponential integrators for semilinear Klein–Gordon equations with solutions of low-regularity

    S. Buchholz, B. Dörich, and M. Hochbruck. “On averaged exponential integrators for semilinear Klein–Gordon equations with solutions of low-regularity”. In:SN Partial Differ. Equ. Appl.2 (2021), pp. 2662–2963

  4. [4]

    Closing the gap between trigonometric integrators and splitting methods for highly oscillatory differential equations

    S. Buchholz, L. Gauckler, V. Grimm, M. Hochbruck, and T. Jahnke. “Closing the gap between trigonometric integrators and splitting methods for highly oscillatory differential equations”. In:IMA J. Numer. Anal.38.1 (2018), pp. 57–74

  5. [5]

    Random data Cauchy theory for supercritical wave equations I: local theory

    N. Burq and N. Tzvetkov. “Random data Cauchy theory for supercritical wave equations I: local theory”. In:Invent. math.173 (2008), pp. 449–475

  6. [6]

    Numerical approximation of discontinuous solutions of the semilinear wave equation

    J. Cao, B. Li, J. Lin, and F. Yao. “Numerical approximation of discontinuous solutions of the semilinear wave equation”. In:SIAM J. Numer. Anal.63.1 (2025), pp. 214–238

  7. [7]

    On the splitting method for the nonlinear Schrödinger equation with initial data inH 1

    W. Choi and Y. Koh. “On the splitting method for the nonlinear Schrödinger equation with initial data inH 1”. In:Discrete Contin. Dyn. Syst.41 (2021), pp. 3837–3867

  8. [8]

    Long-time-step methods for oscilla- tory differential equations

    B. Garcia-Archilla, J. Sanz-Serna, and R. Skeel. “Long-time-step methods for oscilla- tory differential equations”. In:SIAM J. Sci. Comput.20 (1999), pp. 930–963

  9. [9]

    Error analysis of trigonometric integrators for semilinear Klein-Gordon equations

    L. Gauckler. “Error analysis of trigonometric integrators for semilinear Klein-Gordon equations”. In:SIAM J. Numer. Anal.53 (2015), pp. 1082–1106

  10. [10]

    Grafakos

    L. Grafakos. Classical Fourier Analysis. Springer, New York, 2014

  11. [11]

    Error analysis of exponential integrators for oscillatory second-order differential equations

    V. Grimm and M. Hochbruck. “Error analysis of exponential integrators for oscillatory second-order differential equations”. In:J. Phys. A.39 (2006), pp. 5495–5507

  12. [12]

    Guo.Spectral Methods and Their Applications

    B. Guo.Spectral Methods and Their Applications. World Scientific, Singapore, 1998

  13. [13]

    Hairer, C

    E. Hairer, C. Lubich, and G. Wanner.Geometric Numerical Integration: Structure- Preserving Algorithms for Ordinary Differential Equations. Springer, Berlin Heidel- berg, 2006

  14. [14]

    Best onesided approximation and mean approximation by interpolation polynomials of periodic functions

    V. Hristov. “Best onesided approximation and mean approximation by interpolation polynomials of periodic functions”. In:Math. Balkanica, New Series3. Fasc. 3–4 (1989), pp. 418–429

  15. [15]

    A splitting method for the nonlinear Schrödinger equation

    L. Ignat. “A splitting method for the nonlinear Schrödinger equation”. In:J. Differ- ential Equations250 (2011), pp. 3022–3046

  16. [16]

    Numerical dispersive schemes for the nonlinear Schrödinger equation

    L. Ignat and E. Zuazua. “Numerical dispersive schemes for the nonlinear Schrödinger equation”. In:SIAM J. Numer. Anal.47.2 (2009), pp. 1366–1390

  17. [17]

    Low Regularity Full Error Estimates for the Cubic Nonlinear Schrödinger Equation

    L. Ji, A. Ostermann, F. Rousset, and K. Schratz. “Low Regularity Full Error Estimates for the Cubic Nonlinear Schrödinger Equation”. In:SIAM J. Numer. Anal. 62.5 (2024), pp. 2071–2086

  18. [18]

    Global Solutions to Maxwell Equations in a Ferromagnetic Medium

    J. Joly, G. Métivier, and J. Rauch. “Global Solutions to Maxwell Equations in a Ferromagnetic Medium”. In:Ann. Henri Poincaré1 (2000), pp. 307–340

  19. [19]

    The Cauchy problem for a semilinear wave equation I

    L. Kapitanskii. “The Cauchy problem for a semilinear wave equation I”. In:J. Sov. Math. 49 (1990), pp. 1166–1186

  20. [20]

    A second-order low-regularity correction of Lie splitting for the semilinear Klein–Gordon equation

    B. Li, K. Schratz, and F. Zivcovich. “A second-order low-regularity correction of Lie splitting for the semilinear Klein–Gordon equation”. In:ESAIM Math. Model. Numer. Anal.57.2 (2023), pp. 899–919

  21. [21]

    On splitting methods for Schrödinger-Poisson and cubic nonlinear Schrödinger equations

    C. Lubich. “On splitting methods for Schrödinger-Poisson and cubic nonlinear Schrödinger equations”. In:Math. Comp.77 (2008), pp. 2141–2153

  22. [22]

    Error estimates of a Fourier integrator for the cubic Schrödinger equation at low regularity

    A. Ostermann, F. Rousset, and K. Schratz. “Error estimates of a Fourier integrator for the cubic Schrödinger equation at low regularity”. In:Found. Comput. Math.21 (2021), pp. 725–765

  23. [23]

    Fourier integrator for periodic NLS: low regularity estimates via discrete Bourgain spaces

    A. Ostermann, F. Rousset, and K. Schratz. “Fourier integrator for periodic NLS: low regularity estimates via discrete Bourgain spaces”. In:J. Eur. Math. Soc.25.10 (2023), pp. 3913–3952

  24. [24]

    On uniqueness for semilinear wave equations

    F. Planchon. “On uniqueness for semilinear wave equations”. In:Math. Z.244 (2003), pp. 587–599

  25. [25]

    Averaged moduli of smoothness for functions of several variables and function spaces generated by them

    V. Popov and V. Hristov. “Averaged moduli of smoothness for functions of several variables and function spaces generated by them”. In:Trudy Mat. Inst. Steklov164 (1983), pp. 136–141

  26. [26]

    A general framework of low regularity integrators

    F. Rousset and K. Schratz. “A general framework of low regularity integrators.” In: SIAM J. Numer. Anal.59.3 (2021), pp. 1735–1768. REFERENCES 39

  27. [27]

    Error analysis of the Lie splitting for semilinear wave equations with finite-energy solutions

    M. Ruff and R. Schnaubelt. “Error analysis of the Lie splitting for semilinear wave equations with finite-energy solutions”. In:Discrete Contin. Dyn. Syst.(2025). Early Access. doi: 10.3934/dcds.2025009

  28. [28]

    Sogge.Lectures on Non-Linear Wave Equations

    C. Sogge.Lectures on Non-Linear Wave Equations. International Press, Somerville, 2013

  29. [29]

    A counterexample to an endpoint bilinear Strichartz inequality

    T. Tao. “A counterexample to an endpoint bilinear Strichartz inequality”. In:Electron. J. Differ. Equ.151 (2006), pp. 1–6

  30. [30]

    Tao.Nonlinear Dispersive Equations: Local and Global Analysis

    T. Tao.Nonlinear Dispersive Equations: Local and Global Analysis. American Math- ematical Society, Providence, RI, 2006. Karlsruhe Institute of Technology, Department of Mathematics, Engler- straße 2, 76131 Karlsruhe, Germany Email address: maximilian.ruff@kit.edu