Advanced University Physics

A Calculus-Based Exploration of Gravitation

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$.

$$ \vec{F}_{12} = - G \frac{m_1 m_2}{r^2} \hat{r}_{12} $$

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.

$$ \vec{g}(\vec{r}) = -G \frac{M}{r^2} \hat{r} $$

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.

$$ \vec{g}_{total} = \sum_i \vec{g}_i = -G \sum_i \frac{M_i}{r_i^2} \hat{r}_i $$

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:

$$ W(h) = F_g = G \frac{M_E m}{(R_E + h)^2} = mg(h) $$

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$:

$$ g(h) \approx g_0 \left(1 - 2\frac{h}{R_E}\right) $$

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:

$$ g(r) = \left(\frac{GM_E}{R_E^3}\right) r = g_0 \frac{r}{R_E} $$

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$.

$$ \Delta U = U_f - U_i = - \int_{r_i}^{r_f} \vec{F}_g \cdot d\vec{r} $$

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:

$$ U(r) = - G \frac{m_1 m_2}{r} $$

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:

$$ E = K + U = \frac{GMm}{2r} - \frac{GMm}{r} = - \frac{GMm}{2r} $$

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 $$
$$ v_{esc} = \sqrt{\frac{2GM}{R}} $$

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:

$$ \frac{dA}{dt} = \frac{1}{2} |\vec{r} \times \vec{v}| = \frac{|\vec{L}|}{2m} = \text{constant} $$

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:

$$ T^2 = \left(\frac{4\pi^2}{GM}\right) a^3 \implies T^2 \propto a^3 $$

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.

$$ \vec{F}_{net} = \vec{F}_g + \vec{F}_N = m\vec{a}_c $$

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.

Problem Solutions

Intermediate Problem Solutions

1. Hohmann Transfer Orbit Analysis

The total change in velocity is the sum of two burns. The first burn moves the satellite from the initial circular orbit to the elliptical transfer orbit. The second circularizes the orbit at the new radius.
$$ \Delta v_{total} = \sqrt{\frac{GM}{r_1}} \left(\sqrt{\frac{2r_2}{r_1+r_2}} - 1\right) + \sqrt{\frac{GM}{r_2}} \left(1 - \sqrt{\frac{2r_1}{r_1+r_2}}\right) $$

2. Oscillatory Motion in a Gravity Tunnel

The component of gravitational force along the tunnel is $F_x = -(mg_0/R_E)x$, which is a linear restoring force, characteristic of SHM. The period is:
$$ T = 2\pi \sqrt{\frac{R_E}{g_0}} $$

3. The Sun-Earth L1 Lagrange Point

In the inertial frame, the net force on the third body provides the centripetal force for its orbit: $-\frac{GM_S m}{(R-r)^2} + \frac{GM_E m}{r^2} = -m\omega^2(R-r)$. We use Earth's orbital angular velocity, $\omega^2 \approx GM_S/R^3$. Using the binomial approximation $(1-x)^{-2} \approx 1+2x$ for $x=r/R$, the force equation simplifies to:
$$ r \approx R \left( \frac{M_E}{3M_S} \right)^{1/3} $$

4. Energy Requirement for Orbital Insertion

The total energy in a circular orbit is $E = -GMm/(2r)$. The work required is the difference in energy between the two orbits. Let $r_1=R_E+h_1$ and $r_2=R_E+h_2$.
$$ W = \Delta E = E_2 - E_1 = \left(-\frac{GM_E m}{2r_2}\right) - \left(-\frac{GM_E m}{2r_1}\right) = \frac{GM_E m}{2} \left(\frac{1}{r_1} - \frac{1}{r_2}\right) $$

5. Dynamics of a Binary Star System

Let star A orbit at radius $r_A$ and star B at $r_B$. The gravitational force provides the centripetal force for each: $G M_A M_B / d^2 = M_A v_A^2/r_A$. Using $v_A = 2\pi r_A/T$ and the center of mass relation $M_A r_A = M_B r_B$ with $d=r_A+r_B$, we can solve for T.
$$ T = 2\pi \sqrt{\frac{d^3}{G(M_A+M_B)}} $$

Hard Problem Solutions

6. Gravitational Self-Energy of a Sphere

By integrating the work done to add successive shells of mass, $dW = V(r)dm$, from $r=0$ to $r=R$, we find the total self-energy.
$$ U_g = -\frac{3}{5}\frac{GM^2}{R} $$

7. The Gravitational Slingshot

Transform to planet's frame: $\vec{v}'_i = \vec{v}_i - \vec{V}$. After elastic interaction: $|\vec{v}'_f|=|\vec{v}'_i|$. For head-on, $\vec{v}'_f = -\vec{v}'_i$. Transform back to Sun's frame: $\vec{v}_f = \vec{v}'_f + \vec{V} = -(\vec{v}_i - \vec{V}) + \vec{V} = -\vec{v}_i + 2\vec{V}$.
$$ \Delta K = \frac{1}{2}m(v_f^2 - v_i^2) = \frac{1}{2}m((|-v_i+2V|)^2 - v_i^2) = 2m V (v_i + V) $$

8. Gravitational Field of a Finite Rod

Let $\lambda = M/L$. The field from an element $dx$ at position $x$ is $d\vec{g}$. By symmetry, x-components cancel. The y-component is $dg_y = -G(\lambda dx) \cos\theta / r^2$, where $r=\sqrt{x^2+y^2}$ and $\cos\theta=y/r$. Integrating from $x=-L/2$ to $x=L/2$:
$$ \vec{g}(0,y) = -\frac{2GM}{y\sqrt{L^2+4y^2}} \hat{j} $$

9. Orbital Stability with a Perturbation Potential

The effective potential is $U_{eff}(r) = U(r) + L^2/(2mr^2)$. For a stable orbit, we need $d^2U_{eff}/dr^2 > 0$ at the radius of the circular orbit (where $dU_{eff}/dr=0$). This analysis yields the condition:
$$ L^2 > 6m|\beta| \implies \text{The orbit must have sufficient angular momentum.} $$

10. Tidal Stress on a Radially Infalling Rod

The tensile force at the midpoint is the integral of the differential gravitational force on the outer half of the rod. $F_{tidal} = \int_0^{L/2} (F_g(r+x) - F_g(r)) dm$. Approximating the force difference as $dF/dr \cdot x$ and integrating yields the tension $T$. Stress $\sigma = T/A$.
$$ \sigma = \frac{T}{A} = \frac{G M_P m L}{4 A r^3} $$

Advanced (Irodov-like) Problem Solutions

11. Escape Velocity from a Rotating Body

In the inertial frame, the object's initial velocity is $\vec{v}_{launch} = \vec{v}_{rel} + \vec{v}_{rot}$, where $v_{rot} = \omega R$. To escape, total energy must be zero: $\frac{1}{2}m v_{launch}^2 - GMm/R = 0$. This requires $v_{launch} = \sqrt{2GM/R} = v_{esc,static}$. To minimize $v_{rel}$, we must maximize the contribution from rotation by launching prograde (in the direction of rotation).
$$ v_{rel,min} = v_{esc,static} - v_{rot} = \sqrt{\frac{2GM}{R}} - \omega R $$

12. Bertrand's Theorem and Orbital Stability

The condition for closed orbits is that the angular frequency of radial perturbations, $\omega_r$, is a rational multiple of the orbital angular frequency, $\omega_\phi$. For stability for all initial conditions, we require $\omega_r = \omega_\phi$. Analyzing the effective potential for $F(r)=-k/r^n$ shows this condition is only met for $n=2$ (gravity) and $n=-1$ (Hooke's Law, harmonic oscillator).
$$ n = -1 \quad (F(r) = -kr, \text{ a linear restoring force}) $$

13. Field Inside a Spherical Cavity

Using superposition, the field inside is the field of the full sphere ($\vec{g}_{full}$) minus the field of the mass that would have filled the cavity ($\vec{g}_{cavity\_mass}$). Let $\vec{r}$ be a point from the large sphere's center and $\vec{b}$ be the cavity's center. The field is uniform.
$$ \vec{g}_{cavity} = \frac{4}{3}\pi G \rho \vec{b} $$

14. Derivation of the Roche Limit

Consider a point mass on the satellite's surface closest to the planet. The planet's gravity pulls it, while the satellite's gravity and the centripetal force (in the satellite's frame) push it away. The tidal force is the difference in the planet's gravity across the satellite. The limit is where the tidal force exceeds the satellite's self-gravity. $F_{tidal} \approx \frac{2GM_p r_m}{d^3} > \frac{G m_m}{r_m^2}$. Substituting densities yields the result.
$$ d \approx R_p \left(2 \frac{\rho_p}{\rho_m}\right)^{1/3} $$

15. Apsidal Precession from a Perturbation

The angular frequencies of radial ($\omega_r$) and orbital ($\omega_\phi$) motion are found from the effective potential $U_{eff}$. $\omega_r^2 = \frac{1}{m}\frac{d^2U_{eff}}{dr^2}|_{r_0}$ and $\omega_\phi = L/(mr_0^2)$. The precession angle per orbit is $\Delta\phi = 2\pi(\omega_\phi/\omega_r - 1)$. For the given potential, this evaluates to:
$$ \Delta\phi \approx \frac{6\pi h m}{L^2} $$