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

книги из ГПНТБ / Каули, Ч. Теория звездных спектров

.pdf
Скачиваний:
21
Добавлен:
21.10.2023
Размер:
9.13 Mб
Скачать

30

ГЛАВА 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 у) 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). Во-вторых, существенный вклад в допплеровские ши­ рины многих спектральных линий дают турбулентные скорости. Мы полагали, что турбулентные скорости имеют гауссов про­ филь, но это делалось только потому, что мы не знали истин­ ного профиля. По всей вероятности, негауссов профиль турбу­ лентных скоростей дает отклонение истинного контура коэффи­ циента поглощения от функции Я (а, и).

Соседние файлы в папке книги из ГПНТБ