Stabilized finite element method for the radial Dirac equation

A challenging difficulty in solving the radial Dirac eigenvalue problem numerically is the presence of spurious (unphysical) eigenvalues, among the genuine ones, that are neither related to mathematical interpretations nor to physical explanations. Many attempts have been made and several numerical methods have been applied to solve the problem using the finite element method (FEM), the finite difference method, or other numerical schemes. Unfortunately most of these attempts failed to overcome the difficulty. As a FEM approach, this work can be regarded as a first promising scheme to solve the spuriosity problem completely. Our approach is based on an appropriate choice of trial and test function spaces. We develop a Streamline Upwind Petrov-Galerkin method to the equation and derive an explicit stability parameter.

[1]  V. Shabaev,et al.  Spurious states of the dirac equation in a finite basis set , 2008 .

[2]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[3]  Mario A. Storti,et al.  A Petrov-Galerkin formulation for advection-reaction-diffusion problems , 1996 .

[4]  S. Salomonson,et al.  Relativistic all-order pair functions from a discretized single-particle Dirac Hamiltonian. , 1989, Physical review. A, General physics.

[5]  Regina C. Almeida,et al.  A stable Petrov-Galerkin method for convection-dominated problems , 1997 .

[6]  M. Griesemer,et al.  ACCUMULATION OF DISCRETE EIGENVALUES OF THE RADIAL DIRAC OPERATOR , 1999 .

[7]  R. W. Douglass,et al.  The Origin and Nature of Spurious Eigenvalues in the Spectral Tau Method , 1998 .

[8]  Shan Zhao,et al.  On the spurious solutions in the high-order finite difference methods for eigenvalue problems , 2007 .

[9]  Ioan E. Lager,et al.  On the Causes of Spurious Solutions in Electromagnetics , 2002 .

[10]  V M Shabaev,et al.  Dual kinetic balance approach to basis-set expansions for the dirac equation. , 2004, Physical review letters.

[11]  I. Lindgren,et al.  The covariant-evolution-operator method in bound-state QED , 2004 .

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

[13]  Charlotte Froese Fischer,et al.  A B-spline Galerkin method for the Dirac equation , 2008, Comput. Phys. Commun..

[14]  I. Wolff,et al.  The origin of spurious modes in numerical solutions of electromagnetic field eigenvalue problems , 1994 .

[15]  P. E. Long,et al.  Some physical and numerical aspects of boundary layer modeling , 1975 .

[16]  P. A. B. De Sampaio,et al.  A Petrov–Galerkin/modified operator formulation for convection–diffusion problems , 1990 .

[17]  A. Brooks,et al.  A Petrov-Galerkin Finite Element Formulation for Convection Dominated Flows , 1981 .

[18]  Rosenberg.  Virtual-pair effects in atomic structure theory. , 1989, Physical review. A, General physics.

[19]  Werner Scheid,et al.  FINITE ELEMENT FORMULATION OF THE DIRAC EQUATION AND THE PROBLEM OF FERMION DOUBLING , 1998 .

[20]  M. J. P. Cullen.  A Finite Element Method for a Non-linear Initial Value Problem , 1974 .

[21]  Edward Ackad,et al.  Numerical solution of the Dirac equation by a mapped Fourier grid method , 2005 .