pith. sign in

arxiv: 1906.09413 · v1 · pith:Q4ANX3IDnew · submitted 2019-06-22 · 🧮 math.NA · cs.NA

Low-regularity integrators for nonlinear Dirac equations

Katharina Schratz , Yan Wang , Xiaofei Zhao This is my paper

Pith reviewed 2026-05-25 18:21 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords nonlinear Dirac equationlow-regularity integratorDirac-Poisson systemnumerical convergencerough initial datatime integrationSobolev regularity
0
0 comments X

The pith

The ultra low-regularity integrator achieves optimal first-order time convergence in H^r for solutions in H^r without extra regularity.

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

Classical integrators for the nonlinear Dirac equation and Dirac-Poisson system lose derivatives numerically when initial data are rough, forcing the use of smoother solutions to obtain convergence. This paper introduces an ultra low-regularity integrator that works directly with solutions whose regularity is exactly H^r. It proves that the new scheme delivers first-order accuracy in time measured in the same H^r norm, and it supplies an extension to second order. The result matters for physical regimes where solutions naturally sit at low regularity and classical schemes become inefficient or inaccurate. Numerical tests illustrate the improved error behaviour at the lowest admissible regularity.

Core claim

The ultra low-regularity integrator (ULI) for the nonlinear Dirac equation and the Dirac-Poisson system enables optimal first-order time convergence in H^r for solutions belonging to H^r, without requiring any additional regularity on the solution. In contrast to classical methods, ULI overcomes the numerical loss of derivatives and is therefore more efficient and accurate for approximating low-regular solutions. Convergence theorems and the extension of ULI to second order are established.

What carries the argument

The ultra low-regularity integrator (ULI), a time-stepping scheme built on the mild formulation that integrates the nonlinearity while matching the solution's native Sobolev regularity.

If this is right

  • Optimal first-order convergence holds for arbitrary initial data in H^r.
  • The scheme extends to second-order accuracy while retaining the low-regularity property.
  • Numerical experiments confirm the theoretical rates and demonstrate smaller errors than classical integrators at low regularity.
  • The method applies equally to the nonlinear Dirac equation and the Dirac-Poisson system.

Where Pith is reading between the lines

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

  • The same construction could be tested on other semilinear dispersive systems that suffer derivative loss under standard discretizations.
  • Long-time simulations with rough data may become feasible if the integrator also controls growth of constants independent of regularity.
  • Implementation in existing codes would require only a change in the treatment of the Duhamel term rather than a full redesign of the spatial discretisation.

Load-bearing premise

The mild solution formulation of the nonlinear Dirac equation remains well-defined and the nonlinearity satisfies the necessary estimates in H^r without invoking higher regularity.

What would settle it

A computation on an exactly H^r solution that shows the observed convergence rate of ULI dropping below first order while a classical scheme fails even more severely.

Figures

Figures reproduced from arXiv: 1906.09413 by Katharina Schratz, Xiaofei Zhao, Yan Wang.

Figure 1
Figure 1. Figure 1: Convergence of the first-order methods for NDE with external V = Ve (left) and for Dirac-Poisson system (right) under H2.4 -initial data: error = (Φ(tn)− Φ n)/kΦ(tn)kH2 for tn = T = 1. To test the first-order methods, we construct initial data Φ0(x) ∈ H2.4 (T) as in (5.1), and the convergence results of ULI1 (3.10), FD1 (2.2), Lie splitting (2.8) and EI1 (2.5) for solving NDE (2.1) with external V = Ve are… view at source ↗
Figure 2
Figure 2. Figure 2: Convergence of the second-order methods for NDE with external V = Ve (left) and for Dirac-Poisson system (right) under H2.2 -initial data: error = (Φ(tn) − Φ n)/kΦ(tn)kH1 for tn = T = 1. 10-2 10-1 10-2 10-1 100 ULI2 Strang FD2 EI2 O( ) 10-2 10-1 10-2 10-1 100 ULI2 Strang FD2 EI2 O( ) [PITH_FULL_IMAGE:figures/full_fig_p017_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Convergence of the second-order methods for NDE with external V = Ve (left) and for Dirac-Poisson system (right) under H1.4 -initial data: error = (Φ(tn) − Φ n)/kΦ(tn)kH1 for tn = T = 1. methods all suffer from order reduction if the solution does not satisfy the critical regularity require￾ment. We propose a new class of ultra low-regularity integrators (ULI) for solving the NDEs. The great advantage of t… view at source ↗
read the original abstract

In this work, we consider the numerical integration of the nonlinear Dirac equation and the Dirac-Poisson system (NDEs) under rough initial data. We propose a ultra low-regularity integrator (ULI) for solving the NDEs which enables optimal first-order time convergence in $H^r$ for solutions in $H^{r}$, i.e., without requiring any additional regularity on the solution. In contrast to classical methods, ULI overcomes the numerical loss of derivatives and is therefore more efficient and accurate for approximating low regular solutions. Convergence theorems and the extension of ULI to second order are established. Numerical experiments confirm the theoretical results and underline the favourable error behaviour of the new method at low regularity compared to classical integration schemes.

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 manuscript introduces an ultra low-regularity integrator (ULI) for the nonlinear Dirac equation and Dirac-Poisson system. The central claim is that the ULI achieves optimal first-order time convergence in H^r for solutions possessing only H^r regularity, without requiring additional smoothness on the data or solution. Convergence theorems are established, the integrator is extended to second order, and numerical experiments are presented to confirm the rates and demonstrate improved accuracy relative to classical schemes at low regularity.

Significance. If the stated convergence theorems hold under the mild-solution estimates in H^r, the result would be a meaningful advance for low-regularity numerical methods for nonlinear dispersive equations. The avoidance of derivative loss is directly relevant to applications involving rough data in relativistic systems. Credit is due for the explicit statement of the H^r convergence claim, the second-order extension, and the accompanying numerical validation.

minor comments (3)
  1. [Abstract] The abstract and introduction would benefit from a short, explicit statement of the precise form of the nonlinearity and the precise range of r for which the H^r estimates close.
  2. Notation for the Dirac operator and the mild formulation should be collected in a single preliminary subsection to improve readability before the convergence analysis begins.
  3. Figure captions for the numerical experiments should include the precise values of r and the mesh parameters used in each test to allow direct comparison with the theorem statements.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the careful reading of the manuscript, the positive assessment of its significance, and the recommendation for minor revision. No specific major comments appear in the report, so we have no points requiring point-by-point rebuttal or revision at this stage. We will incorporate any minor suggestions that may arise during the revision process.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper constructs the ULI explicitly from the mild formulation of the NDE and derives convergence estimates under the assumption that the mild solution and nonlinearity estimates remain valid in H^r. No load-bearing step equates the claimed first-order convergence in H^r to a fitted quantity, a self-citation chain, or a definitional renaming; the analysis proceeds from the integral equation and standard estimates without reducing the result to its inputs by construction. The derivation is therefore self-contained against the stated premises.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Review performed on abstract alone; the central claim rests on standard functional-analytic assumptions for the Dirac operator and nonlinearity that are not enumerated here.

axioms (1)
  • domain assumption Mild solutions to the nonlinear Dirac equation exist in H^r and satisfy the integral equation without additional regularity
    Required for the convergence statement to apply directly to H^r data

pith-pipeline@v0.9.0 · 5648 in / 1107 out tokens · 26016 ms · 2026-05-25T18:21:55.464590+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

55 extracted references · 55 canonical work pages

  1. [1]

    Alvarez, B

    A. Alvarez, B. Carreras, Interaction dynamics for the solitary waves of a nonlinear Dirac model, Phys. Lett., 86 (1981) pp. 327-332

    work page 1981

  2. [2]

    W. Bao, Y. Cai, X. Jia, Q. Tang , Numerical methods and comparison for the Dirac equation in the nonrela- tivistic limit regime, J. Sci. Comput., 71 (2017), pp. 1094-1134

    work page 2017

  3. [3]

    W. Bao, Y. Cai, X. Jia, J. Yin , Error estimates of numerical methods for the nonlinear Dirac equation in the nonrelativistic limit regime, Sci. China Math., 59 (2016) pp. 1461-1494

    work page 2016

  4. [4]

    W. Bao, X. Li , An efficient and stable numerical method for the Maxwell-Dirac system, J. Comp. Phys., 199 (2004), pp. 663-687

    work page 2004

  5. [5]

    Beals, M

    R. Beals, M. B ´ezard, Low regularity local solutions for field equations, Comm. Partial Differential Equations, 21 (1996), pp. 79-124

    work page 1996

  6. [6]

    Bechouche, N

    P. Bechouche, N. Mauser, F. Poupaud, (Semi)-nonrelativistic limits of the Dirac eqaution with external time- dependent electromagnetic field, Commun. Math. Phys., 197 (1998), pp. 405-425

    work page 1998

  7. [7]

    Bournaveas, Local well-posedness for a nonlinear Dirac equation in spaces of almost critical dimension, Discrete Contin

    N. Bournaveas, Local well-posedness for a nonlinear Dirac equation in spaces of almost critical dimension, Discrete Contin. Dyn. Syst., 20 (2008), pp. 605-616

    work page 2008

  8. [8]

    Bournaveas, Local and global solutions for a nonlinear Dirac system, Adv

    N. Bournaveas, Local and global solutions for a nonlinear Dirac system, Adv. Differential Equations, 9 (2004), pp. 677-698

    work page 2004

  9. [9]

    Bournaveas, D

    N. Bournaveas, D. Gibbeson, Low regularity global solutions of the Dirac-Klein-Gordon equations in one space dimension, Differential Integral Equations, 19 (2006), pp. 211-222

    work page 2006

  10. [10]

    Bournaveas, T

    N. Bournaveas, T. Candy, S. Machihara, Local and global well-posedness for the Chern-Simons-Dirac system in one dimension, Differential Integral Equations, 25 (2012), pp. 699-718

    work page 2012

  11. [11]

    Brinkman, C

    D. Brinkman, C. Heitzinger, P.A. Markowich, A convergent 2D finite-difference scheme for the Dirac-Poisson system and the simulation of graphene, J. Comput. Phys., 257 (2014), pp. 318-332

    work page 2014

  12. [12]

    Y. Cai, Y. Wang, A uniformly accurate (UA) multiscale time integrator pseudospectral method for the nonlinear Dirac equation in the nonrelativistic limit regime, Esaim Math. Model. Numer. Anal., 52 (2018), pp. 543-566

    work page 2018

  13. [13]

    Y. Cai, Y. Wang, Uniformly accurate nested Picard iterative integrators for the Dirac equation in the nonrela- tivistic limit regime, to appear on SIAM J. Numer. Anal. 2019

    work page 2019

  14. [14]

    Cacciafesta, Dispersive properties of the Dirac equation: Strichartz estimates and the nonlinear problem, LAP LAMBERT Academic Publishing, 2012

    F. Cacciafesta, Dispersive properties of the Dirac equation: Strichartz estimates and the nonlinear problem, LAP LAMBERT Academic Publishing, 2012

    work page 2012

  15. [15]

    Candy , Global existence for an L2 critical nonlinear Dirac equation in one dimension, Adv

    T. Candy , Global existence for an L2 critical nonlinear Dirac equation in one dimension, Adv. Differential Equations, 16 (2011), pp. 643-666

    work page 2011

  16. [16]

    D’Ancona, D

    P. D’Ancona, D. Foschi, S. Selberg , Null structure and almost optimal local regularity for the Dirac-Klein- Gordon system, J. Eur. Math. Soc., 9 (2007), pp. 877-899

    work page 2007

  17. [17]

    D’Ancona, D

    P. D’Ancona, D. Foschi, S. Selberg, Local well-posedness below the charge norm for the Dirac-Klein-Gordon system in two space dimensions, J. Hyperbolic Differ. Equ., 4 (2007), pp. 295-330

    work page 2007

  18. [18]

    D’Ancona, D

    P. D’Ancona, D. Foschi, S. Selberg , Null structure and almost optimal local well-posedness of the Maxwell- Dirac system, Amer. J. Math., 132 (2010), pp. 771-839

    work page 2010

  19. [19]

    Delgado, Global solution of the Cauchy problem for the (classical) coupled Maxwell-Dirac and other nonlinear Dirac equations in one space dimension, Proc

    V. Delgado, Global solution of the Cauchy problem for the (classical) coupled Maxwell-Dirac and other nonlinear Dirac equations in one space dimension, Proc. Amer. Math. Soc., 69 (1978), pp. 289-206

    work page 1978

  20. [20]

    Escobedo, L

    M. Escobedo, L. Vega, A semilinear Dirac equation in H s(R3) for s > 1, SIAM J. Math. Anal., 28 (1997), pp. 338-362

    work page 1997

  21. [21]

    Fang, Low regularity solutions for Dirac-Klein-Gordon equations in one space dimension, Electr

    Y. Fang, Low regularity solutions for Dirac-Klein-Gordon equations in one space dimension, Electr. J. Diff. Equations, 102 (2001) pp. 1-19

    work page 2001

  22. [22]

    Faou, Geometric Numerical Integration and Schr¨ odinger Equations, European Math

    E. Faou, Geometric Numerical Integration and Schr¨ odinger Equations, European Math. Soc. Publishing House, Z¨ urich 2012

    work page 2012

  23. [23]

    Fefferman, M.I

    C.L. Fefferman, M.I. Weistein, Honeycomb lattice potentials and Dirac points, J. Amer. Math. Soc., 25 (2012) pp. 1169-1220

    work page 2012

  24. [24]

    Frutos, J.M

    J.D. Frutos, J.M. Sanz-Serna , Split-Step Spectral schemes for nonlinear Dirac systems, J. Comp. Phys., 83 (1989) pp. 407-423

    work page 1989

  25. [25]

    Gautschi , Numerical integration of ordinary differential equations based on trigonometric polynomials, Numer

    W. Gautschi , Numerical integration of ordinary differential equations based on trigonometric polynomials, Numer. Math., 3 (1961), pp. 381-397. LOW-REGULARITY INTEGRATORS FOR NONLINEAR DIRAC EQUATIONS 21

    work page 1961

  26. [26]

    Haddad, L.D

    L.H. Haddad, L.D. Carr , The nonlinear Dirac equation in Bose-Einstein condensates: Foundation and sym- metries. Phys. D, 238 (2009) pp. 1413-1421

    work page 2009

  27. [27]

    Hairer, Ch

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

    work page 2006

  28. [28]

    Hammer, W

    R. Hammer, W. P¨otz, A. Arnold, Single-cone real-space finite difference scheme for the time-dependent Dirac equation J. Comput. Phys., 265 (2014), pp. 50-70

    work page 2014

  29. [29]

    Hasegawa, Bound states of the one-dimensional Dirac equation for scalar and vector double square-well potentials, Physica E 59 (2014) pp

    H. Hasegawa, Bound states of the one-dimensional Dirac equation for scalar and vector double square-well potentials, Physica E 59 (2014) pp. 192-201

    work page 2014

  30. [30]

    Hochbruck, A

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

    work page 2010

  31. [31]

    Hofmanov´a, K

    M. Hofmanov´a, K. Schratz, An exponential-type integrator for the KdV equation, Numer. Math., 136 (2017) pp. 1117-1137

    work page 2017

  32. [32]

    Huang, S

    Z. Huang, S. Jin, P. A. Markowich, C. Sparber, C. Zheng, A time-splitting spectral scheme for the Maxwel- Dirac system, J. Comp. Phys., 208 (2005) pp. 761-789

    work page 2005

  33. [33]

    H. Huh, B. Moon , Low regularity well-posedness for Gross-Neveu equations, Commun. Pure Appl. Anal., 14 (2015) pp. 1903-1913

    work page 2015

  34. [34]

    Kn¨oller, A

    M. Kn¨oller, A. Ostermann, K. Schratz , A Fourier integrator for the cubic nonlinear Schr¨ odinger equation with rough initial data, preprint 2019

    work page 2019

  35. [35]

    Lemou, F

    M. Lemou, F. M´ehats, X. Zhao, Uniformly accurate numerical schemes for the nonlinear Dirac equation in the nonrelativistic limit regime, Commun. Math. Sci., 15 (2017) pp. 1107-1128

    work page 2017

  36. [36]

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

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

    work page 2008

  37. [37]

    Machiharaa, M

    S. Machiharaa, M. Nakamurab, K. Nakanishi, T. Ozawa, Endpoint Strichartz estimates and global solutions for the nonlinear Dirac equation, J. Funct. Anal., 219 (2005) pp. 1-20

    work page 2005

  38. [38]

    Machihara, One dimensional Dirac equation with quadratic nonlinearities

    S. Machihara, One dimensional Dirac equation with quadratic nonlinearities. Discrete Contin. Dyn. Syst., 13 (2005), pp. 277-290

    work page 2005

  39. [39]

    McLachlan, G.R.W

    R.I. McLachlan, G.R.W. Quispel , Splitting methods, Acta Numer., 11 (2002) pp. 341-434

    work page 2002

  40. [40]

    Mocken, C.H

    G.R. Mocken, C.H. Keitel, FFT-split-operator code for solving the Dirac equation in 2+1 dimensions, Comput. Phys. Commun., 178 (2008), pp. 868-882

    work page 2008

  41. [41]

    Ostermann, K

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

    work page 2018

  42. [42]

    Ostermann, F

    A. Ostermann, F. Rousset, K. Schratz , Error estimates of a Fourier integrator for the cubic Schr¨ odinger equation at low regularity, preprint, 2019

    work page 2019

  43. [43]

    Pecher, Low regularity well-posedness for the one-dimensional Dirac-Klein-Gordon system, Electron

    H. Pecher, Low regularity well-posedness for the one-dimensional Dirac-Klein-Gordon system, Electron. J. Diff. Eqns., 150 (2006), pp. 1-13

    work page 2006

  44. [44]

    Pecher, Local well-posedness for the nonlinear Dirac equation in two space dimensions, Commun

    H. Pecher, Local well-posedness for the nonlinear Dirac equation in two space dimensions, Commun. Pure Appl. Anal., 13 (2014) pp. 673-685

    work page 2014

  45. [45]

    P¨otz, Single-cone finite-difference schemes for the (2 + 1)-dimensional Dirac equation in general electromag- netic textures, Phys

    W. P¨otz, Single-cone finite-difference schemes for the (2 + 1)-dimensional Dirac equation in general electromag- netic textures, Phys. Rev. E, 96 (2017) pp. 053312

    work page 2017

  46. [46]

    Thirring, A soluble relativistic field theory, Ann

    W.E. Thirring, A soluble relativistic field theory, Ann. Physics, 3 (1958), pp. 91-112

    work page 1958

  47. [47]

    Selberg, A

    S. Selberg, A. Tesfahun , Low regularity well-posedness of the Dirac-Klein-Gordon equations in one space dimension, Commun. Contemp. Math., 10 (2008) pp. 181-194

    work page 2008

  48. [48]

    Selberg, A

    S. Selberg, A. Tesfahun , Low regularity well-posedness for some nonlinear Dirac equations in one space dimension, Differential Integral Equations, 23 (2010) pp. 265-278

    work page 2010

  49. [49]

    S. Shao, H. Tang , Higher-order accurate Runge-Kutta discontinuous Galerkin methods for a nonlinear Dirac model, Discrete Cont. Dyn. Syst. B, 6 (2006), pp. 623-640

    work page 2006

  50. [50]

    J. Shen, T. Tang, L. Wang, Spectral Methods: Algorithms, Analysis and Applications, Springer, 2011

    work page 2011

  51. [51]

    Soler, Classical, stable, nonlinear spinor field with positive rest energy

    M. Soler, Classical, stable, nonlinear spinor field with positive rest energy. Phys. Rev. D, 1 (1970) pp. 2766-2769

    work page 1970

  52. [52]

    Tao, Nonlinear Dispersive Equations

    T. Tao, Nonlinear Dispersive Equations. Local and Global Analysis. Amer. Math. Soc. Providence 2006

    work page 2006

  53. [53]

    W. Yi, Y. Cai , Optimal error estimates of finite difference time domain methods for the Klein-Gordon-Dirac system, https://doi.org/10.1093/imanum/dry084, IMA J. Numer. Anal., to appear, 2019

    work page doi:10.1093/imanum/dry084 2019

  54. [54]

    W. Yi, X. Ruan, C. Su , Optimal resolution methods for the Klein-Gordon-Dirac system in the nonrelativistic limit regime, J. Sci. Comput., 79 (2019) pp. 1907-1935

    work page 2019

  55. [55]

    J. Xu, S. Shao, H. Tang , Numerical methods for nonlinear Dirac equation, J. Comput. Phys., 245 (2013) pp. 131-149. K. Schratz: Fakult ¨at f ¨ur Mathematik, Karlsruhe Institute of Technology, Englerstr.2, 76131 Karl- sruhe, Germany E-mail address : katharina.schratz@kit.edu Y. Wang: School of Mathematics and Statistics, Central China Normal University, 43...

    work page 2013