An efficient and stable numerical method for the Maxwell-Dirac system

In this paper, we present an explicit, unconditionally stable and accurate numerical method for the Maxwell-Dirac system (MD) and use it to study dynamics of MD. As preparatory steps, we take the three-dimensional (3D) Maxwell Dirac system, scale it to obtain a two-parameter model and review plane wave solution of free MD. Then we present a time-splitting spectral method (TSSP) for MD. The key point in the numerical method is based on a time-splitting discretization of the Dirac system, and to discretize nonlinear wave-type equations by pseudospectral method for spatial derivatives, and then solving the ordinary differential equations (ODEs) in phase space analytically under appropriate chosen transmission conditions between different time intervals. The method is explicit, unconditionally stable, time reversible, time transverse invariant, and of spectral-order accuracy in space and second-order accuracy in time. Moreover, it conserves the particle density exactly in discretized level and gives exact results for plane wave solution of free MD. Extensive numerical tests are presented to confirm the above properties of the numerical method. Furthermore, the tests also suggest the following meshing strategy (or e-resolution) is admissible in the 'nonrelativistic' limit regime (0 < e ≪ 1): spatial mesh size h = O(e) and time step Δt = O(e2), where the parameter e is inversely proportional to the speed of light.

[1]  Anadijiban Das,et al.  A class of exact plane wave solutions of the Maxwell–Dirac equations , 1989 .

[2]  N. Mauser,et al.  The Selfconsistent Pauli Equation , 2001 .

[3]  Semiclassical asymptotics for the Maxwell–Dirac system , 2003, math-ph/0305009.

[4]  P. Dirac.  Principles of Quantum Mechanics , 1982 .

[5]  A. Lisi.  A solitary wave solution of the Maxwell-Dirac equations , 1995 .

[6]  Simonetta Abenda,et al.  Solitary waves for Maxwell-Dirac and Coulomb-Dirac models , 1998 .

[7]  Clotilde Fermanian Kammerer,et al.  Semi-Classical Analysis of a Dirac Equation Without Adiabatic Decoupling , 2004 .

[8]  B. Najman,et al.  The nonrelativistic limit of the nonlinear Dirac equation , 1992 .

[9]  P. Markowich,et al.  Numerical solution of the Gross--Pitaevskii equation for Bose--Einstein condensation , 2003, cond-mat/0303239.

[10]  P D Jarvis,et al.  Algebraic solution for the vector potential in the Dirac equation , 2001 .

[11]  Anadijiban Das,et al.  General solutions of Maxwell-Dirac equations in 1+1-dimensional space-time and a spatially confined solution , 1993 .

[12]  Weizhu Bao,et al.  An Explicit Unconditionally Stable Numerical Method for Solving Damped Nonlinear Schrödinger Equations with a Focusing Nonlinearity , 2003, SIAM Journal on Numerical Analysis.

[13]  Maria J. Esteban,et al.  An overview on linear and nonlinear Dirac equations , 2002 .

[14]  V. Gueorguiev.  Small amplitude solutions of the Maxwell - Dirac equations , 1991 .

[15]  J. Simon,et al.  Asymptotic Completeness, Global Existence and the Infrared Problem for the Maxwell-Dirac Equations , 1995, hep-th/9502061.

[16]  Stefan Keppeler,et al.  A Semiclassical approach to the Dirac equation , 1999 .

[17]  P. Markowich,et al.  On time-splitting spectral approximations for the Schrödinger equation in the semiclassical regime , 2002 .

[18]  Nikolaos Bournaveas,et al.  Local existence for the Maxwell-Dirac equations in three space dimensions , 1996 .

[19]  Herbert Spohn.  Semiclassical Limit of the Dirac Equation and Spin Precession , 2000 .

[20]  J. Chadam.  Global solutions of the Cauchy problem for the (classical) coupled Maxwell-Dirac equations in one space dimension , 1973 .

[21]  Fangfang Sun,et al.  Efficient and Stable Numerical Methods for the Generalized and Vector Zakharov System , 2005, SIAM J. Sci. Comput..

[22]  N. Mauser,et al.  (Semi)-Nonrelativistic Limits of the Dirac Equation with External Time-Dependent Electromagnetic Field , 1998 .

[23]  Shi Jin,et al.  Numerical Study of Time-Splitting Spectral Discretizations of Nonlinear Schrödinger Equations in the Semiclassical Regimes , 2003, SIAM J. Sci. Comput..

[24]  Bernd Thaller,et al.  The Dirac Equation , 1992 .

[25]  A. Das.  An ongoing big bang model in the special relativistic Maxwell–Dirac equations , 1996 .

[26]  Walter Hunziker,et al.  On the nonrelativistic limit of the Dirac theory , 1975 .

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