CC BY
73
УДК 524.35-852/74; 520.22
СПЕКТРОСКОПИЯ ЗВЕЗД В НАЗЕМНОМ УЛЬТРАФИОЛЕТЕ.
II. АТМОСФЕРЫ
2010 В. Е. Панчук, В. Г. Клочкова, М. В. Юшкин
Специальная астрофизическая обсерватория, Нижний Архыз, 369167 Россия Поступила в редакцию 20 января 2009 г.; принята в печать 30 июля 2009 г.
В историческом плане рассмотрены некоторые задачи спектральных исследований звездных атмосфер в наземном ультрафиолете. В качестве примеров приведены фрагменты спектров, полученных авторами на 6-м телескопе БТА.
Ключевые слова: звезды — переменные и пекулярные; методы астрономических наблюдений, приборы и инструменты
1. ВВЕДЕНИЕ
Необходимость спектроскопических наблюдений звезд в УФ-диапазоне объясняется хотя бы тем, что у горячих звезд максимум распределения энергии в спектре лежит именно в этом диапазоне.
Наземный ультрафиолет (3000—3800 iA) чрезвычайно богат линиями атомов и ионов. Высокая плотность линий создает трудности их идентификации и проблему теоретического описания многочисленных спектральных бленд. При переходе в заатмосферный ультрафиолет плотность линий в спектрах холодных звезд возрастает настолько, что указанный диапазон становится малоинформативным.
Задачи звездной спектроскопии в наземном ультрафиолете можно разделить на три группы. В первой группе используется то обстоятельство, что условия формирования линий в звездной атмосфере существенно различаются для разных ультрафиолетовых поддиапазонов. Ко второй группе можно отнести все задачи, где необходимо наблюдать линии определенных химических элементов и/или избранных мультиплетов, не представленных в других диапазонах длин волн. К третьей группе относим задачи обнаружения и диагностики хромосфер холодных звезд, а также задачи спектроскопии околозвездных оболочек. Некоторые из задач перечислены в [1]. В данной работе рассмотрим задачи, принадлежащие первым двум группам.
2. ПРЕДЕЛ БАЛЬМЕРОВСКОЙ СЕРИИ ВОДОРОДА
Привлекательность использования характеристик абсорбционных линий высоких членов баль-
меровской серии состоит в том, что область формирования соответствующего излучения залегает относительно глубоко в атмосфере, т.е. слабо связана с нестабильностями внешних слоев, но хорошо отражает изменения эффективной температуры. Предполагалось и подтверждено расчетами, что в условиях частичной ионизации водорода влияние содержания других доноров электронов является несущественным. Таким образом, область бальмеровского скачка предоставляет нам “чистый инструмент” определения фундаментальных параметров звездных атмосфер, работающий однозначно в достаточно широком диапазоне концентраций металлов, легких элементов и гелия.
2.1. Классификация спектров
В дополнение к гарвардской и йеркской системам классификации спектров, в середине ХХ века была создана французская классификация, основанная на изучении непрерывного спектра в ближней ультрафиолетовой области. Первоначально измерялись два параметра: величина скачка (О) вблизи предела бальмеровской серии водорода и длина волны Аь на которой падение интенсивности достигает половины О. Основным достоинством классификации является ее нечувствительность к химическому составу. К достоинствам также отнесем использование диапазона, где угловая дисперсия призм быстро возрастает, к недостаткам — труднодоступность используемого диапазона для большинства спектрографов из-за использования в них тяжелых сортов стекол и низкого размещения над уровнем моря. При продвижении от ранних подклассов В к подклассу А7 величина бальмеровского скачка растет, затем снижается,
поэтому двумерная классификация является неоднозначной. Для устранения этой неоднозначности в трехмерной французской классификации (IAP) в качестве третьего параметра добавлен спектрофотометрический градиент в области 3800—4600 A,
при этом нечувствительность классификации к содержанию металлов была сохранена.
2.2. Электронные плотности
Привлекательность использования бальмеров-ского предела состояла еще и в том, что, кроме параметров, описывающих характер распределения энергии в спектре, удавалось определить среднюю электронную плотность. В формуле Инглиса-Теллера logpe = 7.4 — 7.5log n + log T электронное давление pe связано с температурой T и номером n последней различимой линии бальмеровской серии. Здесь уже требовалось, по крайней мере, получать спектры со средним спектральным разрешением и вводить поправки за эффект осевого вращения [2]. Методом кривых роста электронное давление удавалось определить, используя линии соседних стадий ионизации данного элемента, и возможность независимого определения электронной плотности ne методом Инглиса-Теллера была принципиальной. Номер линии n — целое число, и с учетом коэффициента 7.5 при log n точность определения электронной плотности была невысокой. На Рис. 1 приводим примеры, иллюстрирующие возможности метода Инглиса-Теллера при использовании высокого спектрального разрешения.
2.3. Метод бальмеровских индексов
Переход от средних характеристик атмосфер к одномерным моделям атмосфер на первый взгляд исключил метод Инглиса-Теллера из практики спектроскопистов. Однако критерий Инглиса-Теллера еще долгое время неявно использовался, например, в алгоритмах расчета моделей звездных атмосфер [3] при определении номера последней учитываемой линии бальмеровской серии. На кильской диаграмме (температура Te — ускорение силы тяжести log g) все изолинии равенства наблюденных и теоретических характеристик водородных линий в диапазоне температур 22000—12000K расположены наклонно, и только изолинии показателей цвета (U — B)0, (B — V)о и бальмеровского скачка D расположены вертикально, см., например, [4]. Это означает, что при неопределенностях в оценке величины межзвездного или околозвездного покраснения единственным параметром для В-звезд, слабо чувствительным к log g, остается величина бальмеровского скачка. Измерение величины скачка
О на спектрограммах высокого разрешения представляет определенные трудности, связанные с экстраполяцией континуума. Поэтому в [4] был предложен и реализован метод бальмеровских индексов, позволяющий расширить число температурных критериев при определении параметров атмосфер В- и А-звезд. При этом использовался достаточно узкий спектральный интервал (200 ), т.е. метод практически свободен от эффекта покраснения. По наблюдениям 73 звезд с дисперсией 28 Л/мм выяснилось, что эффективные температуры, определенные путем сравнения наблюденных и теоретических бальмеровских индексов, систематически завышены относительно определений в общепринятой шкале в среднем на 1000 К, а эффективные ускорения силы тяжести занижены относительно определений в [2] в среднем на Дк^д = 0.5. Было сделано предположение, что это расхождение связано с приближенным характером расчета области бальмеровского скачка в моделях [3], где для определения номера последней учитываемой линии п использовался критерий
— 2/15
Инглиса-Теллера п = 1100пе . Итак, если
ранее по номеру последней наблюдаемой линии оценивалась электронная плотность, то в методе моделей атмосфер [3] ситуация стала обратной: электронная плотность пе на каждом уровне модели вычислялась, а номер п определялся только из условия равенства квазистатической ширины линии расстоянию до следующей линии. Вопрос о последней линии в серии является задачей определения самого высокого связанного состояния атома в плазме: внутриплазменное электрическое микрополе размывает электронные уровни и снижает потенциал ионизации атома. В [5] показано, что последнее реализуемое энергетическое состояние атома определяется номером
п = 3800пе 1/6, учет туннельного эффекта снижает п на единицу, а влиянием электронных соударений при плотностях атмосфер В- и А-звезд можно пренебречь. С учетом этих обстоятельств в [5] был выполнен перерасчет моделей [3] и сетки бальмеровских индексов. Используя полученные на 0.6-м и 6-м телескопах 118 спектрограмм для 57 звезд, авторы [6] пришли к выводу, что новая сетка бальмеровских индексов соответствует температурной калибровке [7]. Внутренняя точность фотографического метода оказалась высокой (ДТе = 200 К, Дк^д = 0.1), и возрастает с переходом на фотоэлектрические приемники. Мы считаем, что, если использовать характеристики только водородного спектра, то метод бальме-ровских индексов является основным свободным от эффектов покраснения спектроскопическим методом определения температур В-звезд. При
Wavelength, А
Wavelength, А
Рис. 1. (а) Бальмеровская серия водорода в спектрах звезд разной светимости А2 1а (вверху) и А0 Уа (внизу). В приводимом участке спектра сверхгиганта видны линии от Ніб 3712 А до Н27 3666 А, в спектре карлика — до Н21 3679 А. По ординате указано число отсчетов на канал (1л). (Ь) Область предела бальмеровской серии водорода в спектрах звезд А2 1а (вверху) и А0 Уа (внизу). В спектре сверхгиганта видны линии до Нзз 3658 А. Спектры получены на НЭС БТА.
СПЕКТРОСКОПИЯ ЗВЕЗД В НАЗЕМНОМ УЛЬТРАФИОЛЕТЕ I/Icont
Рис. 2. Наблюдения ШЕ (тонкая линия) и модельный расчет резонансного дублета С IV 1548 А и 1550 А в спектре 27 СМа, В3 ^е. Приняты величины угла наклона г = 90° и угловой скорости осевого вращения ш = 0.2 (в долях критической скорости). Наблюдаемая асимметрия крыльев дублета не может быть объяснена фотосферными эффектами.
использовании данных, аналогичных приведенным на Рис. 1, точность метода возрастет.
3. ВЛИЯНИЕ ОСЕВОГО ВРАЩЕНИЯ НА СПЕКТР ЗВЕЗДЫ
Основа “широкодиапазонного” подхода к проблеме осевого вращения заложена в простой модели вращающейся сферы однородной яркости [8]. Уже в таком приближении величина уширения изменяется вдоль спектра, т.к. полуширина линии пропорциональна длине волны. Следовательно, при переходе от синего диапазона к ультрафиолетовому эффект может быть обнаружен. Попытки вычислить профили линий, учитывая сплюснутость вращающейся звезды и зависимость физических условий от координат на поверхности звезды, начались сразу, как только это позволило состояние вычислительных средств. В [9] были вычислены профили линий HeI 4471 A, MgII 4481 iA и Call 3934 A с учетом дифференциального вращения. Для выявления эффекта точность фотографических наблюдений оказалась недостаточной. В [10] проведен расчет влияния вращения на линии HeI 4471 A, FeI 4476 A, MgII 4481 A в диапазоне спектральных классов 09—F8. Профили первых двух линий оказались чувствительными к наклону оси вращения, третьей — нет.
Ситуация оживилась после миссии Коперник: в [11] отмечено, что в заатмосферном УФ (1170—1400 A) проекция скорости осевого вращения Z Tau составляет 100 км/с, тогда как по спектрам в видимой области — 300 км/с. Условие ионизационного баланса как для SiIII и SiIV, так и для CIII и CIV дает параметры атмосферы (Te = 27500 K, log g = 4.0), т.е. температуру гораздо более высокую, чем следует из классификации по видимому диапазону (B4 III, что соответствует температуре 17500 K). Практически одновременно были опубликованы модельные расчеты [12], выполненные при следующих упрощениях: сферическая звезда, планковское (чернотельное) распределение энергии, закон гравитационного потемнения по Цейпелю, доплеровский профиль, отсутствие дифференциального вращения. Расчеты выполнены для 4500 iA, 1200 iA и 1000 iA. Оказалось, что в УФ-диапазоне профили имеют более резкие ядра и неглубокие крылья, что интерпретировалось как результат большего градиента потока в УФ, поэтому для холодных В-звезд эффект выражен сильнее. Отношение величин V sin i, определенных по спектру в видимом и УФ-диапазоне, не сильно зависит от угла наклона i и уменьшается, когда уменьшается угловая скорость и, но увеличивается при понижении температуры. Затем были использованы [13] оригинальные наблюдения с КА Коперник в диапазоне
1000—1300 A с разрешением 0.2 A. Для линий
S IV 1062 A, S IV 1073 A, C II 1092 A, Si III 1108 A, III 1110 A, Si IV 1122 A, C III 1248 A проведена процедура разделения V и i. Исследования были продолжены [14] с привлечением новых наземных наблюдений. В частности, получено, что Ве-звезды не вращаются с критической скоростью и имеют случайное распределение осей вращения.
Многие линии в ультрафиолете являются резонансными, их использование для изучения эффекта уширения вращением представляется проблематичным. В качестве примера на Рис. 2 приводим фрагмент рисунка из работы [15], в которой по наблюдениям на IUE выполнено определение скорости вращения и угла наклона i для 7 звезд. Использованные при этом резонансные линии зачастую имеют протяженные асимметричные коротковолновые крылья. Моделирование эффектов вращения для нерезонансных линий: триплета Si III 1113.17 A, 1113.20 A, 1113.23 A и линии MgII 4481 A, выполнено в [16], где для спектральных подклассов В1 и В5 показано, что изменения с длиной волны вращательной полуширины имеют высокую чувствительность к массе звезды. В [17] по линиям HeI 4144 A, HeI 4471 A, триплета
SiII 4128-4131 A, FeI 4476 A, MgII 4481 A исследовано влияние вращения на результаты спектральной классификации. Теоретические расчеты для линий HeI 4471 iA и MgII 4481 iA [18] и результаты наблюдений с дисперсией 20 iA/мм были использованы для определения углов наклона и скорости вращения 69 звезд В0-В3 [19].
В результате моделирования атмосфер быстро-вращающихся звед был сделан вывод о том, что наблюдаемые скорости осевого вращения горячих звезд должны зависеть от длины волны используемой линии. Причина различий V sin i следующая. У быстровращающейся звезды в непрерывном спектре преобладает излучение, выходящее из более горячих околополярных областей. При этом доля вклада излучения, формирующегося в око-лополярных областях, в суммарный спектр звезды возрастает с уменьшением длины волны. Наблюдаемый профиль спектральной линии является суперпозицией профилей линии, формирующихся в различных участках поверхности звезды, обращенной к наблюдателю. Вклад околополярных областей, где дисперсия скорости меньше, в наблюдаемый профиль линии также будет увеличиваться с уменьшением длины волны. Следовательно, профили спектральных линий становятся уже при продвижении в УФ. Физической основой методов разделения V и i является различие в гравитационном потемнении в континууме между далеким УФ (1100 A) и видимым диапазоном. Подчеркнем,
что моделирование эффекта вращения основано на предположении, что сильные линии, наблюдаемые в заатмосферном ультрафиолете, формируются в атмосфере звезды и околозвездная оболочка не вносит заметного вклада в профили этих линий высокой стадии ионизации. Спектроскописты надеялись развернуть работы по массовому определению пространственной ориентации осей вращения. Однако, вскоре появились доказательства оболочечного происхождения ряда линий, или, по крайней мере, наличия существенного вклада оболочки в профили фотосферных линий (см. Рис. 2). Уже в [16] было показано, что влияние вращения слабее, чем это наблюдается [11, 20], или чем это предсказано в упрощенных моделях [12]. Расчеты, проведенные для бланкетированных моделей, в сравнении с данными IUE показали [21], что различия в скорости осевого вращения, полученной в УФ и видимой области для В-звезд, не так велики, как указывалось в пионерских работах, а в ряде случаев находятся в пределах ошибок измерений. Вычисления синтетических спектров показывают, что в заатмосферном ультрафиолете сочетание высокой плотности линий с быстрым вращением приводит к ситуации, когда в широких спектральных интервалах отсутствуют участки континуума (см., например, [22], Рис. 1). Блендирующий эффект большого числа линий, в сочетании со свойствами регистрирующей аппаратуры (сложная зависимость концентрации света в коротких спектральных порядках эшелле-спектрометров IUE), затрудняет определение уровня непрерывного спектра и, следовательно, вносит ошибки в “ультрафиолетовые” величины V sin i.
Мы считаем, что для задачи разделения V и i оптимальным является сочетание наблюдений в видимой области с наблюдениями в наземном ультрафиолете. В наземном ультрафиолете неопределенность проведения уровня непрерывного спектра меньше и точность выделения профилей линий выше. Набор используемых спектральных линий следует дополнять относительно слабыми линиями, нечувствительными к эффектам нестабильности внешних слоев атмосферы. На Рис. 3 по данным наших наблюдений сравнены профили линий, расположенных до и после бальмеровского скачка, в спектре Веги — быстровращающейся звезды, наблюдаемой почти с полюса (V = 175 км/с,
i = 7° [23]). Видно, что профили линий одинаковой интенсивности, сформированных на различной глубине (до и после бальмеровского скачка) по-разному реагируют на осевое вращение. Таким образом, задача определения угла наклона оси вращения, поставленная ранее для сильных линий, может быть решена по слабым линиям. Однако это требует регистрации спектров звезд с высоким отношением сигнал/шум.
г
л
-40 -20 0 20 40
V, кт/Б
Рис. 3. Профили линий одинаковой остаточной интенсивности в спектре а Ьуг, расположенные до (сплошная линия, около 4300 Л) и после (пунктир, около 3300 Л) предела бальмеровской серии. Спектр получен на НЭС БТА.
4. ХИМИЧЕСКИЙ СОСТАВ АТМОСФЕР
Исследования в заатмосферном ультрафиолете дали много новой информации о содержании химических элементов в межзвездной среде и относительно мало о химическом составе звездных атмосфер. Были обнаружены линии ионов, для объяснения природы которых требовались необычно горячие атмосферы. Так появилась концепция горячей околозвездной оболочки [24], которая искажает сильные резонансные линии ветровыми и оболочечными эффектами; значительную роль также играет блендирование. Для задачи определения содержания элементов в звездных атмосферах диапазон 1200—2900 А оказывается менее пригодным и более дорогостоящим, чем наземный ультрафиолет. Исключение составляют элементы, не представленные линиями в наземном ультрафиолете, например, бор (резонансный дублет В1 2496.7 А и 2497.7 А).
4.1. Содержание элементов CNO-группы и металлов
В процессе классификационных обзоров были выделены ОВ-звезды с пекулярными соотношениями интенсивностей линий элементов СЫО [25]. В дальнейшем была поставлена наблюдательная программа, целью которой являлось как изучение химического состава атмосфер таких звезд, так и поиск неоднородностей химического состава при исследовании горячих звезд в различных
скоплениях и группировках. Был построен специализированный спектрограф [26], параметры которого подобраны так, чтобы линии, используемые для определения химического состава (CII 4267 A,
NII 3995 A, 4630 A, OII 4638 A, 4641 A, 4649 A), находились в центральных частях порядков эшель-ного спектра, а параметры решетки эшелле выбраны так, чтобы изменение концентрации света в порядках было менее крутым. Проблему составляет определение содержания азота, т.к. линии NII слабы (W < 100м^1), и при больших Vsini требуется высокое отношение сигнал/шум (S/N) для их обнаружения. Программа наблюдений была выполнена только для ближайших скоплений, т.о., задача не потеряла актуальности до сих пор. В дальнейшем систематика спектров ОВ-звезд [27] проводилась как по видимой области (фотосфер-ные линии), так и в далеком ультрафиолете (ветровые профили), а наблюдения в наземном ультрафиолете по-прежнему практически не проводились.
Радиальный градиент химического состава в Галактике традиционно измерялся по HII-областям
и, в последнее время, по спектрам цефеид, но в ряде случаев наблюдениям доступны только яркие массивные звезды верхней части ГП. Видимая область спектров горячих звезд бедна линиями ионов металлов. Смещение наблюдаемого диапазона из синей области в фиолетовую и ультрафиолетовую расширяет число доступных линий в несколько раз и делает определения химического состава более уверенными. Разработанный нами вариант на-
блюдений на НЭС БТА позволяет регистрировать диапазон 3000—4400 А с величиной спектрального разрешения К = 60000 с перекрытием соседних порядков во всем диапазоне (см. [1], Рис. 2). На Рис. 4 приведен небольшой фрагмент спектра горячей звезды с линиями металлов.
4.2. Содержание бериллия
Резонансный дублет Ве11 3130.42 А и 3131.06 А находится в области, где коэффициент пропускания земной атмосферы в шесть раз ниже, чем на длине волны 4000 А. На Рис. 5 приведен полученный на НЭС БТА фрагмент спектрального порядка, содержащий линии Ве11. Отметим три задачи, связанные с определением содержания бериллия в атмосферах звезд разной металличности (возраста), подробнее см. в [28].
Задача 1. Изучение конвективного перемешивания и оценка дозвездного содержания Be. Триада литий/бериллий/бор интересна тем, что их содержание (по отношению к водороду) в космических лучах и в звездных атмосферах различается на несколько порядков. Это обстоятельство явилось эмпирической основой сценария [29], в котором изотопы Li, Ве, В образуются в реакциях скалывания при столкновениях галактических космических лучей (ГКЛ) с ядрами элементов СЫО-группы, выброшенными в межзвездную среду в результате эволюции звезд предыдущих поколений. Литий наиболее доступен наблюдениям, но и самый хрупкий. При конвективном перемешивании оболочки звезды ядра лития разрушаются, если внутренняя граница конвективной зоны достигает слоев с температурой 2.5 млн К. В спектре Солнца линия LiI 6707 А не видна, следовательно, конвективное перемешивание достигает слоев с указанной температурой. Ядра Ве разрушаются при температуре 3.5 млн К. В спектре Солнца и звезд, подобных Солнцу, линии бериллия наблюдаются (Рис. 5). В этом случае, по-видимому, нижняя граница конвективной зоны не достигает слоев с температурой, достаточной для разрушения бериллия. Известны объекты, где содержание Li уже уменьшено, а содержание Ве еще осталось дозвездным. Эти факты послужили основой для наблюдательных задач по тестированию механизмов перемешивания в оболочках холодных звезд разного возраста (конвективный перенос или перемешивание, индуцированное вращением, см., например, [30]). При понижении эффективной температуры от 6000 К до 5000 К содержание лития в атмосферах звезд 1-го типа населения понижается на три порядка [31, 32], тогда как содержание бериллия остается постоянным (в пределах точности первых наблюдений, [33]) и равным содержанию, наблюдаемому в атмосферах
горячих звезд и Солнца. Выбор механизма истощения легких ядер связан и с проблемой “литиевого провала” — в узком интервале температур наблюдается резкое понижение содержания лития, сменяемое столь же резким возвратом к прежнему значению [34]. В области литиевого провала обнаружена корреляция содержания бериллия и лития [35], интерпретируемая в рамках моделей перемешивания, индуцированного вращением. Эффект исследован для трех рассеянных скоплений, обнаружены различия от скопления к скоплению, поэтому предполагается, что степень истощения лития и бериллия может быть связана с возрастом и металличностью [36].
Задача 2. Оценка количества бериллия космологического происхождения. Если наблюдать горячие маломассивные субкарлики (Те > 5500 К) со слабой внешней конвективной зоной, то по содержанию лития и бериллия в их атмосферах можно определить дозвездное содержание этих элементов. Как только состоялся переход на низ-кошумящие фотоэлектрические приемники, было обнаружено, что у маломассивных субкарликов концентрация ядер лития сначала понижается с
понижением металличности ^е/Н] (т.е. с увеличением возраста), а затем остается неизменной, формируя т.н. литиевое плато [37]. Такое остаточное содержание лития, рассматриваемое как синтезированное в эпоху Большого Взрыва и сохранившееся в атмосферах медленно эволюционирующих малометалличных звезд, соответствовало оценке барионной плотности (в стандартной модели), совпадающей с оценкой, сделанной по содержанию дейтерия и гелия. Были предприняты попытки [38—40] наблюдений резонансного дублета Ве11 в спектрах звезд все более низкой метал-личности в надежде выйти на плато, аналогичное литиевому. Оказалось, что в отличие от Li при низких металличностях содержание Ве не выходит на плато, а монотонно падает с увеличением возраста звезды. Итак, в атмосфере одного и того же субкарлика (возраст которого, в пределах точности, совпадает с возрастом Вселенной) наблюдается “космологическое” содержание лития, а содержание бериллия 1сg{N(Бе)/М(Н)} = —13 на три порядка выше, чем дает стандартная космологическая модель. Попытки “свести концы с концами” предпринимались неоднократно. Во-первых, совершенствовались модели синтеза элементов во время Большого Взрыва. Например, расширение набора реакций первичного нуклеосинтеза (учет столкновений ядер лития с тритием, [41]) дает выход 1сg{Бe/Н} = —13, что на три порядка выше выхода в стандартной однородной модели. Предпринимались также попытки получить значимый
Wavelength, А
Рис. 4. Фрагмент спектра горячей звезды в Ori, B8 Ia, содержащий линии ионизованных металлов. Спектр получен на НЭС БТА.
Wavelength, А
Рис. Б. Фрагмент спектрального порядка (G2V, солнечный химсостав), содержащий линии резонансного дублета Be II (отмечены черточками). Спектр получен на НЭС БТА.
выход бериллия и бора с использованием неоднородной модели [42]. Во-вторых, уточнялись модели синтеза лития, бериллия и бора в реакциях на ГКЛ.
Задача 3. Определение спектра ГКЛ в разные эпохи синтеза легких элементов. По данным наблюдений 16-ти звезд с разрешением
К = 11000 и в предположении, что весь бериллий образован в реакциях скалывания на ГКЛ, в [38] получено, что таким образом можно синтезировать не более 10% лития, наблюдаемого в звездах гало. После того, как спектральное разрешение было
увеличено до К = 20000, в [39] был сделан вывод, что в реакциях на ГКЛ может быть синтезировано не более 1 % ядер лития. Итак, наблюдения бериллия укрепили мнение о том, что у субкарликов наблюдается литий космологического происхождения. Оставалось уточнить модель синтеза бериллия в ГКЛ. В [40] сделано предположение, что для производства наблюдаемого количества ядер бериллия в реакциях скалывания, требуется спектр ГКЛ в ранней Галактике, отличный от современного в области высоких энергий: много больше чем 100 МэВ на нуклон (синтез Ве в реакциях скалывания идет, начиная с 10 МэВ на нуклон).
Наблюдения звезд с экстремальным дефицитом элементов железного пика показывают, что содержание бериллия в них продолжает снижаться [43]: при ^е/Н]= —3.30 [Ве/н]< —13.4. По результатам определения содержания Ве в диапазоне трех порядков металличности можно уточнить, по какой схеме проходят реакции скалывания. Предположим, что реакции скалывания происходят преимущественно в окрестностях сверхновых, т.е. выход ядер Ве определяется концентрацией мишеней (ядер кислорода) и набегающих протонов. В этом случае выход Ве не зависит от содержания тех ядер кислорода, которые синтезированы в предыдущей волне звездной эволюции и уже находятся в межзвездной среде (МЗС). Содержание Ве, образовавшегося в окрестностях сверхновых, должно быть пропорционально содержанию свежесинте-зированного кислорода. Если реакции скалывания идут преимущественно в МЗС, то выход Ве зависит от концентрации ядер кислорода и плотности долгоживущих ГКЛ, т.е. от величин, определяемых числом сверхновых в единицу времени. Поэтому во втором случае выход Ве пропорционален квадрату концентрации ядер кислорода. Вначале из наблюдений получалось, что характер соотношения концентраций Ве и О изменяется с возрастом. Затем, с увеличением точности определения содержания, оказалось, что содержание ядер Ве линейно зависит от содержания кислорода. Здесь важную роль сыграло увеличение числа определений содержания кислорода и повышение их точности [44].
Сценарии образования легких элементов в эпоху ранней Галактики можно уточнить, добавив к определениям содержания бериллия и лития в спектрах звезд разной металличности определения содержания бора В. Считалось, что образование В и Ве в реакциях скалывания на ГКЛ модельно независимо, Ы(В)/Ы(Ве) ^7 [45]. Наблюдения [46,47] линии бора (В1 2496.7 А) на HST дали соотношения содержаний ядер Ы(В)/Ы(Ве), согласующиеся с модельными.
4.3. Звезды с дефицитом металлов
Расчеты химической эволюции ранней Галактики дают упрощенно следующую шкалу соотношения металличности и возраста: величина металлич-ности ^е/Н]= —3 достигается при возрасте Галактики 107 лет, ^е/Н]= —2 — при возрасте 108 лет, ^е/Н]= —1 — при возрасте 109 лет. Предметом проверки этого соотношения являются маломассивные, очень медленно эволюционирующие звезды, “заморозившие” в своих атмосферах химический состав дозвездного вещества. С уменьшением содержания элементов железного пика атмосфера звезды становится все более прозрачной в ультрафиолетовом диапазоне. На Рис. 6 сопоставлены фрагменты спектров четырех звезд разной ме-талличности, исследованных в рамках программы спектроскопии старых звезд в широком диапазоне длин волн [48], выполненной на спектрографе НЭС 6-метрового телескопа БТА. По мере того, как ультрафиолетовые спектры освобождаются от линий элементов железного пика, становятся доступными линии химических элементов, содержание которых изменяется с возрастом иначе, чем содержание железа ^е/Н]. Приведем несколько примеров.
В спектре красного гиганта CD—38°245 экстремально низкой металличности (^е/Н]= —4.5) в области 3356—3376 А обнаружены полосы ЫН, что позволило определить содержание азота [49]. В работе [50] по 20 линиям ОН в области 3138—3155 А определено содержание кислорода в атмосферах 9 малометалличных звезд. Этот подход получил развитие в [51], но после того, как была показана возможность использования ИК-триплета кислорода [44], приоритет в данной задаче был отдан более доступному ИК-диапазону. В спектрах малометалличных звезд при определениях химического состава удается использовать детали полосы СЫ в области 3880—3884 А [52]. В [53] выполнено первое определение содержания элемента из интервала атомных номеров 41<Z<55 для звезд гало: по резонансным линиям AgI 3280.7 А и 3382.9 А определено содержание серебра для 7 звезд гало.
Отношение содержаний ядер, синтезируемых преимущественно в процессах быстрой (г) и медленной (б) нейтронизации, можно уточнить за счет расширения числа используемых химических элементов, если привлечь наблюдения в фиолетовой и ультрафиолетовой области. Кроме того, появляется возможность датировки вещества методом, использующим разные характерные времена распада радиоактивных изотопов. Линия ТЫ1 4019.13 А изотопа 232^ (период полураспада 14 млрд. лет) в комбинации с линией стабильных ядер ШП 4018.82 А впервые использована для определения возраста в [54], при этом предполагалось,
Рис. 6. Фрагмент атласа [48] спектров звезд разной металличности: HD115444, |Те/И]= —2.91; 037—26, |Те/И]= —2.04; HD188510, ^е/Н]= —1.52; и 027—44, ^е/Н]= —0.60. Спектры получены на НЭС БТА.
что темп производства ядер Th и Nd неизменен во времени. Возраст, получаемый при таком предположении, оказывался зависимым от принятой модели химической эволюции Галактики и принятого отношения Th/Eu. Дополнительные наблюдения линий тория и европия [55] ослабили позиции Th/Nd-метода. Следующий заход на проблему был выполнен с применением более мощных наблюдательных средств [56, 57], поддержанных детальными расчетами г-процесса [57]. Такие наблюдения позволяют разделить индивидуальные события г-процесса, наблюдаемые на шкале более короткой, чем время перемешивания в Галактике. Вид классической пары линий Th/Nd меняется от звезды к звезде [58], что заставляет искать дополнительные линии, блендирующие линию тория. На Рис. 7 приведены фрагменты спектров в окрестностях линий ThII 4086.4/А, 4019.1A, 3929.7A и 3539.6A по данным [59]. Линия урана UII 3859.5 А может быть отделена от блендирующих линий CN, но из-за относительно небольшого периода полураспада 238U (4.47 млрд. лет) должна быть сильно ослаблена (или не наблюдаема) в атмосферах самых старых
звезд [59]. Таким образом, в задаче радиоактивной хронометрии торий остается наиболее удобным элементом.
5. ГРАДИЕНТЫ В АТМОСФЕРЕ
На Рис. 8 приведен пример расчетов [60] изменения коэффициента непрерывного поглощения и рассеяния с длиной волны для набора значений обратной температуры 0 = 5040/Т. Из рисунка следует, что в области наземного ультрафиолета геометрическая глубина формирования непрерывного излучения должна быть различной для разных длин волн. Следовательно, формирование абсорбций в разных длинах волн происходит также на разных геометрических глубинах. Оценки, выполненные для несерой модели атмосферы с Те = 10000 К, дают разность глубин формирования спектральных линий, находящихся до и после бальмеровско-го скачка, порядка 1000 км. На разности уровней формирования различных спектральных линий основаны попытки поиска градиентов различных характеристик в звездной атмосфере. Наблюдая линии, находящиеся в разных участках наземного ультрафиолета (например, до и после предела
бальмеровской серии водорода), можно получить некоторые представления относительно градиента поля скоростей, вертикального градиента магнитного поля, градиента концентрации различных частиц (в предположении диффузии элементов в радиальном направлении). Первая попытка поиска изменения продольной компоненты магнитного поля с геометрической глубиной в атмосфере была предпринята в [61], но точности фотографической спектрофотометрии оказалось недостаточно. Результаты измерений изменения магнитного поля с глубиной, выполненные в районе бальмеровского скачка с использованием современных средств, изложены в [62].
Другим способом построения эмпирического разреза атмосферы является исследование линий, формирующихся на различных участках крыльев сильных резонансных линий. В УФ диапазоне такой случай предоставляет дублет H и K Call.
Об использовании этого дублета будет упомянуто в следующей публикации.
6. СПЕКТРАЛЬНАЯ ПЕРЕМЕННОСТЬ
На Рис. 8 видно, что при понижении температуры коэффициент непрерывного поглощения в видимом и близком ИК-диапазонах, где у холодных звезд происходит перенос основной доли лучистой энергии, понижается. Объясняется это тем, что коэффициент поглощения на отрицательных ионах водорода Н_ пропорционален электронной плотности, а при низких температурах основной вклад в электронную плотность дают только легко ионизуемые щелочные металлы. При низких значениях непрозрачности резко возрастает газовое давление. Например, при переходе от АО к 05 газовое давление на уровне фотосферы увеличивается на два порядка, от 103 до 105 дин/см2. Поэтому в переносе энергии повышается роль аэродинамических механизмов, в том числе и крупномасштабной конвекции. По мере дальнейшего понижения эффективной температуры основным окном прозрачности становится область 1.7—2 мкм, переходная между связанно-свободными и свободносвободными переходами на ионах Н_.
Диапазон 3000—3600 А, где в фотосфере на пер-
Wavelength, m
Рис. 8. Изменение коэффициента непрерывного поглощения и рассеяния с длиной волны для разных температур 0 = 5040/T и газового давления log Pg = 3.0 из [60].
выи план выходит релеевское рассеяние на атомах и молекулах водорода, а в атмосфере излучение блокируется сотнями тысяч абсорбционных линий, для холодной звезды является второстепенным в смысле переноса излучения. Проявления фотометрической и спектральной переменности в этом диапазоне могут и не иметь прямой связи с механизмами физической переменности, воздействующими на фотосферу. Хорошим подтверждением является фотометрическая и спектральная переменность звезд типа Миры Кита, где уже более 30 лет известны изменения, происходящие на коротких временах. Здесь, скорее, наблюдаются процессы, происходящие в хромосферах и оболочках [63]. Интерпретация изменений показателя цвета (и — В) затруднена взаимодействием эмиссионных и абсорбционных деталей, а систематический спектральный мониторинг холодных переменных звезд в наземном ультрафиолете не проводился.
7. ВЫВОДЫ Практически все задачи, перечисленные выше, не могут быть решены без повышения точности спектроскопических наблюдений в наземном ультрафиолете. Кроме того, существует ряд задач, требующих систематического спектрального мониторинга в этом диапазоне. Спектры, полученные нами для различных объектов, показали пригодность спектрографа НЭС для большинства перечисленных задач.
БЛАГОДАРНОСТИ Работа поддержана РФФИ (проект 07—02— 00247—а) и Программой Отделения физических наук РАН.
СПИСОК ЛИТЕРАТУРЫ
1. В. Г. Клочкова, В. Е. Панчук и М. В. Юшкин, в Ультрафиолетовая Вселенная. II., ред. Б. М. Шустов, М. Е. Сачков, и Е. Ю. Кильпио (Янус-К, Москва, 2008), с. 46.
2. И. М. Копылов, Изв. КрАО 26,232(1961).
3. R. L. Kurucz, Astrophys. J. Suppl. 40, 1 (1979).
4. С. А. Гуляев, В. Е. Панчук и др., Астрофиз. Исслед. (Изв.САО) 22,3(1986).
5. С. А. Гуляев, О. Г. Патимова и В. В. Плешаков, Астрофиз. Исслед. (Изв.САО) 25, 55 (1987).
6. С. А. Гуляев, В. Е. Панчук и С. Г. Пяткес, Астрофиз. Исслед. (Изв.САО) 25, 60 (1987).
7. И. М. Копылов, Изв. КрАО 30, 69 (1963).
8. G. Shajn and O. Struve, Monthly Notices Roy. Astronom. Soc. 89, 222 (1929).
9. T. R. Stoeckley, Monthly Notices Roy. Astronom. Soc. 140, 121 (1968).
10. G.W. Collins II, Astrophys. J. 191,157(1974).
11. S. R. Heap, Be and Shell Stars, IAU Symp No. 70, Ed. byA.E. Sletteback, p. 165(1976).
12. J. B. Hutchings, Publ. Astronom. Soc. Pacific 88, 5
(1976).
13. J. B. Hutchings and T. R. Stoeckley, Publ. Astronom. Soc. Pacific 89, 19(1977).
14. J. B. Hutchings, J. M. Nemec, and J. Cassidy, Publ. Astronom. Soc. Pacific 89, 19 (1977).
15. M. Ruusalepp, A. Sapar and L. Sapar, Tartu publ. 50, 152(1984).
16. G. H. Sonneborn and G.W. Collins II, Astrophys. J. 213,787(1977).
17. A. Slettebak, T. J. Kuzma, and G.W. Collins II, Astrophys. J. 242, 171 (1980).
18. M. Ruusalepp, Tartu publ. 48, 81 (1981).
19. M. Ruusalepp, Tartu publ. 51, 84 (1986).
20. S. R. Heap, Astrophys. J. 217, 90 (1977).
21. K. G. Carpenter, A. Slettebak, and G. Sonneborn, Astrophys. J. 28б, 741 (1984).
22. R. L. Kurucz, Astrophys. J. 1SS, L21 (1974).
23. Y. Takeda, S. Kawanomoto, and N. Ohishi, Astrophys. J. б78, 446 (2QQ8).
24. J. M. Marlborough, Be Stars, IAU Symp No. 98, Ed. by M. Jaschek and H.-G. Groth, p. 361 (1982).
25. N. R. Walborn, Astrophys. J. 205, 419 (1976).
26. C. D. McKeith, P. L. Dufton and L. Kane, Observatory 98,263(1978).
27. N. R. Walborn, Rev. Mex. Astron. Astrofis. 33, 5 (2QQ8).
28. В. Г. Клочкова и В. Е. Панчук, УФН 1б4,657 (1994).
29. H. Reeves, W. A. Fowler, and F. Hoyle, Nature 22б, 727 (197Q).
3Q. C. P. Deliyannis and M. H. Pinsonneaut, Astrophys. J. Зб5, L67 (199Q).
31. J. A. Thorburn, L. M. Hobbs et al., Astrophys. J. 425, 15Q (1993).
32. В. Г. Клочкова и В. Е. Панчук, Астрон. ж. 40, 829 (1996).
33. R. J. Garcia Lopez, R. Rebolo, and M. R. Perez de Taoro, Astronom. and Astrophys. 302, 184 (1995).
34. A. M. Boesgaard and M. J. Tripicco, Astrophys. J. 302, L49 (1986).
35. C. P. Deliyannis et al., Astrophys. J. 49S, L147 (1998).
36. A. M. Boesgaard, E. Armengaud, and J. R. King, Astrophys. J. 5S2, 41Q (2QQ3).
37. M. Spite and F. Spite, Astronom. and Astrophys. 115, 357 (1982).
38. R. Rebolo et al., Astronom. and Astrophys. 193, 193
(1988).
39. S. G. Ryan et al., Astrophys. J. 34S, L57 (199Q).
4Q. G. Gilmore et al., Nature 357, 379 (1992).
41. R. N. Boyd and T. Kajino, Astrophys. J. ЗЗб, L55
(1989).
42. D. Thomas et al., Astrophys. J. 40б, 569 (1993).
43. F. Primas et al., Astronom. and Astrophys. Зб2, 666 (2QQQ).
44. T. V. Mishenina et al., Astronom. and Astrophys. 353, 978 (2QQQ).
45. T. P. Walker etal., Astrophys. J. 413, 562 (1993).
46. D. K. Duncan, D. L. Lambert, and M. Lemke, Astrophys. J. 401,584(1992).
47. B. Edvardsson et al., Astronom. and Astrophys. 290, 176(1994).
48. V. G. Klochkova et al., Chin. J. Astron. Astrophys. б, 579 (2QQ6).
49. M. S. Bessel and J. Norris, Astronom. J. 2S5, 622 (1984).
5Q. P. E. Nissen et al., Astronom. and Astrophys. 2S5,44Q (1994).
51. A. M. Boesgaard et al., Astronom. J. 117,492(1999).
52. J. Norris, S. G. Ryan, and T. C. Beers, Astrophys. J. 4SS, 35Q (1997).
53. J. L. Crawford et al., Astronom. J. 11б, 2489 (1998).
54. H. R. Butcher, Nature 32S, 127(1987).
55. P. Francois, M. Spite, and F. Spite, Astronom. and Astrophys. 274,821 (1993).
56. C. Sneden et al., Astrophys. J. 4б7, 819 (1996).
57. J. J. Cowan et al., Astrophys. J. 521, 194 (1999).
58. В. Е. Панчук и др., Препринт № 169 (САО РАН, Нижний Архыз, 2QQ2).
59. C. Sneden and J. J. Cowan, Rev. Mex. Astron. Astrofis. 10,221 (2QQ1).
6Q. В. Е. Панчук, Астрометрия и астрофизика 25, 2Q (1975).
61. S.C. Wolff, Publ. Astronom. Soc. Pacific 90, 412 (1978).
62. И. И. Романюк и др., Астрофизический бюллетень б2, 32 (2QQ7).
63. К. В. Бычков и В. Е. Панчук, Астрон. ж. 54, 34Q
(1977).
SPECTROSCOPY OF STARS IN GROUND-BASED ULTRAVIOLET. II. ATMOSPHERES
V.E. Panchuk, V.G. Klochkova, M.V. Yushkin
Some of the tasks of spectroscopic studies of stellar atmospheres in the ground-based ultraviolet are reviewed in historical perspective. Examples used include fragments of spectra that we obtained with the 6-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences.
Key words: stars—variable and peculiar; methods of astronomical observations, equipment and instruments
