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

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

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

402

§48.1. Уравнения гравитации в пространстве Минковского

следует ожидать «схлопывания» всех нейтронных звезд галактического ядра в одну черную дыру, как это часто предполагают. Точно также массивное атомное ядро из множества нуклонов не сливается в одну ядерную каплю. Активное ядро из нейтрон­ ных звезд объясняет не только энергетику квазаров, но и возникновение регулярного магнитного поля за счет выстраивания магнитных моментов отдельных звезд. Образу­ ющаяся осевая конфигурация магнитного поля позволяет коллимировать потоки ре­ лятивистской плазмы, что часто и наблюдается в виде джетов и выбросов из ядер активных галактик.

В заключение рассмотрим переменность нетеплового излучения во всех диапазо­ нах частот от сейфертовских галактик и квазаров. Ядра большинства галактик неод­ нородны и могут состоять из нескольких ядрышек или сгущений звезд. В таком случае при их вращении будут наблюдаться долговременные и квазипериодические изменения блеска. Если размер ядра R ~ 1013 метра, скорость вращения звезд поряд­ ка С = 220 км/с, то период вращения будет Т~ 2 nR/C = 9 лет. И действительно, наи­ большие характерные периоды колебаний блеска активных галактик составляют несколько лет (у квазара ЗС 273 — порядка 13 лет). В ядре нашей Галактики размеры точечных радиоисточников менее 1012 метров (например, у IRS 16), в то же время наблюдаются колебания длинноволнового гамма-излучения с периодом около одно­ го месяца. В течении нескольких месяцев оптически изменяется на 1 — Зт TV-галактикаХ Comae. Кромедолгопериодической переменности, связанной с враще­ нием ядер галактик, иногда проявляется и быстрая изменчивость в течении 1 — 2 суток. Если разделить размер ядра галактики R = 1013 метра на период 24 часа, найдем скорость распространения возмущения вдоль ядра:

—--------- 10е м/с, что почти равно скорости света. 24часа

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

Сравнительно недавно внутри большой эллиптической галактики М87 нашли ядро с массой 3• 109 Мси скоростью вращения газового диска вокруг него до 500 км/с [4]. Газовый диск имеет спиральную структуру, как в спиральных галактиках, кроме этого наблюдается мощная газовая струя. Изучение шарового скопления M l5, расположенного в созвездии Пегаса на расстоянии 42000 световых лет, показало, что в его ядре нет черной дыры, как полагали ранее, но в самом центре в радиусе 0,4 светового года находятся тысячи тесно расположенных звезд [38].

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

§48. Теория тяготения

§48.1. Уравнения гравитации в пространстве Минковского

Теория всемирного тяготения Ньютона, подобно электростатике, не является лоренц-инвариантной и потому не может быть согласована со специальной теорией относительности. В свое время объединение электричества и магнетизма уравнения­ ми Максвелла позволило осознать новое понятие — электромагнитное поле, законы которого автоматически удовлетворяют преобразованиям Лоренца. В § 45 шла речь о противоположности и взаимодополнительности гравитационных и электромагнит­ ных полей, так что в самом общем случае их надо было бы описывать как единый

404

§48.1. Уравнения гравитации в пространстве Минковского

 

divG = -

4я >7>,

divfl = 0,

 

(504)

 

rotC =

dQ

c2rotfl = - 4 n y J

+

dG

 

 

d t

 

 

 

dt 9

 

/ = **Le =

P »V

p =

PQ

 

dtds

-

v 2/c2

J l

~

^ 2/c 2

здесь/э — плотность движущегося вещества (количество массы в единичном объеме),

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

/— плотность тока массы в направлении единичного вектора е через сечение ds9 перпендикулярное вектору е,

сскорость распространения гравитации, которую далее будем считать равной скорости света,

у— гравитационная постоянная.

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

Fj? = mG + m[V х 0],

(505)

где V—скорость движения тела, т — масса покоя.

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

Векторы гравитационного ускорения и кручения можно выразить через

скалярный гр и векторный D потенциалы:

 

 

G = - gradip -

Q = rotD.

(506)

Потенциалы в свою очередь определяются плотностью вещества р и плотностью

тока массы / через объемные интегралы:

 

 

¥ 0 ,0

= -

 

(507)

D M

= - С

л п

 

здесь подразумевается, что потенциалы в точке 1 в момент времени

t нужно

находить путем интегрирования распределения масс (массовых токов) в области v2 в более ранний момент tl - t — Rn /с ; Rn — текущее расстояние между точкой 1

и элементами объема dv2 в момент/', с — скорость света.

и D :

Подставляя (506) в (504), получим волновые уравнения для потенциалов

~ А л у р ’

(508)

1 дгР _

4 я у ,

 

~ dt1 ~

~

 

при условии калибровки: divD = -

с2 dt

§48.1. Уравнения гравитации в пространстве Минковского

405

Здесь Л = —-

+ — 2 + — 2 “ оператор Лапласа.

 

дх

ду

dz

 

д2

 

Используя оператор Д’Аламбера: □

1

 

= А — j — ,

 

можно переписать (508) в виде:

 

 

 

 

 

 

= 4я у р ,

 

(509)

 

 

U

D ^ J .

 

 

 

 

 

С

 

 

Для завершения картины необходимо добавить следующие соотношения:

Закон сохранения массы:

d i v / = - — .

 

(510)

Плотность энергии гравитационного поля:

и = -----— (Cr2 + c2Q2).

(511)

 

 

 

8л у

 

Вектор плотности потока энергии поля: Sr

-----------[G х О ].

(512)

 

 

 

4 лу

 

Плотность импульса поля:

g =

Sr .

 

(513)

Размерности используемых величин таковы:

кручение [£?] = с-1, гравитационное ускорение [<7] = м с"2, скалярный потенци­ ал [у] = м2*с-2 , векторный потенциал [D] = м*с“‘, плотность массы |р] = кг*м"3,

плотность тока массы [/] = кг-с^-м"2, плотность энергии [и] = кг с"2-м_| , плот­

ность потока энергии [ ^ ] = кг • с"3, плотность импульса [g] = к г-с'1' м"2

Далее предполагается, что масса и ее плотность однополярны и положительны, а гравитационное поле состоит из двух компонент — ускорения G и кручения Q. Заметим, что кручение имеет размерность как у угловой скорости, а векторный потенциал — как у линейной скорости.

Используя (507) и (506), кручение от бесконечного прямолинейного тока массы можно записать через векторное произведение:

2R[R dt

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

R ■— кратчайший вектор, соединящий ток массы с точкой, где определяется кручение,

dM

------- вектор тока массы. dt

В центре тока массы кольцеобразной формы с радиусом R кручение равно:

Q = ^ п У dM c2R dt

По аналогии с магнитным моментом введем гравитационный момент Рс как произведение массового тока на обхватываемую им площадь:

Ш -nR2

dt

406

§48.1. Уравнения гравитации в пространстве Минковского

Для расчета кручения в центре вращающегося кольца имеем:

М _

М

 

Д

_ M VR _ I

dt

Г

Г

’ с

2

2 ’

 

 

 

О = -

У

 

 

 

 

 

с2Л3Л

 

 

 

 

 

 

 

 

здесь М масса кольца,

 

 

 

 

 

 

Г — период вращения кольца,

 

 

 

V— линейная скорость вращения кольца,

 

 

/ — спин (момент импульса) кольца.

 

 

 

На больших расстояниях

г

от вращающегося кольца кручение представляется

дипольной формулой:

 

 

 

 

 

 

 

о

 

= J L ( L

г5

 

 

 

 

 

2с2 V

 

 

Для дальнейшего изложения напомним некоторые положения алгебры четырех­ мерных векторов, используемые в теории относительности (смотри например [113]). Вектор местоположения частицы или точки в системе отсчета с временем г, координатами х, у, г и трехмерным радиус-вектором г записывается так:

х * = (ctyx>y, z) = (с/, г) то есть 3-вектор г является составной частью

4-векгорах1, задавая в нем три последние компоненты х, у,

z.

Четырехмерная скорость частицы, движущейся с трехмерной скоростью V :

ni

_ (

с

 

V

г).

 

~

VI - Уг/ ' г '

 

 

V1 - уг/ с

 

Четырехмерный импульс:

 

 

 

 

 

р = ти

= (

т с

 

тУ

г).

- У2/с 2

’ Vl

 

 

 

VI

- ^ г/с :

здесь т — масса покоя частицы,

 

 

 

 

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

 

 

 

 

 

 

Четырехмерная сила:

 

 

 

 

 

 

 

 

dE

 

 

 

F‘ = (— rJL --------

Vl -

г).

(515)

 

 

cVl -

У2 /с 2

У2/с

 

где £ — полная энергия частицы,

F — обычная сила в трехмерном пространстве,

V— скорость движения частицы, на которую действует сила.

Индексы /, к, т, п и другие латинские буквы пробегают значения 0, 1, 2, 3, первые компоненты четырехмерных векторов всегда скалярные величины, а осталь­ ные три компоненты являются компонентами трехмерных векторов. Если символ вектора (или тензора) находится вверху, то этот вектор называется контравариантным, нижний символ у вектора обозначает ковариантный вектор. При одинаковых индексах, встречающихся одновременно в верхних и нижних символах в произведе­ ниях векторов или тензоров, будет подразумеваться необходимость суммирования со­ ответствующих компонент. Переход от ковариантных индексов к контравариантным и наоборот осуществляется с помощью метрического тензора gik:

§48.1. Уравнения гравитации в пространстве Минковского

407

A i = gikA \ A, - g lkAk.

 

В псевдоевклидовом пространстве (пространство-время Минковского) с учетом использумой нами системы так называемых действительных координат метрический тензор сводится к тензору rjik:

 

 

fl

0

0

0 '

 

 

 

0

-1

0

0

 

_

ik _

 

 

 

(516)

Vile

v

-

0

 

 

 

 

0

-1

0

 

 

 

,0

0

0

- h

 

Для тензора второго ранга,

если рассматривать его

как матрицу, индексы

i, к = 0, 1, 2, 3 соответствуют значениям индексов матриц

1, 2, 3,4; первый индекс

слева означает номер строки, другой правый индекс — номер столбца соответствующей матрицы.

Тогда если 4-вектор Ак =

°, А 1, А 2, А 3), то

А,

= {А0 ,Л , ,Л 2

, А ,)= 4ikAk = ( A \ - Л ' , - А 2, - А }), и

 

А0 = А °, А { = А \ А2 = —А 2, А 3 = —Л 3,

или i-я

компонента вектора Л,, есть скалярное произведение /-й строки

тензора г\1к на вектор — строку Ак

Скалярное произведение двух векторов имеет вид:

 

А'В = А ХВ 1 = А 0В° + Л ^ 1 + А 2В 7 + А 3В г =

 

= ПхкАк В1 = Л°^° - А1В1 - А2В2 - АЪВ \

Длина любого вектора является инвариантной величиной при смене системы от­ счета, это же относится и к квадрату его длины, который можно найти из скалярного произведения вектора самого на себя. Для вектора местоположения имеем:

х(х* = (ct)2- x 2- y 2- z2= ( c t f - r \

Заметим, что четырехмерное скалярное произведение в отличие от трехмерного скалярного произведения имеет знаки минус при перемножении компонент 1, 2, 3. В системе отсчета с временем /' и трехмерным радиус-вектором г' длина вектора не должна измениться:

х] = XfX1 = (ctУ - г'2 = s '2 = (с/)2 - г2 = s2.

Величина s носит название интервал и для каждой точки является инвариантом. Найдем длину вектора сдвига (перемещения) dxl = (cdt, dr):

dx] = ds2 = c2(dt)2 - {dr)2 = c V 0 2[l -

.

L

cdt

Учтем теперь, что — = V:

408

§48.1. Уравнения гравитации в пространстве Минковского

 

ds = c d tj 1 -

V2/с 2

 

Величины ds и ds/c = dt0 являются

инвариантами, тогда и величина

<//0

= -у/1 - V2/с 2 dt является инвариантом при преобразованиях системы отсчета.

Трехмерная скорость обычно находится дифференцированием вектора положения г :

 

г = * .

 

 

dt

 

Для четырехмерной скорости имеем аналогично:

 

и = —

cdt

dr,

= (-

) =

Л .

Vl - Г 7 с 2Л ’ Vl -

^ 2/с 2 dt

СV

~’ Vi - к 2/с 2>

Выражения для трехмерных и четырехмерных сил через производные от импульсов по времени:

 

 

F

_

Ф

f

I _ dp‘

_ ф

и

)

_

_d_(m dx_^

 

(517)

 

 

 

 

dt

 

dtn

dt,

 

 

 

Л ,

Л .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'*0

 

 

 

 

 

 

 

 

 

 

 

 

 

d{

me

 

 

 

mV

 

 

 

 

 

 

 

 

V _

 

 

 

 

 

 

 

 

 

 

 

 

pi

-

-Jl

-

КУсг)

 

^ 1

-

K2/c ;)

 

 

 

 

 

 

 

 

Лд/l

-

K2/c 2

Лд/ l “

KV cJ

 

 

Из сравнения с (515) для векторной компоненты силы

 

получается определе­

ние трехмерной силы в теории относительности:

 

 

 

 

 

 

 

 

 

 

 

mV

,

vrd ,

 

т

ч

.

 

 

 

 

Л

 

 

-

 

г) = * 4 (

 

- ^ 2/с г)

+ V1

~

 

Л

 

 

 

V1 /с1'

"

 

 

 

Сопоставляя скалярные компоненты (515) и (517), находим:

 

 

 

 

 

 

 

(/ ,

т с 2

.

_

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-

KJ/c 2>

~

dt

 

 

 

 

 

где£ =

т с

 

 

— полная энергия частицы массы m с учетом энергии

Vl - Кг/с :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

покоя.

 

 

 

 

 

 

 

 

 

 

 

Ускорение —

тела с переменной массой (без учета импульса теряющейся или

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

привносимой массы) вычисляется по следующей формуле:

 

 

 

 

 

 

 

т

 

dV

= (1 -

V \„

 

п77ГГл

dm

V

(518)

 

Vl

-

V '/c1

dt

~ ) F

- V f -

Уг/с г

^

 

 

 

 

 

 

 

 

 

 

dt

 

 

Это соотношение доказывается следующим путем. С одной стороны можно записать:

С другой стороны имеем:

 

§48.1. Уравнения гравитации в пространстве Минковского

409

P V

m V2

тс

- К2/с 2,

=

- m c ^ l - К2/с 3 = £ -

 

•Jl - V2/c 2

~ у1/ с3

 

 

 

</£

*//я

 

 

 

"л ‘

Сравнивая два выражения для —(p'V), подставляя — для случая переменной

dt

dt

 

массы и сокращая результат на скорость

V, приходим к (518). Если масса т

не

просто изменяется со временем, но и переносит за время dt импульс (dm s)

с

абсолютной скоростью s, то (518) необходимо уточнить:

 

dV = (1 - -V)#- + Vi - K7 cJ dm y - V ) .

(5180

д/l - v 2/c 2 Л

 

В обычном случае, например при взлете ракеты, s > F, —- < 0 (масса ракеты

 

 

 

dt

уменьшается), и ускорение — направлено против скорости s.

 

dt

 

 

Если масса т постоянна во

времени,

то тогда dE/dt = V F и (515), (517)

принимают вид:

 

 

 

F '

F V

(519)

= (

к 2/с г

 

cV 1 -

V1 - К2/с 2

При переходе из одной инерциальной системы отсчета в другую все четырехмер­ ные векторы преобразуются с помощью соотношений типа (386) и (388) из теории относительности. Если в неподвижной системе отсчета К координаты точки опреде­ ляются вектором х 1 = ( с /, г), то в системе отсчета К\ движущейся со скоростью К0

вдоль оси X относительно К,

вектор этой же точки равен

х’1 = (с/', г% причем

выполняются соотношения между координатами:

 

t , _ / - У р Х / с 2

#

х - К t

Z = Z.

 

 

у = у ,

V1 - П 7 с 2 ’

V1 - >о2/с 2 '

 

Вообще при частном преобразовании Лоренца в движущуюся систему отсчета компоненты любого контравариантного вектора будут выглядеть так:

АЛ = А0 - У0 А 1/с

г = А 1 -

К А 0/С

Ап = А \ А'3 = А 3. (520)

4 ^ У о / с 2

л/l -

К /С 2

 

Для ковариантных векторов с учетом изменения знака пространственных компо­

нент получается:

 

 

 

Ао + V ,A J c

А\ =

4*

AQ/ C

Л'„ =

VI -

—^2 » -^3 “ ^З*

Vi - у з/с г

*?/«* ’

Применяя преобразования типа (520) к скалярной компоненте четырехмерной скорости и 1 как к скалярной (временной) компоненте вектора х \ можно найти преобразование лоренцева фактора:

410

§48.1. Уравнения гравитации в пространстве Минковского

 

 

 

с

УрУ*

 

с

Vl -

Уг/сг

c j l -

Уг/с 1

 

Vl - Vй/с г

VI -

К /с 2

или:

 

 

 

Vl -

V'2fc 2

V> - к2/сг Vi -

r02/c2

 

1 -

 

(521)

 

 

 

^ , / c 2

 

здесь V и У' — скорости частицы в системах отсчета К к К’ соответственно,

 

VQ— скорость движения системы отсчета К' относительной,

 

Vx проекция скорости частицы

V на ось X системы отсчета К.

Преобразование х-компоненты четырехмерной скорости и 1 проводится также, как для х-компоненты вектора х 1:

Уж__________%____

rx

_ Vi - v \/c 2

V1 - ^ / с 2 _

V 1 -

vx -

к,_______

V l- Г */с’

V 1

-^ 0/с 2

^

2 / c V*2 Vl1 ‘ -

Подставляя в это выражение (521), получаем:

У х - у ,

^ =

1 - w c 2’

аэто есть релятивистская формула сложения скоростей, обратная по отношению

к(392). Из (520) следует, в частности, что у = у ' для вектора х 1, также будет и для

у-компоненты вектора и1, аналогично для z-компонент:

V l-

К

-

у г

г 2/с2 V 1

-У2/с2

Используя (521), приходим к

следующему:

 

г

 

_

/с*

 

r

I - У * У х / с 2 '

Подобно выглядит формула и для z-компоненты скорости И' в системе отсчета К*(сравни с (392)). В (393) и (394) фактически приведен пример обратных преобра­ зований типа (387) и (389) от четырехмерного вектора импульса р’1к вектору р1

В качестве еще одного примера приведем выражение для четырехмерного

волнового вектора произвольной волны:

 

Kt =(й),сХЛ К1 =(й>,-сЛ У ,

 

здесь а)— угловая частота колебаний,

 

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

 

К — волновой вектор, К - 2 п /1 , Я — длина волны.

 

Нетрудно убедиться непосредственной подстановкой, что квадрат вектора

а

также скалярное произведение К(х 1являются инвариантами при преобразованиях Лоренца:

К] = KjК' = 0 г - сгК \ К р = с (су/ - я - г ) .

§48.1. Уравнения гравитации в пространстве Минковского

411

Используя (520) для ковариантного вектора, получим преобразованные компоненты вектора Ki в другой инерциальной системе, двигающейся со скоростью К0 относительно оси X:

о) + VQKX

К' ~ Кх + Улш/сг К', = KY, К'г = Кг .

л/1 - К / с 1

V* - к 1/с2

здесь Кх , Ку, К2— проекции волнового вектора К на оси координат.

Пусть в исходной системе волна распространяется в плоскости XOYиз начала ко­ ординат, тогда Кх = К cos <p,KY= К sin <р, где у> - угол между осью X и вектором К. Если еще волна электромагнитная и монохроматическая, то ее скорость в вакууме равна скорости света и справедливо соотношение со = сК. Тогда угловая частота колебаний d будет равна:

ш(1 + — cosy>)

л/1 -У о /с 1

что целиком совпадает с (376) для выражения релятивистского эффекта Допплера. В исходной неподвижной системе отсчета также имеем: tg tp = KY/KX. Так как Ку = KY - A'siny?, то в движущейся системе отсчета для * л /2 и при К « с можно записать:

tg <р' = Б

 

 

K yjl

- V j/c1

_

siny-y/l

-

К02/с1

 

 

 

Кх +

К0« /с 2

 

 

cosp

+

V„

к*

 

 

 

 

 

 

 

 

tgy>

~ tg#)(l

).

 

 

 

 

 

К

 

 

 

1 +

 

 

 

 

CCOS(p

 

 

 

CCOS(p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если положить <p' = <p

 

+ 6<p, то так как

 

 

 

 

 

 

,

 

,

 

,

sin(y?'

-

<p)

dtp

 

tg

<P ~

 

tg ip = — ^

— —

cos2y>*

 

 

 

 

 

cosy) cosy?

то для отклонения волнового вектора в движущейся системе отсчета получаем зависимость от скорости и направления движения VQ:

 

 

л

vo

 

 

 

dip = —

-sin у?.

 

При = JT/ 2 находим:

tgy>' = -£-^1

-

V*/с2 ~

,

 

 

К

 

 

 

с

/

С

7t

К

Уп

д<р =

-

tp = arctg—--------= - arctg —

---------

 

 

Уо

2

с

с

Следовательно, при малых скоростях общей формулой для отклонения является

у

dtp ------ s in у), которая, как известно, описывает аберрацию света от звезд,

с

открытую Д. Брадлеем в 1725 году.

Теперь мы вновь обращаемся к гравитационным величинам. Четырехмерный

вектор плотности импульса будет иметь такой вид:

 

 

 

J

= /)„ « ' = (

ер*

Уро

=) =

(c p ,J ),

(522)

V1 - V1/с 1

т/1 - V1!с