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 миллиардов лет.