Спотлайт #69: Орбітальна механіка — задача трьох тіл, точки Лагранжа та порожнини Роша

I. Обмежена задача трьох тіл в обертовій системі відліку

Повна задача трьох тіл — три точкові маси, що взаємодіють через ньютонівську гравітацію — загалом не має розв’язку в замкненому вигляді. Анрі Пуанкаре показав ще 1890 року, що рух, як правило, хаотичний: близькі траєкторії розходяться експоненційно, що робить довгострокове передбачення принципово неможливим. Але один особливий випадок виявляється точно розв’язним і надзвичайно корисним: обмежена задача трьох тіл (R3BP), у якій два масивні тіла (первинні тіла, маси M і m, причому M ≫ m) обертаються навколо спільного центра мас по колових орбітах, а третє тіло знехтувано малої маси рухається в їхньому гравітаційному полі.

Ключове спрощення — перейти до системи відліку, що обертається разом із тілами (синодичної системи відліку), яка обертається з тією самою кутовою швидкістю ω, що й первинні тіла. У цій системі первинні тіла нерухомі. Безмасова тестова частинка відчуває не лише гравітаційне притягання обох первинних тіл, а й дві фіктивні сили, що виникають через обертання: відцентрову силу (що штовхає назовні від осі обертання) та силу Коріоліса (що відхиляє рухомі частинки перпендикулярно до їхньої швидкості).

Рівняння руху тестової частинки в цій системі відліку мають вигляд:

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

де ефективний потенціал U* поєднує гравітаційні та відцентрові члени:

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

r1 = відстань від тестової частинки до первинного тіла M
r2 = відстань від тестової частинки до первинного тіла m
omega = sqrt(G*(M + m) / a^3)    (кутова орбітальна швидкість, 3-й закон Кеплера)
a    = відстань між первинними тілами

Члени Коріоліса (-2ωy' і 2ωx') у лівій частині зв’язують рівняння для x і y, роблячи систему неконсервативною у звичайному сенсі. Однак існує одна збережена величина: інтеграл Якобі, який часто позначають C і який відіграє роль повної енергії в обертовій системі відліку.

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

де v^2 = x'^2 + y'^2 — швидкість в обертовій системі відліку.

Частинка в спокої (v = 0) задовольняє C = 2*U*(x, y).
Це визначає криві нульової швидкості: контури U*, що дорівнюють C/2.
Частинка не може перетнути криву нульової швидкості (це вимагало б v^2 < 0).

Криві нульової швидкості (також звані кривими Гілла) — це межі областей, які частинка може досягти для заданого значення C. Високі значення C щільно обмежують частинку поблизу одного з первинних тіл; коли C зменшується, заборонені області скорочуються, і відкриваються нові — врешті-решт дозволяючи частинці вільно циркулювати навколо обох первинних тіл.

II. П’ять точок рівноваги Лагранжа

Точки рівноваги ефективного потенціалу — де градієнт U* дорівнює нулю, і, отже, нерухома частинка відчуває нульову результівну силу — це точки Лагранжа. Їх рівно п’ять, позначених від L1 до L5. Жозеф-Луї Лагранж знайшов L4 і L5 аналітично 1772 року; колінеарні точки L1–L3 були знайдені раніше Ейлером.

Колінеарні точки: L1, L2, L3

Ці три точки лежать на прямій, що з’єднує два первинні тіла. Усі три є сідловими точками U*: стабільні вздовж лінії первинних тіл, нестабільні перпендикулярно до неї. Оскільки вони перебувають у сідлових точках, будь-яке мале збурення змушує частинку експоненційно віддалятися — вони є нестійкими за Ляпуновим. Космічні апарати, розташовані в L1 чи L2, потребують періодичних маневрів утримання позиції кожні кілька тижнів, щоб залишатися на місці.

L1: між первинними тілами, ближче до меншої маси m.
    Наближене положення (відношення мас mu = m/(M+m) << 1):
    x_L1 ~ a * (1 - (mu/3)^(1/3))    від більшого первинного тіла

L2: з протилежного від M боку m.
    x_L2 ~ a * (1 + (mu/3)^(1/3))

L3: з протилежного від m боку M (майже навпроти m).
    x_L3 ~ -a * (1 + 5*mu/12)

    Усі три: нестабільні (сідлові точки U*).
    Характерний час нестабільності ~ орбітальний період / (3*mu)^(1/3)

L1 — це брама між двома гравітаційними ямами: речовина, що протікає через L1 (наприклад, у подвійній зоряній системі), переходить зі сфери впливу одного тіла до іншого. Саме так відбувається перенесення маси в тісних подвійних зорях і катаклізмічних змінних. L2, розташована з протилежного боку меншого первинного тіла, — це місце, де перебуває космічний телескоп Джеймса Вебба — детальніше про це нижче.

Трикутні точки: L4 і L5

L4 і L5 розташовані у вершинах рівносторонніх трикутників, утворених з двома первинними тілами. Вони є максимумами U* (локальними енергетичними пагорбами), а не мінімумами, що наївно наводить на думку про нестабільність — але сила Коріоліса рятує ситуацію. Коли відношення мас задовольняє M/m > 24.96 (приблизно), відхилення Коріоліса від будь-якого малого збурення створює замкнену епіциклічну орбіту навколо L4 чи L5, а не траєкторію втечі. Обидві системи Сонце–Земля і Сонце–Юпітер із запасом задовольняють цей критерій.

L4: на 60 градусів попереду меншого первинного тіла m на його орбіті
L5: на 60 градусів позаду m

Критерій стабільності (умова Рауса):
    27 * mu * (1 - mu) < 1
    тобто mu < 0.0385 (частка маси меншого первинного тіла)

    Сонце-Юпітер: mu ~ 0.00095   --> стабільно
    Земля-Місяць: mu ~ 0.012     --> стабільно
    Плутон-Харон: mu ~ 0.11      --> НЕСТАБІЛЬНО (троянців не спостерігається)

Стабільність L4 і L5 примітна: це не мінімуми потенціалу, але сила Коріоліса створює те, що фактично є відновлювальною силою для малих зміщень. Частинки, захоплені поблизу L4/L5, виконують повільні епіциклічні («пуголовкові») орбіти навколо точки рівноваги з періодом приблизно в 1/sqrt(27*mu) разів більшим за орбітальний період первинних тіл.

III. Порожнини Роша, перенесення маси й криві нульової швидкості на практиці

Порожнина Роша кожного первинного тіла — це область простору, у межах якої речовина гравітаційно пов’язана саме з цим тілом, а не з іншим. Вона визначається як область, обмежена кривою нульової швидкості, що проходить через внутрішню точку Лагранжа L1. Каплеподібна межа, що проходить через L1, називається поверхнею Роша.

У подвійних зоряних системах поняття порожнини Роша є центральним для розуміння перенесення маси. Якщо одна із зір розширюється (наприклад, на стадії червоного гіганта), поки не заповнить свою порожнину Роша, речовина в точці L1 відчуває однакове гравітаційне притягання від обох зір і починає перетікати до супутника. Цей процес живить одні з найенергетичніших явищ у Всесвіті: нові зорі, рентгенівські подвійні системи та наднові типу Ia.

Радіус порожнини Роша (наближення Еґлтона, 1983):
    R_L / a = 0.49 * q^(2/3) / (0.6 * q^(2/3) + ln(1 + q^(1/3)))

    q = m / M  (відношення мас, вторинного тіла до первинного)
    a = орбітальна відстань

Приклад: зоря сонячного типу + білий карлик, a = 1 сонячний радіус
    q ~ 0.5 --> R_L ~ 0.38 * a ~ 0.38 сонячного радіуса
    Коли зоря-донор сходить з головної послідовності й розбухає,
    щойно R_star > R_L, починається перенесення маси, часто нестримне.

Форма кривих нульової швидкості драматично змінюється зі зменшенням константи Якобі C. Для дуже великих C два окремі замкнені овали оточують кожне первинне тіло (порожнини Роша), із забороненою областю між ними й назовні. Коли C падає до значення в L1, внутрішня заборонена область стискається, і дві порожнини Роша стикаються в L1. За ще нижчих значень C крива нульової швидкості проходить через L2 і L3 по черзі, відкриваючи шляхи для виходу речовини із системи повністю — механізм, причетний до формування подвійних систем зі спільною оболонкою.

IV. Реальні приклади: JWST у точці L2 і троянські астероїди

Космічний телескоп Джеймса Вебба в точці L2 системи Сонце–Земля

Точка L2 системи Сонце–Земля розташована приблизно за 1.5 мільйона кілометрів від Землі, у напрямку прямо протилежному Сонцю. Оскільки L2 обертається разом із Землею, апарат там бачить Сонце, Землю й Місяць у межах 50-градусного конуса з одного боку — уся інша півсфера неба перебуває в постійній тіні. Це робить L2 ідеальним місцем для інфрачервоної астрономії: сонцезахисний екран телескопа може одночасно блокувати тепло від Сонця, Землі й Місяця за допомогою єдиної п’ятишарової фольгованої конструкції.

JWST не перебуває точно в L2 (точна точка не дає жодної переваги для утримання позиції й періодично затінювалася б Землею). Натомість він обертається навколо L2 по великій гало-орбіті з амплітудою приблизно 500 000 км, нахиленій так, щоб Земля й Місяць завжди залишалися поза зоною виключення сонцезахисного екрана. Орбітальний період становить близько 6 місяців. Маневри утримання позиції на рівні приблизно 2–4 м/с на рік необхідні, щоб протидіяти нестабільності колінеарної точки Лагранжа.

Троянські астероїди в точках L4 і L5 системи Юпітер–Сонце

Система Сонце–Юпітер задовольняє критерій стабільності Рауса зі значним запасом, і результат вражаючий: понад 12 000 відомих астероїдів (юпітеріанські троянці) лібрують навколо L4 («грецький табір», за 60° попереду Юпітера) і L5 («троянський табір», за 60° позаду). Вони захоплені в «пуголовкові» орбіти навколо точок рівноваги, лібруючи з амплітудами до 30 градусів і періодами в десятки років.

Місія NASA Lucy, запущена в жовтні 2021 року, — перший космічний апарат для дослідження троянських роїв. Вона пролетить повз вісім троянських астероїдів у період між 2027 і 2033 роками, використовуючи гравітаційні маневри біля Землі тричі, щоб досягти необхідної траєкторії. Вважається, що троянці — це первісні залишки ранньої Сонячної системи, можливо захоплені під час періоду міграції гігантських планет, описаного моделлю Ніцци — що робить їх капсулами часу, які зберігають умови 4.5 мільярда років тому.

Власні точки L4 і L5 Землі не порожні: астероїд 2010 TK7 лібрує навколо земної L4, і кілька супутних об’єктів знайдено в L5. Марс також має кілька троянців. Навіть Сатурн і Уран мають невеликі популяції, хоча в Нептуна їх особливо багато — понад 30 відомих, а моделі свідчать, що справжня популяція може конкурувати з юпітеріанською за чисельністю.

Спробуйте самі

Підсумкова думка

Обмежена задача трьох тіл — одна з найпродуктивніших «нерозв’язних» задач у фізиці. Не маючи змоги записати загальні розв’язки, математики натомість картографували структуру простору розв’язків: топологію кривих нульової швидкості, стабільність точок рівноваги, існування періодичних орбіт та інваріантних торів. Саме це структурне знання інженери використовують, коли відправляють JWST до L2 чи Lucy до троянців. Урок ширший за орбітальну механіку: розуміння того, чому система поводиться саме так, як поводиться — які збережені величини існують, які області досяжні, які точки рівноваги стабільні — часто є потужнішим за будь-який конкретний розв’язок. Рівняння — це карта; фізика — це територія.