---
title: Darboux III Oscillator
url: https://www.emergentmind.com/topics/darboux-iii-oscillator
type: topic
---

# Darboux III Oscillator

The Darboux III oscillator is an exactly solvable nonlinear deformation of the isotropic harmonic oscillator defined on an \(N\)-dimensional radially symmetric space of non-constant curvature. In its standard form, it is described by the Hamiltonian
\[
{\cal H}(q,p)=\frac{p^2}{2\,(1+\lambda q^2)}+\frac{\omega^2 q^2}{2\,(1+\lambda q^2)},
\qquad q,p\in\mathbb R^N,
\]
with deformation parameter \(\lambda\) and oscillator frequency \(\omega\). The model admits two complementary interpretations: as an intrinsic oscillator on the Darboux III manifold, and as a position-dependent-mass system with mass function \(m(q)=1+\lambda q^2\). For \(\lambda>0\), the underlying scalar curvature is negative and nonconstant; in the flat limit \(\lambda\to0\), the usual Euclidean harmonic oscillator is recovered [1411.7569, 2209.05293, 2212.13575]. In the literature on two-dimensional Darboux spaces, the label “Darboux III oscillator” is also used for an oscillator-type superintegrable potential on the Darboux III surface with quadratic symmetry algebra and deformed-oscillator realization [2301.03810].

## 1. Geometric setting and defining Hamiltonians

The \(N\)-dimensional Darboux III manifold is conformally flat, with line element
\[
\mathrm d s^2=(1+\lambda q^2)\,\mathrm d q\cdot\mathrm d q
=(1+\lambda q^2)\sum_{i=1}^N \mathrm d q_i^2,
\]
or, in radial coordinates,
\[
\mathrm d s^2=(1+\lambda r^2)\,\bigl(\mathrm d r^2+r^2\,\mathrm d\Omega_{N-1}^2\bigr),
\qquad r=|q|.
\]
Its scalar curvature is
\[
\mathcal R(r)
=-\,\lambda\,
\frac{(N-1)\bigl[2N+3(N-2)\,\lambda\,r^2\bigr]}
{(1+\lambda\,r^2)^3},
\]
so \(\lambda>0\) produces a negative, \(r\)-dependent curvature [2212.13575]. In two dimensions, the same geometry may be written as
\[
ds^2=(1+\lambda\,(x^2+y^2))\,(dx^2+dy^2),
\]
with scalar curvature
\[
R(r)=-\,\frac{4\,\lambda}{(1+\lambda r^2)^3},
\]
which tends to \(0\) as \(r\to\infty\) [2604.26466].

The standard classical Darboux III oscillator is
\[
{\cal H}(q,p)=\frac{p^2}{2(1+\lambda q^2)}+\frac{\omega^2 q^2}{2(1+\lambda q^2)}.
\]
In hyperspherical coordinates,
\[
{\cal H}
=\frac{p_r^2}{2(1+\lambda r^2)}
+\frac{\mathbf L^2}{2\,r^2(1+\lambda r^2)}
+\frac{\omega^2 r^2}{2(1+\lambda r^2)},
\]
where \(\mathbf L^2\) is the squared total angular momentum [2209.05293]. The same Hamiltonian is routinely interpreted as a flat-space system with position-dependent mass \(m(q)=1+\lambda q^2\) [2209.05293, 2212.13575].

The flat limit is smooth:
\[
\lim_{\lambda\to0}{\cal H}
=\frac{p^2}{2}+\frac{\omega^2 q^2}{2}.
\]
This makes the Darboux III oscillator a one-parameter deformation of the Euclidean oscillator rather than an unrelated curved-space model [1411.7569, 1511.08908].

## 2. Superintegrability, integrals of motion, and classical algebraic structure

The Darboux III oscillator is maximally superintegrable. In \(N\) dimensions it has \(2N-1\) functionally independent integrals of motion. One convenient set consists of the angular integrals
\[
C^{(m)}=\sum_{1\le i<j\le m}(q_i p_j-q_j p_i)^2,\qquad
C_{(m)}=\sum_{N-m<i<j\le N}(q_i p_j-q_j p_i)^2,
\]
for \(m=2,\dots,N\), together with the deformed Demkov–Fradkin integrals
\[
I_i=p_i^2-\bigl(2\,{\cal H}-\omega^2\bigr)\,q_i^2,
\qquad i=1,\dots,N,
\]
or, equivalently, the tensorial form
\[
I_{\lambda,ij}=p_i p_j-\bigl(2\,{\cal H}_{\lambda}-\omega^2\bigr)q_i q_j.
\]
These reduce to the standard Fradkin integrals when \(\lambda\to0\) [2209.05293, 1411.7569].

The model also admits a radial factorization and a non-linear Spectrum Generating Algebra. For fixed angular momentum \(\ell^2\), the reduced radial Hamiltonian can be written as
\[
H(r,p_r)=\frac{1}{1+\lambda r^2}\,H_0(r,p_r),
\]
with \(H_0=T_0+V_0\) the undeformed radial oscillator. The factorization introduces ladder functions \(A^\pm\) such that
\[
\{H,A^\pm\}=\mp i\,\alpha(H)\,A^\pm,\qquad
\{A^-,A^+\}=i\,\beta(H),
\]
where \(\alpha(H)\) and \(\beta(H)\) are explicit functions of the energy [1511.08908]. This non-linear Poisson algebra is the classical Spectrum Generating Algebra of the system.

Within that framework, the deformed frequency is energy-dependent:
\[
\Omega_E=\sqrt{\omega^2-\tfrac{2\lambda}{m}E}.
\]
The same structure is used to integrate the radial motion in closed form and to exhibit the deformation of the flat oscillator’s \(\mathfrak{sl}(2)\)-type ladder structure [1511.08908]. The Stäckel-transform perspective places the Darboux III oscillator in the broader family of maximally superintegrable curved systems obtained from the Euclidean oscillator, with the curved Fradkin tensor playing the rôle of the hidden symmetry generator [1103.4554].

## 3. Quantization, exact solvability, and spectral data

Several quantization prescriptions appear in the literature. A central one is the conformal Laplace–Beltrami prescription,
\[
\hat{\cal H}_{c,\lambda}
=-\frac{\hbar^2}{2}\,\Delta_{\rm conf}
+\frac{\omega^2 q^2}{2(1+\lambda q^2)},
\qquad
\Delta_{\rm conf}
=\Delta_{\rm LB}-\frac{N-2}{4(N-1)}\,R(q),
\]
which preserves maximal superintegrability [1411.7569]. For entropy calculations in the usual \(L^2(\mathbb R^N)\), a self-adjoint transformed PDM operator is used instead, supplemented by explicit \(\hbar\)-dependent ordering terms [2209.05293]. In the one-dimensional entropy analysis, the “left-ordering” kinetic term and an equivalent Hermitian “TPDM” form are both written explicitly [2510.06221]. These formulations are distinct at the operator level but are all used to study the same curved nonlinear oscillator.

The exact discrete spectrum is a smooth deformation of the harmonic-oscillator spectrum. In one common \(N\)-dimensional form,
\[
E_{n_1,\dots,n_N}(\lambda)
=
-\,\hbar^2\lambda\Bigl(n+\tfrac{N}{2}\Bigr)^2
+\hbar\Bigl(n+\tfrac{N}{2}\Bigr)
\sqrt{\hbar^2\lambda^2\Bigl(n+\tfrac{N}{2}\Bigr)^2+\omega^2},
\qquad
n=\sum_{i=1}^N n_i,
\]
while in radial-angular notation the same dependence is written through \(2n+\ell\) [2212.13575, 2209.05293]. The flat limit yields
\[
\lim_{\lambda\to0}E_{\lambda,n}
=\hbar\omega\Bigl(n+\tfrac{N}{2}\Bigr).
\]
The degeneracy remains exactly that of the flat \(N\)-dimensional oscillator,
\[
D(E_{\lambda,n})=\binom{n+N-1}{n},
\]
or equivalently \(D(N,k)=\binom{N+k-1}{k}\) for levels labeled by \(k=2n+\ell\) [1411.7569, 2209.05293]. For \(\lambda>0\), the discrete spectrum is bounded by
\[
E<E_{\rm max}=\frac{\omega^2}{2\lambda},
\]
with a continuous spectrum for \(E\ge \omega^2/(2\lambda)\) [2209.05293].

The bound-state eigenfunctions are also explicit. In hyperspherical variables,
\[
\psi_{n,\ell,\bm\mu}^\lambda(r,\bm\theta)
=
R_{n,\ell}^\lambda(r)\,Y_{\ell,\bm\mu}(\bm\theta),
\]
with radial factor
\[
R_{n,\ell}^\lambda(r)
=
N_{n,\ell}^\lambda\,
\sqrt{1+\lambda r^2}\;
r^\ell\,
e^{-\tfrac12\,\Omega_{n,\ell}^\lambda r^2}\;
L_n^{\left(\ell+\frac N2-1\right)}\!\bigl(\Omega_{n,\ell}^\lambda r^2\bigr),
\qquad
\Omega_{n,\ell}^\lambda=\sqrt{\omega^2-2\lambda E_{n,\ell}^\lambda},
\]
and \(Y_{\ell,\bm\mu}\) the standard hyperspherical harmonics [2209.05293]. In Cartesian separation, the wave functions may be written as products of Gaussian-Hermite factors with an energy-dependent width \(\beta=\sqrt{\Omega(E)/\hbar}\) [2212.13575]. Under conformal Laplace–Beltrami quantization, the bound-state wave functions carry the characteristic prefactor \((1+\lambda q^2)^{(2-N)/4}\) in \(L^2({\cal M}_\lambda)\) [1411.7569].

## 4. One- and two-dimensional reductions and information-theoretic analyses

In one dimension,
\[
H(x,p)=\frac{p^2}{2(1+\lambda x^2)}+\frac{\omega^2 x^2}{2(1+\lambda x^2)},
\]
and the quantum eigenvalues are
\[
E_n(\lambda)
=
-\lambda\,(n+\tfrac12)^2
+(n+\tfrac12)\sqrt{\lambda^2\,(n+\tfrac12)^2+\omega^2},
\qquad n=0,1,2,\dots.
\]
The normalized eigenfunctions are
\[
\psi_n(x;\lambda)
=
\mathcal N_n(\lambda)\,\sqrt{1+\lambda x^2}\,
e^{-\Omega_n(\lambda)x^2/2}\,
H_n\!\bigl(\sqrt{\Omega_n(\lambda)}\,x\bigr),
\qquad
\Omega_n(\lambda)=\sqrt{\omega^2-2\lambda E_n(\lambda)},
\]
with \(\mathcal N_n(\lambda)\) given explicitly in terms of \(\Omega_n\), \(n\), and \(\lambda\) [2510.06221]. For \(\lambda\to0\), both the energies and the Hermite-Gaussian wave functions reduce to the standard harmonic-oscillator ones [2510.06221].

A substantial recent literature studies the information-theoretic structure of these states. For Shannon entropy, analytical position-space results were obtained in arbitrary dimension, while momentum-space quantities were computed numerically because the Fourier transform of the Darboux III wave functions cannot be written in closed form [2209.05293]. In one dimension, analytical formulas were further derived for the position-space entropic moments and for the Rényi and Tsallis entropies; momentum-space Rényi and Tsallis entropies were again evaluated numerically from the Fourier-transformed density \(\gamma_n(p)\) [2510.06221].

The consistent pattern across these studies is that increasing \(\lambda\) broadens the position-space density and narrows the momentum-space density. In the Shannon analysis, \(S_\rho\) increases with \(\lambda\), \(S_\gamma\) decreases with \(\lambda\), and the Bialynicki–Birula–Mycielski uncertainty bound remains valid [2209.05293]. In the Rényi/Tsallis analysis, position-space entropies generally increase with \(\lambda\), momentum-space entropies decrease with \(\lambda\), and large \(\lambda\) or highly excited states can induce non-monotonic dependence on \(n\) [2510.06221].

For large \(\lambda\) and/or large \(n\), a specific approximation becomes effective. Writing
\[
\rho_n(x)=\rho_n^{(0)}+\lambda\,\rho_n^{(2)},
\]
the term \(\lambda\rho_n^{(2)}\) dominates, and the wave function is approximated by
\[
\phi_n(x;\lambda)=\sqrt{\lambda}\,\mathcal N_n\,|x|\,e^{-n x^2/2}H_n(\sqrt n\,x).
\]
This approximation has a closed-form Fourier transform involving incomplete gamma functions, and for \(n=0,1,2,3\) it can be written in terms of Dawson’s \(F(x)\); comparison with the exact numerical momentum density shows excellent agreement for \(\lambda\gg1\) [2510.06221].

## 5. Two-dimensional Darboux-space oscillator models and polynomial symmetry algebras

In the algebraic theory of superintegrable systems on two-dimensional Darboux spaces, one Darboux III family can be chosen to play the rôle of an isotropic “oscillator.” In the notation of the classification used in the 2D Darboux-space literature, this is the potential labeled \(V_1\), with the two cross couplings set to zero so that only the radial parameter remains [2301.03810]. The corresponding Hamiltonian is
\[
H=-\frac{\hbar^2}{2m}\,\Delta_g+\frac{\omega^2}{4e^u+1},
\]
and it admits two independent quadratic integrals,
\[
Q_1=\partial_v^2+\frac{\omega^2}{4e^u+1},
\]
and
\[
Q_2
=
e^{-u}
\Bigl[
\sin v\,\partial_u\partial_v-\cos v\,\partial_v^2
\Bigr]
-\omega^2\,\frac{e^u\sin v}{4e^u+1}.
\]
With \(A=Q_1\), \(B=Q_2\), the quadratic algebra is
\[
[A,B]=C,\qquad
[A,C]=-8\,\omega^2\,B,\qquad
[B,C]=2\,A^2-4\,H\,A+12\,\omega^2,
\]
and its Casimir is
\[
K=C^2-2\{A^2,B\}+4H\{A,B\}-4\omega^2B^2-8\omega^2A
=16\,\omega^4-8\,\omega^2H.
\]
A deformed-oscillator realization with
\[
[N,b^\dagger]=b^\dagger,\qquad [N,b]=-b,\qquad
bb^\dagger=\Phi(N+1),\qquad b^\dagger b=\Phi(N)
\]
leads to the structure function
\[
\Phi(N)
=
4\omega^2\,
(N+\eta)\,
\bigl(N+\eta-p-1\bigr)\,
\bigl(E-2\omega(N+\eta)\bigr),
\]
and the finite-dimensional unirreps yield the equally spaced spectrum
\[
E_n=2\,\omega\,n,\qquad n=1,2,\dots,
\]
with \(n\)-fold degeneracy [2301.03810].

This two-dimensional “oscillator” should be distinguished from the \(N\)-dimensional radially symmetric nonlinear oscillator discussed above. The former belongs to the Darboux-space classification of second-order superintegrable systems and is organized by quadratic symmetry algebras; the latter is the curved isotropic oscillator on the conformally flat manifold \(ds^2=(1+\lambda q^2)\,dq^2\). The common terminology reflects the shared Darboux III background and oscillator-type spectral behavior, but the realizations are not identical.

A recurrent misconception concerns the 2023 review of polynomial symmetry algebras in Darboux spaces. That paper explicitly reviews three approaches—deformed-oscillator constructions for finite-dimensional representations, induced-module-type constructions for infinite-dimensional representations, and commutant methods for discovering new models—but its only fully worked example is the Darboux II oscillator. Although Darboux III is mentioned among the twelve second-order superintegrable models, the paper does not write down, for Darboux III, the quantum Hamiltonian, the quadratic integrals, the commutation relations, the deformed-oscillator mapping, the closed-form structure function, the finite-dimensional unirreps, the explicit energy levels, or the corresponding eigenfunctions. For a complete Darboux III algebraic solution, one must therefore consult earlier specialized treatments rather than that review [2309.04928].

## 6. Magnetic couplings, deformations with reflections, and broader developments

The Darboux III oscillator has served as the base model for several exactly solvable extensions. One is the Dunkl-Darboux III oscillator, introduced as a \(\lambda\)-deformation of the \(N\)-dimensional Dunkl oscillator. In this formulation the deformation can again be interpreted either as non-constant curvature or as a position-dependent mass function, and the resulting quantum model remains exactly solvable in arbitrary dimension [2212.13575].

In two dimensions, the oscillator can be coupled to a constant magnetic field. Within the Dunkl-Darboux III framework this is achieved by the symmetric-gauge substitution
\[
\mathbf p\to \mathbf p-\frac{e}{c}\mathbf A,
\]
leading to an exactly solvable Hamiltonian whose spectrum depends on \(\lambda\), the oscillator parameters, and the magnetic field [2212.13575]. A related development is the Fock-Darwin-Darboux system, defined as the generalization of the Fock-Darwin system to a particle moving on the Darboux III space. When the magnetic field vanishes, it reduces to the two-dimensional Darboux III oscillator on the plane [2604.26466].

For the two-dimensional Laplace–Beltrami realization with vanishing magnetic field, separation
\[
\Psi(r,\theta)=R_{n,l}(r)\,\frac{e^{i l\theta}}{\sqrt{2\pi}},
\qquad l\in\mathbb Z,
\]
gives discrete levels
\[
E_{n,l}
=
-\,\lambda\,(2n+|l|+1)^2
+(2n+|l|+1)\sqrt{\omega^2+\lambda^2(2n+|l|+1)^2}.
\]
Writing \(c=2n+|l|+1\), one has \(E_{n,l}=E_c\), so the level labeled by the principal quantum number \(c\) has degeneracy \(c\), exactly as in the flat oscillator. The normalized eigenfunctions are Laguerre-type modes deformed by the conformal factor \(\sqrt{\mu(r)}\), with \(\mu(r)=1+\lambda r^2\) and effective frequency \(\Omega_{n,l}=\sqrt{\omega^2-2\lambda E_{n,l}}\) [2604.26466].

The magnetic generalization also changes the Landau problem in a specific way: on the Darboux III space, the Landau system has no infinitely degenerate Landau levels [2604.26466]. This is one manifestation of how curvature deforms flat-space spectral structures while preserving exact solvability.

Current directions stated in the recent entropy literature include the study of the higher-dimensional Darboux III oscillator, the case of negative \(\lambda\), local information measures, complexity, and quantum-gravitational position-dependent-mass models [2510.06221]. Earlier classical analyses also note that for \(\lambda<0\) the same formulas continue to hold on the open ball \(0<r<r_s=(|\lambda|)^{-1/2}\), and that for \(E>m\omega^2/(2\lambda)\) with \(\lambda>0\), the bounded-motion trigonometric structure is replaced by a hyperbolic one corresponding to unbounded trajectories [1511.08908]. These developments place the Darboux III oscillator at the intersection of superintegrability, curved-space quantization, polynomial symmetry algebras, and information-theoretic analysis.

Source: https://www.emergentmind.com/topics/darboux-iii-oscillator