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

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

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

80

ГЛАВА 1

где dtо — элемент телесного угла, а интегралы берутся по сфере. Используя оптическую глубину х\ —(к + а) т / х , получим урав­

нение переноса в виде

cos Ѳdljdx к = ІХ— xß/(x + оу) — оу4/(х + оу). (2.10.3)

Умножим уравнение (2.10.3) на daJAn и проинтегрировав по сфере, найдем

dtljdxk= X (JK-

В)/(х + од.

(2.10.4)

Умножение (2.10.3) наcos 0dio/4n

и интегрирование

приведут к

dKJdxk — Hk.

(2.10.5)

Приближение Эддингтона допускает, что поле излучения

изотропно всюду, кроме границы т = 0, где

существует только

выходящий

поток, изотропно

распределенный по

полусфере.

Если мы вынесем /% из-под знака

интеграла для /

и Н, а инте­

грирование

проведем только

по

полусфере

(0< ^ Ѳ ^ я /2), то

получим условие на границе

 

 

 

 

 

4 (0) =

2 # х (0).

 

(2.10.6)

Поскольку среднее по сфере значение cos20 равно 7з, то до­ пущение об изотропности поля излучения приводит к соотно­

шению

 

4 а= 4 / 3.

(2.10.7)

Используя соотношение (2.10.7), получим из (2.10.4) и (2.10.5)

d2(4 - B)/dxt = Зх (4 - ß)/(x + оу),

(2.10.8)

поскольку вторая производная от В равна 0. Интегрируя

(2.10.8), найдем

4 В = А ехр (— qxk),

(2.10.9)

где

(2 . 10. 10)

q = Y 3 x /( x + ok),

а А — постоянная интегрирования. Экспонента с положитель­ ным показателем исключена, так как не имеет физического смысла.

Теперь можно записать функцию источника

5 (ту) = В + ахА ехр (— qxK)/(x + оу).

(2.10.11)

Если функцию Планка выразить как функцию оптической глу­ бины та,, то

4 М :

о. -f- b —л-----ту ----- 4 — л ехр (— qxK)

X

 

X + а.

Л ' X + я.

 

 

Х ехр(— ту у) dxk!\x,

(2.10.12)

СПЕКТРАЛЬНЫЕ ЛИНИИ В ЗВЕЗДНЫХ АТМОСФЕРАХ

81

где ji = cos Ѳ. После интегрирования получаем

 

= а +

+

— (<?+ - ) •

(2.10.13)

Из соотношений (2.10.5), (2.10.6) и (2.10.1) находим постоянную интегрирования

'4 = ( І ,’ Т Т Т : - ' І) / ( 1 + Н

(2.10.14)

Легко видеть, что интенсивность в континууме равна

/ (р) = а + &р,

(2.10.15)

поэтому глубина линии записывается через удельную интенсив­ ность следующим образом:

 

q2bn

Р

 

bq2а 1+ ^ q

а +

+ X +

(X \ (?(Х +

1

Г%(р) = 1

 

а +

Ьц

(2.10.16)

 

 

 

 

Поучительны два предельных случая. Первый случай для ли­

нии

произвольной силы,

но <7 —>• 0, <тх/(х + ах)—>

1. Тогда для

всех

р = cos Ѳ дробь в

правой части (2.10.16)

обращается в

нуль. Следовательно, для достаточно сильных линий рассеяния глубина линии стремится к единице. Аналогичный результат можно получить для потока [158].

Второй полезный для рассмотрения случай, когда ок/к <С 1, т. е. случай очень слабых линий. Выражение для глубины линий -упрощается, и мы находим

 

2

,

 

а - 3-6

/*(р)

и ( а + Ьр)(|Сз |х +

(2.10.17)

 

1)(1 + 2 / ^ 3 )

Следовательно, глубины слабых линий прямо пропорцио­ нальны <зА/х. Сравним этот результат с результатом для случая чистого поглощения в модели Шустера — Шварцшильда. При малых оптических глубинах т можно разложить в ряд экспо­ ненту в уравнении (2.8.3) и получить

гі (р) = г0(р) Ѵ р .

(2.10.18)

В разд. 2.8 не было необходимости использовать индекс, а здесь это полезно сделать. Таким образом, тА = хлЯ [уравнение (2.8.5)], где — коэффициент поглощения в линии. Какую величину брать в качестве Я? Поскольку Я есть высота

82 ГЛАВА 2

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

т — оптическая глубина

в континууме, а

х — коэффициент не­

прерывного поглощения,

то Н « ( т— 1)/х.

Следовательно,

г*,(р) ~ (игЛ)Ыр)/р).

(2.10.19)

Ясно, что для обеих схематических моделей x j x и a j x играют сходную роль.

2.11ЗАМЕЧАНИЯ ОБ УПРОЩЕННЫХ МОДЕЛЯХ АТМОСФЕР

Вконце 40-х и начале 50-х годов было написано много ра­ бот, основанных на моделях атмосферы Шустера — Шварцшильда или Милна — Эддингтона. В эти исследования вклады­ валось столь много усилий, что иногда забывали, насколько упрощенными были эти модели. При этом многие исследова­ тели стремились сосредоточиться на модели Милна — Эддинг­ тона, забывая про модель Шустера — Шварцшильда.

Нередко механизм чистого поглощения критиковался потому,

что реальные линии Фраунгофера не исчезают на лимбе. Этот аргумент, в сущности, базируется на случае чистого поглощения для модели Милна — Эддингтона, когда, как легко показать, вы­ полняется соотношение

 

а +

b\ix/(x + kJ

(2 . 11. 1)

гх (е) = 1

а + 6р

которое приводит к

(0) = 0, т.

е. глубина линии

на лимбе

равна нулю.

Шустера — Шварцшильда предсказывает

Напротив, модель

обратно пропорциональную зависимость глубины слабой линии от р. Глубины слабых линий постепенно возрастают; глубины

сильных линий (в модели Шустера — Шверцшильда)

на лимбе

достигают значения г0 = 1— В (0) / (р) ф 0.

изменение

В действительности свойственное многим линиям

при переходе от центра к лимбу можно объяснить механизмом чистого поглощения [79] на основе модели атмосферы, учиты­ вающей зависимость параметров атмосферы от глубины. Такая модель сочетает многие свойства моделей Шустера — Шварц­ шильда и Милна — Эддингтона и обнаруживает свойства, кото­ рыми не обладает ни одна из этих моделей.

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

СПЕКТРАЛЬНЫЕ ЛИНИИ В ЗВЕЗДНЫХ АТМОСФЕРАХ

83

2.12КРИВЫЕ РОСТА ИЗ МОДЕЛИ АТМОСФЕРЫ.

ФУНКЦИИ ВКЛАДА

Когда в распоряжении исследователя имеются надежные данные о звезде, целесообразно выполнить вычисления, осно­ ванные на полной модели атмосферы.

Будем считать, что в нашем распоряжении есть подходящая

модель атмосферы, т.

е. имеются таблицы величин lg P g, lg Ре,

0 = 5040/7 и

для

различных т(Я = 5000 Ä) = то. Подробное

описание методов расчета

таких моделей приводится в обзор­

ной статье Михаласа

[110]

и учебниках [2, 158]. Много полез­

ных работ было выполнено с моделями, в которых распреде­ ление температуры было заимствовано у других звезд или моделей. Одной из основных статей, посвященных этому методу, является статья [23], а более поздние приложения даны Греем [59].

В настоящем рассмотрении обозначим функцию источника через Bh, так как в большинстве вычислений используется пред­ положение о ЛТР. Если берется функция источника в более общем виде, то применимы те же самые соотношения с заменой

Вк на Sx.

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

Глубину линии для Солнца можно записать так:

оо

d x \ /р

оо

 

J Вк ехр( - т 1 / ц )

- IВ %ехр(- х ф ) d x %^

 

ГьМ = 2-----------------то----------------2---------------------------- •

(2-12.1)

J В \

€Хр(-

% i l n ) d x l / l l

 

о

Знаменатель есть просто удельная интенсивность в континууме Я вблизи линии, он постоянен. Два интеграла в числителе можно свести к одному, обозначая

йтл = [ К + ^ )/Ч ] dxl

Легко также заменить переменную интегрирования т£ на То. Наконец, было показано [50], что удобнее интегрировать по lgto,

84 ГЛАВА 2

а не по самому то. Если мы выполним все эти операции, то по­ лучим

 

оо

ехР (— Я /р ) [1 — ехР (— ГІ/Е) К +

 

 

 

Гк(ц) =

—по

КІ) К \ X

 

 

А,

 

 

 

 

 

 

Х К / х 0)(т0/0,43429) r f lg v

 

(2 .12.2 )

Постоянная

0,43429 =

lg е иногда будет обозначаться

как Mod.

Можно назвать подынтегральную функцию в уравнении

(2.12.2) функцией вклада в глубину линии с£(т0).

Широко

ис­

пользуются

и

другие

определения функций вклада,

так

что

не нужно думать, что все функции вклада идентичны.

 

 

Чтобы вычислить

глубины линий посредством

выражения

(2.12.2) , нужно

1) найти коэффициент поглощения в линии в каждой точке атмосферы (ср. гл. 1);

2) определить оптические глубины в линии хх1 и конти­

нууме %кс:

 

 

(2.12.3)

о

о

3)вычислить интеграл для Я;

4)взять интеграл (2.12.2).

Эквивалентные ширины получаются интегрированием по всему контуру линии.

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

линии соотношением

 

00

оо

 

(2.12.4)

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

X. + А. 1

т0 d lgT0, (2.12.5)

 

х0 Mod

— оо

где F i— поток в континууме вблизи линии. Как и в случае

С П Е К Т Р А Л Ь Н Ы Е ЛИ Н И И В З В Е З Д Н Ы Х АТМОСФЕРАХ

85

удельной интенсивности, определим функцию вклада

C \(lgt0)

б глубину линии Rx как подынтегральную функцию выражения (2.12.5). Среднюю функцию вклада C(lgx0) можно определить

соотношением

00

СО

С = \ c xRxd ДА/ J RxdAX.

(2.12.6)

— со

— оо

 

Соотношения (2.12.2) и (2.12.5) универсальны. Они приме­ нимы как к сильным, так и к слабым линиям и простым спосо­ бом выражают зависимость от глубины таких параметров, как допплеровская ширина и постоянная затухания.

Рис. 2.12.1. Средние функции вклада для линий Nil в центре диска Солнца. Кривая для линии с потенциалом возбуждения 5,50 эВ смещена к более глубоким слоям по отношению к кри­ вой для линии с потенциалом возбуждения 0 эВ. [Astophys J., 143,352(1966).]

Рис. 2.12.1 иллюстрирует использование средних функций вклада для установления относительных глубин образования линий. Вычисления сделаны по удельной интенсивности двух солнечных линий Nil с одинаковыми эквивалентными шири­ нами, но разными потенциалами возбуждения: 0,00 и 5,50 эВ. Ясно видно, что функция вклада для высоковозбужденной ли­ нии смещена к большим оптическим глубинам.

86

ГЛАВА 2

Функции вкладов в форме, использованной в (2.12.2) или (2.12.5), при самых малых оптических глубинах отрицательны. Обычно

К + Ч Ж > !>

поэтому функция вклада отрицательна при малых тф. Когда тф

становится

порядка единицы, то

вклад

самопоглощения

ехр(—

или £’2(та+ та) быстро

спадает

к нулю, умень­

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

Интеграл функции вклада от 0 до т* определяет глубину линии, которая наблюдалась бы при отсутствии атмосферы от т* до оо. Отрицательная часть функции вклада показывает, что при отсутствии более глубоких слоев атмосферы мы наблюдали бы эмиссионные линии (отрицательные глубины линии)! Это не должно удивлять, так как в этом и заключается один из за­ конов спектроскопии — закон Кирхгофа, гласящий, что в спектре горячего газа наблюдаются эмиссионные линии. И только если за этим газом находится более горячий источник непрерывного спектра, то наблюдаются линии поглощения.

Г Л А В А 3

Статистическая механика

3.1.ВВЕДЕНИЕ

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

доля атомов, находящихся в объеме газа, способна поглощать в данной спектральной линии. Остальные атомы будут в других стадиях возбуждения или ионизации, а при низких температурах войдут в состав молекул.

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

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

Рассмотрим большое число частиц N, помещенных в полость и пришедших в состояние равновесия. Будем считать, что каж­ дая частица может находиться в любом из большого числа энергетических состояний Еп. Число частиц в энергетическом состоянии Еп обозначим через Nn. Наша задача — найти Nn как функцию Еп.

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

88 ГЛАВА 3

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

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

Рассмотрим описание таких свойств газа, как, например, среднее число атомов в 1 см3 в данном энергетическом состоя­ нии. Для многих задач подобного рода можно допустить диск­ ретность энергетических состояний без привлечения других квантовомеханических концепций и выводить соотношения, ко­ торые выполняются также и при более общих допущениях. Для статистического анализа важны два квантовомеханических принципа: принцип Паули и принцип неразличимости одина­ ковых элементарных частиц. Статистическое описание явлений без учета этих принципов дается классической статистикой, или статистикой Больцмана. Если принимается во внимание прин­ цип неразличимости частиц, то приходят к статистике Бозе Эйнштейна (или статистике Бозе), а если используются оба принципа, то имеют дело со статистикой Ферми Дирака (или статистикой Ферми).

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

Принцип Паули, хорошо известный из элементарных учеб­ ников атомной физики (см., например, [74, р. 122]), гласит, что два электрона в атоме не могут одновременно иметь четыре одинаковых квантовых числа. Эти четыре квантовых числа обеспечивают полное описание элементарного состояния элек­ трона, так что можно также сказать, что принцип Паули запре­

Статис тиче ска я механика

89

щает пребывание двух электронов в одинаковых элементарных состояниях*). Полное описание элементарного состояния сво­ бодной частицы, вероятно, менее знакомо читателю, так как квантовые представления для описания свободных частиц тре­ буются реже. Свободная частица находится в элементарном состоянии, когда ей приписана ячейка фазового пространства (положение х и импульс р) объема /г3 и задана ориентация ее спина. Некоторое понимание причины, почему ячейки размера /г3 играют важную роль в определении состояния свободной частицы, можно получить при обращении к принципу неопре­ деленности, который гласит, что можно «локализовать» частицу в фазовом пространстве только с неопределенностью АхАрх ~ h (см., например, [74, р. 36]). (См. также рассмотрение поведения частицы в ящике в разд. 3.4.)

Неразличимость частиц — свойство, которое проявляется только тогда, когда, например, две частицы близко располо­ жены в пространстве и имеют почти одинаковые импульсы. Если положения или импульсы двух идентичных частиц полностью разделены, то мы можем использовать эти свойства для иден­ тификации частиц: частица А — это та, которая находится в данном месте и с данным импульсом, а частица В — та, которая находится в другом месте и с другим импульсом. В идеальном газе нет скопления частиц в фазовом пространстве, т. е. частиц с одинаковыми положением и импульсом, и поэтому статистика Больцмана правильно предсказывает свойства идеального газа.

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

3.2. КВАНТОВАЯ СТАТИСТИКА И СТАТИСТИКА БОЛЬЦМАНА

Обратимся теперь к нашим N частицам, помещенным в по­ лость. Когда каждой частице приписана определенная энергия, мы получаем «состояние» для системы в , целом. Мы хотим узнать, какое распределение энергий частиц является наибо­ лее вероятным среди бесконечного множества возможных со­ стояний.

Попытаемся вычислить вероятность состояния системы, в ко­ торой имеется Nі частиц с энергией Еи N2 частиц с энергией Е2

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

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