Uniformly accurate numerical schemes for the nonlinear dirac equation in the nonrelativistic limit regime

Mohammed Lemou (1, 2) , Florian Méhats (2, 1) , Xiaofei Zhao (2)

1 IPSO - Invariant Preserving SOlvers
2 IRMAR - Institut de Recherche Mathématique de Rennes

We apply the two-scale formulation approach to propose uniformly accurate (UA) schemes for solving the nonlinear Dirac equation in the nonrelativistic limit regime. The nonlinear Dirac equation involves two small scales \epsilon and \epsilon^2 with epsilon \to 0 in the nonrelativistic limit regime. The small parameter causes high oscillations in time which brings severe numerical burden for classical numerical methods. We transform our original problem as a two-scale formulation and present a general strategy to tackle a class of highly oscillatory problems involving the two small scales \epsilon and\epsilon^2. Suitable initial data for the two-scale formulation is derived to bound the time derivatives of the augmented solution. Numerical schemes with uniform (with respect to \epsilon\in (0;1]) spectral accuracy in space and uniform first order or second order accuracy in time are proposed. Numerical experiments are done to confirm the UA property.

Mots clés: nonrelativistic limit; uniform accuracy; nonlinear Dirac equation; two-scale formulation; highly oscillatory equations

Domaines: Mathématiques [math]

Mohammed Lemou, Florian Méhats, Xiaofei Zhao. Uniformly accurate numerical schemes for the nonlinear dirac equation in the nonrelativistic limit regime. Communications in Mathematical Sciences, 2017, 15 (4), pp.1107-1128. ⟨10.4310/CMS.2017.v15.n4.a9⟩. ⟨hal-01313976v2⟩