2024 ВЕСТНИК САНКТ-ПЕТЕРБУРГСКОГО УНИВЕРСИТЕТА Т. 20. Вып. 1
ПРИКЛАДНАЯ МАТЕМАТИКА. ИНФОРМАТИКА. ПРОЦЕССЫ УПРАВЛЕНИЯ
ПРИКЛАДНАЯ МАТЕМАТИКА
УДК 517.958:531.111 МБС 83С05, 83С25
Решение краевой задачи для уравнения Эйнштейна внутри шара с однородной плотностью
О. И. Дривотин
Санкт-Петербургский государственный университет,
Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9
Для цитирования: Дривотин О. И. Решение краевой задачи для уравнения Эйнштейна внутри шара с однородной плотностью // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2024. Т. 20. Вып. 1. С. 4-9. https://doi.org/10.21638/11701/spbu10.2024.101
Сформулирована краевая задача для уравнения Эйнштейна, описывающая гравитационное поле шара с однородной плотностью распределения массы. При формулировке краевых условий на поверхности шара используется известное решение Шварцшиль-да для пустого пространства-времени, характеризующее метрический тензор вне шара. Координаты, в которых записано решение Шварцшильда для внешней области, отличаются от вводящихся при записи уравнений для компонент метрического тензора внутри шара. Найдена связь между этими координатами, позволяющая использовать решение Шварцшильда при постановке краевых условий для решения внутри шара. Получено решение данной краевой задачи для случая слабого поля. Ключевые слова: уравнение Эйнштейна, метрический тензор, шар с однородной плотностью.
1. Введение. Известны различные решения уравнения Эйнштейна [1], описывающие гравитационное поле в пустом пространстве времени, в пространстве-времени с электромагнитным полем и в шаре, заполненном жидкостью, которая находится под давлением. Настоящая работа посвящена гравитационному полю, создаваемому некоторым распределением массы. Рассмотрен один из простейших случаев, когда масса равномерно распределена внутри некоторого шара и не движется. Сформулирована соответствующая краевая задача для уравнения Эйнштейна. При формулировке краевых условий на поверхности шара учитывается, что вне шара метрический тензор представляет собой решение Шварцшильда для пустого пространства-времени. При этом одна из координат, которые вводятся внутри шара, отличается от соответствующей координаты, используемой при записи решения Шварцшильда. Исходя
(¡5 Санкт-Петербургский государственный университет, 2024
из условия непрерывности и непрерывной дифференцируемости компонент метрического тензора на поверхности шара, найдена связь между данными координатами, что позволяет задать краевые условия на поверхности шара для поля внутри шара. Получено решение краевой задачи для случая слабого поля.
2. Постановка задачи. Введем такое распределение массы, которое можно рассматривать как шар с однородной плотностью. Предположим, что в некоторой локально лоренцевой системе отсчета [2, 3] функция распределения массы в конфигурационном пространстве [4], ассоциированном с этой системой отсчета, равна
Í во, r < R, в \0, r > R,
и масса находится в покое в каждой точке конфигурационного пространства. Здесь во = const, r — расстояние от центра шара. Координата r и две угловые координаты д, ф образуют систему координат в конфигурационном пространстве, которую можно рассматривать как сферическую систему координат с началом в центре шара.
Задача, которую нужно решить, состоит в том, чтобы найти метрический тензор g при следующих условиях: он должен удовлетворять уравнению Эйнштейна [5, 6]
R - 1 gR = aT (1)
внутри и снаружи шара, компоненты gik должны быть непрерывно дифференцируемыми функциями координат на поверхности шара и g должен быть метрическим тензором плоского пространства-времени при r ^ ж. В этом уравнении R обозначает тензор Риччи, T — тензор энергии-импульса, R = gikRik — скалярную кривизну пространства-времени, a — некоторую постоянную.
В разных исследованиях тензор T вводится различным образом. Здесь будем следовать работам [7-9], где T определяется так, что T00 есть функция распределения энергии в конфигурационном пространстве [4], представляющем собой поверхность x0 = const. Таким образом, внутри шара Too = во с2. Здесь с — фундаментальная константа [3], известная как скорость света в вакууме. Все остальные компоненты тензора T равны нулю, поскольку масса покоится.
3. Решение внутри шара. Предположим, что rg ^ R, где rg = 2GM/c2 — гравитационный радиус, G — гравитационная постоянная, M = 4^R3eo/3 — масса шара.
Будем искать внутреннее решение с точностью до членов первого порядка малости по rg/r в виде
\\gik\\ = diag(l + ar2 + b, -1, -r2(1 + dr2), -r2(1 + dr2) sin2 d)
в координатах x0, r, д, ф. Здесь a,b,d — некоторые постоянные, такие, что величины ar2,b,dr2 имеют первый порядок малости по rg/r. Малость добавленных членов к компонентам метрического тензора плоского пространства-времени означает, что рассматривается случай слабого поля.
Вычисляя символы Кристоффеля второго рода Г- с точностью до членов первого порядка малости, получаем, что
pr
г 00
r 1 а
ar.
„2
Г0
r r0
r(1 + 2dr2), Г% = 1 (1 + dr2),
Гфф = - r(1 + 2dr2) sin2 в, Гфф = 1 (1 + dr2),
Гфф = — sin в cos в,
Гфф = cot в.
Все остальные символы Кристоффеля, отличающиеся от символов (2) не только перестановкой нижних индексов, равны нулю.
Компоненты тензора Риччи могут быть записаны в виде [5]
Rij =
дт% 1 d2 ln \g\ , 1 m д ln \g\
+ -Г"
i
dxk 2 dxidxj 2 ij dx
_ pk pm
r mjr ik,
где g = det ||gik||. В рассматриваемом случае
g = — r4[1 + (a + 2d)r2 + b] sin2 в.
Тогда имеем, что
R R
фф
dr дГг дг
R
_ дГг 00 = ,2
00 1 lp r Г 2Г 00"
д Нд
dr
— 2Г000Г0г = 3a,
1 д ln \g\ _ 2Р0 pr _ 2Рв рв _ 2РФ РФ
2 дг2 21 rr1 0r 21 гв1 гв 21 ГфГ Гф
—a — 6d,
l д2 ln\g\
+
2 дв2 дгф
+ 1 Pr + 2 Г в!
д ln \g\ ' dr
— 2Г^в Гвв — (Гфф)2 = (—a — 6d)r2,
+ lpr д ln \ g \ + l рв д ln \ g \ 2pr рф _ 2Рв рФ = + 21 фф дr + 21 фф д в 21 ФФГ гф 21 ффг вф
( a 6d)r2 sin2 в,
R = 6a + 18d.
(3)
ar
Внедиагональные компоненты тензора Риччи равны нулю.
Подставляя величины (3) в уравнение Эйнштейна (1), получим следующие равенства:
—9d = a.Q0<?, (4)
—a — 6d +2(6a +18d)=2a + 3d = 0 (5)
для компонент с индексами 00 и для остальных компонент соответственно. Следовательно, если 6a = a^0c2, уравнение Эйнштейна удовлетворяется.
4. Внешнее решение. Решение снаружи шара есть решение Шварцшильда
Ц~дкГ11 = diag(1 — Ц-, —1 — r, — f2, — f2 sin2 в), (6)
которое удовлетворяет уравнению Эйнштейна и условию на бесконечности. Символ «тильда» означает, что компонента берется в координатах x0 = ct, r, в, ф. Координаты x0 в, ф те же, что и ранее, а координата f — та, которая используется в выражении для решения Шварцшильда. Поскольку rg ^ r, а на поверхности шара r и r отличаются на малую величину, то, как увидим далее, rg/f также мало. В выражении (6) и далее мы удерживаем только члены не выше первого порядка малости по rg/r.
При формулировке краевой задачи для внутренней области можно использовать решение Шварцшильда для внешней области, но следует иметь в виду, что координата r отличается от координаты r и ее значение на поверхности шара заранее
неизвестно. Тем не менее можно определить краевые условия для внутренней задачи, используя производную одной координаты по другой.
Примем во внимание, что координата г есть расстояние от центра шара. Тогда легко увидеть, что дО1 = —1 и дГ/дг = 1 — гд/(2г). Для других компонент получим следующие выражения: д^ = Го°ои\ дОи = д°1\ дфФ = Гф^-
Нетривиальные производные этих компонент по r равны
dgo°out dr
dg§U dr_
dr dr
-g 1 -
2r
dgOU = dgout dr_ = -2_ / - rg
dr dr dr V 2r
с точностью до величин первого порядка малости по rg/r.
Компоненты с индексами фф отличаются от компонент с индексами вв только множителем sin2 в, и учет их дает тот же результат, как для компонент с индексами вв. По этой причине эти компоненты далее не рассматриваем.
5. Краевые условия на поверхности шара. Непрерывная дифференцируе-мость решения означает, что компоненты решения снаружи и внутри шара и их производные по r должны быть равны на поверхности шара r = R, которую обозначим через S:
gools = go°outU,
dgoo ^^ = dgoo dr dr
-U, gee S =
°ut
he ls,
dgee, _ dgt
dr
ls
Эти равенства дают, что
1 + aR2 + b = 1 - -í, R
2aR = ^ (1- \
R2
2R
R2(1 + dR2) = R2,
R + 2dR3 = R ( 1 --^í 2R
(7)
(8) (9)
(10)
Здесь Д обозначает величину координаты Г на поверхности шара. Из уравнения (9) находим, что
Д = Д[ 1 + 1 dД2
Подставляя R в уравнение (10), получим, что
d =
3R3
8п Ggo
9 c
2
(11)
Уравнение (8) дает, что
2Д3
Эти значения а и d удовлетворяют условию (5). Используя равенства (4) и (11), можно видеть, что
8пС
г4 '
°ut
r
g
а
Величина а найдена А. Эйнштейном [5] и используется обычно в его уравнении. Из (7) можно получить, что
6 .
2R
6. Заключение. Найденное решение согласуется с классическим выражением для гравитационного потенциала Ф шара с однородной плотностью массы
) \2пСв0(г2 - R2)/3 - GM/R, r < R, (r)=j - GM/r, r>R,
при условии Ф(то) = 0, если принять во внимание, что для слабого гравитационного поля [5]
Ф
goo = 1 + 2-2.
с2
Полученное решение можно использовать в качестве примера при анализе законов сохранения для гравитационного поля, в которых взаимодействие массы с создаваемым ею полем дает вклад в энергию и импульс гравитационного поля [10].
Литература
1. Stephani H., Kramer D., MacCallum M. A. H., Hoenselaers C., Herlt E. Exact solutions of Einstein's field equations. Ed. 2. Cambridge: Cambridge University Press, 2003. 701 p.
2. Drivotin O. I. Rigorous definition of the reference frame // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2014. Т. 10. Вып. 4. С. 25-36.
3. Дривотин О. И. Об определении геометрии пространства-времени // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2022. Т. 18. Вып. 3. С. 316-327. https://www.doi.org/10.21638/11701/spbu10.2022.302
4. Дривотин О. И. Covariant description of phase space distributions // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2016. Т. 12. Вып. 3. С. 39-52. https://www.doi.org/10.21638/11701/spbu10.2016.304
5. Einstein A. Die grundlage der allgemeinen relativitatsteorie // Ann. d. Phys. 1916. Vol. 49. P. 769-822.
6. Дривотин О. И. Математические основы теории поля. СПб.: Изд-во Санкт-Петербургского университета, 2010. 168 с.
7. Ландау Л. Д., Лифшиц Е. М. Теория поля. М.: Наука, 1973. 504 с.
8. Wang C.-C. Mathematical principles of mechanics and electromagnetism. Pt B. Electromagnetism and gravitation. New York: Plenum, 1979. 386 p.
9. Gron O., Hervik S. Einstein's general theory of relativity. New York: Springer Verlag, 2007. 558 p.
10. Дривотин О. И. О плотности потока импульса гравитационного поля // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2021. Т. 17. Вып. 2. С. 137-147. https://www.doi.org/10.21638/11701/spbu10.2021.204
Статья поступила в редакцию 24 октября 2023 г. Статья принята к печати 26 декабря 2023 г.
Контактная информация:
Дривотин Олег Игоревич — д-р физ.-мат. наук; o.drivotin@spbu.ru
The Einstein equation solution inside a ball with uniform density
O. I. Drivotin
St. Petersburg State University, 7-9, Universitetskaya nab., St. Petersburg, 199034, Russian Federation
For citation: Drivotin O. I. The Einstein equation solution inside a ball with uniform density.
Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2024, vol. 20, iss. 1, pp. 4-9. https://doi.org/10.21638/11701/spbu10.2024.101 (In Russian)
A great number of solutions of the Einstein field equation are known. They describe the gravitational field in the empty spacetime, in the spacetime with electromagnetic field and for a ball filled with a liquid under pressure. The present work is devoted to gravitational field generated by some mass distribution. One of the simplest cases is considered, when mass is uniformely distributed inside a ball and is not moving. The boundary problem for the Einstein equation is formulated. Solution outside the ball is the Schwartzschild solution in vacuum. The coordinates at which the Schwartzschild solution is written are different from the coordinates used in equations for components of the metric tensor inside the ball. Relations between internal and external coordinates are found on the ball surface. They allow to use the Schwartzschild solution for formulation of boundary conditions for internal solution. The solution of the boundary problem is found for the case of weak field. This solution can be used as an example in the analysis of laws of conservation for the gravitational field, in which interaction of mass with field generated by the mass gives a contribution to momentum and energy of the gravitational field.
Keywords: the Einstein equation, metric tensor, ball with uniform mass density. References
1. Stephani H., Kramer D., MacCallum M. A. H., Hoenselaers C., Herlt E. Exact solutions of Einstein's field equations. Ed. 2. Cambridge, Cambridge University Press, 2003, 701 p.
2. Drivotin O. I. Rigorous definition of the reference frame. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2014, vol. 10, iss. 4, pp. 25—36.
3. Drivotin O. I. Ob opredelenii geometrii prostranstva-vremeni [On the determination of spacetime geometry]. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2022, vol. 18, iss. 3, pp. 316-327. https://www.doi.org/10.21638/11701/spbu10.2022.302 (In Russian)
4. Drivotin O. I. Covariant description of phase space distributions. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2016, vol. 12, iss. 3, pp. 39-52. https://www.doi.org/10.21638/11701/spbu10.2016.304
5. Einstein A. Die grundlage der allgemeinen relativitatsteorie. Ann. d. Phys., 1916, vol. 49, pp. 769822.
6. Drivotin O.I. Matematicheskiye osnovy teorii polya [Mathematical foundations of the field theory ]. St. Petersburg, St. Petersburg University Press, 2010, 168 p. (In Russian)
7. Landau L. D., Lifshitz E. M. Teoriya polya [Field theory ]. Moscow, Nauka Publ., 1973, 504 p. (In Russian)
8. Wang C.-C. Mathematical principles of mechanics and electromagnetism. Pt B. Electromagnetism and gravitation. New York, Plenum Publ., 1979, 386 p.
9. Gron O., Hervik S. Einstein's general theory of relativity. New York, Springer Verlag Publ., 2007, 558 p.
10. Drivotin O. I. O plotnosti potoka impulsa gravitatsionnogo polya [On momentum flow density of the gravitational field]. Vestnik of Saint Petersburg University. Applied Mathematics. Computer Science. Control Processes, 2021, vol. 17, iss. 2, pp. 137-147. https://www.doi.org/10.21638/11701/spbu10.2021.204 (In Russian)
Received: October 24, 2023. Accepted: December 26, 2023.
Author's information:
Oleg I. Drivotin — Dr. Sci. in Physics and Mathematics; o.drivotin@spbu.ru