pith. sign in

arxiv: 2504.03346 · v2 · submitted 2025-04-04 · 🧮 math.NA · cs.NA

Error estimates of an exponential wave integrator for the nonlinear Schr\"odinger equation with singular potential

Weizhu Bao , Chushan Wang This is my paper

Pith reviewed 2026-05-22 21:38 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords nonlinear Schrödinger equationexponential wave integratorsingular potentialerror estimatesStrichartz estimatesCoulomb potentialnumerical methodsquantum mechanics
0
0 comments X

The pith

A first-order exponential wave integrator for the nonlinear Schrödinger equation with singular potential achieves first-order L2 convergence in 3D for the Coulomb potential.

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

This paper establishes error bounds for a first-order exponential wave integrator solving the nonlinear Schrödinger equation when the potential is singular and only locally square-integrable rather than bounded. It proves that with initial data in H2, the L2-norm error is first-order in one and two dimensions but three-quarters order in three dimensions under an L2-potential assumption. With a stronger Lp integrability condition on the potential in three dimensions, the rate improves to nearly first order for p larger than 12/5, and specifically reaches first order for the Coulomb potential. These bounds matter because singular potentials appear in standard models of atoms and molecules, so the analysis supports direct numerical use of the integrator without artificial smoothing of the potential.

Core claim

Under the assumption of an L2-potential and H2-initial data, the L2-norm convergence of the EWI is roughly first-order in 1D and 2D, and 3/4-order in 3D. Under stronger Lp-potential (p>2) in 3D the rate increases to almost 3/4 + 3(1/2 - 1/p) if p≤12/5 and becomes first-order if p>12/5. In particular, first-order L2-norm convergence can be achieved for the Coulomb potential in 3D. The key advancements are the use of discrete Strichartz estimates to handle the loss of integrability due to the singular potential that does not belong to L^∞, and the more favorable local truncation error of the EWI, which requires no spatial smoothness of the potential.

What carries the argument

Discrete (in time) Strichartz estimates combined with the local truncation error of the EWI, which compensates for the potential lacking L^∞ regularity and requires no spatial smoothness of the potential.

If this is right

  • The EWI can be used directly on NLSE problems with unbounded potentials while retaining first-order accuracy in one and two dimensions.
  • In three dimensions the method attains first-order L2 convergence for the Coulomb potential under the stated Lp condition.
  • Numerical simulations of quantum systems involving singular potentials can use these explicit error bounds to control accuracy without regularization.
  • The analysis shows that the choice of exponential integrator avoids extra smoothness requirements on the potential that other schemes might impose.

Where Pith is reading between the lines

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

  • The discrete Strichartz technique might extend to error analysis of other time-stepping methods for dispersive equations with singular coefficients.
  • Similar integrability thresholds could appear when analyzing spatial discretizations of the same equation.
  • The results suggest that three-dimensional quantum simulations with Coulomb potentials need not refine the time step more aggressively than first-order accuracy demands.
  • Extensions to higher-order exponential integrators could follow the same truncation-error advantage.

Load-bearing premise

The potential is locally in L2 (or Lp for p>2 in 3D) and the initial data is in H2, with the proof depending on discrete Strichartz estimates to manage the lack of boundedness.

What would settle it

Compute the observed L2-norm convergence rate for the three-dimensional NLSE with Coulomb potential under the EWI as the time step is successively halved; the rate should approach 1 if the central claim holds.

read the original abstract

We analyze a first-order exponential wave integrator (EWI) for the nonlinear Schr\"odinger equation (NLSE) with a singular potential that is locally in $L^2$, which might be locally unbounded. A typical example is the inverse power potential such as the Coulomb potential, which is the most fundamental potential in quantum physics and chemistry. We prove that, under the assumption of $L^2$-potential and $H^2$-initial data, the $L^2$-norm convergence of the EWI is, roughly, first-order in one dimension (1D) and two dimensions (2D), and $\frac{3}{4}$-order in three dimensions (3D). In addition, under a stronger integrability assumption of $L^p$-potential for some $p>2$ in 3D, the $L^2$-norm convergence increases to almost ${\frac{3}{4}} + 3(\frac{1}{2} - \frac{1}{p})$ order if $p \leq \frac{12}{5}$ and becomes first-order if $p > \frac{12}{5}$. In particular, our results show, to the best of our knowledge for the first time, that first-order $L^2$-norm convergence can be achieved when solving the NLSE with the Coulomb potential in 3D. The key advancements are the use of discrete (in time) Strichartz estimates, which allow us to handle the loss of integrability due to the singular potential that does not belong to $L^\infty$, and the more favorable local truncation error of the EWI, which requires no (spatial) smoothness of the potential. Extensive numerical results in 1D, 2D, and 3D are reported to confirm our error estimates and to show the sharpness of our assumptions on the regularity of the singular potentials.

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

0 major / 3 minor

Summary. The paper analyzes error estimates for a first-order exponential wave integrator (EWI) applied to the nonlinear Schrödinger equation with singular potentials that are locally in L² (but not necessarily L^∞), such as the Coulomb potential. Under L²-potential and H² initial data, it proves L²-norm convergence rates of order 1 in 1D/2D and 3/4 in 3D; under stronger L^p assumptions (p>2) in 3D the rate improves to nearly 3/4 + 3(1/2 - 1/p) for p ≤ 12/5 and reaches order 1 for p > 12/5, including first-order convergence for the Coulomb potential. The proofs rely on discrete (in time) Strichartz estimates to handle the loss of integrability and on the EWI's local truncation error that requires no spatial smoothness of the potential. Numerical experiments in 1D–3D confirm the rates and the sharpness of the regularity assumptions.

Significance. If the results hold, the work is significant for providing the first rigorous L²-error analysis achieving first-order convergence for the EWI on the 3D NLSE with Coulomb potential, a setting central to quantum physics and chemistry. The key technical contributions—discrete Strichartz estimates compensating for the potential's lack of L^∞ regularity and the EWI truncation error avoiding spatial smoothness requirements—enable dimension-dependent rates that are precisely characterized and shown to be sharp. The numerical validation of both the rates and the necessity of the assumptions adds substantial value. This advances the applicability of exponential integrators to singular-potential problems without artificial regularity assumptions.

minor comments (3)
  1. Abstract: the phrase 'roughly first-order' in 1D/2D should be replaced by the precise order (including any logarithmic factors or constants) that appears in the main theorem statement.
  2. The transition value p = 12/5 in 3D is stated without derivation in the abstract; a brief indication of how this threshold arises from the discrete Strichartz estimates would improve readability.
  3. Numerical section: the reported convergence tables should include the exact observed orders (with computed slopes) alongside the theoretical predictions to allow direct comparison.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive and accurate summary of our manuscript, as well as for recognizing its significance in providing the first rigorous L²-error analysis achieving first-order convergence for the EWI on the 3D NLSE with Coulomb potential. We appreciate the recommendation for minor revision. No major comments were raised in the report.

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

This is a standard mathematical error-analysis paper deriving L2-norm convergence rates for the EWI scheme under explicit assumptions on the potential (locally L2 or Lp for p>2 in 3D) and initial data (H2). The central claims follow from discrete Strichartz estimates and local truncation error bounds that are established directly from the stated regularity; no fitted parameters are renamed as predictions, no self-definitional reductions appear, and no load-bearing step collapses to a self-citation chain or ansatz smuggled from prior author work. The derivation is self-contained against the paper's own assumptions and external benchmarks, with numerics serving only as confirmation rather than input.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The analysis rests on standard tools from dispersive PDE theory adapted to the discrete setting; no free parameters are fitted, no new entities are postulated, and the assumptions are domain-standard rather than ad-hoc inventions of this paper.

axioms (2)
  • domain assumption Discrete (in time) Strichartz estimates hold for the exponential wave integrator
    Invoked to control the nonlinear term when the potential lacks L^∞ regularity
  • domain assumption The local truncation error of the EWI requires no spatial smoothness of the potential
    Stated as a key property enabling the result for singular potentials

pith-pipeline@v0.9.0 · 5887 in / 1623 out tokens · 55929 ms · 2026-05-22T21:38:37.369354+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 2 Pith papers

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

  1. Optimal error bounds on the exponential wave integrator for nonlinear Schr\"odinger equations with highly singular potential

    math.NA 2026-05 unverdicted novelty 7.0

    Optimal first-order L2 convergence for the first-order EWI on NLSE is established for L^p_loc potentials with p > d/2, reaching the well-posedness threshold for the first time in 3D.

  2. Optimal error bounds on the exponential wave integrator for nonlinear Schr\"odinger equations with highly singular potential

    math.NA 2026-05 unverdicted novelty 7.0

    Optimal first-order L2 convergence is proven for the exponential wave integrator on NLSE with L^p_loc potentials down to the well-posedness threshold, with reduced orders for more singular cases.

Reference graph

Works this paper leans on

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

  1. [1]

    Alama Bronsard, Y

    Y. Alama Bronsard, Y. Bruned, and K. Schratz , Low regularity integrators via decorated trees, 2022, arXiv:2202.01171

    work page arXiv 2022

  2. [2]

    P. W. Anderson , Absence of diffusion in certain random lattices , Phys. Rev., 109 (1958), pp. 1492–1505

    work page 1958

  3. [3]

    Antoine, W

    X. Antoine, W. Bao, and C. Besse, Computational methods for the dynamics of the nonlinear Schr¨ odinger/Gross-Pitaevskii equations, Comput. Phys. Commun., 184 (2013), pp. 2621– 2633

    work page 2013

  4. [4]

    Bao and Y

    W. Bao and Y. Cai , Mathematical theory and numerical methods for Bose-Einstein conden- sation, Kinet. Relat. Models, 6 (2013), pp. 1–135

    work page 2013

  5. [5]

    Bao and Y

    W. Bao and Y. Cai , Optimal error estimates of finite difference methods for the Gross- Pitaevskii equation with angular momentum rotation, Math. Comp., 82 (2013), pp. 99–128

    work page 2013

  6. [6]

    W. Bao, D. Jaksch, and P. A. Markowich , Numerical solution of the Gross-Pitaevskii equation for Bose-Einstein condensation , J. Comput. Phys., 187 (2003), pp. 318–342

    work page 2003

  7. [7]

    W. Bao, B. Lin, Y. Ma, and C. Wang , An extended Fourier pseudospectral method for the Gross-Pitaevskii equation with low regularity potential , East Asian J. Appl. Math., 14 (2024), pp. 530–550

    work page 2024

  8. [8]

    W. Bao, Y. Ma, and C. Wang , Optimal error bounds on time-splitting methods for the non- linear Schr¨ odinger equation with low regularity potential and nonlinearity, Math. Models Methods Appl. Sci., 34 (2024), pp. 803–844

    work page 2024

  9. [9]

    Bao and C

    W. Bao and C. Wang , Error estimates of the time-splitting methods for the nonlinear Schr¨ odinger equation with semi-smooth nonlinearity, Math. Comp., 93 (2024), pp. 1599– 1631

    work page 2024

  10. [10]

    Bao and C

    W. Bao and C. Wang, An explicit and symmetric exponential wave integrator for the nonlinear Schr¨ odinger equation with low regularity potential and nonlinearity, SIAM J. Numer. Anal., 62 (2024), pp. 1901–1928

    work page 2024

  11. [11]

    Bao and C

    W. Bao and C. Wang , Optimal error bounds on the exponential wave integrator for the nonlinear Schr¨ odinger equation with low regularity potential and nonlinearity , SIAM J. Numer. Anal., 62 (2024), pp. 93–118

    work page 2024

  12. [12]

    Becker, N

    S. Becker, N. Galke, R. Salzmann, and L. van Luijk , Convergence rates for the Trotter- Kato splitting, 2024, arXiv:2407.04045

    work page arXiv 2024

  13. [13]

    Besse, B

    C. Besse, B. Bid ´egaray, and S. Descombes , Order estimates in time of splitting methods for the nonlinear Schr¨ odinger equation, SIAM J. Numer. Anal., 40 (2002), pp. 26–40

    work page 2002

  14. [14]

    Bruned and K

    Y. Bruned and K. Schratz , Resonance-based schemes for dispersive equations via decorated trees, Forum Math. Pi, 10 (2022), e2, 76

    work page 2022

  15. [15]

    Burgarth, P

    D. Burgarth, P. Facchi, A. Hahn, M. Johnsson, and K. Yuasa , Strong error bounds for Trotter and Strang-splittings and their implications for quantum chemistry , Phys. Rev. Res., 6 (2024), p. 043155

    work page 2024

  16. [16]

    Carles and C

    R. Carles and C. Su , Scattering and uniform in time error estimates for splitting method in NLS, Found. Comput. Math., 24 (2024), pp. 683–722

    work page 2024

  17. [17]

    K. M. Case, Singular potentials, Phys. Rev., 80 (1950), pp. 797–806

    work page 1950

  18. [18]

    Cazenave, Semilinear Schr¨ odinger Equations, Cour

    T. Cazenave, Semilinear Schr¨ odinger Equations, Cour. Lect. Notes Math. 10, Courant Insti- tute of Mathematical Sciences, New York, 2003

    work page 2003

  19. [19]

    Celledoni, D

    E. Celledoni, D. Cohen, and B. Owren , Symmetric exponential integrators with an appli- cation to the cubic Schr¨ odinger equation, Found. Comput. Math., 8 (2008), pp. 303–317

    work page 2008

  20. [20]

    H. J. Choi, S. Kim, and Y. Koh , Time splitting method for nonlinear Schr¨ odinger equation with rough initial data in L2, J. Differential Equations, 417 (2025), pp. 164–190

    work page 2025

  21. [21]

    Choi and Y

    W. Choi and Y. Koh , On the splitting method for the nonlinear Schr¨ odinger equation with initial data in H 1, Discrete Contin. Dyn. Syst., 41 (2021), pp. 3837–3867

    work page 2021

  22. [22]

    Erd˝os, B

    L. Erd˝os, B. Schlein, and H.-T. Yau, Derivation of the cubic non-linear Schr¨ odinger equation from quantum dynamics of many-body systems , Invent. Math., 167 (2007), pp. 515–614

    work page 2007

  23. [23]

    Fibich , The Nonlinear Schr¨ odinger Equation: Singular Solutions and Optical Collapse , Springer, Cham, 2015

    G. Fibich , The Nonlinear Schr¨ odinger Equation: Singular Solutions and Optical Collapse , Springer, Cham, 2015

    work page 2015

  24. [24]

    W. M. Frank, D. J. Land, and R. M. Spector , Singular potentials, Rev. Mod. Phys., 43 (1971), pp. 36–98

    work page 1971

  25. [25]

    Gross and I

    C. Gross and I. Bloch, Quantum simulations with ultracold atoms in optical lattices , Science, 357 (2017), pp. 995–1001. ERROR ESTIMATES FOR NLSE WITH SINGULAR POTENTIAL 21

    work page 2017

  26. [26]

    Henning and D

    P. Henning and D. Peterseim , Crank-Nicolson Galerkin approximations to nonlinear Schr¨ odinger equations with rough potentials, Math. Models Methods Appl. Sci., 27 (2017), pp. 2147–2184

    work page 2017

  27. [27]

    Hochbruck and A

    M. Hochbruck and A. Ostermann, Exponential integrators, Acta Numer., 19 (2010), pp. 209– 286

    work page 2010

  28. [28]

    L. I. Ignat, A splitting method for the nonlinear Schr¨ odinger equation, J. Differential Equa- tions, 250 (2011), pp. 3022–3046

    work page 2011

  29. [29]

    Jahnke and C

    T. Jahnke and C. Lubich , Error bounds for exponential operator splittings , BIT, 40 (2000), pp. 735–744

    work page 2000

  30. [30]

    L. Ji, A. Ostermann, F. Rousset, and K. Schratz, Low regularity full error estimates for the cubic nonlinear Schr¨ odinger equation, SIAM J. Numer. Anal., 62 (2024), pp. 2071–2086

    work page 2024

  31. [31]

    Lubich , On splitting methods for Schr¨ odinger-Poisson and cubic nonlinear Schr¨ odinger equations, Math

    C. Lubich , On splitting methods for Schr¨ odinger-Poisson and cubic nonlinear Schr¨ odinger equations, Math. Comp., 77 (2008), pp. 2141–2153

    work page 2008

  32. [32]

    N. J. Mauser, Y. Wu, and X. Zhao , The cubic nonlinear Schr¨ odinger equation with rough potential, 2024, arXiv:2403.16772

    work page arXiv 2024

  33. [33]

    Meetz, Singular potentials in nonrelativistic quantum mechanics , Il Nuovo Cimento (1955- 1965), 34 (1964), pp

    K. Meetz, Singular potentials in nonrelativistic quantum mechanics , Il Nuovo Cimento (1955- 1965), 34 (1964), pp. 690–708

    work page 1955

  34. [34]

    Nisoli and A

    C. Nisoli and A. R. Bishop , Attractive inverse square potential, u(1) gauge, and winding transitions, Phys. Rev. Lett., 112 (2014), p. 070401

    work page 2014

  35. [35]

    Ostermann, F

    A. Ostermann, F. Rousset, and K. Schratz , Error estimates at low regularity of splitting schemes for NLS , Math. Comp., 91 (2021), pp. 169–182

    work page 2021

  36. [36]

    Ostermann, F

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

    work page 2021

  37. [37]

    Ostermann, F

    A. Ostermann, F. Rousset, and K. Schratz, Fourier integrator for periodic NLS: low regu- larity estimates via discrete Bourgain spaces, J. Eur. Math. Soc., 25 (2023), pp. 3913–3952

    work page 2023

  38. [38]

    Ostermann and K

    A. Ostermann and K. Schratz , Low regularity exponential-type integrators for semilinear Schr¨ odinger equations, Found. Comput. Math., 18 (2018), pp. 731–755

    work page 2018

  39. [39]

    Pitaevskii and S

    L. Pitaevskii and S. Stringari , Bose-Einstein Condensation and Superfluidity , Oxford Uni- versity Press, 2016

    work page 2016

  40. [40]

    Rousset and K

    F. Rousset and K. Schratz , Resonances as a computational tool , Found. Comput. Math., DOI:s10208-024-09665-8

    work page

  41. [41]

    J. M. Sanz-Serna, Methods for the numerical solution of the nonlinear Schr¨ odinger equation, Math. Comp., 43 (1984), pp. 21–27

    work page 1984

  42. [42]

    Schreiber, S

    M. Schreiber, S. S. Hodgman, P. Bordia, H. P. L¨uschen, M. H. Fischer, R. Vosk, E. Alt- man, U. Schneider, and I. Bloch , Observation of many-body localization of interacting fermions in a quasirandom optical lattice , Science, 349 (2015), pp. 842–845

    work page 2015

  43. [43]

    Schr ¨odinger, Quantisierung als eigenwertproblem , Annalen der Physik, 384 (1926), pp

    E. Schr ¨odinger, Quantisierung als eigenwertproblem , Annalen der Physik, 384 (1926), pp. 361–376

    work page 1926

  44. [44]

    Shao and L

    S. Shao and L. Su , Nonlocalization of singular potentials in quantum dynamics , J. Comput. Electron., 22 (2023), pp. 930–945

    work page 2023

  45. [45]

    Sulem and P.-L

    C. Sulem and P.-L. Sulem , The Nonlinear Schr¨ odinger Equation: Self-focusing and Wave Collapse, Springer, 1999

    work page 1999

  46. [46]

    Yukawa, On the interaction of elementary particles , Prog

    H. Yukawa, On the interaction of elementary particles , Prog. Theor. Phys. Suppl., 1 (1955), pp. 1–10

    work page 1955

  47. [47]

    X. Zhou, Y. Cai, X. Tang, and G. Xu , Accurate and efficient numerical methods for the nonlinear Schr¨ odinger equation with Dirac delta potential, Calcolo, 60 (2023), Paper No. 57

    work page 2023