Британский физик-теоретик Поль Дирак (1902–1984), один из основоположников квантовой механики, однажды написал:
«Было трудно примирить ньютоновскую теорию гравитации, в соответствии с которой гравитационное взаимодействие распространяется мгновенно, с требованиями специальной теории относительности, и Эйнштейн, работавший над этой проблемой, пришёл к обобщению своей теории относительности — вероятно, величайшему из когда-либо сделанных научному открытию».
Общую теорию относительности принято считать теорией исключительной красоты. Эксперименты и наблюдения, проводившиеся на протяжении многих лет, подтвердили непротиворечивость этой теории. Я опишу одно из испытаний, давшее корректное объяснение «аномального» смещения перигелия Меркурия, которое не подчинялось ньютоновскому закону всемирного тяготения.
Проблема классической теории тяготения Ньютона
Прецессия (или смещение) перигелия (ближайшей к Солнцу точки орбиты планеты) имеет множество причин. Вот две из них:
- наличие других планет, взаимовлияние которых вызывает медленное изменение параметров их орбит. Это главная причина;
- сплюснутость солнечного диска (см. рисунок), которое оказывает значительно меньшее влияние.
Скорость смещения перигелия Меркурия не согласуется с предсказанием ньютоновской теории тяготения. Эту аномалию обнаружил французский астроном и математик Урбен Леверье. В 1882 году были опубликованы результаты расчётов Саймона Ньюкома, согласно которым фактическая скорость смещения расходится с предсказанием ньютоновской теории на 43 угловых секунды. Для объяснения этого феномена выдвигались разные гипотезы, но все они оказывались несостоятельными.
Решение, полностью объясняющее это дополнительное смещение, было найдено в рамках общей теории относительности Эйнштейна.
Использование общей теории относительности для расчёта смещения перигелия Меркурия
Решение Шварцшильда — это решение полевых уравнений Эйнштейна, описывающих геометрию вакуумного пространства-времени вокруг Солнца. Другими словами, метрика Шварцшильда — это метрика Солнечной системы, обусловленная искривлением пространства-времени вблизи Солнца. Такая метрика справедлива, если:
- Солнце рассматривается как невращающийся объект;
- гравитационным полем других планет Солнечной системы можно пренебречь.
Решение Шварцшильда имеет следующий линейный элемент:
где r, θ и φ — это сферические координаты (наглядно показанные на рис. 3), а параметр R = 2M известен как радиус Шварцшильда.
Заметим, что из изотропии метрики мы всегда имеем θ = π/2 (орбиты ограничены экваториальным местом). На самом деле, согласно задаче двух тел (в нашем случае это Солнце и планета), движение тела с учётом центральной силы всегда будет происходить в плоскости. На рис. 4 и рис. 5 показаны два типа систем двух тел, движущихся по орбите. Движение, ограниченное плоскостью, актуально и для закона всемирного тяготения Ньютона, и для теории тяготения Эйнштейна. Следовательно, в нашем анализе будет достаточно рассмотреть только геодезические, находящиеся в этой плоскости.
Третье условие справедливости этого анализа состоит в том, что радиальная координата r должна быть намного больше радиуса Солнца. Это не проблема, так как радиус Шварцшильда для Солнца гораздо меньше радиуса самого Солнца. Для сравнения: радиус Шварцшильда для Солнца составляет примерно 2,95×10³ м, тогда как радиус Солнца приблизительно равен 6,96×10⁸ м.
Симметрии в данном пространстве-времени связаны с сохраняющимися величинами для движущихся в нём частиц и фотонов. Поскольку метрика g решения Шварцшильда независима от времени (или инвариантна относительно сдвига начала отсчёта по времени) и сферически симметрична, энергия массивных частиц и энергия фотона сохраняются. Математически можно представить это следующим образом.
В пространстве-времени с метрикой g свободно падающая материальная частица или фотон удовлетворяет условиям геодезического уравнения (обобщения понятия «прямая» для искривлённых пространств), связанного с этим пространством-временем, которое имеет следующий вид (см. Schutz):
Обратите внимание: так как фотоны тоже будут приниматься в расчёт, то параметр λ не может быть собственным временем τ. Геодезическое уравнение может быть записано и в таком виде:
Заметьте, что:
Уравнения 3 и 4 предполагают, что:
Затем определяем следующее:
Символ ~ в обозначении энергии массивной частицы используется для указания на то, что это энергия на единицу массы (см. Schutz). Аналогично вследствие независимости g от φ сохраняется момент импульса. Мы определяем:
где в левой части момент импульса на единицу массы частицы, а в правой — момент импульса фотона. Теперь нам нужны уравнения для орбиты. Тремя составляющими импульса массивной частицы являются:
Импульсы для фотона:
Теперь используем только что полученные компоненты импульса, подставляем их в уравнения |p|=-m² для частицы и фотона и решаем для dr/dλ. Уравнения для dr/dλ принимают такой вид:
Теперь интуиция подсказывает нам переписать эти уравнения с использованием эффективных потенциалов, а именно:
Эти потенциалы показаны на рис. 7. Обратите внимание: поскольку левая сторона обоих уравнений положительна, эффективный потенциал должен быть меньше энергии. Рис. 7 показывает эффективный потенциал (заметьте, что E и V на рис. 7 указывают те же величины с символом ~) для массивных и безмассовых частиц. На рисунке также показаны точки поворота, где dr/dλ=0, запрещённые зоны, где E < V, и круговые орбиты (устойчивые и неустойчивые), где dV²/dr=0.
Смещение перигелия Меркурия
Поскольку наша цель — вычислить смещение перигелия Меркурия, далее будем рассматривать только движение массивных объектов.
Устойчивая круговая орбита возникает на минимуме эффективного потенциала. Пусть M будет массой Солнца. Дифференцируя эффективный потенциал, приводя результат к нулю и решая для r, получаем радиус устойчивой круговой орбиты:
В ньютоновской механике полная орбита планеты на круговой орбите возвращается к своему исходному φ. Используя тот факт, что круговые орбиты имеют E²=V², и задействуя полученные нами выражения, можно вычислить время, требующееся планете, чтобы достичь Δφ=2π, то есть период P:
Итак, в общей теории относительности вращающаяся планета не возвращается в свою начальную точку. Если релятивистские эффекты незначительны, то у нас должен быть эллипс, который медленно вращается вокруг своего центра. И теперь мы можем изучить движение перигелия орбиты. Для этого выполняем три быстрых вычисления (см. Schutz):
- выведем выражение для dφ/dλ, используя момент импульса;
- выведем выражение для dt/dλ, используя энергию на единицу массы;
- определим новую переменную u ≡ 1/r.
Подставляя их в уравнение 10, получаем:
А теперь определим y, то есть отклонение от круглости:
Для ньютоновской орбиты y=0. Чтобы получить релятивистское выражение, подставим уравнение 15 в уравнение 14 и отбросим переменные y³ третьего порядка. Получаем следующее уравнение для почти круговой орбиты:
Решение принимает вид:
где B зависит от начальных условий. Из аргумента косинуса заключаем, что орбита возвращается к тому же радиусу, когда Δ(kφ) = 2π. Полученный результат отличается от ньютоновского наличием k, не равного 1! Если релятивистские эффекты незначительны, можно сделать ещё несколько простых приближений, чтобы получить:
В частности для Меркурия получается сдвиг на 43 угловых секунды в год — значение, которое (как уже упоминалось в начале статьи) было подтверждено экспериментально.
Похоже, даже Эйнштейн был ошеломлен результатом своих вычислений. Ошеломлён настолько, что не мог дальше работать несколько дней. По его собственным словам, он был «вне себя от радости».
Читайте также:
- Почему 0,99999… равно 1
- Плотность вероятности - это не сама вероятность
- 3 случая, когда линейная модель может ошибаться
Перевод статьи Marco Tavora: Einstein, and the Most Beautiful of All Theories
