courseshub.world · Lives in Mathematics and Physics

Charles Messier

Charles Messier, Astronomer of the French Navy

A course in two strands: the life of the clerk from Lorraine whom Louis XV called the ferret of comets, and the astronomy behind his work, from the orbits of Halley’s and Lexell’s comets to the 110 star clusters, nebulae and galaxies that still carry his initial.

Born 26 June 1730, Badonviller Died 12 April 1817, Paris Catalogue 103 objects, 1781 Level advanced undergraduate
Schematic. A nebula keeps its place among the stars; a comet moves from night to night. Messier listed the fixed glows, beginning with M1 near ζ Tauri, so that comet hunters would not be fooled by them.
Lesson 00

How to use this course

Charles Messier is remembered for a list he made to help himself. He wanted comets, and the sky kept offering him faint glows that looked like comets but did not move. He recorded them so that neither he nor anyone else would waste a night on them again. The list grew to 103 objects in 1781, and seven more were added in the twentieth century from his and Pierre Méchain’s notes. It contains the Crab Nebula, the Andromeda galaxy, the Orion Nebula, the Pleiades and the galaxy M87, whose central black hole was imaged in 2019.

Messier himself would have pointed to his comets. He observed more than forty, discovered about a dozen and a half, and recovered Halley’s comet on its first predicted return. He was not a theorist. His skill was patient, accurate observation, and the mathematics of this course is the mathematics that turned observations like his into orbits and physics.

Prerequisites for Part II
  • Trigonometry, logarithms and vectors in three dimensions.
  • Newtonian gravitation and Kepler’s laws.
  • Basic calculus, including small-angle approximations.
  • For Lessons 09 and 10, the magnitude scale and the idea of a parsec.

Each technical lesson opens with an italic paragraph for the general reader, then gives precise statements, derivations and exercises with folded solutions. The companion modules on Caroline Herschel and William Herschel treat sweeping, sidereal time, precession and the surface brightness of nebulae in detail; this course refers to them rather than repeating them.

Reference

Timeline

1730Born 26 June in Badonviller, in the principality of Salm in Lorraine, tenth of twelve children of Nicolas Messier, a court usher, and Françoise Grandblaise.
1741Death of his father.
1744Sees the great comet with six tails (C/1743 X1).
174825 July: annular eclipse of the Sun visible from Badonviller.
1751Goes to Paris and enters the service of Joseph-Nicolas Delisle as draughtsman and copyist.
17536 May: first recorded observation, a transit of Mercury, from the Hôtel de Cluny.
1754Employed as a clerk of the Navy’s map depot.
1757–58Searches for Halley’s comet using Delisle’s chart. In August and September 1758 he finds and measures the nebula in Taurus later numbered M1.
175921 January: recovers Halley’s comet, about four weeks after Johann Georg Palitzsch. Delisle announces the observations on 1 April.
1760sA run of comet discoveries; Louis XV calls him the ferret of comets.
1764Fellow of the Royal Society of London. Finds the Dumbbell Nebula, M27.
1769Foreign member of the Royal Swedish Academy of Sciences. Discovers the great comet of 1769 on 8 August.
177014 June: discovers the comet later named after Anders Johan Lexell, which passes 0.015 AU from Earth on 1 July. Elected to the Paris Academy of Sciences, 30 June. Marries Marie-Françoise de Vermauchampt.
177116 February: presents his catalogue of 45 nebulae and star clusters to the Academy. Becomes Astronomer of the Navy.
c. 1772Death of his wife and, days later, of their newborn son.
1774The 45-object catalogue appears in the Academy’s Mémoires for 1771. Méchain begins to work with him.
1780Second list, 68 objects, in the Connaissance des Temps.
1781Final list of 103 objects, printed in the Connaissance des Temps for 1784. Around this time he is badly injured in a fall into an ice cellar.
1790sThe Revolution ends his pension and salary.
1796Appointed to the Bureau des Longitudes; member of the new Institut.
1806Receives the cross of the Legion of Honour from Napoleon.
1808Publishes a memoir on the great comet of 1769, the year of Napoleon’s birth.
1815Partly paralysed by a stroke.
1817Dies 12 April in Paris, aged 86. Buried in the Père Lachaise cemetery.
Part I · Formation
Lesson 01

Badonviller, 1730–1751

Charles Messier was born on 26 June 1730 in Badonviller, a small town on the western slopes of the Vosges. It then belonged to the principality of Salm, a pocket of the Holy Roman Empire inside Lorraine. He was the tenth of twelve children of Nicolas Messier, an usher of the princely court, and Françoise Grandblaise. Six of the children died young, and his father died in 1741, when Charles was eleven. He grew up with his elder brother Hyacinthe, who held an administrative post in the service of the princes of Salm.

Two events turned the boy to the sky. In 1744 a spectacular comet appeared, the great comet of 1743–44, whose fan of six tails rose before dawn in early March. On 25 July 1748 an annular eclipse of the Sun was visible from Badonviller. Messier had no formal training in mathematics or astronomy. What he had was good eyesight, patience and a neat hand, which proved to be his introduction to Paris.

In 1751 the principality of Salm was partitioned and the family’s position changed. Hyacinthe moved to Senones, and Charles, now 21, went to Paris to look for work.

Lesson 02

Delisle’s clerk and Halley’s comet, 1751–1759

In Paris Messier found a place with Joseph-Nicolas Delisle, astronomer of the Navy and professor at the Collège de France. Delisle had spent more than twenty years in St Petersburg and had brought back a large collection of maps and observations. He needed a copyist, and Messier’s first work was to copy maps, including a map of Beijing and of the Great Wall of China. From 1754 Messier was paid as a clerk of the Navy’s depot of charts.

Delisle’s observatory was a small tower on the roof of the Hôtel de Cluny, the medieval town house of the abbots of Cluny near the Sorbonne. Messier learned to observe there. His first recorded observation is the transit of Mercury of 6 May 1753. He would live and observe at the Hôtel de Cluny for most of his life.

Waiting for the comet

Edmond Halley had predicted that the comet of 1682 would return around 1758. Alexis-Claude Clairaut, with Joseph-Jérôme Lalande and Nicole-Reine Lepaute, refined the prediction by computing the perturbations of Jupiter and Saturn. Delisle drew up a chart of the path the comet should follow across the sky, and from 1757 Messier searched along it. The chart was in error, and for months he searched in the wrong places. In August 1758, while following a different comet, he came upon a faint nebula near the star ζ Tauri that looked like a comet but did not move. He measured its position on 12 September. It would become the first entry in his catalogue, and is now known as the Crab Nebula.

On 21 January 1759 Messier found Halley’s comet. Johann Georg Palitzsch, a farmer and amateur astronomer near Dresden, had seen it on Christmas night 1758, so Messier’s sighting was not the first, but it was independent and well measured. Delisle did not allow him to announce it, and the observations were made public only on 1 April 1759, when the comet had passed perihelion and was moving away. The delay, which was not Messier’s choice, led some to doubt the observations. Lesson 06 computes the orbit he was searching for.

Lesson 03

The ferret of comets, 1759–1771

After 1759 Messier devoted his nights to comets. Hunting them meant sweeping the evening and morning skies near the Sun, where comets brighten, again and again. In the 1760s he found a series of new comets, and his record became famous at court. Louis XV called him the furet des comètes, the ferret of comets. Today 13 comets are credited to him as discoverer, among them the comets of 1763 and 1764, and he observed many more. By Lalande’s count he observed 41 in his career.

Recognition came from abroad first. He was elected a Fellow of the Royal Society of London in 1764 and a foreign member of the Royal Swedish Academy of Sciences in 1769. He was later elected to academies in Germany, Russia and elsewhere. The Paris Academy of Sciences admitted him, after some difficulty, on 30 June 1770.

A catalogue of false alarms

The same sweeps that found comets turned up nebulae and clusters. In 1764 he found the Dumbbell Nebula in Vulpecula (M27) and observed the nebula in Andromeda (M31), described in the tenth century by the Persian astronomer Abd al-Rahman al-Sufi. He recorded each object with its position and a short description. On 16 February 1771 he presented a list of 45 to the Academy. It appeared in print in 1774, in the Mémoires for 1771. Of the 45, he had discovered 17 or 18 himself; the rest were known from earlier observers, including Edmond Halley, Jean-Philippe Loys de Chéseaux and Nicolas-Louis de Lacaille, whose positions he checked.

In 1770, at the age of forty, he married Marie-Françoise de Vermauchampt. Their son was born in 1772. His wife died, and the child died eleven days after his birth. A well-known anecdote says that Messier missed a comet, found by a rival, while he nursed his wife. In 1771 he succeeded Delisle as Astronomer of the Navy.

Part II · The astronomy
Lesson 04

Small telescopes and faint light

For the general reader. Messier worked with instruments that a modern amateur would find poor. His telescopes had apertures of a few inches, and many had single lenses whose images were fringed with colour, or metal mirrors that reflected barely half the light. Yet he found comets and nebulae at the limit of what such instruments could show. This lesson explains what an aperture buys, why eighteenth-century telescopes were so long, and how faint a star such an instrument could reveal.

Goals
  • Compute the light grasp and limiting magnitude of a telescope.
  • Understand the diffraction limit.
  • Derive the condition for an achromatic doublet and see why long focal lengths were used before it.

Messier’s instruments

He used a succession of instruments at the Hôtel de Cluny: refractors with focal lengths up to several metres, and Gregorian reflectors with speculum-metal mirrors up to about 20 cm across. His usual comet-seeking telescopes had apertures of about 10 cm. The reflectors, though larger, delivered less light than their size suggests, because speculum metal reflects only about 60 per cent of the light falling on it, and a Gregorian uses two mirrors.

Proposition (light grasp and limiting magnitude)

A telescope of aperture \(D\) with transmission \(\tau\) collects \(\tau(D/p)^2\) times as much light from a star as an eye of pupil \(p\). If the naked-eye limit is magnitude \(m_0\), the telescopic limit is approximately

\[ m_{\lim}\approx m_0+5\log_{10}\frac{D}{p}+2.5\log_{10}\tau. \]

The magnitude scale is logarithmic: a factor \(F\) in flux is \(2.5\log_{10}F\) magnitudes. The flux gain from the area is \((D/p)^2\), which gives the \(5\log_{10}(D/p)\) term. A star stays a point at any magnification, so, unlike a nebula (Lesson 10 of the Caroline Herschel module), its brightness grows with aperture.

Diffraction limit

A circular aperture of diameter \(D\) images a point source as an Airy pattern whose first dark ring has angular radius \(\theta=1.22\,\lambda/D\). Two equal stars separated by \(\theta\) are just resolved (Rayleigh’s criterion).

Colour and the long telescope

A single lens of refractive index \(n\) has power \(P=1/f=(n-1)K\), where \(K\) depends only on the curvatures of its surfaces. Because \(n\) depends on wavelength, so does \(f\). Write \(V=(n_d-1)/(n_F-n_C)\) for the Abbe number of the glass, which measures how weakly it disperses. Then the spread in power between the blue F line and red C line is \(\Delta P=P/V\), so the focal lengths differ by about \(\Delta f\approx f/V\).

Proposition (chromatic blur and the achromat)

(a) For a single lens of aperture \(D\) and focal length \(f\), the angular diameter of the colour blur at best focus is roughly

\[ \beta\approx\frac{D}{2Vf}. \]

(b) Two thin lenses in contact with powers \(P_1,P_2\) and Abbe numbers \(V_1,V_2\) have the same focal length in the F and C lines if and only if

\[ \frac{P_1}{V_1}+\frac{P_2}{V_2}=0. \]

Sketch. (a) Red and blue rays focus \(\Delta f\approx f/V\) apart. Midway between the foci, a cone of rays of diameter \(D\) converging over a length \(f\) has diameter about \(D\,\Delta f/(2f)\), which subtends \(D/(2Vf)\) at the lens. (b) Powers of thin lenses in contact add, so the change in total power is \(P_1/V_1+P_2/V_2\), which must vanish. \(\square\)

Part (a) explains the absurdly long refractors of the seventeenth and eighteenth centuries: the blur in angle falls as \(1/f\). Part (b) is the principle of the achromatic doublet, a converging crown lens paired with a weaker diverging flint lens, patented by John Dollond in 1758. Achromats spread slowly, and Messier used both kinds of refractor.

Exercise 4.1Take a naked-eye limit of magnitude 6, a pupil of 7 mm, and a 100 mm refractor with \(\tau=0.8\). Estimate the limiting magnitude and the Rayleigh resolution at \(\lambda=550\) nm.

Show solution
\(5\log_{10}(100/7)=5.77\) and \(2.5\log_{10}0.8=-0.24\), so \(m_{\lim}\approx11.5\). Resolution: \(\theta=1.22\times550\times10^{-9}/0.1=6.7\times10^{-6}\) rad \(\approx1.4''\). Such a telescope shows about a hundred times as many stars as the naked eye, and could not resolve most globular clusters into stars, which is why Messier called many of them nebulae without stars.

Exercise 4.2A single crown-glass lens has \(V=59\) and \(D=100\) mm. Estimate the colour blur for \(f=1\) m and for \(f=7\) m. What focal length would bring the blur down to the diffraction limit of Exercise 4.1?

Show solution
For \(f=1\) m, \(\beta\approx0.1/(2\times59\times1)=8.5\times10^{-4}\) rad \(\approx175''\), about three arcminutes. For \(f=7\) m, \(\beta\approx25''\). To reach \(6.7\times10^{-6}\) rad one would need \(f\approx0.1/(2\times59\times6.7\times10^{-6})\approx130\) m. This is why single-lens telescopes were made long, and why even long ones never reached their diffraction limit.

Exercise 4.3A doublet of focal length 1 m uses crown with \(V_1=59\) and flint with \(V_2=36\). Find the focal lengths of the two lenses.

Show solution
\(P_1+P_2=1\) D and \(P_1/59+P_2/36=0\), so \(P_2=-\tfrac{36}{59}P_1\). Then \(P_1(1-36/59)=1\), giving \(P_1=2.565\) D and \(P_2=-1.565\) D: a crown lens of \(f_1=0.39\) m and a flint lens of \(f_2=-0.64\) m.
Lesson 05

Measuring a comet’s place

For the general reader. A comet’s orbit is computed from its positions on the sky over several nights. Messier did not measure those positions from scratch. He compared the comet with a nearby star whose position was already catalogued, measuring the small differences in right ascension and declination. He then had to correct for the bending of starlight by the atmosphere, which lifts objects near the horizon, where comets close to the Sun are usually found.

Goals
  • Convert transit-time and declination differences into positions.
  • Derive the refraction formula for a plane-parallel atmosphere.
  • Estimate the size of refraction errors for low comets.

Differential positions

With a fixed telescope and cross-wires, an observer times the passage of the comet and of a comparison star across a wire perpendicular to the diurnal motion, and measures their separation in declination with a micrometer. If a clock keeping mean solar time shows an interval \(\Delta t\) seconds between the two transits, the difference in right ascension is \(\Delta\alpha=1.0027379\,\Delta t\) seconds of sidereal time, and the angular separation along the parallel of declination is

\[ s=15''\,\Delta\alpha\,\cos\delta, \qquad\text{equivalently}\qquad s=15.04''\,\Delta t\,\cos\delta. \]

The advantage of differential measurement is that most errors of the instrument and the clock cancel, because they affect star and comet alike. Its weakness is that the result is only as good as the star’s catalogue position, which for Messier often came from Lacaille or Flamsteed.

Proposition (refraction in a plane-parallel atmosphere)

If the atmosphere is a stack of horizontal layers and the refractive index at the ground is \(n_0\), a star at true zenith distance \(z_t\) is seen at apparent zenith distance \(z\) with

\[ n_0\sin z=\sin z_t,\qquad R=z_t-z\approx(n_0-1)\tan z. \]

At sea level \(n_0-1\approx2.8\times10^{-4}\), so \(R\approx58''\tan z\).

Proof. At each horizontal interface Snell’s law gives \(n_k\sin z_k=n_{k+1}\sin z_{k+1}\), so \(n\sin z\) is the same in every layer. Above the atmosphere \(n=1\) and the angle is \(z_t\); at the ground it is \(n_0\sin z\). Writing \(z_t=z+R\) with \(R\) small, \(\sin z+R\cos z\approx n_0\sin z\), so \(R\approx(n_0-1)\tan z\). \(\square\)

The result is independent of the structure of the atmosphere, which is why a single constant serves for most of the sky. It fails near the horizon, where the Earth’s curvature matters and \(\tan z\) diverges; the true refraction at the horizon is about \(35'\), and it varies with temperature and pressure. Refraction raises objects without changing their azimuth, so it shifts both right ascension and declination.

Exercise 5.1A comet transits 12.0 s of mean time after a comparison star at declination \(+22^\circ\), and lies \(4'10''\) north of it. Find the separation of the two objects in arcseconds.

Show solution
Along the parallel: \(s=15.04\times12.0\times\cos22^\circ\approx167''\). In declination: \(250''\). For a small field the separation is \(\sqrt{167^2+250^2}\approx301''\), about five arcminutes.

Exercise 5.2Estimate the refraction at apparent zenith distances \(45^\circ\) and \(70^\circ\). Then estimate the differential refraction between a comet and a star \(0.5^\circ\) apart in altitude at \(z=70^\circ\).

Show solution
\(R(45^\circ)\approx58''\) and \(R(70^\circ)\approx58''\tan70^\circ\approx160''\). Differentiating, \(dR/dz=58''\sec^2z\), which at \(70^\circ\) is \(58''\times8.55=496''\) per radian, or about \(8.7''\) per degree. Half a degree gives a differential error of about \(4''\): small compared with the absolute refraction, which is why differential measurement against a nearby star works well even low in the sky.
Lesson 06

Halley’s comet on its ellipse

For the general reader. Halley’s comet spends almost all of its 76-year orbit far beyond the planets, too faint to see. It is visible only for a few months near perihelion, its closest approach to the Sun. In 1758 observers knew roughly when it would come back but not exactly where to look. This lesson computes where the comet was on the night Messier found it, using the equation Kepler devised to locate a body on an ellipse at a given time.

Goals
  • State and derive Kepler’s equation.
  • Solve it numerically by Newton’s method.
  • Compute how long a long-period comet spends near the Sun.
Theorem (Kepler’s equation)

Let a body move on an ellipse with semi-major axis \(a\), eccentricity \(e\) and period \(P\) about the Sun at a focus, passing perihelion at time \(T\). Define the mean anomaly \(M=2\pi(t-T)/P\). The eccentric anomaly \(E\) satisfies

\[ M=E-e\sin E, \]

and the distance and true anomaly are

\[ r=a(1-e\cos E),\qquad \tan\frac{\nu}{2}=\sqrt{\frac{1+e}{1-e}}\,\tan\frac{E}{2}. \]

Proof. Put the ellipse’s centre at the origin with the Sun at \((ae,0)\), and parametrize by \(x=a\cos E\), \(y=b\sin E\) with \(b=a\sqrt{1-e^2}\). The area swept from perihelion, seen from the Sun, is the area of the elliptic sector from the centre, \(\tfrac12abE\), minus the triangle with vertices at the centre, the Sun and the body, \(\tfrac12\,ae\,b\sin E\). By Kepler’s second law this area is the fraction \((t-T)/P\) of the total \(\pi ab\):

\[ \frac{ab}{2}(E-e\sin E)=\pi ab\,\frac{t-T}{P}\quad\Longrightarrow\quad E-e\sin E=M. \]

The distance is \(r^2=(a\cos E-ae)^2+b^2\sin^2E=a^2(1-e\cos E)^2\). The half-angle formula follows from \(r\cos\nu=a(\cos E-e)\) and \(r\sin\nu=b\sin E\). \(\square\)

Solving the equation

Kepler’s equation has no closed-form solution, but Newton’s method converges rapidly: starting from \(E_0=M\) (or \(E_0=\pi\) for large \(e\) near aphelion), iterate

\[ E_{k+1}=E_k-\frac{E_k-e\sin E_k-M}{1-e\cos E_k}. \]

Because \(f(E)=E-e\sin E-M\) is increasing with \(f'(E)=1-e\cos E\ge1-e>0\), the root is unique. For \(e\) close to 1 and small \(M\), the equation is nearly cubic: \(M\approx(1-e)E+eE^3/6\), which is the elliptic counterpart of Barker’s equation for parabolas (Lesson 05 of the Caroline Herschel module).

1P/Halley (modern values): \(a\approx17.9\) AU, \(e\approx0.967\), \(P\approx76\) yr, \(q\approx0.59\) AU. At the 1759 return, perihelion fell on 13 March 1759.

Exercise 6.1Messier saw the comet on 21 January 1759, 51 days before perihelion. Using \(a=17.9\) AU and \(e=0.967\), find its distance from the Sun and its true anomaly that night.

Show solution
\(P=17.9^{3/2}=75.7\) yr \(=27{,}660\) d, so \(M=-2\pi\times51/27{,}660=-0.01158\) rad. Newton’s method from \(E_0=-0.1\) converges to \(E\approx-0.263\) rad. Then \(r=17.9(1-0.967\cos0.263)\approx1.18\) AU and \(\nu\approx-91^\circ\). The comet was just outside the Earth’s orbit and approaching the Sun.

Exercise 6.2For how many days in each orbit is Halley’s comet within 2 AU of the Sun? What fraction of its period is that?

Show solution
\(r=2\) when \(\cos E=(1-2/17.9)/0.967=0.919\), so \(E=0.407\) rad and \(M=E-e\sin E=0.407-0.967\times0.396=0.0243\). The interval is \(2M/2\pi\times27{,}660\approx213\) days, about 0.8 per cent of the period. A searcher has a window of a few months every 76 years.

Exercise 6.3Show that for \(e\) close to 1 and small \(E\), \(M\approx(1-e)E+E^3/6\). Use it to estimate \(E\) in Exercise 6.1 without iteration.

Show solution
\(\sin E\approx E-E^3/6\), so \(E-e\sin E\approx(1-e)E+eE^3/6\approx(1-e)E+E^3/6\). With \(1-e=0.033\) and \(M=-0.01158\), try \(E=-0.26\): \(0.033\times(-0.26)+(-0.26)^3/6=-0.00858-0.00293=-0.01151\). Close to \(M\), so \(E\approx-0.26\), in agreement with Newton’s method.
Lesson 07

Lexell’s comet and Jupiter

For the general reader. On 14 June 1770 Messier found a faint comet in Sagittarius. On 1 July it passed the Earth at about six times the Moon’s distance, the closest approach of any comet in recorded history, and crossed the sky so fast that its motion could be seen in an hour. Anders Johan Lexell, working in St Petersburg, showed that it moved on a 5.6-year ellipse and argued that Jupiter had thrown it onto that orbit in 1767. In 1779 Jupiter threw it off again, and it was never seen after.

Goals
  • Estimate the parallax and apparent motion of a very close comet.
  • Derive the Tisserand parameter from the circular restricted three-body problem.
  • Use it to predict the comet’s orbit after a Jupiter encounter.

D/1770 L1 (Lexell), orbit of 1770: \(a=3.15\) AU, \(e=0.786\), \(q=0.675\) AU, inclination \(1.55^\circ\), \(P=5.58\) yr. Closest approach to Earth on 1 July 1770 at about 0.015 AU (2.2 million km). After the 1779 encounter with Jupiter, Johann Karl Burckhardt found a perihelion distance of 3.33 AU.

Parallax

A body at distance \(d\) is seen from points on the Earth’s surface in directions that differ by up to the horizontal parallax \(\pi_h=R_\oplus/d\). For the Moon this is about a degree; for the Sun about \(8.8''\). A comet at 0.015 AU has a parallax of almost ten arcminutes, so observations from different places, and at different hours from the same place, give visibly different positions. Lexell had to correct for this before an orbit would fit.

The Tisserand parameter

Consider a comet of negligible mass moving under the Sun and a planet on a circular orbit of radius \(a_J\). In a frame rotating with the planet at angular speed \(n_J\), the Jacobi integral

\[ C_J=2\Big(\frac{GM_\odot}{r_1}+\frac{GM_J}{r_2}\Big)+n_J^2(x^2+y^2)-v_{\text{rot}}^2 \]

is exactly conserved, where \(r_1,r_2\) are the distances to the Sun and planet and \(v_{\text{rot}}\) the speed in the rotating frame.

Proposition (Tisserand’s criterion)

Far from the planet, the Jacobi integral reduces to a function of the comet’s heliocentric elements:

\[ T_J=\frac{a_J}{a}+2\cos i\,\sqrt{\frac{a}{a_J}(1-e^2)} \]

takes approximately the same value before and after a close encounter with the planet. Here \(i\) is the inclination to the planet’s orbital plane.

Sketch. Far from the planet, drop the \(GM_J/r_2\) term. The inertial velocity is \(\mathbf v=\mathbf v_{\text{rot}}+n_J\hat z\times\mathbf r\), so \(v_{\text{rot}}^2=v^2-2n_J(\hat z\cdot\mathbf r\times\mathbf v)+n_J^2(x^2+y^2)\). Substituting, \(C_J\approx 2GM_\odot/r-v^2+2n_Jh_z\), with \(h_z=\sqrt{GM_\odot a(1-e^2)}\cos i\). The vis-viva equation gives \(2GM_\odot/r-v^2=GM_\odot/a\). With \(n_J^2a_J^3=GM_\odot\), dividing by \(GM_\odot/a_J\) gives \(T_J\). \(\square\)

Tisserand’s parameter is how astronomers decide whether a newly found comet might be an old one in disguise: if two orbits have very different \(T_J\), no encounter with Jupiter can connect them.

Exercise 7.1Find the horizontal parallax of Lexell’s comet at 0.015 AU. Then, modelling both bodies as coplanar, with the comet at 1 AU on its 1770 orbit and the Earth on a circular orbit at 29.8 km/s, estimate the comet’s speed relative to Earth and its apparent angular motion per day at closest approach.

Show solution
\(d=2.24\times10^6\) km, so \(\pi_h=6378/2.24\times10^6=2.84\times10^{-3}\) rad \(\approx9.8'\). Vis-viva at \(r=1\): \(v=29.8\sqrt{2-1/3.15}=38.6\) km/s. The transverse component is \(h/r=29.8\sqrt{a(1-e^2)}=29.8\sqrt{1.205}=32.7\) km/s, so the radial component is \(\sqrt{38.6^2-32.7^2}=20.6\) km/s. Relative to Earth: \(\sqrt{(32.7-29.8)^2+20.6^2}\approx20.8\) km/s. If this velocity is perpendicular to the line of sight at closest approach, the angular rate is \(20.8/2.24\times10^6\) rad/s, about \(46^\circ\) per day, or nearly two degrees an hour.

Exercise 7.2Compute \(T_J\) for the 1770 orbit, with \(a_J=5.20\) AU. Assuming \(T_J\) and the low inclination are unchanged, find \(a\), \(e\) and \(P\) after the 1779 encounter, given Burckhardt’s \(q=3.33\) AU.

Show solution
\(T_J=5.20/3.15+2\sqrt{(3.15/5.20)(1-0.786^2)}=1.653+2\times0.481=2.615\). After: write \(a(1-e^2)=q(1+e)=q(2-q/a)\), and solve \(5.20/a+2\sqrt{q(2-q/a)/5.20}=2.615\). Bisection gives \(a\approx9.1\) AU, \(e=1-q/a\approx0.64\), \(P\approx28\) yr. With perihelion beyond Mars, the comet would be far too faint to see with eighteenth-century telescopes, as Lexell and Laplace concluded.
Lesson 08

How bright is a comet?

For the general reader. A planet shines by reflecting sunlight, so its brightness depends on its distances from the Sun and from us. A comet does more: as it nears the Sun its ices evaporate and the gas and dust spread into a coma thousands of times larger than the nucleus. Its brightness therefore rises steeply as it approaches the Sun. This is why comet hunters, Messier among them, searched the twilight skies, low in the west after sunset and in the east before dawn.

Goals
  • Derive the magnitude law for a sunlit body.
  • Use the empirical comet law with activity index \(n\).
  • Compute where in the sky comets near the Sun appear.
Proposition (magnitude of a sunlit body)

A body at distance \(r\) AU from the Sun and \(\Delta\) AU from the Earth, of fixed size and reflectivity, has apparent magnitude

\[ m=H+5\log_{10}r+5\log_{10}\Delta+\phi(\alpha), \]

where \(H\) is its magnitude at \(r=\Delta=1\) AU and zero phase angle and \(\phi\) is a phase correction. For comets the empirical law is

\[ m=H_0+5\log_{10}\Delta+2.5\,n\log_{10}r, \]

with activity index \(n\) typically between 2 and 6, often taken as 4.

Derivation. The flux received by the body is proportional to \(r^{-2}\), and the reflected flux reaching us to \(r^{-2}\Delta^{-2}\). In magnitudes this is \(2.5\log_{10}(r^2\Delta^2)=5\log_{10}r+5\log_{10}\Delta\). For a comet the amount of reflecting dust and fluorescing gas itself grows as the comet nears the Sun, roughly as a power of \(r\); writing the total as \(r^{-n}\) gives the second form, with \(n=2\) for an inert body. \(\square\)

Elongation

The elongation \(\varepsilon\) of a body from the Sun, seen from Earth at 1 AU, satisfies \(r^2=1+\Delta^2-2\Delta\cos\varepsilon\). A comet with \(r<1\) can never appear at elongation greater than \(\arcsin r\).

The second statement follows from the triangle Sun–Earth–comet: the elongation is largest when the line of sight is tangent to the sphere of radius \(r\) about the Sun. Mercury, at \(r\approx0.39\), never strays more than about \(23^\circ\) to \(28^\circ\) from the Sun. A comet at perihelion inside Mercury’s orbit is similarly confined to the twilight sky.

Exercise 8.1A comet has \(H_0=6\) and \(n=4\). Find its magnitude (a) at \(r=1.5\), \(\Delta=1.0\) and (b) at \(r=1.2\), \(\Delta=0.5\). Could Messier’s refractor of Exercise 4.1 show it in each case?

Show solution
(a) \(m=6+0+10\log_{10}1.5=7.8\). (b) \(m=6+5\log_{10}0.5+10\log_{10}1.2=6-1.5+0.8=5.3\). Both are within reach of an 11th-magnitude telescope, though a comet’s light is spread over a coma, so the practical limit for a diffuse comet is several magnitudes brighter than for a star.

Exercise 8.2A comet is at \(r=0.5\) AU. What is its greatest possible elongation from the Sun? If it is at \(\Delta=1.2\) AU, what is its elongation?

Show solution
Maximum: \(\arcsin0.5=30^\circ\). At \(\Delta=1.2\): \(\cos\varepsilon=(1+1.44-0.25)/2.4=0.9125\), so \(\varepsilon\approx24^\circ\). The comet would be seen only low in a bright twilight sky.
Lesson 09

What the Messier objects are

For the general reader. To Messier, all his objects were the same kind of thing: faint patches that were not comets. We now know that his list samples almost every kind of object in the universe that a small telescope can show. It contains clusters of young stars in our Galaxy, ancient globular clusters orbiting it, clouds of glowing gas where stars form, shells thrown off by dying stars, the debris of a supernova, and about forty galaxies millions of light years away.

Goals
  • Classify the 110 objects and estimate their distances.
  • Convert angular size to physical size.
  • Use the virial theorem to weigh a globular cluster.

The 110 objects by type (modern classification): about 40 galaxies; 29 globular clusters; 26 open clusters; 4 planetary nebulae; 7 diffuse nebulae, most with embedded clusters; 1 supernova remnant (M1); 1 Milky Way star cloud (M24); 1 asterism (M73) and 1 double star (M40). M102 is disputed: it may be a duplicate of M101 or the galaxy NGC 5866.

Distance modulus and size

An object of absolute magnitude \(M\) at distance \(d\) parsecs has \(m-M=5\log_{10}(d/10\,\text{pc})\), neglecting absorption. An object of angular size \(\theta\) radians at distance \(d\) has linear size \(\ell=\theta d\).

Theorem (virial theorem)

For a gravitationally bound system of point masses in a steady state, the time averages of the total kinetic energy \(K\) and gravitational potential energy \(W\) satisfy \(2\langle K\rangle+\langle W\rangle=0\).

Proof. Let \(G=\sum_i\mathbf p_i\cdot\mathbf r_i\). Then \(\dot G=\sum m_iv_i^2+\sum\mathbf F_i\cdot\mathbf r_i=2K+W\), where the second sum equals \(W\) because the gravitational potential is homogeneous of degree \(-1\) in the positions (Euler’s theorem on homogeneous functions). For a bound system \(G\) stays bounded, so the time average of \(\dot G\) over a long interval tends to zero. \(\square\)

For a cluster of mass \(M\), one-dimensional velocity dispersion \(\sigma\) and characteristic radius \(R\), \(K=\tfrac32M\sigma^2\) and \(W=-\alpha GM^2/R\) with \(\alpha\) of order 1, so \(M\approx3\sigma^2R/(\alpha G)\). The precise factor depends on the cluster’s density profile and on which radius is used, so such masses are good to a factor of about two.

Exercise 9.1The Pleiades (M45) lie at about 136 pc and the cluster’s brightest stars span about \(2^\circ\). Find the distance modulus and the cluster’s linear size. The brightest star, Alcyone, has \(m\approx2.9\). What is its absolute magnitude?

Show solution
\(5\log_{10}13.6=5.67\). Size: \(136\times(2\pi/180)\approx4.7\) pc. Alcyone: \(M\approx2.9-5.7=-2.8\), roughly a thousand times the Sun’s luminosity (\(M_\odot\approx4.8\)).

Exercise 9.2Take a globular cluster such as M13 with \(\sigma\approx7\) km/s and \(R\approx15\) pc, and \(\alpha=1\). Estimate its mass in solar masses.

Show solution
\(M\approx3\sigma^2R/G=3\times(7\times10^3)^2\times(15\times3.086\times10^{16})/(6.674\times10^{-11})\approx1.0\times10^{36}\) kg \(\approx5\times10^5\,M_\odot\). Half a million stars, in a patch that Messier’s telescope could not resolve into a single one.

Exercise 9.3In 2019 the Event Horizon Telescope imaged the shadow of the black hole in M87, of mass about \(6.5\times10^9\,M_\odot\) at 16.8 Mpc. For a non-rotating black hole the shadow has diameter \(2\sqrt{27}\,GM/c^2\). Find its angular diameter.

Show solution
\(GM_\odot/c^2=1.477\) km, so \(GM/c^2=9.6\times10^9\) km. The diameter is \(2\times5.196\times9.6\times10^9\approx1.0\times10^{11}\) km. At \(16.8\times3.086\times10^{19}\) km this subtends \(1.9\times10^{-10}\) rad \(\approx40\) microarcseconds, close to the measured ring of about 42 microarcseconds. Messier’s M87 contains the first black hole ever imaged.
Lesson 10

M1: from nebula to pulsar

For the general reader. The first object in Messier’s catalogue, found in 1758 while he waited for Halley’s comet, has become one of the most studied objects in the sky. John Bevis had seen it in 1731. In the 1840s Lord Rosse drew it with crab-like arms and gave it its name. In the twentieth century astronomers found it was expanding, traced the expansion back to a star recorded by Chinese astronomers in 1054, and in 1968 found at its centre a neutron star spinning thirty times a second. It is the remnant of a supernova.

Goals
  • Date the explosion from the nebula’s expansion.
  • Derive the spin-down luminosity and characteristic age of a pulsar.
  • Estimate a pulsar’s magnetic field.

Expansion age

A feature moving outward at speed \(v\) perpendicular to the line of sight, at distance \(d\), has proper motion \(\mu=v/(4.74\,d)\) arcseconds per year, with \(v\) in km/s and \(d\) in parsecs. If the nebula has expanded uniformly from a point, a feature at angular distance \(\theta\) from the centre is \(\theta/\mu\) years old. In 1921 John Charles Duncan compared photographs taken years apart and found the nebula expanding; in the same year Knut Lundmark noted that it lay near the position of the guest star of 1054. Careful modern measurements trace the expansion back to several decades after 1054, which implies that the nebula has been accelerated, pushed from inside by the pulsar wind.

The pulsar

The Crab pulsar has period \(P\approx33.7\) ms, increasing at \(\dot P\approx4.2\times10^{-13}\) s per second.

Proposition (spin-down)

A rotating neutron star with moment of inertia \(I\) and period \(P\) loses rotational energy at the rate

\[ L_{\text{sd}}=-\frac{d}{dt}\Big(\frac12I\Omega^2\Big)=\frac{4\pi^2I\dot P}{P^3}. \]

If it spins down as a magnetic dipole, \(\dot\Omega\propto-\Omega^3\), and if its initial period was much shorter than \(P\), its age is \(\tau_c=P/(2\dot P)\).

Proof. With \(\Omega=2\pi/P\), \(\dot\Omega=-2\pi\dot P/P^2\), so \(-I\Omega\dot\Omega=I(2\pi/P)(2\pi\dot P/P^2)=4\pi^2I\dot P/P^3\). For the age: \(\dot\Omega=-k\Omega^3\) means \(P\dot P=2\pi\cdot\)const, so \(d(P^2)/dt=2P\dot P\) is constant. Integrating from \(P_0\) to \(P\), \(t=(P^2-P_0^2)/(2P\dot P)\approx P/(2\dot P)\). \(\square\)

The dipole formula for spin-down also gives the surface magnetic field at the equator, \(B\approx3.2\times10^{19}\sqrt{P\dot P}\) gauss, with \(P\) in seconds, for \(I=10^{45}\) g cm\(^2\) and radius 10 km.

Exercise 10.1On photographs taken around 1950, a filament \(180''\) from the centre of the nebula moves outward at \(0.20''\) per year. Find the expansion age and the year of origin. At a distance of 2.0 kpc, how fast is the filament moving?

Show solution
Age \(=180/0.20=900\) years, giving an origin around 1050, close to the supernova of 1054. Speed: \(v=4.74\,\mu d=4.74\times0.20\times2000\approx1{,}900\) km/s. Averaging over many filaments gives dates somewhat later than 1054, the sign of acceleration described above.

Exercise 10.2With \(I=10^{38}\) kg m\(^2\), \(P=0.0337\) s and \(\dot P=4.21\times10^{-13}\), compute \(L_{\text{sd}}\), \(\tau_c\) and \(B\). Compare \(\tau_c\) with the true age and \(L_{\text{sd}}\) with the Sun’s luminosity, \(3.8\times10^{26}\) W.

Show solution
\(L_{\text{sd}}=4\pi^2\times10^{38}\times4.21\times10^{-13}/0.0337^3\approx4.3\times10^{31}\) W, about \(10^5\) solar luminosities: enough to power the nebula’s glow. \(\tau_c=0.0337/(8.42\times10^{-13})=4.0\times10^{10}\) s \(\approx1{,}270\) years, against a true age of about 970 years in 2024; the difference means the initial period was not negligible, or the braking is not purely dipolar. \(B\approx3.2\times10^{19}\sqrt{1.42\times10^{-14}}\approx3.8\times10^{12}\) G.
Part III · Life and legacy
Lesson 11

Academy, catalogue, Revolution, 1771–1817

As Astronomer of the Navy from 1771 Messier had a salary, a royal pension and the observatory on the Hôtel de Cluny. He kept observing comets, eclipses, occultations and transits, and contributed papers regularly to the Academy’s Mémoires. In 1775 Lalande gave his name to a small constellation he drew on a celestial globe, Custos Messium, the guardian of the harvests, a pun on Messier’s name. It did not survive into modern star atlases.

From about 1774 he worked with Pierre Méchain, a younger astronomer and later director of the Paris Observatory, who found many new nebulae and passed them to Messier. The catalogue grew to 68 objects in the Connaissance des Temps published in 1780, and to 103 in the edition for 1784, completed in 1781. The last entries, many of them Méchain’s, were added hastily and some were not checked. Around 1781 Messier fell about eight metres into an ice cellar and was badly injured. He recovered only slowly, and the catalogue stopped growing.

Revolution and Empire

The Revolution swept away the institutions that paid him. His pension and salary ended, the Academy was abolished in 1793, and he was left poor. He continued to observe. When the Republic rebuilt its scientific institutions, he was appointed to the Bureau des Longitudes in 1796 and became a member of the Institut. In 1806 Napoleon gave him the cross of the Legion of Honour. Two years later Messier published, at his own expense, a memoir on the great comet of 1769, connecting it with the year of the Emperor’s birth. The flattery shows how much he still needed patronage in his late seventies.

His eyesight and health failed in his last years. In 1815 a stroke left him partly paralysed. He died at his home in Paris on 12 April 1817, aged 86, and was buried in the Père Lachaise cemetery. Jean-Baptiste Delambre wrote an account of his life for the Institut.

DiscussionMessier valued his comets far above his catalogue, and his contemporaries agreed. Today the catalogue is what everyone knows. What does this say about how the value of scientific work is judged, and by whom? Are there modern equivalents, data sets compiled as by-products that outlive the research they served?

Lesson 12

The catalogue after Messier

Soon after 1781 William Herschel obtained a copy of Messier’s final catalogue. Its hundred-odd objects showed him that the nebulae were a field worth surveying, and with Caroline Herschel he swept the sky from Slough for twenty years and catalogued about 2,500. Their work, extended by John Herschel to the southern sky, became the General Catalogue of 1864 and the New General Catalogue of J. L. E. Dreyer in 1888, which gives every Messier object an NGC number. Messier’s list nevertheless kept its place, because it gathers the brightest and most striking objects visible from northern latitudes.

The list also grew. In the twentieth century historians found, in Messier’s and Méchain’s notes, observations of objects that belonged in it. Camille Flammarion added M104 in 1921; Helen Sawyer Hogg added M105 to M107 in 1947; Owen Gingerich added M108 and M109 in 1960; and Kenneth Glyn Jones added M110, a companion of the Andromeda galaxy, in the 1960s. The standard list has 110 entries. In early spring the Sun lies in a part of the sky with no Messier objects, and amateur astronomers attempt a Messier marathon, observing all of them between dusk and dawn.

Legacy

  • Comets: about 13 discoveries now credited to him, a careful recovery of Halley’s comet in 1759, and the discovery of Lexell’s comet, the closest comet ever recorded.
  • The catalogue: 110 objects, the first systematic list of non-stellar objects, now a standard observing list and a sampler of astrophysics, from star formation to supermassive black holes.
  • Honours: Fellow of the Royal Society (1764), member of the Paris Academy (1770), of the Bureau des Longitudes (1796) and of the Legion of Honour (1806). The lunar crater Messier and the asteroid 7359 Messier are named after him.

Final projectChoose one Messier object. Write a two-page account for a strong undergraduate that gives Messier’s own description, explains what the object is now known to be, and computes two quantities from Part II (for example a distance modulus, a linear size, a virial mass, an expansion age or a limiting magnitude).

Reference

Sources and further reading

Biography

  1. Charles Messier, Wikipedia (with the list of comets credited to him).
  2. Observatoire de Paris, Bibnum, Charles Messier (1730–1817).
  3. Charles Messier and Charles Messier, French and German Wikipedia.
  4. J. J. O’Connor and E. F. Robertson, Charles Messier, MacTutor History of Mathematics.
  5. J.-B. J. Delambre, Histoire de l’astronomie au dix-huitième siècle, Paris, 1827 (contains his account of Messier).
  6. K. G. Jones, Messier’s Nebulae and Star Clusters, 2nd ed., Cambridge University Press, 1991.
  7. Encyclopaedia Britannica, 11th ed. (1911), article “Messier, Charles”.

Astronomy

  1. C. Messier, Catalogue des nébuleuses et des amas d’étoiles, Mémoires de l’Académie royale des sciences for 1771 (published 1774); final version in Connaissance des Temps for 1784.
  2. Messier object, Wikipedia (editions, additions, classification).
  3. Lexell’s Comet and Crab Nebula, Wikipedia.
  4. Event Horizon Telescope Collaboration, First M87 Event Horizon Telescope results. I. The shadow of the supermassive black hole, arXiv:1906.11238.
  5. J. M. A. Danby, Fundamentals of Celestial Mechanics, 2nd ed., Willmann-Bell, 1988 (Kepler’s equation, Tisserand’s criterion).
  6. C. D. Murray and S. F. Dermott, Solar System Dynamics, Cambridge University Press, 1999 (Jacobi integral and Tisserand parameter).
  7. D. R. Lorimer and M. Kramer, Handbook of Pulsar Astronomy, Cambridge University Press, 2005 (spin-down and characteristic age).