Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.astro.spbu.ru/staff/viva/Book/ch4L/node9.html
Дата изменения: Fri Nov 19 19:22:03 2010
Дата индексирования: Tue Oct 2 02:24:29 2012
Кодировка: koi8-r

Поисковые слова: п п п п п п п п
Решения: Излучение



previous up next
Next: Спектры Солнца и звезд Up: Задачи Previous: Звездные величины

Решения

8. Излучение


gif 8.1 Функция Планка в шкале длин волн имеет вид
 displaymath5213
а приближение Вина дается формулой
displaymath5215
Поэтому
eqnarray1717
Значит, относительная погрешность приближения Вина tex2html_wrap_inline3263 = tex2html_wrap_inline5221. С другой стороны, исследуя функцию Вина на максимум, легко найти, что
displaymath5223
Этот очень полезный результат есть закон смещения Вина. Обычно ограничиваются тем, что отмечают постоянство произведения tex2html_wrap_inline5225. Однако очень важно и численное значение tex2html_wrap_inline5227, точнее, то, что это число заметно превосходит 1. Действительно, мы имеем tex2html_wrap_inline5229, так что tex2html_wrap_inline3263 мало при tex2html_wrap_inline5233. Интенсивность в максимуме, которую мы (теперь обоснованно) вычисляем в приближении Вина, есть
displaymath5235

Если не пользоваться приближением Вина, а работать с точной планковской функцией, то оказывается, что tex2html_wrap_inline5237 равно не точно 5, а 4.965 (проверьте!). Все остальное, включая заключение, что tex2html_wrap_inline5239, остается в силе.

Длинноволновое приближение Рэлея-Джинса, противоположное приближению Вина, обеспечивает относительную погрешность tex2html_wrap_inline5241 лишь при tex2html_wrap_inline5243 (проверьте). Скажем, точность в 10% приближение Рэлея-Джинса дает лишь при tex2html_wrap_inline5245, превосходящих tex2html_wrap_inline5247 в 25 раз!

В итоге оказывается, что вся планковская кривая tex2html_wrap_inline5249, как ее видит глаз на обычном графике, неотличима от виновской кривой. Большинство же студентов (да, пожалуй, и экс-студентов тоже) ошибочно полагает, что приближение Вина для tex2html_wrap_inline5249 применимо только слева от максимума, приближение Рэлея-Джинса начинает работать слегка правее него, а сам максимум хорошо описывается лишь точной формулой Планка. На самом деле, как мы убедились, все совсем не так.


gif 8.2 Планковские кривые, соответствующие разным T, не пересекаются, поскольку tex2html_wrap_inline5255. Отсюда следует, что при любом (фиксированном) tex2html_wrap_inline5245 значения tex2html_wrap_inline5249 монотонно возрастают с T. Далее, высота максимума планковской кривой, т.е. максимальное значение интенсивности, пропорциональна tex2html_wrap_inline5263 для tex2html_wrap_inline5249 и tex2html_wrap_inline5267 для tex2html_wrap_inline5269. Это легко показать, исследуя на максимум соответствующие функции Планка (см. задачу gif). Площадь же под обеими кривыми, и tex2html_wrap_inline5249 и tex2html_wrap_inline5269, растет tex2html_wrap_inline5275 (закон Стефана-Больцмана). Поэтому с ростом температуры кривая Планка в шкале длин волн "заостряется", а в шкале частот -- "притупляется".


gif 8.3 Пусть f(x) дифференцируема в точке tex2html_wrap_inline5279. Для простоты считаем, что tex2html_wrap_inline5281 и tex2html_wrap_inline5283. Для получения степенной аппроксимации f(x) в окрестности tex2html_wrap_inline5279, т.е. представления f(x) вида
displaymath5291
поступим следующим образом. Будем рассматривать tex2html_wrap_inline5293 как функцию tex2html_wrap_inline5295. Тогда имеем обычную линеаризацию в окрестности tex2html_wrap_inline5279:
displaymath5299
Отсюда, потенцируя, получаем вышеприведенную степенную аппроксимацию f(x), причем обнаруживается, что
displaymath5303
Подобные степенные аппроксимации используются в физике (и, в частности, в астрофизике) буквально на каждом шагу. К сожалению, ни в одном известном авторам курсе математического анализа об этом нет ни слова -- хотя учить этому следовало бы всех, даже изучающих анализ не слишком глубоко. Видимо, считается, что студент сам все это сообразит, когда немного "подрастет". Мы решили нарушить традицию и не ждать, когда это случится.

Приведенная в условии задачи степенная аппроксимация зависимости чернотельной интенсивности tex2html_wrap_inline5269 от T в окрестности tex2html_wrap_inline5309 получается только что описанным стандартным способом. Выкладку предоставляем читателю.

Степенная аппроксимация функции Планка, которую почему-то не отыщешь ни в одном учебнике, позволяет понять многие качественные особенности солнечного и звездных спектров. См., в частности, задачу gif.


gif 8.4 А почему, собственно, она должна равняться (3/2)kT? Ведь фотон -- не классическая частица, движущаяся с нерелятивистской скоростью. А только к таким частицам и применима классическая формула (3/2)kT.

Чтобы найти среднюю энергию одного чернотельного фотона tex2html_wrap_inline5315, надо объемную плотность энергии поля излучения
displaymath5317
поделить на число фотонов в единице объема
displaymath5319
Сделав в обоих интегралах одну и ту же замену tex2html_wrap_inline5321, обнаруживаем, что
displaymath5323
где
displaymath5325
Для оценки A можно воспользоваться приближением Вина, т.е. пренебречь 1 по сравнению с tex2html_wrap_inline5329 в двух последних интегралах (ср. с обсуждением в задаче gif). Тогда немедленно получим, что tex2html_wrap_inline5331, поскольку (советуем это запомнить)
displaymath5333
Последнее легко доказывается интегрированием по частям. Итак, tex2html_wrap_inline5335; вычислив интегралы точно, мы нашли бы, что A=2.70, так что окончательно
displaymath5339

Таким образом, средняя энергия одного чернотельного фотона без малого вдвое больше средней энергии теплового движения нерелятивистской частицы. Однако вклад каждого фотона в давление почти в точности такой же, как и каждой частицы: давление излучения tex2html_wrap_inline5341, газовое же давление P = n k T, где tex2html_wrap_inline5345 и n -- концентрации фотонов и частиц, соответственно. (Как вы думаете, почему так получается? Впрочем, это уже скорее физика, чем астрономия. Но ведь решаемся же мы нет-нет да и приучить вас к "физической математике", так почему же не поучить чуть-чуть и "астрономической физике"?) Из только что сказанного следует, что отношение концентраций фотонов и частиц есть одновременно (с точностью tex2html_wrap_inline5349) и отношение давления излучения к газовому (см. задачу gif).


gif 8.5 Выражение для tex2html_wrap_inline5345 уже появлялось в решении предыдущей задачи:
displaymath5353
Подстановка tex2html_wrap_inline5321 приводит его к виду
displaymath5357
где
displaymath5359
Таким образом, tex2html_wrap_inline5361. Чтобы найти точное значение коэффициента пропорциональности, надо получить C. Как мы знаем (см. решение предыдущей задачи), приближенно можно считать, что C=2. Точное же значение C получается так:
eqnarray1798
где tex2html_wrap_inline5369 -- дзета-функция Римана:
displaymath5371
Число tex2html_wrap_inline5373 не выражается через какие-либо "стандартные" постоянные (tex2html_wrap_inline5375, e, постоянную Эйлера и т.п.). Оно равно tex2html_wrap_inline5377.

После подстановки всех постоянных в полученное выше выражение для tex2html_wrap_inline5345 находим, что
displaymath5381

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

А вот как формулу tex2html_wrap_inline5383 можно получить совсем просто, комбинируя другие известные результаты. Плотность лучистой энергии равновесного излучения равна
displaymath5385
где a -- постоянная плотности излучения:
displaymath5389
а средняя энергия одного фотона tex2html_wrap_inline5391 (см. предыдущую задачу). Поэтому
displaymath5393
что сразу же и дает коэффициент 20 при tex2html_wrap_inline5395. На самом деле, конечно, ничего принципиально нового в таком способе расчета нет -- просто мы использовали готовое численное значение постоянной плотности излучения a (оно есть, например, у Аллена [1]).


gif 8.6 Частота фотона, испускаемого при переходе атома водорода с уровня m на уровень n, дается известной формулой
displaymath5403
где tex2html_wrap_inline5405 -- частота предела ионизации с первого уровня, tex2html_wrap_inline5407 Гц.

Перейдем от частот к длинам волн:
displaymath5409
где tex2html_wrap_inline5411 Angstrem.

Искомый переход найдем перебором. Вначале положим n=1. Тогда для m=2 получаем tex2html_wrap_inline5417Angstrem. Это знаменитая линия лайман-альфа, или tex2html_wrap_inline5419. Ясно, что для m>2 будем иметь tex2html_wrap_inline5423Angstrem, т.е. переходы на первый уровень нам не подходят. Возьмем n=2. Тогда для m=3 получим tex2html_wrap_inline5429Angstrem. Это и есть искомый переход (линия Htex2html_wrap_inline5431). Различие в длине волны в четвертом знаке (6 вместо 3) нас смущать не должно, так как при расчете мы использовали значение tex2html_wrap_inline4859 лишь с тремя значащими цифрами.


gif 8.7 По формуле из решения предыдущей задачи находим
displaymath5435
Таким образом, линия межзвездного водорода H tex2html_wrap_inline5437 лежит в субмиллиметровом диапазоне. Излучение в нем поглощается земной атмосферой. Наземные наблюдения линии невозможны.


gif 8.8 Линия H tex2html_wrap_inline5439 возникает при переходе tex2html_wrap_inline5441 в атоме водорода. Из общей сериальной формулы для водорода (см. задачу gif)
displaymath5443
полагая m=n+1 и считая, что tex2html_wrap_inline5447, находим
displaymath5449
Согласно этой формуле, линия H tex2html_wrap_inline5451, например, имеет длину волны около 5 см.

Подобные радиолинии, возникающие при переходах между близкими высокорасположенными уровнями, давно уже наблюдаются в туманностях. Как вы думаете, есть ли надежда обнаружить их также в радиоизлучении Солнца? (Ср. задачу gif.)


gif 8.9 Исходим из сериальной формулы для водорода в шкале частот:
displaymath5453
Полагая в ней tex2html_wrap_inline5455 и считая, что tex2html_wrap_inline5447 и tex2html_wrap_inline5459, получим
displaymath5461
Отсюда видно, что, действительно, при увеличении tex2html_wrap_inline5463 на единицу частота соответствующего перехода возрастает на одну и ту же величину tex2html_wrap_inline5465 которая есть не что иное как частота линии H tex2html_wrap_inline5439.


gif 8.10 Ионизация атомов водорода с n-го уровня может производиться фотонами с длиной волны короче, чем та, которую имеет излучение, образующееся при переходе атома водорода с уровня tex2html_wrap_inline5471 на уровень n. По формуле из решения задачи gif находим, что эта длина волны равна tex2html_wrap_inline5475 Angstrem. При n = 2 имеем 3648 Angstrem (на самом деле 3646 Angstrem). Таким образом, излучение видимого диапазона не способно ионизовать атомы водорода со второго уровня. Это очень важное заключение.


gif 8.11 Причина этого -- различие плотности. В атмосферах белых карликов она значительно выше, чем в солнечной хромосфере (почему?). Поэтому средние расстояния между атомами в хромосфере гораздо больше, чем в атмосферах белых карликов. Но радиус n-й боровской орбиты tex2html_wrap_inline5481 быстро растет с n, именно, tex2html_wrap_inline5485. Понятно, что он не может быть больше среднего расстояния между атомами -- иначе станет непонятно, какому именно атому принадлежит электрон, находящийся на этом уровне, т.е. произойдет его "обобществление". Поэтому, чем выше плотность, тем меньшее число уровней реализуется, а потому и тем меньшее число бальмеровских линий может возникать.

Удивительно, но факт: просто подсчитывая число бальмеровских линий, которые видны в спектре той или иной звезды, можно оценить плотность ее атмосферы!



ї В.В.Иванов, А.В.Кривов, П.А.Денисенков
HTML by Igor Drozdovsky
Последнее обновление:

previous up next
Next: Спектры Солнца и звезд Up: Задачи Previous: Звездные величины