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

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

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

130 ГЛАВА 4

полезно сравнить ее измерения с опубликованными эквивалент­ ными ширинами для сходной звезды. Две такие кривые сравне­ ния были описаны в конце разд. 2.1.

Для многих звезд сравнение средних постоянных затухания позволяет оценить среднее газовое давление Pg. Из теории уши-

рения ван дер Ваальса

(п = 6) можно вывести

(ср. гл. 5 и 6)

lgY =

lgPff — 0,3 lg Т + const.

(4.5.4)

Если известно, что уширение ван дер Ваальса преобладает (бо­ лее холодные звезды), то это уравнение можно использовать в форме

М = [**] — 0,3 [Г].

(4.5.5)

В (4.5.5) некоторые систематические ошибки

в эквивалентных

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

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

И первоначальная кривая роста для звезды сравнения, и кривая роста для исследуемой звезды строились в предположе­ нии, что lgx£ не меняется заметным образом для всех иссле­

дуемых линий.

Строя отдельные кривые роста для разных областей длин волн, можно оценить зависимости lgx£ от X.

Когда lgx£ сильно меняется с X, как, например, вблизи бальмеровского скачка в звездах класса А, следует ожидать раз­ личия средних атмосферных параметров в двух интервалах длин волн. В принципе можно выполнить дифференциальный анализ отдельно в различных спектральных областях, но не следует ожидать высокой точности от такого анализа, поскольку в каж­ дом отдельном интервале мало данных и средние физические условия хорошо определяться не будут. Дифференциальный анализ содержания элементов для звезд класса А проводить на­ много труднее, чем для других спектральных классов.

4.6. СРЕДНЕЕ ЭЛЕКТРОННОЕ ДАВЛЕНИЕ

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

К О Л И Ч ЕС ТВЕН Н Ы Й ХИ М И ЧЕСКИ Й АНАЛИЗ ЗВ Е ЗД Н Ы Х АТМ ОСФЕР 13)

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

Большинство методов определения [Ре] связано с уравнением

Саха, которое записывается так:

 

 

lg (iV+/A/0) = - l g P e + 2 ,5 1 g r - e x +

lg(«+/»0) - 0 ,1 8 . (4.6.1)

В дифференциальной форме оно проще:

 

 

[N+I№\ - [u+lu°\ = - [Ре] +

2,5 [Г] - ДѲх.

(4.6.2)

Измерения звездных спектров дают возможность оценить ле­ вую часть уравнения (4.6.2). По отдельным средним кривым роста для звезды сравнения и исследуемой звезды можно по­ лучить lg X = \g(W'/h). Затем для линий ионов используется формула

[JV+] + [С/и+кІ1 = [Х+] + ДѲх».

(4.6.3)

полученная из (4.4.5) при переходе к разности логарифмов, а для нейтральных линий аналогичная формула

 

 

[ Л + [С /Л £ Ы * ° ] + ДѲх„.

(4.6.4)

В

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

и (4.6.4)

стоят известные

величины

если

мы

можем считать,

что ДѲвозб

известно из кривой роста.

Вообще

довольно хорошая оценка ДѲвозб получается лишь для

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

стут с увеличением ДѲ.

___________

 

Определив

 

а+ =

[Д"+] +

ДѲхп,

(4.6.5)

а0=

[> ] +

ДѲхя,

(4.6.6)

где усреднение проводится по всем ионизованным и нейтраль­ ным линиям данного элемента соответственно, и вычтя среднее значение (4.6.4) из среднего (4.6.3), получим*)

[N+In °] + [uju+\ =

а+ а0= -

[Рв] — 2,5 [Ѳ] — ДѲх. (4.6.7)

При выводе (4.6.7) считалось, что

одинаково для

нейтраль­

ных и ионизованных

атомов. Поскольку величины а

известны,

*) Здесь мы заменяем средние логарифмов логарифмами средних. Это допущение должно быть хорошим, поскольку усредняемые числа должны быть близки.

5*

132

ГЛАВА 4

величина [Ре] получается сразу, как только определены А0 и [Ѳ]. Все этапы определения (4.6.7) можно выполнить графически. Описанный метод идеально подходит для машинных вычисле­ ний. Когда вычислительная машина выполнит численные ре­ дукции, измеренные значения 1g(W%/%) обычно выдаются на пер­ фокартах. Легко ввести в машину затабулированные средние кривые роста и непосредственно вычислить а+ и а0. Если это сделано, то разумно вывести на печать некоторые промежуточ­ ные результаты или графики, так как машина включает в усред­ нения также величины, содержащие грубые и очевидные ошибки.

[n*/№]

Рис. 4.6.1 Различие в степени ионизации для звезды HR 774 и звезды сравнения я6 Огі [35].

Чаще всего для определения [Ре] оценивают ДѲ и затем про­ водят вычисления по формуле (4.6.7) для Fel, Fell; Til, Till и т. д. Величины АѲ можно получить из фотометрических изме­ рений, или из АѲвозб, или даже из точной спектральной класси­ фикации, если это возможно. Фотометрические измерения, ве­ роятно, дают более точную оценку АѲ, чем можно получить из предположения, что АѲ = АѲП0.зб. Но АѲвозб не подвержено эф­

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

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

Если несколько элементов хорошо наблюдается в двух

ста­

диях

ионизации, то возможно одновременное решение для

[Ре]

и для

АѲ. Строится график величины а+а0 в функции /

 

для

разных элементов. Наклон этой кривой дает АѲ, и отсюда можно определить Ѳ и \Ре\

К О Л И Ч ЕС ТВЕН Н Ы Й ХИ М И ЧЕСКИ Й АНАЛИЗ ЗВ Е ЗД Н Ы Х

АТМ ОСФЕР

133

Рис. 4.6.1 иллюстрирует кривую такого рода [35]. Ординатой

рисунка фактически является а+— а0, поскольку

температуры

обеих звезд были настолько близки, что [и°/м+] =

0. Точки

для

Sc и Zr не были определены надежно и их разбросу не припи­ сывалось физического смысла.

Можно извлечь и другую информацию из диаграммы на рис. 4.6.1. Сравним значения [Эл/Н], полученные по линиям различных нейтральных и ионизованных элементов. Если все эти значения одинаковы, то все точки на диаграммах типа рис. 4.6.1 будут попадать на одну прямую линию.

Например, различия в определении [Эл/Н] по линиям ней­ трального и ионизованного циркония проявятся в том, что точки не лягут на принятую линию.

4.7. ОЦЕНКА ЭЛЕКТРОННЫХ КОНЦЕНТРАЦИЙ ПО ЛИНИЯМ ВОДОРОДА

Когда водород находится главным образом в ионизованном состоянии, число No,2 атомов, которые способны поглощать бальмеровские фотоны, зависит как от температуры, так и от электронного давления. Учитывая, что + N+ та N+, получим

lg (No. Ж ) = 3,39Ѳ + lg Ре- 2,5 lg Т + 0,78.

(4.7.1)

Поскольку Ѳ.<; 0,5, величина lg Л, играет важную роль при оп­ ределении относительного содержания N0,2.

«

Когда

водород преимущественно нейтрален, т. е.

+ N+ л;

lg (N0j2/N0) = 0,30 — 10,200.

(4.7.2)

~

Здесь Ѳ>, 0,5, а электронное давление вовсе не входит в фор­ мулу. Поскольку при величине Ѳ стоит большой численный ко­ эффициент, относительное содержание /Ѵ0,ч сильно зависит от температуры.

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

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

Самый простой метод оценки электронного давления по ли­ ниям Бальмера можно осуществить по формуле Инглиса — Теллера (см. приложение I). Если пщ — главное квантовое

134 ГЛАВА 4

число последней бальмеровской линии, которая еще отчетливо

видна, то

(4.7.3)

lg 2Ne 23,26 — 7,5 lg пт.

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

Крылья более низких членов бальмеровской серии также дают возможность оценивать lgPe, если есть оценка N0i2H, где Н — толщина обращающего слоя.

Приведем остроумный метод опре­ деления N0'2Н> предложенный Ун-

 

зольдом [156].

 

члены серии

 

Самые высокие

 

Бальмера настолько сильно ушире­

 

ны штарковским эффектом, что их

 

глубины в центре малы. На языке

 

модели

Шустера — Шварцшильда:

Рис. 4.7.1. Определение N0,2H

обращающий слой

оптически

то­

(схематический график).

нок на

длинах

волн

наиболее

вы­

 

соких

членов

серии Бальмера.

Поэтому эти линии можно рассматривать как слабые, откуда следует

Г л = (леУтс2) l 2N0,2f2nH.

(4.7.4)

Силы осцилляторов для линий водорода хорошо известны, поэтому измерения Wy, давали бы N0t2H для самых высоких чле­ нов серии Бальмера, если бы крылья не накладывались друг на друга. Унзольд следующим образом попытался учесть эти наложения. Наносим на график величины lgNQ2H, определен­ ные из соотношения (4,7.4), в функции п для верхнего уровня линии Бальмера. lg N0t2H сначала будет расти, так как линии на самом деле становятся слабыми; затем наступит максимум, после которого lg M0t2H будет падать, поскольку измеряются слишком малые и происходит наложение крыльев линий. Если максимум достаточно плоский, то можно провести экстра­ поляцию к п — оо и определить N0t2H. Такая кривая схематиче­ ски показана на рис. 4.7.1.

Теперь предположим, что контуры линий водорода можно описать формулой

где

Ч = Ч Н>

(4.7.6)

КОЛИЧЕСТВЕННЫЙ ХИМИЧЕСКИЙ АНАЛИЗ ЗВЕЗДНЫХ АТМОСФЕР 135

а максимальная глубина в центре берется из наблюдений. Интерполяционная формула (4.7.5) широко использовалась (ср. [158, § 103]) и может быть пригодной как для контуров, образо­ ванных чистым поглощением, так и для контуров в случае рас­

сеяния. Для

<

1

результат R% ~ х% следует

из модели

чи­

стого

поглощения

и

Шустера — Шварцшильда,

описанной

в

разд.

2.8 при ц = 1

— Rx- Теперь

 

 

 

 

 

 

xx = axN0i2H,

(4.7.7)

где ал — коэффициент поглощения на частицу. Форма коэффи­ циента поглощения на частицу исследована в гл. 5 и 6. В дан­ ном случае достаточно ограничиться простым приближением «самой близкой возмущающей частицы» (разд. 5.4) и получить для крыльев линии*)

ах = (3/2) W l m c 2) llf (А Л /г (Ак )~ \

(4.7.8)

Здесь AA.0 — смещение спектральной линии в электрическом поле

с напряженностью <§ь

 

&0 = е (4яЛ73)ѵ\

(4.7.9)

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

Для ДЯо мы используем приближение (ср. разд. 6.5 или [158])

\ l Q= (3hl8n2cme)l2nk&0 = 0,0l92l2nkg Q (ед. СГС), (4.7.10)

где я* — разность между

квантовыми числами верхнего и ниж­

него состояний, дающих

штарковские компоненты, а йа— сред­

нее значение для чисел «а одной водородной линии, йа равно

2,24; 5,96; 11,8 и 15,9 соответственно для На, Hß,

Ну и Н6.

После подстановки численных значений получим в системе

СГС

 

ссх= 1,56 • 10-2%1(пк)ъ (2Ые) A r t

(4.7.11)

Унзольд [158, табл. 51] приводит удобную таблицу сил осцилля­ торов для линий водорода, усредненных по интенсивностям от­ дельных компонент.

*) Нормировка контура поглощения для «наиближайшей возмущающей частицы» выбрана так, чтобы

О

(* (*х dX = (яе2/отс2) Agf.

о

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

ß = АЯ/ДЯо.

136

ГЛАВА 4

 

Если объединить уравнения (4.7.11) и (4.7.7) с результатами

экстраполяции для

отыскания

то можно получить при­

ближенную электронную концентрацию Ne.

Представленная здесь схема позволяет определить абсолют­ ную величину Ne. Очевидно, если в распоряжении имеются дан­ ные для двух звезд, то можно найти [А7е]. Нет нужды приводить подробные формулы.

4.8. ОПРЕДЕЛЕНИЕ ОТНОСИТЕЛЬНОГО СОДЕРЖАНИЯ ЭЛЕМЕНТОВ

При помощи уравнения (4.4.5) можно записать

 

[#] = [*] + ДѲх„ + Ы + К ]

(4.8.1)

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

Содержание элемента можно найти для каждой наблюдае­ мой стадии ионизации по формуле Саха. Пусть x —N\!{N\-\-N0), а отношение (А71 + N0)/N,n обозначено через (Эл/Н), где Na — концентрация атомов и ионов водорода. Тогда из (4.8.1) для нейтральных атомов получаем

[Эл/Н] = [X] + ДѲ*„ + [и] + [KifNH] - [ 1 - х ] ,

(4.8.2)

а для линий, образованных ионами, имеем

 

[Эл/Н] = [X] + ДѲХ„ + [и] + [xi/NH] - [х].

(4.8.3)

Обобщение формул на элементы в r-й и + 1)-й стадиях иони­ зации очевидно, и нет нужды приводить эти формулы отдельно.

Чтобы вычислить X для обеих звезд, необходимо задаться

абсолютными величинами Ѳ и 1gP e для одной звезды. Затем, используя разности логарифмов ДѲ и А lg Ре, можно определить

соответствующие

величины для

второй звезды. Вообще [я] и

[1— X] не будут

чувствительны к малым ошибкам в оценках Ѳ

и lg Ре, и таким образом дифференциальный характер

анализа

сохраняется.

 

делят на Nn потому,

что для

Коэффициент поглощения

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

 

< =

^ (Н _)а (Н _) +

іѴна(Н),

(4.8.4)

где а(Н~)

и а (Н )— коэффициенты

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

на частицу

для Н~

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

водорода

соответственно.

КОЛИЧЕСТВЕННЫЙ ХИМИЧЕСКИЙ АНАЛИЗ ЗВЕЗДНЫХ АТМОСФЕР 137

Концентрация отрицательных ионов водорода просто связана с числом Nn посредством уравнения Саха:

NHNe/N (H -)^K (T ),

(4.8.5)

и из (4.8.5) и (4.8.4) следует

 

 

nx/NH~f{P e’

атомные параметры).

(4.8.6)

Табличные значения коэффициента непрерывного поглоще­ ния в звездах со стандартным содержанием элементов вычис­ лены Боде [15, 16]. По этим таблицам и приближенным значе­ ниям Ѳ и Ре легко вычислить KcxfNn. Обычно при анализе ис­

пользуется одна средняя длина волны, но в [15, 16] приведены вычисления для ряда длин волн, так что в случае необходимо­ сти можно рассмотреть и зависимость от длины волны.

Таблицами непрозрачности можно пользоваться и в тех слу­ чаях, когда отношение содержаний водорода и металлов А от­ личается от стандартного, до тех пор пока ( ^ / jVh) не зависит

от ІѴн. Если основными источниками непрозрачности являются Н (связанно-свободные переходы), Н- и рэлеевское рассеяние на нейтральном водороде, то изменения А отражаются только

на электронном давлении. В таком случае отношение

в

бедных металлами звездах будет таким же, как и

в звездах

с высоким содержанием металлов при одинаковых электронных давлениях.

Величина kcx/N h зависит от Nн, когда существенна непро­ зрачность, вызванная металлами, или когда имеет место рэлеев­ ское рассеяние на молекулах Нг. Такие случаи должны рас­ сматриваться отдельно. Вероятно, лучший способ определения

— это взять одну из программ вычисления непрозрачности,

применяемых при расчете атмосферных моделей. Такие про­ граммы имеются почти во всех современных астрономических учреждениях.

До сих пор мы ничего не сказали о гелии. В большинстве звезд вклад гелия в непрозрачность атмосферы незначителен. Что же произойдет со спектром нормальной звезды при замене части атмосферного водорода гелием? Для конкретности возь­ мем солнечный спектр и предположим, что существующее от­ ношение Не/Н = В равно 0,15 по числу атомов. Каким бы стал солнечный спектр, если бы В увеличилось от 0,15 до 0,30? Пусть, для простоты, атомы водорода замещены атомами гелия таким образом, чтобы изменение давления было незначительным.

Поскольку гелий не вносит вклад в непрерывное поглощение, а интенсивность слабых линий зависит от отношения

то интенсивность линий поглощения должна возрасти. Такое из­ менение мы восприняли бы как признак понижения отношения

138 ГЛАВА 4

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

Общим результатом изменения отношения В на Солнце с 0,15 до 0,30 было бы усиление металлических линий и возраста­ ние среднего газового давления. К сожалению, изменения га­ зового давления не так легко обнаружить, как усиление метал­ лических линий в спектре. Мы не можем измерить [Pg] с такой

высокой точностью, как

или л/Н]. Таким образом,

большинство определений

относительного содержания дают

Эл/Н, и мы мало что можем сказать об отношении Эл/(Н + Не).

4.9. СПЕКТРАЛЬНЫЙ СИНТЕЗ: ОБЩИЕ ЗАМЕЧАНИЯ

Изучение звездных спектров с привлечением полных моде­ лей атмосфер называется тонким или детальным анализом. Унзольд [158] использовал слово feinanalyse (нем. — тонкий или детальный анализ) для описания детального рассмотрения на основании полной (учитывающей зависимость от глубины) мо­ дели атмосферы в противоположность grobanalyse (нем. — грубый анализ), когда применяется схематическая модель Шу­ стера — Шварцшильда или подобные ей модели. Метод спек­ трального синтеза основан на полных моделях тонкого анализа, но, кроме того, он должен удовлетворять и дополнительным требованиям.

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

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

КОЛИЧЕСТВЕННЫЙ ХИМИЧЕСКИЙ АНАЛИЗ ЗВЕЗДНЫХ АТМОСФЕР 139

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

Молекулярные особенности являются важными индикаторами как температуры, так и ускорения силы тяжести на поверхно­ сти звезды. Например, при низкой дисперсии появление полосы G — особенно чувствительный индикатор температуры, и было бы весьма полезно развить теоретический метод, который давал бы количественные значения свойства, описываемого в спект­ ральной классификации как появление полосы G. Ясно, что ис­ следователи, работающие с малой дисперсией, имеют ряд чув­ ствительных критериев температуры и яркости, которые пока не могут применяться исследователями, работающими с высо­ кими дисперсиями, и теоретиками.

Метод спектрального синтеза только начинает применяться (главным образом к Солнцу). Много полезных исследований, смягчающих некоторые требования для полного спектрального синтеза, уже сделано и может быть еще сделано. Но мы не должны считать анализ какого-либо спектра законченным, пока не выполнен полный спектральный синтез.

4.10.МОДЕЛЬ АТМОСФЕРЫ

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

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

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

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

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