The nonlinear Dirac equation is an important model in quantum physics with a set of conservation laws and a multi-symplectic formulation. In this paper, we propose energy-preserving and multi-symplectic wavelet algorithms for this model. Meanwhile, we evidently improve the efficiency of these algorithms in computations via splitting technique and explicit strategy. Numerical experiments are conducted during long-term simulations to show the excellent performances of the proposed algorithms and verify our theoretical analysis.

In this paper, we propose two new explicit multi-symplectic splitting methods for the nonlinear Dirac (NLD) equation. Based on its multi-symplectic formulation, the NLD equation is split into one linear multi-symplectic system and one nonlinear infinite Hamiltonian system. Then multi-symplectic Fourier pseudospectral method and multi-symplectic Preissmann scheme are employed to discretize the linear subproblem, respectively. And the nonlinear subsystem is solved by a symplectic scheme. Finally, a composition method is applied to obtain the final schemes for the NLD equation. We find that the two proposed schemes preserve the total symplecticity and can be solved explicitly. Numerical experiments are presented to show the effectiveness of the proposed methods.

We study the spectral stability of solitary wave solutions to the nonlinear Diracequation in one dimension. We focus on the Dirac equation with cubic nonlinearity, knownas the Soler model in (1+1) dimensions and also as the massive Gross-Neveu model.Presented numerical computations of the spectrum of linearization at a solitary wave showthat the solitary waves are spectrally stable. We corroborate our results by findingexplicit expressions for several of the eigenfunctions. Some of the analytic results holdfor the nonlinear Dirac equation with generic nonlinearity.