Login

Concepts

Read through concept notes with worked explanations and diagrams, chapter by chapter.

Concept 06 / 08Gravitation

Kepler's Three Laws of Planetary Motion and Their Dynamical Origin

Planets move on ellipses with the Sun at a focus, sweep equal areas in equal times because angular momentum is conserved, and obey T squared proportional to a cubed, which follows directly from the inverse-square law.

Figure 1
Kepler extracted three empirical laws from Tycho Brahe's naked-eye observations of Mars between 1609 and 1619, decades before Newton. Newton then showed that all three are consequences of the inverse-square law, which is why they are the sharpest evidence for it. Law 1, the law of orbits. Every planet moves in an ellipse with the Sun at one focus. The other focus is empty. In polar coordinates centred on the Sun, the orbit is r=l1+ecos⁡θr = \frac{l}{1+e\cos\theta}, where ee is the eccentricity and l=a(1−e2)=b2al = a(1-e^2) = \frac{b^2}{a} is the semi-latus rectum. The nearest point (perihelion for the Sun, perigee for the Earth) is at rmin=a(1−e)r_{min} = a(1-e) and the farthest (aphelion, apogee) at rmax=a(1+e)r_{max} = a(1+e). Adding these, rmin+rmax=2ar_{min}+r_{max} = 2a: the semi-major axis is the arithmetic mean of the extreme distances. A circle is the special case e=0e=0. Most planetary orbits are nearly circular: the Earth has e=0.017e = 0.017, so its distance from the Sun varies by only about 3 percent over the year, which is why the seasons cannot be caused by it. Comets have ee close to 1. Law 2, the law of areas. The line joining a planet to the Sun sweeps out equal areas in equal intervals of time. In the diagram above the two shaded sectors have equal area, so the planet covers the fat short arc near the Sun in the same time as the long thin arc far away; hence it moves fastest at perihelion. This law is simply angular momentum conservation in disguise, and it holds for any central force whatsoever, not just an inverse-square one. The proof: the area swept in time dtdt by the position vector is the area of the thin triangle, dA=12∣r⃗×dr⃗∣=12∣r⃗×v⃗∣dtdA = \frac{1}{2}|\vec{r}\times d\vec{r}| = \frac{1}{2}|\vec{r}\times\vec{v}|dt. Therefore dAdt=12∣r⃗×v⃗∣=∣L⃗∣2m\frac{dA}{dt} = \frac{1}{2}|\vec{r}\times\vec{v}| = \frac{|\vec{L}|}{2m}. Since the gravitational force is central, its torque about the Sun, τ⃗=r⃗×F⃗\vec{\tau} = \vec{r}\times\vec{F}, is zero because F⃗\vec{F} is parallel to r⃗\vec{r}. So L⃗\vec{L} is constant, and the areal velocity dAdt\frac{dA}{dt} is constant. The constancy of the direction of L⃗\vec{L} also tells us the orbit is planar. An immediate corollary: at perihelion and aphelion the velocity is perpendicular to the radius, so L=mvprp=mvaraL = mv_pr_p = mv_ar_a, giving vprp=varav_pr_p = v_ar_a. Law 3, the law of periods. The square of the period is proportional to the cube of the semi-major axis: T2∝a3T^2 \propto a^3, with the same constant of proportionality for every body orbiting the same central mass. Derivation for a circular orbit (this is the derivation JEE expects). For a planet of mass mm in a circular orbit of radius rr about a mass MM, gravity supplies the centripetal force: GMmr2=mω2r=m4π2T2r\frac{GMm}{r^2} = m\omega^2r = m\frac{4\pi^2}{T^2}r. Cancel mm and rearrange: T2=4π2GMr3T^2 = \frac{4\pi^2}{GM}r^3. The full result for an ellipse replaces rr by the semi-major axis aa: T2=4π2a3GMT^2 = \frac{4\pi^2 a^3}{GM}. The constant 4π2GM\frac{4\pi^2}{GM} depends only on the central mass, which is why every planet, asteroid and comet in the Solar System lies on one straight line in a plot of log⁡T\log T against log⁡a\log a, and why the Moon and every artificial satellite lie on a different line set by the Earth's mass. Strictly the two-body result is T2=4π2a3G(M+m)T^2 = \frac{4\pi^2a^3}{G(M+m)}, where aa is now the semi-major axis of the relative orbit. For planets m≪Mm \ll M so the correction is tiny (for Jupiter it is one part in a thousand), but for binary stars, where the masses are comparable, it is the whole point: measuring TT and aa gives the total mass of the system. Two energy results tie the laws together and are worth memorising, because they turn many hard problems into one-liners. The total energy of an elliptical orbit depends only on the semi-major axis: E=−GMm2aE = -\frac{GMm}{2a}. Circular orbits are the case a=ra = r. The vis-viva equation, obtained by combining that with E=12mv2−GMmrE = \frac{1}{2}mv^2-\frac{GMm}{r}: v2=GM(2r−1a)v^2 = GM\left(\frac{2}{r}-\frac{1}{a}\right). This single formula gives the speed anywhere on any orbit. At r=ar=a it gives v=GMav = \sqrt{\frac{GM}{a}}, the circular speed; for a→∞a \to \infty (parabolic escape) it gives v=2GMrv = \sqrt{\frac{2GM}{r}}, the escape speed. Gotchas. The Sun sits at a focus, not at the centre. The quantity in Law 3 is the semi-major axis, not the average of the distances over time (the time-average of rr is actually a(1+e22)a(1+\frac{e^2}{2})). Law 2 is about areal speed, which is constant; the linear speed and angular speed are certainly not constant. Kepler's laws apply to any inverse-square two-body system, including satellites around the Earth and electrons in the Bohr model. And the laws fail slightly when other planets perturb the orbit or when general relativity matters, which is how the precession of Mercury's perihelion was discovered. Key Takeaways: - Law 1: ellipse, Sun at a focus; rmin=a(1−e)r_{min}=a(1-e), rmax=a(1+e)r_{max}=a(1+e), rmin+rmax=2ar_{min}+r_{max}=2a. - Law 2: dAdt=L2m\frac{dA}{dt} = \frac{L}{2m} is constant, because a central force exerts no torque; valid for any central force. Corollary vprp=varav_pr_p = v_ar_a. - Law 3: T2=4π2a3GMT^2 = \frac{4\pi^2a^3}{GM}; two-body form T2=4π2a3G(M+m)T^2 = \frac{4\pi^2a^3}{G(M+m)}. - Energy: E=−GMm2aE = -\frac{GMm}{2a} (depends only on aa). Vis-viva: v2=GM(2r−1a)v^2 = GM\left(\frac{2}{r}-\frac{1}{a}\right). - Angular momentum: L=mGMa(1−e2)L = m\sqrt{GMa(1-e^2)}, depends on both aa and ee; energy depends only on aa. - Common mistakes: putting the Sun at the centre; using ravgr_{avg} instead of aa; assuming angular velocity is constant because areal velocity is; forgetting that a body on a given orbit has the same period regardless of its own mass. Brain-Teaser 1 (a comet's aphelion): A comet passes perihelion at a distance rpr_p from the Sun with speed vp=1.6 GMrpv_p = \sqrt{\frac{1.6\,GM}{r_p}}. Find its aphelion distance, the ratio of its speeds at the two ends, and its period in terms of rpr_p. Solution: First get aa from the energy, which needs no angular momentum bookkeeping. Em=12vp2−GMrp=0.8GMrp−GMrp=−0.2GMrp\frac{E}{m} = \frac{1}{2}v_p^2-\frac{GM}{r_p} = \frac{0.8GM}{r_p}-\frac{GM}{r_p} = -\frac{0.2GM}{r_p}. Since Em=−GM2a\frac{E}{m} = -\frac{GM}{2a}, we get GM2a=0.2GMrp\frac{GM}{2a} = \frac{0.2GM}{r_p}, so a=2.5rpa = 2.5r_p. Then ra=2a−rp=5rp−rp=4rpr_a = 2a - r_p = 5r_p - r_p = 4r_p. Speeds: angular momentum conservation gives vpva=rarp=4\frac{v_p}{v_a} = \frac{r_a}{r_p} = 4. Verify independently with vis-viva at aphelion: va2=GM(24rp−12.5rp)=GMrp(0.5−0.4)=0.1GMrpv_a^2 = GM\left(\frac{2}{4r_p}-\frac{1}{2.5r_p}\right) = \frac{GM}{r_p}(0.5-0.4) = \frac{0.1GM}{r_p}. Since vp2=1.6GMrpv_p^2 = \frac{1.6GM}{r_p}, the ratio vp2va2=16\frac{v_p^2}{v_a^2} = 16, so vpva=4\frac{v_p}{v_a} = 4. The two methods agree. Period: T=2πa3GM=2π15.625 rp3GM=2.52.5⋅2πrp3GMT = 2\pi\sqrt{\frac{a^3}{GM}} = 2\pi\sqrt{\frac{15.625\,r_p^3}{GM}} = 2.5\sqrt{2.5}\cdot 2\pi\sqrt{\frac{r_p^3}{GM}}, i.e. about 3.953.95 times the period of a circular orbit of radius rpr_p. Eccentricity, for completeness: rp=a(1−e)r_p = a(1-e) gives e=1−rp2.5rp=0.6e = 1-\frac{r_p}{2.5r_p} = 0.6. Brain-Teaser 2 (how long to fall into the Sun): If the Earth were suddenly stopped in its orbit, how long would it take to fall into the Sun? Answer without doing any integral. Solution: A body dropped radially is not a special case at all, it is a degenerate ellipse of eccentricity e→1e \to 1, squashed flat. Its aphelion is the starting distance rr, and its perihelion is essentially zero (the centre of the Sun). So 2a′=rmax+rmin=r+0⇒a′=r22a' = r_{max}+r_{min} = r + 0 \Rightarrow a' = \frac{r}{2}. Kepler's third law compares this degenerate orbit with the Earth's actual circular orbit of radius rr and period T=1T = 1 year: T′T=(a′a)3/2=(12)3/2=122\frac{T'}{T} = \left(\frac{a'}{a}\right)^{3/2} = \left(\frac{1}{2}\right)^{3/2} = \frac{1}{2\sqrt{2}}. The fall to the centre is half of that degenerate orbit (down and back would be the full "orbit"), so the falling time is t=T′2=T42=365.255.657=64.6t = \frac{T'}{2} = \frac{T}{4\sqrt{2}} = \frac{365.25}{5.657} = 64.6 days. A three-line answer to a problem that looks like it needs a nasty integral. The same trick answers "how long does a satellite take to fall to the Earth's centre if stopped", and it is a perennial favourite of JEE Advanced. Brain-Teaser 3 (the speed at the end of the minor axis): A planet moves on an ellipse of semi-major axis aa and eccentricity ee. Show that its speed at either end of the minor axis equals the geometric mean of its perihelion and aphelion speeds, and that it equals the speed of a circular orbit of radius aa. Solution: From vis-viva, at perihelion r=a(1−e)r = a(1-e): vp2=GM(2a(1−e)−1a)=GMa⋅2−(1−e)1−e=GMa⋅1+e1−ev_p^2 = GM\left(\frac{2}{a(1-e)}-\frac{1}{a}\right) = \frac{GM}{a}\cdot\frac{2-(1-e)}{1-e} = \frac{GM}{a}\cdot\frac{1+e}{1-e}. Similarly at aphelion r=a(1+e)r = a(1+e): va2=GMa⋅1−e1+ev_a^2 = \frac{GM}{a}\cdot\frac{1-e}{1+e}. Multiplying, vp2va2=(GMa)2v_p^2v_a^2 = \left(\frac{GM}{a}\right)^2, and since both speeds are positive, vpva=GMav_pv_a = \frac{GM}{a}, so the geometric mean is vpva=GMa\sqrt{v_pv_a} = \sqrt{\frac{GM}{a}}. Now the end of the minor axis. The distance from a focus to that point is exactly aa, since b2+c2=a2(1−e2)+a2e2=a\sqrt{b^2+c^2} = \sqrt{a^2(1-e^2)+a^2e^2} = a. Vis-viva there gives v2=GM(2a−1a)=GMav^2 = GM\left(\frac{2}{a}-\frac{1}{a}\right) = \frac{GM}{a}. So v=GMa=vpvav = \sqrt{\frac{GM}{a}} = \sqrt{v_pv_a}, which is both the geometric mean of the extreme speeds and the speed on a circular orbit of radius aa. The conceptual payoff: the total energy of the elliptical orbit depends only on aa, so at the one place where rr happens to equal aa the kinetic energy must match that of the circular orbit of radius aa. Energy knows only about aa; the eccentricity is stored entirely in the angular momentum, L=mGMa(1−e2)L = m\sqrt{GMa(1-e^2)}.