книги из ГПНТБ / Каули, Ч. Теория звездных спектров
.pdf130 ГЛАВА 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 и химического состава. Исследователь химического со става, начинающий изучение звезды методом спектрального син теза, будет иметь примерное представление о величине этих
*) Конвективные модели могут вычисляться для нескольких предполагае мых значений длины пути перемешивания. Солнечные модели, учитывающие движения холодных и горячих элементов, имеют другие переменные пара метры. Эти лишние степени свободы редко используются в исследованиях звезд.
