$\mathrm{L}^{2}$--convergence of the time-splitting scheme for nonlinear Dirac equation in 1+1 dimensions

TL;DR

This paper proves that the time-splitting numerical scheme for the nonlinear Dirac equation in 1+1 dimensions converges strongly in L^2 to the true solution, ensuring reliable long-term approximations.

Contribution

It establishes strong L^2 convergence of the time-splitting scheme for the nonlinear Dirac equation, using pointwise estimates and a modified Glimm functional for stability.

Findings

Strong L^2 convergence of the scheme to the global solution

Uniform boundedness of the Glimm-type functional over time

Precompactness of the solution set in C([0,T];L^2)

Abstract

We study the time-splitting scheme for approximating solutions to the Cauchy problem of the nonlinear Dirac equation in 1+1 dimensions. Under the assumption that the initial data for the scheme are convergent in $\mathrm{L}^{2}(\mathbb{R})$, we prove that the approximate solutions constructed by the corresponding time-splitting scheme are strongly convergent in $\mathrm{C}([0,\infty);\mathrm{L}^{2}(\mathbb{R}))$ to the global strong solution of the nonlinear Dirac equation. To achieve this, we first establish the pointwise estimates for time-splitting solutions. Based on these estimates, a modified Glimm-type functional is carefully designed to show that it is uniformly bounded in time, which yields $\mathrm{L}^2$ stability estimates for the scheme. Furthermore, we prove that the set of time-splitting solutions is precompact in $\mathrm{C}([0,T];\mathrm{L}^{2}(\mathbb{R}))$ for any $T>0$. Finally, we show that the limit of any subsequence of the time-splitting solutions is the unique strong solution to the Cauchy problem of the nonlinear Dirac equation.