Dirac Oscillator: Exact Relativistic Model

• The Dirac oscillator is an exactly solvable relativistic model defined by a non-minimal momentum substitution that produces harmonic oscillator behavior with strong spin–orbit coupling.
• It exhibits distinctive positive- and negative-energy branches with spectral asymmetries that provide deep insights into relativistic scattering, nuclear structure, and spin dynamics.
• The model’s mapping to Jaynes–Cummings systems and its experimental analogs in photonic and microwave platforms enable practical quantum simulations and applications in quantum optics.

The Dirac oscillator is an exactly solvable relativistic bound-state model obtained from the free Dirac equation by a non-minimal substitution of the momentum operator, conventionally written as $$\mathbf{p}\to \mathbf{p}-i m\omega\,\beta\,\mathbf{r}$$. It is linear in both position and momentum, and in the nonrelativistic limit it yields a harmonic oscillator supplemented by a strong spin–orbit coupling. Across the literature, it appears as a paradigm of relativistic quantum mechanics, a bridge to quantum optics through Jaynes–Cummings-type mappings, a basis for nuclear and scattering calculations, and an experimentally accessible analog system in photonic and microwave platforms (Quimbay et al., 2012, Tsutsui et al., 21 Dec 2025, Franco-Villafañe et al., 2013).

1. Definition and operator structure

The free Dirac Hamiltonian is written as

$$H_{\mathrm{free}}=\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m,$$

or, in $(3+1)$ dimensions,

$$i\hbar \frac{\partial}{\partial t}\Psi(\mathbf{r},t)=\left[c\,\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m c^2\right]\Psi(\mathbf{r},t).$$

The Dirac oscillator is defined by the non-minimal replacement

$$\mathbf{p}\longrightarrow \mathbf{p}-i m\omega\,\beta\,\mathbf{r},$$

which gives

$$H_{\mathrm{DO}}=\boldsymbol{\alpha}\cdot\left(\mathbf{p}-i m\omega\,\beta\,\mathbf{r}\right)+\beta m$$

in units $\hbar=c=1$, or the equivalent $c$-dependent form in standard units (Quimbay et al., 2012, Boumali et al., 25 Feb 2026).

A central reason for the model’s importance is that decoupling the Dirac equation into large and small components produces a three-dimensional isotropic harmonic-oscillator operator together with a strong spin–orbit term. One explicit form is

$$\left(E^2-m^2c^4\right)\phi= \left[ c^2\left(\mathbf{p}^{\,2}+m^2\omega^2 r^2\right)-3mc^2\hbar\omega-\frac{4mc^2\omega}{\hbar}\,\mathbf{L}\cdot\mathbf{S} \right]\phi,$$

so the model is not merely a relativistic oscillator in name; its squared Hamiltonian contains the oscillator and spin–orbit structures exactly (Boumali et al., 25 Feb 2026).

In one spatial dimension, equivalent formulations use Pauli matrices directly. A representative Hamiltonian is

$$\mathcal{H}=\sigma_y\,(p-i\sigma_z m\omega x)+\sigma_z m,$$

while in the photonic realization a Dirac-like equation is written as

$$H_{\mathrm{free}}=\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m,$$0

with the Dirac-oscillator choice

$$H_{\mathrm{free}}=\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m,$$1

This realizes the oscillator as a constant mass plus a linear non-minimal coupling (Cao et al., 2019, Longhi, 2010).

2. Exact solvability and spectral structure

The one-dimensional Dirac oscillator has a discrete relativistic spectrum with separate positive- and negative-energy branches. In the formulation used for the fiber-Bragg-grating analog, the stationary-state problem reduces to harmonic-oscillator equations for each spinor component,

$$H_{\mathrm{free}}=\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m,$$2

with eigenvalues

$$H_{\mathrm{free}}=\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m,$$3

A notable feature is spectral asymmetry: the positron branch contains the special level $$H_{\mathrm{free}}=\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m,$$4, which has no counterpart in the electron branch (Longhi, 2010).

A complementary analytic treatment of the $$H_{\mathrm{free}}=\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m,$$5 model organizes the eigenvalues as

$$H_{\mathrm{free}}=\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m,$$6

distinguishing “conventional” and “unconventional” branches. In that framework, the absence of a conventional negative-energy state with quantum number $$H_{\mathrm{free}}=\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m,$$7 is tied to the fact that $$H_{\mathrm{free}}=\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m,$$8 is independent of $$H_{\mathrm{free}}=\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m,$$9, so the relevant branch point in the complex-$(3+1)$0 plane disappears (Cao et al., 2019).

In three dimensions, spherical symmetry leads to two spin–orbit families. A compact representation introduces

$(3+1)$1

$(3+1)$2

and the undeformed spectrum

$(3+1)$3

This makes explicit that the oscillator scale and the spin–orbit splitting are intertwined already at the exact relativistic level (Boumali et al., 25 Feb 2026).

The nonrelativistic limit does not reduce to the plain oscillator without residue. In one dimension, the Dirac oscillator behaves like a harmonic oscillator plus a spin term in a resonant field of frequency $(3+1)$4 (Cao et al., 2019). In three dimensions, the oscillator survives together with the large spin–orbit interaction (Boumali et al., 25 Feb 2026). A persistent misconception is therefore to regard the model as only a relativistic reparametrization of the Schrödinger oscillator; its relativistic branch structure and spin dependence are essential, not peripheral.

3. Algebraic, analytic, and field-theoretic viewpoints

The Dirac oscillator admits several algebraic reformulations. In $(3+1)$5 dimensions it can be written as

$(3+1)$6

with spectrum

$(3+1)$7

and in the massless limit

$(3+1)$8

This form makes transparent the coupling between a two-level degree of freedom and oscillator ladder operators (Franco-Villafañe et al., 2013).

In $(3+1)$9 and $$i\hbar \frac{\partial}{\partial t}\Psi(\mathbf{r},t)=\left[c\,\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m c^2\right]\Psi(\mathbf{r},t).$$0 dimensions the system can be mapped onto the Jaynes–Cummings and anti-Jaynes–Cummings models. For time-dependent frequency, the optical counterpart yields modified angular-momentum dynamics, spin-orbit entanglement, and noticeable changes in the Zitterbewegung. One time dependence produces aperiodic evolution of the observables, whereas another admits analytical solutions (Tsutsui et al., 21 Dec 2025). In a related supersymmetric extension, the generalized Dirac oscillator in $$i\hbar \frac{\partial}{\partial t}\Psi(\mathbf{r},t)=\left[c\,\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m c^2\right]\Psi(\mathbf{r},t).$$1 dimensions maps to an anti–Jaynes–Cummings-like Hamiltonian in which the spin operators couple with the supercharges (Ghosh et al., 22 Feb 2025).

Supersymmetric quantum mechanics also supplies a natural route to generalized one-dimensional models. Writing

$$i\hbar \frac{\partial}{\partial t}\Psi(\mathbf{r},t)=\left[c\,\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m c^2\right]\Psi(\mathbf{r},t).$$2

the upper component satisfies

$$i\hbar \frac{\partial}{\partial t}\Psi(\mathbf{r},t)=\left[c\,\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m c^2\right]\Psi(\mathbf{r},t).$$3

For the isotonic choice

$$i\hbar \frac{\partial}{\partial t}\Psi(\mathbf{r},t)=\left[c\,\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m c^2\right]\Psi(\mathbf{r},t).$$4

the effective potential becomes harmonic plus singular, and the nonrelativistic limit reproduces the isotonic oscillator rather than the ordinary harmonic oscillator (Ghosh et al., 22 Feb 2025).

At the field-theoretic level, canonical quantization is possible because the Dirac oscillator is characterized by the absence of the Klein paradox and by the completeness of its eigenfunctions. In that formulation, the field is constituted by infinite degrees of freedom identified as decoupled quantum linear harmonic oscillators, with quanta given by relativistic oscillator energies rather than free-particle energies (Quimbay et al., 2012).

4. Experimental and analog realizations

The Dirac oscillator is not only a formal model. A photonic realization was proposed in fiber Bragg gratings, where the refractive-index modulation

$$i\hbar \frac{\partial}{\partial t}\Psi(\mathbf{r},t)=\left[c\,\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m c^2\right]\Psi(\mathbf{r},t).$$5

induces coupled-mode equations that can be recast in Dirac form. With

$$i\hbar \frac{\partial}{\partial t}\Psi(\mathbf{r},t)=\left[c\,\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m c^2\right]\Psi(\mathbf{r},t).$$6

the choice

$$i\hbar \frac{\partial}{\partial t}\Psi(\mathbf{r},t)=\left[c\,\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m c^2\right]\Psi(\mathbf{r},t).$$7

implements the one-dimensional Dirac oscillator exactly in the idealized infinite-grating limit (Longhi, 2010).

In that photonic system, bound states appear as narrow transmission peaks inside the stop band. Positive-energy, electron-like states lie above the Bragg frequency, and negative-energy, positron-like states lie below it. The transmission spectrum displays the asymmetry of the Dirac-oscillator branches, including the additional negative-detuning resonance associated with the extra bound state at $$i\hbar \frac{\partial}{\partial t}\Psi(\mathbf{r},t)=\left[c\,\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m c^2\right]\Psi(\mathbf{r},t).$$8 (Longhi, 2010). For typical telecom parameters $$i\hbar \frac{\partial}{\partial t}\Psi(\mathbf{r},t)=\left[c\,\boldsymbol{\alpha}\cdot\mathbf{p}+\beta m c^2\right]\Psi(\mathbf{r},t).$$9, $$\mathbf{p}\longrightarrow \mathbf{p}-i m\omega\,\beta\,\mathbf{r},$$0, and $$\mathbf{p}\longrightarrow \mathbf{p}-i m\omega\,\beta\,\mathbf{r},$$1 nm, the characteristic scales are $$\mathbf{p}\longrightarrow \mathbf{p}-i m\omega\,\beta\,\mathbf{r},$$2 mm and $$\mathbf{p}\longrightarrow \mathbf{p}-i m\omega\,\beta\,\mathbf{r},$$3 ps, with an example grating length $$\mathbf{p}\longrightarrow \mathbf{p}-i m\omega\,\beta\,\mathbf{r},$$4 cm (Longhi, 2010).

The first experimental realization of the one-dimensional Dirac oscillator was achieved with microwaves in a chain of coupled dielectric disks. The tight-binding representation uses dimers with constant intra-dimer coupling $$\mathbf{p}\longrightarrow \mathbf{p}-i m\omega\,\beta\,\mathbf{r},$$5 and inter-dimer couplings

$$\mathbf{p}\longrightarrow \mathbf{p}-i m\omega\,\beta\,\mathbf{r},$$6

which reproduce the spectrum

$$\mathbf{p}\longrightarrow \mathbf{p}-i m\omega\,\beta\,\mathbf{r},$$7

In the massless case, the measured resonances follow

$$\mathbf{p}\longrightarrow \mathbf{p}-i m\omega\,\beta\,\mathbf{r},$$8

and in the effective massive case

$$\mathbf{p}\longrightarrow \mathbf{p}-i m\omega\,\beta\,\mathbf{r},$$9

The experiment used identical dielectric disks with height $$H_{\mathrm{DO}}=\boldsymbol{\alpha}\cdot\left(\mathbf{p}-i m\omega\,\beta\,\mathbf{r}\right)+\beta m$$0 mm, radius $$H_{\mathrm{DO}}=\boldsymbol{\alpha}\cdot\left(\mathbf{p}-i m\omega\,\beta\,\mathbf{r}\right)+\beta m$$1 mm, refractive index $$H_{\mathrm{DO}}=\boldsymbol{\alpha}\cdot\left(\mathbf{p}-i m\omega\,\beta\,\mathbf{r}\right)+\beta m$$2, and isolated TE resonance $$H_{\mathrm{DO}}=\boldsymbol{\alpha}\cdot\left(\mathbf{p}-i m\omega\,\beta\,\mathbf{r}\right)+\beta m$$3 GHz (Franco-Villafañe et al., 2013).

An important clarification concerns the mass term in that microwave realization. Because all disks are identical, an intrinsic on-site-energy asymmetry is absent; the “massive DO” is instead produced by finite-size effects through a distorted coupling pattern, and the gap vanishes as the number of sites $$H_{\mathrm{DO}}=\boldsymbol{\alpha}\cdot\left(\mathbf{p}-i m\omega\,\beta\,\mathbf{r}\right)+\beta m$$4 due to chiral symmetry (Franco-Villafañe et al., 2013). This distinguishes an effective spectral gap from a fundamental mass parameter in the underlying continuum model.

5. External fields, geometry, and deformations

The Dirac oscillator remains exactly tractable in several nontrivial backgrounds. In a uniform axial magnetic field, the full $$H_{\mathrm{DO}}=\boldsymbol{\alpha}\cdot\left(\mathbf{p}-i m\omega\,\beta\,\mathbf{r}\right)+\beta m$$5-dimensional problem separates into a two-dimensional oscillator in the plane perpendicular to the field and a one-dimensional oscillator along the field direction. The transverse levels involve the combination of oscillator and Larmor frequencies through $$H_{\mathrm{DO}}=\boldsymbol{\alpha}\cdot\left(\mathbf{p}-i m\omega\,\beta\,\mathbf{r}\right)+\beta m$$6, while the spectrum also contains spin–orbit and spin–Larmor couplings (Moniruzzaman et al., 2015).

Curved and topologically nontrivial backgrounds deform the angular structure of the bound states. In a spinning cosmic string spacetime with curvature and torsion, exact spectra depend on effective angular quantum numbers involving the oscillator frequency, the angular deficit parameter $$H_{\mathrm{DO}}=\boldsymbol{\alpha}\cdot\left(\mathbf{p}-i m\omega\,\beta\,\mathbf{r}\right)+\beta m$$7, the torsional parameters $$H_{\mathrm{DO}}=\boldsymbol{\alpha}\cdot\left(\mathbf{p}-i m\omega\,\beta\,\mathbf{r}\right)+\beta m$$8 and $$H_{\mathrm{DO}}=\boldsymbol{\alpha}\cdot\left(\mathbf{p}-i m\omega\,\beta\,\mathbf{r}\right)+\beta m$$9, and the longitudinal momentum $\hbar=c=1$0. The flat-space Moshinsky spectrum is recovered when curvature and torsion vanish (Boumali, 19 Sep 2025). In a spinning cosmic string spacetime, the oscillator produces an implicit energy equation because the rotation parameter enters the effective angular term through $\hbar=c=1$1 (Hosseinpour et al., 2019). In gravity’s-rainbow versions of the cosmic-string background, the rainbow functions can either preserve or distort the symmetry between positive and negative energy branches, depending on the chosen scenario (Bakke et al., 2018).

Quantum-group and high-energy deformations modify the spectral algebra rather than the oscillator-spinor structure alone. In the $\hbar=c=1$2-Dirac oscillator, the deformation preserves parity while breaking charge conjugation and time reversal symmetries, and the deformation parameter breaks the infinite degeneracy of the Dirac oscillator; for $\hbar=c=1$3, the standard energy eigenvalues and eigenfunctions are recovered (Andrade et al., 2013). In doubly special relativity, the three-dimensional modified Dirac oscillator retains the oscillator-spinor eigenfunctions dictated by spherical symmetry, while the relation between $\hbar=c=1$4 and the energy is deformed in branch-dependent ways; the undeformed limit $\hbar=c=1$5 or $\hbar=c=1$6 is recovered smoothly (Boumali et al., 25 Feb 2026).

Recent extensions also include non-Abelian gauge structure and non-Hermitian variants. A non-Abelian extension introduces an additional internal-space term in the effective momentum substitution,

$\hbar=c=1$7

while preserving total angular momentum $\hbar=c=1$8 (Guerah et al., 11 Apr 2025). By contrast, the inverted Dirac oscillator replaces the usual non-Hermitian momentum modification with a Hermitian one, $\hbar=c=1$9, producing a non-Hermitian Hamiltonian with an inverted potential, continuous spectrum, and eigenfunctions that fail to be square-integrable; the model is pseudo-$c$0-symmetric and is related to the usual Dirac oscillator by an unbounded, non-unitary transformation (Maamache, 13 Jun 2026).

6. Applications, methodological roles, and interpretive issues

Beyond its status as an exactly solvable model, the Dirac oscillator has become a practical computational tool. In covariant density functional theory, a Dirac-oscillator basis has been used for self-consistent nuclear structure calculations. For a selected set of doubly-magic nuclei, binding energies and ground-state densities obtained in this basis reproduced with high accuracy those derived using the Runge–Kutta method, and the results suggest a path toward generalization to systems with axial symmetry (Yang et al., 2020).

The same basis is useful in relativistic scattering theory. Dirac oscillators provide an excellent expansion basis for relativistic $c$1-matrix techniques, allowing reaction and bound-state problems to be treated within a common formalism. In relativistic impulse approximation calculations, the combination of the Dirac oscillator and the $c$2-matrix approach permits exact calculation of exchange terms and construction of scattering waves orthogonal to bound-state wave functions (Grineviciute et al., 2014).

The model also occupies a stable place in quantum optics and analog simulation. The exact Jaynes–Cummings and anti–Jaynes–Cummings correspondences in low dimensions explain why the Dirac oscillator appears naturally in trapped-ion, microwave, and photonic platforms (Tsutsui et al., 21 Dec 2025, Franco-Villafañe et al., 2013). This suggests that its role is dual: it is both a relativistic quantum system and a transfer principle linking relativistic bound-state dynamics to experimentally tunable boson–spin couplings.

Several recurrent misunderstandings can be stated precisely. The Dirac oscillator is not simply the harmonic oscillator with relativistic notation; its spectrum can be asymmetric between particle and antiparticle branches, and special low-lying states may lack naive partners (Longhi, 2010, Cao et al., 2019). It is not only a textbook model; direct microwave and photonic analog realizations exist and resolve its spectral features experimentally (Franco-Villafañe et al., 2013, Longhi, 2010). Nor is every “mass term” in an analog realization fundamental; in the microwave experiment, the observed gap can be an effective finite-size feature rather than an intrinsic onsite asymmetry (Franco-Villafañe et al., 2013).

Taken together, these developments place the Dirac oscillator at the intersection of relativistic spectral theory, supersymmetric quantum mechanics, quantum simulation, and computational many-body physics. Its defining substitution is elementary, but the resulting structure is unusually rich: exact solvability coexists with nontrivial spin–orbit dynamics, asymmetric relativistic branches, analytic continuations, experimental analogs, and controlled generalizations to magnetic, curved, deformed, and non-Abelian settings.