1. Введение: Проклятие третьей звезды в часовом механизме Вселенной

В 1687 году Исаак Ньютон опубликовал свои монументальные «Математические начала натуральной философии» (Principia), навсегда изменив представление человечества о космосе. Доказав, что падение яблока и движение Луны подчиняются одному и тому же закону всемирного тяготения, Ньютон полностью и точно решил задачу двух тел (Two-Body Problem): траектории двух материальных точек, притягивающихся по закону обратных квадратов, вечно остаются коническими сечениями (эллипсами, параболами или гиперболами), выражаемыми через элементарные функции, что стало фундаментом детерминизма «демона Лапласа».

Однако стоит добавить в систему Солнце–Земля всего одно третье небесное тело — Луну, перейдя от \(N=2\) к \(N=3\), то есть к задаче трёх тел (Three-Body Problem), — как аналитическая магия рушится. Почему система из двух звёзд решается методами первокурсника, а три звезды ставили в тупик величайших математиков более трёх столетий?

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

2. Полная интегрируемость задачи двух тел и размерностная стена задачи трёх тел

2.1 Гамильтонова формулировка и теорема Лиувилля — Арнольда

Для \(N\) точечных масс \(m_i\) с радиус-векторами \(\mathbf{r}_i \in \mathbb{R}^3\) уравнения движения Ньютона имеют вид:

В \(6N\)-мерном фазовом пространстве с каноническими координатами \(\mathbf{q} = (\mathbf{r}_1, \dots, \mathbf{r}_N)\) и импульсами \(\mathbf{p} = (m_1\dot{\mathbf{r}}_1, \dots, m_N\dot{\mathbf{r}}_N)\) система задаётся функцией Гамильтона \(H(\mathbf{q}, \mathbf{p})\):

Согласно теореме Лиувилля — Арнольда (1855 / 1963), автономная гамильтонова система с \(n\) степенями свободы (\(2n\)-мерное фазовое пространство) вполне интегрируема в квадратурах тогда и только тогда, когда существуют \(n\) функционально независимых первых интегралов \(F_1 (=H), F_2, \dots, F_n\), находящихся в инволюции (скобки Пуассона попарно равны нулю: \(\{F_i, F_j\} = 0\)). Тогда компактные поверхности уровня диффеоморфны \(n\)-мерным инвариантным торам \(\mathbb{T}^n\), а движение по ним условно-периодично.

По теореме Нётер из симметрий группы Галилея для любой системы \(N\) тел следуют ровно 10 классических алгебраических интегралов движения:

1. Полный импульс \(\mathbf{P}\) (3) и равномерное движение центра масс \(\mathbf{G}_{\text{CM}} = M\mathbf{R}_{\text{CM}} - \mathbf{P}t\) (3): Безусловно понижают размерность фазового пространства на 6 единиц (\(6N \to 6N - 6\)).
2. Полный момент импульса \(\mathbf{L} = \sum \mathbf{r}_i \times \mathbf{p}_i\) (3): В силу коммутационных соотношений \(\{L_x, L_y\} = L_z \ne 0\) из них можно выделить ровно 2 интеграла в инволюции (\(\|\mathbf{L}\|^2\) и \(L_z\)), что с помощью исключения узлов Якоби (1842) понижает размерность фазового пространства на 4.
3. Полная механическая энергия \(H = E\) (1).

2.2 Баланс степеней свободы: почему \(N=2\) интегрируется, а \(N=3\) — нет

  flowchart TD
    subgraph TwoBody["Задача двух тел (N = 2): 12D фазовое пространство (6 степеней свободы)"]
        TB1["Исходное фазовое пространство: 6 × 2 = 12 измерений"] --> TB2["Отделение центра масс (6 интегралов)<br/>12 - 6 = 6D относительное пространство (n_rel = 3)"]
        TB2 --> TB3["Есть 3 коммутирующих интеграла: H_rel, |L|^2, L_z<br/>Условие Лиувилля — Арнольда выполнено (3 = 3)!"]
        TB3 --> TB4["Вектор Лапласа — Рунге — Ленца A -> Суперинтегрируемость<br/>Точное решение: конические сечения Кеплера"]
    end

    subgraph ThreeBody["Задача трёх тел (N = 3): 18D фазовое пространство (9 степеней свободы)"]
        THB1["Исходное фазовое пространство: 6 × 3 = 18 измерений"] --> THB2["Координаты Якоби / Центр масс (6 интегралов)<br/>18 - 6 = 12D относительное пространство (n_rel = 6)"]
        THB2 --> THB3["Для n_rel = 6 степеней свободы есть лишь 3 интеграла в инволюции:<br/>H_rel, |L|^2, L_z (Дефицит 3 интегралов!)"]
        THB3 --> THB4["Остаётся неприводимое 7D–8D пространство!<br/>Теоремы Брунса (1887) и Пуанкаре (1889): Неинтегрируемость"]
    end

• При \(N=2\): После отделения центра масс относительное движение (\(\mathbf{r} = \mathbf{r}_1 - \mathbf{r}_2\), приведённая масса \(\mu_{\text{red}} = \frac{m_1 m_2}{m_1+m_2}\)) имеет \(n_{\text{rel}} = 3\) степени свободы и ровно 3 коммутирующих интеграла (\(H_{\text{rel}}, \|\mathbf{L}\|^2, L_z\)). Более того, скрытая симметрия \(SO(4)\) потенциала \(1/r\) сохраняет вектор Лапласа — Рунге — Ленца (LRL): \(\mathbf{A} = \mathbf{p} \times \mathbf{L} - \mu_{\text{red}} G m_1 m_2 \frac{\mathbf{r}}{\|\mathbf{r}\|}\) делая систему максимально суперинтегрируемой с замкнутыми эллиптическими орбитами \(r(\theta) = \frac{p}{1 + e\cos(\theta - \varpi)}\).
• При \(N=3\): Размерность фазового пространства возрастает до \(6 \times 3 = 18\). После исключения центра масс остаётся 12 измерений (\(n_{\text{rel}} = 6\) степеней свободы), но по-прежнему лишь 3 интеграла в инволюции! Даже после исключения узлов Якоби и фиксации уровня энергии остаются 8 симплектических измерений (или 7-мерное многообразие энергии).
  • Теорема Брунса (1887) и Теорема Пуанкаре о неинтегрируемости (1889) доказали отсутствие любых других независимых алгебраических или аналитических первых интегралов.
  • Ловушка ряда Зундмана (1912): Хотя финский математик Карл Зундман построил сходящийся при всех \(t\) степенной ряд по \(t^{1/3}\) (при \(\mathbf{L} \ne \mathbf{0}\)), Давид Белорицкий (1930) показал, что для вычисления всего одного оборота требуется просуммировать не менее \(10^{8\,000\,000}\) членов ряда, что делает его физически бесполезным!

3. Круговая ограниченная задача трёх тел (CR3BP) и пять точек Лагранжа (\(L_1\)–\(L_5\))

3.1 Уравнения во вращающейся синодической системе и интеграл Якоби

В CR3BP два массивных тела \(m_1 \ge m_2\) обращаются по круговым орбитам (\(a=1, \omega=1\)) вокруг общего барицентра, а третье тело пренебрежимо малой массы (\(m_3 \to 0\)) движется в их поле. С безразмерным параметром масс \(\mu = \frac{m_2}{m_1 + m_2} \in (0, 1/2]\) основные тела неподвижны во вращающейся системе координат \((x,y,z)\) в точках \((-\mu, 0, 0)\) и \((1-\mu, 0, 0)\). Вводя эффективный потенциал \(\Omega(x, y, z)\):

получаем уравнения движения с учётом силы Кориолиса (\(-2\boldsymbol{\omega}\times\mathbf{v}\)):

Умножая на \(\dot{x}, \dot{y}, \dot{z}\) и интегрируя, находим единственный закон сохранения CR3BP — интеграл Якоби (постоянную Якоби \(C_J\)):

Так как \(v_{\text{rot}}^2 = 2\Omega - C_J \ge 0\), движение ограничено областью Хилла \(2\Omega(x,y,z) \ge C_J\), ограниченной поверхностями нулевой скорости (Zero-Velocity Surfaces) \(2\Omega = C_J\).

3.2 Точки равновесия \(L_1\)–\(L_5\): Телескоп «Джеймс Уэбб» и критерий Рауса

Условие \(\nabla \Omega = \mathbf{0}\) в плоскости \(z=0\) даёт пять точек либрации:

1. Три коллинеарные точки Эйлера (\(L_1, L_2, L_3\), 1767): Лежат на прямой \(y=0\) и определяются уравнением пятой степени \(x - \frac{(1-\mu)(x+\mu)}{|x+\mu|^3} - \frac{\mu(x-1+\mu)}{|x-1+\mu|^3} = 0\). При \(\mu \ll 1\) точки \(L_1\) и \(L_2\) отстоят от \(m_2\) на радиус сферы Хилла \(r_{\text{Hill}} \approx a(\mu/3)^{1/3}\) (около 1,5 млн км от Земли в системе Солнце–Земля).
2. Две треугольные точки Лагранжа (\(L_4, L_5\), 1772): Записывая \(\Omega(r_1, r_2) = \frac{1-\mu}{2}(r_1^2 + \frac{2}{r_1}) + \frac{\mu}{2}(r_2^2 + \frac{2}{r_2}) - \frac{1}{2}\mu(1-\mu)\), из \(\frac{\partial \Omega}{\partial r_1} = \frac{\partial \Omega}{\partial r_2} = 0\) мгновенно получаем \(r_1 = r_2 = 1\) — вершины двух равносторонних треугольников \(L_{4,5} = \left(\frac{1}{2}-\mu, \pm\frac{\sqrt{3}}{2}, 0\right)\).

Линейный анализ устойчивости (\(\lambda^4 + (4 - \Omega_{xx} - \Omega_{yy})\lambda^2 + (\Omega_{xx}\Omega_{yy} - \Omega_{xy}^2) = 0\)) показывает:

• Седловая неустойчивость \(L_1, L_2, L_3\) и гало-орбита телескопа «Джеймс Уэбб» (JWST): Так как \(\Omega_{xx}\Omega_{yy} < 0\), корни имеют вид \(\pm\sigma, \pm i\omega_{xy}, \pm i\omega_z\) (\(\text{Saddle} \times \text{Center} \times \text{Center}\)), поэтому \(L_1, L_2, L_3\) всегда неустойчивы. Однако на 4-мерном центральном многообразии существуют трёхмерные периодические гало-орбиты (Halo Orbits). Космический телескоп «Джеймс Уэбб» (JWST) обращается вокруг точки \(L_2\) по гало-орбите: он никогда не заходит в тень Земли, закрывается единым пятислойным тепловым экраном сразу от Солнца, Земли и Луны (охлаждаясь до \(<40\text{ K}\)) и тратит на удержание орбиты всего \(\Delta v \approx 2\text{–}4\text{ м/с}\) в год!
• Гироскопическая стабилизация \(L_4, L_5\) силой Кориолиса и критическое отношение масс Рауса: Хотя \(L_4\) и \(L_5\) являются максимумами (вершинами) потенциала \(\Omega\), сила Кориолиса (\(-2\boldsymbol{\omega}\times\mathbf{v}\)) заворачивает скатывающуюся частицу в замкнутую орбиту вокруг вершины! Характеристическое уравнение \(\lambda^4 + \lambda^2 + \frac{27}{4}\mu(1 - \mu) = 0\) даёт чисто мнимые корни (линейную устойчивость) тогда и только тогда, когда дискриминант положителен (критерий Рауса, 1875): \(1 - 27\mu(1 - \mu) > 0 \iff \mu < \mu_{\text{crit}} = \frac{1 - \sqrt{23/27}}{2} \approx 0.0385209 \quad \left(\frac{m_1}{m_2} > 24.96\right)\) В системе Солнце–Юпитер (\(\mu \approx 0.000954 \ll 0.03852\)) эта устойчивость удерживает более 12 000 троянских астероидов в окрестностях \(L_4\) и \(L_5\)!

4. Премия короля Оскара II и ошибка Анри Пуанкаре (1889), породившая теорию хаоса

В 1889 году Анри Пуанкаре выиграл международный конкурс, учреждённый королём Швеции и Норвегии Оскаром II. Во время печати мемуара в Acta Mathematica молодой редактор Ларс Эдвард Фрагмен заметил пробел в доказательстве. Перепроверив выкладки, Пуанкаре с ужасом обнаружил, что устойчивое многообразие \(W^s(p_0)\) и неустойчивое многообразие \(W^u(p_0)\) гиперболической периодической орбиты не сливаются гладко, а пересекаются трансверсально (под ненулевым углом)!

  flowchart TD
    A["Исходное допущение Пуанкаре (Рукопись 1888 г.)<br/>W^s и W^u гладко сливаются в замкнутую сепаратрису"] --> B["Вопрос редактора Фрагмена (Июль 1889 г.)<br/>Пуанкаре перепроверяет ряды теории возмущений"]
    B --> C["Открытие трансверсальной гомоклинической точки q_0!<br/>W^s и W^u пересекаются под ненулевым углом"]
    C --> D["Инвариантность + Сохранение площади (Лиувилль)<br/>Одно пересечение влечёт бесконечное число пересечений P^k(q_0)"]
    D --> E["Гомоклиническое сплетение (Пуанкаре / Подкова Смейла)<br/>Экспоненциальная чувствительность -> Рождение теории хаоса (1890)!"]

В силу инвариантности относительно отображения Пуанкаре \(P\) одно трансверсальное пересечение в точке \(q_0\) порождает бесконечное множество гомоклинических точек \(P^k(q_0)\) (\(k \in \mathbb{Z}\)), а сохранение фазовой площади вынуждает петли многообразий бесконечно растягиваться и складываться, образуя гомоклиническое сплетение (Homoclinic Tangle / подкова Смейла, 1967). Пуанкаре за свой счёт заплатил 3 581 крону (больше самой премии в 2 500 крон!), чтобы изъять и уничтожить отпечатанный тираж, и в 1890 году опубликовал исправленный труд, положивший начало теории динамического хаоса.

5. Современные чудеса: Орбита-восьмёрка и КАМ-теория в Солнечной системе

5.1 Хореография «восьмёрки» (Мур 1993; Шенсине и Монтгомери 2000)

Спустя более двух веков после Эйлера и Лагранжа Кристофер Мур (1993) численно открыл, а Ален Шенсине и Ричард Монтгомери (2000) строго доказали в Annals of Mathematics с помощью принципа наименьшего действия (вариационного исчисления) существование периодической орбиты в форме восьмёрки (Figure-Eight Orbit) для трёх равных масс (\(m_1=m_2=m_3=1\)) с нулевым моментом импульса (\(\mathbf{L}=\mathbf{0}\)):

Минимизируя функционал действия \(\mathcal{A}[\mathbf{q}] = \int_0^T (T - V)\,dt\) с \(D_6\)-симметрией на сфере форм (Shape Sphere \(S^2\)), они доказали, что минимум действия строго меньше действия любой траектории со столкновением (\(\mathcal{A}_{\text{eight}} < \mathcal{A}_{\text{coll}}\)). Впоследствии Шуваков и Дмитрашинович (2013) и Ли и Ляо (2017) открыли более 2 000 новых семейств периодических орбит.

5.2 Теорема КАМ, диффузия Арнольда и 5-миллионолетний горизонт Ляпунова (Ласкар)

В близкой к интегрируемой гамильтоновой системе \(H(\mathbf{I}, \boldsymbol{\theta}) = H_0(\mathbf{I}) + \epsilon H_1(\mathbf{I}, \boldsymbol{\theta})\) резонансы порождают проблему малых знаменателей \(\mathbf{k}\cdot\boldsymbol{\omega}(\mathbf{I}) \approx 0\). Теорема Колмогорова — Арнольда — Мозера (КАМ, 1954–1963) доказала, что при \(\epsilon \ll 1\) большинство инвариантных торов сохраняется, если их частоты удовлетворяют диофантову условию:

тогда как при \(n \ge 3\) степенях свободы разрушенные резонансы сливаются в единую сеть, допускающую медленную диффузию Арнольда.

Наконец, многомиллиарднолетние численные интегрирования Жака Ласкара (1989 / 2009) показали, что внутренняя Солнечная система хаотична с временем Ляпунова всего \(\tau_{\text{L}} \approx 5\text{ миллионов лет}\): начальная погрешность положения Земли в 15 метров вырастает до \(1\text{ а.е.}\) за 100 миллионов лет, делая невозможным точный прогноз орбитальной фазы за этим пределом, а вековой резонанс между Меркурием и Юпитером (\(g_1 \approx g_5\)) создаёт \(\sim 1\%\)-ную вероятность катастрофического роста эксцентриситета Меркурия (\(e > 0.6\)) и столкновения внутренних планет в ближайшие 5 миллиардов лет.