Скаляр Риччи для диагонального метрического тензора

Мне было интересно, существует ли общая формула для вычисления скаляра Риччи для любой диагонали. н × н метрический тензор?

Для некомпактного многообразия используйте уравнение потока Риччи:
0 г я Дж "=" 2 р я Дж
Для компактного многообразия используйте нормализованное уравнение потока Риччи:
0 г я Дж "=" 2 р я Дж + 2 н р а в г г я Дж
Здесь, н означает н измерения пространства-времени.

Ответы (3)

Предыдущий ответ правильный, но не дает практического алгоритма для людей, потому что вычислять тензор кривизны — это кошмар. Вам нужен хороший алгоритм раздачи, или же вам нужен пакет символьных манипуляций. Я предпочитаю ручные расчеты для симметричных Анзацев, потому что они всегда показательны.

Традиционный упрощенный метод заключается в использовании форм кривизны, и этот метод описан у Misner Thorne and Wheeler. Это необходимо для понимания Kerr Solutions. Также стоит изучить формализм Ньюмена-Пенроуза, потому что он дает физическое понимание. Я предпочитаю использовать свой собственный математически неизящный самодельный метод, потому что большая часть упрощения в продвинутых методах на самом деле связана только с использованием того, что в литературе по информатике называется «вычислениями с разреженными матрицами».

Если у вас есть матрица, состоящая в основном из нулей, например кривизна по диагонали, вам не следует записывать ее в матричной форме, если только вы не хотите накачать сильные мышцы руки для письма. Ввести нековариантные базисные тензоры "l_{ij}", которые отличны от нуля в позиции i,j. Затем запишите метрику в формате «в основном плюс» следующим образом:

г мю ν "=" А л 00 + Б л 11 + С л 22
г мю ν "=" 1 А л 00 + 1 Б л 11 + 1 С л 22

Где я перешел к трем измерениям, чтобы избежать безнадежно длинного ответа, и где я надеюсь, что обозначения ясны. Для теоретической элегантности базисный тензор л 00 действительно должно быть написано как л мю ν 00 если вы хотите, чтобы это соответствовало обычным соглашениям об указателях, но поскольку цель состоит в том, чтобы мышцы письма были как можно более дряблыми, не делайте этого. Поскольку l смехотворно зависят от координат, вы можете выразить смехотворно нетензорные объекты, такие как коэффициенты связи и тензор энергии псевдонапряжения.

Расчет соединения

Есть приемы расчета связи, вроде вывода геодезического уравнения, но я не буду их использовать. Если вы используете базис-тензоры, получить коэффициенты связи совсем не составит труда, а с практикой вы сможете выполнять большую часть работы в уме для более простых анзацев.

Во-первых, дифференцировать метрику. Поскольку «диагональ» не является большим упрощением, я буду предполагать «диагональность и зависимость только от x1 и x0». Я буду использовать штрих для дифференцирования относительно Икс 1 , и точка для дифференцирования по Икс 0 :

г мю ν , α "=" А ˙ л 000 А л 001 + Б ˙ л 110 + Б л 111 + С ˙ л 220 + С л 221

Урок в том, что дифференциация тривиальна. Обратите внимание, что это симметрично по первым двум индексам и ничего особенного по третьему индексу. Символы Кристоффеля симметричны по двум последним индексам и ничего особенного по первому. Вы передаете симметрию между позициями индекса следующим образом:

п я | Дж к "=" Вопрос я Дж | к + Вопрос я к | Дж Вопрос Дж к | я

Где P симметричен в двух последних позициях, а Q симметричен в первых двух. Привыкайте к этому, потому что это часто всплывает. Первый член имеет тот же порядок индексов только благодаря правильному соглашению о порядке индексов, второй член принудительно симметрирует вторую и третью позиции, а последний член требуется, чтобы P сохранял всю информацию в Q. Вы можете выполнить эту процедуру на l автоматически, просто заменив л 001 с л 001 + л 010 л 100 , и так далее. Это то, что вы делаете, чтобы получить Г от производной g.

Вот формула для всех нижних индексов Г (так никогда не писалось, потому что Г это не тензор, но я делаю это здесь, просто чтобы не печатать рукой).

Г мю | ν о "=" 1 2 ( А ˙ л 000 А ( л 001 + л 010 л 100 )
+ Б ˙ ( л 110 + л 101 л 011 ) + Б л 111
+ С ˙ ( л 220 + л 202 л 022 ) + С ( л 221 + л 212 л 122 ) )

Чтобы поднять индексы, когда метрика является диагональной, тривиально, вы просто повышаете индекс на l и делите на соответствующую диагональную запись:

Г ν о мю "=" А ˙ 2 А л 00 0 А 2 А ( л 01 0 + л 10 0 ) А 2 Б л 00 1 + Б ˙ 2 Б ( л 10 1 + л 01 1 ) Б ˙ 2 А л 11 0 + Б 2 Б л 111 + . . .

Где остальное должно быть очевидно. С практикой это занимает минуту, чтобы сделать это вручную.

Расчет кривизны Риччи.

Чтобы рассчитать кривизну, важно отслеживать по мере продвижения , потому что это вдвое сокращает работу. В кривизне Римана всегда есть куча хлама Вейля, который вас в основном не волнует.

Я всегда пишу формулу для тензора Римана так:

р ν λ о мю "=" Г ν λ , о мю ( А С ) + Г ν с мю Г λ о с ( А С )

«AS» означает вычесть одно и то же выражение с ν и о поменяны местами. Эта форма обладает тем свойством, что она антисимметрична по нижнему первому и третьему индексам, поэтому это не обычное соглашение для тензора Римана, который антисимметричен по последним двум индексам. Но это легко исправить в конце. Поверь мне, это лучшая конвенция, ты все исправишь в конце.

След Риччи в этом соглашении находится на первых двух индексах:

р мю ν "=" р α мю ν α

Это важно, потому что каждый член, который вы получаете в тензоре Римана, имеет букву «l», и если верхнее число l не совпадает с крайним левым нижним числом, то этот член не вносит вклад в тензор Риччи. Чтобы получить тензор Риччи, вы просто игнорируете все л б с г а с а б и запиши л с г вместо тех л которые имеют а "=" б .

Теперь продифференцируем выражение для Г , прикрепляя индекс в конце. Я продемонстрирую с одним из вкладов, от принятия Икс 1 производная от первого члена:

Г ν λ , о мю ( А С ) "=" . . + ( А ˙ 2 А ) ( л 001 0 л 100 0 )

Где второй член является антисимметризатором. Но теперь взгляните на две буквы l --- есть ли у одной из них совпадающий индекс слева вверху и слева внизу? Да! Так что сотрите верхнее и нижнее левое число из l, и у вас останется вклад в тензор Риччи:

( А ˙ 2 А ) л 01

Все это можно проделать в своей голове, термин за термином. Если вы получите л который л 221 0 это не может способствовать Риччи, потому что нижний первый и третий индексы не совпадают с верхними, если вы получаете л 121 0 потом его убивают антисимметризацией и т.д. и т.п., это все очевидно.

Далее нужно умножить Г с собой и след. Здесь вы отрабатываете все условия. Но это конечный расчет, а не безнадежный. Вы не получаете вклад, если только крайний левый нижний индекс не совпадает с верхним индексом, так что остается только несколько терминов, и, кроме того, вы получаете ноль на некоторых l-членах после антисимметризации и трассировки (в вашей голове). Тогда остается всего несколько оставшихся слагаемых, и это реальные вклады в кривизну, так что избежать их вычисления невозможно.

Первый раз занимает некоторое время, но с практикой это занимает всего несколько минут для более простых анзацев и несколько часов для самых сложных.

Комментарий по поводу эффективности --- в этом методе все еще есть некоторые "чудесные отмены" терминов, поэтому он не принимает во внимание всю присутствующую математическую структуру. Но на практике эти отмены составляют менее 50% оставшихся условий, которые вы получаете, поэтому вы выполняете только двойную работу. Таким образом, не так много возможностей для улучшения в изучении лучших методов. Я считаю, что отслеживать знаки — самая сложная часть, и иногда полезно использовать для этой цели метрику «все плюс», имея в виду, что A должен быть отрицательным.
большое спасибо за четкий и очень полезный ответ! Я практиковал ваш метод на простом случае. И я действительно доволен. Однако я заметил, что если я хочу получить правильный результат для Г , мне нужно разделить на запись, которая соответствует верхнему индексу Г . Например, в вашем уравнении для Г ν о мю , Я кладу
. . . А 2 А л 11 0 . . . Б ˙ 2 А л 11 0 . . .
и я получаю правильные результаты. Я делаю это правильно?
Это повышение индекса--- Я сделал это и в ответе (и сделал глупые ошибки, теперь я вижу, я их исправил--- спасибо): он исходит из г мю ν спереди. Будьте осторожны при переходе к недиагональным метрикам для умножения/отслеживания на г мю ν вместо.

Скаляр Риччи - это просто след тензора Риччи, который, в свою очередь, представляет собой тензорное сокращение тензора кривизны Римана, который может быть выражен в символах Кристоффеля, определяемых локальной метрикой.

р "=" р я я "=" г я Дж р я Дж

р я Дж "=" р к я к Дж

р р о мю ν "=" мю Г ν о р ν Г мю о р + Г мю λ р Г ν о λ Г ν λ р Г мю о λ

Г к л я "=" 1 2 г я м ( л г м к + к г м л м г к л )

как вы видите, когда у вас есть метрика, вы можете рассчитать все эти количества.

Обратите внимание, что нет гарантии ЛЮБОГО упрощения тензора или скаляра Риччи даже в случае диагональной метрики.
Ну, есть сокращение с точки зрения количества необходимых вычислений. Если г является диагональным, поэтому г 1 , поэтому нам нужно вычислить только диагональные элементы Риччи. Кроме того, единственными ненулевыми членами символов Кристоффеля являются те, у которых два из я , к , л то же самое: сокращение количества терминов от О ( Н 3 ) ( Н будучи размерностью многообразия) к О ( Н 2 ) .
Конечно. Но вы не можете ничего ожидать от формы вашего окончательного ответа ни для тензора Риччи, ни для скаляра Риччи. Но ОП только просит р и не надо вычислять р а б для а б в диагональном случае, не говоря уже об упрощении в Г с.

Несмотря на то, что этот вопрос старый, но да, такие формулы есть, вы можете проверить книгу Ландау о классических полях (том 2, параграф 92) или проверить эту работу (там формулы не так хороши, как у Ландау)

Тензор Риччи диагональной метрики

Обновлять:

Так как книга русская и старая, подытожу здесь результат (с небольшими изменениями):

Если предположить, что диагональная метрика имеет вид:

г мю мю ты мю е 2 л мю : ты мю "=" знак г мю мю
затем определите:
л мю ν р л р , ν ( л ν л р ) , мю л р , мю ν
тогда тензор Римана имеет вид:
р р мю р ν "=" { г мю мю ( л мю ν р + л мю , ν л р , мю ) : р мю ν г мю мю л мю мю р + г ν ν л р р мю г р р г мю мю λ р , мю л мю , λ л р , λ г λ λ : р мю "=" ν
и тензор Риччи:
р мю ν "=" { р мю , ν ( л мю ν р + л мю , ν л р , мю ) : мю ν р мю [ л мю мю р + г мю мю г р р ( л р р мю л мю , р λ мю , р л λ , р ) ] : мю "=" ν

Остальные компоненты тензора равны нулю.

Я не смог найти такую ​​формулу в пункте 92, не могли бы вы указать номер страницы? Какое издание вы используете?
@ user41404: Пожалуйста, посмотрите обновление.