Study of instability of the Fourier split-step method for the massive Gross-Neveu model

Stability properties of the well-known Fourier split-step method used to simulate a soliton and similar solutions of the nonlinear Dirac equations, known as the Gross--Neveu model, are studied numerically and analytically. Three distinct types of numerical instability that can occur in this case, are revealed and explained. While one of these types can be viewed as being related to the numerical instability occurring in simulations of the nonlinear Schrodinger equation, the other two types have not been studied or identified before, to the best of our knowledge. These two types of instability are {\em unconditional}, i.e. occur for arbitrarily small values of the time step. They also persist in the continuum limit, i.e. for arbitrarily fine spatial discretization. Moreover, one of them persists in the limit of an infinitely large computational domain. It is further demonstrated that similar instabilities also occur for other numerical methods applied to the Gross--Neveu soliton, as well as to certain solitons of another relativistic field theory model, the massive Thirring.

[1]  E. V. Zemlyanaya,et al.  Vibrations and Oscillatory Instabilities of Gap Solitons , 1998 .

[2]  Sihong Shao,et al.  Numerical methods for nonlinear Dirac equation , 2012, J. Comput. Phys..

[3]  A. Khare,et al.  Approximate analytic solutions to coupled nonlinear Dirac equations , 2016, 1603.08043.

[4]  Thiab R. Taha,et al.  Analytical and numerical aspects of certain nonlinear evolution equations. I. Analytical , 1984 .

[5]  T. Lakoba.  Numerical study of solitary wave stability in cubic nonlinear Dirac equations in 1D , 2018 .

[6]  A. Newell,et al.  On the coleman correspondence and the solution of the massive thirring model , 1977 .

[7]  Erwan Faou,et al.  Birkhoff normal form for splitting methods applied to semilinear Hamiltonian PDEs. Part I. Finite-dimensional discretization , 2009, Numerische Mathematik.

[8]  S. Orfanidis,et al.  Soliton solutions of the massive thirring model , 1975 .

[9]  S. D. Ellis,et al.  Chiral Confinement: An Exact Solution of the Massive Thirring Model , 1975 .

[10]  G. Knorr,et al.  The integration of the vlasov equation in configuration space , 1976 .

[11]  P. Markowich,et al.  Numerical simulation of a generalized Zakharov system , 2004 .

[12]  Katharina Schratz,et al.  Splitting methods for nonlinear Dirac equations with Thirring type interaction in the nonrelativistic limit regime , 2021, J. Comput. Appl. Math..

[13]  Ludwig Gauckler,et al.  Splitting Integrators for Nonlinear Schrödinger Equations Over Long Times , 2010, Found. Comput. Math..

[14]  A. Saxena,et al.  Solitary waves in the Nonlinear Dirac Equation , 2017, 1707.01946.

[15]  T. Taha,et al.  Analytical and numerical aspects of certain nonlinear evolution equations. II. Numerical, nonlinear Schrödinger equation , 1984 .

[16]  Guo-Wei Wei,et al.  Numerical methods for the generalized Zakharov system , 2003 .

[17]  S. Y. Lee,et al.  Exact localized solutions of two-dimensional field theories of massive fermions with Fermi interactions , 1975 .

[18]  J. Atai,et al.  Bragg solitons in systems with separated nonuniform Bragg grating and nonlinearity. , 2017, Physical review. E.

[19]  R. H. Hardin.  Application of the split-step Fourier method to the numerical solution of nonlinear and variable coefficient wave equations , 1973 .

[20]  H. Yoshida.  Construction of higher order symplectic integrators , 1990 .

[21]  B. Malomed,et al.  Slowdown and splitting of gap solitons in apodized bragg gratings , 2004, nlin/0405036.

[22]  E. Kuznetsov,et al.  On the complete integrability of the two-dimensional classical Thirring model , 1977 .

[23]  B. Herbst,et al.  Split-step methods for the solution of the nonlinear Schro¨dinger equation , 1986 .

[24]  J. M. Sanz-Serna,et al.  Split-step spectral schemes for nonlinear Dirac systems , 1989 .

[25]  Fernando Vadillo,et al.  An integrating factor for nonlinear Dirac equations , 2010, Comput. Phys. Commun..

[26]  Yan Wang,et al.  A uniformly accurate (UA) multiscale time integrator pseudospectral method for the nonlinear Dirac equation in the nonrelativistic limit regime , 2018 .

[27]  W. Bao,et al.  Error estimates of numerical methods for the nonlinear Dirac equation in the nonrelativistic limit regime , 2015, 1511.01192.

[28]  Hanquan Wang,et al.  An efficient and spectrally accurate numerical method for computing dynamics of rotating Bose-Einstein condensates , 2006, J. Comput. Phys..

[29]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[30]  Herbert G. Winful,et al.  Pulse compression in optical fiber filters , 1985 .

[31]  Stability of solitary waves in the nonlinear Dirac equation with arbitrary nonlinearity. , 2014, Physical review. E, Statistical, nonlinear, and soft matter physics.

[32]  I. V. Barashenkov,et al.  Spinor solitons and their PT-symmetric offspring , 2018, Annals of Physics.

[33]  T. Lakoba.  Instability of the finite‐difference split‐step method applied to the nonlinear Schrödinger equation. II. moving soliton , 2016 .

[34]  J. Atai,et al.  Moving Bragg grating solitons in a semilinear dual-core system with dispersive reflectivity , 2017, Scientific Reports.

[35]  Xianggui Li,et al.  Time‐splitting methods with charge conservation for the nonlinear Dirac equation , 2017 .

[36]  V. Zakharov.  Collapse of Langmuir Waves , 1972 .

[37]  B. Carreras,et al.  Interaction dynamics for the solitary waves of a nonlinear Dirac model , 1981 .

[38]  M. Suzuki,et al.  Fractal decomposition of exponential operators with applications to many-body theories and Monte Carlo simulations , 1990 .

[39]  Xiaowei Jia,et al.  Numerical Methods and Comparison for the Dirac Equation in the Nonrelativistic Limit Regime , 2015, J. Sci. Comput..

[40]  G. Stewart.  Introduction to matrix computations , 1973 .

[41]  B. Hermansson,et al.  High-order generalized propagation techniques , 1991 .

[42]  David J. Kaup,et al.  The squared eigenfunctions of the massive Thirring model in laboratory coordinates , 1996 .

[43]  M. Ablowitz,et al.  Analytical and Numerical Aspects of Certain Nonlinear Evolution Equations , 1984 .

[44]  X Antoine,et al.  A simple pseudospectral method for the computation of the time-dependent Dirac equation with Perfectly Matched Layers , 2018, J. Comput. Phys..

[45]  S. Cox,et al.  Exponential Time Differencing for Stiff Systems , 2002 .

[46]  Xiaofei Zhao,et al.  Uniformly accurate numerical schemes for the nonlinear Dirac equation in the nonrelativistic limit regime , 2016, 1605.02475.

[47]  André D. Bandrauk,et al.  Exponential split operator methods for solving coupled time-dependent Schrödinger equations , 1993 .

[48]  Natalia M. Litchinitser,et al.  Pulse compression using fiber gratings as highly dispersive nonlinear elements , 1998 .

[49]  Weizhu Bao,et al.  An efficient and stable numerical method for the Maxwell-Dirac system , 2004 .

[50]  G. Berkolaiko,et al.  On Spectral Stability of Solitary Waves of Nonlinear Dirac Equation in 1D , 2009, 0910.0917.