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

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

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

100

ГЛАВА 3

отдельных атомов, а Еп— то же энергетическое состояние, отне­ сенное к основному (низшему энергетическому) состоянию моле­ кулы, то

El = En + D,

(3.4.23)

и сумма по состояниям для молекулы есть

 

QBH= exp (— D/kT) 2 gn exp (— EJkT).

(3.4.24)

П

 

Подставляя формулы (3.4.24), (3.4.14) и (3.4.21) в (3.4.12), ПО-

лучим

( N a N b \ i N C

( 2 я ц £ 7 У / а

Q a ( b h ) Q b ( b h )

exp

(— D/kT).

(3.4.25)

 

Vc

Ä*

Q^Jbh)

 

 

 

 

Здесь

Qc ( b h ) — сумма в

формуле (3.4.24).

Заменим

теперь

N a / V a

- + Na (число частиц типа А на единицу объема)

и т. д. и

введем

р = тА-тв/{тА + тв) как «приведенную» массу одной

из частиц.

3.5.СТАТИСТИЧЕСКИЕ ВЕСА И СУММЫ ПО СОСТОЯНИЯМ

Вастрономической литературе символом и обычно пользуют­ ся для обозначения суммы по состояниям для атомов и ионов. Символ б часто используется для обозначения статистических весов, но мы будем применять символ gn.

При вычислении сумм по состояниям, конечно, следует ис­ пользовать статистические веса, которые не пропорциональны объему. Пропорциональность объему имеет место для функции распределения кинетической энергии поступательного движения,

имы использовали эту пропорциональность в уравнении равно­ весия (3.4.25), записывая его через плотность числа частиц, а не через полное число частиц в воображаемой полости.

Согласно разд. 3.2, статистический вес gn есть число дискрет­ ных элементарных состояний с энергией Еп. Мы будем считать, что читатель встречался с понятием статистических весов в атом­ ной физике (ср. [74]). Некоторое число gn элементарных атомных состояний может иметь одну и ту же энергию Еп. Такой уровень называют ^„-кратно вырожденным, или имеющим статистиче­ ский вес gn. На языке квантовой механики gn — это число ли­ нейно независимых собственных функций ф„, которые соответ­ ствуют собственному значению энергии Еп (ср. [122, § 14]).

В атомах степень вырождения уровня с полным угловым моментом / равна 2/ + 1, это число и есть статистический вес та­ кого уровня. Мур [112] составила таблицы атомных энергетиче­ ских уровней, модернизированные затем в последующих публи­ кациях (см., например, [114]). В этих таблицах приводится пол­

СТАТИСТИЧЕСКАЯ МЕХАНИКА

101

ный угловой момент электронов /: орбитальный угловой момент плюс спин, а спиновый угловой момент ядра не учитывается. Следовательно, приводимая величина J не есть полный угловой момент, который требуется при вычислении статистических весов.

Если / — спиновый угловой момент ядра, то полный статисти­ ческий вес уровня равен (21 + 1)(2/ + 1) [74, р. 193]. Таким об­ разом, спин ядра вносит постоянный множитель в статистические веса всех энергетических уровней атома. Этот множитель сокра­ щается, когда берутся отношения чисел заселенностей при со­ ставлении уравнения Больцмана или ионизации. По этой при­ чине в астрономии обычно пренебрегают ядерным спином. В не­ которых атомах взаимодействие ядерного спина с электронной структурой достаточно велико, тогда желательно определить за­ селенность энергетического уровня с заданным значением пол­ ного электронного плюс ядерного углового момента F = J + I. Статистический вес такого уровня равен 2 Е + 1 , а уравнение Больцмана для него имеет вид

N (aJF) = N - Ш - exp [ _ %(aJF)lkT\,

(3.5.1)

где aJF означает уровень атома с полным угловым моментом F, с суммой орбитального и спинового угловых моментов электро­ нов /, а а — символ, описывающий электронную конфигурацию. Сумма по состояниям и должна задаваться выражением

и = (2І + 1 )2 (2 / + 1)ехр[ — x(aJ)lkT}.

(3.5.2)

аJ

 

Заметим, что энергетические уровни атома по-прежнему обо­ значаются посредством аJ в выражении для суммы по состоя­ ниям. Это соответствует обозначению таблиц [112], которые не учитывают сверхтонкой структуры. Соответствующие обозначе­ ния уровня и его расщепление приходится брать из рассеянных

влитературе ссылок или вычислять из теории. Краткое изложе­ ние вопросов по сверхтонкой структуре и некоторые библиогра­ фические ссылки приводятся в разд. 5.2.

Полное рассмотрение статистических весов и сумм по состоя­ ниям для двухатомных молекул заставило бы нас довольно по; дробно знакомиться со структурой молекул. Мы же приведем элементарное вводное изложение, которое послужит посылкой для многих вычислений. Дальнейшие подробности можно найти

в[132, 144]. При использовании уравнения Больцмана или урав­ нения диссоциации для гетероядерных двухатомных молекул (состоящих из разных атомов) ядерным спином можно прене­ бречь. Поэтому мы пока ограничимся рассмотрением таких моле­ кул, а позднее дадим необходимые уточнения для гомоядерных

двухатомных молекул (состоящих из одинаковых атомов).

102

ГЛАВА 3

Для наших целей будет использована простейшая аппрокси­ мация двухатомной молекулы колеблющимся ротатором. Слу­ чаи, когда это приближение не верно, рассматривались в [26, 40]. Вычисления в этих случаях в принципе совершенно просты и легко проводятся на электронно-вычислительных машинах. Но на практике могут быть трудности из-за недостатка физи­

ческих данных.

Если предположить, что электронная, вращательная и коле­ бательная энергии независимы, то полная сумма по состояниям для двухатомной молекулы есть

Q = QeQvibQroU

(3.5.3)

где Qe, QVib и Qrot — электронная, колебательная

и вращатель­

ная суммы по состояниям соответственно.

 

Рассмотрим теперь электронную сумму по состояниям. Со­ стояния электронов классифицируются в соответствии со значе­ нием числа А, которое является абсолютной величиной Л — (век­ торной) компоненты электронного углового момента L вдоль междуядерной оси. В отличие от случая зеемановского расщеп­ ления, при наличии двух ядер электроны подвергаются дей­ ствию электрического поля с энергией взаимодействия, которая одинакова для состояний с одинаковой проекцией L на междуядерную ось. Таким образом, состояния с одинаковым Л яв­ ляются двукратно вырожденными, так как + Л и —Л приводят к одному и тому же значению Л = |Л |. Если Л = 0, то двукрат­ ного вырождения состояния не возникает. Вырождение, возни­

кающее из-за углового момента электронов, можно

записать

в виде

 

(2-бол ).

(3.5.4)

Если интерпретировать набор возможных состояний как мо­ лекулярные орбитали, то можно видеть, что при заданном А формула (3.5.4) объясняет все остальное орбитальное вырожде­ ние. Кроме орбитального вырождения, имеется 2S -)- 1-кратное спиновое вырождение, как и у атома. Таким образом, статисти­ ческий вес ge для данного электронного состояния равен

ge — (2 — б0л) (2S -ф- 1).

(3.5.5)

Для ряда двухатомных молекул можно определить энергию электронных термов (энергия в см-1) Тп и получить электрон­ ную сумму по состояниям

Qe = 2 (2 — бол) (25 + I) ехр (— hcTJkT).

(3.5.6)

П

 

Молекулярные спектры изучены значительно хуже, чем атом­ ные, и не всегда можно быть уверенным, что все нужные зна­ чения термов для вычисления (3.5.6) известны.

СТАТИСТИЧЕСКАЯ МЕХАНИКА

ЮЗ

Колебательные состояния не вырождены, следовательно, их статистические веса равны единице. Принято включать энергию нулевого колебательного состояния в величину электронного терма. Если ѵ — колебательное квантовое число, то значения терма есть G0(v)= vhvosJc, где vosc— частота колебаний. Отме­ тим, что это соотношение, не учитывающее энергии нулевого колебательного состояния, является просто планковским посту­ латом относительно энергии рассматриваемых им осцилляторов. Величину Vosc можно определить из спектроскопических дан­ ных. Колебательная сумма по состояниям становится равной

Qvtb^= 1 + ех р [— hv0JkT]-\-exv[— 2hvosc/kT] + . .. .

(3.5.7)

При X = exp [h v o sJk T ]

получаем

 

Qvib 1

X-f X2 -f . . . = 1 /( 1 — х),

(3.5.8)

откуда следует приближение

 

Qvib ~

[1 — ехр (— hvosc/kГ)]-1.

(3.5.9)

Наибольший вклад в полную сумму по состояниям двух­ атомной молекулы дает вращение. В магнитном поле вектор углового момента J может иметь 2/ -К 1 ориентацию, следова­ тельно, статистические веса равны 2 / - f l . Для обычного рота­ тора энергия уровня с квантовым числом углового момента I равна

Е = / (7 + 1) h2/2I,

(3.5.10)

где / — момент инерции молекулы *), который можно оценить из спектроскопических данных. Для реальных молекул вводится константа вращения Вѵ, определенная таким образом, чтобы

У (/+ 1 )Я Ѵ

(3.5.11)

давало наилучшее приближение к истинным значениям вра­ щательных термов колеблющегося ротатора, выраженных в см-1. Тогда

 

СО

 

 

 

 

Q,0/= 2 ( 2 /+

1)

ехр[— / ( / + l)Bvhc/kT]>

(3.5.12)

 

/=о

 

 

 

Сумму (3.5.12) можно аппроксимировать интегралом

 

 

оо

 

 

 

Qro,

{21+ 1) ехр [— / ( /

+

l)Bvhc/kT]dJ^kT/hcBv.

(3.5.13)

 

о

 

 

 

Вернемся теперь к вопросу о двухатомных молекулах, со­ стоящих из одинаковых атомов. У таких систем имеется центр

*) Символ I также используется для обозначения спинового углового момента ядер.

104 ГЛАВА 3

симметрии, налагающий некоторые ограничения чисто квантово­

механического происхождения на число

возможных состояний

с заданной энергией [75, § III-2, Ѵ-2;

102,

§ 75, 83;- 122, § 43, 48].

Полный статистический вес вращательного уровня гетеро­

ядерной молекулы в (21а + 1) (2/в +

1)

раз больше, чем дает

вычисление по формуле (3.5.13), где

іа

и ів — ядерные спины.

Полный статистический вес гомоядерной молекулы с учетом ядерного спина можно получить умножением вращательного статистического веса 2/ + 1 на множитель, который, следуя [144], мы назовем «ядерным спиновым статистическим весом».

Этот множитель зависит, во-первых, от класса симметрии

(s

или а) ядерной собственной функции и, во-вторых, от того,

яв­

ляется ли спиновый угловой момент ядра целым или полуцелым. Ядерные спиновые статистические веса представлены в табл. 3.5.1.

Таблица 3.5.1

Ядерные спиновые статистические веса гомоядерных молекул

 

Классификация уровня

S

 

а

/

полуцелое

(21 + 1) /

( 2 / +

1)(/ + 1)

I

целое

( 2 / + ! ) ( / + 1)

(2/ +

1)7

Для гетероядерных молекул ядерный спиновый статистиче­ ский вес всегда равен (21А + 1)(2Ів -\- 1). Сравнивая это зна­ чение с данными табл. 3.5.1 и учитывая примерно одинаковую встречаемость уровней класса симметрии s и а, получим, что ядерные статистические веса гомоядерных молекул вдвое меньше весов гетероядерных молекул. Таким образом, можно вычислить сумму по состояниям для гомоядерных молекул из соотношения

п

ЬТ

(2/+1)2

(3.5.14)

Чг°‘

hcBv

2

 

Если (3.5.14) используется в уравнении диссоциации, то спино­ вый угловой момент ядра нужно также включить в атомную сумму по состояниям. Поэтому часто при вычислении диссоциа­ ции гомоядерных молекул опускают множитель ( 2 /+ 1 ) 2 как

в(3.5.14), так и в атомной сумме по состояниям.

3.6.УРАВНЕНИЕ ИОНИЗАЦИИ САХА

ИУРАВНЕНИЕ ДИССОЦИАЦИИ

Равновесие между атомом и ионизованным атомом и элек­ троном описывается уравнением (3.4.25). Необходимо только правильно определить сумму по состояниям и поправку нуль-

 

СТАТИСТИЧЕСКАЯ МЕХАНИКА

105

пункта энергии D. Из этого уравнения мы имеем

 

N AN „

(2mnkT)3^ 2и.

(3.6.1)

N-

■exp(— Xi/kT).

hA

 

Здесь уже учтено, что электронная сумма по состояниям равна 2 в соответствии с двумя возможными ориентациями спина и что вместо приведенной массы ц можно использовать массу электрона т. Энергетические уровни иона отсчитываются от его основного состояния, которое выше основного состояния атома на величину %і (эВ). Следовательно, D = %і.

Перейдя от плотностей чисел частиц к парциальным давле­ ниям Рі, Ро и Ре для иона, нейтрального атома и электрона со­ ответственно, получим уравнение Саха, или, как его принято

называть в астрономии, уравнение ионизации

 

 

lg ^7 = - lg ^

+ 4 lg г + lg

- Ѳх, -

° ’48-

<3-6-2)

Здесь %і нужно брать в электрон-вольтах,

Ре— в динах на 1 см2,

а Ѳ= 5040/7’.

и колебательные

суммы

по состояниям

Если вращательные

аппроксимируются формулами (3.5.9) и (3.5.13) соответственно, то уравнение диссоциации молекул приобретает вид

РАРв

» АиВ

(2np,kT)sP

hcBv X

Рс

Q-symQe(^)

h3

 

X [1 — exp(— Av0JC/fer)]exp(— D/kT). (3.6.3)

Это соотношение приближенно выполняется для гомоядерных молекул, если Qsym = ѴгДля гетероядерных молекул Qsym — 1.

Построение теории, необходимой для вывода уравнений ионизации и диссоциации, было проведено химиком ван’т Хоф­ фом. Впервые эти идеи были применены к рассмотрению иони­ зации атома не Саха, а Эггертом [47], и в физике плазмы урав­ нение ионизации часто называется уравнением Эггерта — Саха. Интересы Эггерта были сосредоточены на веществе в недрах звезд. Саха [128] применил ионизационное уравнение к звезд­ ным атмосферам и тотчас же установил, что оно дает ключ к пониманию спектральной последовательности звезд. Вероятно, исключительная важность этой работы для астрономии привела к тому, что обычно в астрономической литературе уравнение ионизации и называется уравнением Саха.

3.7. ТЕОРИЯ ДЕБАЯ И ПЛАЗМЕННЫЕ ВЗАИМОДЕЙСТВИЯ

Дебай и Хюккель нашли интересное приложение формулы Больцмана к теории ионных взаимодействий. Их работа свя­ зана с подвижностью ионов в электролитических (жидких)

106 ГЛАВА 3

растворах. Было установлено, что скорость дрейфа ионов в рас­ творе по направлению к электродам несколько меньше, чем это следовало из вычислений. Причина такого расхождения объяс­ нена в теории Дебая и Хюккеля.

Согласно теории Дебая — Хюккеля, каждый ион окружен «облаком» противоположно заряженных ионов. Поэтому дрейф иона через раствор происходит несколько медленнее, так как ему приходится тащить за собой свою «свиту»*). Понятия тео­ рии Дебая — Хюккеля широко используются в физике плазмы, где ее называют просто теорией Дебая.

Пусть ф(г) — работа, необходимая для переноса единичного положительного заряда из бесконечности в точку с радиусом-

вектором

г в плазме.

По

определению ф(г)— потенциальная

энергия

единичного заряда

в точке г.

Потенциальная

энергия

положительного иона

равна

Z+e\р (г),

а отрицательного

иона —

Z~e\р (г). В этих обозначениях Z+ Z~ и е — положительные числа. Запишем формулу Больцмана в виде

N+ (г) — А ехр [— Z+ety {г)/кТ],

(3.7.1)

где А — постоянная, которую еще нужно определить, и огра­ ничимся случаем, когда показатель экспоненты (3.7.1) мал. Тогда

N+ (r) = A { \ ~ \ Z +e^(r)lkT\).

(3.7.2)

Такое приближение должно выполняться в большинстве ситуа­ ций для плазмы. Правда, если средние энергии кулоновского взаимодействия значительно больше тепловых энергий, то мож­ но ожидать и кристаллизацию!

Проинтегрируем (3.7.2) по объему нейтральной плазмы, за­

писав элемент объема dxdydz как dr:

 

I N+ dr = AV.

(3.7.3)

V

 

Интеграл по объему члена в квадратных скобках в (3.7.2) об­ ратился в нуль, поскольку среднее значение ф(г) для нейтраль­ ной плазмы равно нулю. Из (3.7.3) следует, что А = No, где No — средняя плотность числа положительных зарядов.

Те же аргументы можно применить и к N~(r)— плотности

числа отрицательных зарядов, следовательно,

 

N+ (г) = jVo*ехр [— Z+eip )/kT\,

 

N~(r) — Nö ехр[+ Z~ety {г)/кТ}.

(3.7.4)

*) Гласстоун [58] применил теорию Дебая — Хюккеля к проводимости электролитов.

СТАТИСТИЧЕСКАЯ МЕХАНИКА

107

Рассмотрим фиксированный отрицательный заряд. Вблизи него

ф(г) отрицательно, таким образом, из уравнений (3.7.4)

сле­

дует N+(r) > /Ѵ~(г), т. е. вокруг отрицательного заряда

суще­

ствует облако положительных зарядов. У читателя может воз­ никнуть вопрос, как объяснить с точки зрения геометрии, что в нейтральной плазме положительные заряды экранируют от­ рицательный и наоборот. Для иллюстрации на рис. 3.7.1 пока­

зана

решетка обычной поваренной

 

 

 

соли.

Конечно, кристаллы

значи­

 

 

 

тельно отличаются от плазмы, но

 

 

 

геометрия конфигурации легче вос­

 

 

 

принимается, если мы рассматри­

 

 

 

ваем

неподвижные частицы.

Ионы

 

 

 

Na+ и С1~ представлены черными и

 

 

 

белыми кружками.

На шести

гра­

 

 

 

нях, которые окружают централь­

 

 

 

ный

белый кружок

(С1~),

имеется

 

 

 

14 черных кружков

(Na+) и 12 бе­

 

 

 

лых. Таким образом, облако вокруг

 

 

 

иона СИ положительное, а кристалл

 

 

 

в целом электрически нейтрален.

Рис. 3.7.1. Часть

кристалли­

Теория Дебая утверждает, что ана­

ческой

решетки

хлористого

логичная ситуация

осуществляется

 

натрия.

 

в среднем в плазме.

потенциала ф(г)

можно

вывести

на основе

Выражение для

уравнения Пуассона

Ѵ2ф =

4лре,

 

(3.7.5)

 

 

 

где ре — плотность электрического заряда. В общем случае мо­ жет быть несколько сортов положительных и (для общности мы также будем считать) отрицательных ионов, так что результи­ рующая плотность числа зарядов равна

Ре = е ( S Zf N t - 2 ZJN jy

(3.7.6)

Используя два первых члена разложения (3.7.4), получим

ре = - е2 [2 ( z t f N t i + 2 (Zf)2Wol] т т .

(3.7.7)

Введем обозначение

(3.7.8)

и будем предполагать сферическую симметрию. Тогда уравне­ ние Пуассона примет вид

Ц24> (г)

2 Цф (г)

(3.7.9)

d r 2

+ г dr = Р£>2Ф(7),

108

ГЛАВА 3

 

где считается,

что ф является только функцией скаляра

г. При

г —►оо имеем соотношение ф" =

const X Ф. свидетельствующее

об экспоненциальности решения,

а при г —> 0 мы знаем,

что по­

тенциал должен меняться по закону 1/г (зависимость для то­ чечного заряда). Легко видеть, что соотношение

ф (г) = ±

Z±e exp ( г/рд)/г

(3.7.10)

является решением (3.7.9).

Знаки выбираются так,

чтобы

при

г —►0 получить истинный потенциал: для положительного

иона

берется + Z +, а для отрицательного —Z~.

В большинстве интересных для астрофизики случаев можно считать Z+ = Z~ = 1, откуда для радиуса Дебая получаем

Ро = (£778ле2АдЧ

(3.7.11)

где Ne — плотность числа электронов.

Дебаевское экранирование довольно существенно при интер­ претации ряда плазменных. явлений. Из формулы (3.7.10) мы видим, что потенциальная энергия' быстро спадает к нулю на расстояниях, значительно превосходящих pD- В плазме действие иона не проявляется на расстояниях значительно больших, чем ро- Мы воспользуемся этим обстоятельством в гл. 6 при рас­ смотрении уширения линии вследствие кулоновских взаимодей­ ствий. Сила кулоновского взаимодействия спадает по закону 1/г2, а число частиц возрастает как г3. Уширяющее действие мно­ гих зарядов на излучатель в данной точке обратилось бы в бес­ конечность, если бы не дебаевское экранирование.

3.8. АТОМНЫЕ СУММЫ ПО СОСТОЯНИЯМ. ПРОБЛЕМА СХОДИМОСТИ

В большинстве задач достаточно вычислять атомные суммы по состояниям суммированием по известным энергетическим уровням атома [112]. На практике необходимость включать все известные уровни возникает редко. Для примера рассмотрим Fel, для которого в циркуляре Национального бюро стандар­ тов № 467 приведен список, содержащий около 500 уровней. При температуре 5040 К вклад от последних 400 уровней в сумму по состояниям составляет меньше 1%. При более высоких тем­ пературах относительный вклад сильно возбужденных состояний несколько возрастает, но тогда большая часть железа в звезд­ ных атмосферах однократно ионизована.

Рассмотрим для простоты концентрацию нейтральных ато­ мов железа в s-м состоянии возбуждения. При вычислении силы линии, возникшей при переходе с этого уровня, нам потребуется знать отношение N0, s/ (Д70 -+- N\), где N0 и Л/)— полные числа нейтральных и однократно ионизованных атомов. Если железо

СТАТИСТИЧЕСКАЯ МЕХАНИКА

109

преимущественно однократно ионизовано, то N0 + Ni « Nu сочетание уравнений Больцмана и Саха дает соотношение

N о, і

g 0, s N e

h2

'h

Xo - %o, .

(3.8.1)

/V,

1

2nmkT

exp

kT

 

 

которое не зависит от суммы по состояниям для Fel.

До сих пор принимались во внимание только известные уров­ ни. Ясно, что важно рассмотреть уровни, которые должны суще­ ствовать, на пока еще не известны. Для изолированного атома число теоретических энергетических уровней бесконечно. Для водорода и ионизованного гелия теория Бора предсказывает

уровни с главными квантовыми числами п и энергиями

(отсчи­

тываемыми от основного состояния)

 

=

Ry(l — Д )

(3.8.2)

для всех значений я. Здесь Ry — постоянная Ридберга, a Z — за­ ряд ядра, воздействующий на электрон. Можно ожидать, что для других атомов и ионов справедлива та же самая формула с заменой я на «эффективное» квантовое число *)

«’ = £ [R y /ta ,-x J ]‘/2

(3-8.3)

и с учетом мультиплетности (2S1 + 1) (2L1 + 1) исходного терма.

Теперь сумма по состояниям должна включать высоколежащие уровни и, следовательно, иметь слагаемые вида

(2S' + 1) (2U + 1 ) ^ 2 (л*)2ехр{ - [и, - д ^ г ] Д Д , (3.8.4)

где знак штрих при суммировании показывает, что в этой фор­ муле учитываются только высоколежащие водородоподобные со­ стояния.

При любой заданной температуре ряд в формуле (3.8.4) при я * оо не будет сходиться, так как множитель Больцмана при­

ближается к постоянной величине ехр(—%0/kT). Ясно, что здесь что-то неверно, так как если, например, суммы по состояниям расходятся, то отношения

Nr,s/Nr = (gr. Jur) exp (— Xr, JkT)

(3.8.5)

должны стремиться к нулю, и мы не сможем наблюдать никаких спектральных линий! Нетрудно видеть, в чем источник затруд­ нений: атомный радиус возрастает пропорционально я2, таким образом, если сумма по состояниям неограниченно растет, то и размеры атома также должны расти неограниченно.

*) Если исходный терм не является основным состоянием следующей стадии ионизации, то не будет обычным потенциалом ионизации [164].

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