Codex @codex  0

An equatorial coordinate system labels a viewing direction on the celestial sphere by right ascension $\alpha$, increasing eastward around the reference celestial equator, and declination $\delta$, increasing northward from that celestial equator. These two angles locate a direction, without specifying the object's Euclidean distance. In the traditional convention, zero right ascension points toward the reference equinox, the northward crossing of the Sun's annual path through the celestial equator. An hour circle through the direction supplies the north-south arc on which declination is measured. ^[ref].

For a concrete construction, choose Cartesian coordinates with $x$ pointing toward the reference equinox, $z$ toward the north reference celestial pole, and positive $y$ given by the cross product of the positive $z$ and $x$ directions. The direction's unit vector isThis is a spherical coordinate system with polar angle $\pi /2-\delta$ and no radial coordinate. The trigonometric functions require consistent angular units: convert right ascension in hours to radians before evaluating them. The synthetic direction $\alpha =4^{\mathrm{h}}=60^{\circ }$, $\delta =30^{\circ }$ gives $\mathbf{n}=(\sqrt{3}/4,3/4,1/2)$; its squared norm is $3/16+9/16+1/4=1$. A change of reference axes applies a rotation to this unit vector before recovering the two angles. Away from the celestial poles, recover right ascension using the quadrant-aware $atan2(n_y,n_x)$ reduced modulo $2\pi$, and declination using $\arcsin n_z$; a plain arctangent of $n_y/n_x$ loses the quadrant. Adding a common offset to every right ascension describes only the special case of a rotation about the same polar axis.

The reference specification matters. For an equatorial coordinate system tied to the Earth, precession changes the reference celestial equator and reference equinox, and nutation distinguishes their mean and true versions. The date specifying those axes is distinct from the astronomical epoch at which a catalogue's source directions apply: advancing a star's proper motion and changing reference axes are separate operations. A modern catalogue may instead use fixed reference frame axes, aligned approximately with the traditional J2000.0 axes, without a moving reference equinox. Thus a pair of numbers labelled merely "J2000" does not settle every convention. State the reference frame, the astronomical epoch, and whether a direction is catalogued or apparent. ^[ref].

During one night's elementary fixed-direction approximation, a star's right ascension and declination remain constant while its altitude, azimuth and hour angle change because the Earth turns. In angular units, with $\Theta$ the local sidereal time, west-positive hour angle $H$ obeysAt upper meridian transit, $H=0$ and the local sidereal time equals the transiting direction's right ascension. Pair apparent sidereal time with apparent right ascension referred to the true celestial equator and reference equinox of date; mean sidereal time refers to the mean reference equinox. Raw fixed-frame catalogue coordinates require conversion before precision pointing. ^[ref].

This separation makes the equatorial coordinate system useful: a catalogue can identify a target independently of the observer's current astronomical horizon, while local sidereal time tells the observer which hour circle is on the celestial meridian. Precision observations also account for proper motion, stellar parallax, aberration of light and atmospheric refraction; the fixed-direction construction above omits those effects. A reference-axis rotation preserves angles between directions, whereas changing the observer or correcting atmospheric refraction is a different operation. ^[ref].

Now take zero torque and $\mathbf{\omega }(0)=(1,0,1)$. The axial angular velocity is ${\omega }_3=1$ and the complex transverse component obeys $\dot{z}=i\nu z$ with $z(0)=1$ and $\nu =(I_3-I_1)/I_1$. Hence the angular momentum components in the principal body frame areThe transverse projection therefore traces a circle with signed angular frequency $\nu$. Its magnitude is constant, $L=\sqrt{I_1^2+I_3^2}$.

In the inertial frame, zero torque means that $\mathbf{L}$ is fixed. Usingand the derivative of a body-fixed basis vector, the symmetry axis satisfiesThus there is uniform torque-free precession of a symmetric top about the fixed angular momentum, at a constant cone angle $\beta$ and space angular frequencyThe cone angle stays constant, with no nutation. The body-frame and space-frame frequencies of free symmetric-top precession differ: $\nu$ describes rotating components, whereas ${\Omega }_{\mathrm{space}}$ describes the symmetry axis in space. The formulas also cover $I_3=I_1$, where the angular momentum components in the principal body frame are constant.

Noether theorem associates a conserved quantity ${\sum }_i{p}_i\delta q_i$ to a continuous symmetry of the Lagrangian (with the usual boundary-term correction if the Lagrangian changes by a total derivative). The angles $\psi$ and $\varphi$ are cyclic coordinates, so for the heavy symmetric topTime-translation symmetry gives conserved energy, whose Hamiltonian expression isIt is quadratic in the momenta; $p_{\theta }=I_1\dot{\theta }$. Eliminating the cyclic velocities and putting $\mu =\cos \theta$, so ${\dot{\mu }}^2={\sin }^2\theta {\dot{\theta }}^2$, yieldsThe first term has degree three and the second degree two, displaying the cubic governing nutation.

A rapidly spinning gyroscope approximately has angular momentum $\mathbf{L}=I_3{\omega }_3\hat{\mathbf{e}}$ along its axle. With no torque in an inertial frame, this vector retains its direction. A transverse torque changes its direction according to $d\mathbf{L}/dt=\mathbf{\tau }$; it does not initially push the axle in the direction of the applied force.

For a fixed pivot, centre of mass $l\hat{\mathbf{e}}$ from it and gravity $-Mg\hat{\mathbf{z}}$, the torque is $Mgl\hat{\mathbf{z}}\times \hat{\mathbf{e}}$. In slow steady precession about the vertical, $d\mathbf{L}/dt=\Omega \hat{\mathbf{z}}\times \mathbf{L}$, givingThe factor $\sin \theta$ appears both in torque and in the radius traced by the angular-momentum tip, so it cancels for a nonvertical axle. Reversing spin reverses this precession. The full symmetric-top steady condition is $I_1\cos \theta {\Omega }^2-I_3{\omega }_3\Omega +Mgl=0$; the approximation neglects its quadratic term. Sudden release can also produce nutation, so large spin alone does not guarantee a steady sweep.

A wheel-on-pivot experiment can compare precession periods at different measured spin speeds and reverse the spin. Keep the pivot and centre of mass fixed; bearing friction, spin decay, an initial push and support motion can mimic departures from the predicted inverse-spin relation. A freely gimballed direction reference and a gravity-loaded precessing wheel test different torque conditions.^[ref]

For torque-free rotation, angular momentum is fixed in inertial space, but angular velocity need not point along it. In principal body axes,and cyclic permutations. A symmetric top with transverse moment $A$ and axial moment $C$ has constant ${\omega }_3$, while ${\omega }_1+i{\omega }_2$ rotates in body coordinates at signed frequency $(C-A){\omega }_3/A$. The symmetry axis precesses around the fixed angular momentum in space at rate $L/A$. These body-frame and space-frame rates describe different motions and need not coincide.

Poinsot construction expresses the same torque-free motion geometrically. The energy ellipsoid $\sum I_i{\omega }_i^2=2T$ is fixed in the body. Its normal at the angular-velocity point is parallel to $\mathbf{L}$, and $\mathbf{L}\cdot \mathbf{\omega }=2T$ places that point in a fixed plane perpendicular to $\mathbf{L}$. The ellipsoid rolls on this invariable plane without slipping: the instantaneous velocity of its contact point is proportional to $\mathbf{\omega }\times \mathbf{\omega }=0$. The contact curve is the polhode on the body ellipsoid and the herpolhode on the plane.^[ref]

For $I_1<I_2<I_3$, linearization about spin $\Omega$ on the third axis givesThe extreme-moment axes are stable; the intermediate one is unstable. The major-axis relaxation theorem adds weak internal dissipation while retaining negligible external torque. At fixed $L$,Energy loss favours rotation about the axis of largest principal moment. Here “major axis” means maximum inertia, not necessarily the body's longest geometric direction. A perfectly rigid torque-free body does not dissipatively settle; its energy is conserved. Distinguishing stability from relaxation also explains why minimum-inertia rotation can be stable in the ideal rigid model.

Forced precession occurs when an external torque changes $\mathbf{L}$, with $\dot{\mathbf{L}}=\mathbf{\tau }$. For steady precession, $\mathbf{\tau }={\mathbf{\Omega }}_p\times \mathbf{L}$. A rapidly spinning, pivoted heavy gyroscope has $L\simeq C{\omega }_3$ along its symmetry axis, givingThe torque is perpendicular to the nearly axial angular momentum, so it mainly changes its direction. For a steady inclined symmetric top, the more complete condition isIt exhibits slow and fast precession branches when real solutions exist. A generic release excites nutation as well: inclination oscillates instead of staying fixed. Axis constraints and bearing forces must be included when specifying which angular momentum and pivot are being used.

For slow gyroscopic precession, the gravitational torque turns the spin angular momentum: $Mg\ell \sin \theta =\Omega I_3\omega \sin \theta$, with $I_3=M{R}^2/2$. ThusThe inclination cancels at this order. More precisely, uniform precession of a heavy symmetric top satisfiesIts slow root gives $T_p\simeq 1.57\mathrm{s}$. About $1.6\mathrm{s}$ is the slow-precession answer; arbitrary initial conditions can also produce nutation rather than exact steady precession.

Inertia is a tensor. For a rigid body, distances between material points are fixed. Write each position relative to the centre of mass as $\mathbf{r}$, so its rotational velocity is $\mathbf{\omega }\times \mathbf{r}$. Integrating $\mathbf{r}\times (\mathbf{\omega }\times \mathbf{r})$ gives angular momentumThe rotational kinetic energy isThe inertia tensor is real and symmetric, so orthogonal principal axes diagonalize it. Angular momentum is parallel to angular velocity only along a principal axis, or within a subspace of equal principal moments. For translation plus rotation about the centre of mass, $T=M{v}_{\mathrm{CM}}^2/2+T_{\mathrm{rot}}$. Moving the reference to a parallel axis adds $M(d^{2}1-\mathbf{d}{\mathbf{d}}^{\mathsf{T}})$ to the tensor, the tensor form of the parallel-axis theorem.

Euler's equations. The external torque about the centre of mass obeys $\mathbf{\tau }=(d\mathbf{L}/dt{)}_{\mathrm{space}}$. In body axes, the transport rule givesThe principal moments are constant in that frame, henceThese are Euler's equations. With zero torque, spatial angular momentum and kinetic energy are conserved, although the body components generally vary. For distinct moments, small perturbations about rotation around the largest or smallest principal moment oscillate, whereas those about the intermediate moment grow. This is the intermediate-axis instability.

A lamina. Choose the lamina's plane as $z=0$. Then $I_{xz}=I_{yz}=0$, so the normal is a principal axis, andThis perpendicular-axis theorem applies to a planar mass distribution; it is not a general identity for a three-dimensional body. Rotate the in-plane axes to diagonalize the remaining symmetric two-by-two block. For a uniform rectangular lamina with sides $a,b$ and mass $M$, its central principal moments are $I_x=M{b}^2/12$, $I_y=M{a}^2/12$ and $I_z=M(a^2+b^2)/12$. If $a>b$, then $I_x<I_y<I_z$: torque-free rotation about the long in-plane axis or the normal is stable, while that about the short in-plane axis is unstable. Coincident moments in a square require the degenerate, rather than strict intermediate-axis, case. Under a constraint to rotate only about the normal, the motion reduces to the scalar equation ${\tau }_z=I_z{\dot{\omega }}_z$.

A symmetric top without torque. Let $I_1=I_2=I_{\perp }$ and $I_3=I_{\parallel }$. Euler's equations reduce toThus ${\omega }_1+i{\omega }_2$ rotates at signed body-frame frequency $\nu$, while the axial component stays fixed. The symmetry axis has constant angle to fixed spatial angular momentum. Indeed $\mathbf{L}=I_{\perp }\mathbf{\omega }+(I_{\parallel }-I_{\perp }){\omega }_3{\mathbf{e}}_3$, soIt precesses in space about $\mathbf{L}$ with angular speed $L/I_{\perp }$. This spatial precession and the body-frame frequency $\nu$ describe different motions and should not be identified.

A heavy top. If a symmetric top is pivoted a distance $\ell$ from its centre of mass, gravity gives torque ${\mathbf{r}}_{\mathrm{CM}}\times M\mathbf{g}$. Here the relevant moments are about the fixed pivot. The axial angular momentum and the vertical angular-momentum component are conserved, because gravity supplies neither an axial torque nor a torque about the vertical. Using tilt $\theta$, precession angle $\varphi$ and spin angle $\psi$, its Lagrangian isEliminating the two conserved cyclic momenta leaves one-dimensional motion in an effective tilt potential; bounded oscillation in tilt is nutation. Steady precession at $\dot{\varphi }=\Omega$ and fixed tilt satisfiesFor a rapidly spinning top on the slow branch, the first term is small andGravity then changes the direction of angular momentum rather than simply tipping the top straight down. The approximation fails near conditions where nutation or the fast-precession branch matters; the full equations retain those possibilities.

For torque-free precession of a symmetric top, take principal moments of inertia $I_1=I_2=I_{\perp }$ and $I_3=I_{\parallel }$, with no external torque. In an inertial frame moving with its centre of mass, its angular momentum $\mathbf{L}$ is constant in an inertial frame, and its rotational kinetic energy is conserved.

In the principal body axes, Euler equations for a torque-free rigid body reduce toTherefore ${\omega }_3$ is constant andThus the transverse angular velocity rotates uniformly in the body frame; its signed rate reverses between oblate and prolate inertia tensors. This body-frame wobble rate should be distinguished from the space-frame precession of the symmetry axis.

Let ${\mathbf{e}}_3$ be the body's symmetry axis. Sinceits angle $\theta$ to $\mathbf{L}$ is constant, with $\cos \theta =I_{\parallel }{\omega }_3/L$, andHence the symmetry axis describes a cone about fixed $\mathbf{L}$ at the space-frame precession rateThere is no nutation of this angle for torque-free motion of an exactly symmetric top. In the usual Euler angles, ${\omega }_3=\dot{\psi }+\dot{\varphi }\cos \theta$, so $\dot{\psi }=-{\Omega }_b$. Both rates describe the same motion in different axes. A freely tossed axisymmetric disc illustrates the effect; a top supported at a point under gravity is a different, externally torqued problem.

Directions. Choose $+x$ along the horizontal spindle from the pivot to the disk, and $+z$ vertically upward. With lever arm $\mathbf{r}=l\hat{\mathbf{x}}$ and weight ${\mathbf{F}}_g=-mg\hat{\mathbf{z}}$, the torque about the pivot isClockwise precession viewed from above has $\mathbf{\Omega }=-\Omega \hat{\mathbf{z}}$. The relation $\mathbf{\tau }=\mathbf{\Omega }\times \mathbf{L}$ therefore requires the spin angular momentum to point along $-x$, from the disk toward the pivot. The disk spins clockwise when viewed from its outer face looking along the spindle toward the pivot. Looking from the pivot reverses that apparent rotation. The right-hand rule fixes these signs unambiguously.

The drawing shows an instantaneous top view; the disk's angular momentum and the horizontal spindle turn clockwise together. The pivot force has zero torque about the pivot.

Rate of precession. For steady horizontal motion the spin contribution has magnitude $L=I\omega$ and turns through $d\varphi =\Omega dt$ in time $dt$. Its change is perpendicular to itself and has magnitude $dL=Ld\varphi$. HenceThis is the usual slow steady gyroscope precession relation. Any constant vertical angular momentum associated with steady precession has zero time derivative; the horizontal spin contribution supplies the changing component here. Transient nutation is excluded by the prescribed horizontal steady motion.

Clock design. Convert both rotation rates to angular frequencies:For the moment of inertia of a uniform disc, $I=m{R}^2/2$, so the mass cancels:

Spin friction. At fixed geometry, reducing $\omega$ increases $\Omega \propto 1/\omega$. The precession period, $T_p=2\pi /\Omega$, therefore decreases. This trend holds while slow steady precession remains a suitable description; it cannot be extrapolated to a stopped disk.

Counterweight. To oppose the disk's gravitational torque, place the counterweight on the other side of the pivot, at $-l_c\hat{\mathbf{x}}$. The net torque isA horizontal axle can remain fixed when this vanishes, givingFor the fixed spin direction already found, counterbalanced gyroscope precession has signed rateMoving the counterweight outward from the pivot first slows the clockwise motion, eliminates its driving torque at balance, and then reverses the steady motion to counterclockwise. A counterweight on the disk's side would add to the gravitational torque. These statements neglect axial spin moment of inertia of the counterweight and assume the stated horizontal motion, without a separately imposed precession at balance.

Eighteenth-century support for the Newtonian system accumulated through successful predictions, the resolution of apparent discrepancies and increasingly precise terrestrial tests. The relevant system included laws of motion, universal gravitation and a programme of mathematical celestial mechanics; accepting one element did not automatically settle every optical or metaphysical dispute.

Edmond Halley identified earlier comet appearances with the same orbit and used Newtonian gravitation to predict a return. Its recovery in 1758, with the subsequent perihelion in 1759 and perturbative corrections for planetary attractions, made an especially persuasive test of a forecast made long before the event. It also brought comets into the same quantitative dynamics as planets rather than treating them as exceptional signs.^[ref]

Geodetic expeditions compared meridian arcs at different latitudes. The observed polar flattening agreed with the qualitative Newtonian expectation for a rotating self-gravitating Earth, against competing prolate interpretations. The inference required both accurate surveying and a theory of how rotation, gravity and internal mass distribution affect shape; exact flattening could not be predicted from rotation alone without an internal-structure model.^[ref]

Lunar and planetary perturbation theory provided a sustained, more demanding test than unperturbed Keplerian ellipses. Discrepancies in the Moon's apsidal motion initially appeared troublesome, but improved approximations within inverse-square gravitation recovered the missing contribution. Agreement achieved by a better controlled calculation supported the theory while demonstrating why a failure of an early approximation was not automatically a failure of the force law.

James Bradley's annual stellar aberration furnished evidence for Earth's orbital motion and finite light propagation; later nutation measurements also connected careful positional astronomy to the dynamics of Earth's axis. Aberration supports the moving-Earth framework, though it is not by itself a unique test of an inverse-square force.^[ref]

Finally, Henry Cavendish's torsion-balance experiment in 1798 measured attraction between laboratory masses and inferred Earth's mean density. It extended gravitation from celestial fitting to a quantitatively controlled terrestrial attraction; describing the original aim as determining Earth's density is historically more accurate than saying he set out to measure a modern separately defined $G$. Together these tests supplied converging support across scales. They did not prove Newtonian mechanics exactly or establish every part of Newton's natural philosophy, but they made its gravitational research programme substantially more credible.

Model the assembly as a heavy symmetric top with a fixed frictionless support, a uniform gravitational field, and centre of mass on its symmetry axis. Assume constant spin, negligible damping, and initial conditions prepared for uniform precession of a heavy symmetric top. Large spin means the transverse angular momentum from precession is small: $I_1|\Omega |\sin \theta \ll |I_3{\omega }_3|$. Also require $|\Omega |\ll |{\omega }_3|$. A large spin rate alone does not select steady motion after an arbitrary release.

For an original illustrative calculation, stipulateThe transverse principal moments of inertia $I_1,I_2$ are about the pivot, including the parallel axis theorem contribution. These inputs giveThe ratio ${\Omega }_0/{\omega }_3=0.00368$ is small; the more direct angular momentum ratio $I_1{\Omega }_0\sin 60^{\circ }/L_s=0.00358$ is small too. These are model calculations, not experimental results.

Check the approximation against the existing exact uniform precession of a heavy symmetric top condition. Its total angular momentum isAt constant $\theta ,\Omega ,L_s$, differentiating this expression and equating it to the gravitational torque givesFor $L_s>0$ the slow solution can be evaluated without subtracting two nearly equal numbers:It exceeds ${\Omega }_0$ by about $0.208\%$. The other solution is about $177.409rad/s$ and violates the assumed slow precession; the quadratic's two roots must not be conflated. Its discriminant must be nonnegative for these steady solutions to exist.

The original curves hold $I_1,I_3,\theta$ fixed and vary spin and the imposed torque coefficient $K$ in ${\mathbf{\tau }}_O=K\hat{\mathbf{z}}\times \hat{\mathbf{e}}$. For gravity alone, $K=Mgd$; the other curves are hypothetical changes in the transverse torque, not observations of this wheel. Solid curves use the slow quadratic root; dashed curves use $K/(I_3{\omega }_3)$. The second panel displays their relative difference. Increasing spin improves this approximation for the parameters shown. It does not remove nutation caused by inappropriate initial conditions.

For a spherical pendulum with nonzero conserved vertical angular momentum $h$, eliminate the azimuthal angular velocity using $h=m{l}^2{\sin }^2\theta \dot{\varphi }$. The polar effective potential is $V_{\mathrm{eff}}=h^2/(2m{l}^2{\sin }^2\theta )-mgl\cos \theta$. A conical pendulum at $0<{\theta }_0<\pi /2$ is a minimum satisfying $h^2=m^{2}g{l}^3{\sin }^4{\theta }_0/\cos {\theta }_0$. Linearization at fixed $h$ yields polar angular frequency ${\omega }_{\theta }^2=V_{\mathrm{eff}}^{\prime \prime }({\theta }_0)/(m{l}^2)={\Omega }^2(1+3{\cos }^2{\theta }_0)$, where ${\Omega }^2=g/(l\cos {\theta }_0)$ is the conical revolution frequency squared. Holding $\dot{\varphi }$ fixed during the perturbation would describe a different, externally driven constraint. The polar perturbation is a nutation; the exactly planar $h=0$ limit requires separate treatment.

Purdue's Bicycle Wheel Gyroscope demonstration uses a bicycle wheel on a long axle, with a ball at one end seated in a freely rotating socket on a heavy stand. Spin the wheel by hand or with the supplied motor, seat the ball, and release the nearly horizontal axle. The readout is the axle sweeping around the vertical while the wheel continues spinning. Compare with the stopped wheel, which falls when one end is released; reverse the spin with the same support and the precession reverses. These are the page's documented apparatus, presentation, and reversal checks. ^[ref]

An alternative in Physics Demonstrations by Julien Clinton Sprott (section 1.3, Materials and Procedure) suspends one axle end by stainless steel wire and a fishing-line swivel. The swivel allows precession without winding up the suspension. Sprott describes the axle bobbing as nutation and precession becoming faster as spin decreases. A suspended support can move, whereas the calculation here assumes a fixed $O$; use the socket version for the closest realization. ^[ref]

A proposed quantitative extension is to record an overhead view of an axle marker against a fixed angle scale, together with a side view for $\theta (t)$. Time successive overhead passages to estimate $\Omega =2\pi /T$; count passages of a rim mark to estimate $|{\omega }_3|=2\pi f$. Record the spin sense too. Measure $M$ and the centre of mass offset $d$ of the whole moving assembly, and determine its axial moment of inertia $I_3$ independently. Wheel, axle, support, release angle, and observation interval should stay the same between spin trials. Measuring both $\Omega$ and $\theta (t)$ distinguishes azimuthal precession from nutation; an overhead movie alone can miss the latter. These camera readouts and parameter checks are proposed here, not reported measurements from the demonstration pages.

For slow precession, compare $\Omega I_3{\omega }_3$ with $Mgd$, and test whether reversing ${\omega }_3$ reverses $\Omega$ with nearly unchanged magnitude. A brief interval minimizes spin loss from friction; repeat at several starting spin rates rather than interpreting an entire slowing run as one steady state. A suspension restoring torque governed by its torsion constant, socket friction, support motion, a misplaced centre of mass, and a push at release can all supply extra torque or change the initial conditions. The swivel reduces the twist that activates the suspension's torsion constant; it does not demonstrate zero torque.

Steady precession keeps $\theta$ constant; nutation changes $\theta$ while the axle may also precess. Sudden release generally excites both. A small initial drop is consistent with conservation of energy and conserved vertical angular momentum: some gravitational potential energy supplies the motion around the vertical, and the axial contribution to vertical angular momentum adjusts. Bearing friction may damp the bobbing, but an ideal undamped heavy symmetric top need not settle into steady precession. Chapter 20: Rotation in space explicitly distinguishes uniform precession from release transients and nutation (section 20-3, paragraphs following equation 20.15).^[ref]