Solving the radial Dirac equations: a numerical odyssey
Abstract
We discuss, in a pedagogical way, how to solve for relativistic wave functions from the radial Dirac equations. After an brief introduction, in Section II we solve the equations for a linear Lorentz scalar potential, Vs(r), that provides for confinement of a quark. The case of massless u and d quarks is treated first, as these are necessarily quite relativistic. We use an iterative procedure to find the eigenenergies and the upper and lower component wave functions for the ground state and then, later, some excited states. Solutions for the massive quarks (s, c, and b) are also presented. In Section III we solve for the case of a Coulomb potential, which is a time-like component of a Lorentz vector potential, Vv(r). We re-derive, numerically, the (analytically well-known) relativistic hydrogen atom eigenenergies and wave functions, and later extend that to the cases of heavier one-electron atoms and muonic atoms. Finally, Section IV finds solutions for a combination of the Vs and Vv potentials. We treat two cases. The first is one in which Vs is the linear potential used in Sec. II and Vv is Coulombic, as in Sec. III. The other is when both Vs and Vv are linearly confining, and we establish when these potentials give a vanishing spin-orbit interaction (as has been shown to be the case in quark models of the hadronic spectrum).
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.