Cyclic Coordinates and Conserved Momenta

Write the central-force problem in Cartesian coordinates and the Hamiltonian is an honest mess: xx and yy both appear, tangled inside V(x2+y2)V(\sqrt{x^2+y^2}), and nothing is obviously conserved. Write the same problem in polar coordinates and the angle θ\theta simply vanishes from HH. Hamilton's equation p˙θ=−∂H/∂θ\dot p_\theta=-\partial H/\partial\theta then reads p˙θ=0\dot p_\theta=0 before you have solved anything: angular momentum is conserved, and a two-dimensional problem collapses to one radial equation plus a final integral. A coordinate that is absent from the Hamiltonian in this way is called cyclic (or ignorable), and it is the single cheapest conservation law in physics — you get it by inspection. This page is the Hamiltonian half of the story, and the half with the sharper payoff. You have already seen (in Cyclic Coordinates and Conservation of Conjugate Momenta) that a coordinate missing from the Lagrangian conserves its conjugate momentum; here the same fact becomes a genuine reduction: the conserved momentum turns into a mere number in the reduced Hamiltonian, the cyclic pair drops out of the active dynamics entirely, and Routh's hybrid construction packages the whole procedure into an effective one-dimensional problem. But the deepest lesson is a provocation. Whether a coordinate is cyclic is a property of your coordinate choice, not of the system — the same orbit hides its conservation law in Cartesian variables and displays it in polar ones. So why stop at one good coordinate? The dream of making every coordinate cyclic is precisely what launches the theory of canonical transformations, the chapter this article opens the door to.

Key takeaways
  • A coordinate qkq_k is cyclic (ignorable) if it does not appear explicitly in the Hamiltonian — even though its conjugate momentum pkp_k usually does.
  • Hamilton's equation p˙k=−∂H/∂qk\dot p_k=-\partial H/\partial q_k then gives p˙k=0\dot p_k=0 immediately: pkp_k is a constant of the motion, fixed by the initial conditions.
  • The Hamiltonian reduction is complete: the conserved pk=αp_k=\alpha enters the reduced Hamiltonian as a parameter, the active phase space shrinks from 2n2n to 2(n−1)2(n-1) dimensions, and qk(t)q_k(t) is recovered at the end by one quadrature.
  • On the Lagrangian side the reduction is not automatic — the cyclic velocity q˙k\dot q_k still appears in LL; the Routhian R=∑cycpjq˙j−LR=\sum_{\text{cyc}}p_j\dot q_j-L fixes this by Legendre-transforming only the cyclic variables.
  • Cyclicity depends on the choice of coordinates: central forces have no cyclic Cartesian coordinate but a cyclic polar angle. Choosing — or transforming to — coordinates adapted to the symmetry is the whole game.
  • Key idea: If ∂H/∂qk=0\partial H/\partial q_k = 0 then pkp_k is conserved; the cyclic coordinate qkq_k drops out of the active dynamics and is recovered at the end by one integration.

Where this lives today

The purest industrial descendant of the cyclic coordinate is guiding-center and gyrokinetic theory in magnetized plasmas. A charged particle in a tokamak gyrates around a field line millions of times per transit; the gyrophase is (to excellent approximation) an ignorable angle, and its conjugate momentum — the magnetic moment μ\mu — is the frozen quantity that makes the simulation tractable. Gyrokinetic codes used to design fusion experiments like ITER are, at bottom, Routh's maneuver executed asymptotically: eliminate the fast cyclic angle, keep its conserved momentum as a parameter, and evolve only the slow remaining dynamics — a million-fold saving in time resolution. The same reduction, in its modern geometric dress (symplectic reduction), runs through spacecraft attitude control, the falling-cat problem and robot locomotion — where the "cyclic" directions are overall rotations whose conserved angular momentum constrains what shape-changes can achieve — and through molecular-dynamics codes, which exploit conserved total momentum to factor out the center-of-mass drift of every simulation box. Whenever a simulation "works in the co-moving/co-rotating frame with fixed angular momentum," a cyclic coordinate has been ignored, in Routh's precise sense of the word.

The physical picture

Think of a cyclic coordinate as a dial you can turn without changing the machine's internal state. Nothing in the energy bookkeeping notices the dial's setting, so whatever momentum is associated with turning it has nothing to push against — it stays frozen at its initial value. A frictionless puck on an infinite featureless table has two such dials (xx and yy: linear momentum conserved); a spinning axisymmetric top has its spin angle (pψp_\psi conserved); a planet orbiting the Sun has the polar angle (pθ=ℓp_\theta=\ell conserved). The crucial subtlety is that the same physical system shows or hides its dials depending on the coordinates you draw. That is the hero picture:

Figure
One orbit, two coordinate charts — Cartesian coordinates hide the conservation law that polar coordinates hand you by inspection
Cartesian chart (x, y) H depends on x AND y no coordinate is cyclic — nothing frozen by inspection polar chart (r, θ) r θ H = p²ᵣ/2m + p²θ/2mr² + V(r) θ never appears — p_θ frozen at its initial value same orbit, same physics — the conservation law is only VISIBLE in the chart adapted to the symmetry a cyclic coordinate is a property of the chart, not of the system

Geometrically, each conserved pkp_k traps the trajectory on a level sheet pk=αp_k=\alpha inside the 2n2n-dimensional phase space. The motion never leaves that sheet, so the active dynamics — the part you must actually integrate — lives in two fewer dimensions. Stack up one cyclic coordinate after another and the problem cascades downward:

Figure
The reduction cascade — every cyclic pair freezes onto a sheet of phase space, leaving fewer active equations and one quadrature IOU
full phase space: 2n dimensions sheets p_k = const slice the space; the trajectory is trapped on ONE sheet p_k = α reduced problem H(q₁…q_n₋₁, p₁…p_n₋₁; α) 2(n−1) active dimensions α enters as a fixed parameter, like a mass or a spring constant the quadrature IOU, paid last q_k(t) = q_k(0) + ∫ (∂H/∂p_k) dt′

The bookkeeping is hierarchical, and worth internalizing as a recipe: solve the reduced problem for the non-cyclic variables first, treating each conserved momentum as a number; only then integrate q˙k=∂H/∂pk\dot q_k=\partial H/\partial p_k — now an explicit function of time — to recover where the ignorable coordinate actually went. The ignorable coordinate is not unphysical; it is simply deferrable.

Cyclic coordinates in the Hamiltonian language

Everything starts from a definition you can check by inspection. Given a Hamiltonian H(q,p,t)H(q,p,t) — built, as always, by the Legendre transformation from LL — scan its argument list:

∂H∂qk=0  ⟺  qk is cyclic(1)\frac{\partial H}{\partial q_k}=0 \;\Longleftrightarrow\; q_k\ \text{is cyclic} \qquad (1)

where: HH — the Hamiltonian H(q,p,t)H(q,p,t); qkq_k — the coordinate under test; pkp_k — its conjugate momentum, which may (and usually does) still appear in HH. Read it physically: the energy landscape is flat along the qkq_k direction. Note what is not required — pkp_k absent, q˙k\dot q_k absent, or HH time-independent. Only the coordinate itself must be missing. The conservation law is then a one-line consequence of Hamilton's canonical equations (see Hamilton's Canonical Equations — we use them here, not re-derive them). The momentum equation p˙k=−∂H/∂qk\dot p_k=-\partial H/\partial q_k collapses:

p˙k=−∂H∂qk=0    ⟹    pk=α=const(2)\dot p_k=-\frac{\partial H}{\partial q_k}=0 \;\;\Longrightarrow\;\; p_k=\alpha=\text{const} \qquad (2)

where: α\alpha — the frozen value of pkp_k, set once by the initial conditions. This is the entire proof — no clever integrating factor, no Noether machinery. The Hamiltonian formalism was built so that conservation laws would read off like this: one flat direction in HH, one frozen momentum. Whether that frozen number is a linear momentum, an angular momentum, or something less familiar (the canonical pϕp_\phi of a charged particle includes a field term) depends on the coordinate, not on the theorem.

Derivation ACyclic implies conserved, and the problem genuinely drops a dimension

We are proving the two claims that make cyclic coordinates the workhorse of Hamiltonian problem-solving: that pkp_k is conserved the moment qkq_k is absent from HH, and — the part the Lagrangian formulation cannot match — that the cyclic pair (qk,pk)(q_k,p_k) then leaves the active dynamics entirely, reducing an nn-degree-of-freedom problem to n−1n-1 degrees of freedom plus one final integration. The strategy is to read Hamilton's equations with the flat direction imposed, check that cyclicity is a property of the physics and not of the formalism by carrying it through the Legendre transform, and then count what is left. Watch the contrast with the Lagrangian side: there the conserved quantity exists but the cyclic velocity q˙k\dot q_k still contaminates the remaining equations.

Step 1. Impose the flat direction on the momentum equation. For every index, Hamilton's equations read q˙i=∂H/∂pi\dot q_i=\partial H/\partial p_i and p˙i=−∂H/∂qi\dot p_i=-\partial H/\partial q_i. If ∂H/∂qk=0\partial H/\partial q_k=0, then p˙k=0\dot p_k=0, so pk(t)=pk(0)≡αp_k(t)=p_k(0)\equiv\alpha for all time.

Step 2. Check that cyclicity transfers through the Legendre transformation. Differentiating H=∑jpjq˙j−LH=\sum_j p_j\dot q_j-L at fixed (p,t)(p,t), with q˙j=q˙j(q,p,t)\dot q_j=\dot q_j(q,p,t) implicit, the chain rule gives

∂H∂qk=∑jpj∂q˙j∂qk−∂L∂qk−∑j∂L∂q˙j∂q˙j∂qk=−∂L∂qk\frac{\partial H}{\partial q_k}=\sum_j p_j\frac{\partial \dot q_j}{\partial q_k}-\frac{\partial L}{\partial q_k}-\sum_j\frac{\partial L}{\partial \dot q_j}\frac{\partial \dot q_j}{\partial q_k}=-\frac{\partial L}{\partial q_k}

because pj=∂L/∂q˙jp_j=\partial L/\partial\dot q_j makes the two sums cancel identically. Why this matters: qkq_k is cyclic in HH exactly when it is cyclic in LL — the property belongs to the system-plus-chart, not to which formalism you opened. The conserved quantity here is the same one the Lagrangian sibling article finds via d(∂L/∂q˙k)/dt=0\mathrm{d}(\partial L/\partial\dot q_k)/\mathrm{d}t=0.

Step 3. Substitute pk=αp_k=\alpha into the Hamiltonian. Since qkq_k never appears and pkp_k appears only as the constant α\alpha, the function

Hred=H(q1,…,qn−1, p1,…,pn−1; α, t)H_{\text{red}}=H(q_1,\dots,q_{n-1},\,p_1,\dots,p_{n-1};\,\alpha,\,t)

generates the complete dynamics of the remaining n−1n-1 pairs: their 2(n−1)2(n-1) canonical equations close among themselves, with α\alpha playing the role of a fixed parameter like a mass or a coupling constant. Why the Lagrangian cannot do this directly: setting ∂L/∂qk=0\partial L/\partial q_k=0 conserves pk=∂L/∂q˙kp_k=\partial L/\partial\dot q_k, but this is a relation among the velocities — q˙k\dot q_k itself still appears in LL and hence in every remaining Euler–Lagrange equation. To exploit the constant you must solve for q˙k\dot q_k and eliminate it by hand, and doing that carelessly inside LL produces wrong equations of motion. The cure is Routh's construction, derived next.

Step 4. Recover the ignorable coordinate by quadrature. Along the solved reduced trajectory, q˙k=∂H/∂pk∣pk=α\dot q_k=\partial H/\partial p_k\big|_{p_k=\alpha} is a known function of time, so

qk(t)=qk(0)+∫0t∂H∂pk dt′q_k(t)=q_k(0)+\int_0^t \frac{\partial H}{\partial p_k}\,\mathrm{d}t'

is a single integration — no differential equation remains for qkq_k.

Which statement when

you needreach fornotes
spot a conserved momentum by inspectioneq. (1)scan HH for missing coordinates
the conservation law itselfeq. (2)pk=αp_k=\alpha, one line from Hamilton's equations
the reduced problemHred(q,p;α,t)H_{\text{red}}(q,p;\alpha,t)cyclic pair gone; α\alpha a parameter
the ignorable coordinate's historyquadrature, Step 4 abovedo it last
reduction on the Lagrangian sideRouthian, eq. (3)half Legendre transform
effective 1D dynamics for central forcesDerivation CVeff=ℓ2/2mr2+V(r)V_{\text{eff}}=\ell^2/2mr^2+V(r)

The Routhian: reduce like a Hamiltonian, solve like a Lagrangian

Sometimes you want the cyclic coordinates eliminated and the surviving equations in familiar Euler–Lagrange form — for instance to reuse small-oscillation machinery about a steady motion. Routh's construction does exactly this: perform the Legendre transformation only over the cyclic variables, leaving the honest coordinates in Lagrangian dress. Suppose qs+1,…,qnq_{s+1},\dots,q_n are cyclic. Define

R(q1,…,qs, q˙1,…,q˙s, ps+1,…,pn, t)  =  ∑j=s+1npj q˙j  −  L(3)R(q_1,\dots,q_s,\,\dot q_1,\dots,\dot q_s,\,p_{s+1},\dots,p_n,\,t) \;=\; \sum_{j=s+1}^{n} p_j\,\dot q_j \;-\; L \qquad (3)

where: RR — the Routhian, a hybrid of HH and −L-L; the sum runs over the cyclic coordinates only, whose velocities are eliminated in favor of the conserved momenta pj=βjp_j=\beta_j; the non-cyclic q˙i\dot q_i survive untouched. Think of RR as a partial Legendre transform (the full-transform machinery is in Legendre Transformation): Hamiltonian where that pays (the frozen directions), Lagrangian where that is comfortable (the live ones). The dynamics generated by RR splits cleanly along the same line:

ddt∂R∂q˙i−∂R∂qi=0    (i≤s),q˙j=∂R∂pj,    p˙j=−∂R∂qj=0    (j>s)(4)\frac{\mathrm{d}}{\mathrm{d}t}\frac{\partial R}{\partial \dot q_i}-\frac{\partial R}{\partial q_i}=0 \;\;(i\le s), \qquad \dot q_j=\frac{\partial R}{\partial p_j},\;\; \dot p_j=-\frac{\partial R}{\partial q_j}=0 \;\;(j>s) \qquad (4)

where: the first set are Euler–Lagrange equations with RR standing in for the Lagrangian of the reduced system; the second set are Hamilton's equations for the cyclic pairs, the momentum ones confirming pj=βjp_j=\beta_j. In practice you use only the first set — with every βj\beta_j inserted as a number — and pay the quadrature IOU for the qjq_j at the end.

Go deeper · Derivation B · 1 minRouth's hybrid function and its split equations of motion
Derivation BRouth's hybrid function and its split equations of motion

The goal is to justify equation (4): that after the partial Legendre transform (3), the non-cyclic coordinates obey Euler–Lagrange equations generated by RR, while the cyclic pairs obey Hamilton's equations — so the conserved momenta can be legally substituted as constants before deriving the reduced equations of motion. The strategy is the standard differential bookkeeping of any Legendre transform: expand dR\mathrm{d}R from the definition, cancel the conjugate-pair terms, and read off partial derivatives. The point to watch is which variables RR actually depends on — that is the entire content of the construction.

Step 1. Take the total differential of definition (3), using dL=∑i(∂L∂qidqi+∂L∂q˙idq˙i)+∂L∂tdt\mathrm{d}L=\sum_i\big(\tfrac{\partial L}{\partial q_i}\mathrm{d}q_i+\tfrac{\partial L}{\partial\dot q_i}\mathrm{d}\dot q_i\big)+\tfrac{\partial L}{\partial t}\mathrm{d}t over all nn coordinates:

dR=∑j>s(q˙j dpj+pj dq˙j)−∑i∂L∂qi dqi−∑i∂L∂q˙i dq˙i−∂L∂t dt\mathrm{d}R=\sum_{j>s}\big(\dot q_j\,\mathrm{d}p_j+p_j\,\mathrm{d}\dot q_j\big)-\sum_{i}\frac{\partial L}{\partial q_i}\,\mathrm{d}q_i-\sum_{i}\frac{\partial L}{\partial \dot q_i}\,\mathrm{d}\dot q_i-\frac{\partial L}{\partial t}\,\mathrm{d}t

Step 2. Cancel the cyclic velocity differentials. For j>sj>s, the definition of momentum gives ∂L/∂q˙j=pj\partial L/\partial\dot q_j=p_j, so the +pj dq˙j+p_j\,\mathrm{d}\dot q_j and −pj dq˙j-p_j\,\mathrm{d}\dot q_j terms annihilate. Moreover ∂L/∂qj=0\partial L/\partial q_j=0 because those coordinates are cyclic. What survives is

dR=∑j>sq˙j dpj−∑i≤s(∂L∂qi dqi+∂L∂q˙i dq˙i)−∂L∂t dt\mathrm{d}R=\sum_{j>s}\dot q_j\,\mathrm{d}p_j-\sum_{i\le s}\Big(\frac{\partial L}{\partial q_i}\,\mathrm{d}q_i+\frac{\partial L}{\partial \dot q_i}\,\mathrm{d}\dot q_i\Big)-\frac{\partial L}{\partial t}\,\mathrm{d}t

Why the cancellation is the construction: no dq˙j\mathrm{d}\dot q_j survives for the cyclic directions, so RR genuinely depends on (qi,q˙i)(q_i,\dot q_i) for i≤si\le s and on pjp_j for j>sj>s — precisely the hybrid argument list claimed in (3).

Step 3. Read off coefficients: ∂R/∂pj=q˙j\partial R/\partial p_j=\dot q_j and ∂R/∂qj=0=−p˙j\partial R/\partial q_j=0=-\dot p_j for the cyclic pairs (Hamilton's equations, conservation included), and ∂R/∂qi=−∂L/∂qi\partial R/\partial q_i=-\partial L/\partial q_i, ∂R/∂q˙i=−∂L/∂q˙i\partial R/\partial\dot q_i=-\partial L/\partial\dot q_i for the live coordinates. Substituting the last two into the Euler–Lagrange equation for qiq_i flips both signs at once, leaving its form intact with RR in place of LL — which is (4). Why −R-R is the honest effective Lagrangian: since both partials of RR carry a minus sign relative to LL, the function Leff=−RL_{\text{eff}}=-R has the standard signs. For the central force this is where the centrifugal term famously flips: the kinetic piece ℓ22mr2\tfrac{\ell^2}{2mr^2} of LL reappears inside −R-R as a repulsive potential +ℓ22mr2+\tfrac{\ell^2}{2mr^2} in VeffV_{\text{eff}} — forget the flip and steady radii come out wrong.

Reduction in practice: the central force, start to finish

The planar central-force problem is the canonical demonstration — two degrees of freedom in, one out, and every step of the cascade visible. It is also the computation behind Kepler's second law and behind Worked Example 1 below.

Go deeper · Derivation C · 1 minCentral force reduced to one radial degree of freedom by the cyclic angle
Derivation CCentral force reduced to one radial degree of freedom by the cyclic angle

We are executing the full reduction pipeline on a particle of mass mm in a central potential V(r)V(r): identify the cyclic coordinate, freeze its momentum, build the reduced one-dimensional Hamiltonian with an effective potential, and end with the quadrature that recovers the angle. The strategy is simply Derivation A run on a concrete system, so each abstract step acquires a face. Watch how the conserved ℓ\ell migrates from "momentum" to "parameter": the term pθ2/2mr2p_\theta^2/2mr^2, kinetic in origin, becomes a repulsive ℓ2/2mr2\ell^2/2mr^2 barrier in the radial problem.

Step 1. Lagrangian and momenta in plane polar coordinates. From L=12m(r˙2+r2θ˙2)−V(r)L=\tfrac12 m(\dot r^2+r^2\dot\theta^2)-V(r), the conjugate momenta are pr=mr˙p_r=m\dot r and pθ=mr2θ˙p_\theta=mr^2\dot\theta — the angular momentum about the force center.

Step 2. Build the Hamiltonian by the (full) Legendre transform and scan it:

H=pr22m+pθ22mr2+V(r)H=\frac{p_r^2}{2m}+\frac{p_\theta^2}{2mr^2}+V(r)

The angle θ\theta is nowhere in sight: it is cyclic, by inspection, exactly as promised by the hero figure.

Step 3. Freeze the conserved momentum: p˙θ=−∂H/∂θ=0\dot p_\theta=-\partial H/\partial\theta=0, so pθ=ℓp_\theta=\ell, a number fixed by the launch conditions. Why this is rotational symmetry speaking: θ\theta absent from HH is the statement that the physics is invariant under rotations about the center — the symmetry-side of this coin is told properly in the Lagrangian sibling article and in Noether's theorem, and we do not retell it here.

Step 4. Read off the reduced one-dimensional problem. With ℓ\ell inserted, the radial pair obeys the Hamiltonian

Hred=pr22m+ℓ22mr2+V(r)⏟Veff(r)H_{\text{red}}=\frac{p_r^2}{2m}+\underbrace{\frac{\ell^2}{2mr^2}+V(r)}_{V_{\text{eff}}(r)}

a bona fide 1D system: r˙=pr/m\dot r=p_r/m, p˙r=−Veff′(r)\dot p_r=-V_{\text{eff}}'(r). Turning points, stability of circular orbits, oscillation frequencies — all now read off a one-dimensional potential. Why the barrier term is physical: ℓ2/2mr2\ell^2/2mr^2 is the kinetic energy the particle is obliged to keep because angular momentum cannot leak away; shrinking rr at fixed ℓ\ell forces faster rotation. It is not a new force — it is the old conservation law wearing a potential's clothes.

Step 5. Pay the quadrature IOU. Once r(t)r(t) is solved, the angle follows from one integration:

θ(t)=θ0+∫0tℓm r(t′)2 dt′\theta(t)=\theta_0+\int_0^t \frac{\ell}{m\,r(t')^2}\,\mathrm{d}t'

and the areal-velocity statement 12r2θ˙=ℓ/2m=const\tfrac12 r^2\dot\theta=\ell/2m=\text{const} is Kepler's second law, free of charge.

The canonical-transformation dream

Here is the thought that turns a bookkeeping trick into a research program. Cyclicity is chart-dependent: the central force has zero cyclic Cartesian coordinates but one cyclic polar coordinate. Choosing better coordinates created a conserved momentum out of thin air. Push the greed to its logical end: could there be coordinates in which every coordinate is cyclic? Then H=H(P1,…,Pn)H=H(P_1,\dots,P_n) alone, every PiP_i is constant, every QiQ_i evolves linearly in time — the system is solved by pure algebra:

H=H(P)    ⟹    P˙i=0,Q˙i=∂H∂Pi=const(5)H=H(P) \;\;\Longrightarrow\;\; \dot P_i=0,\qquad \dot Q_i=\frac{\partial H}{\partial P_i}=\text{const} \qquad (5)

where: (Q,P)(Q,P) — a hypothetical set of phase-space coordinates adapted so perfectly to the dynamics that nothing accelerates. Reaching them requires transformations that mix coordinates and momenta wholesale — far beyond the point transformations Q=Q(q,t)Q=Q(q,t) of Lagrangian mechanics — while preserving the canonical form of Hamilton's equations. Which changes of variables are allowed, and how to build them, is precisely the subject of the next article, Definition of Canonical Transformations; that the goal is achievable for a large class of systems is the triumph of Hamilton–Jacobi theory and action-angle variables.

Figure
From zero cyclic coordinates to all of them — the road that motivates the canonical-transformations chapter
generic chart 0 cyclic coordinates 2n coupled ODEs conservation laws hidden choose wisely adapted chart some coordinates cyclic reduced problem + quadratures this article canonical transformation perfect chart ALL coordinates cyclic H = H(P) only solved by algebra mixing q and p is essential — point transformations of q alone cannot get there the rules of the road: Definition of Canonical Transformations, then Hamilton–Jacobi theory reduction by cyclic coordinates, taken to its logical conclusion

Limiting cases

CaseResultMeaning
free particle, Cartesian chartH=(px2+py2+pz2)/2mH=(p_x^2+p_y^2+p_z^2)/2mall three coordinates cyclic; full momentum vector frozen — maximal reduction, motion solved by inspection
central force, polar chartH=pr2/2m+pθ2/2mr2+V(r)H=p_r^2/2m+p_\theta^2/2mr^2+V(r)one cyclic angle; pθ=ℓp_\theta=\ell frozen, radial 1D problem with VeffV_{\text{eff}} remains
central force, Cartesian chartno coordinate cyclicsame conserved ℓ=xpy−ypx\ell=xp_y-yp_x exists but is invisible to inspection — cyclicity is chart-dependent
heavy symmetric topϕ,ψ\phi,\psi cyclic in H(θ,pθ,pϕ,pψ)H(\theta,p_\theta,p_\phi,p_\psi)two frozen momenta; only nutation θ\theta is dynamical — a 3-DOF problem collapses to 1
HH explicitly time-dependent, qkq_k absentpkp_k still conservedcyclic conservation survives when energy conservation dies
all coordinates cyclicH=H(P)H=H(P): P˙=0\dot P=0, Q˙=const\dot Q=\text{const}the canonical-transformation dream (5); realized by Hamilton–Jacobi / action-angle variables

Worked examples

1 · Satellite in Earth's gravity: freeze the angle, read the orbit (confidence builder)

A satellite of mass m=800 kgm=800\ \mathrm{kg} is at r0=7000 kmr_0=7000\ \mathrm{km} from Earth's center, moving purely tangentially at v0=7.50 km/sv_0=7.50\ \mathrm{km/s}. Take GM⊕=3.986×1014 m3/s2GM_\oplus=3.986\times10^{14}\ \mathrm{m^3/s^2}. Find the conserved pθp_\theta, the effective potential picture at r0r_0, and the other apsis of the orbit.

Step 1. Identify the cyclic coordinate and freeze its momentum. In the plane of motion, H=pr2/2m+pθ2/2mr2−GMm/rH=p_r^2/2m+p_\theta^2/2mr^2-GMm/r has θ\theta cyclic, so pθ=ℓ=mr0v0=800×(7.0×106)×7500=4.20×1013 kg m2/sp_\theta=\ell=mr_0v_0=800\times(7.0\times10^{6})\times7500=\mathbf{4.20\times10^{13}\ kg\,m^2/s} — locked for the whole flight.

Step 2. Evaluate the reduced 1D energy. The launch is tangential, so pr=0p_r=0 and r0r_0 is an apsis. Centrifugal term: ℓ2/2mr02=12mv02=12×800×(7500)2=2.25×1010 J\ell^2/2mr_0^2=\tfrac12 mv_0^2=\tfrac12\times800\times(7500)^2=2.25\times10^{10}\ \mathrm{J}. Gravity: −GMm/r0=−3.1888×1017/7.0×106=−4.555×1010 J-GMm/r_0=-3.1888\times10^{17}/7.0\times10^{6}=-4.555\times10^{10}\ \mathrm{J}. Total: E=Veff(r0)=−2.305×1010 J<0E=V_{\text{eff}}(r_0)=-2.305\times10^{10}\ \mathrm{J}<0 — bound orbit.

Step 3. Which way does rr move first? Veff′(r0)=−ℓ2/mr03+GMm/r02=−6429+6508=+79 NV_{\text{eff}}'(r_0)=-\ell^2/mr_0^3+GMm/r_0^2=-6429+6508=+79\ \mathrm{N}, so the radial force −Veff′=−79 N-V_{\text{eff}}'=-79\ \mathrm{N} points inward: the satellite falls toward Earth initially. (Cross-check: circular speed at r0r_0 is GM/r0=7.546 km/s\sqrt{GM/r_0}=7.546\ \mathrm{km/s}; we launched slower, so r0r_0 must be the apogee.)

Step 4. Find the perigee from the frozen quantities alone. Semi-major axis from the vis-viva energy: a=−GMm/2E=3.1888×1017/(4.611×1010)=6.916×106 ma=-GMm/2E=3.1888\times10^{17}/(4.611\times10^{10})=6.916\times10^{6}\ \mathrm{m}, so rp=2a−r0=6.832×106 mr_p=2a-r_0=6.832\times10^{6}\ \mathrm{m} (perigee altitude ≈454 km\approx 454\ \mathrm{km}). Perigee speed from the conserved pθp_\theta: vp=ℓ/mrp=4.20×1013/(800×6.832×106)=7.685 km/sv_p=\ell/mr_p=4.20\times10^{13}/(800\times6.832\times10^{6})=\mathbf{7.685\ km/s}. Energy check there: 12mvp2−GMm/rp=2.362×1010−4.667×1010=−2.305×1010 J\tfrac12mv_p^2-GMm/r_p=2.362\times10^{10}-4.667\times10^{10}=-2.305\times10^{10}\ \mathrm{J} ✓ — same EE, as conservation demands.

2 · Spherical pendulum: one cyclic angle, one live angle (exam level)

A pendulum bob of mass m=0.20 kgm=0.20\ \mathrm{kg} swings on a rigid massless rod of length l=1.0 ml=1.0\ \mathrm{m}. At t=0t=0 it is at polar angle θ0=60∘\theta_0=60^\circ from the downward vertical with θ˙0=0\dot\theta_0=0 and azimuthal rate ϕ˙0=3.0 rad/s\dot\phi_0=3.0\ \mathrm{rad/s}. Take g=9.8 m/s2g=9.8\ \mathrm{m/s^2}. Find the conserved momentum, decide whether the bob initially rises or falls, and compute the initial θ¨\ddot\theta.

Step 1. Scan for the cyclic coordinate. With L=12ml2(θ˙2+sin⁡2 ⁣θ ϕ˙2)+mglcos⁡θL=\tfrac12ml^2(\dot\theta^2+\sin^2\!\theta\,\dot\phi^2)+mgl\cos\theta, the azimuth ϕ\phi is absent: cyclic. Its frozen momentum is pϕ=ml2sin⁡2 ⁣θ ϕ˙=0.20×1.0×(0.8660)2×3.0=0.20×0.75×3.0=0.450 kg m2/sp_\phi=ml^2\sin^2\!\theta\,\dot\phi=0.20\times1.0\times(0.8660)^2\times3.0=0.20\times0.75\times3.0=\mathbf{0.450\ kg\,m^2/s}.

Step 2. Build the reduced problem for the live angle. Eliminating ϕ˙\dot\phi via the Routhian (Leff=−RL_{\text{eff}}=-R) gives a 1D system in θ\theta with Ueff(θ)=pϕ2/(2ml2sin⁡2 ⁣θ)−mglcos⁡θU_{\text{eff}}(\theta)=p_\phi^2/(2ml^2\sin^2\!\theta)-mgl\cos\theta — note the sign flip that turns the azimuthal kinetic term into a repulsive barrier around the vertical.

Step 3. Compare against steady (conical) motion. A horizontal circle at θ0=60∘\theta_0=60^\circ would need ϕ˙c=g/(lcos⁡θ0)=9.8/0.50=4.43 rad/s\dot\phi_c=\sqrt{g/(l\cos\theta_0)}=\sqrt{9.8/0.50}=4.43\ \mathrm{rad/s}. Our bob spins at only 3.0 rad/s3.0\ \mathrm{rad/s} — too slow for the centrifugal barrier to hold it out at 60∘60^\circ, so gravity wins and the bob initially falls toward the vertical (θ\theta decreases... careful: θ\theta is measured from the downward vertical, so falling toward the vertical means θ\theta decreases).

Step 4. Compute the initial angular acceleration from ml2θ¨=−Ueff′(θ)ml^2\ddot\theta=-U_{\text{eff}}'(\theta). The barrier term: pϕ2cos⁡θ0/(ml2sin⁡3 ⁣θ0)=(0.450)2×0.50/(0.20×0.6495)=0.1013/0.1299=0.779 N mp_\phi^2\cos\theta_0/(ml^2\sin^3\!\theta_0)=(0.450)^2\times0.50/(0.20\times0.6495)=0.1013/0.1299=0.779\ \mathrm{N\,m}. The gravity term: mglsin⁡θ0=0.20×9.8×0.8660=1.697 N mmgl\sin\theta_0=0.20\times9.8\times0.8660=1.697\ \mathrm{N\,m}. So ml2θ¨=0.779−1.697=−0.918 N mml^2\ddot\theta=0.779-1.697=-0.918\ \mathrm{N\,m} and θ¨0=−0.918/0.20=−4.59 rad/s2\ddot\theta_0=-0.918/0.20=\mathbf{-4.59\ rad/s^2} — consistent with Step 3's verdict. The motion oscillates in θ\theta between two turning circles while ϕ\phi advances by the quadrature ϕ˙=pϕ/ml2sin⁡2 ⁣θ\dot\phi=p_\phi/ml^2\sin^2\!\theta, fastest when the bob dips closest to the vertical — the same 1/ ⁣sin⁡2θ1/\!\sin^2\theta speedup a Kepler orbit shows at perihelion via 1/r21/r^2.

Origins · 1 minEdward Routh and how this was discovered
ER
Edward Routh
1831–1907 · Cambridge

Edward Routh was Victorian Cambridge's supreme coach: over three decades he trained more than 600 students for the Mathematical Tripos, among them 27 Senior Wranglers — J. J. Thomson, Rayleigh, and Larmor passed through his drilling. Beating the young James Clerk Maxwell to Senior Wrangler in 1854 made his name; the mathematics that kept it was his 1877 Adams Prize essay on the stability of motion, which introduced both the Routh–Hurwitz stability criterion of modern control theory and the "ignoration of coordinates" — the Routhian — that eliminates cyclic variables before a stability analysis. That combination was no accident: to decide whether a spinning system (a governor, a top, a bicycle) is stable, one first strips away the uniform spinning via the cyclic coordinates, then examines small oscillations of what remains. Routh's procedure is still how one studies steady motions, from satellite attitude dynamics to the gaits of walking robots.

History
1788
— Lagrange's Mécanique analytique shows coordinates absent from the dynamics yield first integrals
1834
— Hamilton's canonical equations make cyclic-implies-conserved a one-line inspection
1877
— Routh's Adams Prize essay introduces the Routhian and the ignoration of coordinates for stability analysis
1904
— Whittaker's Analytical Dynamics systematizes reduction by ignorable coordinates and quadratures
1918
— Noether's theorem generalizes the coordinate-shift symmetry behind cyclicity to arbitrary continuous symmetries
1974
— Marsden and Weinstein recast the cyclic-pair dimension drop as symplectic reduction, the modern geometric form

Going further

References. Goldstein, Poole & Safko, Classical Mechanics, 3rd ed., §2.6, §8.2–8.3 · Landau & Lifshitz, Mechanics, 3rd ed., §14, §40–41 · Morin, Introduction to Classical Mechanics, Ch. 6 · Taylor, Classical Mechanics, §7.8 and §13.4 · Whittaker, Analytical Dynamics, Ch. III (ignoration of coordinates) The key identities used throughout are: the Legendre-transform relation ∂H/∂qk=−∂L/∂qk\partial H/\partial q_k=-\partial L/\partial q_k (derived by the chain-rule argument in the derivation above); the central-force momenta pθ=mr2θ˙p_\theta=mr^2\dot\theta and effective potential ℓ2/2mr2+V(r)\ell^2/2mr^2+V(r); and the symmetric-top Hamiltonian with its two cyclic angles. The worked-example numbers are illustrative arithmetic computed directly from the displayed formulas: satellite (ℓ=4.20×1013 kg m2/s\ell=4.20\times10^{13}\ \mathrm{kg\,m^2/s}, E=−2.305×1010 JE=-2.305\times10^{10}\ \mathrm{J}, rp=6.832×106 mr_p=6.832\times10^{6}\ \mathrm{m}, vp=7.685 km/sv_p=7.685\ \mathrm{km/s}); spherical pendulum (pϕ=0.450 kg m2/sp_\phi=0.450\ \mathrm{kg\,m^2/s}, ϕ˙c=4.43 rad/s\dot\phi_c=4.43\ \mathrm{rad/s}, θ¨0=−4.59 rad/s2\ddot\theta_0=-4.59\ \mathrm{rad/s^2}).

Builds on: The Hamiltonian Function · Generalized Momenta · Hamilton's Canonical Equations · Legendre Transformation Leads to: Definition of Canonical Transformations · Phase Space · Hamiltonian Flow and the Phase Fluid · Equivalence of Lagrangian and Hamiltonian Formulations

If conservation of angular momentum is a physical fact, how can it depend on my coordinate choice whether a coordinate is cyclic?

The conservation law never depends on the chart — only its visibility does. In polar coordinates the constant is the momentum conjugate to a coordinate that happens to be absent from HH, so it is found by inspection; in Cartesian coordinates the very same constant exists as the combination xpy−ypxxp_y-yp_x, but no Cartesian coordinate is cyclic, so inspection finds nothing and you would need Poisson-bracket or Noether machinery to detect it. "Cyclic" is a statement about the pair (system, chart). The practical moral: invest in coordinates aligned with the expected symmetry before building HH — and the theoretical moral, developed from the next article onward, is that transforming to charts with more cyclic coordinates is a solution method in itself.

Why exactly is the Lagrangian reduction incomplete if the momentum is conserved either way?

Conserved is not the same as eliminated. On the Lagrangian side, ∂L/∂qk=0\partial L/\partial q_k=0 gives you the constant pk=∂L/∂q˙kp_k=\partial L/\partial\dot q_k — but this is an equation relating velocities, and q˙k\dot q_k itself still appears inside LL, hence inside every remaining Euler–Lagrange equation. To use the constant you must solve for q˙k\dot q_k and substitute — and substituting into LL directly is illegal (it changes the variational problem and yields wrong equations; the sign of the centrifugal term comes out backwards). The Hamiltonian never has this problem because velocities were already traded for momenta by the Legendre transformation: setting pk=αp_k=\alpha in HH is a legitimate substitution of a constant, and the remaining canonical equations close among themselves. The Routhian is precisely the device that makes the elimination legal in Lagrangian dress.

Is the Routhian a Hamiltonian or a Lagrangian? The signs keep confusing me.

It is deliberately both, split by coordinate type: Hamiltonian in the cyclic directions (which it eliminates), Lagrangian in the live ones (which it keeps). The sign subtlety is that R=∑cycpjq˙j−LR=\sum_{\text{cyc}}p_j\dot q_j-L has both ∂R/∂qi\partial R/\partial q_i and ∂R/∂q˙i\partial R/\partial\dot q_i equal to minus the corresponding derivatives of LL, so the Euler–Lagrange equation keeps its form with RR in place of LL — but the honest effective Lagrangian is Leff=−RL_{\text{eff}}=-R. The place this bites is the central force: the azimuthal kinetic energy +ℓ2/2mr2+\ell^2/2mr^2 in LL reappears in LeffL_{\text{eff}} as a potential term −ℓ2/2mr2-\ell^2/2mr^2 subtracted, i.e. a repulsive barrier +ℓ2/2mr2+\ell^2/2mr^2 inside UeffU_{\text{eff}}. If your steady-orbit radii come out wrong by a sign, this flip is almost always the culprit.

Does a cyclic coordinate still give a conserved momentum when H depends explicitly on time?

Yes. The momentum equation p˙k=−∂H/∂qk\dot p_k=-\partial H/\partial q_k involves only the qkq_k-dependence of HH; explicit time dependence is irrelevant to it. So a bead on a wire being cranked with a prescribed time-dependent profile can perfectly well conserve some pkp_k while its energy is being pumped. The two conservation questions — is qkq_k absent (conserves pkp_k), is tt absent (conserves HH) — are independent switches. This independence is also why cyclic-coordinate reduction remains useful for driven systems where energy methods fail.

How is this related to Noether's theorem — is a cyclic coordinate just a special case?

Yes, the simplest one. A coordinate absent from HH (equivalently LL) means the dynamics is invariant under rigid shifts qk→qk+ϵq_k\to q_k+\epsilon; the conserved Noether charge of that one-parameter symmetry is exactly pkp_k. Noether's theorem generalizes in two directions cyclic coordinates cannot reach: symmetries that mix several coordinates (a rotation in Cartesian variables, where no single coordinate shifts), and symmetries acting on time itself. The Lagrangian-side account — including the theorem and its charge formula — lives in Cyclic Coordinates and Conservation of Conjugate Momenta and the Noether article; this page's contribution is the Hamiltonian mechanism (p˙k=−∂H/∂qk\dot p_k=-\partial H/\partial q_k) and the reduction that the conserved charge finances.

Can every system be brought to a chart where all coordinates are cyclic?

No — and the boundary is one of the deepest facts of mechanics. Systems for which it is possible (given nn independent conserved quantities in involution) are called integrable in the Liouville sense: the central force, the heavy symmetric top, all small-oscillation problems. For them, Hamilton–Jacobi theory and action-angle variables construct the all-cyclic chart explicitly, and motion winds uniformly on invariant tori. Generic systems with n≥2n\ge2 degrees of freedom are not integrable — the three-body problem being the classic casualty — and no chart, however clever, makes all their coordinates cyclic; chaos lives in the obstruction. The dream of equation (5) is thus a precise diagnostic: how close a system can get to it measures how solvable it is.

The conserved p for a charged particle came out different from mv in my problem set. Did I make an error?

Probably not — you rediscovered that the canonical momentum, not the kinetic one, is what cyclicity conserves. For a particle in a magnetic field the Lagrangian contains q v⋅Aq\,\mathbf{v}\cdot\mathbf{A}, so the momentum conjugate to an ignorable azimuth is pϕ=mr2ϕ˙+q rAϕp_\phi=mr^2\dot\phi+q\,rA_\phi — kinetic angular momentum plus a field contribution. It is this combined object that stays constant in an axisymmetric field (the basis of magnetic-mirror and tokamak confinement), while mr2ϕ˙mr^2\dot\phi alone drifts. The rule: always derive pk=∂L/∂q˙kp_k=\partial L/\partial\dot q_k honestly before freezing it; whenever velocity-dependent potentials are in play, conserved canonical momenta and mvm\mathbf{v} part company.

When should I use the Routhian instead of just working with the reduced Hamiltonian?

Use the reduced Hamiltonian when you want phase-space tools: portraits in the (r,pr)(r,p_r) plane, Liouville-volume arguments (see Liouville's Theorem), or a launchpad for canonical transformations. Use the Routhian when the remaining analysis is naturally Lagrangian — above all, small oscillations about steady motion. Routh's own use case is still the standard one: a steady state with uniformly advancing cyclic coordinates (a precessing top, an orbiting satellite, a spinning governor) appears in the Routhian description as an equilibrium of the live coordinates, so the entire small-oscillation toolkit (expand LeffL_{\text{eff}} to second order, read off frequencies) applies verbatim. The two routes always agree — Derivations B and C produce the same radial equation — so the choice is purely about which toolbox you want open.

Before this

Where this leads