книги из ГПНТБ / Каули, Ч. Теория звездных спектров
.pdf30 |
ГЛАВА f |
Чтобы связать эту величину с коэффициентами Эйнштейна А и В, нужно перейти от монохроматического плоскополяризованного электромагнитного излучения к функции распределения плотности энергии, которую мы возьмем равной функции рас пределения иѵ излучения абсолютно черного тела.
Плотность энергии излучения абсолютно черного тела в ин тервале частот (ѵ, dv) с центром на частоте излучения (1.5.17)
составляет (сравните с уравнением (1.3.4): заменено на 2А)
uv dv = ЗЛ2/2л. |
(1.5.22) |
Подставляя А2 из (1.5.22) в уравнение (1.5.21), получим |
|
атат = -у -«ѵі (т\ р х I«) Pf (ѵ, t)dv, |
(1.5.23) |
где через f(v, t) обозначено довольно громоздкое выражение в фигурных скобках (1.5.21).
Теперь мы должны проинтегрировать это выражение по всем V, поскольку, как и в классическом случае, атом способен ре
агировать на некоторую малую область частот около частоты ѵ,
определяемой соотношением |
(1.5.20). Полагая |
|
|
|
|
/ = |
Jоо f(v,t)dv, |
(1.5.24) |
|
получим |
—оо |
|
|
|
|
оо |
|
|
|
|
|
|
|
|
|
I ~ ( 2t/hh) J (sin2xlx2) dx. |
(1.5.25) |
||
Интеграл, |
входящий в (1.5.25), является табличным |
и |
равен |
|
л. Чтобы |
получить полное |
выражение для а*тат, нужно |
доба |
|
вить еще два слагаемых, учитывающих направления у и г. Окон чательный результат
а‘т(іт = |
| ( т | Р |а) р uvt, |
(1.5.26) |
где |
|
|
I (m| Р |п> Р = ! (т IРх Iп) f + |
1(т\Р„ |п) |2 + | (щ| Р2 |п) |2, |
(1.5.27) |
является матричным элементом полного электрического диполь ного момента атомов.
Мы находим, что а*тат пропорционально времени, |
следова |
тельно, |
|
Впт — (8л3/3/г2) |(m| Р |я) |2. |
(1.5.28) |
ТЕОРИЯ КОЭФФИЦИЕНТА ПОГЛОЩЕНИЯ В ЛИНИИ |
31 |
Можно связать теперь коэффициенты Эйнштейна А с матрич ными элементами дипольного момента при условии, что состоя ния m и л не вырождены (g m = g n — 1):
Атп= (64я4ѵ3/3/гс3) I(m| P |n) f. |
(1.5.29) |
Таким образом, скорость спонтанного излучения энергии кван товомеханической системой в расчете на один атом есть
(4со4/Зс3) \{т \ Р \п) I2. |
(1.5.30) |
Сравнение с соответствующим классическим выражением, кото рое можно получить из уравнений (1.1.4) и (1.1.5), показывает, что классическому дипольному моменту соответствует удвоен ный квантовомеханический момент, т. е. 2 | ( т | Р | « ) | .
Волновые функции Т т и Ч'п должны рассматриваться как функции, дающие наиболее полное из возможных описание атомных состояний т и п . Это значит, что оптический электрон характеризуется четырьмя квантовыми числами: п, I, / и или п, I, ті и т„, или, другими словами, (т\Р\п) выражает диполь ный момент для переходов между двумя элементарными зеемановскими состояниями, или просто состояниями. С другой сто роны, наши соотношения между коэффициентами А, В и f учи тывают переходы между вырожденными уровнями, имеющими статические веса g.
Следуя Кондону и Шортли [29], рассмотрим переход между двумя элементарными состояниями а и 4, где а — верхнее, а<5 — нижнее состояние. Пусть состояния а и 4 являются одним из gM- и состояний, образующих вырожденные уровни s i и Я.
Коэффициенты Эйнштейна А для перехода между состояниями а и 4 обозначим через А (а, 4). Подобным же образом обозна чим аналогичный коэффициент для перехода между двумя уров нями s i и Я: A (si,È). Энергия І(а,4), спонтанно излученная одной (зеемановской) компонентой, равна
I (а, 4 ) ~ dW/dt — NaA (а, |
4) hv — Na(4со4/Зс3)| (а | Р \ 4)f . (1.5.31) |
||
Величина A(si, |
Я) |
позволяет выразить общую энергию в спек |
|
тральной линии |
( |
s |
i которая получается суммированием |
интенсивностей всех (зеемановских) компонент (а->4). Соответ
ствующая интенсивность в линии I (si, Я) |
дается формулой |
I (si, Ш) = ЫЛА (si, ЯІ hv = Na(4co4/3c3) 2 ' | |
<«| P [*) |2. (1.5.32) |
ab |
|
Штрих при знаке суммирования означает, что сумма берется по всем возможным переходам между состояниями а и 4. Эту сумму Кондон и Шортли [29] назвали силой линии и обозначили
S(si, ^) = 2 l ( a | Р \4)\\ |
(1.5.33) |
ab
32 |
ГЛАВА 1 |
Здесь мы опустим штрих при знаке суммирования и заметим, что число слагаемых в сумме равно числу зеемановских компо нент в линии.
Поскольку Ы л = gjiNa, |
соотношения (1.5.32) |
и (1.5.33) дают |
при со = 2яѵ |
|
|
А (st-, Я) = |
- J - - ^ т - 5 (st, Я). |
(1.5.34) |
Аналогичные соотношения для коэффициентов Эйнштейна В, сил осцилляторов f и сил линий S можно получить из уравнений
(1.5.34), (1.4.9) и (1.4.10):
f {ß, |
= |
S (st, |
Я), |
(1.5.35) |
В [Я, |
s4) = -±- = |
^ - S ( s t , |
Я). |
(1.5.36) |
Преимущество понятия «силы линии» заключается в том, что сила линии одинакова для излучения и поглощения, т. е.
S(st, Я) = S(M, st), как легко видеть из определения |
(1.5.33). |
А соотношение (1.5.35) показывает, что величина |
st) |
также симметрична относительно уровней s t и Я. При переста новке s t и Я в уравнении (1.5.35) мы получим выражение для эмиссионной силы осциллятора *):
f(st, |
= |
(1.5.37) |
1.6. ЕСТЕСТВЕННАЯ ШИРИНА СПЕКТРАЛЬНОЙ ЛИНИИ
Мы видели, как классическое выражение для коэффициента поглощения можно ввести в квантовомеханическое выражение путем подстановки N —*■Nnfnm. Однако вопрос о том, как интер претировать в квантовой теории постоянную затухания, остается все еще открытым. Согласно классической теории, ширина спек тральной линии задается постоянной которая определяет уменьшение энергии классического осциллятора (уравнение
( 1. 1.8) )
W (t) = W0exp ( — yt). |
(1.6.1) |
Конечно, это излучение энергии непрерывно, тогда как в случае атома энергия освобождается дискретными порциями, или кван тами, равными hv. Чтобы найти аналогию с классическим слу чаем, мы должны рассмотреть ансамбль N осцилляторов. Число
*) Иногда в формулу (1.5.37) вводят знак минус. Делается это по тем же соображениям, по которым иногда вводят при рассмотрении линий погло щения знак минус для частоты линии излучения.
іевіЛЯЯС1' >.
ТЕОРИЯ КОЭФФИЦИЕНТА ПОГЛОЩЕНИЯ В ЛИНИИ |
33 |
классических осцилляторов Nmn, колеблющихся на частоте ѵтп, будет
N mn = W mn/ h \ mn, |
( 1. 6. 2) |
где Wmn — общая энергия всех осцилляторов (начальная). Для этих осцилляторов из (1.6.1) получим
dNmn/dt = — \N mn. |
(1.6.3) |
Согласно квантовой теории, скорость уменьшения заселенности уровня т из-за переходов с уровня т на п есть
(dN/dt)m^ n = - N mAmn. |
(1.6.4) |
Сравнивая (1.6.3) и (1.6.4), мы видим, что Атп играет роль постоянной затухания. Полная скорость уменьшения заселенно сти уровня т получается при суммировании по всем более низ ким, чем т, уровням п, так что мы считаем
Ут = 2 Атп. |
(1.6.5) |
п < т |
|
В квантовой теории постоянная затухания выражается сум мой обратных времен жизни как верхнего, так и нижнего со стояния, т. е. если
У п = 1 > А п1, |
(1.6.6) |
1< Ѣ |
|
Т О |
|
Утп = Ут + Уп. |
(1.6.7) |
Физический смысл того, что нижнее состояние расширяет спект ральную линию так же, как и верхнее, заключается в том, что верхнее и нижнее состояния размыты. Каждое состояние имеет конечную ширину Д1К, которая прямо связана со временем жиз ни состояния принципом неопределенности Гейзенберга, т. е. ДІЕ-Г^г/г, где Т — время жизни состояния. Если приписать уширение спектральных линий конечной ширине состояний, то мы должны пользоваться (1.6.7).
Если плотность излучения высока, то при вычислении вре мени жизни уровня нужно учесть также и переходы, индуциро ванные излучением. По мнению Слэтера [136], квантовая теория линии поглощения еще не полностью разработана. Для случая, когда поле излучения влияет на время жизни состояний, Миннарт и Мулдерс [111] предложили следующее выражение:
Ут = 2 Атп + |
^ иѵ (т -* п) Втп + |
uv {tn-*l)Bml (1.6.8) |
n<in |
п<т |
l> т2 |
2 Ч . К ау л и
34 |
Глава і |
и аналогичное |
выражение для уп■ Если в качестве иѵ можно |
взять функцию Планка, то, используя соотношения между коэф фициентами А и В и принцип детального равновесия, нетрудно выразить утп как функцию параметров атома и температуры.
1.7. УШИРЕНИЕ СПЕКТРАЛЬНЫХ ЛИНИЙ
Только в тщательно контролируемых лабораторных экспери ментах можно наблюдать линии с шириной, примерно равной естественной ширине. В большинстве же случаев спектральные линии будут уширяться в результате: а) теплового эффекта Допплера; б) турбулентных движений излучающего газа; в) эф фекта Зеемана и сверхтонкой структуры; г) многочисленных ти пов атомных, электронных или ионных столкновений (уширение давлением).
Здесь мы рассмотрим лишь механизмы уширения, приводя щие к профилям Гаусса или профилям Лоренца, и обсудим ре зультат одновременного действия нескольких таких механизмов. Общепринято, что турбулентное движение газа приводит к гаус сову профилю, тогда как многие типы уширения вследствие давления дают профиль Лоренца, т. е. профиль той же формы, что и при естественном уширении линии. Поэтому мы рассмот рим сумму двух гауссовых профилей, сумму двух профилей Лоренца и, наконец, сочетание профиля Лоренца с профилем Гаусса.
Вначале покажем, что тепловой эффект Допплера приводит к гауссову контуру. Объединим хорошо известную формулу Доп плера
АХ/Х = и/с |
(1.7.1) |
е известным фактом, что распределение скоростей ѵ по лучу зрения (одна компонента) является максвелловским. Если P(v)dv — вероятность того, что скорость частицы массы т будет заключена в интервале (и, dv), то
Р (V) dv = С, ехр ( — mv2/2kT) dv, |
(1.7.2) |
где постоянная Сі определяется условием |
|
оо |
* |
J P { v ) d v = 1. |
|
—00 |
|
Вероятность P(AX)dAX того, что спектральная линия будет сме щена в интервал (ДА, dAX)*), находится из (1.7.1) и (1.7.2):
Р (ДА) d АА = С2ехр [ — (ДА/ДА^)2] d ДА, |
(1.7.4) |
*) То есть интервал от ДА до ДА + d ДА.
ТЕОРИЯ КОЭФФИЦИЕНТА ПОГЛОЩЕНИЯ В ЛИНИИ |
35 |
где величина АXd , называемая допплеровской шириной, дается
формулой
AXD = (X/с) У 2kT/tn |
(1.7.5) |
и Сг — другая нормирующая постоянная.
Аналогично если мы предположим, что турбулентное движе ние газа сообщает излучающим атомам такое распределение по лучевым скоростям £, что вероятность P(l)dl попадания скоро
сти в интервал (|, dl) есть |
|
|
Р ( ^ )^ = |
С3е х р [ - Ш 2]. |
(1-7.6) |
то контуром линии будет |
|
|
Р (ДА) d ДА, = |
С4 ехр [ — (ДАс/А^)2]. |
(1.7.7) |
Постоянные С3 и С4 — нормирующие постоянные; It равно умно женной на У 2 среднеквадратичной турбулентной скорости вдоль луча зрения.
Отметим, наконец, что в результате столкновения или близ кого прохождения возмущающей частицы и поглощающего ато ма поглощенный фотон может быть смещен от естественной длины волны А0 на величину ДА. Рассмотрим контур Лоренца
(1.4.3) как функцию, |
пропорциональную вероятности |
P(AX)dAX |
|||
того, |
что поглощается |
фотон в интервале длин волн (ДА, dAX). |
|||
Тогда |
|
Csd АХ |
|
||
|
Р (ДА) d ДА |
(1.7.8) |
|||
|
ДА2 + |
(устА2/4 яс)2 |
|||
|
|
|
|
||
где |
С5 — нормирующая постоянная, |
а уст — параметр, |
который |
||
определяет ширину контура поглощения, обусловленную давле нием (гл. 5).
Четыре контура, с которыми мы будем иметь дело, даются формулами (1.4.3), (1.7.4), (1.7.7) и (1.7.8).
1.8. СВЕРТКА ДВУХ КОНТУРОВ
Теперь исследуем форму спектральной линии, которая уши ряется совместным действием двух механизмов. Результирую щий контур будет сверткой двух контуров. Ниже мы определим постоянные множители двух контуров, исходя из физических требований нашей задачи.
Если действует лишь один механизм уширения, то энергия в спектральной линии распределяется по контуру, который обо значим через fi(x), где х заменяет ДА. Энергия, содержащаяся в каждом элементарном прямоугольнике с основанием dx и вы сотой fi(x), перераспределяется вторым механизмом уширения по контуру fz(y), но с центром в точке х. Энергия, попадающая
2
36 |
ГЛАВА I |
|
|
б |
интервал [у, dy) в результате перераспределения |
энергии |
|
fi(x)dx, равна (рис. 1.8.1) |
|
|
|
|
f1{x)dxf2{y — x)dy. |
■ |
(1.8.1) |
Ее можно приравнять энергии элемента комбинированного кон тура d(fc(y)dy). Комбинированный контур получается суммиро ванием вкладов от всех таких элементарных полос, т. е.
оо |
|
fc (у) — J h ( x ) h ( y — x)dx. |
(1.8.2) |
—со |
|
Этот результат есть хорошо известное выражение для свертки двух функций распределений.
Рассмотрим три следующих случая:
1)/і и і2— контуры Гаусса,
2)fi и f2— контуры Лоренца,
3)fi — контур Лоренца, а f2— контур Гаусса.
Для двух первых случаев результаты можно получить при помощи теоремы о свертке, которую мы сейчас докажем.
Рис. 1.8.1. Свертка двух контуров.
Рассмотрим преобразование Фурье для функции f3(x), кото рая является сверткой функций fi(x) и f2(x):
со
/з (* )= J fЛу) f2 ІХ — у) dy. |
(1.8.3) |
—оо |
|
Предположим, что интеграл (1.8.3) существует и что существуют преобразования Фурье для fu f2, f3, которые обозначим через
Fi(k), F2(k), F3(k). Тогда
оо оо
F3( k ) = J exp ( — ikx) dx |
j fi (y)f2(x — y)dy. |
(1.8.4) |
—09 |
r - c o |
|
\
ТЕОРИЯ КОЭФФИЦИЕНТА ПОГЛОЩЕНИЯ В ЛИНИИ |
37 |
Проведем интегрирование сначала по переменной я. Для этого, переписывая
ехр ( — ikx) = exp [ — ik(x — z/)]exp( — iky), |
(1.8.5) |
||
получим |
|
|
|
оо |
оо |
|
|
Fa(k)= J fl (y)exp( — iky)dy |
J f2(x — y ) X |
|
|
—оо |
—-оо |
|
|
Хехр[ |
— ik{x —y)\d{x — у), |
(1.8.6) |
|
где. во втором интеграле вместо |
dx записано d(x — у), |
потому |
|
что во втором интеграле у постоянно. Следовательно, |
|
||
F3(k) = |
F1(k)F2(k), |
(1.8.7) |
|
или преобразование Фурье свертки равно произведению преоб разований Фурье функций, входящих в свертку. Иногда проще перемножить преобразования Фурье и найти обратное преобра зование Фурье от произведения, чем прямо вычислять свертку.
Теперь рассмотрим случай, когда f\ и f2— гауссовы контуры:
fi (х) = |
Аехр[ — {х/а)2] |
(а), |
U{х) = |
Вехр [ — (х/Ь)2} |
( 1.8 .8 ) |
(б ). |
Соответствующие преобразования Фурье равны
F, {k) — Аа У п ехр[— (а/г/2)2] |
(а), |
F2 (А) = ВЪ Ѵл ехр { — {Ыг/2)2) |
(1.8.9) |
(б). |
А преобразование Фурье F%(k) свертки /1 и f2 равно
F3(k) = АВаЬлехр[ — (а2 + b2)(k/2)2]. |
(1.8.10) |
Записывая с2 = а2+ Ь2, мы можем найти обратное преобразо вание Фурье Fa(k), т. е. /3, исходя из соотношения для fi{x) и Fi(k):
fa(x) — (ABab/c ] / л)ехр[ — (х/с)2]. |
(1.8.11) |
Мы можем не заботиться о постоянном множителе в этом выра-' жении, поскольку конечное выражение для коэффициента по-: глощения должно нормироваться так, чтобы выполнялось соот ношение (1.4.4). Наш результат можно резюмировать следую щим образом: Свертка двух функций Гаусса является гауссовой функцией, ширина которой равна квадратному корню из суммы квадратов ширин функций свертки.
Теперь положим |
|
і |
fi (#) = Aj{x2+ а2) |
(а), |
|
Ш = ВІ(х2 + Ь2) |
(б). |
и ’8Л ; |
38 |
ГЛАВА I |
|
|
|
Их преобразования Фурье соответственно равны |
|
|||
Fl (k) = |
A ^ e x p [ — a \k \] |
(а), |
(1.8.13) |
|
F2( k ) ^ B ^ e x p l - b \ k \ ] |
(б). |
|||
|
||||
Следовательно, |
|
|
|
|
F3l(k) = - ^ - n 2exp[— (a + b ) [k |], |
(1.8.14) |
|||
отсюда |
|
|
|
|
h (*) = АВя -т |
‘ |
(1.8.15) |
||
|
хг + (а + Ь)2 |
|
||
Итак, свертка двух контуров Лоренца является контуром Ло ренца, ширина которого равна сумме ширин контуров свертки.
Легко теперь скомбинировать гауссовы контуры теплового и
турбулентного уширений. Можно записать |
|
fD(ДА) = Сд ехр [ — (ДА/ДАд)2], |
(1.8.16) |
где теперь |
|
|
(1.8.17) |
и CD— нормирующая постоянная. |
совпадает с |
Комбинированный контур поглощения (1.8.15) |
выражением (1.4.3), с тем лишь исключением, что теперь мы понимаем под у сумму у от естественного уширения и разных типов уширений давлением. Вспомним, что интеграл контура поглощения не зависит от ширины. Функция хл(ДА) (1.4.3) нор мирована таким образом, чтобы дать результат (1.4.4), поэтому нам нужно выбрать такой нормирующий множитель CD, чтобы контур Допплера был нормирован к единице. Физически это означает, что полное или интегральное поглощение в элементар ном объеме не меняется с шириной профиля коэффициента по глощения.
Теперь мы в состоянии выразить коэффициент поглощения для затухания совместно с допплеровским уши'рением как сверт ку двух эффектов, действующих независимо. Найдем свертку коэффициента поглощения, выраженного как функция разности
длины волны ДА от центра линии, с контуром Допплера |
(1.8.16). |
На основании уравнения (1.4.3) можно записать |
|
fy (АЛ) = (e2l 2/2mc2) yKNJ |
|
где |
(1.8.19) |
Ѵх = у К І 2лс |
есть постоянная затухания в сантиметрах.
ТЕОРИЯ КОЭФФИЦИЕНТА ПОГЛОЩЕНИЯ в л и н и и |
39 |
|||
Нормированный контур Допплера есть |
|
|||
!° (АЛ) = |
у - - дя ехр [ - |
(ЛЛ/ДХ0П |
(1.8.20) |
|
Взяв свертку от (1.8.18) и |
(1.8.20), |
можно записать |
результи |
|
рующий контур |
|
|
|
|
f4D = ~ |
r - ^ |
N nfnmH{а, ДЛ/ДЛ0). |
(1.8.21) |
|
Параметр затухания а определяется соотношением |
|
|||
a |
2а = Уя/АЛ0, |
(1.8.22) |
||
|
оо |
|
|
|
|
|
|
|
|
|
|
— Jоо |
-<)■•• |
(1'8'23) |
Если существенно лишь допплеровское уширение, то форма коэффициента поглощения должна совпадать с (1.8.20), а нор мировка выбрана такой, чтобы удовлетворялось соотношение (1.4.4). Окончательное выражение будет
(АЛ) = |
KlNJnm |
ехр [ — (ДЛ/ДЛ0)*], |
(1.8.24) |
где |
|
|
|
Ѵпе21тс2 = |
4,99 • ІО-13 |
(ед. СГС). |
(1.8.25) |
Приложение IV знакомит с трансцендентной функцией Н(а,ѵ), определяемой соотношением (1.8.23). Оценка значений этой функции в нескольких точках контура при каждом шаге рас чета модели атмосферы может быть одним из самых трудоем ких этапов вычисления профиля линии. Вследствие этого изо бреталось множество алгоритмов для оценки Я (а, о).
Вычисление функции Я (а, ѵ), несомненно, является задачей, которая представляет возможность специалистам по приклад ной математике испытать свои силы. Но полезно отметить, что в большинстве физических ситуаций функцию Я (а, ѵ) нужно рассматривать лишь как некоторую аппроксимацию истинного профиля коэффициента поглощения. Во-первых, многие профили, уширенные давлением, не лоренцевы, а большинство других должно рассматриваться только как близкие к лоренцевым (гл. 6). Во-вторых, существенный вклад в допплеровские ши рины многих спектральных линий дают турбулентные скорости. Мы полагали, что турбулентные скорости имеют гауссов про филь, но это делалось только потому, что мы не знали истин ного профиля. По всей вероятности, негауссов профиль турбу лентных скоростей дает отклонение истинного контура коэффи циента поглощения от функции Я (а, и).
