---
title: Generalized Dirac Oscillator
url: https://www.emergentmind.com/topics/generalized-dirac-oscillator
type: topic
---

# Generalized Dirac Oscillator

Searching arXiv for recent and foundational papers on the generalized Dirac oscillator.
The generalized Dirac oscillator is a relativistic oscillator model obtained by replacing the linear non-minimal coupling of the standard Dirac oscillator with a more general interaction function. In \((1+1)\) dimensions this is commonly written as the replacement \(m\omega x\to W(x)\), leading to
\[
\bigl(c\sigma_x(p_x-iW(x)\sigma_z)+mc^2\sigma_z\bigr)\psi=E\psi,
\]
or equivalently, in the notation \(f(x)\),
\[
H_{\text{GDO}}=\sigma_x p+\sigma_y f(x)+\sigma_z m
\]
with \(c=1\) in the latter convention. In this sense, the standard Dirac oscillator is the special case \(W(x)=m\omega x\) or \(f(x)=m\omega x\), whereas the generalized model admits arbitrary local couplings, complex interactions, nonlocal kernels, curved-space versions, and deformed kinematics, while still often retaining an exactly solvable supersymmetric structure [2607.04348] [2603.03572] [1301.2035].

## 1. Defining equations and canonical formulations

The standard starting point is the Moshinsky–Szczepaniak substitution
\[
\mathbf p\to \mathbf p-i m\omega \beta \mathbf r,
\]
which preserves linearity in momentum and coordinate. The generalized Dirac oscillator replaces the linear term by a function. In \((1+1)\) dimensions this appears as
\[
p \longrightarrow p-i\,\beta\, f(x),
\]
or, in the equivalent \(W(x)\) notation,
\[
m\omega x \to W(x).
\]
For real \(f(x)\), the Hamiltonian
\[
H_{GDO}=c\,\sigma_x\left(p_x-i\sigma_z f(x)\right)+\beta mc^2
\]
is Hermitian; for complex \(f(x)\), the same formal structure supports pseudo-Hermitian and \(\mathcal{PT}\)-symmetric realizations [1301.2035].

The generalized formulation is not confined to one dimension. In \((2+1)\) dimensions one finds radial substitutions such as
\[
p_x \to p_x-i\beta\,x\,f(r),\qquad p_y \to p_y-i\beta\,y\,f(r),
\]
and in curved backgrounds the ordinary replacement
\[
p_\mu\to p_\mu+m\omega\beta x_\mu
\]
is extended to
\[
p_\mu \to p_\mu + m\omega\,\beta\,f_\mu(x_\mu),
\]
with the cosmic-string literature often activating only the radial component \(f_\rho(\rho)\) [1811.03916]. In the recent three-dimensional DSR treatment, the undeformed spatial operator remains the ordinary Dirac oscillator, while the deformation acts on the energy reconstruction rather than on the oscillator eigenfunctions themselves [2603.15632].

## 2. Supersymmetric factorization and solvable sectors

A central structural fact is that the generalized Dirac oscillator decouples into supersymmetric partner equations. In the local \((1+1)\)-dimensional case, one introduces first-order operators
\[
\hat B=\hat p-iW(x),\qquad \hat B^+=\hat p+iW(x),
\]
or equivalently
\[
\mathcal A=p-if(x),\qquad \mathcal A^\dagger=p+i f^*(x),
\]
and obtains
\[
H_-=\hat B^+\hat B,\qquad H_+=\hat B\hat B^+,
\]
with partner potentials
\[
V_\pm(x)=f^2(x)\pm \hbar f'(x).
\]
The two spinor components satisfy Schrödinger-type equations with the same positive-\(\epsilon\) spectrum, differing only possibly in the zero mode. In exact SUSY, the zero mode satisfies
\[
\hat B\phi_0=0,\qquad \phi_0=Ce^{-\int dx\,W(x)},
\]
and the corresponding Dirac state has \(\chi_0=0\), so the ground state is separable between spin and spatial degree of freedom [2607.04348].

This factorization is the basis of exact solvability. The DSR–Morse analysis recalls the standard shape-invariance condition
\[
V_+(x;a_0)=V_-(x;a_1)+R(a_0),\qquad a_1=f(a_0),
\]
from which
\[
\epsilon_n=\sum_{j=0}^{n-1}R(a_j)
\]
follows algebraically. For the pseudo-Hermitian complexified Morse interaction
\[
f(x)=D-(A+iB)e^{-\alpha x},
\]
the partner potentials are shape invariant and the spatial spectrum is
\[
\epsilon_n=D^2-(D-n\hbar\alpha)^2,\qquad n<\frac{D}{\hbar\alpha},
\]
so the Morse sector has only finitely many bound states [2603.03572].

A distinct exactly solvable direction is the isotonic generalization. In \((1+1)\) dimensions, choosing
\[
W(x)=ax+\frac{b}{x},\qquad x>0,
\]
yields an effective isotonic potential, and the exact relativistic spectrum is
\[
E^2=(mc^2)^2\left[1+\frac{a}{m^2c^2}\Big(4n+2b+1+\sqrt{1+4b(b+1)}\Big)\right].
\]
In the non-relativistic limit this reduces to an isotonic oscillator Hamiltonian on the half-line. In \((2+1)\) dimensions, the analogous choice
\[
\mathbf W(\mathbf r)=a\mathbf r+\frac{b}{r^2}\mathbf r
\]
produces a radial isotonic equation and an anti-Jaynes–Cummings-like Hamiltonian in which the spin operators couple with the supercharges [2502.16165].

## 3. Pseudo-Hermiticity, \(\mathcal{PT}\) symmetry, and complex interactions

Generalized Dirac oscillators with complex interactions form a major branch of the subject. The basic criterion is \(\eta\)-pseudo-Hermiticity,
\[
H^\dagger=\eta H\eta^{-1},
\]
with the translation-type metric
\[
\eta=e^{-\theta p_x}.
\]
Because \(\eta\) implements an imaginary shift,
\[
\eta x\eta^{-1}=x+i\hbar\theta,
\]
the condition
\[
f(x+i\hbar\theta)=f^*(x)
\]
guarantees pseudo-Hermiticity of the local GDO. For such systems one may define \(\rho=\sqrt{\eta}\) and obtain a Hermitian counterpart \(h=\rho H\rho^{-1}\). The 2013 analysis gave explicit complex Morse-type and periodic Rosen–Morse-type interactions satisfying this shifted-conjugation condition and showed that the non-Hermitian GDO can be isospectral to a Hermitian Dirac Hamiltonian with a real interaction function [1301.2035].

The 2026 DSR review extends this framework by distinguishing the ordinary adjoint from the metric adjoint
\[
\mathcal A^\#=\eta^{-1}\mathcal A^\dagger \eta,
\]
so that the supersymmetric partner Hamiltonians become
\[
H_-=\mathcal A^\#\mathcal A,\qquad H_+=\mathcal A\mathcal A^\#.
\]
The same paper also formulates the \(\mathcal{PT}\)-symmetry criterion
\[
V(x)=V^*(-x)
\]
for the local partner Schrödinger operators and emphasizes that pseudo-Hermiticity or \(\mathcal{PT}\) symmetry secures a real spatial spectrum \(\{\epsilon_n\}\), while DSR modifies only the final map \(\epsilon_n\mapsto E_n\) [2603.03572].

These complex extensions also admit model correspondences. The generalized Dirac oscillator is exactly identified with a generalized anti-Jaynes–Cummings Hamiltonian
\[
H_{GAJC}=\Omega(\sigma_+A^\# + \sigma_-A)+\delta\sigma_z
\]
under the identification \(\Omega=c\), \(\delta=mc^2\). Under the sign reversal \(f(x)\to -f(x)\), which interchanges \(A\leftrightarrow A^\#\), the model becomes a generalized Jaynes–Cummings system. This places the GDO inside a broader algebraic family of spin-boson-like models [1301.2035].

## 4. External fields, position-dependent mass, and curved geometry

One widely studied extension couples the generalized Dirac oscillator to a scalar electric potential. In \((1+1)\) dimensions,
\[
\left(\sigma_x(p_x-iW(x)\sigma_z)+m\sigma_z\right)\psi=(E-U(x))\psi,
\]
and for the special choice
\[
U(x)=\kappa W(x)
\]
the second-order equation reduces to a SUSY form with the energy-dependent superpotential
\[
\widetilde W(x)=\sqrt{1-\kappa^2}\left(W(x)+\frac{\kappa E}{1-\kappa^2}\right).
\]
This yields exact solutions whenever \(\widetilde W(x)\) is shape invariant. The same analysis shows that sufficiently strong electric fields destroy the bounded eigenstates: bound states require
\[
|\kappa|<1,
\]
whereas for \(|\kappa|>1\) the effective confining term changes sign and discrete states disappear [1804.06091].

A related extension introduces a position-dependent mass. The \((1+1)\)-dimensional equation
\[
\big[\sigma_x(p_x-i\sigma_z f(x))+\sigma_z m(x)+V(x)-E\big]\psi(x)=0
\]
becomes exactly reducible when
\[
f(x)=\kappa_f W(x),\qquad m(x)=\kappa_m W(x),\qquad V(x)=\kappa_v W(x).
\]
The resulting scalar equation has supersymmetric form with
\[
\tilde W(x)= \sqrt{\kappa_f^2+\kappa_m^2-\kappa_v^2}\,W(x) +\frac{\kappa_v E}{\sqrt{\kappa_f^2+\kappa_m^2-\kappa_v^2}},
\]
and bound states exist only below the critical electric-field strength
\[
|\kappa_v|<\sqrt{\kappa_f^2+\kappa_m^2}.
\]
The same paper also constructs exact zero-energy states, including a step-profile solution that reduces to the Jackiw–Rebbi mode when the oscillator term vanishes [1808.03962].

Curved-space realizations are likewise standard. In cosmic-string space-time,
\[
ds^2=-dt^2+d\rho^2+\alpha^2\rho^2 d\varphi^2+dz^2,
\]
the conical defect enters through the effective angular parameter
\[
\lambda=\left(\frac{l+\frac12}{\alpha}\mp\frac12\right),
\]
and radial couplings such as the Cornell choice
\[
f_\rho(\rho)=a\rho-\frac{b}{\rho}
\]
lead to confluent-hypergeometric solutions, whereas Yukawa-, Hulthén-, generalized-Morse-, and singular-type couplings lead to hypergeometric or power-series solutions [1811.03916]. In \((1+2)\)-dimensional cosmic-string space-time with Aharonov–Casher coupling, the Coulomb-type generalized oscillator
\[
f(r)=\frac{N_1}{r}
\]
produces exact radial Laguerre solutions, and the relativistic energy levels depend explicitly on the Coulomb strength \(N_1\), the AC phase, the AC frequency, and the deficit parameter \(\alpha\); the degeneracy is broken by the Coulomb-strength parameter under the combined influence of curvature and the Aharonov–Casher effect [2012.08764].

## 5. Nonlocal kernels and deformed relativistic kinematics

A recent development replaces the local multiplicative interaction by a nonlocal integral operator
\[
(\hat F\phi)(x)=\int_{\mathbb R}dx'\,f(x,x')\phi(x').
\]
The resulting nonlocal generalized Dirac oscillator preserves factorization with
\[
A=p_x-i\hat F,\qquad A^\#=p_x+i\hat F,
\]
and the two spinor components satisfy nonlocal Schrödinger-type equations with partner kernels
\[
V_{1,2}(x,x')=(f\star f)(x,x')\mp \hbar(\partial_x+\partial_{x'})f(x,x'),
\]
where
\[
(f\star f)(x,x')=\int dy\,f(x,y)f(y,x').
\]
The same work extends pseudo-Hermiticity to the kernel level through
\[
f(x+i\hbar\theta,x'+i\hbar\theta)=f^*(x',x),
\]
and adapts current-based localization to obtain energy-dependent equivalent local potentials together with multiplicative Perey damping factors. The localization breaks down at current zeros, which diagnose spurious solutions of the nonlocal Schrödinger problem [2603.06717].

Another major direction is DSR. In one dimension, the spatial spectral problem is first solved as
\[
H_\pm\psi=\epsilon\psi,
\]
and only afterward is the relativistic energy reconstructed. In the undeformed theory,
\[
\epsilon=E^2-m^2,
\]
but in the Magueijo–Smolin prescription the map becomes
\[
E^2-m^2\left(1-\frac{E}{k}\right)^2=\epsilon_n,
\]
whereas in the Amelino–Camelia prescription
\[
\frac{E^2-m^2}{\left(1+\frac{E}{2k}\right)^2}=\epsilon_n.
\]
The AC relation imposes the admissibility condition
\[
\epsilon_n<4k^2,
\]
so DSR can remove otherwise acceptable bound states solely through the deformed energy map; in the massless limit, the MS deformation collapses to the undeformed relation while the AC deformation remains nontrivial [2603.03572].

The same logic extends to the three-dimensional modified Dirac oscillator. There the bound-state eigenfunctions retain the oscillator-spinor structure dictated by spherical symmetry, while DSR deforms the algebraic relation between \((N,j,\ell)\) and the relativistic energy. In the generalized first-order Planck-length expansion, the master equation becomes
\[
\frac{E^{2}}{c^{2}}-\Lambda_{N j}^{(\pm)}-2l_p\alpha_2\frac{E^{3}}{c^{3}}+2l_p(\alpha_3-\alpha_1)\frac{E}{c}\Lambda_{N j}^{(\pm)} = m^{2}c^{2},
\]
with
\[
\Lambda_{N j}^{(-)} = m\hbar\omega\left[2(N-j)+1\right],\qquad
\Lambda_{N j}^{(+)} = m\hbar\omega\left[2(N+j)+3\right].
\]
This shows directly that the deformation signal increases with excitation through the oscillator scale and the spin–orbit splitting [2603.15632].

## 6. Entanglement, equivalent models, and conceptual significance

The most recent conceptual extension treats the generalized Dirac oscillator as a bipartite relativistic system with Hilbert-space factorization
\[
\mathcal H=\mathcal H_{\rm spin}\otimes \mathcal H_{\rm cont}.
\]
Because the upper and lower spinor components are SUSY partners, the spin amplitudes are fixed algebraically by the relativistic energy:
\[
\gamma=\sqrt{\frac{E/c-mc}{E/c+mc}},\qquad
a=\frac{1}{\sqrt{1+\gamma^2}},\qquad
b=\frac{\gamma}{\sqrt{1+\gamma^2}}.
\]
For excited states,
\[
|\psi\rangle=a\left(|\uparrow\rangle\tilde\phi(x)+\gamma|\downarrow\rangle\tilde\chi(x)\right),
\]
and the reduced spin density matrix has eigenvalues
\[
\lambda_{1,2}=\frac12\left(1\pm\sqrt{\frac{m^2c^4}{E^2}+\left(1-\frac{m^2c^4}{E^2}\right)|Z|^2}\right),
\]
where
\[
Z=\int dx\,\tilde\phi^*(x)\tilde\chi(x).
\]
The von Neumann entropy
\[
S=-\lambda_1\log_2\lambda_1-\lambda_2\log_2\lambda_2
\]
therefore depends only on the relativistic energy \(E\) and the SUSY-partner overlap \(Z\). For odd superpotentials, \(W(-x)=-W(x)\), one has \(Z=0\), the reduced density matrix is diagonal, the entanglement vanishes in the nonrelativistic limit, and it reaches the maximal qubit value \(S=1\) as \(E\to\infty\) [2607.04348].

This entanglement result sharpens a broader pattern in the GDO literature. Hidden supersymmetric structure is not only a spectral device; it also controls reduced states, model equivalences, and effective descriptions. The generalized Dirac oscillator has been mapped to generalized anti-Jaynes–Cummings and Jaynes–Cummings Hamiltonians, to isotonic and Morse families, to magnetic-field problems in two dimensions through the exact identity
\[
\sigma_x\left[p_x-i\sigma_z f(x)\right]+\sigma_y p_y
=
\sigma_xp_x+\sigma_y\left[p_y-f(x)\right],
\]
and, in the nonlocal case, to equivalent local problems with energy-dependent potentials and damping factors [1301.2035] [2102.13561] [2603.06717].

Taken together, these developments define the generalized Dirac oscillator less as a single potential than as a framework. Its unifying feature is the replacement of the linear Dirac-oscillator coupling by a structured interaction—local or nonlocal, real or complex, flat-space or curved-space, undeformed or DSR-deformed—while preserving enough first-order factorization to make relativistic spectra, spinor structure, and, in some cases, quantum-information observables accessible in closed form.

Source: https://www.emergentmind.com/topics/generalized-dirac-oscillator