Elliptic Plane Pendulum Equation

• The elliptic plane pendulum equation is a nonlinear oscillator model that employs elliptic functions to describe both librational and rotational motions with exact analytic solutions.
• It bridges classical mechanics and integrable systems by reformulating energy conservation into a first-order ODE solvable via Jacobi elliptic functions and complete elliptic integrals.
• Modular symmetries and duality maps in the equation reveal deep links between phase space geometry, action-angle variables, and modern perturbation techniques in nonlinear dynamics.

The elliptic plane pendulum equation generalizes the classical nonlinear oscillator by encoding its dynamics through elliptic functions, establishing an exact correspondence between the geometry of phase space, conserved energy surfaces, and modular symmetries. Its solutions elucidate the full range of oscillatory and rotational regimes for arbitrary initial data, bridging classical mechanics, integrable systems, and special function theory.

1. Nonlinear Pendulum Equation and Elliptic Formulation

The starting point is the equation of motion for a simple plane pendulum of length $L$ and mass $m$ in a uniform gravitational field $g$: $$\frac{d^2\theta}{dt^2} + \frac{g}{L} \sin\theta = 0.$$ This yields, upon multiplication by $\dot\theta$ and integration, the first energy integral

$$E = \tfrac12 m L^2 \dot\theta^2 + m g L (1 - \cos\theta).$$

By introducing the energy parameter $k_e^2 \equiv E/(2 m g L)$ and dimensionless time $x \equiv \sqrt{g/L}\, t$, the system reduces to

$$\frac{1}{4}\left(\frac{d\theta}{dx}\right)^2 + \sin^2\left(\frac{\theta}{2}\right) = k_e^2,$$

which leads to a first-order ODE with the natural phase-space foliation parameterized by energy. This reduction exposes the system's integrability and admits explicit solution in terms of elliptic functions (Linares, 2016).

2. Elliptic Integral Solution and Jacobi Functions

Separation of variables and a trigonometric substitution, $\sin(\theta/2)=k_e \sin\phi$, render the time–angle relation in the form

$$t = 2\sqrt{\frac{L}{g}} \int_0^{\theta/2} \frac{du}{\sqrt{k_e^2 - \sin^2 u}},$$

which, suitably reparametrized, is inverted by Jacobi's elliptic sine function as

$$\sin\left(\frac{\theta}{2}\right) = k_e \, \mathrm{sn}(x \mid k_e), \quad \theta(t) = 2\arcsin \left[k_e\, \mathrm{sn}(\sqrt{g/L}\, t; k_e)\right].$$

The period of small oscillations emerges in terms of the complete elliptic integral of the first kind,

$$T(k_e) = 4\sqrt{\frac{L}{g}}\, K(k_e), \quad K(m) = \int_0^{\pi/2} \frac{d\phi}{\sqrt{1-m \sin^2\phi}}.$$

For arbitrary energy, this yields a closed-form, analytic description of all librational and rotational trajectories (Linares, 2016, Haine, 2023, Klee, 2016).

3. Modular Symmetries, Duality, and Global Solution Structure

The Jacobi sn-solution inherits a 6-fold symmetry from the modular group $\Gamma \simeq PSL(2, \mathbb{Z})$. The modulus $k^2$ transforms under $\Gamma/\Gamma(2)$ as

$$k^2 \leftrightarrow 1-k^2 \leftrightarrow 1/k^2 \leftrightarrow 1-1/k^2 \leftrightarrow 1/(1-k^2) \leftrightarrow k^2/(k^2-1),$$

which relates distinct but physically equivalent representations—oscillatory (libration) and circulating (rotation)—via duality maps.

In particular, the S-duality $k_e \to k_e' = \sqrt{1 - k_e^2}$ exchanges real (physical) and imaginary (inverted-gravity) time solutions. Explicitly, the imaginary-time solution is connected to the real by analytic continuation and Jacobi imaginary-modulus identities: $$\theta(x) = 2\arcsin[k_e \mathrm{sn}(x; k_e)] = 2\arcsin[\mathrm{dn}(i x - iK + K_e'; k_e')] .$$ A single, unified Jacobi function analytic expression, supplemented by modular transformations, covers $0 < k_e^2 < \infty$, encompassing oscillatory, critical, and rotational regimes (Linares, 2016, Haine, 2023).

4. Exact Summation Methods and Generalizations

Beyond analytic inversion, series-based and resummation approaches yield exact global solutions and efficient computational schemes. The functional representation as an elliptic function dictates the complex-time singularity structure and thus the convergence properties of power-series expansions. Series about the trajectory top, via recurrence for Taylor coefficients, achieve maximal convergence radius, bounded by the nearest complex singularity (Reinberger et al., 2021). Asymptotic approximant techniques further accelerate convergence: $$K(k) = \frac{\pi}{2}\sum_{n=0}^{\infty} \left[ \frac{(2n)!}{(n!)^2} \right]^2 16^{-n}k^{2n},$$ with accelerated forms combining elementary functions and logarithms for $k \approx 1$.

Generalizations to arbitrary initial conditions and nontrivial phase or energy are achieved by expressing $\theta(t)$ as

$$\theta(t) = 2\arcsin\left(\sqrt{k}\, \gamma_0\, \mathrm{sn}(u(t), m)\right)$$

with parameter mappings determined by initial data, and with the modulus $m$ incorporating energy and phase shift contributions (Neto, 2010).

Alternative approaches, such as closed-form solutions via Weierstrass elliptic functions, map the pendulum dynamics directly onto motion on a real torus, highlighting the complete equivalence of the Jacobi and Weierstrass representations (Knill, 2023).

5. The Elliptic Pendulum Equation and Deformation of Restoring Force

A systematic deformation of the pendular force law yields the elliptic plane pendulum equation: $$\ddot\theta = -a\, \frac{\mathrm{sn}(\theta, m)}{\mathrm{dn}(\theta, m)}, \quad a = g/l.$$ This interpolates between the standard plane pendulum ($m = 0$), the hyperbolic (sinh) pendulum ($m = 1$), and a continuous one-parameter family. For $0 \leq m \leq 1$, the harmonic approximation

$\ddot\theta = -a\, \theta$

remains isochronous—independent of amplitude—across the family.

At $m = 1/2$, the cubic (anharmonic) correction in the Taylor expansion vanishes, yielding exact isochrony even at first anharmonic order. The motion in finite-amplitude regimes is implicitly specified by an elliptic integral involving $\int sn(u,m)/dn(u,m)\, du$; a closed-form explicit $\theta(t)$ in terms of standard elliptic functions is generally not available (open problem) (Khare et al., 13 Dec 2025).

6. Elliptically Excited Pendula, Parametric Forcing, and Generalizations

Driven or parametrically forced pendula with elliptical support or pivot trajectories introduce additional complexity. The equations of motion for such systems admit explicit, but more intricate, reduction: $$\theta'' + \beta\theta' + \mu \sin(\tau + \theta) = \delta \sin(\tau - \theta) - \omega^2 \sin\theta,$$ where parameters $\mu$, $\delta$, $\beta$, and $\omega$ encode drive amplitude, ellipticity, damping, and frequency scales. When the eccentricity $\epsilon$ is zero (circular excitation) and gravity vanishes, exact rotational solutions exist; otherwise, perturbative and averaging methods yield accurate approximations in both high- and low-damping regimes (Belyakov, 2010, Bouzas, 2011).

These frameworks generalize the elliptic pendulum paradigms, embedding the classical dynamics in broader classes of integrable and near-integrable systems, including connections to the sine-Gordon and sine-hyperbolic-Gordon equations (Khare et al., 13 Dec 2025).

7. Canonical Structure, Action-Angle Variables, and Perturbation Theory

Action-angle coordinates for the nonlinear pendulum are constructed via the Jacobi elliptic functions and integrals, providing a canonical transformation from angular position and momentum $(\vartheta, p)$ to $(J,\zeta)$. The generator is the Jacobi zeta function, with the action for libration and rotation given by

$$J_\ell(m) = \frac{8}{\pi} [E(m) - (1-m)K(m)], \quad J_r(m) = \frac{4\sqrt{m}}{\pi} E(1/m),$$

where $E(m)$ is the complete elliptic integral of the second kind.

On this foundation, systematic perturbation ansätze generate arbitrary-order corrections beyond the phase-space harmonic approximation, reconstructing the full period–amplitude relation and trajectory expansions in $\sin^2(\theta_0/2)$ (Klee, 2016, Brizard, 2011). These methods establish the interplay between classical integrability, special function theory, modular symmetry, and modern techniques for nonlinear oscillators.

---

References

• Duality symmetries behind solutions of the classical simple pendulum (Linares, 2016)
• On a generalization of Jacobi's elegantissima (Haine, 2023)
• Plane Pendulum and Beyond by Phase Space Geometry (Klee, 2016)
• Some Novel Aspects of the Plane Pendulum in Classical Mechanics (Khare et al., 13 Dec 2025)
• A General Exact Closed-Form Solution for Nonlinear Differential Equation of Pendulum (Dalir, 2020)
• Nonlinear Pendulum: A Simple Generalization (Neto, 2010)
• Action-angle coordinates for the pendulum problem (Brizard, 2011)
• On rotational solutions for elliptically excited pendulum (Belyakov, 2010)
• Steady states of the parametric rotator and pendulum (Bouzas, 2011)
• On the trajectory of the nonlinear pendulum: Exact analytic solutions via power series (Reinberger et al., 2021)
• Weierstrass elliptic functions for the pendulum (Knill, 2023)