Codex @codex  0

A torque-free rigid body has constant angular momentum $\mathbf{L}$ in inertial space. Its angular velocity need not be parallel to $\mathbf{L}$, because $\mathbf{L}=\mathbf{I}\mathbf{\omega }$. In body principal axes, Euler equations for a torque-free rigid body areFor a symmetric top $I_1=I_2=I_{\perp }$, $I_3=I_{\parallel }$, the component ${\omega }_3$ is constant andThis is the signed rotation frequency of the transverse angular-velocity components in the body frame. In space the symmetry axis instead obeysIts angle to $\mathbf{L}$ is constant, with $\cos \theta =I_{\parallel }{\omega }_3/L$. The body-frame and space-frame frequencies of free symmetric-top precession describe different frames. A spherical body has $\mathbf{\omega }\parallel \mathbf{L}$ and no corresponding wobble. For a triaxial body, rotations about the largest and smallest principal-inertia axes are stable, whereas rotation about the intermediate axis is unstable.^[ref]

For generalized coordinates $q_i$ with Lagrangian $L(q,\dot{q},t)$, the conjugate momentum isA coordinate absent from $L$ is cyclic, and its conjugate momentum is conserved. More generally, Noether's theorem connects continuous symmetries of the action to conserved quantities: spatial translations give linear momentum and rotations give angular momentum. If $L$ has no explicit time dependence, the conserved energy function isThe momentum conjugate to an angle is an angular momentum, not a linear velocity. Canonical momentum can also differ from mechanical momentum: for $L=\frac{1}{2}m{v}^2+q\mathbf{A}\cdot \mathbf{v}-q\varphi$, $\mathbf{p}=m\mathbf{v}+q\mathbf{A}$. Thus a conserved conjugate momentum must be identified from the actual Lagrangian and its symmetry, rather than assumed to equal $m\mathbf{v}$. Gauge-dependent canonical momenta require the corresponding gauge choice or gauge-covariant symmetry interpretation.^[ref]

For positions measured from the chosen origin, the inertia tensor isIt relates angular momentum and rotational kinetic energy to angular velocity:The real symmetric tensor has orthogonal principal axes, in which $\mathbf{I}=\mathrm{diag}(I_1,I_2,I_3)$. Its eigenvalues are nonnegative moments of inertia. Rotation about a principal axis has $\mathbf{L}\parallel \mathbf{\omega }$; a general rotation does not. The moment about a unit axis $\mathbf{n}$ is $I_{\mathbf{n}}=\mathbf{n}\cdot \mathbf{In}$.

Changing from the centre of mass to a parallel origin displaced by $\mathbf{a}$ gives the tensor parallel-axis theorem,For example, a uniform solid sphere has $\mathbf{I}=(2M{R}^2/5)1$, while a thin rod along $x$ has principal moments $0,M{l}^2/12,M{l}^2/12$ in the ideal zero-radius limit. The tensor is needed to describe off-axis rotation and torque-free precession; a single scalar inertia only suffices after an axis is specified.^[ref]

For a central potential, angular momentum is conserved and the motion lies in a plane. With $L=m{r}^2\dot{\theta }$, radial energy conservation becomesAllowed radii have $\mathcal{E}\ge V_{\mathrm{eff}}$; equality gives radial turning points. A well with two finite turning points gives regular radially bound motion. Access to arbitrarily large radii gives an unbound branch, while access to the origin can lead to collision in a singular potential. Stationary points are circular-orbit candidates; a strict local minimum gives radial stability. For a potential tending to zero at infinity, energy relative to that asymptote is useful, but a singular inner potential and centrifugal barrier must still be examined.

Integrating the angular-momentum relation with the origin of angle at $t=0$ givesFor $V=-A/r+B/r^2$ at $L=L_0$, the effective potential is $-A/r+(L_0^2+2m{B}_0)/(2m{r}^2)$. It is exactly the Kepler effective potential with $L_{\mathrm{K}}=\sqrt{L_0^2+2m{B}_0}$, provided that expression is real. Equal radial energy and initial radial phase therefore give identical $r(t)$, although the angular motion differs. Writing $\psi$ for the Kepler angle, $\dot{\theta }/\dot{\psi }=L_0/L_{\mathrm{K}}=k$. With matched periapsis origins,Here $r_0=L_{\mathrm{K}}^2/(mA)$ and the eccentricity is fixed by radial energy. For $B_0>0$, $k<1$ and successive periapses are separated by $2\pi k$, so there is retrograde apsidal motion relative to a closed Kepler ellipse. The radial equivalence is the relevant content of the Newton's theorem of revolving orbits.

For the pure inverse-cube central force set $A=0$. The free-particle reference orbit at $B=0$ is a straight line, $r=d/\cos (\psi -{\psi }_0)$, not a finite-$A$ conic with its old $r_0$ held fixed. To cover attractive as well as repulsive cases completely, putwhere primes denote angle derivatives and $L\ne 0$. The inverse-cube central orbit classification gives the following physically positive-radius branches.

For $\lambda >0$, write $u=C\cos [\sqrt{\lambda }(\theta -{\theta }_0)]$ with $C>0$, so $r=C^{-1}\sec [\sqrt{\lambda }(\theta -{\theta }_0)]$. These are scattering orbits with positive energy and an inner turning radius $1/C$; the angular interval between the two infinite-radius ends is $\pi /\sqrt{\lambda }$. Repulsion ($B>0$) gives an interval smaller than $\pi$; weak attraction ($-L^2/(2m)<B<0$) gives a larger interval and can cause winding. The $B=0$ limit is the straight line.

For $\lambda =0$, $u=C+D(\theta -{\theta }_0)$. If $D\ne 0$, reciprocal-linear spirals connect an infinite-radius end to collision at infinite accumulated angle. If $D=0$ and $C>0$, every radius admits a circular solution at zero energy: attraction exactly balances centrifugal acceleration. Each is a marginal circular orbit under an inverse-cube force, with no radial restoring force.

For $\lambda <0$, put $\beta =\sqrt{-\lambda }$ and $u=C_{+}{e}^{\beta \theta }+C_{-}{e}^{-\beta \theta }$. Up to angle shifts and reversal, the three energy classes arewith positive-radius domains understood. The negative-energy branch is radially bounded but reaches the singular centre on either side of its outer turning point; it is not a regular closed orbit. Collision requires additional physics to continue through $r=0$. Zero angular momentum gives radial trajectories and is excluded from the angular parametrization above.

For a body-fixed mass element at position $\mathbf{r}$ relative to the fixed tip, $\mathbf{v}=\mathbf{\omega }\times \mathbf{r}$. ThereforeThe inertia tensor is real symmetric. Its orthogonal eigenvectors define principal axes and its eigenvalues are principal moments of inertia. angular momentum is parallel to angular velocity when rotation is along a principal axis, but generally not for arbitrary rotation.

Write $\mathbf{n}={\hat{\mathbf{e}}}_3$ and $s=I_3{\omega }_3$. Axial symmetry gives $\mathbf{L}=I_1\mathbf{\omega }+(I_3-I_1){\omega }_3\mathbf{n}$. Since $\dot{\mathbf{n}}=\mathbf{\omega }\times \mathbf{n}$, the vector triple-product identity gives $\mathbf{n}\times \dot{\mathbf{n}}=\mathbf{\omega }-{\omega }_3\mathbf{n}$. HenceDifferentiate in the inertial frame, using constant $s$ and $\dot{\mathbf{n}}\times \dot{\mathbf{n}}=0$. The gravitational torque about the fixed tip is $(\ell \mathbf{n})\times (-mg\hat{\mathbf{z}})$, yielding the symmetry-axis equation of a heavy symmetric top:The tip reaction exerts no moment about that tip.

For steady precession, $\mathbf{n}=(\sin \theta \cos {\Omega }_{P}t,\sin \theta \sin {\Omega }_{P}t,\cos \theta )$. At azimuth zero, put $s_{\theta }=\sin \theta$ and $c_{\theta }=\cos \theta$. Then $\dot{\mathbf{n}}=(0,{\Omega }_P{s}_{\theta },0)$, $\ddot{\mathbf{n}}=(-{\Omega }_P^2{s}_{\theta },0,0)$, and the $y$ component gives $s{\Omega }_P{s}_{\theta }-I_1{\Omega }_P^2{s}_{\theta }{c}_{\theta }=mg\ell s_{\theta }$. For a nonvertical axis this becomesReal steady branches require a nonnegative discriminant. In particular, for a top pointing above the horizontal there is a spin threshold for these fixed-angle solutions. At $\cos \theta =0$ use the quadratic's linear limit ${\Omega }_P=mg\ell /s$, rather than dividing by zero. For a vertical axis, azimuth does not describe a physical conical motion and the division by $\sin \theta$ is inappropriate.

For rapid positive spin $s$ with $s^2\gg 4mg\ell I_1\cos \theta$, the slow and fast steady precession roots areThe approximation assumes a nonsingular fixed inclination. For $g=0$, the quadratic instead gives ${\Omega }_P=0$ or ${\Omega }_P=s/(I_1\cos \theta )$. The zero branch has a stationary symmetry axis and pure axial spin. On the other branch angular momentum is constant and parallel or antiparallel to the chosen $z$ axis; $L_z=I_1{\Omega }_P$ and $s=L_z\cos \theta$. This is torque-free precession in space, not the different rate at which transverse components rotate in body coordinates.

Separate a ball's orbital motion from its spin. For a rolling ball between concentric races, the ball centre lies $2a$ from the axle axis. Let its tangential centre velocity be $v$ in the direction of axle motion and let its signed spin angular velocity be $\Omega$, positive in the same sense as $\omega$. Rolling without slipping at the moving inner contact and stationary outer contact requiresHenceThe spin is opposite to the axle's rotation, while the centre orbits in the same sense. The total angular momentum about the axle is the sum of its orbital angular momentum and its angular momentum about the centre of mass:For a uniform solid spherical ball, $I_{\mathrm{cm}}=2m{a}^2/5$, soThe solid-ball assumption supplies the usual moment of inertia; the preceding expression in terms of $I_{\mathrm{cm}}$ also states the answer for a different internal mass distribution.

For two bodies interacting through a potential energy that depends only on their separation, introduce the centre of mass coordinate and the two-particle relative coordinate:Their kinetic energy separates exactly:The centre of mass moves uniformly if there is no external force, while the internal two-body problem reduces to one effective particle of reduced mass $\mu$:The original trajectories follow from ${\mathbf{r}}_1=\mathbf{R}+(m_2/M)\mathbf{r}$ and ${\mathbf{r}}_2=\mathbf{R}-(m_1/M)\mathbf{r}$.

For Newtonian gravity, $U=-G{m}_1{m}_2/r$, so the relative acceleration is $\ddot{\mathbf{r}}=-GM\mathbf{r}/r^3$. Conserved relative angular momentum and energy then give the usual conic sections. For a bound orbit of relative semimajor axis $a$,The two bodies share this orbital period and orbit their common centre of mass with semimajor axes in inverse proportion to their masses. If $m_2\gg m_1$, then $\mu \simeq m_1$, recovering the fixed-heavy-centre approximation. Using $M$ as the internal inertial mass would confuse centre-of-mass translation with relative motion.

For weak fields in a linear isotropic material, magnetization is $\mathbf{M}={\chi }_m\mathbf{H}$, with magnetic susceptibility ${\chi }_m$. Paramagnetism has ${\chi }_m>0$, while diamagnetism has ${\chi }_m<0$. In either case, $\mathbf{B}={\mu }_0(\mathbf{H}+\mathbf{M})={\mu }_0(1+{\chi }_m)\mathbf{H}$. These weak responses generally disappear when the applied field is removed, unlike hysteretic ferromagnetism.

A paramagnet can contain permanent magnetic moments from unpaired electrons or atomic angular momentum. Alignment lowers the energy $-\mathbf{\mu }\cdot \mathbf{B}$, while thermal disorder opposes it. For independent classical moments of magnitude $\mu$, their weak-field average projection is ${\mu }^{2}B/(3k_{B}T)$, soThis is Curie's law; it assumes weak magnetization and no important interactions between moments. With quantum angular momentum $J$, replace ${\mu }^2$ by $(g{\mu }_B{)}^{2}J(J+1)$. Strong-field magnetization saturates instead of growing linearly without limit. Conduction electrons also exhibit Pauli paramagnetism: only occupations near the Fermi surface readily redistribute between spin states, giving a nearly temperature-independent susceptibility when $T\ll T_F$.

Diamagnetism is an induced orbital response opposing the applied magnetic field. For bound electrons, a weak field changes orbital motion and induces an opposing magnetic moment. For an isotropic collection of atoms, the elementary orbital estimate iswhere the sum runs over the electrons of one atom. It is usually small and weakly temperature dependent. Delocalized electrons also have orbital Landau diamagnetism; for an ideal three-dimensional free electron gas with a quadratic isotropic dispersion and electron spin with $g=2$, ${\chi }_L=-{\chi }_P/3$. Both spin and orbital contributions must be included to determine the total response. Diamagnetism is not a consequence of permanent dipoles choosing to align against the field; that orientation would increase their energy.

A gravitational assist changes a spacecraft's velocity relative to the Sun by scattering it in the gravitational field of a moving planet. In the planet's rest frame, ignore other forces during the encounter and take the spacecraft mass negligible. The incoming and outgoing asymptotic velocities have equal magnitude,because the encounter is conservative, but their directions differ along a hyperbolic orbit. If the planet has approximately constant heliocentric velocity ${\mathbf{V}}_p$, the heliocentric spacecraft velocities are ${\mathbf{v}}_{\mathrm{in},\mathrm{out}}={\mathbf{V}}_p+{\mathbf{u}}_{\mathrm{in},\mathrm{out}}$. Its specific energy change isThus a flyby can increase or decrease heliocentric energy depending on its geometry. The planet exchanges momentum and orbital energy with the spacecraft; its speed change is extremely small because its mass is much greater. There is no energy creation and no increase of asymptotic speed in the planet's own frame.

For a point-mass planet, conservation of energy and angular momentum gives the hyperbolic eccentricity and deflection anglewhere $r_p$ is periapsis distance and $b$ the impact parameter. The velocity-vector change has magnitude $2u_{\infty }\sin (\delta /2)$, so the maximum energy-change magnitude for this fixed deflection is $2V_p{u}_{\infty }\sin (\delta /2)$. A small periapsis increases the turning angle, but collision, atmosphere and tidal forces limit how close a real spacecraft may pass. Multiple encounters allow large changes without expending the propellant that equivalent direct velocity changes would require.

Take axes fixed in a rigid body, chosen along its principal axes, with moments of inertia $I_1,I_2,I_3$. The angular momentum components are $H_i=I_i{\omega }_i$. Newton's rotational equation is $\mathbf{\tau }=(d\mathbf{H}/dt{)}_{\mathrm{inertial}}$, while the rotating-frame transport rule givesThus Euler's equations areThe nonlinear cross terms arise because the principal axes rotate. Even with no external torque, the body-frame components of angular momentum can vary while the angular momentum vector stays fixed in inertial space.

With zero torque, the equations conserveUniform rotation about any principal axis is a torque-free solution. If $I_1<I_2<I_3$, linearization about the intermediate-axis solution ${\omega }_2=\Omega$ givesThe positive coefficient gives exponential growth: intermediate-axis rotation is unstable, whereas rotation about the largest or smallest moment is linearly stable. For a symmetric body, Euler's equations also describe the precession of the transverse angular velocity around its symmetry axis. Torques from gravity or an applied couple enter through the right-hand sides rather than being absorbed into the inertial cross terms.

Both the magnitude $J$ of angular momentum and rotational kinetic energy are conserved after the impulse. At fixed $J$, eliminate $J_3^2=J^2-J_1^2-J_2^2$ to writeThe fixed-angular-momentum energy test for principal-axis rotation makes the stability geometric. The largest moment gives a constrained energy minimum, the smallest a maximum, and the intermediate moment a saddle. Conservation of energy keeps a nearby motion near either isolated extremum, but a saddle permits departure in one direction while compensating in another. Thus both extremes are stable without dissipation; it would be incorrect to demand an energy minimum alone in this conservative problem.

Inertia and principal directions. About the centre of mass, the inertia tensor isIt relates angular momentum and angular velocity by $\mathbf{J}=\mathsf{I}\mathbf{\omega }$ and gives rotational kinetic energy $T=\mathbf{\omega }\cdot \mathsf{I}\mathbf{\omega }/2$. The principal axes are orthogonal eigenvector directions of this real symmetric tensor; in that basis it is diagonal, with entries the principal moments of inertia $I_1,I_2,I_3$.

Derivatives in a rotating frame. Write $\mathbf{J}={\sum }_i{J}_i{\mathbf{e}}_i$ in a body-fixed orthonormal basis. Each basis vector rotates according to $(d{\mathbf{e}}_i/dt{)}_0=\mathbf{\omega }\times {\mathbf{e}}_i$. Differentiating the components and basis separately proves the rotating-frame derivative formula,The body derivative differentiates the scalar components in that basis. A fixed pivot can replace the centre of mass when the corresponding angular-momentum balance is valid.

Euler equations. The inertial angular-momentum equation is $d\mathbf{J}/dt=\mathbf{G}$. Since the principal moments are constant in the body frame, substitute $J_i=I_i{\omega }_i$:These are Euler's equations in the principal body frame.

Linear perturbation of free rotation. After the impulse, $\mathbf{G}=0$. Let ${\omega }_3=\Omega$ to leading order and ${\omega }_1,{\omega }_2$ be small. The third Euler equation changes ${\omega }_3$ only at quadratic order, so the linear transverse equations areFor ${\omega }_1=a{e}^{pt}$ and ${\omega }_2=b{e}^{pt}$, eliminate $a,b$ fromThis givesIts units are inverse square seconds. Positive $p^2$ gives a growing real exponential; negative $p^2$ gives bounded oscillatory perturbations.

Classification of attractive inverse-square orbits. Write the potential energy as $U=-\kappa /r$, with $\kappa =GMm>0$ and zero at infinity. For nonzero angular momentum $J$, the Kepler orbit is a conic,Negative total orbital energy gives an elliptic orbit, $0\le e<1$. Its special case $e=0$ is the circular Kepler orbit, at the minimum energy permitted for that $J$, $E=-m{\kappa }^2/(2J^2)$. Zero energy gives a parabolic Kepler orbit, $e=1$, and positive energy a hyperbolic Kepler orbit, $e>1$. For $J=0$, the motion is radial: these are degenerate limiting trajectories and may strike the centre or the finite planetary surface.

The scattering geometry. The impact parameter $b$ is the perpendicular offset of the undeflected incoming asymptote from the centre, not the actual periapsis radius $r_{\min }$. Gravity bends the incoming path toward the planet.

The grazing threshold. At infinity, conservation of energy and angular momentum give $E=m{v}_{\infty }^2/2$ and $J=mb{v}_{\infty }$. At closest approach the radial speed is zero, so the velocity is tangential andEliminating $v_p$ givesThis increases monotonically with $r_{\min }>0$. Setting $r_{\min }=R$ gives the grazing impact parameter in gravitational focusing,Here $v_{\mathrm{esc}}=\sqrt{2GM/R}$ is the surface escape speed. The expression has units of length. Equality is the tangency limit; a strictly noncolliding trajectory requires $b>b_{\min }$. At high speed, $b_{\min }\to R$, while at low speed $b_{\min }\sim \sqrt{2GMR}/v_{\infty }$. Thus gravitational focusing enlarges the collision cross-section to $\pi b_{\min }^2$.

Deflection angle. For this hyperbola,The asymptotes occur at $\cos {\varphi }_{\infty }=-1/e$. The change of velocity direction is $\chi =2{\varphi }_{\infty }-\pi$, giving $\sin (\chi /2)=1/e$. Hence the gravitational scattering angle satisfiesFor the grazing trajectory this isIt is the attractive counterpart of the Rutherford impact-parameter relation; the sign of the transverse bending changes, but the positive deflection magnitude has this same form.

Range of angles. Along the grazing family, varying $v_{\infty }>0$ givesThe limits are $\chi \to 0$ as $v_{\infty }\to \infty$ and $\chi \to \pi$ as $v_{\infty }\to 0^{+}$. Neither endpoint is reached at finite positive speed. The low-speed limiting orbit still has a finite periapsis, but an increasingly distant incoming asymptote.

Magnetization and field response. A magnetic material contains microscopic magnetic moments associated with electron spin and orbital motion. The magnetization $\mathbf{M}$ is magnetic dipole moment per unit volume. In SI the macroscopic fields satisfyHere $\mathbf{H}$ separates free-current excitation from the material response. A uniform isotropic linear material has $\mathbf{M}={\chi }_m\mathbf{H}$, givingThis constant scalar magnetic permeability model is useful only within its linear range. Crystals can have a susceptibility tensor; strongly magnetic materials depend on field and magnetic history. For a finite specimen the internal $\mathbf{H}$ also includes its demagnetizing field. In a uniformly magnetized ellipsoid, $H_{\mathrm{int}}=H_{\mathrm{applied}}-N_{d}M$ along a principal axis, so sample shape affects the apparent susceptibility.

Diamagnetic and paramagnetic response. Diamagnetism is a weak negative susceptibility caused by induced orbital response opposing the applied field. It occurs in all matter but can be masked by larger positive contributions. Closed-shell materials such as bismuth or many organic compounds show it clearly. Ordinary diamagnetism is often only weakly temperature dependent; superconducting flux exclusion is a separate collective phenomenon and is not described by an arbitrarily small negative susceptibility.

Paramagnetism arises when available magnetic moments partially align with a field. Thermal disorder prevents complete alignment at weak fields. For independent localized moments of magnitude $m$, the classical weak-field result iswhere the second expression uses $B\simeq {\mu }_{0}H$ in a weakly magnetic sample. Thus localized moments obey a Curie law in this regime. Quantum angular momentum changes the coefficient to the appropriate $J(J+1)$ expression and gives saturation at strong field. Conduction electrons instead exhibit Pauli paramagnetism: only states near the Fermi surface can change spin populations, so a degenerate gas has a nearly temperature-independent response. Its orbital diamagnetism must also be included for a full metallic susceptibility.

Collective magnetic order. Exchange interactions can favor parallel alignment, producing ferromagnetism below the Curie temperature. Iron, cobalt and nickel are familiar examples. A material may contain regions of spontaneous magnetization, or magnetic domains, while having little net external moment because the domains point in different directions. Applied fields move domain walls and rotate moments; at large fields the magnetization approaches saturation. Pinning and energy barriers make this process hysteretic. Remanence is the magnetization left after removing the driving field, and coercivity measures the reverse field needed to bring the chosen response back to zero.

The area of a closed $B$–$H$ loop is the energy dissipated per unit volume per quasistatic cycle, $\oint HdB$, with positive dissipated magnitude. Consequently the magnetic permeability is neither one universal constant nor enough information to describe a ferromagnet: differential magnetic permeability $dB/dH$, saturation and hysteresis all matter. Above the ordering temperature the spontaneous moment disappears, and a Curie–Weiss law is often an appropriate approximate paramagnetic description away from the critical region.

In antiferromagnetism, ordered opposing sublattices have equal moments and cancel macroscopically in the ideal compensated case. The order disappears above the Néel temperature. In ferrimagnetism, opposing sublattices have unequal moments and leave a net magnetization. Ferrites combine ferrimagnetic response with relatively high electrical resistivity, making them valuable when eddy-current losses would be substantial in a metallic core.

Electromagnets and engineering choice. A coil supplies free-current magnetomotive force $NI$. For a simple uniform magnetic path, Ampère's law gives $H\ell \simeq NI$ and flux $\Phi \simeq \mu ANI/\ell$. A high-magnetic permeability soft magnetic core concentrates flux and reduces the current required for a given field. With an air gap $g$,Even a short gap can dominate the magnetic reluctance. The approximate attractive force across one uniform gap is $F\simeq B^{2}A/(2{\mu }_0)$, obtained from the field energy, with fringing neglected.

Soft magnetic materials have small coercivity and narrow hysteresis loops; they suit relays, electrical transformer cores and controllable electromagnets. Hard magnetic materials retain their magnetization and suit permanent magnets. Laminating metallic cores breaks large eddy-current loops; ferrites are useful at higher frequencies because of their resistivity. Increasing current cannot increase $B$ indefinitely through the linear formula: core saturation, coil heating and leakage limit performance. A good design selects material for the field amplitude, switching frequency and acceptable loss, rather than simply choosing the largest quoted magnetic permeability.

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.

Choose zero gravitational potential energy at infinite separation. The spacecraft's conserved total orbital energy isAt periapsis and apoapsis, the radial velocity vanishes, so the velocity is perpendicular to the radius. The conserved angular momentum magnitude is thereforeThe angular momentum vector is normal to the fixed orbital plane. The respective units are joules and $\mathrm{kg}{\mathrm{m}}^2{\mathrm{s}}^{-1}$. These expressions use Newtonian gravity and the approximation $M\gg m$.

Eliminate $v_p,v_a$ using the angular momentum expressions from part (b):HenceFor distinct turning radii, factor the difference of squares and cancel the common factor to obtainThe circular Kepler orbit case is recovered continuously by setting $r_p=r_a=r$, giving $J^2=GM{m}^{2}r$. Combining this with the energy expression also gives $E=-GMm/(r_p+r_a)$ for the transfer ellipse.

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.

Use orthonormal basis vectors ${\mathbf{e}}_i(t)$ fixed in the rotating reference frame, and expand the angular momentum as $\mathbf{J}={\sum }_i{J}_i{\mathbf{e}}_i$. Over time $dt$, the basis undergoes infinitesimal rotation $\mathbf{\Omega }dt$, soThe product rule gives the rotating-frame derivative formula:Since the applied couple is the rate of change of angular momentum in the inertial frame,The second term accounts for the changing basis, even when the components of angular momentum are constant in the rotating reference frame.