Skip to content

Gravitational Fields

Info: Board Coverage AQA Paper 2 | Edexcel CP3 | OCR (A) Paper 2 | CIE P4

1. Newton”s Law of Universal Gravitation

Section titled “1. Newton”s Law of Universal Gravitation”

Newton’s Law. Every point mass attracts every other point mass with a force directed along the line Joining them, whose magnitude is:

F=Gm1m2r2\boxed{F = \frac{Gm_1 m_2}{r^2}}

Where G=6.67×1011G = 6.67 \times 10^{-11} N m2^2 kg2^{-2} is the gravitational constant.

Geometric Origin of the Inverse Square Law

Section titled “Geometric Origin of the Inverse Square Law”

Consider a point mass emitting gravitational flux uniformly in all directions. The flux through a Sphere of radius rr is constant (by Gauss’s law for gravity):

gdA=4πGMenc\oint \mathbf{g} \cdot d\mathbf{A} = -4\pi G M_{\mathrm{enc}}

Since the surface area is 4πr24\pi r^2The flux density (field strength) must be g=GM/r2g = GM/r^2. The Force on a test mass mm is then F=mg=GMm/r2F = mg = GMm/r^2. This inverse square law is a direct geometric Consequence of flux conservation in three-dimensional space.

  • Always attractive (no negative mass exists).
  • Infinite range — the force extends to arbitrarily large distances, weakening as 1/r21/r^2.
  • Acts on all objects with mass — it is the universal force binding large-scale structures.
  • Extremely weak compared to electromagnetism: Fe/Fg1036F_e/F_g \sim 10^{36} for elementary particles.

The weakness of gravity means it is only dominant at macroscopic scales where the near-perfect Cancellation of positive and negative electric charges renders electromagnetic forces negligible.

Definition. The gravitational field strength g\mathbf{g} at a point is the force per unit mass on a Small test mass placed at that point:

g=Fm\boxed{\mathbf{g} = \frac{\mathbf{F}}{m}}

SI units: N kg1^{-1}Equivalent to m s2^{-2}.

For a point mass MM at distance rr:

g=GMr2\boxed{g = \frac{GM}{r^2}}

Directed radially inward towards MM.

Proof. From Newton’s law: F=GMm/r2F = GMm/r^2. Dividing by mm: g=F/m=GM/r2g = F/m = GM/r^2. \square

Theorem. A uniform spherical shell of mass MM and radius RR produces: (a) the same field as a Point mass MM at its centre for all external points (r>Rr \gt R), and (b) zero field at all internal Points (r<Rr \lt R).

Proof (external, r>Rr \gt R). Consider a test mass mm at distance rr from the centre. Divide the Shell into thin annular rings perpendicular to the line from the centre to mm. A ring at polar angle θ\theta has radius RsinθR\sin\thetaWidth RdθR\,d\thetaAnd mass:

dM=M4πR22πR2sinθdθ=M2sinθdθdM = \frac{M}{4\pi R^2} \cdot 2\pi R^2 \sin\theta\,d\theta = \frac{M}{2}\sin\theta\,d\theta

Every element of the ring is at distance s=r2+R22rRcosθs = \sqrt{r^2 + R^2 - 2rR\cos\theta} from mm. By the Axial symmetry of the ring, the transverse components of force cancel, leaving only the component Along the axis. The angle α\alpha between the force direction and the axis satisfies:

cosα=rRcosθs\cos\alpha = \frac{r - R\cos\theta}{s}

The axial force contribution from the ring is:

dF=GmdMs2cosα=GmM2(rRcosθ)sinθdθ(r2+R22rRcosθ)3/2dF = \frac{Gm\,dM}{s^2}\cos\alpha = \frac{GmM}{2}\cdot\frac{(r - R\cos\theta)\sin\theta\,d\theta}{(r^2 + R^2 - 2rR\cos\theta)^{3/2}}

Substitute u = \cos\theta$$du = -\sin\theta\,d\theta. When \theta = 0$$u = 1; when \theta = \pi$$u = -1:

F=GmM211(rRu)(du)(r2+R22rRu)3/2F = \frac{GmM}{2}\int_{1}^{-1}\frac{(r - Ru)(-du)}{(r^2 + R^2 - 2rRu)^{3/2}}

The integral evaluates to 2rr2R21r2(r2R2)1r\frac{2r}{r^2 - R^2} \cdot \frac{1}{r^2} \cdot (r^2 - R^2) \cdot \frac{1}{r} After careful algebra, giving:

F=GMmr2F = \frac{GMm}{r^2}

This is identical to the field of a point mass MM at the centre. \square

Proof (internal, r<Rr \lt R). The same integral with r<Rr \lt R evaluates to zero. Physically, for Every mass element pulling the test mass in one direction, there is a compensating element on the Opposite side. The nearer element pulls more strongly (shorter distance) but is subtended by a smaller Solid angle, and these two effects cancel exactly. \square

### Field Strength at Altitude

At height hh above a planet of radius RR and surface field g0g_0:

g=g0(RR+h)2\boxed{g = g_0\left(\frac{R}{R + h}\right)^2}

For hRh \ll RThe binomial approximation gives gg0(12h/R)g \approx g_0(1 - 2h/R).

Definition. The gravitational potential VV at a point is the work done per unit mass in bringing a Small test mass from infinity to that point:

V=GMr\boxed{V = -\frac{GM}{r}}

SI units: J kg1^{-1}.

V=Wm=1mrGMmr2dr=GMrdrr2=GM[1r]r=GMrV = \frac{W}{m} = \frac{1}{m}\int_{\infty}^{r} \frac{GMm}{r'^2}\,dr' = GM\int_{\infty}^{r}\frac{dr'}{r'^2} = GM\left[-\frac{1}{r'}\right]_{\infty}^{r} = -\frac{GM}{r}

\square

Why negative? We define V=0V = 0 at infinity. As the test mass moves inward, gravity does positive Work, so the potential decreases. The potential at any finite rr is therefore negative, becoming more Negative as rr decreases. An external agent must supply energy GMm/rGMm/r to move the mass from rr back To infinity.

g=dVdr\boxed{g = -\frac{dV}{dr}}

Proof. V=GM/r=GMr1V = -GM/r = -GM r^{-1}. Then:

dVdr=GMr2=GMr2\frac{dV}{dr} = GM\,r^{-2} = \frac{GM}{r^2}

The field (directed inward, i.e. In the negative radial direction) is:

gr=GMr2=dVdrg_r = -\frac{GM}{r^2} = -\frac{dV}{dr}

\square

The minus sign confirms that the field points in the direction of decreasing potential.

For two masses MM and mm separated by rr:

Ep=GMmr\boxed{E_p = -\frac{GMm}{r}}

Connection to Ep=mghE_p = mgh. For height hREh \ll R_EThe Taylor expansion gives:

ΔEp=GMmRE+h+GMmRE=GMmhRE(RE+h)GMmhRE2=mgh\Delta E_p = -\frac{GMm}{R_E + h} + \frac{GMm}{R_E} = \frac{GMmh}{R_E(R_E + h)} \approx \frac{GMmh}{R_E^2} = mgh

Since g=GM/RE2g = GM/R_E^2. The linear formula mghmgh is the first-order approximation of the full Gravitational potential energy.

Definition. The escape velocity vev_e is the minimum launch speed for an object to reach infinity With zero residual speed from the surface of a body of mass MM and radius RR.

At launch: E_k = \frac{1}{2}mv_e^2$$E_p = -GMm/R. At infinity: E_k = 0$$E_p = 0. By energy Conservation:

12mve2GMmR=0\frac{1}{2}mv_e^2 - \frac{GMm}{R} = 0

ve=2GMR\boxed{v_e = \sqrt{\frac{2GM}{R}}}

\square

The circular orbital speed at radius rr is vorb=GM/rv_{\mathrm{orb}} = \sqrt{GM/r}. Therefore:

ve=2vorbv_e = \sqrt{2}\,v_{\mathrm{orb}}

Escape requires exactly twice the kinetic energy of a circular orbit: 12mve2=2×12mvorb2\frac{1}{2}mv_e^2 = 2 \times \frac{1}{2}mv_{\mathrm{orb}}^2. A spacecraft in circular orbit needs a speed increase of (21)×100%41.4%(\sqrt{2} - 1) \times 100\% \approx 41.4\% to escape.

Worked Example: Escape from MarsCalculate the escape velocity from Mars ($M = 6.42 \times 10^{23}$ kg, $R = 3.39 \times 10^6$ m).

Answer. ve=2×6.67×1011×6.42×10233.39×106=8.56×10133.39×106=2.53×107=5020v_e = \sqrt{\frac{2 \times 6.67 \times 10^{-11} \times 6.42 \times 10^{23}}{3.39 \times 10^6}} = \sqrt{\frac{8.56 \times 10^{13}}{3.39 \times 10^6}} = \sqrt{2.53 \times 10^7} = 5020 m s1^{-1} =5.02= 5.02 km s1^{-1}.

## 5. Orbital Mechanics

For a satellite of mass mm in a circular orbit of radius rr around mass MM:

QuantityExpression
Orbital speedv=GM/rv = \sqrt{GM/r}
Orbital periodT=2πr3/(GM)T = 2\pi\sqrt{r^3/(GM)}
Centripetal accelerationa=GM/r2=ga = GM/r^2 = g
Kinetic energyEk=GMm/(2r)E_k = GMm/(2r)
Potential energyEp=GMm/rE_p = -GMm/r
Total energyEtotal=GMm/(2r)E_{\mathrm{total}} = -GMm/(2r)

Proof of energy relations. From GMmr2=mv2r\frac{GMm}{r^2} = \frac{mv^2}{r}: v2=GM/rv^2 = GM/r.

Ek=12mv2=GMm2rE_k = \frac{1}{2}mv^2 = \frac{GMm}{2r}. Ep=GMm/rE_p = -GMm/r. Total: Ek+Ep=GMm2rE_k + E_p = -\frac{GMm}{2r}.

Note: Ek=12EpE_k = \frac{1}{2}|E_p| and Etotal=EkE_{\mathrm{total}} = E_k. This is the virial theorem for bound Gravitational systems: 2Ek+Ep=02E_k + E_p = 0. \square

Key insight. The total energy is negative — the satellite is gravitationally bound. To move to a Higher orbit, energy must be added (the orbit becomes less negative). The kinetic energy decreases With increasing rrBut the total energy increases (potential energy increase dominates).

For any Keplerian orbit (circular or elliptical) with semi-major axis aa:

v2=GM(2r1a)\boxed{v^2 = GM\left(\frac{2}{r} - \frac{1}{a}\right)}

Setting a=ra = r recovers the circular orbit result. For a parabolic escape trajectory (aa \to \infty): v2=2GM/rv^2 = 2GM/rGiving the escape speed.

v=GMrv = \sqrt{\frac{GM}{r}}

The satellite’s mass mm cancels. This is the same reason that all objects fall at the same rate in a Gravitational field (equivalence principle).

Worked Example: Satellite Orbit ChangeA 500 kg satellite is in a circular orbit of radius $7.0 \times 10^6$ m. Calculate the energy Required to move it to a circular orbit of radius $1.4 \times 10^7$ m.

Answer. GM=3.98×1014GM = 3.98 \times 10^{14} N m2^2 kg1^{-1}.

E1=GMm2r1=3.98×1014×5002×7.0×106=1.42×1010E_1 = -\frac{GMm}{2r_1} = -\frac{3.98 \times 10^{14} \times 500}{2 \times 7.0 \times 10^6} = -1.42 \times 10^{10} J.

E2=3.98×1014×5002×1.4×107=7.11×109E_2 = -\frac{3.98 \times 10^{14} \times 500}{2 \times 1.4 \times 10^7} = -7.11 \times 10^9 J.

ΔE=E2E1=7.11×109(1.42×1010)=7.1×109\Delta E = E_2 - E_1 = -7.11 \times 10^9 - (-1.42 \times 10^{10}) = 7.1 \times 10^9 J =7.1= 7.1 GJ.

Every planet moves in an elliptical orbit with the Sun at one focus.

Proof sketch. Starting from F=GMmr2r^\mathbf{F} = -\frac{GMm}{r^2}\hat{\mathbf{r}} (central force), the orbit Equation in polar coordinates is:

r=a(1e2)1+ecosθr = \frac{a(1 - e^2)}{1 + e\cos\theta}

Where aa is the semi-major axis and ee is the eccentricity. For E<0E \lt 0 (bound orbit), e<1e \lt 1 And the orbit is an ellipse with the central mass at one focus. \square

A line joining a planet to the Sun sweeps out equal areas in equal times.

Proof. The gravitational force is central (Fr\mathbf{F} \parallel \mathbf{r}), so torque τ=r×F=0\boldsymbol{\tau} = \mathbf{r} \times \mathbf{F} = 0. Angular momentum L=mrvL = mrv_\perp is constant.

Area swept in time dtdt: dA=12rvdt=L2mdtdA = \frac{1}{2}r \cdot v_\perp\,dt = \frac{L}{2m}\,dt.

dAdt=L2m=const\frac{dA}{dt} = \frac{L}{2m} = \mathrm{const}

\square

The planet moves fastest at perihelion (closest approach) and slowest at aphelion (farthest point), Consistent with conservation of angular momentum: small rr requires large vv_\perp.

T2=4π2GMa3\boxed{T^2 = \frac{4\pi^2}{GM}\,a^3}

Proof for circular orbits. Equating gravitational and centripetal force:

GMmr2=mv2r    v=GMr\frac{GMm}{r^2} = \frac{mv^2}{r} \implies v = \sqrt{\frac{GM}{r}}

Since T=2πr/vT = 2\pi r/v:

T=2πrGM/r=2πr3GMT = \frac{2\pi r}{\sqrt{GM/r}} = 2\pi\sqrt{\frac{r^3}{GM}}

T2=4π2r3GMT^2 = \frac{4\pi^2 r^3}{GM}

For elliptical orbits, replace rr with the semi-major axis aa. \square

The proportionality constant 4π2/(GM)4\pi^2/(GM) depends only on the central body, not on the orbiting Object. This is how Kepler determined the relative distances of the planets from the Sun using only Their observed periods.

Definition. A geostationary orbit is a circular, prograde, equatorial orbit with period equal to One sidereal day (T=86164T = 86164 s), causing the satellite to remain fixed above a point on the equator.

r3=GMT24π2=3.98×1014×(86164)24π2=7.54×1022r^3 = \frac{GMT^2}{4\pi^2} = \frac{3.98 \times 10^{14} \times (86164)^2}{4\pi^2} = 7.54 \times 10^{22}

r=4.22×107 m=42200 km\boxed{r = 4.22 \times 10^7\ \mathrm{m} = 42\,200\ \mathrm{km}}

Altitude above Earth’s surface: h=422006370=35830h = 42200 - 6370 = 35\,830 km.

v=GMr=2πrT=3070 ms1v = \sqrt{\frac{GM}{r}} = \frac{2\pi r}{T} = 3070\ \mathrm{m\,s}^{-1}

  1. Correct radius: r42200r \approx 42\,200 km (from Kepler’s third law with T=86164T = 86164 s).
  2. Equatorial plane: Any inclination causes the satellite to trace a figure-eight (anisotropic pattern) as seen from the ground.
  3. Prograde rotation: The satellite must orbit west to east, matching Earth’s rotation.

From Newton’s apple to quantum particles, physics explains how the world works through measurable quantities and testable laws. Forces cause acceleration, energy transforms between forms but is never lost, and waves carry information across vast distances. These concepts form the foundation for engineering, astronomy, and our understanding of reality itself.