Хаос задачи трех тел преодолен? Математики научились находить точные орбиты по точкам выравнивания

Взаимодействие двух космических тел — например, одиночной звезды и обращающейся вокруг нее планеты — описывается классическими законами Кеплера. Траектория планеты в такой изолированной системе неизменна во времени: это правильный эллипс, круг или гипербола. Параметры движения рассчитываются точно, а будущее системы можно предсказать на миллионы лет вперед.

Однако в реальном космосе изолированных систем не существует. Стоит добавить к гравитационному взаимодействию третье тело, как математическая простота исчезает. Общая задача трех тел не имеет решения в виде универсальных формул. Гравитационное поле каждого участника непрерывно меняет траектории остальных, порождая сложную динамику. Малейшее отклонение в начальной скорости любого из компонентов со временем может привести к непредсказуемым последствиям: вытягиванию орбит, столкновению или гравитационному выбросу одного из участников в межзвездное пространство.

При этом в нашей планетной системе наблюдаются зоны удивительного долголетия. В окрестностях точек Лагранжа системы Солнце — Юпитер (в зонах гравитационного равновесия, которые опережают планету-гигант на шестьдесят градусов по ее орбите и отстают от нее на такое же расстояние) обращаются тысячи так называемых троянских астероидов. Гравитация массивного Юпитера и Солнца уравновешивает инерцию этих малых тел, удерживая их в относительно стабильных коридорах на протяжении миллиардов лет.

При этом в аналогичных зонах системы Солнце — Сатурн крупных популяций астероидов практически нет. Причина этой разницы кроется именно в гравитации третьего участника: Юпитер своей массой дестабилизирует орбитальные области Сатурна, постепенно выталкивая из них любые малые объекты.

Чтобы рассчитывать движение космических аппаратов, поведение спутников или орбиты астероидов в таких системах, астрономам требуется анализировать так называемую ограниченную задачу четырех тел. В ней участвуют три массивных объекта (например, звезда и две планеты) и одно малое тело (зонд или астероид), масса которого ничтожна и не способна повлиять на гигантов.

Однако корректно рассчитать движение малого тела невозможно без точного понимания того, как движутся сами три массивных объекта.

Пределы стандартных расчетов

На протяжении десятилетий небесная механика использовала приближенную схему, известную как бициркулярная модель. Чтобы обойти неразрешимость задачи трех тел, исследователи прибегали к искусственному допущению: двум главным телам принудительно задавали строго круговое движение вокруг общего центра масс, а третье массивное тело заставляли обращаться вокруг них по правильной внешней окружности.

Подобный расчет называют некогерентным. Его главный недостаток заключается в нарушении фундаментальной физики: заданные круговые орбиты не удовлетворяют законам ньютоновской гравитации. В реальности массивные тела непрерывно деформируют движение друг друга. Если планеты массивны или сближаются слишком тесно, погрешность некогерентных моделей накапливается лавинообразно, приводя к ошибочным выводам о стабильности всей системы.

Чтобы расчеты совпадали с физической реальностью, астродинамике требуются когерентные модели. Это означает, что траектории всех массивных участников должны в каждый момент времени строго подчиняться уравнениям тяготения.

Команда математиков из Барселонского университета и Университета Сантьяго-де-Компостела разработала численный метод, позволяющий системно находить и непрерывно отслеживать целые семейства таких точных решений. В центре их внимания оказались так называемые относительные периодические орбиты.

Геометрия вращающегося мира

В классической механике периодической называют орбиту, на которой все тела через строго определенный отрезок времени — период обращения — возвращаются в исходные точки пространства с теми же самыми скоростями. Но в задаче трех тел такие абсолютно периодические орбиты изолированы. Они представляют собой редкие математические исключения, которые разрушаются при минимальном изменении параметров системы.

Относительные периодические орбиты устроены более гибко. С точки зрения неподвижного наблюдателя траектории трех тел не замыкаются в простые кольца: тела выписывают сложные узоры, которые кажутся хаотическими. Однако во вращающейся системе координат это движение становится строго замкнутым.

Через определенный промежуток времени вся система трех тел полностью воспроизводит свое взаимное расположение и внутренние скорости, но сама плоскость этой конфигурации за время пути поворачивается на постоянный угол.

Принцип относительной периодичностиТела не обязаны возвращаться в исходные точки неподвижного пространства. Достаточно, чтобы через один цикл взаимные расстояния между ними и их относительные скорости в точности повторились, а вся система лишь сместилась на определенный угол поворота.

Если этот угол поворота составляет рациональную долю от полной окружности (например, ровно половину, треть или четверть круга), то спустя несколько полных циклов орбита замкнется и в неподвижном пространстве, превратившись в абсолютно периодическую. Если же угол поворота иррационален, траектория никогда не замкнется, но форма взаимного движения тел останется стабильной неограниченно долго.

Главная ценность таких орбит в том, что они образуют непрерывные семейства: плавно изменяя всего один физический параметр (например, расстояние между планетами), исследователи могут проследить непрерывную эволюцию устойчивых орбит от простых форм до крайне сложных.

Метод сизигий: упрощение уравнений через прямую линию

Основная сложность задачи трех тел заключается в количестве неизвестных параметров. Чтобы описать движение трех материальных точек на плоскости, требуется двенадцать дифференциальных уравнений: по две координаты и по две составляющие скорости для каждого тела. Расчет устойчивости в двенадцатимерном пространстве параметров требует колоссальных вычислительных мощностей.

Математики решили эту проблему с помощью геометрического явления, известного в астрономии как сизигия.

Сизигией называют момент времени, когда три или более небесных тела выстраиваются строго вдоль одной прямой линии. В Солнечной системе примером сизигии служат затмения, когда Солнце, Луна и Земля оказываются на одной геометрической оси.

В математической модели исследователи использовали сизигию как строгий фильтр:

Условие расчетной сизигии

1. В начальный момент времени все три тела располагаются строго на горизонтальной координатной оси.
2. Векторы скоростей всех трех тел направлены строго вертикально (перпендикулярно линии выравнивания).
3. Общий центр масс фиксируется в нулевой точке координат, что исключает из расчетов движение центра тяжести системы.
4. Единица длины выбирается так, чтобы расстояние от центра до первого тела всегда равнялось единице.

Благодаря законам сохранения импульса, энергии и масштаба двенадцатимерная система уравнений сжимается всего до трех независимых переменных: начального расстояния до второго тела и двух вертикальных скоростей.

Компьютерный алгоритм запускает движение тел из этой прямой линии и рассчитывает их гравитационное взаимодействие до тех пор, пока система снова не выстроится вдоль одной прямой. Если во время этого второго выравнивания расстояния между телами и величины их скоростей в точности совпали с исходными, а направления скоростей вновь перпендикулярны линии (пусть сама эта линия и повернулась на определенный угол), то относительная периодическая орбита успешно зафиксирована.

После этого метод численного продолжения позволяет делать плавные шаги: изменяя расстояние между телами, алгоритм автоматически подбирает новые скорости, выстраивая непрерывную кривую решений.

Отсечение математического шума: проверка устойчивости

Мало просто найти замкнутую орбиту — необходимо выяснить, сможет ли реальная система существовать на ней, не распадаясь от малейшего толчка.

В теории дифференциальных уравнений стабильность периодического движения оценивают с помощью матрицы возмущений (матрицы монодромии). Она показывает, как случайные отклонения от заданной траектории меняются за время одного периода. Если эти отклонения со временем затухают или остаются в фиксированных границах, орбита устойчива. Если хотя бы одно возмущение начинает расти по экспоненте, орбита физически нестабильна.

Однако в задаче об относительных периодических орбитах стандартная матрица монодромии содержит системную погрешность. Фундаментальные физические законы — сохранение полной механической энергии, сохранение момента вращения, симметрия времени и симметрия поворота в пространстве — порождают в уравнениях паразитные коэффициенты, в точности равные единице. Еще два паразитных коэффициента возникают из-за самого поворота координатной сетки между выравниваниями.

Эти нейтральные математические величины не отражают реальную физическую прочность орбиты. Однако при расчетах на компьютере они накапливают ошибки округления и искажают результаты.

Авторы применили строгий метод матричной проекции (дефляцию на основе сингулярного разложения). Алгоритм рассчитывает направления, отвечающие за физические законы сохранения, и проецирует матрицу возмущений в очищенное четырехмерное пространство параметров.

В результате удаляются все фиктивные коэффициенты. Если оставшиеся четыре показателя лежат в пределах допустимых комплексных значений, орбита признается строго устойчивой. Если показатели выходят в опасную вещественную зону, траектория считается нестабильной.

Три сценария: планеты, луны и двойные звезды

Используя разработанный вычислительный аппарат, авторы проанализировали три классические астрономические структуры, получив точные карты их поведения.

1. Планетные системы (семейства Пуанкаре)

В этой схеме рассматривается тяжелая звезда и две планеты, обращающиеся вокруг нее в одном направлении. Одна планета по массе сопоставима с Юпитером (одна тысячная массы звезды), а масса второй варьируется от земной (одна миллионная) до юпитерианской.

Пока планеты находятся далеко друг от друга, их орбиты почти идеально круговые, а вся система находится в устойчивом равновесии. Однако при постепенном сокращении дистанции между ними алгоритм зафиксировал критические точки:

• Орбитальный резонанс три к одному: внутреннее тело совершает ровно три оборота за то время, пока внешнее делает один. В этой конфигурации гравитационные толчки планет синхронизируются. Если внешняя планета достаточно массивна, круговые орбиты теряют устойчивость.
• Точка поворота: при дальнейшем сближении система достигает барьера, где взаимное расстояние между телами на прямой линии минимально, а общий момент вращения падает практически до нуля. Тела оказываются на пороге критического сближения.
• Переход к вытянутым орбитам: классические теории прошлого века утверждали, что за этой границей регулярные орбиты перестают существовать. Однако новый метод показал, что семейство решений продолжается. Орбиты претерпевают кардинальную трансформацию: почти круговые траектории вытягиваются в удлиненные эллипсы. Более того, после прохождения зоны дестабилизации сильно вытянутые эллиптические орбиты вновь восстанавливают устойчивость: планеты могут двигаться по вытянутым путям миллиарды лет, избегая опасных сближений.

2. Спутниковые системы (семейства Хилла)

Второй сценарий воспроизводит систему из центральной звезды, массивной планеты и обращающегося вокруг нее легкого спутника.

В астрономии граница гравитационного влияния планеты называется сферой Хилла. Внутри этой области притяжение планеты преобладает над приливным притяжением центральной звезды. Для выбранных в исследовании параметров расчетный радиус сферы Хилла составил около тридцати двух тысячных долей от расстояния между планетой и звездой.

Алгоритм начал расчет со спутника, расположенного в непосредственной близости от планеты, а затем стал постепенно увеличивать его дистанцию:

Эволюция спутниковой орбиты

1. Вблизи планеты спутник стабильно удерживается ее притяжением, а энергия их взаимного движения отрицательна.
2. При удалении спутника на расстояние около сорока пяти тысячных долей единицы энергия взаимной связи меняет знак на положительный. Спутник покидает сферу Хилла.
3. Вместо падения на звезду или хаотического выброса система переходит в новое устойчивое состояние: спутник становится самостоятельной планетой. Он переходит на околозвездную орбиту, которая регулярно пересекает траекторию родительской планеты, но сохраняет периодический ритм движения.

Этот расчет математически воспроизводит феномен временного захвата малых тел. Именно по такой траектории околоземные астероиды (например, объект 2020 CD3) на несколько месяцев или лет переходят на орбиту вокруг Земли, становясь временными спутниками, а затем снова сходят с нее и возвращаются на стабильный путь вокруг Солнца.

3. Двойные звездные системы

В третьей конфигурации исследователи смоделировали поведение планеты малой массы в системе из двух одинаковых звезд, обращающихся вокруг общего центра тяжести.

Динамика третьего тела полностью зависела от его начальной скорости при пересечении линии сизигии:

• Циркумбинарные орбиты: при высокой начальной скорости планета обращается вокруг обеих звезд сразу по широкой замкнутой внешней дуге. Вся система сохраняет стабильность на протяжении сотен расчетных оборотов.
• Ретроградные орбиты: при уменьшении начальной скорости система претерпевает инверсию. Скорость планеты меняет знак на противоположный: она начинает двигаться в направлении, обратном вращению звездных компонентов (по часовой стрелке против общего вращения системы против часовой стрелки). При сближении со звездами планета формирует сложные многоконтурные замкнутые траектории с центральным сужением, оставаясь гравитационно привязанной только к одной из звезд.

Фундамент для проекта TROY и космической навигации

Построенные математиками орбиты не просто описывают теоретические траектории на плоскости. В физике нелинейных систем доказано фундаментальное правило: вокруг любой устойчивой периодической траектории в пространстве параметров всегда формируются объемные коридоры стабильности (инвариантные торы).

Это означает, что физические системы, параметры которых слегка отклоняются от расчетного идеала — например, если планеты имеют небольшие взаимные наклоны орбит, а их массы определены с погрешностью, — все равно будут следовать вдоль этих траекторий, сохраняя стабильность в течение астрономических промежутков времени.

Это открытие имеет прямое прикладное значение для проекта TROY — международной наблюдательной программы по поиску троянских экзопланет. Астрофизики ведут поиск землеподобных каменистых миров, которые скрываются в точках Лагранжа массивных газовых гигантов у других звезд.

Обнаружить сравнительно небольшую каменистую планету на расстоянии в сотни световых лет напрямую с помощью телескопов крайне трудно. Однако если в звездной системе присутствует крупный газовый гигант, стабильность его лагранжевых зон можно заранее рассчитать математически.

Прежние некогерентные расчеты нередко признавали такие зоны нестабильными там, где они в реальности сохраняли устойчивость, либо требовали месяцев работы суперкомпьютеров для перебора миллионов случайных траекторий.

Новый метод последовательных выравниваний дает наблюдателям точный координатный каркас. Относительные периодические орбиты трех тел задают абсолютно непротиворечивое гравитационное поле. Добавляя в эту готовую систему пробную частицу с нулевой массой, исследователи получают возможность быстро и безошибочно рассчитывать границы зон стабильности, находить безопасные орбиты для космических аппаратов и с высокой точностью определять сектора звездных систем, где могут существовать пригодные для жизни миры.

Источник: arXiv