Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Физика и философия подобия от преонов до метагалактик

..pdf
Скачиваний:
7
Добавлен:
19.11.2023
Размер:
41.03 Mб
Скачать

442

§48.4. Метрика внутри однородного шара. Квазигалилеев подход

Приблизительно в таком же виде выражаются остальные компоненты тензора Rik . С целью их упрощения наложим на метрический тензор gik так называемые ус­ ловия изотермичности:

что дает 4 соотношения между компонентами тензора hjk . Вновь учитывая, что

стоит использовать только компоненты g n

при g 00 = 1 и g pp = — 1 = 1, 2 , 3),

для 5 = 0 получаем одно из соотношений в виде:

 

 

1 д

+ h22 + Л 33)

dhn.

дИт

dhm

+

= -J1L +

+

^ 2 1

2 cdt

 

дх

ду

Ъъ

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

где □ = А — 1=--8j - — оператор Д’Аламбера,

Сdt

акомпоненты тензора hjk выражены через гармонические координаты.

Вследствие малости hjk скалярная кривизна R равна:

R ~ 8

~ *оо "

Кц ~ -^зз =

^?hss)'

 

 

1

5 = 1

Запишем уравнение для метрики Эйнштейна — Гильберта в дважды

ковариантных индексах:

 

4 ~ \ s lkR - -Zf - T ik.

(573)

В качестве тензора плотности энергии-импульса материи Tik

возьмем только

тензор tik (537) для несжимаемой жидкости (при L = const согласно (546), мы же примем здесь 1 = 0 ) без учета других видов энергии:

'ft “ (Л

-S< kP»’

(574)

 

где р 0— плотность покоящегося вещества, Р0— изотропное давление, с — скорость света,

Uj — ковариантная 4-скорость.

В нашем приближении можно считать, что в (574) gik = rjjk,

ui

= S iku

с

V

 

~ Vik(-

У*/сг ’

) ~ (с , - У),

 

 

- ур,

J l -

- V1 /с г

 

U„ = c, Up -

р= 1,2,3.

 

 

 

Подставляя Rik' 4 k

в (573), получаем следующие уравнения:

 

 

 

^00 т(^оо

^fiss)

_ 16з гу р .

(575)

 

 

 

§48.4. Метрика внутри однородного шара. Квазкгалилеев подход

443

v =

- ^

т г ( р о + 4 ) ^ .

 

(576)

 

 

 

с

 

с

 

 

v + ^оо -

2 * J

=

L

+ % ^

р + л ]

(577)

^

^=>i

J

С

С

J

 

□ V

=

^ (

Р

о

+ 4 )Р *^9 , ^ .

(578)

здесь р, q — 1, 2 , 3.

Мы видим, что компоненты добавочного метрического тензора hjk удовлетворя­ ют волновым уравнениям, признаком которых является даламбертиан. Тем самым проявляются колебания метрики пространства-времени как отклонения от стацио­ нарного значения, фиксируемого тензором rjjk.

Рассмотрим уравнения (575) — (578) в статическом случае, когда произведениями

скоростей движения

вещества

VpVq можно

пренебречь из-за их малости, а

даламбертианы превращаются в лапласианы. Уравнения (577), (578) примут вид:

 

 

hpp +

 

 

 

 

 

_ 16яуР0

(579)

 

 

 

 

 

 

 

 

 

 

 

 

 

Ahpq = 0.

 

 

(580)

Уравнение (580) имеет общее решение в виде:

 

 

 

 

 

h

 

= г

+ ^£2-

 

 

 

 

npq

^pq

т

г

,

 

где Cpq, Cpq — некоторые константы,

 

 

 

 

г —текущий радиус.

 

 

 

 

 

 

 

 

Считая давление

Р0 постоянным внутри однородного по плотности шара, из

(579) согласно (567) находим:

 

 

 

 

 

 

 

 

ь »

+ f (*00

-

 

2*® >

=

 

+ СРР>

(581)

где Срр— некоторые константы.

 

 

 

 

 

 

Аналогично при

постоянной

 

плотности

р 0 в пренебрежении

временной

производной в даламбертиане решением (575) будет:

 

 

*00

- \(h m -

f h n )

= ^

+ Ст ,

(582)

 

 

1

 

5=1

 

 

 

 

где CQQ— некоторая константа.

 

 

 

 

 

 

Просуммируем решение (581) по индексу

р

от 1 до 3 и подставим в (582),

3

 

3

 

 

 

 

 

 

 

учитывая, что ^ h pp =

2 *ss:

 

 

 

 

 

 

 

р = 1

5=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8 яу Р 0гг

 

 

 

 

 

 

/7= 1

 

 

 

,

_

Алур„гг ^ 4 я у /,„/^

х

 

 

*0. -

 

+ — —

+

2

 

444 §48.4. Метрика внутри однородного шара. Квазигалилеев подход

Из (581) находим hpp :

_ 4л урУ _ 4луРУ

+ см_ I Vr

№ ~ Зсг

Зс4

рр 7 1 Z d W

Произведем дальнейшие упрощения. Из решения для hpp видно, что эти величи­ ны совпадают с точностью до константы, поэтому в симметричном случае можно по­ ложить: hpp = Л„ = Л22 = Л33,

Си = ^22 = С33 = Срр = Т ^ р р р '

1

Используем далее решение Шварцильда (566) для метрики вокруг гравитирую­ щей массы в сферических координатах:

dsl{r) = gw(r)c2dt2 + gu{r)dr2 + g22(r)(dQ2 + sin2f?dp2),

SooW = 1 - 7 ^ ' *uW = ---------= - г2-

1 " 7 7 “

 

В гармонических координатах в общем виде интервал имеет вид:

 

ds2 = gifc dxl dxk,

(583)

причем в статическом случае при нулевых скоростях движения вещества основной вклад в интервал согласно (575) — (578) должны вносить только диагональные члены метрического тензора. Свяжем компоненты метрического тензора gik(r) решения Шварцильда в сферических координатах с диагональными компонентами метрического тензора в (583) в гармонических координатах с помощью преобразования пространственных координат:

х = rsinQcosp, у = rsinQsinp, z = rcosQ.

Подставляя в (583) дифференциалы dx, dy, dz, выраженные в сферических коор­ динатах, и считая, что#п = g22 =gn = gpp , находим равенство интервала Шварцильда

и (583) при следующих условиях:

 

 

 

 

 

*ooW = Як»

*ц(г) = gpp

 

Так как мы приняли соотношение gjk = rjik

+ hik , откуда

= 1 + Aw ,

gpp = - l + hpp , то для границы сферического тела при радиусе R,

когда метрика

Шварцильда стыкуется с внутренней метрикой тела, можно записать:

hJR)

 

2 уМ

 

 

 

Rc2

 

 

 

 

 

 

hpp{R) ” 1 4* gpp =

1 -

 

1

2 уМ

 

 

2 yM

Rc2

 

 

 

 

 

 

 

 

Rc2

 

 

Подставляя hw{R) = hpp(R) в решения для

hm и hpp при r= R

и учитывая, что

М = р ^ я Л 3, находим константы и полное решение внутри сферической массы:

8 яу Р 0Д2

г

 

вуМ

 

3

с4

00

Rc2 '

 

§48.5. Определение вклада энергии гравитационного поля в метрику шара

445

L

_

4 я у р ,,г 2

4 яуР ^Я 1 ~ Г2)

ЪуМ

 

П00

 

1 2

4

п 2

 

 

 

Зс'

с

R e

 

,

_ А л у р У

А лyP0(R* - г2)

ЗуМ

 

рр ~

Зс2 +

3 ?

Т Г

( )

С учетом (568) компоненты йм и hpp внутри тела оказываются функциями давления и скалярного гравитационного потенциала у>в:

 

ь

_

2ip„

 

АлуР0(Я2 -

г2)

 

 

00

~

с2

 

?

 

t

 

_

2f„

,

4*у/>02 - г2)

Подставляя эти значения

 

и

Ирр в (575), в статическом случае приходим к

уравнению Пуассона для гравитационного потенциала внутри тела:

 

 

 

 

ДV'a =4яур0-

 

Из (576) видно, что если обозначить:

 

 

V

и

= -

~Dp,

(Ро + V '

= JP>

 

 

 

С

 

с

 

то (576) будет эквивалентно уравнению (508) для векторного гравитационного потенциала D и плотности тока массы / с учетом вклада массы-энергии от давле­ ния PQ. Поэтому величины hQp следует считать соответствующими функциями ком­ понент векторного потенциала Dp. Из вышеизложенного вытекает, что метрику непосредственно нельзя отождествлять с 1равитационным полем — метрика есть сложная функция от состояния движения масс (либо функция соответствующих гра­ витационных потенциалов с учетом состояния вещества). Из (584) находим в центре сферического однородного тела при г = 0 :

*о,(0) = ~ 4* yR2(P» /

, /т

АяуЯ\р, -

/>/ (З с 2) ]

(585)

V<°>

— г

 

Если масса тела М и его радиус R или плотность вещества р 0 и давление Р0 стремятся одновременно к нулю, то hw и Ирр также стремятся к нулю и метрика становится галилеевой, а пространство-время плоским.

§ 48 .5 . Определение вклада энергии гравитационного поля в метрику однородного шара

Поскольку гравитационное поле обладает плотностью энергии (511) и описыва­

ется тензором плотности энергии-импульса Uik (533), то это необходимо учитывать в уравнении (573) в тензоре материи Т[к. Как будет показано далее, это приведет к поправкам в следующем порядке малости по отношению к решению Шварцильда. Ре­ шение задачи начнем аналогично [117J. Используя сферические координаты JC° = ct, х1 = г , х2= Q , JC3 = , будем искать неизвестные метрические коэффициенты gik в статическом интервале риманова пространства-времени для тела, обладающего сферической симметрией:

446 §48.5. Определение вклада энергии гравитационного поля в метрику шара

= gwc2dt2 + g tldr2 + gndQ2 + gnd<p2.

Обозначим gm = В , gH= — К , ga = — E , g}i = — E sin2 Q и будем считать, что величины By Ку £ зависят только от радиуса г. Метрический тензор и его детерминант выглядят так:

 

0

0

0

'

'1/ £

0

0

0

'

0

- л :

0

0

 

0

- 1

0

0

 

 

 

 

 

,

 

 

 

 

 

0

0

- Е

0

 

0

0

- 1

а

 

,0

0

0

- £ sin 2G;

, 0

0

0

- \/(Е sin2Q)y

 

 

 

g = detgik = -

BKE2sin2Q.

 

 

 

Запишем (573) для тензоров со смешанными индексами:

 

 

 

 

 

 

V

-1<5?Д

= ^

г Д

 

 

(586)

 

 

 

 

2

с

 

 

 

 

гдеR k = Rin gnk — тензор Риччи со смешанными индексами, его можно найти, зная Rin из (571) и метрический тензор gnk,

6k — символ Кронекера, то есть тензор, у которого главные диагональные

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

По формулам (570), (571) и (572) находим ненулевые коэффициенты Кристоффеля, компоненты тензора Риччи и скалярную кривизну 2?, учитывая

равенство нулю производных по времени и по углам Q, <р:

 

 

г0

= г>0

J L

г. = J L

 

г > = A L

г ' = - А 1

 

■ *01

1 10

 

 

2В ’ 00

 

2К ’

 

11

2 К '

22

2 К ’

 

 

 

 

 

 

 

_

E'sin2Q

г 2

г2

г*3

_ гЗ

Е'

 

 

33

2 К

1

12 *“

1

21

*" 1

13

1

31

2Е ’

 

 

Гh = - sinficosG,

 

Г%} = Г пг

 

=

ctgQ,

 

 

 

 

В"

 

В Т

 

 

(В')2

 

t

 

В'Е'

 

 

 

я ; = 2 ВК

 

4 В К 1

4 ВгК

 

2 В К Е ’

 

 

В"

В'К'

 

(ВГ

 

+ А 1 . .

(E 'f

К 'Е 1

 

 

2ВК

4ВК 2

 

4В1 К + К Е

 

 

2 К Е 1

2К гЕ

 

 

д 2 _ д ? _

Е"

 

К'Е'

В'Е'

1

 

 

2

3

2 £ £

4 £ 2£

4 5 £ £

£ ’

 

5 " 5 ' Г

( £ ') 2

+

2 £ "

 

 

(Е')2

 

* ' £ '

В'Е'

2

ВК

2ВК1

2ВгК

 

КЕ

 

 

2КЕг

 

КгЕ + ВКЕ

Е'

здесь один штрих означает первую производную по координате г, а два штриха - вторую производную.

Тензор 7}k в (586) состоит из двух частей: T-k = U k + к,

§48.5. Определение вклада энергии гравитационного поля в метрику шара 447

где Uiк —тензор плотности энергии-импульса гравитационного поля,

t k — аналогичный тензор вещества.

Используя определение тензора Ulk из (533), запишем его в смешанных коорди­ натах:

и -: = &ти пк = 4л у (~Ф1тФтк + ^ Ф 1тФш \

здесь Фш = ds Dm - дт Ds — тензор гравитационного поля (529),

Dm — ковариантный 4-векгор потенциала.

В сферических координатах скалярный потенциал поля f зависит только от радиуса, а для статических масс векторный потенциал D равен нулю. В результате за пределами тела имеем:

Di = (3£ , - D) = ( -

, 0 , 0 , 0 ) ,

СГ С

Ф = - ф

= - У—

0 ns Ш8* ф

g* ф 01 = - ф ш = уМ

*10

г2с

9

ВК г2с

 

 

уМ 2

уМ 2

 

 

йлВКг*

и? = и 2

 

 

8 лВ К г4 9

где у — гравитационная постоянная,

М— гравитационная масса шара,

г— текущий радиус,

с — скорость света.

Внутри шара надо использовать потенциал грв из (568):

VB =

(г1 - з д 2),

 

3

где р 0 — однородная внутри шара плотность вещества, R — радиус шара.

Компоненты тензора гравитационного поля примут вид:

ф

= - ф

 

dipв _ _ 4 я у р 0Г

*01

*10

дг

3с

 

 

 

 

ф »'

= _ ф«» = 4?т

г

3сВК '

Гравитационный тензор плотности энергии-импульса внутри тела внешне похож на этот же тензор вне тела с той разницей, что нужно учитывать только массу вещества М{г) внутри текущего радиуса г :

UQ° - и х

У М \ Г )

2

_

уМ \г)

ВлВКг*

2

5

iitB K r* '

 

В качестве тензора вещества возьмем тензор /*

(537)

для несжимаемой жидкости

при L = 0 , тогда можно записать:

 

 

 

 

= gin*"* = (Ро +

 

-

W i ’

где р 0 и Р0 — соответственно плотность и давление в элементе покоящегося вещества,

448 §48.5. Определение вклада энергии гравитационного поля в метрику шара

iff — ковариантная 4-скорость.

 

 

В статическом случае вещество неподвижно,

и

если вектор положения

х 1 = (ct , г , Q , <р), то вектор смещения вещества dxl

= (сdt , 0 , 0 , 0 ) ,

dxк = gkj dx‘ =(Bcdt , 0 , 0

,

0),

интервал:

2 = dx

dxk - Bc2dt2,

ds = 4 B cdt, dt0 = — = -/fi A,

 

fc

 

 

 

 

 

 

 

 

 

c

4-скорость:

и' -

= (-£= , 0 , 0 , 0 ) ,

 

uf = {4B c

, 0 , 0 , 0 ) ,

 

dt0

-JB

 

 

 

 

 

 

 

 

 

ненулевые компонента тензора вещества:

/0° = р 0 с2,

 

 

 

 

 

= , 22

= / 33

=

 

- р 0.

 

 

Уравнения (586) принимают вид:

 

 

 

 

 

 

 

 

Е"

(Е')2

 

К'Е'

1

_

 

у2М \г)

8 я у р 0

 

КЕ

4 К Е г

1К2Е

 

Е

 

 

В К г4с4

с2

 

(Е')г

В'Е'

 

| 1

_

 

у2М \г)

%луР„

 

~ 4 К Е 2 2ВК Е

+ Е

 

 

В К г4с*

 

с4

_ Л . _ + в 'к ' +

W

_

 

£ "

+ ( £ '>2

+

_

 

25АГ

4BK 2

4ВгК

 

2К Е

4 К Е 1

4 К 2Е

 

 

В'Е'

_ у2М \г)

 

8 луР й

 

(589)

 

 

4ВК Е

В К г4с4

 

с4

 

 

 

 

 

 

Вычитая уравнение (588) из (587), находим:

 

 

 

 

 

Е'

d \

( B K E S

 

8 лу(р„с2 + Pt)

 

2КЕ dr[

\(E ')2 j

 

 

 

4

 

 

 

 

 

 

 

 

 

Если рассматривать область вне гравитирующего тела, где плотность вещества р„ = 0 и давление Р0= 0, то получится:

2

In ВКЕ = const, ВК = С, (Е')г

Постоянную С, можно найти из условия для В, К, Е на бесконечности, где метрика должна быть псевдоевклидовой:

В(оо) = 1, К(*>) = 1, Е(оо) = г\ отсюда С,= 1.

Если положить В К - 1 и г2 = Е только для правой части уравнений (587), (588), 589), то можно быть уверенным, что вносимая тем самым ошибка будет незначитель­ на, и уравнения (587), (588) можно переписать так:

!^

V

1 _

y2M \r)

8 я у р 0

1 d

1

r d 4 I dr

К

Е

c4E 2

c2

dr

 

 

 

 

1

d-jE d Ы в Щ

+ 1

= - tlfW L -

(591)

Y j E

dr Tr 1

1 E

C 4E 2

C4

§48.5. Определение вклада энергии гравитационного поля в метрику шара 449

Решим эти уравнения вначале вне вещества, затем внутри гравитирующего тела и сравним оба решения на границе сферического тела при его радиусе R. Интеграл от (590) вне вещества при р 0 = 0 и M(r) = М дает:

к

= •

VI

(592)

угМ 2

 

V I -

dr ) '

 

с4VI +

с г

где С2 — некоторая константа.

Решение уравнения (591) при Р0 = 0 и величине К из (592) приводит к

следующему значению В:

 

 

V I -

угМг + с 2

 

В = с 3

с*VI

(593)

VI

Из условия на бесконечности следует, что ВК= 1 и постоянная С3 = 1. Нетрудно проверить, что решение уравнения (589) вне вещества (когда M(r) = М) сводится к уже полученным решениям (592) и (593) при сделанных выше предположениях.

Однако строгое решение уравнений (587), (588), (589) без упрощающих условий ВК = 1 и г2 = Е для правых частей этих уравнений могло бы в принципе дать решение для Ву К, Е как функций от радиуса г, в то время как по (592), (593) К и В являются функциями от Е, а зависимость Е от радиуса г не определена.

Перейдем теперь к решению уравнения (590) внутри гравитирующего тела. Для упрощения интегрирования предположим, что гравитационную массу М{г) можно представить в виде:

М(г) =

4 * p „(V I)J

»

 

.

 

где р 0 — некоторая средняя плотность вещества тела.

 

В результате получим:

 

 

 

К(Е) = ______________VI_____________

(594)

Г= _ 2 уМ(Е)

16л У р о (У1)5

 

с2

45с4

4

 

гдъМ (Е) = f 4jrp0 (VZ)2*/VZ — некоторая текущая масса гравитирующего тела, До)

С4 — некоторая постоянная.

Сравнивая (592) и (594) при значении лГЕ, соответствующем радиусу тела Л, "можно оценить С2 :

с _ 2 уМ а 6 угМ г

C l ~ ~ ^ ~ + ~ S ? R +C i’

(595)

 

гдеЛ/ 0 = f 4 n p 0( jE ) 2d->[E—масса гравитирующего тела при r - R

(596)

Д (0)

450 §48.5. Определение вклада энергии гравитационного поля в метрику шара

и при J E = -JE(R), причем в (595) в знаменателе второго члена справа вместо

•JE(R) подставлено R}что не вносит большой ошибки.

При движении внутрь однородной по плотности сферической массы по направ­ лению к ее центру гравитационный потенциал, как известно, стремится к константе и согласно (568) зависит только от текущего радиуса г :

_

2 л у р йгг

_ 3уМ

W

3

2 R '

здесь M,R — гравитационная масса и радиус всего тела.

Воспользуемся результатами предыдущего раздела о метрике в центре сфериче­

ского тела, где согласно (585) получается:

 

4лгуД2[/з0 -

/(Зс2)]

h№(0) =

'С"

Так какв данном случае £,,(г) = gpp или - К - — 1 +Ирр , то имеем:

4яуЯ2\р0 - Р0/(Зс2)1

К(г = 0) = 1 - ЛдДО) = 1 + -

- 1 с2

К

(597)

здесь £п(г) — компонента метрического тензора в сферических координатах, gpp — компонента метрического тензора в гармонических координатах.

В конце этого раздела будет показано, что 4 е = г + const,

jE (r = 0) = JE(Q) ~

и отсюда

= 1.

 

с

аг

Тогда из (594) следует, что в пренебрежении степенями малой величины ~JE(0)

величина К(Е) при г = 0 равна:

 

 

 

 

1

-

^4

 

 

 

+ с 4

 

4т'

 

 

 

1

+

 

 

 

 

 

 

 

 

Сравнивая данное значение с (597), находим С4 через ^

 

( 0):

 

_

i7tYR24 m [ p 0 - P 0 /(3 C2)]

.

(598)

ц

-

 

 

 

-2

 

С учетом значения -/Е(О) постоянная С4 оказывается величиной второго поряд-

1алости по отношению к квадрату скорости света. Интегрируем (591) внутри вещества с учетом (594):

г

 

\6ягу2р 0Е2 г

%яуР0Е

d[\n( В Щ _

1 +

 

'

+ ~

7 ~

rfVF

r= 2

 

 

 

= F{E),

уМ{Е) 16л2 у 'р Ц Л у

 

VЕ ------^

+ ------^

------ + С<

 

 

с-

 

45с"

 

или В(Е) =

С, exp[f F (E)d-fE\

Преобразуем неопределенный интеграл от Е(Е):

§48.5. Определение вклада энергии гравитационного поля в метрику шара 451

f F ( E ) d j E

(

с2

45с4

*]

= /

 

16л2у2р 1 ( Л ) 5 + С<

 

4 1 - 2уЩ Е)

 

 

 

45с4

 

 

 

£ЧРо + PJc2)d4I

 

 

Г F

2 у * ( Я ) ,

К>яг у гр 1 { 4 Ё У

, „

 

 

 

45с4

 

 

'J J _

2уМ(£)

16я2у2р 2(л/£)5

+ С<

*<*) = С5

с2

45с4

 

 

41

 

 

 

 

 

8 * у ^ _________ Е(р0 + P J c 2) d 4 I _________

 

 

 

 

(599)

*2

щ » 4 Ё

- 2уМ2(Е) + 16” У р

) 5 + с4

Под знаком экспоненты стоит выражение, пропорциональное отношению гравитационной энергии к энергии покоя вещества внутри текущего радиуса ^E (r).

Из (593) и (599) с учетом (595) следует, что на границе тела при г = R произведение экспоненты на постоянную С5 должно равняться единице. Постоянную С5 можно оценить с помощью данных предыдущего раздела, где компонента gw(r) метрического тензора в сферических координатах связывалась с компонентой метрического тензора в гармонических координатах:

£ооМ = £оо или в обозначениях данного раздела ^ (г ) = В: В = 1 + hw , где hmсогласно (585). Тогда при г = 0 находим:

В ( г = Щ = \ -

(600)

 

С

При г = 0 экспонента в (599) обращается в единицу и приблизительно получаем:

В(г = 0) = С54 т +с< J m

Сравнивая данное выражение с (600) с учетом (598), находим, что С5 < 1:

_ 4я у Р г(р0 + Рр/с2)

^С*

с%

4жуЛ2[р 0 - Р0/(Зс2)]

1

/»2

Уравнение (589) внутри вещества сводится к релятивистскому уравнению гидро­ статики, связывающему текущую массу М(Е), давление и плотность вещества в зави­ симости от координаты 4 Ё . Если бы зависимость давления и плотности вещества от радиуса г была известной, то можно было бы найти величину Е как функцию от г внутри вещества. За пределами гравитирующего тела давление и плотность вещества равны нулю, а зависимость Е от г оказывается точно не определенной в связи со сложностью совместного решения уравнений (587), (588), (589). Обычно