Примем теперь во внимание взаимное притяжение небесных тел. Масса самой большой планеты, Юпитера, в тысячу с небольшим раз меньше Солнечной. Примерно во столько же раз ускорение каждой планеты, вызванное притяжением других планет, меньше ускорения к Солнцу. Дифференциальное уравнение движения можно записать в форме
w=F0+µF1 (10)
Здесь индексом 0 отмечено основное ускорение, индексом 1 — вызванное притяжением планет друг к другу возмущающее ускорение; малый параметр µ~0,001. Уравнений типа (10) надо написать несколько, по числу планет. Движение при µ=0 нам известно. При истинном малом значении µ траектория чуть-чуть отклоняется от невозмущенной. Допустимо считать, что орбита по-прежнему является эллипсом, но его элементы (большая полуось, эксцентриситет и т.д.) медленно меняются со временем со скоростью порядка µ. Этот прием мы уже рассматривали на примере ИСЗ.
Фундаментальный вопрос: накапливаются ли возмущения со временем или колеблются около некоторого среднего значения? В первом случае мы говорим о вековых возмущениях; за время T0/µ орбиты менялись бы до неузнаваемости. Здесь Т0 — характерный период, 10 лет для Солнечной системы, примерно год Юпитера. Критическое время T0/µ равно всего десяти тысячам лет, совсем немного в истории Земли. К нашему счастью, при отсутствии резонансов возмущения большой полуоси, эксцентриситета и наклона не накапливаются, эти важнейшие для жизни на Земле элементы лишь колеблются в узких пределах.
Надо сказать, что к этому результату математики и астрономы шли три столетия. Очень уж трудно доказать эту теорему, ведь уравнения (10) настолько сложны, что до сих пор не найдено в аналитическом виде их общего решения, пригодного на космогонических временах порядка T0/µ3. Ньютон полагал, что возмущения накапливаются. В образах того времени Великий Часовщик создал часы не абсолютного совершенства, нуждающиеся в ремонте один раз в несколько десятков тысяч лет. Лаплас и Лагранж продлили устойчивость движения планет до миллионов лет. В свое время это вызвало бурный энтузиазм в образованных кругах, результат Лапласа-Лагранжа назвали теоремой об устойчивости Солнечной системы. Забавно, что теорема эта приятна и теистам (Часовщик создал часы высочайшего совершенства), и атеистам (Часовщик не нужен, по Пушкину — и без него все шло своим порядком). Во второй половине XX в. советские математики А. Н. Колмогоров, В.И. Арнольд и независимо их американский коллега Ю. Мозер продлили время устойчивости до миллиардов лет. Их результаты уточняются и сейчас, но главное уже сделано.
Любопытно, что факт устойчивости по крайней мере земной орбиты давно известен геологам и палеонтологам. С человеческой точки зрения климат Земли сильно менялся с геологическими эпохами. Но во всяком случае океаны никогда не покрывались сплошным льдом и никогда не нагревались выше 40°С. В первом случае океаны бы никогда не растаяли, так как бело-голубая Земля отражала бы почти все падающее на нее излучение Солнца в космос. Некому было бы читать эту книгу, да и написать ее. Во втором случае мы бы наблюдали грандиозные вымирания растительного и животного мира, по сравнению с которыми гибель динозавров показалась бы мелкой неприятностью. Следовательно, Земля получала от Солнца примерно столько же энергии в прошлом, сколько она получает сейчас. Последние 3 млрд. лет Солнце обладает почти постоянной светимостью. Значит, большая полуось и эксцентриситет земной орбиты существенно не менялись.
Стоит обратить внимание на слова примерно и почти. Колебания эксцентриситета амплитудой в 0,03—0,04 имеют место; ими, согласно хорошо аргументированной гипотезе югославского ученого М. Миланковича, объясняются ледниковые периоды в плейстоцене.
Итак, орбиты восьми больших планет около 4 млрд. лет назад приняли современный вид и с тех пор оставались примерно такими же. А Плутон? Мы помним о его резонансе с Нептуном. Оказывается, в резонансном случае эволюция сильно зависит от фазы. Орбиты Плутона и Нептуна близки к пересечению, расстояние между ними менее 2 а.е. При сближении Плутон перешел бы на существенно другую орбиту, а после серии сближений упал бы на Нептун, или Солнце, или был бы выброшен за пределы Солнечной системы. Но фазы «подобраны» так, что в точках сближения орбит планеты никогда не бывают вместе, расстояние между ними всегда больше 18 а.е. Доказано, что такое состояние длится многие миллиарды лет, орбита Плутона устойчива и сохраняет резонансность с Нептуном. Желающие могут связать это с мудростью Часовщика; нежелающие — с естественным отбором. Начиная с 1992 г. открыты уже сотни планеток диаметрами порядка сотен километров, двигающиеся по схожим с плутоновой орбитам. Они устойчивы, поэтому мы их и видим. Множество же тел, попавших на неустойчивые орбиты, исчезло так, как описано чуть выше.