Spotlight #69: Mechanika orbitalna — problem trzech ciał, punkty Lagrange'a i płaty Roche'a

I. Ograniczony problem trzech ciał we współobracającym się układzie odniesienia

Pełny problem trzech ciał — trzy punktowe masy oddziałujące grawitacją newtonowską — generalnie nie ma rozwiązania w postaci zamkniętej. Henri Poincaré pokazał w 1890 roku, że ruch jest z reguły chaotyczny: bliskie trajektorie rozbiegają się wykładniczo, co czyni długoterminową predykcję zasadniczo niemożliwą. Jeden szczególny przypadek okazuje się jednak dokładnie rozwiązywalny i niezwykle użyteczny: ograniczony problem trzech ciał (R3BP), w którym dwa masywne ciała (ciała pierwotne, masy M i m, przy czym M ≫ m) krążą wokół wspólnego środka masy po orbitach kołowych, a trzecie ciało o pomijalnej masie porusza się w ich polu grawitacyjnym.

Kluczowe uproszczenie polega na przejściu do współobracającego się (synodycznego) układu odniesienia, który obraca się z tą samą prędkością kątową ω co ciała pierwotne. W tym układzie ciała pierwotne są nieruchome. Bezmasowa cząstka testowa doświadcza nie tylko przyciągania grawitacyjnego obu ciał pierwotnych, ale też dwóch sił pozornych wynikających z obrotu: siły odśrodkowej (odpychającej od osi obrotu) i siły Coriolisa (odchylającej poruszające się cząstki prostopadle do ich prędkości).

Równania ruchu cząstki testowej w tym układzie mają postać:

x'' - 2*omega*y' = dU*/dx
y'' + 2*omega*x' = dU*/dy

gdzie efektywny potencjał U* łączy człony grawitacyjne i odśrodkowe:

U*(x, y) = -G*M/r1 - G*m/r2 - 0.5*omega^2*(x^2 + y^2)

r1 = odległość cząstki testowej od ciała pierwotnego M
r2 = odległość cząstki testowej od ciała pierwotnego m
omega = sqrt(G*(M + m) / a^3)    (prędkość kątowa orbity, III prawo Keplera)
a    = odległość między ciałami pierwotnymi

Człony Coriolisa (-2ωy' i 2ωx') po lewej stronie sprzęgają równania dla x i y, czyniąc układ niezachowawczym w zwykłym sensie. Istnieje jednak jedna wielkość zachowana: całka Jacobiego, często oznaczana C, która pełni rolę całkowitej energii w obracającym się układzie.

C = 2*U*(x, y) - v^2

gdzie v^2 = x'^2 + y'^2 to prędkość w obracającym się układzie odniesienia.

Cząstka w spoczynku (v = 0) spełnia C = 2*U*(x, y).
To definiuje krzywe zerowej prędkości: kontury U* równe C/2.
Cząstka nie może przekroczyć krzywej zerowej prędkości (wymagałoby to v^2 < 0).

Krzywe zerowej prędkości (nazywane też krzywymi Hilla) to granice obszarów, które cząstka może osiągnąć dla danej wartości C. Wysokie wartości C ciasno ograniczają cząstkę w pobliżu jednego z ciał pierwotnych; gdy C maleje, obszary zabronione kurczą się i otwierają się nowe, aż w końcu cząstka może swobodnie krążyć wokół obu ciał pierwotnych.

II. Pięć punktów równowagi Lagrange'a

Punkty równowagi potencjału efektywnego — gdzie gradient U* znika, a więc nieruchoma cząstka doświadcza zerowej siły wypadkowej — to punkty Lagrange'a. Jest ich dokładnie pięć, oznaczonych od L1 do L5. Joseph-Louis Lagrange znalazł L4 i L5 analitycznie w 1772 roku; punkty współliniowe L1–L3 zostały odkryte wcześniej przez Eulera.

Punkty współliniowe: L1, L2, L3

Te trzy punkty leżą na linii łączącej dwa ciała pierwotne. Wszystkie trzy są punktami siodłowymi U*: stabilne wzdłuż linii ciał pierwotnych, niestabilne prostopadle do niej. Ponieważ znajdują się w punktach siodłowych, każde małe zaburzenie powoduje wykładnicze oddalanie się cząstki — są niestabilne w sensie Lapunowa. Statki kosmiczne umieszczone w L1 lub L2 wymagają okresowych manewrów utrzymania pozycji co kilka tygodni, by pozostać na miejscu.

L1: między ciałami pierwotnymi, bliżej mniejszej masy m.
    Przybliżone położenie (stosunek mas mu = m/(M+m) << 1):
    x_L1 ~ a * (1 - (mu/3)^(1/3))    od większego ciała pierwotnego

L2: po dalszej stronie m od M.
    x_L2 ~ a * (1 + (mu/3)^(1/3))

L3: po dalszej stronie M od m (niemal naprzeciw m).
    x_L3 ~ -a * (1 + 5*mu/12)

    Wszystkie trzy: niestabilne (punkty siodłowe U*).
    Charakterystyczna skala czasowa niestabilności ~ okres orbitalny / (3*mu)^(1/3)

L1 to brama między dwiema studniami grawitacyjnymi: materia przepływająca przez L1 (na przykład w układzie gwiazd podwójnych) przechodzi ze sfery wpływu jednego ciała do drugiego. Dokładnie tak działa transfer masy w ciasnych układach gwiazd podwójnych i zmiennych kataklizmicznych. L2, położony po przeciwnej stronie mniejszego ciała pierwotnego, to miejsce, w którym przebywa Kosmiczny Teleskop Jamesa Webba — więcej o tym poniżej.

Punkty trójkątne: L4 i L5

L4 i L5 znajdują się w wierzchołkach trójkątów równobocznych utworzonych z dwoma ciałami pierwotnymi. Są maksimami U* (lokalnymi wzniesieniami energii), a nie minimami, co naiwnie sugerowałoby niestabilność — ale ratuje je siła Coriolisa. Gdy stosunek mas spełnia M/m > 24,96 (w przybliżeniu), odchylenie Coriolisa dowolnego małego zaburzenia tworzy zamkniętą orbitę epicykliczną wokół L4 lub L5 zamiast trajektorii ucieczki. Zarówno układ Słońce–Ziemia, jak i Słońce–Jowisz spełniają ten kryterium z dużym zapasem.

L4: 60 stopni przed mniejszym ciałem pierwotnym m na jego orbicie
L5: 60 stopni za m

Kryterium stabilności (warunek Routha):
    27 * mu * (1 - mu) < 1
    czyli mu < 0,0385 (ułamek masy mniejszego ciała pierwotnego)

    Słońce-Jowisz: mu ~ 0,00095   --> stabilne
    Ziemia-Księżyc: mu ~ 0,012    --> stabilne
    Pluton-Charon: mu ~ 0,11      --> NIESTABILNE (nie obserwuje się trojańczyków)

Stabilność L4 i L5 jest niezwykła: nie są to minima potencjału, ale siła Coriolisa tworzy coś, co skutecznie działa jak siła przywracająca dla małych przemieszczeń. Cząstki uwięzione w pobliżu L4/L5 wykonują powolne epicykliczne (typu „kijanka") orbity wokół punktu równowagi z okresem w przybliżeniu 1/sqrt(27*mu) razy dłuższym niż okres orbitalny ciał pierwotnych.

III. Płaty Roche'a, transfer masy i krzywe zerowej prędkości w praktyce

Płat Roche'a każdego ciała pierwotnego to obszar przestrzeni, w którym materia jest grawitacyjnie związana z tym ciałem, a nie z drugim. Definiuje się go jako obszar zamknięty przez krzywą zerowej prędkości przechodzącą przez wewnętrzny punkt Lagrange'a L1. Łzowaty kontur przechodzący przez L1 nazywa się powierzchnią Roche'a.

W układach gwiazd podwójnych pojęcie płata Roche'a jest kluczowe dla zrozumienia transferu masy. Jeśli jedna gwiazda rozszerza się (jak w fazie czerwonego olbrzyma), aż wypełni swój płat Roche'a, materia w punkcie L1 odczuwa równe przyciąganie grawitacyjne od obu gwiazd i zaczyna przepływać do towarzysza. Ten proces napędza jedne z najbardziej energetycznych zjawisk we Wszechświecie: nowe, podwójne układy rentgenowskie i supernowe typu Ia.

Promień płata Roche'a (przybliżenie Eggletona, 1983):
    R_L / a = 0,49 * q^(2/3) / (0,6 * q^(2/3) + ln(1 + q^(1/3)))

    q = m / M  (stosunek mas, wtórnego do pierwotnego)
    a = odległość orbitalna

Przykład: gwiazda podobna do Słońca + biały karzeł, a = 1 promień słoneczny
    q ~ 0,5 --> R_L ~ 0,38 * a ~ 0,38 promienia słonecznego
    Gdy gwiazda-dawca schodzi z ciągu głównego i puchnie,
    gdy tylko R_gwiazdy > R_L, zaczyna się transfer masy, często lawinowy.

Kształt krzywych zerowej prędkości zmienia się dramatycznie wraz ze spadkiem stałej Jacobiego C. Dla bardzo dużych C dwie oddzielne zamknięte owalne krzywe otaczają każde ciało pierwotne (płaty Roche'a), z obszarem zabronionym między nimi i na zewnątrz. Gdy C spada do wartości w L1, wewnętrzny obszar zabroniony zaciska się, a dwa płaty Roche'a stykają się w L1. Dla jeszcze niższych C krzywa zerowej prędkości przechodzi kolejno przez L2 i L3, otwierając drogi ucieczki materii z układu całkowicie — mechanizm powiązany z formowaniem się układów podwójnych ze wspólną otoczką.

IV. Przykłady z rzeczywistości: JWST w L2 i asteroidy trojańskie

Kosmiczny Teleskop Jamesa Webba w punkcie L2 układu Słońce–Ziemia

Punkt L2 układu Słońce–Ziemia znajduje się około 1,5 miliona kilometrów od Ziemi, dokładnie po przeciwnej stronie od Słońca. Ponieważ L2 obraca się razem z Ziemią, statek kosmiczny widzi tam Słońce, Ziemię i Księżyc w obrębie 50-stopniowego stożka po jednej stronie — cała pozostała półkula nieba jest stale w cieniu. Czyni to L2 idealną lokalizacją dla astronomii w podczerwieni: osłona przeciwsłoneczna teleskopu może jednocześnie blokować ciepło od Słońca, Ziemi i Księżyca za pomocą jednej pięciowarstwowej struktury foliowej.

JWST nie znajduje się dokładnie w L2 (dokładny punkt nie daje żadnej przewagi w utrzymaniu pozycji i byłby okresowo zaćmiewany przez Ziemię). Zamiast tego krąży wokół L2 po dużej orbicie halo o amplitudzie około 500 000 km, nachylonej tak, by Ziemia i Księżyc zawsze pozostawały poza strefą wykluczenia osłony przeciwsłonecznej. Okres orbitalny wynosi około 6 miesięcy. Manewry utrzymania pozycji rzędu 2–4 m/s rocznie są konieczne, by przeciwdziałać niestabilności współliniowego punktu Lagrange'a.

Asteroidy trojańskie w punktach L4 i L5 układu Jowisz–Słońce

Układ Słońce–Jowisz spełnia kryterium stabilności Routha z dużym zapasem, a wynik jest spektakularny: ponad 12 000 znanych asteroid (trojańczycy Jowisza) libruje wokół L4 (obóz grecki, 60° przed Jowiszem) i L5 (obóz trojański, 60° za nim). Są uwięzione na orbitach typu „kijanka" wokół punktów równowagi, librując z amplitudami do 30 stopni i okresami liczonymi w dekadach.

Misja NASA Lucy, wystrzelona w październiku 2021 roku, to pierwszy statek kosmiczny badający roje trojańczyków. Przeleci obok ośmiu asteroid trojańskich w latach 2027–2033, wykorzystując trzykrotnie asysty grawitacyjne Ziemi, by osiągnąć niezbędną trajektorię. Uważa się, że trojańczycy to pierwotne pozostałości wczesnego Układu Słonecznego, prawdopodobnie przechwycone podczas okresu migracji wielkich planet opisanego przez model z Nicei — co czyni je kapsułami czasu z warunków sprzed 4,5 miliarda lat.

Własne punkty L4 i L5 Ziemi nie są puste: asteroida 2010 TK7 libruje wokół ziemskiego L4, a garstkę towarzyszy koorbitalnych znaleziono w L5. Mars ma również kilku trojańczyków. Nawet Saturn i Uran mają niewielkie populacje, choć u Neptuna jest ich szczególnie dużo — ponad 30 znanych, a modele sugerują, że prawdziwa populacja może dorównywać liczebnie tej u Jowisza.

Wypróbuj sam

Myśl na zakończenie

Ograniczony problem trzech ciał to jeden z najbardziej produktywnych „nierozwiązywalnych" problemów w fizyce. Nie mogąc zapisać ogólnych rozwiązań, matematycy zamiast tego skartografowali strukturę przestrzeni rozwiązań: topologię krzywych zerowej prędkości, stabilność punktów równowagi, istnienie orbit okresowych i niezmienniczych torusów. Ta strukturalna wiedza jest dokładnie tym, co inżynierowie wykorzystują, wysyłając JWST do L2 czy Lucy do trojańczyków. Lekcja wykracza poza mechanikę orbitalną: zrozumienie, dlaczego układ zachowuje się tak, a nie inaczej — jakie wielkości są zachowane, jakie obszary są dostępne, które punkty równowagi są stabilne — jest często potężniejsze niż jakiekolwiek konkretne rozwiązanie. Równania są mapą; fizyka jest terytorium.