Past exam of the mathematics course of the University of Cambridge / 2004 / iii / Paper 66 / 4

An isolating integral is a single-valued integral of motion $I(\mathbf{x},\mathbf{v})$ whose fixed value confines an orbit to a lower-dimensional level set of phase space. Along the flow, $dI/dt=0$; for an autonomous Hamiltonian, this means $\{I,H\}=0$. Several isolating integrals are independent when their gradients are linearly independent on the regular region considered. Counting them requires a specified domain: merely local orbit labels in a flow box do not establish global single-valued integrals. The following cases concern regular globally meaningful orbit constraints.

One integral can occur in an autonomous system with a fully chaotic region and no additional continuous orbit constraints. The specific orbital energy is conserved in a time-independent gravitational potential. If a chaotic stellar orbit is dense in a connected region of its energy hypersurface, every continuous integral of motion is constant on that region. Thus it cannot supply a new independent constraint there, and energy is the only independent isolating integral. A generic nonaxisymmetric potential can contain such regions; not every orbit in such a potential is necessarily chaotic.

No global isolating integral is possible in a sufficiently general time-dependent system without surviving symmetries or other invariants. Foreven energy is then not conserved, and nonsymmetric forcing can destroy all global constraints on the orbit. Time dependence alone does not imply this result: a central time-dependent force still conserves angular momentum, and a rigidly rotating potential can conserve a Jacobi integral of a rigidly rotating potential. The case of no integrals requires absence of such additional structure. A fixed point or locally defined flow labels are not the global orbit family considered here.

More than three integrals occur when additional symmetry or degeneracy restricts motion beyond a generic three-frequency invariant torus. In a time-independent central gravitational potential, energy and the three components of specific angular momentum are four generically independent isolating integrals. They fix the orbital plane as well as the radial-motion parameters. A nonclosed planar rosette typically has two-dimensional orbit closure, rather than a three-dimensional invariant torus.

The Kepler orbit is a stronger example. In addition to $E$ and $\mathbf{L}$, the conserved Laplace-Runge-Lenz vectorfixes the apsidal direction. Its constraints $\mathbf{A}\cdot \mathbf{L}=0$ and $A^2={\mu }^2+2E{L}^2$ leave five independent integrals among $E,\mathbf{L},\mathbf{A}$. It is a maximally superintegrable Hamiltonian system, with closed nondegenerate bound orbits.

More than three integrals need not form a set of Poisson-commuting functions; the three-integral limit applies to an independent commuting set, rather than to every conserved orbit constraint. At a nonstationary point in six-dimensional phase space, integral gradients annihilate the nonzero flow vector, so at most five are independent. This explains the isolating integrals and orbital closure dimension count. A resonance on one isolated orbit does not by itself prove an additional smooth global integral on a neighborhood.

Use $M$ for the mass parameter throughout; the PDF uses $m$ in its formulas and $M$ in the subsequent prose for the same quantity. Set $a=r_j>0$. The Jaffe model density simplifies toIts enclosed mass isIt tends to $M$ as $r\to \infty$, verifying the total mass. It tends to zero at the center, so the cusp carries no extra point mass.

Spherical symmetry gives ${\Phi }^{\prime }(r)=GM(r)/r^2=GM/[r(r+a)]$. Integrating with $\Phi (\infty )=0$ givesThis also directly satisfies the Poisson equation for Newtonian gravity: differentiating $r^2{\Phi }^{\prime }=GMr/(r+a)$ gives $4\pi G{r}^2\rho (r)$. The circular speed isFor $r\ll a$, $v_c=\sqrt{GM/a}[1-r/(2a)+O((r/a{)}^2)]$ is approximately constant. For $r\gg a$, $v_c=\sqrt{GM/r}[1-a/(2r)+O((a/r{)}^2)]\propto r^{-1/2}$. Thus the finite-mass model has an inner nearly flat galaxy rotation curve and an outer Keplerian decline.