2025-11-11T15:13:09.350045

Action-angle coordinates of spherical pendulums with symmetric quadratic potentials

Peng, Tang
We study the spherical pendulum system with an arbitrary potential function $V = V (z)$, which is an integrable system with a first integral whose Hamiltonian flow is periodic. We give an explicit solution to this integrable system and then we compute its action-angle coordinates. In the special case where the potential function is symmetric quadratic like $V = z^2$, we represent its action-angle coordinates in terms of elliptic integrals, and calculate the monodromy.
academic

Action-angle coordinates of spherical pendulums with symmetric quadratic potentials

Basic Information

  • Paper ID: 2509.04207
  • Title: Action-angle coordinates of spherical pendulums with symmetric quadratic potentials
  • Authors: Chengle Peng, Xiudi Tang
  • Classification: math.SG (Symplectic Geometry)
  • Publication Date: November 7, 2025
  • Paper Link: https://arxiv.org/abs/2509.04207

Abstract

This paper investigates spherical pendulum systems with arbitrary potential functions V=V(z)V = V(z), which constitute integrable systems with first integrals whose Hamiltonian flows are periodic. The article provides explicit solutions for this integrable system and computes its action-angle coordinates. For the special case where the potential is a symmetric quadratic function V=z2V = z^2, the authors express the action-angle coordinates in terms of elliptic integrals and compute the monodromy.

Research Background and Motivation

Problem Background

  1. Integrable Systems Theory: In Hamiltonian dynamics, Liouville integrable systems represent an important class of dynamical systems possessing the maximum number of functionally independent conserved quantities. Semitoric systems constitute a special subclass exhibiting rotational symmetry in one component.
  2. Significance of Spherical Pendulums: The spherical pendulum is one of the few natural Hamiltonian systems that is both integrable and possesses focus-focus singularities. In linear potential fields (such as gravity on Earth's surface), there exists a unique focus-focus point, which is the classical example discovered by Duistermaat where action-angle coordinates cannot exist globally.
  3. Challenges with General Potentials: When allowing more general potential functions, multiple focus-focus points and other singular behaviors may arise. For the case V(z)=z2V(z) = z^2, there are exactly two focus-focus points on the same fiber of the momentum map.

Research Motivation

  1. Theoretical Completeness: Although the topological properties of spherical pendulums have been studied, explicit computation of their action-angle coordinates remains incomplete.
  2. Technical Challenges: Existing research primarily focuses on Dirac-Poisson structures and topological properties, lacking comprehensive analysis of general spherical pendulums.
  3. Practical Applications: Action-angle coordinates are crucial for quantization and dynamical analysis.

Core Contributions

  1. Explicit Solutions for General Spherical Pendulums: For spherical pendulums with general potential function V(z)V(z), explicit solutions for the Hamiltonian flows JJ and HH are provided (Formula 3.10), with action coordinates expressed as definite integrals (Formula 5.8).
  2. Complete Analysis of Quadratic Potentials: Thorough investigation of the case V(z)=z2V(z) = z^2, providing explicit formulas for the joint flow in terms of elliptic integrals (Formula 4.3) and action-angle coordinates (Formula 5.6).
  3. Monodromy Computation: Calculation of the monodromy of the focus-focus fiber (Formula 5.7), consistent with the theory of focus-focus fibers in integrable systems.
  4. Explicit Computation of Symplectic Invariants: Provides one of the few examples of semitoric systems where symplectic invariants are actually computed.

Methodology Details

Problem Formulation

Study the spherical pendulum system (TS2,ω,F)(T^*S^2, \omega, F), where:

  • TS2T^*S^2 is the cotangent bundle of the two-dimensional sphere
  • ω\omega is the standard symplectic structure
  • F=(J,H)F = (J, H) is the momentum map, with JJ the angular momentum and HH the mechanical energy

The objective is to compute the action-angle coordinates of the system, particularly for the symmetric quadratic potential V(z)=z2V(z) = z^2.

Coordinate System Setup

Embedding Coordinates

Embed TS2T^*S^2 into TR3T^*\mathbb{R}^3 using position coordinates (x,y,z)(x,y,z) and momentum coordinates (u,v,w)(u,v,w), subject to constraints: x2+y2+z2=1,xu+yv+zw=0x^2 + y^2 + z^2 = 1, \quad xu + yv + zw = 0

Local Coordinate Charts

Three coordinate charts are defined:

  1. North Pole Coordinates UNU_N: using (ρ,η,θ,ϕ)(\rho, \eta, \theta, \phi)
  2. South Pole Coordinates USU_S: using (ρ,η,θ,ϕ)(\rho, \eta, \theta, \phi)
  3. Equatorial Coordinates UEU_E: using (z,w,θ,ϕ)(z, w, \theta, \phi)

Hamiltonian System

Momentum Map

J=xvyu,H=12(u2+v2+w2)+V(z)J = xv - yu, \quad H = \frac{1}{2}(u^2 + v^2 + w^2) + V(z)

where {J,H}=0\{J, H\} = 0, ensuring system integrability.

Image of Momentum Map

F(TS2)={(j,h)R2jR,hj22}F(T^*S^2) = \left\{(j,h) \in \mathbb{R}^2 \mid j \in \mathbb{R}, h \geq \frac{j^2}{2}\right\}

Derivation of Explicit Solutions

Equations of Motion

In north pole coordinates, the Hamiltonian equations are: dρdt=ηcosδ,dηdt=(ρη(1ρ2sin2δ)1ρ2+(1ρ2)V~N(ρ)η)ηcosδ\frac{d\rho}{dt} = \eta \cos \delta, \quad \frac{d\eta}{dt} = -\left(\frac{\rho\eta(1-\rho^2\sin^2\delta)}{1-\rho^2} + \frac{(1-\rho^2)\tilde{V}_N'(\rho)}{\eta}\right)\eta\cos\delta

where δ=ϕθ\delta = \phi - \theta, V~N(ρ)=V(1ρ2)\tilde{V}_N(\rho) = V(\sqrt{1-\rho^2}).

Separation of Variables

Using the conserved quantity ρηsinδ=j\rho\eta\sin\delta = j, the equations separate into: dρdt=ϵ(1ρ2)(2ρ2(hV~N(ρ))j2)ρ\frac{d\rho}{dt} = \epsilon\frac{\sqrt{(1-\rho^2)(2\rho^2(h-\tilde{V}_N(\rho))-j^2)}}{\rho}

This is a separable differential equation in ρ\rho and tt.

Elliptic Integral Solutions for Quadratic Potentials

For V(z)=z2V(z) = z^2, introduce parameters: k=1+h(1h)2+2j21+h+(1h)2+2j2,n=1+h(1h)2+2j22k = \sqrt{\frac{1+h-\sqrt{(1-h)^2+2j^2}}{1+h+\sqrt{(1-h)^2+2j^2}}}, \quad n = \frac{1+h-\sqrt{(1-h)^2+2j^2}}{2}

Explicit Expression of Joint Flow

ρ(s,t)=1nsin2γ(t)\rho(s,t) = \sqrt{1-n\sin^2\gamma(t)}θ(s,t)=θ0+ϵ2njkΠ(γ0,n,k)ϵ2njkΠ(γ(t),n,k)+s\theta(s,t) = \theta_0 + \epsilon\sqrt{2n}jk\Pi(\gamma_0,n,k) - \epsilon\sqrt{2n}jk\Pi(\gamma(t),n,k) + s

where γ(t)=am(F(γ0,k)ϵ2nt/k,k)\gamma(t) = \text{am}(F(\gamma_0,k) - \epsilon\sqrt{2n}t/k, k) is the Jacobi amplitude function.

Experimental Setup

Theoretical Verification

Since this is pure mathematical theoretical work, verification is conducted through:

  1. Symmetry Checks: Utilizing system symmetries to verify the correctness of period calculations
  2. Limiting Cases: Examining degenerate cases with special parameter values
  3. Monodromy Verification: Confirming that computed monodromy is consistent with known theory

Parameter Ranges

  • Elliptic modulus: 0<kn<10 < k \leq \sqrt{n} < 1
  • Momentum map values: (j,h)Br(j,h) \in B_r (regular value region)

Experimental Results

Period Lattice Computation

For regular values (j,h)Br(j,h) \in B_r, the period lattice is: Λj,h=(1,0)Z+(22kjnΠ(n,k),22k1nK(k))Z\Lambda_{j,h} = (1,0)\mathbb{Z} + \left(-2\sqrt{2}k\sqrt{\frac{j}{\sqrt{n}}}\Pi(n,k), 2\sqrt{2}k\frac{1}{\sqrt{n}}K(k)\right)\mathbb{Z}

Action Coordinates

A1=JA_1 = JA2=22π(n1nkK(k)+nkE(k)(n1)(k2n)nkΠ(n,k))A_2 = \frac{2\sqrt{2}}{\pi}\left(\frac{n-1}{\sqrt{n}}kK(k) + \frac{\sqrt{n}}{k}E(k) - \frac{(n-1)(k^2-n)}{\sqrt{n}k}\Pi(n,k)\right)

Monodromy Matrix

The monodromy around the focus-focus value (0,1)(0,1) is: dA+(dA)1=(1021)\text{dA}^+ \circ (\text{dA}^-)^{-1} = \begin{pmatrix} 1 & 0 \\ 2 & 1 \end{pmatrix}

This is in complete agreement with theoretical predictions for focus-focus fibers.

Special Case Verification

  1. j=0,h<1j = 0, h < 1: Λj,h=(1,0)Z+(0,22K(h))Z\Lambda_{j,h} = (1,0)\mathbb{Z} + (0, 2\sqrt{2}K(\sqrt{h}))\mathbb{Z}
  2. j=0,h>1j = 0, h > 1: Λj,h=(1,0)Z+(0,22hK(1h))Z\Lambda_{j,h} = (1,0)\mathbb{Z} + (0, \frac{2\sqrt{2}}{\sqrt{h}}K(\frac{1}{\sqrt{h}}))\mathbb{Z}

Classification of Semitoric Systems

  • Álvaro Pelayo and San Vũ Ngọc classified simple semitoric systems
  • Joseph Palmer and colleagues extended the classification to all semitoric integrable systems

History of Spherical Pendulum Research

  • Duistermaat discovered the global non-existence of action-angle coordinates in linear potential fields
  • Efstathiou studied the Dirac-Poisson structure for the case V(z)=z2V(z) = z^2
  • Dullin computed Taylor series invariants for spherical pendulums

Other Explicit Examples

  • Le Floch-Pelayo: coupled angular momenta
  • Alonso-Dullin-Hohloch: coupled spin oscillators

Conclusions and Discussion

Main Conclusions

  1. Complete Solution: First complete explicit expression of action-angle coordinates for spherical pendulums with symmetric quadratic potentials
  2. Elliptic Integral Representation: All key quantities are expressed in terms of elliptic integrals, facilitating numerical computation and theoretical analysis
  3. Monodromy Confirmation: Computational results are in complete agreement with general theory of focus-focus singularities

Limitations

  1. Coordinate Chart Restrictions: Complete computation requires transitions between multiple coordinate charts, increasing complexity
  2. Special Potential Functions: Detailed analysis is limited to quadratic potentials, with general potentials given only integral expressions
  3. Non-strict Semitoric: Since JJ is not a true circle map, the system does not strictly belong to semitoric systems

Future Directions

  1. Taylor Series Invariants: Compute Taylor series and other invariants for general spherical pendulums
  2. Other Potential Functions: Extend explicit computations to additional special potential functions
  3. Quantization Applications: Utilize action-angle coordinates for quantization studies

In-Depth Evaluation

Strengths

  1. Theoretical Completeness: Provides a complete derivation chain from Hamiltonian equations to action-angle coordinates
  2. Computational Precision: All formulas are explicit, facilitating verification and application
  3. Methodological Innovation: Cleverly combines elliptic integral theory with symplectic geometry techniques
  4. Result Significance: Represents one of the few semitoric systems where symplectic invariants are completely computed

Weaknesses

  1. Limited Scope: Complete results apply only to specific quadratic potential functions
  2. Technical Complexity: Requires deep background in elliptic functions and symplectic geometry for complete understanding
  3. Practical Limitations: As pure mathematical theoretical work, direct application scenarios are limited

Impact

  1. Theoretical Contribution: Provides important explicit examples for integrable systems theory
  2. Methodological Value: Demonstrates powerful applications of elliptic integrals in symplectic geometry
  3. Inspirational Significance: Provides methodological guidance for analysis of other integrable systems

Applicable Scenarios

  1. Theoretical Physics: Precise analysis of classical mechanical systems
  2. Mathematical Physics: Research in integrable systems and symplectic geometry
  3. Numerical Computation: Provides exact computational foundation for spherical pendulum systems

References

The paper cites 18 important references covering integrable systems theory, semitoric systems classification, elliptic integral theory, and other key areas, providing a solid theoretical foundation for the research.