Nonlinear Lissajous orbits and particular superintegrability

Abstract: We investigate the geometry of classical trajectories generated by separable two-dimensional polynomial potentials of the form $$V(x,y)=\tfrac{1}{2}\big(x<sup>{2N}+A\,y<sup>{2N}\big)$$, where $N=1,2,\ldots,$ and $A&gt;0$. Special emphasis is placed on the emergence of nonlinear Lissajous figures and on the distinction between global and particular superintegrability in the Liouville sense. In the harmonic case ($N=1$) closed periodic orbits are a consequence of an additional \emph{global} integral of motion whenever the frequency ratio is rational, rendering the system maximally superintegrable. In contrast, for anharmonic oscillators, already in the quartic case ($N=2$), the oscillation frequencies depend on the partial energies, so periodic Lissajous-type trajectories occur only under nonlinear resonance conditions fixed by the initial data. Accordingly, the extra conserved quantities that characterize these closed orbits are not global invariants but \emph{particular} (trajectory-dependent) integrals that emerge only on the resonant trajectories. For higher-degree potentials $N\geq3$, the resonant trajectories are naturally described by hyperelliptic phase constraints rather than by a universal polynomial orbit equation.

• The paper establishes that closed nonlinear Lissajous orbits arise on energy-dependent resonant manifolds, distinguishing particular superintegrability from global superintegrability.
• It solves the quartic case with Jacobi elliptic functions, deriving explicit algebraic orbit curves for 1:n resonances such as n=2 and n=3, while showing that resonance frequencies depend on partial energies.
• For degrees N≥3, it replaces elliptic solutions with hyperelliptic phase variables and demonstrates that phase-based particular integrals remain conserved only on resonant shells, where no universal polynomial orbit equation is generally available.

Overview

This paper by Escobar-Ruiz and Azuaje studies the classical dynamics of two-dimensional separable polynomial potentials of the form $V(x,y)=\tfrac12(x^{2N}+Ay^{2N})$ with $A>0$, with the aim of characterizing the configuration-space geometry of closed trajectories — nonlinear analogues of Lissajous figures — and the conserved quantities associated with them. The central conceptual contribution is a sharp distinction between global superintegrability, in which an additional integral exists throughout phase space, and particular superintegrability, in which trajectory-dependent integrals are conserved only on resonant invariant manifolds selected by nonlinear, energy-dependent resonance conditions.

The analysis proceeds through three regimes: the harmonic case ($N=1$), where closed orbits reflect global superintegrability; the quartic case ($N=2$), where the dynamics is elliptic and resonant orbits can be written as explicit algebraic curves via Jacobi elliptic multiplication formulas; and the general case $N\geq3$, where hyperelliptic phase variables replace polynomial orbit equations.

The harmonic oscillator: global superintegrability

For $N=1$, the frequencies $\omega_x=\sqrt{1/m}$ and $\omega_y=\sqrt{A/m}$ are energy-independent, so closure is controlled solely by the frequency ratio $\omega_x/\omega_y = 1/\sqrt{A}$. When this ratio is rational, $p/q$, the motion is periodic with period $A>0$0, and eliminating time from the amplitude-phase solutions yields the implicit Lissajous relation

$A>0$1

For the phase-locked branch $A>0$2 at $A>0$3, these reduce to explicit algebraic curves (a line, a parabola-like quadratic, and a cubic). Crucially, for each rational ratio there exists a third global integral $A>0$4 — angular momentum $A>0$5 for $A>0$6, a cubic integral for $A>0$7, a quartic integral for $A>0$8 — whose Poisson bracket with the Hamiltonian vanishes identically on all of phase space. The system is therefore maximally superintegrable: three functionally independent integrals for two degrees of freedom, and every bounded trajectory is closed.

The quartic oscillator: elliptic dynamics and nonlinear resonance

For $A>0$9, both separated motions are solved exactly in terms of Jacobi elliptic functions $N=1$0 with fixed modulus $N=1$1. The physical frequencies scale as $N=1$2, giving the nonlinear frequency ratio

$N=1$3

which depends explicitly on the partial energies. This is the structural departure from the harmonic case: closure is no longer a property of the potential alone but a condition on initial data. On the equal-energy shell $N=1$4, choosing $N=1$5 realizes the $N=1$6 resonance, and the phase-locked orbit ($N=1$7) follows from the multiplication theorem for $N=1$8:

$N=1$9

where $N=2$0, $N=2$1, and $N=2$2 are polynomials determined by the modulus. Explicitly, for $N=2$3 the orbit is $N=2$4, and for $N=2$5 a degree-nine curve; note that the polynomial degree grows with resonance order, as in the harmonic case.

The associated conserved quantities are particular integrals: functions $N=2$6 satisfying $N=2$7 on the resonant manifold $N=2$8 but not identically on phase space. A representative example is the angular momentum: for the isotropic quartic oscillator,

$N=2$9

so $N\geq3$0 is not global, yet it vanishes and is conserved on the invariant lines $N\geq3$1. Analogous energy-shell representatives $N\geq3$2, $N\geq3$3 are constructed for the $N\geq3$4 and $N\geq3$5 resonances using the identity $N\geq3$6.

General degree: hyperelliptic regime

For $N\geq3$7 the separated motions are governed by hyperelliptic integrals. The authors introduce a generalized cosine $N\geq3$8, defined as the real $N\geq3$9-periodic solution of $N=1$0, where $N=1$1. The frequency scaling becomes

$N=1$2

so the resonance condition reads $N=1$3, or equivalently $N=1$4 for a $N=1$5 resonance. This makes explicit that nonlinear Lissajous trajectories are selected jointly by the anisotropy parameter and the partial-energy ratio, not by $N=1$6 alone — a point illustrated by resonance curves in the $N=1$7 plane.

On the equal-energy shell with $N=1$8, the resonant solution is $N=1$9, $\omega_x=\sqrt{1/m}$0, and the orbit is given implicitly by

$\omega_x=\sqrt{1/m}$1

or equivalently by a hyperelliptic integral constraint modulo $\omega_x=\sqrt{1/m}$2. The paper emphasizes that no universal polynomial relation $\omega_x=\sqrt{1/m}$3 exists in this regime — a genuine qualitative change from both the harmonic and quartic cases. The corresponding particular integrals are the phase combinations

$\omega_x=\sqrt{1/m}$4

with single-valued representatives built from trigonometric functions of the phase difference; their Poisson brackets with $\omega_x=\sqrt{1/m}$5 vanish only after restriction to the resonant shell.

Nontriviality and the obstruction to globalization

A substantive methodological point concerns whether these constructions are tautological. The authors argue they are not: the particular integrals arise intrinsically from the phase-locking condition $\omega_x=\sqrt{1/m}$6 on resonant invariant manifolds, and hence characterize entire resonant families rather than individual parametrized curves. Their nontrivial content is twofold — constructive (they define the resonant manifolds) and negative (their brackets fail to vanish globally).

The paper also clarifies what such relations do not mean. Restricted to a resonant manifold, functional independence among phase-space functions drops, and additional relations appear; but these should not be read as remnants of harmonic superintegrability. In particular, the authors explicitly warn against representing the $\omega_x=\sqrt{1/m}$7 resonance by conditions such as $\omega_x=\sqrt{1/m}$8 or straight-line orbits $\omega_x=\sqrt{1/m}$9, and against identities relating quartic particular integrals to harmonic partial energies — the latter are merely auxiliary functions under the quartic flow and are not conserved. The obstruction to promoting a particular integral to a global one is dynamical rather than merely algebraic: because the resonance condition itself depends on the partial energies, no fixed phase-space function can encode it globally.

Limitations and open questions

Several restrictions bound the scope of the results. The analysis is confined to separable homogeneous even-power potentials; non-separable perturbations, which would destroy the partial-energy integrals underlying the entire construction, are not treated. The explicit algebraic orbit equations and phase-space representatives are derived only for the equal-energy shell $\omega_y=\sqrt{A/m}$0 and the phase-locked branch $\omega_y=\sqrt{A/m}$1; other relative phases yield distinct curves whose explicit forms are not written out. For $\omega_y=\sqrt{A/m}$2, the claim that no universal polynomial orbit equation exists is asserted rather than proven in full generality, and the possibility that special low-order cases reduce to algebraic relations is acknowledged but not systematically classified. Finally, the treatment is entirely classical; the quantum analogue of particular superintegrability on resonant manifolds — e.g., whether trajectory-dependent integrals induce spectral degeneracies or quasi-exact solvability — remains open.

Conclusion

The paper establishes that in separable anharmonic oscillators, closed Lissajous-type motion survives only on resonant invariant manifolds selected by energy-dependent nonlinear resonance conditions, and that the associated conserved quantities are genuinely particular integrals — conserved on those manifolds but not globally. Integrability is thereby shown to reorganize rather than disappear when global superintegrability is lost: Liouville integrability persists through separability, while the additional structure localizes on resonant sets defined by hyperelliptic (for $\omega_y=\sqrt{A/m}$3) or elliptic-algebraic (for $\omega_y=\sqrt{A/m}$4) phase constraints.