Курс теоретической астрофизики
Шрифт:
4
3
dI
dr
=-
H
,
(5.22)
где H — полный поток излучения в фотосфере и I — средняя полная интенсивность излучения.
Величина , определённая формулой (5.21), есть средний коэффициент поглощения. Вводя соответствующую ему оптическую глубину по формуле
=
r
dr
,
(5.23)
вместо (5.22) имеем
4
3
dI
d
=
H
,
(5.24)
Так как поток излучения H постоянен в фотосфере, то интегрирование (5.24) даёт
I
=
H
2
1
+
3
2
.
(5.25)
Здесь мы воспользовались граничным условием: 2I=H при =0.
Считая, что величина I равна полной интенсивности излучения при термодинамическом равновесии, т.е. I=T/, и выражая полный поток излучения через эффективную температуру Te по формуле
H
=
T
4
e
,
вместо (5.25) находим
T
=
T
4
e
1
2
+
3
4
,
(5.26)
т.е. ранее полученную формулу (4.20).
Таким образом, определяя средний коэффициент поглощения формулой (5.21) и пользуясь приближением Эддингтона, мы приходим к такой же зависимости между температурой и оптической глубиной, как и в случае, когда коэффициент поглощения не зависит от частоты. Однако вычислить точно величину мы не можем, так как в формулу (5.21) входит поток излучения H в реальной фотосфере, в которой коэффициент поглощения зависит от частоты. Поэтому средний коэффициент поглощения приходится вычислять приближённо.
Для приближённого вычисления величины были предложены следующие способы.
1. Будем считать, что поток излучения H равен потоку излучения из абсолютно чёрного тела, т.е. H=B(T) где B(T) — планковская интенсивность при температуре T. Тогда
=
B(T)d
B(T)d
.
(5.27)
2. Возьмём выражение для H, даваемое формулой (5.19). Заменяя в ней I на планковскую интенсивность B(T), находим
H
=-
4
3
1
dB(T)
dT
dT
dr
.
(5.28)
Подстановка (5.28) в (5.21) даёт
=
dB(T)
dT
d
·
1
dB(T)
dT
d
^1
.
(5.29)
Формула (5.29) была предложена Росселандом [2].
3. Примем для H выражение, которое получается в случае, когда коэффициент поглощения не зависит от частоты. Обозначая поток излучения для этого случая через
H
0
,
получаем
=
H
0
H
(5.30)
Формулу (5.30) предложил Чандрасекар [4], табулировавший также величину
H
0
H
.
Мы не будем сравнивать между собой различные способы вычисления величины Отметим только, что вычисления по формулам (5.27) и (5.30) проще, чем по формуле (5.29). Это особенно заметно в случае сложного химического состава, так как в формулы (5.27) и (5.30) члены, соответствующие разным атомам, входят аддитивно. Однако формула (5.29), по-видимому, точнее.
Для примера найдём средний коэффициент поглощения по формуле (5.27) в случае, когда поглощение вызывается атомами водорода.
Пользуясь формулой (5.11) для и формулой (4.2) для B(T), получаем
0
B
(T)
d
=
n
e
n
2^2ekT
33ch(2mkT)^3/^2
2h
c^2
x
x
0
1+2
kT
i=i
1
i^3
exp
kT
exp
–
h
kT
d
.
(5.31)
Здесь для простоты мы положили gi=1 и g=1. Меняя порядок интегрирования и суммирования и производя интегрирование, находим
0
B
(T)
d
=
n
e
n
2^2ekT
33ch(2mkT)^3/^2
2h
c^2
x
x
kT
h
1
+
2,4
kT
.
(5.32)
Кроме того, имеем
0
B
(T)
d
=
2h
c^2
kT
h
0
x^3dx
ex– 1
=
2h
c^2
kT
h
15
.
(5.33)
Подстановка (5.32) и (5.33) в формулу (5.27) даёт
=
40
3
eh
me(2m)^3/^2
1
+
2,4
kT
nen
(kT)/^2
.
(5.34)
Формулу (5.34) мы получили для атома водорода, но она справедлива без изменений и для водородоподобных ионов (так как атомный номер Z входит в ) Приближённо формула (5.34) справедлива и для других атомов.