Фейнмановские лекции по гравитации
Шрифт:
8
GT
=
G
=
R
–
1
2
g
R
.
(14.1.1)
Правая часть этого уравнения есть ”геометрическая” часть, здесь мы подставляем выражения для кривизны через компоненты метрического тензора. Если мы предполагаем статические, сферически симметричные решения, тогда элементы метрического тензора в точности определяются функциями (r) и (r) такими, что
(ds)^2
=
e
(dt)^2
–
e
(dr)^2
–
r^2
sin^2
(d)^2
+
(d)^2
.
(14.1.2)
Левая
G^1
=-
e
–
'
/
r
–
(e
–
– 1)
/
r^2
=-
8Gp
(14.1.3а)
G^2
=-
e
–
''
2
–
''
4
+
^2
4
+
('-')
2r
=-
8Gp
(14.1.3б)
G
=-
e
–
'
/
r
–
(e
–
– 1)
/
r^2
=
8G
(14.1.3в)
Модель, которую мы будем использовать, будет задаваться теми выражениями, которые мы подставим для давления p и плотности . Эти величины представляют давление и плотность, которые могли бы быть действительно измерены наблюдателем, стоящим в какой-либо выделенной точке. Мы не получим правильных решений до тех пор, пока мы не проследим за тем, чтобы наш физический тензор T удовлетворял законам сохранения. Для нашего случая сферической симметрии только радиальная компонента дивергенции тензора имеет значение; мы должны иметь
T^1
r
+
1
2
'
(T^1-T)
+
1
r
(T^1-T^2)
=
0
=-
1
2
'
(p+)-p'
,
(14.1.4)
что по сути дела утверждает то, что давления в радиальном направлении уравновешены, как это и должно быть в нашем статическом решении. Это уравнение (равенства нулю дивергенции) служит тому, чтобы исключить '. Далее мы получаем соотношение для того, чтобы исключить exp(-). Сначала мы перепишем G через новую функцию M(r), как показано в следующих соотношениях
G
=
1
r^2
d
dr
r(1-e
–
)
.
(14.1.5а)
Если мы положим
M(r)
=
1
2
r(1-e
–
)
,
e
–
=
1
–
2M(r)
r
,
(14.1.5б)
тогда
dM
dr
=
4r^2
G
.
(14.1.5в)
Оказывается,
1
–
2M
r
1
r
dp
dr
=-
(p-)
4Gp
+
M
r^3
.
(14.1.6)
Вместе с дифференциальным уравнением для M(r) и с уравнением состояния, связывающим величины p и , мы имеем систему связанных уравнений, которые могут быть в принципе разрешены для функций M(r), p и ; с подходящими граничными условиями они могли бы описывать сверхзвезду в приближении статического решения.
Какого рода уравнение состояния мы возьмём? Масса, образованная из 10 солнечных масс, является очень сильно разреженной, будучи размазанной по области с галактическими размерами; даже при температуре несколько единиц 10 градусов Кельвина, газовое давление является довольно низким. Тем не менее, оказывается, что плотность излучения, которая пропорционально T, даёт существенную часть энергии массы покоя нуклонной плотности. Мы получаем осмысленное приближение, пренебрегая газовым давлением по сравнению с давлением излучения; в том же самом духе, мы пренебрегаем небольшим увеличением массы нуклона, вызываемым их скоростями. В единицах энергии массы покоя нуклона мы имеем тогда, если s - плотность нуклонов, что
=
s+
,
(14.1.7а)
p
=
1
3
.
(14.1.7б)
Эти уравнения связывают p и , но мы всё ещё нуждаемся в том, чтобы в точности определить для того, чтобы иметь уравнение состояния. Мы делаем адиабатическое приближение, которое делает каждый, пытающийся иметь дело с такими проблемами, такое, что распределение температуры является тем же самым, как будто это есть величина, которая падает вместе с первоначально однородным распределением без всякого перемешивания или переноса энергии между различными областями. Если мы сжимаем вещество внутри ящика, все частоты вырастают на один и тот же множитель, обратно пропорциональный длине ящика. Так как энтропия является постоянной для адиабатического процесса, то температура должна увеличиваться таким же образом. Таким образом, плотность нуклонов пропорциональна кубу температуры, и плотность энергии излучения пропорциональна T. На языке температуры, измеренной в единицах 10 градусов, и энергии, в единицах массы покоя нуклона, имеем
=
aT
,
s
=
aT
.
(14.1.8)
Величина amn где mn есть масса нуклона, есть константа, имеющая значение 8.4 г/см^3; - параметр, связанный с не зависящей от радиуса энтропией на барион соотношением (энтропия на барион) = 4/(3). Эти результаты могут быть выведены также из общего условия для адиабатического сжатия, которое может быть выражено как
s^2
d(/s)