Dunkl Creation and Annihilation Operators

• Dunkl creation and annihilation operators are parity-deformed ladder operators that replace standard derivatives with Dunkl differential-difference operators incorporating reflection symmetry.
• They yield modified commutation relations and Fock structures which alter energy spectra, symmetry algebras, and thermodynamic properties across nonrelativistic and relativistic oscillator models.
• Their algebra forms the basis for advanced constructions like Schwinger–Dunkl algebras and su(1,1) realizations, providing insights into superintegrability, quantum dynamics, and many-body extensions.

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 (Genest et al., 2012, Genest et al., 2013, Hamil et al., 2022, Zenkhri et al., 15 Aug 2025, Guillén et al., 18 Feb 2026, Ojeda-Guillén et al., 24 Apr 2026).

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^2=1$0, while the Dunkl derivative is

$R^2=1$1

or, in alternative notation, $R^2=1$2. The deformation enters the canonical structure through a parity operator: one paper writes a Yang–Wigner type relation $R^2=1$3, with $R^2=1$4 and $R^2=1$5 (Hamil et al., 2022). In the oscillator realization, the resulting ladder algebra is

$R^2=1$6

and in multidimensional settings this generalizes componentwise to $R^2=1$7 or $R^2=1$8 (Genest et al., 2012, Zenkhri et al., 15 Aug 2025).

A characteristic feature of this algebra is that the reflection operator anticommutes with the ladders. In the standard one-dimensional realization, $R^2=1$9 and $D_x=\partial_x+\frac{\mu}{x}(1-R)$0; in planar and higher-dimensional models this becomes $D_x=\partial_x+\frac{\mu}{x}(1-R)$1 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, $D_x=\partial_x+\frac{\mu}{x}(1-R)$2 with $D_x=\partial_x+\frac{\mu}{x}(1-R)$3. This parity resolution is not ancillary: it is the mechanism by which Dunkl operators alter spectra, symmetry algebras, and thermodynamic quantities (Genest et al., 2012, Hamil et al., 2022).

The natural inner product also changes. For the one-dimensional and planar oscillator models, the relevant weight is $D_x=\partial_x+\frac{\mu}{x}(1-R)$4 in one dimension and the product weight $D_x=\partial_x+\frac{\mu}{x}(1-R)$5 in several variables; with respect to this measure the Dunkl derivatives are anti-Hermitian and the corresponding Hamiltonians are self-adjoint for $D_x=\partial_x+\frac{\mu}{x}(1-R)$6 (Genest et al., 2012, Genest et al., 2013). In the recent anharmonic construction, the same weighted inner product is used to ensure that $D_x=\partial_x+\frac{\mu}{x}(1-R)$7 is the Hermitian adjoint of $D_x=\partial_x+\frac{\mu}{x}(1-R)$8 (Ojeda-Guillén et al., 24 Apr 2026).

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

$D_x=\partial_x+\frac{\mu}{x}(1-R)$9

which equals $\theta$0 on even states and $\theta$1 on odd states. In the one-dimensional Dunkl-Fock basis,

$\theta$2

and $\theta$3 (Ojeda-Guillén et al., 24 Apr 2026). The same $\theta$4-number structure appears in the planar oscillator, where the Cartesian basis carries the action

$\theta$5

with analogous formulas in the $\theta$6 variable; in normalized states the coefficients are the square roots of the $\theta$7-numbers (Genest et al., 2012).

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

$\theta$8

with analogous formulas for $\theta$9 (Ojeda-Guillén et al., 24 Apr 2026). 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 $$a=\frac{x+D_x}{\sqrt{2}}, \qquad a^\dagger=\frac{x-D_x}{\sqrt{2}},$$0 and $$a=\frac{x+D_x}{\sqrt{2}}, \qquad a^\dagger=\frac{x-D_x}{\sqrt{2}},$$1 satisfying $$a=\frac{x+D_x}{\sqrt{2}}, \qquad a^\dagger=\frac{x-D_x}{\sqrt{2}},$$2 and a modified number operator $$a=\frac{x+D_x}{\sqrt{2}}, \qquad a^\dagger=\frac{x-D_x}{\sqrt{2}},$$3. Its eigenvalues are parity dependent: $$a=\frac{x+D_x}{\sqrt{2}}, \qquad a^\dagger=\frac{x-D_x}{\sqrt{2}},$$4 so $$a=\frac{x+D_x}{\sqrt{2}}, \qquad a^\dagger=\frac{x-D_x}{\sqrt{2}},$$5 for even $$a=\frac{x+D_x}{\sqrt{2}}, \qquad a^\dagger=\frac{x-D_x}{\sqrt{2}},$$6 and $$a=\frac{x+D_x}{\sqrt{2}}, \qquad a^\dagger=\frac{x-D_x}{\sqrt{2}},$$7 for odd $$a=\frac{x+D_x}{\sqrt{2}}, \qquad a^\dagger=\frac{x-D_x}{\sqrt{2}},$$8 (Zenkhri et al., 15 Aug 2025). 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,

$$a=\frac{x+D_x}{\sqrt{2}}, \qquad a^\dagger=\frac{x-D_x}{\sqrt{2}},$$9

with $[a,a^\dagger]=1+2\mu R$0, $[a,a^\dagger]=1+2\mu R$1 and the analogous $[a,a^\dagger]=1+2\mu R$2-operators. The Cartesian eigenstates are products of one-dimensional Dunkl oscillator wavefunctions, and the energy spectrum is

$[a,a^\dagger]=1+2\mu R$3

For fixed $[a,a^\dagger]=1+2\mu R$4, the level $[a,a^\dagger]=1+2\mu R$5 has $[a,a^\dagger]=1+2\mu R$6-fold degeneracy, as in the ordinary isotropic oscillator, but the symmetry algebra is deformed by reflections rather than being the undeformed Schwinger $[a,a^\dagger]=1+2\mu R$7 algebra (Genest et al., 2012). In three dimensions the same construction yields

$[a,a^\dagger]=1+2\mu R$8

with degeneracy $[a,a^\dagger]=1+2\mu R$9 (Genest et al., 2013).

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 (Genest et al., 2012, Genest et al., 2013). 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 $\partial_x$0, $\partial_x$1. The paper introduces Dunkl creation and annihilation operators linear in $\partial_x$2 and $\partial_x$3, with commutator

$\partial_x$4

The resulting relativistic spectrum for both the Dunkl–Klein–Gordon and Dunkl–Dirac oscillators is

$\partial_x$5

which depends explicitly on parity and reduces to the usual relativistic oscillator in the even sector when $\partial_x$6 (Hamil et al., 2022). 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

$\partial_x$7

or equivalently $\partial_x$8, $\partial_x$9, $R f(x)=f(-x)$0, one obtains symmetry generators commuting with the full Hamiltonian. The operator

$R f(x)=f(-x)$1

is the Dunkl angular momentum and controls polar separation. Together with the involutions $R f(x)=f(-x)$2, $R f(x)=f(-x)$3 and the central Hamiltonian, these operators generate the Schwinger–Dunkl algebra $R f(x)=f(-x)$4, an extension of $R f(x)=f(-x)$5 with reflections (Genest et al., 2012).

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 $R f(x)=f(-x)$6 Hahn polynomials, and the same structure is identified with the Clebsch–Gordan problem of $R f(x)=f(-x)$7 (Genest et al., 2012). The follow-up representation-theoretic analysis constructs Cartesian and circular bases from parabosonic creation and annihilation operators. In the Cartesian basis, $R f(x)=f(-x)$8 is diagonal and $R f(x)=f(-x)$9 acts tridiagonally; in the circular basis, $R^2=1$0 becomes block upper-triangular with $R^2=1$1 blocks, while $R^2=1$2 is tridiagonal. The Cartesian–circular expansion coefficients are Krawtchouk polynomials, and in the general anisotropic case the eigenvectors of $R^2=1$3 are generated by Heun polynomials with components expressed in para-Krawtchouk, ანუ complementary Bannai–Ito, polynomials; in the fully isotropic case they reduce to little $R^2=1$4 Jacobi or ordinary Jacobi polynomials (Genest et al., 2013).

In three dimensions, the same paradigm produces a reflection-deformed extension of $R^2=1$5. The one-dimensional Dunkl ladders $R^2=1$6 realize $R^2=1$7 in each coordinate, and bilinears built from them generate the Schwinger–Dunkl algebra $R^2=1$8. The corresponding constants of motion commute with the Hamiltonian, the system is maximally superintegrable, and separation is available in Cartesian, cylindrical, and spherical coordinates (Genest et al., 2013). 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 $R^2=1$9 realizations

A second major theme is that quadratic combinations of Dunkl ladders often generate undeformed $R x=-xR$0 algebras even though the underlying first-order commutator is reflection deformed. In the singular Dunkl oscillator, one-dimensional Hamiltonians

$R x=-xR$1

admit quadratic ladder operators

$R x=-xR$2

satisfying

$R x=-xR$3

These give an $R x=-xR$4 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 $R x=-xR$5 with central involutions (Genest et al., 2013).

The recent parametric-amplifier generalization makes the quadratic closure especially explicit. Defining

$R x=-xR$6

the commutation relations become

$R x=-xR$7

The Casimir acts as

$R x=-xR$8

with parity-dependent Bargmann indices

$R x=-xR$9

The Dunkl parametric-amplifier Hamiltonian

$R^2=1$00

is diagonalized both by an $R^2=1$01 tilting transformation and by a generalized Bogoliubov transformation, with stability condition $R^2=1$02 and exact spectrum

$R^2=1$03

For squeezed vacuum, the Mandel parameter is independent of $R^2=1$04, whereas $R^2=1$05 depends explicitly on $R^2=1$06 and modifies photon bunching (Guillén et al., 18 Feb 2026).

The same $R^2=1$07 mechanism appears in the Dunkl anharmonic, or Kerr, oscillator. There the Hamiltonian

$R^2=1$08

is rewritten in terms of $R^2=1$09 and $R^2=1$10, yielding parity-dependent energies

$R^2=1$11

The quantum dynamics exhibits collapse-and-revival phenomena with fundamental revival period $R^2=1$12 independent of $R^2=1$13, while $R^2=1$14 produces a perfect half-period revival at $R^2=1$15; the quadrature variance also shows interference-induced squeezing around $R^2=1$16 for specific deformation values (Ojeda-Guillén et al., 24 Apr 2026). 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

$R^2=1$17

which is identified as a Dunkl Anti–Jaynes–Cummings model. The same work computes parity-resolved thermal quantities from the spectrum,

$R^2=1$18

and shows that in the high-temperature limit the parity dependence washes out, with $R^2=1$19 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” (Hamil et al., 2022).

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,

$R^2=1$20

so Fermi–Dirac statistics are preserved, but the single-particle operator content is deformed by reflections. The grand partition function factorizes as

$R^2=1$21

equivalently corresponding to a Hamiltonian $R^2=1$22. The average occupation becomes

$R^2=1$23

and the paper imposes $R^2=1$24 to keep $R^2=1$25 for all energies (Zenkhri et al., 15 Aug 2025). 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 $R^2=1$26, 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 $R^2=1$27 and $R^2=1$28 symmetry, while the parametric-amplifier and Kerr constructions recover the standard squeezed-state and Kerr-medium formulas when $R^2=1$29 (Genest et al., 2012, Genest et al., 2013, Guillén et al., 18 Feb 2026, Ojeda-Guillén et al., 24 Apr 2026). 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.