---
title: Dunkl Creation and Annihilation Operators
url: https://www.emergentmind.com/topics/dunkl-creation-and-annihilation-operators
type: topic
---

# Dunkl Creation and Annihilation Operators

Searching arXiv for recent and foundational papers on Dunkl creation and annihilation operators.
arXiv search query: "Dunkl oscillator creation annihilation operators"

Dunkl creation and annihilation operators are ladder operators obtained by replacing the ordinary derivative in harmonic-oscillator constructions with a Dunkl differential-difference operator tied to a reflection symmetry. In one dimension, the reflection operator acts by parity, $R f(x)=f(-x)$, with $R^2=1$, and the Dunkl derivative is typically written as $D_x=\partial_x+\frac{\mu}{x}(1-R)$, although some papers denote the deformation parameter by $\theta$. The corresponding ladder operators,
\[
a=\frac{x+D_x}{\sqrt{2}}, \qquad a^\dagger=\frac{x-D_x}{\sqrt{2}},
\]
or their coordinatewise analogues, satisfy the deformed commutator $[a,a^\dagger]=1+2\mu R$ rather than the ordinary Heisenberg relation. This inserts parity directly into the oscillator algebra and makes Dunkl ladders the basic building blocks for solvable nonrelativistic, relativistic, superintegrable, and statistical-mechanical models [1212.4459] [1312.3877] [2202.02871] [2508.11806] [2602.16743] [2604.22945].

## 1. Algebraic definition and parity deformation

The defining ingredient is the replacement of $\partial_x$ by a reflection-symmetric differential-difference operator. In the one-dimensional formulations used across the literature, the reflection operator obeys $R f(x)=f(-x)$, $R^2=1$, $R x=-xR$, and $R\partial_x=-\partial_xR$, while the Dunkl derivative is
\[
D_x=\partial_x+\frac{\mu}{x}(1-R)
\]
or, in alternative notation, $D_x=\partial_x+\frac{\theta}{x}(1-R)$. The deformation enters the canonical structure through a parity operator: one paper writes a Yang–Wigner type relation $[\hat x,\hat p]=i(1+\mu \hat R)$, with $\hat p=\frac{1}{i}\hat D$ and $\hat R f(x)=f(-x)$ [2202.02871]. In the oscillator realization, the resulting ladder algebra is
\[
[a,a^\dagger]=1+2\mu R,
\]
and in multidimensional settings this generalizes componentwise to $[A_i^-,A_i^+]=1+2\mu_iR_i$ or $[\phi_i,\bar\phi_i]=\delta_{ij}(1+2\theta R_i)$ [1212.4459] [2508.11806].

A characteristic feature of this algebra is that the reflection operator anticommutes with the ladders. In the standard one-dimensional realization, $Ra=-aR$ and $Ra^\dagger=-a^\dagger R$; in planar and higher-dimensional models this becomes $\{A_i^\pm,R_i\}=0$ together with commutation across distinct coordinates. The Hilbert space therefore decomposes into even and odd sectors, and the deformed commutator reduces sectorwise to a scalar, $1+2\mu s$ with $s=\pm1$. This parity resolution is not ancillary: it is the mechanism by which Dunkl operators alter spectra, symmetry algebras, and thermodynamic quantities [1212.4459] [2202.02871].

The natural inner product also changes. For the one-dimensional and planar oscillator models, the relevant weight is $|x|^{2\mu}$ in one dimension and the product weight $\prod_i |x_i|^{2\mu_i}$ in several variables; with respect to this measure the Dunkl derivatives are anti-Hermitian and the corresponding Hamiltonians are self-adjoint for $\mu_i>-1/2$ [1212.4459] [1312.3877]. In the recent anharmonic construction, the same weighted inner product is used to ensure that $a_\mu^\dagger$ is the Hermitian adjoint of $a_\mu$ [2604.22945].

## 2. Fock structure, Dunkl numbers, and parity sectors

The deformed ladder algebra admits a Fock-like basis, but the ladder coefficients are parity sensitive. A recurring device is the Dunkl number
\[
[n]_\mu = n+\mu\bigl(1-(-1)^n\bigr),
\]
which equals $2m$ on even states and $2m+1+2\mu$ on odd states. In the one-dimensional Dunkl-Fock basis,
\[
a_\mu|n\rangle=\sqrt{[n]_\mu}\,|n-1\rangle, \qquad
a_\mu^\dagger|n\rangle=\sqrt{[n+1]_\mu}\,|n+1\rangle,
\]
and $R|n\rangle=(-1)^n|n\rangle$ [2604.22945]. The same $\mu$-number structure appears in the planar oscillator, where the Cartesian basis carries the action
\[
A_x^\dagger |n_x,n_y\rangle=[n_x+1]_{\mu_x}|n_x+1,n_y\rangle,\qquad
A_x|n_x,n_y\rangle=[n_x]_{\mu_x}|n_x-1,n_y\rangle,
\]
with analogous formulas in the $y$ variable; in normalized states the coefficients are the square roots of the $\mu$-numbers [1212.4459].

Because the ladders anticommute with reflection, they flip parity. On even states the lowering coefficient is undeformed, while on odd states the coefficient is shifted by the deformation parameter. In the recent Kerr-type model this is made explicit as
\[
a_\mu|2m\rangle=\sqrt{2m}\,|2m-1\rangle,\qquad
a_\mu|2m+1\rangle=\sqrt{2m+1+2\mu}\,|2m\rangle,
\]
with analogous formulas for $a_\mu^\dagger$ [2604.22945]. This parity flip is the elementary local action from which higher quadratic symmetries are assembled.

A related formulation appears in the Dunkl ideal Fermi-gas paper, which introduces operators $\phi_i$ and $\bar\phi_i$ satisfying $[\phi_i,\bar\phi_i]=\delta_{ij}(1+2\theta R_i)$ and a modified number operator $\mathcal N_i=\bar\phi_i\phi_i$. Its eigenvalues are parity dependent:
\[
\mathcal N_i|n_i\rangle=\bigl[n_i+\theta(1-(-1)^{n_i})\bigr]|n_i\rangle,
\]
so $\mathcal N_i|n_i\rangle=n_i|n_i\rangle$ for even $n_i$ and $(n_i+2\theta)|n_i\rangle$ for odd $n_i$ [2508.11806]. In this sense, Dunkl ladder operators do not merely deform transition amplitudes; they also alter the spectral meaning of occupation-number operators.

## 3. Oscillator Hamiltonians and spectral consequences

For the isotropic Dunkl oscillator, the ladder operators reconstruct the Hamiltonian exactly. In the plane,
\[
H=H_x+H_y,\qquad H_i=\frac12\{A_i,A_i^\dagger\},
\]
with $A_x=(x+D_x)/\sqrt2$, $A_x^\dagger=(x-D_x)/\sqrt2$ and the analogous $y$-operators. The Cartesian eigenstates are products of one-dimensional Dunkl oscillator wavefunctions, and the energy spectrum is
\[
E=n_x+n_y+\mu_x+\mu_y+1.
\]
For fixed $N=n_x+n_y$, the level $E_N=N+\mu_x+\mu_y+1$ has $(N+1)$-fold degeneracy, as in the ordinary isotropic oscillator, but the symmetry algebra is deformed by reflections rather than being the undeformed Schwinger $u(2)$ algebra [1212.4459]. In three dimensions the same construction yields
\[
\mathcal H=\sum_{i=1}^3\left(-\frac12\mathcal D_i^2+\frac12 x_i^2\right),
\qquad
\mathcal E_N=N+\mu_1+\mu_2+\mu_3+\frac32,
\]
with degeneracy $\frac{(N+1)(N+2)}{2}$ [1312.3877].

The ladder framework also governs separated wavefunctions. In Cartesian coordinates, generalized Hermite polynomials appear; in polar, cylindrical, and spherical coordinates the solutions involve Laguerre and Jacobi families, with parity labels entering the separation constants and admissible quantum numbers [1212.4459] [1312.3877]. This is not a secondary analytic artifact: it reflects the fact that the reflection operators commute with the Hamiltonians while anticommute with the basic ladders.

Relativistic Dunkl oscillators sharpen the same mechanism. In the one-dimensional Dunkl–Klein–Gordon and Dunkl–Dirac systems, the Hamiltonians commute with reflection, so eigenfunctions can be chosen with $\hat R\psi^s=s\psi^s$, $s=\pm$. The paper introduces Dunkl creation and annihilation operators linear in $m\omega\hat x$ and $\hbar \hat D$, with commutator
\[
[\hat a_D^\dagger,\hat a_D]=1+2\mu \hat R.
\]
The resulting relativistic spectrum for both the Dunkl–Klein–Gordon and Dunkl–Dirac oscillators is
\[
\frac{E_n^s}{m}=\pm \sqrt{4nr+2r\left(\mu+\frac12\right)(1-s)+1},
\qquad r\equiv \frac{\omega}{m},
\]
which depends explicitly on parity and reduces to the usual relativistic oscillator in the even sector when $\mu\to0$ [2202.02871]. The parity split is therefore intrinsic to the ladder algebra itself.

## 4. Schwinger constructions and reflection-deformed symmetry algebras

In two dimensions, combining two independent Dunkl modes via a Schwinger construction produces the symmetry algebra of the planar oscillator. With
\[
J_+=A_x^\dagger A_y,\qquad J_-=A_y^\dagger A_x,\qquad J_0=\frac{H_x-H_y}{2},
\]
or equivalently $J_1=(J_++J_-)/2$, $J_2=(J_+-J_-)/(2i)$, $J_3=J_0$, one obtains symmetry generators commuting with the full Hamiltonian. The operator
\[
J_2=\frac{i}{2}(xD_y-yD_x)
\]
is the Dunkl angular momentum and controls polar separation. Together with the involutions $R_x$, $R_y$ and the central Hamiltonian, these operators generate the Schwinger–Dunkl algebra $sd(2)$, an extension of $u(2)$ with reflections [1212.4459].

This algebra organizes both the degeneracies and the interbasis transforms. In the planar model, Cartesian–polar overlap coefficients are expressed as linear combinations of dual $-1$ Hahn polynomials, and the same structure is identified with the Clebsch–Gordan problem of $sl_{-1}(2)$ [1212.4459]. The follow-up representation-theoretic analysis constructs Cartesian and circular bases from parabosonic creation and annihilation operators. In the Cartesian basis, $J_3$ is diagonal and $J_2$ acts tridiagonally; in the circular basis, $J_2$ becomes block upper-triangular with $2\times2$ blocks, while $J_3$ is tridiagonal. The Cartesian–circular expansion coefficients are Krawtchouk polynomials, and in the general anisotropic case the eigenvectors of $J_2$ are generated by Heun polynomials with components expressed in para-Krawtchouk, ანუ complementary Bannai–Ito, polynomials; in the fully isotropic case they reduce to little $-1$ Jacobi or ordinary Jacobi polynomials [1302.6142].

In three dimensions, the same paradigm produces a reflection-deformed extension of $u(3)$. The one-dimensional Dunkl ladders $A_\pm^{(i)}=(x_i\mp \mathcal D_i)/\sqrt2$ realize $sl_{-1}(2)$ in each coordinate, and bilinears built from them generate the Schwinger–Dunkl algebra $sd(3)$. The corresponding constants of motion commute with the Hamiltonian, the system is maximally superintegrable, and separation is available in Cartesian, cylindrical, and spherical coordinates [1312.3877]. A plausible implication is that Dunkl creation and annihilation operators play in these models the same structural role that bosonic ladders play in ordinary oscillator-based representation theory, but with involutions replacing continuous rotational symmetry as the primary deformation data.

## 5. Quadratic closures and \(su(1,1)\) realizations

A second major theme is that quadratic combinations of Dunkl ladders often generate undeformed $su(1,1)$ algebras even though the underlying first-order commutator is reflection deformed. In the singular Dunkl oscillator, one-dimensional Hamiltonians
\[
H_i=-\frac12(D_{x_i}^{\mu_i})^2+\frac{x_i^2}{2}+\frac{\alpha_i+\beta_iR_i}{2x_i^2}
\]
admit quadratic ladder operators
\[
A_i^\dagger=(a_i^\dagger)^2-\frac{\alpha_i+\beta_iR_i}{2x_i^2},\qquad
A_i=a_i^2-\frac{\alpha_i+\beta_iR_i}{2x_i^2},
\]
satisfying
\[
[H_i,A_i^\dagger]=2A_i^\dagger,\qquad [H_i,A_i]=-2A_i,\qquad [A_i,A_i^\dagger]=-4H_i.
\]
These give an $su(1,1)$ dynamical algebra in each coordinate. In two dimensions, the resulting symmetry generators form a cubic algebra in the ladder basis and, when rewritten in terms of separation symmetries, a quadratic Hahn algebra, identified as a special case of Askey–Wilson $AW(3)$ with central involutions [1305.2126].

The recent parametric-amplifier generalization makes the quadratic closure especially explicit. Defining
\[
K_+^\mu=\frac12(a_\mu^\dagger)^2,\qquad
K_-^\mu=\frac12 a_\mu^2,\qquad
K_0^\mu=\frac14(a_\mu^\dagger a_\mu+a_\mu a_\mu^\dagger),
\]
the commutation relations become
\[
[K_0^\mu,K_\pm^\mu]=\pm K_\pm^\mu,\qquad [K_-^\mu,K_+^\mu]=2K_0^\mu.
\]
The Casimir acts as
\[
C^\mu|n\rangle=\left(\frac{\mu^2}{4}-\frac{\mu}{4}R-\frac{3}{16}\right)|n\rangle=k(k-1)|n\rangle,
\]
with parity-dependent Bargmann indices
\[
k_+=\frac14+\frac{\mu}{2},\qquad k_-=\frac34+\frac{\mu}{2}.
\]
The Dunkl parametric-amplifier Hamiltonian
\[
H_\mu=\omega a_\mu^\dagger a_\mu+f a_\mu^2+f^*(a_\mu^\dagger)^2
\]
is diagonalized both by an $su(1,1)$ tilting transformation and by a generalized Bogoliubov transformation, with stability condition $\omega>2|f|$ and exact spectrum
\[
E_n=\frac12\Omega_\mu\left(n+\mu+\frac12\right),\qquad
\Omega_\mu=2\sqrt{\omega^2-4|f|^2}.
\]
For squeezed vacuum, the Mandel parameter is independent of $\mu$, whereas $g^{(2)}(0)$ depends explicitly on $\mu$ and modifies photon bunching [2602.16743].

The same $su(1,1)$ mechanism appears in the Dunkl anharmonic, or Kerr, oscillator. There the Hamiltonian
\[
H_\mu=\omega a_\mu^\dagger a_\mu+\frac{\lambda}{2}(a_\mu^\dagger)^2a_\mu^2
\]
is rewritten in terms of $K_0^\mu$ and $K_+^\mu K_-^\mu$, yielding parity-dependent energies
\[
E_{2m}^\mu=2\lambda m^2+m(2\omega+2\lambda\mu-\lambda),\qquad
E_{2m+1}^\mu=2\lambda m^2+m(2\omega+2\lambda\mu+\lambda)+\omega(1+2\mu).
\]
The quantum dynamics exhibits collapse-and-revival phenomena with fundamental revival period $t_R=2\pi/\lambda$ independent of $\mu$, while $\mu=\tfrac12$ produces a perfect half-period revival at $t=\pi/\lambda$; the quadrature variance also shows interference-induced squeezing around $t\approx\pi$ for specific deformation values [2604.22945]. This suggests that the parity deformation survives in the spectrum and dynamics even when the quadratic algebra itself is undeformed.

## 6. Relativistic, thermodynamic, and many-body extensions

Dunkl creation and annihilation operators also support thermodynamic and open-system interpretations. In the relativistic one-dimensional Dunkl–Dirac oscillator, the ladder operators allow the Hamiltonian to be rewritten as
\[
\hat{\mathcal H}_D=g(\sigma^-\hat a_D+\sigma^+\hat a_D^\dagger)+m\sigma_z,
\]
which is identified as a Dunkl Anti–Jaynes–Cummings model. The same work computes parity-resolved thermal quantities from the spectrum,
\[
Z^s=\sum_{n=0} e^{-(E_n^s-E_0^s)/(K_BT)},
\]
and shows that in the high-temperature limit the parity dependence washes out, with $Z^+=Z^-=\tau^2/(2r)$ in the notation of that paper. The model is therefore presented as “an appropriate scenario for the theory of an open quantum system coupled to a thermal bath” [2202.02871].

A many-body extension is developed for the ideal Fermi gas in Dunkl kinematics. There, the second-quantized fermionic operators still satisfy the canonical anticommutation relations,
\[
\{c_i,c_j^\dagger\}=\delta_{ij},\qquad \{c_i,c_j\}=0=\{c_i^\dagger,c_j^\dagger\},
\]
so Fermi–Dirac statistics are preserved, but the single-particle operator content is deformed by reflections. The grand partition function factorizes as
\[
\mathcal Z_{FD}=\prod_i\left[1+e^{-\beta(\epsilon_i-\mu)(1+2\theta)}\right],
\]
equivalently corresponding to a Hamiltonian $H_D=\sum_i \epsilon_i(1+2\theta)c_i^\dagger c_i$. The average occupation becomes
\[
f(\epsilon_i)=\frac{1+2\theta}{e^{\beta(\epsilon_i-\mu)(1+2\theta)}+1},
\]
and the paper imposes $-1/2\le \theta\le 0$ to keep $0\le f\le 1$ for all energies [2508.11806]. In this formulation, Dunkl creation and annihilation operators do not alter the CAR; they modify the effective spectral weights carried by occupied modes.

Across these settings, the vanishing-deformation limit restores the ordinary derivative, the standard commutator $[a,a^\dagger]=1$, the undeformed number operator, and the usual oscillator or many-body thermodynamics. The planar and three-dimensional Dunkl oscillators reduce to the ordinary harmonic oscillator with $u(2)$ and $u(3)$ symmetry, while the parametric-amplifier and Kerr constructions recover the standard squeezed-state and Kerr-medium formulas when $\mu\to0$ [1212.4459] [1312.3877] [2602.16743] [2604.22945]. This continuity is central: Dunkl creation and annihilation operators are best understood not as an unrelated ladder formalism, but as a parity-deformed extension of conventional oscillator methods whose distinctive content is encoded in the reflection operator.

Source: https://www.emergentmind.com/topics/dunkl-creation-and-annihilation-operators