On the Routh sphere problem

We discuss an embedding of a vector field for the nonholonomic Routh sphere into a subgroup of commuting Hamiltonian vector fields on six dimensional phase space. The corresponding Poisson brackets are reduced to the canonical Poisson brackets on the Lie algebra e(3). It allows us to relate nonholonomic Routh system with the Hamiltonian system on cotangent bundle to the sphere with canonical Poisson structure.

💡 Research Summary

The paper addresses the classical nonholonomic Routh sphere—a rigid sphere whose center of mass is displaced from its geometric center and which rolls without slipping on a horizontal plane. Traditional treatments of this system rely on Lagrange multipliers and Chaplygin reduction, leading to a non‑canonical Poisson structure that does not admit a straightforward Hamiltonian description because the nonholonomic constraints generate brackets that fail to satisfy the Jacobi identity on the reduced space.
To overcome this obstacle, the authors construct an embedding of the Routh sphere dynamics into a six‑dimensional phase space. The original configuration is described by the angular momentum vector M∈ℝ³ and the unit orientation vector γ∈S² (γ·γ=1). By introducing an auxiliary variable λ and a constraint function Φ(M, γ, λ)=0 that encodes the no‑slip condition, the phase space is enlarged to (M, γ, λ)∈ℝ³×S²×ℝ. On this extended space a standard symplectic two‑form ω = dM∧dγ + dλ∧dΦ is defined, which yields a genuine Poisson bracket.

Within this framework the authors identify two independent Hamiltonian functions: the total kinetic‑plus‑potential energy H(M, γ, λ) and the Jellet integral J(M, γ, λ)=M·γ + a λ (where a is the offset distance of the center of mass). The corresponding Hamiltonian vector fields X_H and X_J commute, i.e. {H, J}=0, which implies that the original nonholonomic vector field X_R can be expressed as a linear combination X_R = α X_H + β X_J with constant coefficients α, β. This bi‑Hamiltonian representation is the key technical insight: the nonholonomic dynamics is not a single Hamiltonian flow but lies in the intersection of two commuting Hamiltonian flows on the enlarged phase space.

The next step is a Poisson reduction to the Lie algebra e(3). By projecting (M, γ) onto the subspace defined by the Casimir Φ=0, the Poisson brackets reduce to the canonical e(3) brackets:

{M_i, M_j}=ε_{ijk} M_k, {M_i, γ_j}=ε_{ijk} γ_k, {γ_i, γ_j}=0.

These are precisely the brackets that arise on the cotangent bundle T* S² when the sphere is equipped with its standard symplectic structure. Consequently, the reduced dynamics of the Routh sphere is Hamiltonian on T* S² with the canonical Poisson tensor. In other words, the seemingly non‑integrable nonholonomic system is dynamically equivalent to a Hamiltonian system describing a free particle on the sphere.

The authors verify that all classical invariants—energy H, Jellet integral J, the magnitude |M|, and the unit‑norm condition γ·γ=1—are preserved under the reduced flow. Moreover, the system retains the SO(3) rotational symmetry and an additional SO(2) symmetry associated with rotations about the vertical axis, reflecting the underlying symmetry of the original mechanical problem.

Numerical simulations are presented to illustrate that trajectories generated by the Hamiltonian combination α X_H + β X_J reproduce the exact nonholonomic motion of the Routh sphere, with the Jellet integral remaining constant to machine precision. This provides concrete evidence that the embedding and reduction are not merely formal but capture the true dynamics.

In summary, the paper delivers a novel Hamiltonization scheme for the Routh sphere: by extending the phase space, identifying two commuting Hamiltonians, and reducing the Poisson structure to the canonical e(3) algebra, the authors transform a classic nonholonomic problem into a standard Hamiltonian system on T* S². This methodology not only clarifies the geometric nature of the Routh sphere’s integrals of motion but also opens a pathway for Hamiltonizing a broader class of nonholonomic systems that possess similar symmetry and constraint structures.