Когда речь заходит о теории относительности, частенько на ровном месте разрастаются споры, которые были занесены в почву непонимания и обильно удобрены мифами, недосказанностью и недостаточной математической подготовкой. Даже на лекциях от некоторых профессоров можно услышать, что детище гения Эйнштейна не имеет практической пользы, а на робкие попытки пролепетать что-то про спутниковые системы навигации они пренебрежительно отмахиваются, дескать, там все сложно и двояко.  Так что совершенно естественно желание попробовать провести некоторые расчеты самолично, потрогать формулы, покрутить параметры, чтобы постепенно заложить интуицию в столь горячей теме. Вся красота и непостижимая мощь общей теории относительности аккуратно упакована в уравнение гравитационного поля Эйнштейна: Оно говорит, что искривление пространства-времени определяется материей. То есть, гравитация воспринимается не как сила, а как следствие искривленной геометрии царства 4d, а искривления порождены материей (тензор энергии-импульса в правой части уравнения). Движению объекта можно поставить в соответствие траекторию в четырехмерном пространстве-времени - мировую линию. Просто представьте себе четырехмерный брусок в котором статично зависли макаронины: галактики, планеты или молекулы - последние на некоторое мгновение переплетаются, чтобы образовать читателя этой статьи, а затем вновь разносятся в пространстве, очерчивая пути кусочков кожи и турбулентных потоков выдохнутого воздуха. Вся эта лапша подчинена кривизне пространства-времени, и она же, собираясь в массивные сплетения, эту кривизну задает. Для работы с искривленной геометрией у нас есть математические инструменты: метрический тензор, символы Кристоффеля, тензор кривизны Римана, тензор Риччи и скалярная кривизна. Для того, чтобы создать интуитивный образ, вообразите натянутый кусок ткани деформированный массивным телом... Хотя есть куда более умозрительный и математически точный вариант в видео A new way to visualize General Relativity: И давайте тогда сразу приведем список полезных источников для погружения в тему: Вспомнить специальную теорию относительности помогут видео от MinutePhysics  на русском. Вспоминаем сущность линейной алгебры с 3Blue1Brown. Постигаем магию тензорной алгебры с серией статей на хабре. Вывод уравнения Эйнштейна смотрим в Einstein Field Equations - for beginners! Понять математику ОТО поможет серия The Maths of General Relativity от ScienceClic English (настойчиво советую просмотреть весь материал канала − на редкость качественная постановка и визуализация сложных тем), книжки Introduction to general relativity 2010, Gerard ’t Hooft, A first course in general relativity 2009 Schutz, General Relativity, Black Holes, and Cosmology 2020, Andrew J. S. Hamilton   и более щадящая Cosmology, Daniel Baumann   Видео VSauce: вниз - это куда? если вы боитесь формул. На русском книг не предложу, но, думаю, если мы хором позовем @Tyusha, то она может чего посоветовать. Все мы со школьной скамьи приучены к метрике Минковского: решая задачки в евклидовой геометрии, постепенно привыкаешь к плоскому пространству. Но чего нам действительно не хватало, так это построения геодезических линий на глобусе или небесной сфере. Собственно, сферическая симметрия это вотчина метрики Шварцшильда. Гравитационное поле массивного сферического тела, черные дыры - это все сюда. Есть еще разные другие метрики: для заряженных, для вращающихся тел, для расширяющейся вселенной и т. д.  В плоском пространстве-времени все довольно просто: Событие задается вектором с четырьмя компонентами Зависящие от наблюдателя пространственно-временные координаты превращаются в инвариантные линейные элементы с помощью тензора метрики (мы держим в уме соглашение Эйнштейна о суммировании и Лоренц-инвариантность) Тензор метрики Минковского - это диагональная матрица 4х4 Определяем интервал в этом плоском четырехмерьи. Обобщение теоремы Пифагора в искривленное пространство-время. Записываем его через собственное время - часики-то тикают, и у каждого свои. Вводим 4d-скорость. dx/dτ - скорость движения в направлении оси Х, dt/dτ - скорость изменения временнóй компоненты и т.д. Допустим мы хотим решить уравнение Эйнштейна (найти метрический тензор g_{μν} ) для точки на некотором удалении от статичного незаряженного сферического тела. Само по себе решение в лоб трудоемко, но правильный выбор системы координат и учет симметрий чрезвычайно упрощают получение результата. Так что вывод метрики Шварцшильда вполне посильный труд. Получается, расстояние между двумя событиями в пространстве-времени в окрестности массивного сферического тела в вакууме имеет форму Если занулить массу М, получим пустое плоское пространство-время Минковского. Чем ближе мы к массивному телу, тем сильнее ощущаем кривизну. Внутри тела метрика не работает, но если оно компактное, то мы можем найти особое положение, при котором первое слагаемое стремится к нулю, и, соответственно, второе уходит в бесконечность. Как вы догадались, речь идет о радиусе Шварцшильда - горизонте черной дыры. Для разогрева, попробуйте рассчитать в этой метрике замедление времени для искусственного спутника Земли и сравните с результатами какого-нибудь эксперимента. Запишем тензор метрики в нормальных единицах (c=G=1). Кстати, в этой геометрической системе время и масса имеют размерность длины. Ответьте, чему равен ваш возраст в метрах? А сколько километров весит Солнце? Воспользуемся симметрией смещения во времени и вращения вокруг оси z, чтобы ввести сохраняющиеся величины: Это энергия и угловой момент на единицу массы покоя. К слову, они нам еще могут пригодиться, если мы вдруг надумаем поиграть с черными дырами и червоточинами. Итак, сохранение углового момента подразумевает, что орбита лежит в заданной плоскости, что позволяет выбрать для переменной θ конкретное значение, скажем θ = π/2. Чтобы перейти к скоростям, разделим выражение для интервала в метрике Шварцшильда на dτ² и получим выражение для полной энергии тела на орбите Уравнение полной энергии слагается из кинетической и потенциальной, так что у нас есть аналитическое выражение для гравитационного потенциала в релятивистском и классическом случаях.  Давайте их визуализируем! Используем язык Julia (для питонистов в конце тоже будет ссылка).  Окей, как и предполагалось, в классической механике существует только два различных типа орбит: замкнутые: эллипсы и круги (орбиты 2 и 3)  незамкнутые орбиты рассеяния: гиперболы (орбита 1)  В общей теории относительности существует три типа орбит: спиральные: радиальное погружение (орбита 1) незамкнутые (орбита 2) замкнутые: прецессирующие эллипсы (орбита 3) Очевидно, устойчивость орбиты определяется характером экстремума.  Вспомним формулу полной энергии и запишем ее в привычной размерной форме  Заметим, что при малых скоростях четвертое слагаемое дает пренебрежимо малый вклад, и формула будет описывать классическое орбитальное движение. А дальше, выполнив дифференцирование последних двух слагаемых, имеем радиальную силу Опять-таки, в классическом случае, когда скорость света кажется бесконечной, у нас будет обычная формула из школьной физики. Вспомнив второй закон Ньютона, записываем уравнение движения или, на плоскости и с понижением порядка: Здесь появились константы А и В, чтобы можно было переобозначив их легко вернуться к размерным переменным. Четыре дифурки решаем Рунге-Куттой-4: и давайте нарисуем все три типа релятивистских орбит Что будет если стартануть с горизонта событий? С позиции за горизонтом? Что произойдет при достижении r = 0? Каждое столетие перигелий орбиты Меркурия увеличивается на 5300 угловых секунд, но только около 5260 угловых секунд могут быть объяснены ньютоновской механикой. Да, были способы разной степени костыльности, но наша цель - проверить, может ли общая теория относительности объяснить оставшиеся 40 угловых секунд. Запишем уравнение движения в форме удобной для подстановки измеряемых констант: Здесь R = 1 AU расстояние между Землей и Солнцем, T = 1 земной год, m - масса Меркурия, M - масса Солнца. Еще полезно будет знать ряд измеримых параметров: Меркурий Перигелий 0.307 499 АЕ Афелий 0.466 697 АЕ Год 0.240 846 земн. лет Макс. орбитальная скорость 58.98 км/с Мин. орбитальная скорость 38.86 км/с В перигелии скорость максимальна. В афелии она минимальна. В обоих случаях скорость нормальна к вектору расстояния между Солнцем и Меркурием.  Присмотримся к константе B. Эта величина характеризуется небольшим числом, что означает, что общая релятивистская поправка на ньютонову орбиту невелика. На самом деле, если мы продолжим наивно вычислять орбиту Меркурия, интегрируя уравнение движения, мы получим неправильную скорость прецессии из-за ошибок численного округления. Поэтому мы рассмотрим различные значения для B для набора орбит с более высокой скоростью прецессии, выполним линейную регрессию и экстраполируем прецессию Меркурия. Итак, мы выбираем ряд значений В, чтобы потом оценить для них прецессию. Чтобы удостовериться, что диапазон подобран правильно, проверим, что при В = 0 прецессия отличается от оной при В > 0 на несколько порядков: Чтобы найти наклон кривой, характеризующей изменение прецессии, воспользуемся пакетом Optim.jl Найденный параметр P задает наклон прямой, и теперь мы можем оценить релятивистскую поправку для прецессии меркурия! Ожидаемое значение 43 угловые секунды за век, так что очень неплохо. Для остальных планет поправки куда меньше, но их тоже удается воспроизвести, причем даже аналитически. Методики нахождения смещения перигелия смотрим в предложенных ранее книгах и в дополнительных источниках: Two conserved angular momenta in Schwarzschild spacetime geodesics   The Precession of Mercury’s Perihelion, Owen Biesel   A primer to numerical simulations: The perihelion motion of Mercury   Computational physics examples as IPython Notebooks   Теорию опробовали на иных космических телах в центре нашей галактики и за ее пределами, так что любой желающий может найти в открытом доступе наблюдаемые данные, выбрать аналитику или численные методы себе на вкус и удостовериться в эффективности теории относительности. P.S. На всякий случай добавим ссылочки - если вдруг кому очень захочется всех заверить в никчемности релятивизма, то ему следует по пунктам разобрать каждый эксперимент, указать ошибки и предложить способ решить их. https://en.wikipedia.org/wiki/Tests_of_special_relativity https://en.wikipedia.org/wiki/Tests_of_general_relativity What is the experimental basis of Special Relativity? Experimental Tests of General Relativity В любом случае, все эти теории подвергаются действию естественного отбора, так что хорошим тоном будет воздержание от подъема агрессивного бурления в комментах.