Uranus to Neptune: the eye, the chart and the equation

On the night of 13 March 1781, between ten and eleven o’clock, William Herschel was working through faint stars near H Geminorum at Bath. He was looking for stellar parallax, not for planets. One object in the field looked larger than its neighbours. He had been observing at a magnification of 227. He swapped in eyepieces giving 460 and then 932, on the reasoning that a star is a point whose image should not grow in proportion, while a real disk should. The object grew. On 26 April his report was read to the Royal Society under the title Account of a Comet (Herschel 1781).

It was not a comet. It was the first planet found since antiquity, and the first found with a telescope. Herschel did not know that, and the title is an honest record of what he could claim. His instrument was a reflector of his own construction, not the giant telescope he built later. What he had done was to notice an anomaly in morphology, test it with a controlled experiment, and then measure the object’s position against reference stars on later nights until its motion could be handed to people who computed orbits.

This chapter follows the discoveries of the century and a quarter that came after: Uranus, the first asteroids, Neptune, the satellites of the outer planets, and the photographic asteroids of the 1890s. It also follows one failure, the planet Vulcan. The reason to tell these stories in a book about algorithms is that every logical stage of a modern survey is already present. Someone must decide where to look. An observation must be collected. A body must be distinguished from stars and instrumental defects. Its motion must be measured. An inverse problem must be solved, a return predicted, and an independent confirmation obtained. In 1800 each of these stages was a person, a notebook, a printed chart or a table of logarithms. The decisive innovations were often in the connections between the stages, not in a larger telescope.

Uranus: a disk that should not be there

The magnification test

Herschel’s test deserves a modern statement, because it is the first deliberate use of image shape as a detection statistic in this story. A body of diameter $D$ at distance $\Delta$ subtends an angle $\theta_d \simeq D/\Delta$. What the eye sees is that disk convolved with the telescope’s diffraction pattern, the blur from atmospheric seeing, and the optical aberrations. A star is effectively a point, so its image is the blur alone. Raising the magnification enlarges both images on the retina, but a resolved disk keeps its intrinsic size relative to the field while a star’s apparent size is set by the blur and the eye’s response.

In modern terms one would compare two hypotheses for the same pixels $I$: a point source with flux $F$ seen through a point-spread function $P$, and a disk $D_{\theta_d}$ of angular diameter $\theta_d$ seen through the same $P$. The log-likelihood ratio

measures how much better a resolved disk explains the data. Nobody computed this in 1781. It makes explicit what information Herschel’s experiment was after. It also shows the test’s limits. Comets have comae, close double stars blur together, a telescope out of focus makes every star look fat, and bad seeing erases the size information altogether. A better classifier cannot recover a disk that seeing has washed out. Herschel’s report correctly treated shape as a reason to keep looking, not as proof (Herschel 1781).

From comet to planet

The second test was motion. Herschel’s report gives micrometer distances from reference stars on repeated dates, notes the changing apparent diameter, and records the practical difficulty of measuring a moving object. Those tables turned a visual impression into a communicated, testable object. The planetary interpretation came from the orbit computers, not from the first night. Gauss, writing much later in the preface to his Theoria motus, notes that the new body could not be reconciled with a parabolic comet orbit, while a circular hypothesis gave a useful first approximation (Gauss 1809).

The detection lesson is easy to miss. Uranus is bright enough that earlier astronomers had recorded it as a star. Collecting its photons was never the missing capability. Recognising a moving, extended Solar System object, and communicating enough positions to determine what it was, was. Once Uranus was identified, those earlier “star” positions gained a new value as a long astrometric baseline. Sixty years later the same old positions would be part of the evidence for Neptune (Adams 1846).

A rule for where to look

Uranus seemed to confirm a numerical pattern in the planetary distances. Johann Daniel Titius had introduced it in a translation of Charles Bonnet, and Johann Elert Bode promoted it from 1772 (Piazzi 1802). In astronomical units,

The limit $n=-\infty$ gives 0.4 for Mercury. Then $n=0,1,2$ give 0.7, 1.0 and 1.6 for Venus, Earth and Mars, $n=4$ gives 5.2 for Jupiter and $n=5$ gives 10.0 for Saturn. The slot $n=3$, at 2.8, was empty. The slot $n=6$ gives 19.6, and Uranus’s mean distance of about 19.2 fitted it well. The next slot, $n=7$, would be 38.8 (Bovaird & Lineweaver 2013). These are outputs of the formula, not measured distances.

The rule has no force law behind it. It does not follow from Newtonian gravity, and nobody could say why it should hold. In the language of Chapter 1 it was a prior: a probability distribution over semimajor axes concentrated at a few predicted slots. A sharp prior shrinks the search, which is why it was useful. But every conclusion drawn from it is conditional on an unproved architecture. Neptune, at about 30 au, would break the rule. In the meantime it organised a search.

Piazzi’s account describes an association of 24 astronomers who divided the zodiac among themselves to look for the missing planet between Mars and Jupiter (Piazzi 1802). The coordination is associated with Franz Xaver von Zach and is often called the “celestial police.” It is a pleasing story that the police then found Ceres. They did not.

Ceres: lost and found

A star that moved

On 1 January 1801, at Palermo, Giuseppe Piazzi was measuring stars for a catalogue. He was looking for star 87 in Nicolas-Louis de Lacaille’s zodiacal catalogue and also measured a preceding object of about the eighth magnitude. The next evening its transit time and its zenith distance did not agree with the first night. He first suspected an error of his own. The following nights showed systematic motion. In his own account, suspicion became certainty on the third night, and he waited until the fourth before speaking of it (Piazzi 1801).

This was precision transit astronomy, and the detail matters. Piazzi’s instrument was the Palermo Circle, built by Jesse Ramsden. Its five-foot vertical circle was an angle-measuring structure, not a five-foot aperture. The telescope magnified 50 times on the circle and 80 on the transit instrument. An object crossed the field in about two minutes. Right ascension came from the clock time of meridian passage: a clock error $\delta t$ produces a position error $\delta\alpha \simeq \omega_\oplus\,\delta t$, where $\omega_\oplus$ is the sidereal rotation rate. Declination came from the measured zenith distance after correcting for latitude and refraction. Piazzi did not notice an unfamiliar dot. He noticed that two measured numbers had changed. The discovery was a residual against a catalogue (Piazzi 1801; O’Connor & Robertson).

Clouds immediately created the first association problem in this book. Returning to the field after a gap, Piazzi found several similar stars, measured them all, and used the next night’s motion to decide which one was his object. That is a linking problem with three candidates and one extra epoch to break the tie. It is the same problem that Chapter 12 treats with millions of candidates.

He followed the object to 11 February. Then it approached the Sun in the sky and he fell ill. Gauss later summarised what was left: a span of 41 days in which the body had moved through a geocentric arc of only three degrees. That is a baseline and an angle, not a claim of 41 nightly observations (Gauss 1809).

The angles-only problem

Three degrees of arc in six weeks is a short arc. To see why recovering the body was hard, write down what an observation measures. Let $\mathbf r(t)$ be the body’s heliocentric position and $\mathbf R_i$ the observer’s at time $t_i$. An angular measurement gives only the unit vector

where light time has been neglected. Each observation supplies two angles and loses the range $\rho_i$. Three observations supply six angles, and an orbit has six elements, so the count balances. The ranges are the unknowns that the dynamics must supply.

Over a short arc the motion obeys Kepler’s two-body equation, $\ddot{\mathbf r} = -\mu\,\mathbf r/r^3$, with $\mu$ the Sun’s gravitational parameter. Any position on the orbit can be written in terms of the position and velocity at a middle time:

where $f$ and $g$ are the Lagrange coefficients, known functions of the time intervals and the orbit. Eliminating the unknown velocity $\dot{\mathbf r}_2$ gives

Substituting $\mathbf r_i = \mathbf R_i + \rho_i\mathbf u_i$ turns this into a small linear system in the three ranges, once $c_1$ and $c_3$ are known. They are not known exactly, because they depend on the orbit. The standard modern treatment truncates $f$ and $g$ to a few terms, which reduces the problem to an eighth-degree polynomial in the middle heliocentric distance, and then iterates. Several positive roots may need to be tested (Orekit IodGauss). This textbook polynomial belongs to modern implementations. It is not a claim about Gauss’s notebook.

Once a preliminary orbit exists, every observation refines it. With residual vector $\mathbf d = \mathbf y - \mathbf h(\theta_0)$ between measured positions $\mathbf y$ and predicted positions $\mathbf h$, Jacobian $J = \partial\mathbf h/\partial\theta$, and measurement covariance $C$, one solves

whose normal equations are $(J^{\mathsf T}C^{-1}J)\,\delta\theta = J^{\mathsf T}C^{-1}\mathbf d$. The matrix $(J^{\mathsf T}C^{-1}J)^{-1}$ is the local covariance of the fitted elements. On a short arc it is only locally meaningful: several distinct orbits can fit, and the uncertainty region is long and curved. Counting six angles does not guarantee a well-conditioned orbit. Chapter 13 develops this for modern surveys.

Gauss and the recovery

Carl Friedrich Gauss, aged 24, took up the problem in October 1801. His achievement was an ephemeris accurate enough to find an unresolved eighth-magnitude object after months behind the Sun, without imposing a circular orbit as the final answer. Piazzi reports that Gauss’s predicted positions differed from the competing predictions by about nine degrees in longitude (Piazzi 1802). That is a change in where to point the telescope, not a refinement of residuals.

The recovery itself was a sequence, not a moment. Gauss’s preface gives a footnote observation by von Zach on 7 December 1801. Piazzi’s more detailed account explains that the candidate stayed uncertain through bad weather and became convincing only on 31 December. Heinrich Olbers also recovered the body around the turn of the year (Gauss 1809; Piazzi 1802). A single celebrated “rediscovery date” hides observation, suspicion and demonstrated motion.

A second compression should also be resisted: that Gauss invented least squares and used the textbook “Gauss method” to find Ceres. His Theoria motus of 1809 says the methods had changed so much since the original computation that scarcely any resemblance remained. The book sets out a mature orbit-determination framework with a probabilistic treatment of errors and least squares. It is not a record of the 1801 calculation (Gauss 1809).

Charts as detectors

Ceres was followed within a few years by Pallas, found by Olbers on 28 March 1802, Juno by Karl Ludwig Harding on 1 September 1804, and Vesta, again by Olbers, on 29 March 1807 (MPC discovery circumstances). Then nothing for 38 years. The fifth, Astraea, came on 8 December 1845, from Karl Ludwig Hencke’s patient visual work.

The long gap did not mean the population was sparse. It meant that the method was expensive. Finding a star-like object whose position changes requires comparing a field with something: a chart, a catalogue, or the observer’s own memory of the sky. An observer who has spent years learning a region can notice a newcomer. Anyone else needs a map.

Friedrich Wilhelm Bessel proposed the Berliner Akademische Sternkarten, star charts of the zodiac, in 1824. Remchin and Schrimpf’s archival study stresses that the existing catalogues together held roughly 50,000 stars, but did not form a complete catalogue to any uniform magnitude (Remchin & Schrimpf 2019). That matters for how a chart works as a detector. An object missing from the chart could be a moving body, an ordinary star the chart omitted, or a chart error. The false-candidate rate of the method is set by the completeness of the reference.

The figure shows what happened once Hencke’s success and the charts combined. The count of numbered main-belt asteroids, grouped by discovery year, goes from 5 in 1845 to 36 by 1855 and over 200 by 1880. The map had become part of the detector. It externalised memory, standardised the comparison, and let someone who had not spent years learning a field search it. The labour also divided. Chart construction, discovery, position measurement and ephemeris calculation were done by different people. Gauss’s students and other computers specialised in the orbits of particular new bodies, and the processing network mattered as much as the discoverer’s eye (Remchin & Schrimpf 2019).

In modern terms there are two operations here. In catalogue differencing, one looks for detections that have no stationary counterpart in the reference. In repeated-field comparison, one tests whether a detection’s position changes. The first is economical when the reference is complete. The second repairs the reference’s omissions and resolves which of several faint objects is the mover. Matching $N$ detections against $M$ chart entries naively costs $O(NM)$ comparisons. Sorting both into coordinate strips, or using a spatial index, brings it to about $O(M\log M + N\log M + K)$, with $K$ the number of returned neighbours. A printed map, ordered by right ascension in hour-wide sheets, is a human spatial index. The direct descendant of the method is image differencing against a reference (Chapter 10).

Neptune: a planet found by its pull

The trouble with Uranus

By the 1820s Uranus was misbehaving. Alexis Bouvard’s tables of its motion could not reconcile the old pre-discovery positions, the ones recorded by observers who had taken it for a star, with the newer observations. Throwing out the old data did not help, because the new observations then drifted away from the tables too. John Couch Adams’s account gives the size of the problem: after his corrections, the longitude discrepancies of a no-new-planet theory were about +92.7 arcseconds in 1712, −47.6 in 1750, +24.2 in 1804 and −66.6 in 1840 (Adams 1846). A minute and a half of arc is small, a twentieth of the Moon’s diameter, but it was far outside the errors of the measurements.

Two young mathematicians, independently, asked whether an unseen planet beyond Uranus could cause the discrepancy, and where it would have to be.

The physics of the inverse problem

Work in heliocentric coordinates. If Uranus is at $\mathbf r$ and a hypothetical outer planet of mass $m_p$ is at $\mathbf r_p$, the acceleration of Uranus is

where $\mu_U = G(M_\odot + m_U)$ and $\mathbf a_{\rm known}$ collects the perturbations by Jupiter, Saturn and the rest. The bracket has two terms. The first is the direct pull of the new planet on Uranus. The second, the indirect term, is the pull of the new planet on the Sun, which must be subtracted because the coordinate origin is the Sun. Dropping it is a serious physical error.

If the perturber is far away, $r \ll r_p$, expand the direct term in powers of $r/r_p$. The leading (uniform) piece cancels the indirect term exactly, because a uniform field accelerates the Sun and Uranus equally. What is left is a tidal acceleration of order $Gm_p r/r_p^3$. A distant planet does not pull Uranus away from the Sun. It stretches the Sun–Uranus system. This is why perturbation searches constrain combinations like $m_p/r_p^3$ more readily than mass and distance separately, a degeneracy that returns in Chapter 4 and Chapter 8. For Uranus and Neptune the distances are too similar for the leading tidal term alone to reproduce the historical problem (Rodríguez-Moris & Docobo 2024).

A rough scale helps. Using the modern values in the next table, Neptune has about $5.15\times10^{-5}$ of the Sun’s mass. At conjunction it is about $30.1 - 19.2 \approx 11$ au from Uranus, while the Sun is about 19 au away. The ratio of Neptune’s direct pull to the Sun’s is about $5.15\times10^{-5}\times(19.2/10.9)^2 \approx 1.6\times10^{-4}$ (an estimate, using modern distances). A perturbing acceleration of one part in six thousand, acting for decades and modulated by the planets’ changing relative positions, produced the arc-minute residuals.

A second subtlety is that the observed residual is not simply the response to the perturber. Linearise the motion about a reference orbit $\mathbf r_0(t)$:

with $q$ the perturber’s parameters. The solution is a homogeneous part, which is just a small change in Uranus’s own initial orbit, plus a particular part driven by the perturber. Uranus’s orbit is fitted to the same data. So a substantial part of any perturbation is absorbed into a slightly different Uranus orbit and never shows up as a residual. Le Verrier understood this and solved for Uranus’s element corrections and the new planet’s parameters together (Le Verrier 1846). This is also why integrating an orbit with and without Neptune from identical initial conditions does not reproduce the residuals of a refitted no-Neptune ephemeris (Krajnović 2016).

In modern language, let the columns of a matrix $H$ hold the derivatives of the predicted longitudes with respect to Uranus’s own elements (the nuisance directions), and let $S(q)$ be the template signal from a perturber with parameters $q$ and unit mass. After whitening by the measurement covariance $C$, project out the nuisance directions:

With $d = P_\perp C^{-1/2}\,\delta\lambda$ the projected longitude residuals and $s = P_\perp C^{-1/2}S(q)$ the projected template, the best non-negative mass and the improvement in fit are

This is a matched filter in which the noise includes the unknown orbit of the planet being perturbed. Only the part of the perturber’s signal orthogonal to Uranus’s own orbit changes is detectable. Le Verrier did not write this statistic, and a modern search would also need a false-alarm calibration that accounts for maximising over $q$. But it is the structure of what he did (Rodríguez-Moris & Docobo 2024). Chapter 7 develops it.

Adams and Le Verrier

Adams’s own account describes a first calculation in 1843 that assumed a circular orbit at twice Uranus’s mean distance, which he took as 38.4 au. This is close to, but not the same as, the Titius–Bode value of 38.8. Later calculations allowed eccentricity and included further reductions. He formed 21 equations from groups of observations each spanning three years, checked the known perturbations, applied corrections including one to Jupiter’s mass, expanded the perturbations to second order in the eccentricities, and solved by least squares and successive substitution (Adams 1846; Airy 1846). His results were communicated privately in 1845, and a revised hypothesis in 1846. They did not lead to the discovery.

Urbain Le Verrier’s sequence was public and operationally decisive. He presented a reassessment of Uranus’s theory on 5 November 1845, an argument for an exterior planet with an approximate position on 1 June 1846, and refined elements with a search prescription on 31 August 1846 (Krajnović 2016). George Airy’s description of Le Verrier’s final analysis specifies 33 equations of condition in nine unknowns: five for the new planet and four corrections to Uranus (Airy 1846). These are reduced equations and model parameters, not raw observations. Le Verrier also explored which solutions remained admissible rather than publishing only a single best point. Airy describes varying the elements until the observational error became unacceptable, forming an envelope of possible positions. The envelope was conditional on the allowed eccentricities and periods. It was not a calibrated confidence region in the modern sense.

The historical computation used analytic expansions, periodic terms in longitude, and elimination of unknowns by hand. Le Verrier’s printed installment in Astronomische Nachrichten shows numerical equations of condition and an elimination argument in which he first removes Uranus’s corrections using selected epochs, then constrains combinations of the perturber’s mass, longitude, eccentricity and perihelion (Le Verrier 1846). It was not “measure the acceleration and point along it.”

23 September 1846

Le Verrier wrote to Johann Gottfried Galle at the Berlin Observatory on 18 September. The letter arrived on 23 September (Galle 1877). That night Galle and Heinrich d’Arrest went to the observatory’s refractor.

They first looked for what Le Verrier’s mass and distance implied: a resolved disk. Galle’s later account says this gave no certainty. The real disk is only slightly larger than two arcseconds, while the prediction suggested more than three. The morphology test that had worked for Herschel was marginal here. D’Arrest suggested using the Academy’s star charts. They found that Carl Bremiker’s recently completed sheet for Hora XXI covered the region, returned to the refractor, and compared the sky against the chart. An object of about the eighth magnitude was not on the chart. The following night it had moved (Galle 1877).

The search took a single night because three things came together: a localised prediction, a new and good reference chart of exactly that region, and a second epoch to confirm motion. Take away the chart and they would have been back to judging marginal disks one star at a time.

Good pointing, wrong orbit

The prediction was close on the sky. Krajnović tabulates Le Verrier’s final predicted true heliocentric longitude at the discovery epoch as 325°58′, against 326°57′ for Neptune, a difference of 59 arcminutes. (His narrative elsewhere says “one minute of arc,” which conflicts with his own table and should not be repeated.) These are heliocentric longitudes, not the separation seen from Berlin (Krajnović 2016).

The orbit was much worse than the pointing. Airy gives Le Verrier’s predicted semimajor axis as 36.154 au, eccentricity 0.10761, and mass about 1/9300 of the Sun’s (Airy 1846). The comparison below is from Krajnović’s Table 2. Walker’s column is an early orbit fitted to observations of the planet after its discovery.

The predicted orbits put the planet 20% too far away and made it twice as massive as it is. Mass and distance trade against each other in the perturbation, so the residuals constrained the direction to the perturber far better than its distance or mass. The inference was a good localiser and a poor estimator of the global orbit, conditioned on assumptions about distance that the Titius–Bode prior had encouraged.

Adams’s fit shows the same thing from another angle. After including his planet, his report gives a largest residual of 2.35 arcseconds for observations from 1781 to 1840. But the residuals then grew: 12.40 and 10.18 arcseconds in 1845 under his two distance hypotheses (Adams 1846). The wrong orbit fitted the fitted interval and began to fail outside it. Any modern reanalysis has to state which observations were used and how they were weighted, because these numbers depend on those choices.

Vulcan: the same method, a different answer

Le Verrier tried the strategy again, inward. Mercury’s perihelion advanced faster than the known planets could explain. A contemporary review by J. R. Hind gives the unexplained excess as 38 arcseconds per century (Hind 1876). Le Verrier proposed an unseen mass inside Mercury’s orbit, either one planet or a group.

An intramercurial planet would be hard to see in the Sun’s glare. Two observing methods were available. The first was to watch the solar disk for a small dark spot crossing it. This gives up sky coverage for contrast: the planet would be a silhouette against a uniform bright background, visible only during a transit. The second was to search near the Sun during total eclipses, which removes the glare for a few minutes.

Edmond Lescarbault reported a dark spot crossing the Sun on 26 March 1859. This encouraged the planet interpretation, but Hind’s own review shows how hard it became to fit the scattered reports onto one orbit. In March 1877 S. J. Perry kept a careful watch of the Sun with an 8-inch achromatic telescope through predicted transit windows, logging observing hours and interruptions for weather. His conclusion was appropriately narrow. No transit occurred while he was adequately observing. He did not claim to have excluded all interior orbits (Perry 1877; Capron 1877). Perry’s logged gaps are exactly the input a modern analysis needs: a non-detection excludes only the trajectories that would have crossed the observed disk while someone was looking.

During the total eclipse of 29 July 1878 James Craig Watson and Lewis Swift reported unfamiliar objects near the Sun. A contemporary Nature report sets out C. H. F. the objection by Peters that known stars could account for Watson’s objects, and the dispute over the precision of Watson’s temporary setting circles (Nature 1879). An apparent new star near the eclipsed Sun is especially vulnerable to catalogue and pointing errors: the observer has a few minutes, uncalibrated pointing, and no second epoch.

By 1895 Simon Newcomb declined to put a hypothetical group of unknown masses into his planetary tables, and discussed alternatives and systematic errors (Newcomb 1895). He noted that changing the Sun’s shape enough to affect Mercury would also change the nodes and inclinations of other planets. That is a cross-planet consistency test, and it is the right one. A real planet must predict compatible transits, sky positions and effects on other bodies.

The anomaly was real. General relativity later explained it. In relativity the perihelion advances per orbit by

where $a$ and $e$ are Mercury’s semimajor axis and eccentricity. With standard modern constants this is about $5.0\times10^{-7}$ radians per orbit, or about 43 arcseconds per century over Mercury’s roughly 415 orbits a century (my calculation). That is close to, but not the same as, Hind’s 38. The excess attributed to an unknown cause depended on how well the known perturbations had been computed, and the historical figure should not be silently replaced by the modern one.

The lesson is precise. A residual after removing known effects does not identify its cause. An interior planet, a distributed mass, a non-spherical Sun, a change in the law of gravity, or errors in reduction can all produce part of the same signature. Neptune and Vulcan used the same mathematics. One produced a planet within a degree of its predicted place. The other produced a decades-long chase after a body that did not exist. The difference was not in the inversion. It was that Neptune was then seen, with a chart, and moved.

Moons: a prior from the primary, a penalty from glare

Satellite searches started from a better localisation prior than any planet search. A moon must lie near its planet. The difficulty was the opposite one: the planet’s glare.

Galileo recognised Jupiter’s four large moons in 1610 from their changing arrangement around a common primary. Christiaan Huygens found Titan in 1655. Giovanni Domenico Cassini found Iapetus in 1671, Rhea in 1672, and Tethys and Dione in 1684 (JPL satellite discovery circumstances). The method was visual: find a faint point near the planet and follow its position relative to the planet until it was clearly travelling with it.

Herschel found Titania and Oberon around Uranus in 1787, and Enceladus and Mimas around Saturn in 1789. His Saturn report describes how the moons’ visibility near the rings and planet changed, and his use of his large reflector (Herschel 1790). William Lassell found Triton in 1846, just weeks after Neptune itself, and Ariel and Umbriel in 1851. Hyperion was found independently by William and George Bond and by Lassell in 1848.

Write a moon’s measured position relative to its planet as

where $\Pi$ is the projection onto the sky and $\mathbf s$ is the satellite’s orbit around the planet. The planet shrinks the search area to a few arcminutes but adds a bright, structured background that varies across that area, and makes the reference point (the planet’s centre) uncertain. A moon can disappear into glare or lie almost along the line of sight to the planet, so a single failure to see it does not refute it.

Asaph Hall found both moons of Mars in August 1877 with the US Naval Observatory’s 26-inch refractor, during a favourable opposition. His original tables begin the outer satellite on 11 August and the inner on 17 August in the astronomical convention of the day, in which the day began at noon. Civil-date retellings often differ by one day for that reason. Hall warned that his positions still needed corrections for differential refraction and for the fact that the apparent centre of a gibbous Mars is not its true centre (Hall 1877). That warning is an early version of a nuisance parameter in modern astrometry.

On 9 September 1892 Edward Emerson Barnard, with the Lick Observatory refractor, found a fifth moon of Jupiter, now Amalthea. The telegram reporting it gives about thirteenth magnitude and prints an early orbital period with a corrective footnote (Lick Observatory 1892). It was the last Jovian moon found by eye.

Photography reached satellites before the century ended. William H. Pickering recognised Saturn’s moon Phoebe on plates that DeLisle Stewart had exposed at Harvard’s Arequipa station in August 1898. Edward C. Pickering announced it in 1899. The report names four plates, A3227, A3228, A3230 and A3233, exposed for 60, 120, 122 and 120 minutes on 16–18 August, and gives the moon’s displacements relative to Saturn as (0, 0), (−2.4, 0.1), (−10.7, 2.4) and (−22.0, 2.9) arcseconds in its own coordinate convention. The object was about photographic magnitude 15.5, and the short arc allowed very different orbital periods (Pickering 1899). Recording, examining, identifying and announcing were done at different times by different people. A stored plate let the recognition happen long after the exposure. Chapter 27 and Chapter 19 continue this line.

Wolf: motion becomes a mark

The trail

Max Wolf at Heidelberg found the asteroid (323) Brucia on 22 December 1891, the first new asteroid discovered photographically. His own report says the initial work used a five-inch Kranz portrait lens, a fast camera lens made for portraits. The six-inch Voigtländer lens at focal ratio 1:5, often shown in illustrations of the discovery, came into use in February 1892 (Wolf 1895).

The method is simple to state. Guide the camera on the stars through a long exposure. Stars come out as points. An asteroid, moving against the stars, comes out as a short trail. The observer scans the developed plate for short lines among the dots.

The difficulty is that other things make short lines. Close double stars, chains of faint stars, defects in the emulsion and guiding errors all imitate trails. Wolf therefore coupled the shape cue to repeated plates. Auguste Charlois at Nice adopted the method quickly and checked candidates with a visual telescope. Wolf’s report and Edward Holden’s contemporary summary describe a complete observing and validation procedure, not a lone lucky plate (Wolf 1895; Holden 1896).

Wolf specifies exposures of at least one and a half hours, preferably two, with the six-inch lens for objects around magnitude 12–13. Under good conditions, the area available for planets brighter than magnitude 12.6 was more than 70 square degrees. For fainter objects the usable field shrank, so this is not a uniform limiting magnitude over a fixed rectangle (Wolf 1895). Seventy square degrees in two hours was an enormous gain over a visual observer checking stars against a chart one at a time.

Why a longer exposure stops helping

A trail is not a free gift. Consider an idealised linear detector, which a photographic plate is not, and a source of flux $F$ (counts per unit time) moving at angular speed $v$ across a background of $b$ counts per unit time per unit area. Let $w$ be the width of the image, set by seeing and optics. In an exposure of length $T$ the trail has length $vT$. If $vT \gg w$, the source deposits $FT$ counts spread over an area of about $w\,vT$, and that area collects background $b\,w\,vT\cdot T$. A detector that sums along the whole trail gets

The exposure time cancels. Once the trail is much longer than the image width, a longer exposure adds source counts and background in step, and the signal-to-noise ratio stops growing. The constant equals the signal-to-noise of an exposure that lasts just long enough for the object to cross one image width, $t_s = w/v$. To see this, a point source in that time gives $Ft_s/\sqrt{b\,w^2 t_s} = F/\sqrt{b\,w\,v}$.

A worked number shows how far into this regime Wolf was. A main-belt asteroid near opposition moves about 36 arcseconds per hour (Chapter 4 derives this). In two hours it trails about 72 arcseconds, far longer than any star image on the plate. Wolf’s long exposures bought him a long, distinctive mark that the eye could recognise against the stars and defects, but not depth. Depth would have needed many shorter exposures added along the motion, which is shift-and-stack (Chapter 17). With $N$ such exposures, each lasting about $t_s$, the signal-to-noise grows as $\sqrt N$ instead of saturating. This scaling is a modern idealisation, not a fitted property of nineteenth-century emulsions, which respond nonlinearly to light.

Two plates, offset in time

Wolf’s best protocol was a staggered pair of exposures (Wolf 1895). Plate A is exposed from hour 0 to hour 2. Plate B, with the second camera, is exposed from hour 1 to hour 3. Both receive two hours. A real asteroid leaves trails of the same direction and length on both plates, but B’s trail is shifted along the motion by half its length. A close pair of stars appears identically on both. A defect appears on one plate only. The whole sequence takes three hours instead of four for two consecutive exposures.

In modern terms, a trail of unknown position, direction and length on a background-subtracted linear image $I - B$ with pixel variance $V$ is found with the matched statistic

where $T_\gamma$ is a trail template with parameters $\gamma$. Searching over $P$ positions, $Q$ trail shapes and $\ell$ samples per trail costs $O(PQ\ell)$ directly. Line-sum transforms and hierarchical schemes reuse partial sums (Chapter 11). The second plate adds a different kind of evidence: agreement in direction and displacement between independent images, which no single-image statistic can provide. Wolf reported no formal threshold or false-alarm probability. His two-plate test was his false-alarm control.

The yield

Holden’s summary of Wolf’s report gives the objects recorded on the six-inch plates from February 1892 to April 1895, split into those already known and those new (Wolf 1895; Holden 1896).

In 1892 one object in three on Wolf’s plates was new. By 1895 it was about one in seven. That decline is descriptive. It does not prove the belt was becoming complete at his limit, and it is not an efficiency per observing hour. But the figure in the asteroid-cohort chart above shows what photography did overall: after 1891 the discovery rate climbed again, and in the twentieth century photographic surveys would carry the count from hundreds to thousands (Chapter 4). Wolf also made the trade-off between field and image scale explicit in comparing portrait lenses with reflectors (Wolf 1897). A wide, fast lens covers sky. A long reflector gives small images and reaches fainter stars. The tension between the two is the design question of every survey since.

Landmarks

What the nineteenth century could not measure

Historical noise models are incomplete. None of these discoveries comes with a distribution of rejected candidates, a measured completeness, a plate-by-plate calibration, or the variance between observers. Their limiting magnitudes are operational descriptions. Giving Herschel a sigma threshold or Wolf a false-alarm rate would invent a result.

Detection and identity fail independently. A body can be recorded and not recognised (the pre-1781 Uranus observations, Challis’s Neptune). It can be recognised as moving and then lost (Ceres before Gauss). It can be linked to the wrong earlier observation. Phoebe’s arc allowed very different orbits. A useful historical database would record first exposure, first recognition, announcement and secure orbit as separate events.

The reference can dominate the sensitivity. A complete chart turns a hard visual search into checking a short list of missing entries. An incomplete or wrong chart produces false candidates, and an aggressive match against it can erase a real mover. Old star positions, proper motions, clock corrections, field distortion and printing errors are all part of the signal chain.

Short arcs are not identifiable. Many trajectories fit a short sequence of angles, and a single local covariance can hide disconnected solutions. More arithmetic cannot replace a baseline that was never observed. A recent paper on measuring asteroid distances from one observatory in one night, using the observer’s own motion, is a modern continuation of this point and discusses the astrometric systematics that limit it (Fernandes et al. 2025).

Perturbations are conditional evidence. Neptune’s good pointing did not validate its predicted orbit, and Mercury’s residual did not validate Vulcan. Unknown masses, initial conditions, reference-frame errors and the law of gravity have to be tested together. Fits made after the answer is known are particularly vulnerable to choosing data and weights that make the answer look inevitable. Even recent reconstructions contain slips: one 2024 numerical study misdates Galileo’s sightings of Neptune and writes Bode where Bouvard is meant (Rodríguez-Moris & Docobo 2024), and Krajnović’s “one minute of arc” contradicts his own table.

Where the room is

The arithmetic in this chapter is tiny. Le Verrier’s final system had 33 equations and nine unknowns. A laptop solves it in microseconds, and no GPU is needed to reproduce it. If there is room for a new method here, it is not in speeding up these calculations. It is in the information-use lessons they teach, which carry over directly to modern problems with millions of hypotheses.

Better prediction shrinks the search area. Gauss’s orbit moved the search by nine degrees. Le Verrier’s turned a survey into a single night. In both cases the gain came from inference, with the telescope unchanged. In a modern survey the same principle appears as propagating every credible short-arc orbit and choosing the follow-up observation that best separates them, with weather and glare included. The objective is expected reduction in recovery area. This can be a large operational gain without collecting one more photon (Chapter 13).

Separate linear parameters from nonlinear search. In a modern Neptune-like search, factor the nuisance basis (the perturbed planet’s own orbit and the known masses) once, project both residuals and templates into its complement, and solve for each template’s mass analytically. This avoids refitting the same nuisance directions for every trial perturber, and it shows which part of each template is observable at all. This is an algorithmic assessment, not a measured speedup. It matters for the Planet Nine problem in Chapter 8, where the trial bank is large.

Batch only when the ensemble is large. A grid of candidate perturbers, many nuisance models, or many calibration simulations can share background ephemerides and run as batched propagations and small matrix operations on a GPU. A historical-scale problem offers nothing to batch. Any claim of a tenfold gain needs a stated baseline and identical acceptance, sensitivity and calibration.

Throughput from observation to comparison is a detection limit. Challis had Neptune twice. The bottleneck was the comparison, not the telescope. The modern equivalent is a pipeline whose latency or review capacity decides whether a candidate is followed up before it is lost.

Integrate along the motion, do not trail. Wolf’s trails show the limit of single long exposures: once the trail is longer than the image, signal-to-noise stops growing. Splitting the time into short exposures and adding them along trial motions removes that ceiling and leads directly to shift-and-stack (Chapter 17) and its tree and transform accelerations (Chapter 21, Chapter 39). Digitised plate archives could also be searched this way for precoveries, although emulsion nonlinearity and damage need empirical calibration.

Know which limits are hard. Information lost to glare, saturated plates, bad seeing or unobserved intervals cannot be recovered by faster summation. Neither can a unique perturber be inferred when its effect lies almost entirely in the fitted nuisance subspace. Those are physical and informational limits. Chart completeness, scheduling, registration, indexing and repeated calculation are engineering limits, and those are the ones a new method can move. The general case is in Chapter 42.

References

1. Adams, J. C. (1846). An Explanation of the observed Irregularities in the Motion of Uranus, on the Hypothesis of Disturbance caused by a more distant Planet. Monthly Notices of the Royal Astronomical Society 7, 149–152. ADS scan
2. Airy, G. B. (1846). Account of some circumstances historically connected with the discovery of the Planet exterior to Uranus. Monthly Notices of the Royal Astronomical Society 7, 121–144. ADS scan
3. Bovaird, T., Lineweaver, C. H. (2013). Exoplanet Predictions Based on the Generalised Titius–Bode Relation. arXiv preprint. arXiv:1304.3341
4. Capron, J. R. (1877). Report of Examination of the Sun's Disk at Guildown, Guildford, on 21st, 22nd, and 23rd March 1877, for the suspected planet Vulcan. Monthly Notices of the Royal Astronomical Society 37, 348–349. doi:10.1093/mnras/37.6.348
5. Fernandes, M. B., Scolnic, D., Peterson, E. et al. (2025). Measuring the Distances to Asteroids from One Observatory in One Night with Upcoming All-Sky Telescopes. arXiv preprint. arXiv:2502.07881
6. Galle, J. G. (1877). Ein Nachtrag zu den … Berichten über die erste Auffindung des Planeten Neptun. Astronomische Nachrichten 89, 349–351. ADS scan
7. Gauss, C. F. (1809; English translation by C. H. Davis, 1857). Theoria motus corporum coelestium / Theory of the Motion of the Heavenly Bodies Moving about the Sun in Conic Sections. Little, Brown. archive.org
8. Hall, A. (1877). Observations of the Satellites of Mars, Made with the 26-inch Refractor of the U. S. Naval Observatory. Astronomische Nachrichten 91, 11–14. ADS scan
9. Herschel, W. (1781). Account of a Comet. Philosophical Transactions of the Royal Society of London 71, 492–501. archive.org
10. Herschel, W. (1790). Account of the Discovery of a Sixth and Seventh Satellite of the Planet Saturn. Philosophical Transactions of the Royal Society of London 80. doi:10.1098/rstl.1790.0001
11. Hind, J. R. (1876). The Intra-Mercurial Planet or Planets. Nature 14, 469–470. doi:10.1038/014469a0
12. Holden, E. S. (1896). The Photography of Planetoids, by Professor Max Wolf. Publications of the Astronomical Society of the Pacific 8, 23–26. archive.org
13. JPL Solar System Dynamics (accessed 2026). Planetary Satellite Discovery Circumstances. Online register. ssd.jpl.nasa.gov
14. Krajnović, D. (2016). The Contrivance of Neptune. Astronomy & Geophysics 57(5), 5.28–5.34. arXiv:1610.06424
15. Le Verrier, U. J. (1846). Recherches sur les mouvements d'Uranus (Fortsetzung). Astronomische Nachrichten 25, 65–80. ADS scan
16. Lick Observatory (1892). Discovery of a Fifth Satellite of Jupiter by Professor Barnard. Publications of the Astronomical Society of the Pacific 4, 199. ADS scan
17. Minor Planet Center (accessed 2026). Discovery Circumstances: Numbered Minor Planets (1)–(5000). Online catalogue. minorplanetcenter.net
18. Nature (1879). The Intra-Mercurial Planet Question. Nature 20, 597–598. doi:10.1038/020597a0
19. Newcomb, S. (1895). The Elements of the Four Inner Planets and the Fundamental Constants of Astronomy. Astronomical Papers of the American Ephemeris 6. archive.org
20. O'Connor, J. J., Robertson, E. F. (accessed 2026). Giuseppe Piazzi. MacTutor History of Mathematics. mathshistory.st-andrews.ac.uk
21. Orekit contributors (accessed 2026). IodGauss: Gauss angles-only initial orbit determination. Orekit 13.1 documentation. orekit.org
22. Perry, S. J. (1877). The Planet Vulcan. Monthly Notices of the Royal Astronomical Society 37, 347–348. doi:10.1093/mnras/37.6.347
23. Piazzi, G. (1801). Risultati delle osservazioni della nuova stella scoperta il dì 1. Gennaio all'Osservatorio Reale di Palermo. Palermo. archive.org
24. Piazzi, G. (1802). Della scoperta del nuovo pianeta Cerere Ferdinandea. Stamperia Reale, Palermo. archive.org
25. Pickering, E. C. (1899). A New Satellite of Saturn. Astrophysical Journal 9, 274–276. ADS scan
26. Remchin, J., Schrimpf, A. (2019). Research on asteroids of Christian Ludwig Gerling and his students in the nineteenth century. arXiv preprint. arXiv:1903.12544
27. Rodríguez-Moris, G., Docobo, J. A. (2024). The Discovery of Neptune Revisited. arXiv preprint. arXiv:2405.06310
28. Wolf, M. (1895). Die Photographie der Planetoiden. Astronomische Nachrichten 139, 97–110. ADS scan
29. Wolf, M. (1897). Reflector and Portrait Lens in Celestial Photography. Nature 55, 582–584. doi:10.1038/055582a0