Da Ceres forsvandt
Nytårsdag 1801 opdagede Giuseppe Piazzi i Palermo et svagt, vandrende lyspunkt — den første asteroide, Ceres. Han nåede at følge den i seks uger, før den gled ind i Solens skær. Spørgsmålet var nu: hvor dukker den op igen? Datidens metoder krævede lange observationsrækker; Piazzis korte bue på 9° syntes håbløst utilstrækkelig.
Gauss angreb problemet på en helt ny måde. Han indså, at tre retningsbestemmelser — tre par af rektascension og deklination med tilhørende tidspunkter — i princippet indeholder netop de seks tal, en banes seks elementer kræver. I december 1801 fandt Franz von Zach Ceres mindre end en halv grad fra Gauss' forudsagte position. Metoden, han senere udgav i Theoria Motus (1809), bruges i sin kerne stadig i dag, når nyopdagede asteroider og kometer skal have deres første bane.
Jupiter er naturligvis ingen nyopdagelse — men netop derfor er den et perfekt prøvelegeme: Vi kender facit og kan efterprøve, hvor godt metoden virker på en bue, der kun udgør halvanden procent af det 11,86 år lange omløb.
Tre aftener ved teleskopet
De tre billeder er taget med en Unistellar eVscope2 fra observationsstedet ved Marshøj (56° N, 11° Ø). Billedernes tidsstempler giver observationstidspunkterne med sekundpræcision — og tiden er metodens halvdel af målingen. Den anden halvdel er retningen til planeten: rektascension \( \alpha \) og deklination \( \delta \). I planet-tilstand indeholder eVscope-billederne ingen baggrundsstjerner at pladeløse på, så retningerne nedenfor er de astrometriske positioner, der hører til de fotograferede tidspunkter (svarende til almanak- eller plate‑solve-astrometri, korrigeret for lysets løbetid). Det er præcis de tre talpar, Gauss' maskineri skal fodres med.
| nr. | Dato & kl. (UT) | Juliansk dato | Rektascension α | Deklination δ |
|---|---|---|---|---|
| 1 | 14. jan 2026 · 18:33:47 | 2 461 055,27346 | 07ʰ 23ᵐ 03,3ˢ | +22° 20′ 17″ |
| 2 | 14. feb 2026 · 17:11:51 | 2 461 086,21656 | 07ʰ 08ᵐ 12,5ˢ | +22° 50′ 22″ |
| 3 | 22. mar 2026 · 17:47:22 | 2 461 122,24123 | 07ʰ 05ᵐ 01,3ˢ | +22° 57′ 42″ |
Bemærk Jupiters retrograde bevægelse: rektascensionen aftager hen over de tre datoer, fordi Jorden — der var i opposition til Jupiter den 10. januar — overhaler indenom. Netop denne sløjfebevægelse, set fra en platform der selv er i fart, er det, metoden skal vikle ud.
Sigtelinjer fra en bevægelig platform
Hver observation fastlægger en sigtelinje: en halvlinje fra Jorden i retningen \( \hat{\boldsymbol\rho} \), givet ved
hvor \( \mathbf R \) er Jordens heliocentriske stedvektor (kendt fra almanakken), \( \mathbf r \) planetens ukendte heliocentriske stedvektor og \( \rho \) den ukendte afstand fra Jorden til planeten langs sigtelinjen. Tre observationer giver tre ligninger af denne type — ni skalare ligninger med ni ubekendte (\( \mathbf r_1,\mathbf r_2,\mathbf r_3 \) bundet sammen af tyngdeloven og de tre \( \rho_i \)).
Gauss' maskineri
1. Tidsintervaller og Lagranges f- og g-funktioner
Med \( \tau_1=t_1-t_2 \), \( \tau_3=t_3-t_2 \) og \( \tau=\tau_3-\tau_1 \) (her \( \tau_1=-30{,}943 \) døgn, \( \tau_3=+36{,}025 \) døgn, \( \tau=66{,}968 \) døgn) udnytter Gauss, at de tre stedvektorer ligger i samme baneplan, så \( \mathbf r_1 \) og \( \mathbf r_3 \) kan skrives som kombinationer af \( \mathbf r_2 \) og \( \mathbf v_2 \):
hvor rækkeudviklingerne i første omgang afkortes som vist. Heri gemmer tyngdeloven sig: \( \mu=GM_\odot=k^2 \) med Gauss' egen konstant \( k=0{,}01720209895 \) (enheder: au og døgn).
2. Skalartripelprodukter og ottendegradsligningen
Indsættes f- og g-udtrykkene i sigtelinjeligningerne og elimineres \( \rho_1,\rho_3 \) med skalartripelprodukter af enhedsvektorerne (Curtis' \( D \)-determinanter), står man tilbage med én ligning i én ubekendt — planetens heliocentriske afstand ved midterobservationen:
For vores Jupiter-data bliver koefficienterne \( a=-27{,}5796 \), \( b=+38{,}4971 \), \( c=-13{,}6332 \) (au-potenser), og polynomiet har tre positive reelle rødder:
Flertydigheden er ægte og velkendt: tre sigtelinjer alene kan ikke altid skelne en nær, hurtig bane fra en fjern, langsom. Her vælger fysikken for os — et legeme på en næsten jordnær bane ville flytte sig mange grader pr. døgn, ikke Jupiters magelige bueminutter.
3. Fra afstand til tilstandsvektor
Med \( r_2=5{,}2259 \) au følger de tre skrå afstande \( \rho_1,\rho_2,\rho_3 \) direkte, dermed \( \mathbf r_1,\mathbf r_2,\mathbf r_3 \), og hastigheden
Herefter itereres: med \( (\mathbf r_2,\mathbf v_2) \) beregnes f og g eksakt (universelle variable i stedet for de afkortede rækker), de skrå afstande genberegnes — og lysets løbetid \( \rho_i/c \) (33–41 minutter!) trækkes fra observationstiderne. Processen konvergerer her på 6 iterationer til
4. Baneelementerne
Tilstandsvektoren oversættes til keplerelementer via impulsmoment \( \mathbf h=\mathbf r\times\mathbf v \), excentricitetsvektoren \( \mathbf e=(\mathbf v\times\mathbf h)/\mu-\hat{\mathbf r} \) og energien \( \varepsilon=v^2/2-\mu/r=-\mu/2a \). Resultatet står i næste afsnit.
Jupiters bane — bestemt fra Djursland
| Element | Gauss' metode | Reference | Afvigelse |
|---|---|---|---|
| Halve storakse \(a\) | 5,2073 au | 5,2029 au | +0,08 % |
| Excentricitet \(e\) | 0,0483 | 0,0484 | −0,2 % |
| Inklination \(i\) | 1,3038° | 1,3044° | −0,0006° |
| Knudelængde \(\Omega\) | 100,53° | 100,47° | +0,05° |
| Perihellængde \(\varpi\) | 15,80° | 14,73° | +1,07° |
| Omløbstid \(P\) | 11,883 år | 11,862 år | +0,18 % |
Med en bue på blot 1,5 % af omløbet rammes \( a \), \( e \), \( i \) og \( \Omega \) på promilleniveau. Kun perihellængden \( \varpi \) afviger en grads penge — ganske som forventet: med \( e\approx 0{,}05 \) er banen så nær cirkulær, at perihelets retning er den dårligst bestemte størrelse. Det var i øvrigt det samme, Gauss oplevede med Ceres.
Hvorfor virker det så godt?
To gaver fra geometrien: For det første blev første billede taget fire døgn efter opposition, hvor parallaksen fra Jordens egen bevægelse er størst — sigtelinjerne »sakser« kraftigt og fastlægger afstanden skarpt. For det andet er 67 døgn nok til, at Jupiters egen krumning af banen kan mærkes i tredjeordensleddene af f- og g-rækkerne. Gauss' metode er netop bygget til denne mellemvægt: buer for korte til klassisk banetilpasning, men lange nok til at krumningen taler.
Forbehold, der hører med
Beregningen er geocentrisk: fra 56° N flytter topocentrisk parallakse Jupiter under 2 buesekunder, hvilket er uden betydning her. Lysets løbetid er derimod ikke til at ignorere (op til 41 minutter) og indgår i iterationen. Observationstiderne er aflæst af billedernes tidsstempler (CET = UT+1); positionerne er astrometriske J2000-størrelser hørende til netop disse tidspunkter. Vil man gøre øvelsen helt »rå«, fotograferer man planeten i vidfelt-tilstand med baggrundsstjerner og pladeløser selv — metoden nedstrøms er nøjagtig den samme.
Fra Palermo 1801 til Marshøj 2026
Det forbløffende ved Gauss' metode er økonomien: seks målte tal (tre vinkelpar) plus tre aflæsninger af et ur bestemmer samtlige seks baneelementer — fordi Newtons tyngdelov leverer al den manglende information. Hver gang en amatørastronom i dag indrapporterer astrometri af en nyfunden komet til Minor Planet Center, er det denne maskine (i moderniseret form), der spinner i den anden ende.
Og dermed er ringen sluttet: et teleskop på en mark på Djursland, tre klare aftener hen over en vinter, og matematik fra 1801 — det er alt, hvad der skal til for at veje Solsystemets største planet ind på sin plads.