1. Newton’s Law of Gravitation
Newton's Law of Universal Gravitation provides the foundational framework for understanding celestial mechanics. It posits that the gravitational force between two point masses, $m_1$ and $m_2$, is attractive, directly proportional to the product of their masses, and inversely proportional to the square of their separation distance, $r$.
Here, $G$ is the universal gravitational constant, and $\hat{r}_{12}$ is the unit vector pointing from $m_1$ to $m_2$. The negative sign signifies the attractive nature of the force. The associated gravitational field $\vec{g}$ created by a mass $M$ is a vector field representing the force per unit mass at any point in space.
For a collection of masses, the total gravitational field at a point is the vector sum of the fields from each individual mass, an application of the Principle of Superposition.
2. Weight
Weight is the gravitational force exerted on an object by a large astronomical body. For an object of mass $m$ at a height $h$ above the Earth's surface (mass $M_E$, radius $R_E$), its weight is given by:
Proof: Variation of $g$ with Altitude
For altitudes $h \ll R_E$, a first-order approximation can be found using the binomial expansion $(1+x)^n \approx 1+nx$ for small $x$. We can write $g(h)$ as:
$$ g(h) = G \frac{M_E}{R_E^2(1 + h/R_E)^2} = g_0 (1 + h/R_E)^{-2} $$Where $g_0 = GM_E/R_E^2$ is the acceleration at the surface. Using the binomial approximation with $x=h/R_E$ and $n=-2$:
Proof: Variation of $g$ with Depth
Assuming the Earth has a uniform density $\rho$, the gravitational force on a mass $m$ at a distance $r$ from the center ($r < R_E$) is due only to the mass enclosed within the sphere of radius $r$, according to the Shell Theorem. The mass enclosed is $M_{enc} = M_E (r^3/R_E^3)$. The gravitational force is $F_g(r) = G M_{enc}m/r^2$.
Since $F_g = mg(r)$, the gravitational acceleration inside the Earth is:
This shows that $g$ decreases linearly as one approaches the center of the Earth, becoming zero at the center.
3. Gravitational Potential Energy
Since the gravitational force is conservative, we can define a scalar potential energy function $U$ such that the work done by gravity is $W_g = -\Delta U$. The change in potential energy is the negative of the work done by the gravitational force in moving a mass $m$ from an initial position $r_i$ to a final position $r_f$.
Substituting $\vec{F}_g = -G M m/r^2 \hat{r}$ and $d\vec{r} = dr \hat{r}$, the integral becomes:
$$ \Delta U = - \int_{r_i}^{r_f} \left(- G \frac{M m}{r^2}\right) dr = GMm \int_{r_i}^{r_f} \frac{1}{r^2} dr = GMm \left[-\frac{1}{r}\right]_{r_i}^{r_f} = -G M m \left(\frac{1}{r_f} - \frac{1}{r_i}\right) $$By defining the zero point of potential energy at infinite separation ($U(\infty) = 0$), the potential energy of a two-mass system is always negative:
4. The Motion of Satellites
Proof: Total Energy of a Circular Orbit
For a satellite in a circular orbit, the gravitational force provides the centripetal force: $G Mm/r^2 = mv^2/r$. The kinetic energy is $K = \frac{1}{2} mv^2 = GMm/(2r)$. The potential energy is $U = - GMm/r$. The total energy is the sum, a result known as the Virial Theorem for this system:
The negative energy signifies a bound system. An object with $E \ge 0$ is unbound and will follow a parabolic or hyperbolic trajectory.
Proof: Escape Velocity
Escape velocity is the minimum speed an object must have to escape a gravitational field, i.e., to reach $r=\infty$ with zero final speed. We use conservation of energy. The minimum energy for escape is $E_{total}=0$. Initial energy at radius $R$ must equal the final energy at infinity:
$$ E_i = E_f \implies \frac{1}{2}mv_{esc}^2 - \frac{GMm}{R} = 0 + 0 $$5. Kepler’s Laws and the Motion of Planets
Proof of Kepler's Second Law (Law of Areas)
This law is a direct consequence of the conservation of angular momentum. Since gravity is a central force, the torque $\vec{\tau} = \vec{r} \times \vec{F}_g$ is zero, which implies the angular momentum $\vec{L} = \vec{r} \times \vec{p}$ is constant. The infinitesimal area $dA$ swept by the position vector $\vec{r}$ in time $dt$ is $dA = \frac{1}{2} |\vec{r} \times d\vec{r}| = \frac{1}{2} |\vec{r} \times \vec{v}dt|$. The rate of change of area is:
Proof of Kepler's Third Law (Law of Periods)
The total area of an ellipse is $A = \pi ab$, where $a$ is the semi-major axis. The orbital period is $T = A / (dA/dt)$. Using $dA/dt = L/2m$ and expressions for an ellipse relating $L$ and $a$, we find:
$$ T = \frac{2m(\pi a^2 \sqrt{1-\epsilon^2})}{\sqrt{GMm^2 a(1-\epsilon^2)}} = \frac{2\pi a^{3/2}}{\sqrt{GM}} $$Squaring both sides gives the celebrated result:
6. Spherical Mass Distributions
Proof of Shell Theorem via Gauss's Law for Gravity
Gauss's Law for Gravity states that the flux of the gravitational field through any closed surface is proportional to the enclosed mass: $\oint \vec{g} \cdot d\vec{A} = -4\pi G M_{enc}$.
Outside the shell ($r > R$): We choose a spherical Gaussian surface of radius $r$. By symmetry, $\vec{g}$ is radial. So, $g(4\pi r^2) = -4\pi G M_{enc}$. The enclosed mass is the entire shell, $M_{enc}=M$.
$$ g(4\pi r^2) = -4\pi GM \implies g = -G\frac{M}{r^2} $$Inside the shell ($r < R$): The enclosed mass is zero, $M_{enc}=0$, so $g=0$. This elegantly proves both parts of the Shell Theorem.
7. Apparent Weight and the Earth’s Rotation
An object on the surface of the rotating Earth is in a non-inertial reference frame. Its apparent weight, $\vec{W}_{app}$, (the normal force $\vec{F}_N$) differs from its true gravitational weight, $\vec{F}_g$. The net force provides the centripetal force required for its circular path of radius $r_c = R_E \cos\lambda$ around the Earth's axis, where $\lambda$ is the latitude.
The apparent weight is $\vec{W}_{app} = \vec{F}_N = m\vec{g}_{true} - m\vec{a}_c$. The centripetal acceleration, $a_c = \omega^2 R_E \cos\lambda$, is largest at the equator ($\lambda=0$) and zero at the poles ($\lambda=90^\circ$).
Advanced Problems in Gravitation
Intermediate Problems
1. Hohmann Transfer Orbit Analysis
An artificial satellite is in a circular orbit of radius $r_1$ around the Earth. To move it to a coplanar circular orbit of radius $r_2 > r_1$, a Hohmann transfer orbit is employed. Determine the total change in velocity required for the maneuver.
2. Oscillatory Motion in a Gravity Tunnel
A straight, frictionless tunnel is bored through a non-rotating Earth of uniform density. Demonstrate that an object dropped into the tunnel undergoes simple harmonic motion and derive an expression for its period.
3. The Sun-Earth L1 Lagrange Point
Determine the distance $r$ of the L1 Lagrange point from the Earth, assuming the Earth (mass $M_E$) and Sun (mass $M_S$) are in circular orbits at a separation $R$. Provide a linearized solution for $r \ll R$.
4. Energy Requirement for Orbital Insertion
A probe of mass $m$ is in a circular LEO at altitude $h_1$. Calculate the minimum energy required to place it into a circular GEO at altitude $h_2$.
5. Dynamics of a Binary Star System
Two stars of mass $M_A$ and $M_B$ are separated by a constant distance $d$. Derive an expression for the orbital period $T$ of the system.
Hard Problems
6. Gravitational Self-Energy of a Sphere
Calculate the gravitational self-energy $U_g$ of a sphere of uniform mass density, with total mass $M$ and radius $R$.
7. The Gravitational Slingshot
A spacecraft (initial velocity $\vec{v}_i$) performs a hyperbolic fly-by of a planet (velocity $\vec{V}$). For a head-on encounter, determine the spacecraft's final velocity $\vec{v}_f$ and the change in its kinetic energy.
8. Gravitational Field of a Finite Rod
A uniform rod of mass $M$ and length $L$ lies on the x-axis. Derive the gravitational field $\vec{g}$ at a point $(0, y)$.
9. Orbital Stability with a Perturbation Potential
A test mass moves in a potential $U(r) = -k/r - \beta/r^3$. Determine the condition on $\beta$ that permits stable circular motion.
10. Tidal Stress on a Radially Infalling Rod
A uniform rod of length $L$ falls radially towards a planet. Derive the tensile stress at its midpoint.
Advanced (Irodov-like) Problems
11. Escape Velocity from a Rotating Body
An object is at the equator of a rotating asteroid. Determine the minimum launch speed $v_{rel}$ relative to the surface required for escape, and find the optimal launch direction.
12. Bertrand's Theorem and Orbital Stability
For a central force $F(r)=-k/r^n$, only $n=2$ and one other integer value of $n$ result in stable, closed orbits. Find this other value.
13. Field Inside a Spherical Cavity
A uniform sphere has an off-centre spherical cavity. Prove that the gravitational field inside the cavity is uniform and derive its value.
14. Derivation of the Roche Limit
A fluid satellite orbits a planet. Derive the orbital distance (Roche limit) at which tidal forces from the planet overcome the satellite's self-gravity, disintegrating it.
15. Apsidal Precession from a Perturbation
For a potential $U(r) = -k/r + h/r^2$, calculate the angle of perihelion advance per revolution for a nearly circular orbit.