“…Equation (1.1) occurs as a core part of many related models, such as the Gross-Neveu system [33], the Dirac-Klein-Gordon system [16,53], the Maxwell-Dirac system [18,32] and the Chern-Simons-Dirac system [10]. Due to these applications, the numerical solutions of the NDE (1.1) are widely interested in computational physics [4,11,24,32,35,40,45], especially the bound states, the solitary waves and the interaction dynamics of solitons [1,29,55]. These existing works in computational physics or applied mathematics mainly assume that the initial data of NDE is smooth enough for theoretical and/or numerical studies.…”