The so-called Dirac oscillator was proposed as a modification of the free Dirac equation which reproduces many of the properties of the simple harmonic oscillator but accompanied by a strong spin-orbit coupling term. It has yet to be extended successfully to the arbitrary spin S case primarily because of the unwieldiness of general spin Lorentz invariant wave equations. It is shown here using the formalism of totally symmetric multispinors that the Dirac oscillator can, however, be made to accommodate spin by incorporating it into the framework of Galilean relativity. This is done explicitly for spin zero and spin one as special cases of the arbitrary spin result. For the general case it is shown that the coefficient of the spin-orbit term has a 1/S behavior by techniques which are virtually identical to those employed in the derivation of the g-factor carried out over four decades ago.PACS numbers: 03.65.Pm; 03.65.-w
I INTRODUCTIONIn 1989 Moshinski and Szczepaniak [1] proposed a formulation of the simple harmonic oscillator (SHO) which was based on the Dirac equation. They argued that since the first order differential Dirac equation was linear in the momentum p but led to a quadratic dependence on p in the relativistic energy-momentum relation, it was natural to expect that the quadratic dependence on the spatial coordinate in the SHO could emerge by a linear dependence on it in the corresponding Dirac equation. Thus they proposed that the momentum four-vector p µ (µ = 0, 1, 2, 3) be replaced by (p − iMωrβ, p 0 ) where M is the particle mass, ω a frequency parameter, and β the zero component of the Dirac set γ µ . This leads one to consider the modified Dirac equationfor the four component wave function ψ in units in whichh = c = 1. Upon writing ψ in terms of two two-component wave functions ψ 1 and ψ 2 it is readily shown that ψ 1 satisfies a Schrödinger equation with eigenvalues given bywhere n = 0, 1, 2, ..., ℓ is the angular momentum quantum number, L is the angular momentum operator, and σ are the usual Pauli spin matrices. In the low energy limit the left hand side of this relation becomes the nonrelativistic energy E and one has E = ω[2n + (ℓ + 1) 2 − (j + 1 2 ) 2 ]where j = ℓ ± 1 2 is the usual total angular momentum quantum number. .