Comprehensive University Study Notes on Gravitation and Orbital Mechanics

Fundamentals of Gravitation and Universal Gravitation Law

• Historical Discovery and Phenomenological Context:

  • Gravitation represents the universal attractive force existing between any two arbitrary masses in the universe.

  • Discovered by Sir Isaac Newton in 1665 following his observation of an apple falling from a tree toward the Earth's surface.

  • Gravitation is fundamentally the weakest of the four fundamental interactions of nature (gravitational, electromagnetic, weak nuclear, and strong nuclear).

  • Despite its minimal relative magnitude, it is the most dominant large-scale force in the cosmos, responsible for initiating the gravitational collapse that produces stars and governing the large-scale structure and dynamics of galaxies, planetary systems, and the universe.

• Newton's Law of Universal Gravitation:

  • Statement: Every particle in the universe attracts every other particle with a force that is directly proportional to the product of their masses and inversely proportional to the square of the distance separating their centers.

  • Scalar Form: For two point masses $m_1$ and $m_2$ separated by a distance $r$, the magnitude of the mutual gravitational attraction is given by:     $F = \frac{G m_1 m_2}{r^2}$

  • Universal Gravitational Constant ($G$): The proportionality factor $G$ has the value:     $G = 6.674 \times 10^{-11}\,\text{N}\cdot\text{m}^2/\text{kg}^2$

  • Historical Evaluation of $G$: Newton did not explicitly formulate the constant $G$ in his original expressions; the universal gravitational constant was first experimentally measured in a laboratory setting in the late eighteenth century by Henry Cavendish utilizing a sensitive torsional balance.

  • Nature of $G$: The constant $G$ is completely independent of the medium separating the masses, the physical state, temperature, chemical composition, or shape of the interacting bodies.

• Vector Formulation and Third Law Consistency:

  • Defining a unit vector $\hat{\mathbf{r}}_{12}$ directed along the displacement line from particle 1 to particle 2, the gravitational force $\mathbf{F}_{12}$ exerted by particle 1 on particle 2 is formulated as:     $\mathbf{F}_{12} = -\frac{G m_1 m_2}{r^2} \hat{\mathbf{r}}_{12}$

  • The negative sign explicitly denotes that the force is attractive, pulling particle 2 inward toward particle 1.

  • By Newton's third law of motion, the force exerted by particle 2 on particle 1 ($\mathbf{F}_{21}$) constitutes an action-reaction pair with $\mathbf{F}_{12}$:     $\mathbf{F}_{21} = -\mathbf{F}_{12}$     $\mathbf{F}_{21} = \frac{G m_1 m_2}{r^2} \hat{\mathbf{r}}_{12}$

• Applicability, Boundary Constraints, and the Inverse-Square Law:

  • Point-Mass Constraint: The formulation strictly applies to point particles. For arbitrary extended bodies of finite, irregular dimensions, the formula cannot be directly evaluated because no unique, single spatial separation distance $r$ exists across all differential mass elements.

  • Spherical Symmetry Exception: If interacting extended bodies are uniform spheres or spherically symmetric shells, the spatial separation $r$ is defined precisely as the distance between their geometric centers. A uniform mass sphere acts gravitationally on all external points as though its entire mass were concentrated at its geometric center.

  • Weakness in Daily Scenarios: Gravitational attraction between typical everyday objects is extraordinarily small ($\sim 10^{-7}\,\text{N}$ or less) and undetectable without precision instrumentation. Significant gravitational phenomena manifest only when at least one interacting body possesses astronomical mass (e.g., interaction between the human body and the Earth, or between the Earth and the Sun). Gravitational interactions between subatomic particles (such as an electron and proton in an atom) are entirely negligible compared to electromagnetic forces.

• Distinction Between Gravitation and Gravity:

  • Gravitation: The general attractive interaction occurring between any pair of bodies in the universe.

  • Gravity: The specific manifestation of gravitation representing the Earth's attractive pull on a body situated on or near its surface, directed radially inward toward the Earth's center.

  • Weight ($W$): The measure of the force of gravity exerted on a body of mass $m$ by the Earth:     $W = mg$

  • Comparison of Symbols $g$ and $G$:

  • Acceleration due to gravity ($g$): A local vector quantity representing the free-fall acceleration near a planet. At the Earth's surface at mean sea level, its standard average magnitude is $9.8\,\text{m/s}^2$. Its value varies with altitude, depth, latitude, and planetary mass.

  • Universal gravitational constant ($G$): A universal scalar constant having the identical value ($6.674 \times 10^{-11}\,\text{N}\cdot\text{m}^2/\text{kg}^2$) across all spatial regions and reference frames in the universe.

Kepler's Laws of Planetary Motion

• Historical Shift from Geocentric to Heliocentric Paradigms:

  • Johannes Kepler established three empirical laws of planetary motion, published sequentially in 1609 (First and Second Laws) and 1619 (Third Law).

  • Kepler's empirical calculations, founded upon the meticulous observational data of Tycho Brahe, decisively overturned the classical Ptolemaic geocentric model (which posited Earth at the motionless center of the cosmos orbited by the Sun, planets, and celestial spheres).

  • Kepler, alongside Nicolaus Copernicus and Galileo Galilei, established the heliocentric (Sun-centered) architecture of the solar system.

  • Solar System Context: The Sun's planetary system comprises planets orbiting at varying distances: Mercury (innermost, smallest, hottest), Venus (brightest, closest to Earth), Earth, Mars, Jupiter (largest, hosting 12 prominent moons in traditional classical counts), Saturn, Uranus, Neptune, and Pluto.

• Kepler's First Law (The Law of Orbits):

  • Statement: Every planet moves in an elliptical orbit around the Sun, with the Sun situated at one of the two foci of the ellipse.

  • Elliptical Geometry:

  • An ellipse is geometrically defined as the locus of all coplanar points whose distances $r_1$ and $r_2$ from two fixed focal points $F_1$ and $F_2$ sum to a constant value equal to the major axis length:       $r_1 + r_2 = 2a$

  • Major Axis: The longest line passing through both foci and the geometric center, possessing total length $2a$ (where $a$ is the semi-major axis).

  • Minor Axis: The shortest line passing through the geometric center perpendicular to the major axis, possessing total length $2b$ (where $b$ is the semi-minor axis).

  • The Sun resides at focal position $F_1$; the alternate focus $F_2$ is completely empty space.

  • Eccentricity ($e$): Dimensionless parameter measuring orbital elongation, bounded by $0 \le e < 1$. When $e = 0$, the ellipse degenerates to a circle ($a = b$); when $e \to 1$, the ellipse elongates into a line segment.

  • Earth's Orbital Eccentricity: Approximately $e \approx 0.0167$, establishing that Earth's orbit is nearly circular with minimal deviation.

  • Orbital Extremes:

  • Perihelion: The point of closest approach between the planet and the Sun. For an orbit around Earth, this position is designated perigee. The perihelion distance is:       $r_p = a - ea = a(1 - e)$

  • Aphelion: The point of maximum separation between the planet and the Sun. For an orbit around Earth, this position is designated apogee. The aphelion distance is:       $r_a = a + ea = a(1 + e)$

  • Major axis sum relation:       $r_p + r_a = a(1 - e) + a(1 + e) = 2a$

• Kepler's Second Law (The Law of Areas):

  • Statement: The radius vector joining the Sun to a planet sweeps out equal planar areas in equal intervals of time; that is, the areal velocity of the planet is strictly constant.

  • Kinematic Implication: A planet travels with minimum linear speed when positioned farthest from the Sun at aphelion, and accelerates to its maximum linear speed when nearest to the Sun at perihelion.

  • Derivation and Angular Momentum Conservation Proof:

  • Consider a planet of mass $m$ executing planar motion around the Sun at the origin in polar coordinates $(r, \theta)$.

  • In an infinitesimal time interval $dt$, the radius vector rotates through an angular displacement $d\theta$, traversing an arc length $r\,d\theta$.

  • The infinitesimal triangular area $dA$ swept out by the radius vector is:       $dA = \frac{1}{2} (r) (r\,d\theta) = \frac{1}{2} r^2 d\theta$

  • Dividing by $dt$ yields the areal velocity:       $\frac{dA}{dt} = \frac{1}{2} r^2 \frac{d\theta}{dt} = \frac{1}{2} r^2 \omega$

  • The orbital angular momentum $L$ of the planet relative to the Sun is:       $L = m r^2 \frac{d\theta}{dt} = m r^2 \omega$

  • Substituting $r^2 \frac{d\theta}{dt} = \frac{L}{m}$ into the areal velocity expression:       $\frac{dA}{dt} = \frac{L}{2m}$

  • Gravitational force is a central force directed strictly along the radial position vector $\mathbf{r}$. The external torque $\boldsymbol{\tau}$ acting on the planet about the Sun is identically zero:       $\boldsymbol{\tau} = \mathbf{r} \times \mathbf{F} = \mathbf{0}$

  • By rotational dynamics, $\frac{d\mathbf{L}}{dt} = \boldsymbol{\tau} = \mathbf{0}$, proving that the angular momentum $L$ is a conserved constant of motion.

  • Because $L$ and $m$ are strictly constant, the areal velocity is invariant:       $\frac{dA}{dt} = \text{constant}$

• Kepler's Third Law (The Law of Periods):

  • Statement: The square of the orbital period of revolution ($T$) of any planet around the Sun is directly proportional to the cube of the semi-major axis ($a$) of its elliptical orbit:     $T^2 \propto a^3 \implies T^2 = K a^3$

  • Mass Independence of Constant $K$: The constant of proportionality $K$ depends exclusively on the mass of the central attracting body (the Sun) and is independent of the planet's mass. Thus, $K$ is identical for all planets orbiting the Sun, but adopts a different numerical value for satellites orbiting Earth or other primary bodies.

  • Extension to Circular Orbits: For a circular orbit of radius $r$, the semi-major axis equals the radius ($a = r$), maintaining the exact functional relationship $T^2 \propto r^3$.

Dynamic Applications and Problem Solutions of Keplerian Motion

• Orbital Scaling and Planetary Distance Determinations:

  • Orbital Period of Jupiter (Illustration 5.1):

  • Given: Earth orbital period $T_E = 1.0\,\text{year}$, Earth orbital radius $r_E = 1.5 \times 10^{11}\,\text{m}$, Jupiter period $T_J = 11.6\,\text{years}$.

  • From Kepler's Third Law: $\left(\frac{r_J}{r_E}\right)^3 = \left(\frac{T_J}{T_E}\right)^2$.

  • Solving for Jupiter's orbital radius $r_J$:       $r_J = r_E \left(\frac{T_J}{T_E}\right)^{2/3} = (1.5 \times 10^{11}\,\text{m}) (11.6)^{2/3} = 7.68 \times 10^{11}\,\text{m}$

  • Orbital Period of Mars (Illustration 5.2):

  • Given: Mean Sun-Mars distance is $1.524$ times the Sun-Earth distance ($r_M = 1.524\,r_E$), Earth period $T_E = 1\,\text{year}$.

  • Applying period proportionality:       $T_M = T_E \left(\frac{r_M}{r_E}\right)^{3/2} = (1\,\text{year}) (1.524)^{3/2} = 1.88\,\text{years}$

• Orbital Extremes Kinematics (Conservation of Angular Momentum):

  • Speed at Perihelion versus Aphelion (Illustration 5.3):

  • At perihelion ($P$) and aphelion ($A$), the velocity vectors $\mathbf{v}_p$ and $\mathbf{v}_a$ are strictly perpendicular to their respective radius vectors $\mathbf{r}_p$ and $\mathbf{r}_a$.

  • Angular momentum conservation requires:       $L = m v_p r_p = m v_a r_a \implies \frac{v_p}{v_a} = \frac{r_a}{r_p}$

  • Because $r_a > r_p$, the speed at perihelion is greater than at aphelion ($v_p > v_a$).

  • Area Traversals: Given that the bounded orbital area $S_{BAC} > S_{CPB}$, and because areal velocity $\frac{dA}{dt}$ is constant, the planet requires a strictly longer time to traverse arc $BAC$ than arc $CPB$.

  • Speed Ratio for Specified Eccentricity (Illustration 5.4):

  • Consider an elliptical orbit with semi-major axis $a$ and eccentricity $e = \frac{1}{2}$.

  • Perihelion distance: $r_p = a(1 - e) = a\left(1 - \frac{1}{2}\right) = \frac{a}{2}$.

  • Aphelion distance: $r_a = a(1 + e) = a\left(1 + \frac{1}{2}\right) = \frac{3a}{2}$.

  • Given speed at perihelion is $v_1$, let speed at aphelion be $v_2$.

  • Conservation of angular momentum yields:       $m v_1 r_p = m v_2 r_a \implies v_2 \left(\frac{3a}{2}\right) = v_1 \left(\frac{a}{2}\right) \implies v_2 = \frac{v_1}{3}$

• Concept Application Exercise Solutions (Kepler's Laws):

  • Orbit Radii Ratio and Revolutions Count (CAE 5.1, Q1):

  • Orbit radii are $r_A = r$ and $r_B = 4r$.

  • Period ratio: $\frac{T_A}{T_B} = \left(\frac{r}{4r}\right)^{3/2} = \left(\frac{1}{4}\right)^{3/2} = \frac{1}{8}$, meaning $T_B = 8 T_A$.

  • In a quarter year of planet B (time $\Delta t = \frac{T_B}{4}$):       $\text{Number of revolutions of A} = \frac{\Delta t}{T_A} = \frac{T_B / 4}{T_B / 8} = 2$

  • Areal Velocity Equivalence (CAE 5.1, Q2):

  • Areas swept are $A_1 = A$ and $A_2 = 2A$ during intervals $\Delta t_1$ and $\Delta t_2$.

  • Constancy of areal velocity dictates:       $\frac{A}{\Delta t_1} = \frac{2A}{\Delta t_2} \implies \frac{\Delta t_1}{\Delta t_2} = \frac{1}{2}$

  • Orbit Periods from Perihelion/Aphelion Geometry (CAE 5.1, Q3):

  • Semi-major axis for an orbit with extremes $r_p$ and $r_a$ is $a = \frac{r_p + r_a}{2}$.

  • For orbit 1 with extremes $a$ and $b$, $a_1 = \frac{a+b}{2}$; for orbit 2 with extremes $a$ and $c$, $a_2 = \frac{a+c}{2}$.

  • Ratio of periods:       $\frac{T_1}{T_2} = \left(\frac{a_1}{a_2}\right)^{3/2} = \left(\frac{a+b}{a+c}\right)^{3/2}$

Extended Mass Distributions and Superposition Principles

• Principle of Superposition in Gravitation:

  • The net gravitational force exerted on an arbitrary target particle 1 by an ensemble of $N$ distinct external particles is the vector sum of all individual forces acting independently:     $\mathbf{F}_{1,\text{net}} = \sum_{i=2}^{N} \mathbf{F}_{1i} = \mathbf{F}_{12} + \mathbf{F}_{13} + \mathbf{F}_{14} + \dots + \mathbf{F}_{1N}$

  • Pairwise Independence: The presence of third-party masses has zero physical influence on the individual mutual gravitational interaction between any specified pair of masses.

• Superposition Demonstration (Illustration 5.5):

  • Point masses $m$ and $M$ separated by distance $r$ experience mutual attractive force $F = \frac{GMm}{r^2}$.

  • When an additional point mass $3m$ is placed alongside $m$:

  • Force on $M$ due strictly to $m$ remains invariant: $F = \frac{GMm}{r^2}$.

  • Total gravitational force on $M$ becomes:       $F' = \frac{G(m + 3m)M}{r^2} = \frac{4GMm}{r^2} = 4F$

• Centripetal Balance and Solar Mass Determination (Illustration 5.6):

  • Earth mass $m = 6 \times 10^{24}\,\text{kg}$, orbital radius $r = 1.5 \times 10^{11}\,\text{m}$, orbital period $T = 3.2 \times 10^7\,\text{s}$.

  • Angular velocity:     $\omega = \frac{2\pi}{T} = \frac{2\pi}{3.2 \times 10^7\,\text{s}} = 1.96 \times 10^{-7}\,\text{rad/s}$

  • Centripetal force:     $F = m r \omega^2 = (6 \times 10^{24}\,\text{kg}) (1.5 \times 10^{11}\,\text{m}) (1.96 \times 10^{-7}\,\text{rad/s})^2 = 3.46 \times 10^{22}\,\text{N}$

  • Centripetal force is supplied entirely by solar gravity:     $\frac{G M_{\text{Sun}} m}{r^2} = F \implies M_{\text{Sun}} = \frac{F r^2}{G m} = \frac{(3.46 \times 10^{22}\,\text{N}) (1.5 \times 10^{11}\,\text{m})^2}{(6.67 \times 10^{-11}\,\text{N}\cdot\text{m}^2/\text{kg}^2) (6 \times 10^{24}\,\text{kg})} = 1.95 \times 10^{30}\,\text{kg}$

• Derivation of Kepler's Third Law from Newton's Gravity (Illustration 5.7):

  • For a planet of mass $m$ in a circular orbit of radius $r$ around a central mass $M$:     $\frac{m v^2}{r} = \frac{GMm}{r^2} \implies v = \sqrt{\frac{GM}{r}}$

  • The orbital period is:     $T = \frac{2\pi r}{v} = \frac{2\pi r}{\sqrt{GM/r}} = \frac{2\pi r^{3/2}}{\sqrt{GM}}$

  • Squaring both sides yields:     $T^2 = \left(\frac{4\pi^2}{GM}\right) r^3$

  • Because $\frac{4\pi^2}{GM}$ is invariant for any satellite orbiting central mass $M$, $T^2 \propto r^3$.

• Three-Body Dynamic Equilibrium in Equilateral Configuration (Illustration 5.8):

  • Three identical bodies of mass $m$ are positioned at the vertices of an equilateral triangle of side length $r$, revolving in a circumscribed circle about their common center of mass while maintaining relative configuration.

  • Force on vertex $C$ due to vertex $A$: $F_1 = \frac{Gm^2}{r^2}$ directed along $CA$.

  • Force on vertex $C$ due to vertex $B$: $F_2 = \frac{Gm^2}{r^2}$ directed along $CB$.

  • The vectors are inclined at angle $60^\circ$. Net resultant attractive force along the symmetry axis $CD$ is:     $F_{\text{net}} = \sqrt{F_1^2 + F_2^2 + 2F_1 F_2 \cos(60^\circ)} = \sqrt{3} F_1 = \sqrt{3} \frac{Gm^2}{r^2}$

  • The radius of the circumscribed orbit is $R_{\text{circ}} = OC = \frac{r/2}{\cos(30^\circ)} = \frac{r}{\sqrt{3}}$.

  • Centripetal acceleration balance:     $\frac{m v^2}{R_{\text{circ}}} = F_{\text{net}} \implies \frac{m v^2}{r/\sqrt{3}} = \sqrt{3} \frac{Gm^2}{r^2} \implies \frac{\sqrt{3} m v^2}{r} = \sqrt{3} \frac{Gm^2}{r^2} \implies v = \sqrt{\frac{Gm}{r}}$

• Three-Particle Collinear Equilibrium Search (Illustration 5.9):

  • Particles of mass $m$ and $2m$ reside at positions $x = 0$ and $x = d$. A third particle of mass $m$ is positioned at coordinate $x$.

  • Case 1: Net force on the particle at $x = 0$ is directed rightward and equals $\frac{6Gm^2}{d^2}$:     $F_{\text{net}} = \frac{Gm(m)}{x^2} + \frac{Gm(2m)}{d^2} = \frac{6Gm^2}{d^2} \implies \frac{Gm^2}{x^2} + \frac{2Gm^2}{d^2} = \frac{6Gm^2}{d^2}$     $\frac{1}{x^2} = \frac{4}{d^2} \implies x = \frac{d}{2}$

  • Case 2: Net force on the particle at $x = 0$ is directed leftward with magnitude $\frac{6Gm^2}{d^2}$ (placing third mass at negative coordinate $x = -y$):     $\frac{Gm^2}{y^2} - \frac{2Gm^2}{d^2} = \frac{6Gm^2}{d^2} \implies \frac{1}{y^2} = \frac{8}{d^2} \implies y = \frac{d}{2\sqrt{2}}$

• Spherical Cavity Analysis via Negative Mass Superposition (Illustration 5.10):

  • A solid lead sphere of total initial mass $M$ and radius $R$ has a spherical cavity carved within it. The cavity boundary passes through the sphere's origin $O$ and touches the outer edge; its diameter is $R$, so its radius is $R' = R/2$, with center $O'$ at displacement $R/2$.

  • A test mass $m$ is placed at an external point $P$ along the line of centers at distance $r$ from $O$.

  • Intact Sphere Attraction: $F = \frac{GMm}{r^2}$ directed along $PO$.

  • Cavity Removed Mass: The volume of the carved sphere is $\left(\frac{1}{2}\right)^3 = \frac{1}{8}$ of the whole, giving mass $M' = \frac{M}{8}$. Center of cavity is located at distance $r - \frac{R}{2}$ from $P$.

  • Removed Mass Attraction: $F' = \frac{G(M/8)m}{(r - R/2)^2} = \frac{GMm}{8(r - R/2)^2}$ along $PO$.

  • Superposed Net Force:     $F_{\text{cavity}} = F - F' = \frac{GMm}{r^2} - \frac{GMm}{8(r - R/2)^2} = \frac{GMm}{r^2} \left[1 - \frac{1}{8\left(1 - \frac{R}{2r}\right)^2}\right]$

• Gravitational Attraction of a Continuous Thin Rod (Illustration 5.11):

  • Consider a thin uniform rod of total mass $M$ and length $L$. A particle of mass $m$ is placed collinearly at distance $d$ from the near end of the rod.

  • Linear mass density of the rod: $\lambda = \frac{M}{L}$.

  • Consider an infinitesimal element $dx$ at distance $x$ from the test mass $m$, carrying mass $dm = \left(\frac{M}{L}\right) dx$.

  • Infinitesimal force:     $dF = \frac{Gm\,dm}{x^2} = \frac{GMm}{L} \frac{dx}{x^2}$

  • Integrating across the entire rod from $x = d$ to $x = d + L$:     $F = \int_{d}^{d+L} \frac{GMm}{L} \frac{dx}{x^2} = \frac{GMm}{L} \left[-\frac{1}{x}\right]_{d}^{d+L} = \frac{GMm}{L} \left(\frac{1}{d} - \frac{1}{d+L}\right) = \frac{GMm}{L} \left(\frac{L}{d(d+L)}\right) = \frac{GMm}{d(d+L)}$

Gravitational Field, Field Intensity, and Continuous Geometries

• Gravitational Field Concept:

  • A mass modifies the physical properties of the space surrounding it, establishing a gravitational field.

  • When a secondary mass enters this region, it experiences an attractive force mediated directly by the local field.

  • Intensity of Gravitational Field ($\mathbf{E}$): The gravitational force experienced per unit test mass placed at a given point in space:     $\mathbf{E} = \lim_{m_0 \to 0} \frac{\mathbf{F}}{m_0}$

  • Field of a Point Mass $M$ at Distance $x$:     $\mathbf{E} = -\frac{GM}{x^2} \hat{\mathbf{r}}$

  • The negative sign indicates the field is directed radially inward toward the field-generating source mass.

  • Properties of Field Intensity:

  • Equivalence to Acceleration: For a free test mass $m_0$, $F = m_0 a \implies a = \frac{F}{m_0} = E$. The field intensity at any point is numerically and dimensionally equal to the free-fall acceleration experienced by a test particle at that coordinate.

  • Terrestrial Surface Field: At the Earth's surface ($x = R$), $E = \frac{GM}{R^2} = g$.

  • Units and Dimensions: SI unit is $\text{N/kg}$ or $\text{m/s}^2$; CGS unit is $\text{dyne/g}$ or $\text{cm/s}^2$. Dimensional formula: $[\text{M}^0 \text{L} \text{T}^{-2}]$.

• Symmetry Method in Field Intensity Calculations:

  • Equilateral Triangle ($n = 3$): Identical point masses at vertices generate three equal magnitude field vectors at the centroid spaced by $120^\circ$. Vector sum: $\sum \mathbf{E} = \mathbf{0}$.

  • Square ($n = 4$): Opposite pairs at diagonal vertices produce equal and opposite field vectors that cancel identically.

  • Regular Polygons and Polyhedra: For any regular $n$-gon or regular solid with identical masses at all vertices, rotational symmetry ensures the net gravitational field at the geometric center is zero:     $\mathbf{E}_{\text{center}} = \mathbf{0}$

  • Symmetrical Doubling (Illustration 5.14): If masses $m$ are at vertices $B$ and $C$ of an equilateral triangle, and mass $2m$ is at vertex $A$ with centroid distances $AG = BG = CG = 1\,\text{m}$:

  • Decompose vertex $A$ into $m + m$. The symmetric set of three masses $m$ yields zero field at centroid $G$.

  • The net field is due entirely to the extra mass $m$ at $A$: $\mathbf{E}_{G} = \frac{Gm}{(1)^2} \hat{\mathbf{r}}_{GA} = Gm\,\hat{\mathbf{r}}_{GA}$.

  • Net force on a test mass $2m$ at $G$ is $\mathbf{F} = (2m)\mathbf{E}_G = 2Gm^2\,\hat{\mathbf{r}}_{GA}$.

• Newton's Shell Theorem:

  • Exterior Theorem: A uniform spherical shell of mass $M$ attracts an external particle as if the entire mass of the shell were concentrated as a point mass at its geometric center.

  • Interior Theorem: A uniform spherical shell exerts zero net gravitational force on any particle situated anywhere in its interior cavity, regardless of position.

• Field of a Thin Uniform Spherical Shell of Radius $R$:

  • Outer points ($r > R$): $E_{\text{out}} = \frac{GM}{r^2}$

  • Surface points ($r = R$): $E_s = \frac{GM}{R^2}$

  • Inner cavity points ($r < R$): $E_{\text{in}} = 0$

• Field of a Uniform Solid Sphere of Mass $M$ and Radius $R$:

  • External points ($x \ge R$):     $E = \frac{GM}{x^2}$

  • Surface point ($x = R$):     $E_s = \frac{GM}{R^2}$

  • Internal points ($x < R$):

  • Draw a concentric sphere of radius $x$. By the Shell Theorem, the spherical shell region between $x$ and $R$ produces zero net field at $x$.

  • Enclosed mass $m(x)$ of the inner sphere of radius $x$ is:       $m(x) = M \left(\frac{\frac{4}{3}\pi x^3}{\frac{4}{3}\pi R^3}\right) = M \frac{x^3}{R^3}$

  • The interior field is governed solely by this enclosed mass:       $E_{\text{in}} = \frac{G m(x)}{x^2} = \frac{G}{x^2} \left(M \frac{x^3}{R^3}\right) = \frac{GMx}{R^3}$

  • Expressed in terms of uniform density $\rho$:       $E_{\text{in}} = \frac{G \left(\frac{4}{3}\pi x^3 \rho\right)}{x^2} = \frac{4}{3}\pi G \rho x$

  • Inside the sphere, the field scales directly linearly with radial distance ($E \propto x$), achieves its global maximum at the surface ($x = R$), and decays as an inverse-square law outside ($E \propto \frac{1}{x^2}$).

• Gravitational Field of a Uniform Circular Ring on its Axial Line:

  • Ring of mass $M$, radius $R$, centered at the origin in the $yz$-plane. Field evaluated at axial point $P(x, 0, 0)$.

  • Element of arc length $dl$ has mass $dm = \frac{M}{2\pi R} dl$.

  • Distance from element to point $P$ is $\sqrt{x^2 + R^2}$.

  • Elemental field magnitude:     $dE = \frac{G\,dm}{x^2 + R^2} = \frac{G M dl}{2\pi R(x^2 + R^2)}$

  • Symmetry Cancellation: Perpendicular components $dE\sin\theta$ cancel upon integration over the full ring. Axial components $dE\cos\theta$ sum constructively, where $\cos\theta = \frac{x}{\sqrt{x^2 + R^2}}$.

  • Total axial field:     $E = \int dE\cos\theta = \int_0^{2\pi R} \frac{GMx\,dl}{2\pi R(x^2 + R^2)^{3/2}} = \frac{GMx}{2\pi R(x^2 + R^2)^{3/2}} (2\pi R) = \frac{GMx}{(x^2 + R^2)^{3/2}}$

• Gravitational Field of a Uniform Disc on its Axial Line:

  • Disc of mass $M$, radius $R$, area $\pi R^2$, surface density $\sigma = \frac{M}{\pi R^2}$.

  • Decompose into concentric rings of radius $r$ and radial width $dr$. Annular area $dA = 2\pi r\,dr$, mass $dm = \frac{M}{\pi R^2} 2\pi r\,dr = \frac{2Mr\,dr}{R^2}$.

  • Axial field produced by this annular ring:     $dE = \frac{G(dm)x}{(x^2 + r^2)^{3/2}} = \frac{2GMx\,r\,dr}{R^2(x^2 + r^2)^{3/2}}$

  • Integrating across the disc radius from $r = 0$ to $r = R$:     $E = \frac{2GMx}{R^2} \int_0^R \frac{r\,dr}{(x^2 + r^2)^{3/2}}$

  • Using substitution $z = x^2 + r^2$, $dz = 2r\,dr$:     $\int_0^R \frac{r\,dr}{(x^2 + r^2)^{3/2}} = \left[-\frac{1}{\sqrt{x^2 + r^2}}\right]_0^R = \frac{1}{x} - \frac{1}{\sqrt{x^2 + R^2}}$

  • Resulting field magnitude along the axis toward the disc center:     $E = \frac{2GM}{R^2} \left(1 - \frac{x}{\sqrt{x^2 + R^2}}\right) = \frac{2GM}{R^2}(1 - \cos\theta)$     where $\theta$ represents the semi-vertical angle subtended by the disc rim at point $P$.

• Field of a Hollow Open Hemisphere at Center of Base (Illustration 5.15):

  • Open hollow hemisphere of mass $M$, radius $R$. Slice into horizontal rings at angle $\theta$ with angular width $d\theta$.

  • Ring radius is $R\sin\theta$, distance to base center is $R$, mass $dm = \frac{M}{2\pi R^2} (2\pi R\sin\theta)(R\,d\theta) = M\sin\theta\,d\theta$.

  • Axial component at the base center:     $dE = \frac{G\,dm}{R^2} \cos\theta = \frac{GM}{R^2} \sin\theta \cos\theta\,d\theta$

  • Integrating from $\theta = 0$ to $\theta = \pi/2$:     $E = \frac{GM}{R^2} \int_0^{\pi/2} \sin\theta \cos\theta\,d\theta = \frac{GM}{R^2} \left[\frac{\sin^2\theta}{2}\right]_0^{\pi/2} = \frac{GM}{2R^2}$

• Gravitational Field of a Semicircular Rod (Illustration 5.16):

  • Thin rod of mass $M$, length $L$, bent into a semicircle of radius $R$ ($L = \pi R \implies R = L/\pi$). Linear density $\lambda = \frac{M}{L} = \frac{M}{\pi R}$.

  • Element at angle $\theta$ subtends $d\theta$, mass $dm = \lambda R\,d\theta$.

  • Transverse components cancel by bilateral symmetry; radial components along the axis of symmetry sum constructively:     $E = \int_{-\pi/2}^{\pi/2} \frac{G\,dm}{R^2} \cos\theta = \frac{G\lambda}{R} \int_{-\pi/2}^{\pi/2} \cos\theta\,d\theta = \frac{2G\lambda}{R} = \frac{2GM}{\pi R^2} = \frac{2\pi GM}{L^2}$

  • Gravitational force on particle of mass $m$ at center of curvature:     $F = mE = \frac{2\pi GMm}{L^2}$

  • For a complete closed circular ring, symmetric contributions cancel completely at the geometric center, yielding $E_{\text{center}} = 0$ and $F = 0$.

• Uniform Ring Above a Uniform Solid Sphere (Illustration 5.17):

  • Ring of mass $m$ and radius $r$ positioned at vertical height $z = r\sqrt{3}$ directly above a sphere of mass $M$ and radius $r$.

  • Outside the sphere, its field is identical to a point mass $M$ at its center.

  • Distance from sphere center to any point on the ring is $d = \sqrt{(r\sqrt{3})^2 + r^2} = 2r$.

  • Sphere field magnitude at the ring element: $E = \frac{GM}{(2r)^2} = \frac{GM}{4r^2}$.

  • Axial projection factor: $\cos\theta = \frac{r\sqrt{3}}{2r} = \frac{\sqrt{3}}{2}$.

  • Net downward vertical gravitational force on the entire ring of mass $m$:     $F = \int dm\,E\cos\theta = m \left(\frac{GM}{4r^2}\right) \left(\frac{\sqrt{3}}{2}\right) = \frac{\sqrt{3}GMm}{8r^2}$

• Sphere with Linearly Increasing Density Profile (Illustration 5.18):

  • Density profile: $\rho(r) = \rho_0 \left(\frac{r}{a}\right)$ for $r \le a$, where $\rho_0$ is surface density.

  • Enclosed mass within radius $x$ ($x < a$):     $m(x) = \int_0^x \rho(r) 4\pi r^2 dr = \int_0^x \rho_0 \frac{r}{a} 4\pi r^2 dr = \frac{4\pi\rho_0}{a} \int_0^x r^3 dr = \frac{\pi\rho_0 x^4}{a}$

  • Interior field at $x < a$:     $E_{\text{in}}(x) = \frac{G m(x)}{x^2} = \frac{G}{x^2} \left(\frac{\pi\rho_0 x^4}{a}\right) = \frac{\pi G \rho_0 x^2}{a}$

  • Total mass of sphere ($x = a$):     $M = \frac{\pi\rho_0 a^4}{a} = \pi\rho_0 a^3$

  • Exterior field at $x \ge a$:     $E_{\text{out}}(x) = \frac{GM}{x^2} = \frac{\pi G \rho_0 a^3}{x^2}$

Acceleration Due to Gravity and Environmental Variations

• Fundamental Formulation of Acceleration Due to Gravity ($g$):

  • By Newton's second law, $F = mg$. Equating to gravitational surface attraction $F = \frac{GMm}{R^2}$ yields:     $g = \frac{GM}{R^2}$

  • Expression in Terms of Mean Planetary Density ($\rho$):     $M = \frac{4}{3}\pi R^3 \rho \implies g = \frac{G\left(\frac{4}{3}\pi R^3 \rho\right)}{R^2} = \frac{4}{3}\pi \rho G R$

  • Invariance: Acceleration due to gravity is independent of the mass, geometry, or material composition of the falling object; it depends solely on central planetary parameters ($M$, $R$, $\rho$).

  • Standard Terrestrial Sea-Level Values: $g = 9.80\,\text{m/s}^2 = 981\,\text{cm/s}^2 = 32\,\text{ft/s}^2$.

  • Terrestrial Mean Density Calculation (Illustration 5.19):     $\rho_E = \frac{3g}{4\pi G R} = \frac{3(9.80\,\text{m/s}^2)}{4\pi (6.674 \times 10^{-11}\,\text{N}\cdot\text{m}^2/\text{kg}^2)(6.37 \times 10^6\,\text{m})} = 5.50 \times 10^3\,\text{kg/m}^3$

  • Lunar Mass Evaluation (Illustration 5.20):

  • Given $g_m = 1.67\,\text{m/s}^2$, $R_m = 1.74 \times 10^6\,\text{m}$:       $M_m = \frac{g_m R_m^2}{G} = \frac{(1.67\,\text{m/s}^2)(1.74 \times 10^6\,\text{m})^2}{6.67 \times 10^{-11}\,\text{N}\cdot\text{m}^2/\text{kg}^2} = 7.6 \times 10^{22}\,\text{kg}$

• Variation of $g$ with Altitude (Height $h$):

  • At height $h$ above the surface ($r = R + h$):     $g' = \frac{GM}{(R + h)^2} = \frac{GM}{R^2 \left(1 + \frac{h}{R}\right)^2} = g \left(\frac{R}{R + h}\right)^2 = g \left(1 + \frac{h}{R}\right)^{-2}$

  • Small Altitude Approximation ($h \ll R$): Expanding via binomial series:     $g' \approx g \left(1 - \frac{2h}{R}\right)$

  • Fractional and Percentage Decay:

  • Absolute decrease: $\Delta g = g - g' = \frac{2hg}{R}$

  • Fractional decrease: $\frac{\Delta g}{g} = \frac{2h}{R}$

  • Percentage decrease: $\frac{\Delta g}{g} \times 100\% = \frac{2h}{R} \times 100\%$

  • Significant Altitude Example (Illustration 5.21): Height where $g$ drops by $36\%$ ($g' = 0.64g$):     $0.64g = g\left(\frac{R}{R+h}\right)^2 \implies \frac{R}{R+h} = 0.8 \implies 8R + 8h = 10R \implies h = \frac{R}{4} = \frac{6400\,\text{km}}{4} = 1600\,\text{km}$

• Variation of $g$ with Depth ($d$):

  • At depth $d$ below the surface, the effective radius is $R - d$. Assuming uniform planetary density $\rho$:     $g_d = \frac{4}{3}\pi \rho G (R - d)$

  • Taking the ratio with surface gravity $g = \frac{4}{3}\pi \rho G R$:     $g_d = g \left(1 - \frac{d}{R}\right)$

  • Center of the Earth: At $d = R$, $g_{\text{center}} = g\left(1 - \frac{R}{R}\right) = 0$.

  • Fractional and Percentage Decrease:

  • Absolute decrease: $\Delta g = g - g_d = \frac{dg}{R}$

  • Fractional decrease: $\frac{\Delta g}{g} = \frac{d}{R}$

  • Percentage decrease: $\frac{\Delta g}{g} \times 100\% = \frac{d}{R} \times 100\%$

  • Comparison of Gradients: For small vertical excursions ($h, d \ll R$), the rate of decrease of gravity with elevation above the Earth is precisely twice the rate of decrease with depth below the surface:     $\left(\frac{\Delta g}{g}\right)_{\text{altitude}} = 2 \left(\frac{\Delta g}{g}\right)_{\text{depth}}$

  • Mine Shaft Weight Loss (Illustration 5.22): Body taken to depth $d = 16\,\text{km}$ in Earth of radius $R = 6400\,\text{km}$:     $\frac{\Delta W}{W} \times 100\% = \frac{d}{R} \times 100\% = \frac{16\,\text{km}}{6400\,\text{km}} \times 100\% = 0.25\%$     At the center of the Earth ($d = R$), weight decreases by $100\%$, reaching exactly zero.

• Effect of the Earth's Ellipsoidal Shape:

  • The Earth is an oblate spheroid, flattened at the geographic poles and bulging at the equator.

  • Equatorial radius exceeds polar radius: $R_e - R_p \approx 21\,\text{km}$.

  • Surface gravity dependence:     $g = \frac{GM}{R^2} \implies g_p = \frac{GM}{R_p^2}, \quad g_e = \frac{GM}{R_e^2}$

  • Because $R_e > R_p$, acceleration due to gravity is minimum at the equator and maximum at the poles:     $g_p > g_e$

  • Numerical Difference at Sea Level: The polar acceleration exceeds equatorial acceleration by approximately $0.018\,\text{m/s}^2$ ($1.80\,\text{cm/s}^2$) due solely to geometrical oblateness.

• Effect of Earth's Axial Rotation:

  • Consider Earth rotating about its polar axis with angular speed $\omega$. At a location with latitude $\lambda$, a particle rotates in a circular path of radius $r = R\cos\lambda$.

  • In the rotating non-inertial terrestrial frame, the particle experiences a true gravitational pull $\mathbf{F}_g = mg$ toward the Earth center, and a centrifugal pseudo-force directed radially outward from the rotation axis:     $F_c = m r \omega^2 = m (R\cos\lambda) \omega^2$

  • The effective acceleration due to gravity $g'$ is the resultant of the two force fields:     $g' = \sqrt{g^2 + (R\omega^2\cos\lambda)^2 - 2g(R\omega^2\cos\lambda)\cos\lambda}$

  • Since $R\omega^2 \approx 0.034\,\text{m/s}^2 \ll g$, higher-order terms in $\left(\frac{R\omega^2}{g}\right)$ are negligible. Applying a first-order binomial expansion:     $g' = g \left(1 - \frac{2R\omega^2\cos^2\lambda}{g}\right)^{1/2} \approx g - R\omega^2\cos^2\lambda$

  • Latitudinal Extremes:

  • At the Poles ($\lambda = 90^\circ$, $\cos\lambda = 0$):       $g_p' = g - R\omega^2(0)^2 = g$

  • At the Equator ($\lambda = 0^\circ$, $\cos\lambda = 1$):       $g_e' = g - R\omega^2$

  • Pole-to-Equator Rotational Discrepancy:       $\Delta g_{\text{rot}} = g_p' - g_e' = R\omega^2$       With $R = 6.4 \times 10^6\,\text{m}$ and $\omega = \frac{2\pi}{86400}\,\text{rad/s}$, $R\omega^2 \approx 0.034\,\text{m/s}^2$.

  • Combined Effects: Gravity increases continuously from the equator to the poles due to both rotational centrifugal diminution and ellipsoidal flattening.

• Critical Rotational Speeds and Equatorial Weightlessness:

  • Condition for Apparent Weightlessness at the Equator ($g_e' = 0$):     $g - R\omega_c^2 = 0 \implies \omega_c = \sqrt{\frac{g}{R}}$

  • Numerical speed required:     $\omega_c = \sqrt{\frac{9.8\,\text{m/s}^2}{6.4 \times 10^6\,\text{m}}} \approx 1.24 \times 10^{-3}\,\text{rad/s}$

  • Percentage Increase Required (Illustration 5.25):

  • Current angular velocity: $\omega_0 = \frac{2\pi}{86400} \approx 7.27 \times 10^{-5}\,\text{rad/s}$.

  • Ratio: $\frac{\omega_c}{\omega_0} \approx 17$, representing an increase of approximately $1700\%$.

  • Rotational Stability of a Homogeneous Planet (Illustration 5.28):

  • For loose rocks to remain on a planet's equator without flying into space, the required centripetal acceleration cannot exceed gravity:       $\omega^2 R \le g \implies \left(\frac{2\pi}{T}\right)^2 R \le \frac{4}{3}\pi \rho G R$

  • Solving for the critical rotation period $T$:       $\frac{4\pi^2}{T^2} \le \frac{4}{3}\pi \rho G \implies T^2 \ge \frac{3\pi}{G\rho} \implies T \ge \sqrt{\frac{3\pi}{G\rho}}$