CC BY
33
УДК 524.31.01-355
ИССЛЕДОВАНИЕ НЕСТАЦИОНАРНОСТИ АТМОСФЕРЫ а CYG. IV: НЕКОТОРЫЕ ОСОБЕННОСТИ ПЕРЕМЕННОСТИ ПРОФИЛЕЙ ЛИНИЙ
СО ВРЕМЕНЕМ
©2010 А. Х. Рзаев*
Специальная астрофизическая обсерватория, Нижний Архыз, 369167 Россия Поступила в редакцию 6 мая 2009 г.; принята в печать 25 июня 2009 г.
По 240 ПЗС—спектрам, полученным на кудэ-эшелле-спектрографе 2-м телескопа ШАО НАН Азербайджана в 1998—2000 гг., было исследовано изменение со временем профилей и лучевых скоростей линий ионов, Не! и Н^ в спектре сверхгиганта а Cyg. Показано, что эти изменения обусловлены движениями типа пульсаций, существующих в атмосфере этой звезды. Для линий ионов и Не! колебания происходят в основной фундаментальной моде с периодом около 12.0 ± 0.5^ и амплитудой 5.0 ± 0.5 км/с. Их продолжительность для линий ионов составляет около 35 дней.
Далее, вследствие появления разницы между лучевыми скоростями сильных и слабых линий ионов постепенно происходит затухание этих колебаний в течение примерно 5.0 ± 1.0^, после чего процесс повторяется. Для линии Н^ были найдены два значимых периода и две амплитуды, а также три характерных поведения переменности лучевой скорости для синей и красной половин профиля абсорбции: параметры переменности (период Р, амплитуда А) одинаковые; Р и Л одинаковые, но наблюдается сдвиг фазы между колебаниями лучевой скорости синей и красной половин профиля абсорбции; Р и Л различается для двух половин профиля абсорбции. Лучевая скорость центра массы звезды, определенная по 7-скорости составляет —4.5 ± 0.5 км/с. Средняя скорость расширения слоев атмосферы, где формируется линия Н^, около —16.5 ± 0.5 км/с, она меняется со временем с амплитудой около 3.0 км/с.
Ключевые слова: звезды — переменные и пекулярные
Большие значения периода (Р > Р0 14.3й) при-
писывались нерадиальным пульсациям. Иноу [5] обнаружил стратификацию лучевой скорости в атмосфере а Cyg. Линии ионов, предположительно фотосферного происхождения, показывали амплитуду изменения лучевой скорости около 6 км/с. В более высоких слоях, где формируются бальмеров-ские линии с малыми номерами и линия MgII Л 4481 Л, амплитуда достигала 15 км/с. А в самых высоких слоях атмосферы, где формируются линии На, Н и К MgII, амплитуда возрастала от 50 до 150 км/с.
По эшелле-ПЗС спектрам (Н=20000) сверхгигант а Cyg исследован Гейдельбергской группой [6, 7]. Было показано, что переменность лучевой скорости и профилей линий ионов обусловлена движениями типа пульсаций. Одновременно возбуждается несколько мод пульсаций. Характерное время переменности лучевой скорости со временем различается для разных сезонов. Для линий ионов не обнаружено стратификации лучевой скорости. По линиям Н27—Н8 обнаружен
1. ВВЕДЕНИЕ
Сверхгигант a Cyg — одна из ярких звезд на северном небе. Она была исследована в широком интервале спектра от далекого УФ до радиодиапазона [1]. Абт [2] был первым, кто по длительным спектральным наблюдениям исследовал нестационарность сверхгигантов, в том числе и a Cyg. Он обнаружил, что у сверхгигантов наблюдаются изменения лучевой скорости, напоминающие пульсации. В случае a Cyg характерное время переменности лучевой скорости, найденное
по линиям Fell Л 4508, 4515, 4520 и 4522 A, оказалось около 13.8d. Ссылаясь на работу Паддока [3] он принял эту величину равной 11.7d. Люси [4] по 447 измерениям лучевой скорости линий ионов Паддока [3] пришел к заключению, что переменность обусловлена одновременным возбуждением нескольких дискретных пульсаций. Люси нашел 16 значимых периодов в интервале от 6.9 до 100.8d.
E-mail: abid@sao.ru
небольшой Бальмеровский прогресс. Переменность профиля На происходит из-за наложения переменных сине- и красносмещенных эмиссий на стабильный “фотосферный” профиль. Переменность проявляется одновременно и симметрично относительно скорости центра массы звезды. Отмечено, что определить характерное время этих изменений весьма трудно. По изменениям эквивалентной ширины линии На было найдено, что оно близко к величине периода вращения звезды. Авторы сделали заключение, что переменность профиля линии На обусловлена модуляцией вращения.
Наши недавние исследования показали, что традиционный метод изучения изменений со временем лучевой скорости (измеряемой по бисектору абсорбции) приемлем для линий ионов, когда наблюдается радиальная пульсация. Если разделить линии по бисектору на две части и исследовать переменность лучевой скорости отдельно по синей и красной половин профиля линии, то в случае линий Н@ и На картины и параметры переменности лучевой скорости для синей и красной половин профиля могут различаться. Для изучения стратификации лучевой скорости в атмосфере звезды необходимо исследовать переменность лучевой скорости отдельно для синей и красной половин профилей линий на разных уровнях их остаточной интенсивности [8, 9] (далее, соответственно, статья II и III). Параметры переменности (амплитуда и характерное время) также могут различаться для разных уровней остаточной интенсивности г (статья III).
Нами было получено, что для линий ионов изменения лучевой скорости обеих половин профиля абсорбции на всех уровнях ее интенсивности происходят одинаково. Как и в работе Кауфер и др. [7], в наших наблюдениях для конкретного сезона наблюдений для всех линий ионов параметры пульсаций (период, амплитуда) одинаковые. Но для разных сезонов наблюдений периоды различаются.
В случае линии Нв картины и параметры переменности лучевой скорости для синей и красной половин профиля абсорбции различаются. Для красной половины профиля параметры переменности одинаковы для всех уровней остаточной интенсивности линии, а для синей половины отличаются для разных уровней (статья II).
Наши исследования показали, что в случае линии На картина переменности лучевой скорости также обусловлена нерадиальными пульсациями (статья III). Параметры пульсации различаются для синей и красной половин профиля абсорбции. Для красной половины профиля амплитуда и период увеличиваются от ядра в сторону крыла линии. На синей половине профиля дополнительно наблюдаются абсорбционные детали. Их появление
и исчезновение, а также незначительные миграции свидетельствуют о клочковатой структуре оболочки звезды. Прослеживается движение отдельных сгустков в верхних слоях атмосферы звезды, которые проявляются в виде абсорбционных деталей.
Наши исследования также показали, что в спектре a Cyg в течение 5—6 ночей (с 06.09.99 г. по 11.09.99 г.) заметных изменений фотометрических и позиционных параметров линий не обнаруживается, т. е. звезда находилась в состоянии “покоя” (статья II).
В настоящей статье по вышеуказанной методике исследована переменность профилей линий ионов, HeI и Нв в спектре a Cyg. В дополнение к 120 спектрам, полученным в 1998—99 гг. (статья II), в исследование были включены еще 120 эшелле— спектров, полученных в 2000 г., что позволило более подробно исследовать переменность этих линий со временем.
2. НАБЛЮДЕНИЯ И ОБРАБОТКА
Спектры были получены в кудэ спектрографе 2м телескопа ШАО НАН Азербайджана на эшелле-спектрографе CEGS с матрицей 580 х 530 элементов [10]. Размер матрицы ставил ограничение как на число порядков, так и на длину одного порядка. Наклон решетки выбирался таким образом, что в кадр попадал спектральный диапазон ЛЛ 4340— 7000A и одной экспозицией можно было фиксировать линии Ha и Нв [11].
В 1998 г. в течение двух, в 1999 г. в течение полутора и в 2000 г. в течение четырех месяцев были получены по 40, 80 и 120 ПЗС-эшелле-спектров с R = 36000 и 250 < S/N < 500. Во всех трех сезонах каждую ночь получались как минимум по
4 спектра (статья III). Изменения лучевой скорости и профиля линий в течение одной ночи не были обнаружены, поэтому профили усреднялись. Общее количество наблюдательных ночей в течение 3-х лет составило 52, соответственно, по 11, 16 и 25 ночей (точек) в 1998, 1999 и 2000 гг.
Обработка эшелле изображений выполнена стандартной методикой с помощью пакета программ MIDAS. Измерение фотометрических и позиционных параметров линий осуществлено пакетом программ DECH20 [12]. Как и в предыдущих работах, для уменьшения шумов и увеличения точности фотометрических и позиционных измерений мы применили процедуру сглаживания спектра методом гауссианы (статья II). Измерение лучевой скорости производилось следующим образом. Профили исследуемых линий записывались в r — V координатах. Далее отдельно для синей и красной половин профиля на разных уровнях остаточной интенсивности r с шагом ±0.05 (например,
для г = 0.45,0.5 и т.д.) брались отчеты лучевой скорости. Для каждого уровня г (например, для г = 0.55) данные о лучевой скорости усреднялись по 5 соседним уровням с шагом ±0.01. Такой подход значительно уменьшил ошибки измерений, дал возможность детально исследовать картину переменности лучевой скорости со временем. Для сильных линий ионов, линий Нв и для слабых линий ионов и Не! среднеквадратические отклонения лучевой скорости от среднего, соответственно, не превышали ±0.3, ±0.5 и ±1.1 км/с.
3. ПЕРЕМЕННОСТЬ ЛИНИИ
3.1. Переменность лучевой скорости
3.1.1. Линии ионов
Как было отмечено выше, в случае линий ионов переменность лучевой скорости со временем, исследованная по бисектору абсорбции на всех уровнях интенсивности г, реально отражает переменность лучевой скорости по линии в целом (статья II). Данные по лучевым скоростям линий ионов, полученные за 2000 г. также подтвердили это. Поиск возможных периодов изменения лучевой скорости со временем проводился методом гармонического анализа, который был подробно описан в статье III.
Как было отмечено в статье II, картины переменности отдельных профилей относительно среднего профиля для сезонов 1998 и 1999 гг. были одинаковые. Но периоды переменности лучевой скорости различались и равнялись 20 и 12.5й, соответственно, для 1998 и 1999 гг. Для профилей линий ионов, полученных в 2000 г., картины переменности профилей относительно среднего профиля оказались такими же, что для предыдущих годов. А по измерениям лучевой скорости, полученным в 2000 г., мы никаких значимых периодов не нашли. Были получены периоды около 14 и 30 дней с низкой значимостью. Отметим, что в 2000 г. (как и в 1999 г.) в течение 5 ночей (JD—1740-1745) в пределах ошибок измерений также не обнаружено заметных изменений фотометрических и позиционных параметров линий. Наши исследования показали, что обнаружение разных периодов для разных сезонов наблюдений связано с этим фактом. При поиске периодов необходимо исключать постоянные значения лучевой скорости.
Переменность лучевой скорости со временем для разных сезонов наблюдений представлена на Рис. 1. На Рис. 1с, чтобы не менять масштаб по оси абсцисс, данные за JD 1808, 1810, 1866 и 1868 не нанесены. Поиск возможных периодов по данным Рис. 1 выполнялся по следующей схеме. Мы исключили постоянные значения лучевой скорости, т.е. данные JD: 1428—31 и 1740—45 за 1999
- 1999
■ А У (К
1\ У
■т ’ V V і;
1420 1430 1440 1450 1460 1470
Л> 2450000+
Рис. 1. Переменность лучевой скорости со временем для линий ионов в разные годы наблюдений. Вертикальные размеры символов соответствуют ошибкам измерений лучевой скорости. Кривая (синусоида) проведена методом гармонического анализа. По оси абсцисс указаны юлианские даты наблюдений.
и 2000 гг. Далее поиск проводился как по данным отдельных сезонов, так и по всем данным за 3 года с 1998 по 2000 гг. Во всех случаях параметры переменности лучевой скорости (период, амплитуда и 7-скорость) были одинаковыми. Но данные плохо ложились на общую фазовую кривую. Сопоставление фазовых кривых (построенных как для отдельных сезонов, так и для всех данных за 3 года)
Phase
Рис. 2. Фазовая кривая переменности лучевой скорости для линий ионов. Квадратики, треугольники и темные кружки, соответственно, данные за 1998, 1999 и 2000 гг. Начальная фаза — JD 2451746. Кривые на Рис. 1 за 1998 и 1999 гг. сдвинуты по фазе. В правом нижнем углу вертикальный отрезок указывает ошибку измерений лучевой скорости.
показало, что необходимо учитывать сдвиг фазы (возникающий за счет того, что в отдельные периоды наблюдения наблюдаются постоянные значения лучевой скорости) как для отдельных сезонов, так и для отдельных сетов наблюдений. На Рис. 1 для каждого сезона наблюдений приведены кривые переменности лучевой скорости со временем. Кривые построены для параметров переменности лучевой скорости: P = 12.5 ± 1.5d; A = 5.0 ± 0.5 км/с; Y = —4.5 ± 0.5 км/с, найденных методом гармонического анализа. Фазовая кривая для данных, полученных за 3 года, с учетом сдвига фаз для разных сезонов, приведена на Рис. 2.
3.1.2. Линии HeI
В спектре a Cyg линии HeI очень слабые, поэтому они до сих пор не были исследованы. В среднем центральные интенсивности r0 линий HeI Л 4471,
5876 и 6678 A составляют, соответственно, 0.89, 0.88 и 0.95. На спектрах, полученных в 1999 и 2000 гг., отношение S/N было более высоким, чем в
1998 г. и переменность лучевой скорости линий HeI со временем исследовалась для этих двух сезонов по вышеуказанным трем линиям. Как и в случае линий ионов, лучевые скорости измерялись по бисектору абсорбции на уровне половины интенсивности линий. Картины переменности лучевой скорости со временем для линий HeI оказались такими же, как для линий ионов того же сезона. Параметры переменности лучевой скорости для линий HeI, с учетом ошибок измерений, оказались одинаковыми с параметрами для линий ионов: P = 12.5 ± 2.0d; A = 6.0 ± 1.5 км/с; y = —5.5 ± 1.5 км/с. Фазовая
Рис. 3. То же, что на Рис. 2, но для линий НеЕ Квадратики и темные кружки, соответственно, данные за 1999 и 2000 гг.
кривая по данным двух сезонов с учетом сдвига фаз между этими сезонами приведена на Рис. 3.
3.1.3. Линия Нв
Как и в предыдущей работе (статья II), измерение лучевой скорости выполнялось отдельно для синей и красной половин профиля абсорбции Нв на разных уровнях остаточной интенсивности. Центральная остаточная интенсивность линии Нв в среднем составляет г0 = 0.35 ± 0.015. Наши исследования показали, что в области ядра (где т < 0.4) и в крыльях линии (т > 0.7) точность измерений ухудшается, поэтому измерения проводились в интервале
0.4 < т < 0.7 с шагом Ат = 0.05 (т.е. для семи уровней т = 0.4,0.45,0.5,0.55, 0.6,0.65,0.7). Для каждого уровня измерения лучевой скорости выполнялись вышеописанным методом (см. раздел 2). Такой подход показал, что для обеих половин профиля абсорбции на всех вышеуказанных уровнях т параметры переменности лучевой скорости со временем одинаковые. Поэтому кривые на рисунках приведены для среднего уровня т = 0.55 (т.е. средняя лучевая скорость по семи уровням). Ошибки нахождения амплитуды и 7-скорости также были определены как отклонение от среднего значения этих величин, найденных по семи уровням. На Рис. 4—7 вертикальным отрезком отмечены ошибки нахождения амплитуды, а ошибки измерения лучевой скорости приведены в разделе 2.
Поиск возможных периодов переменности лучевой скорости со временем также проводился вышеописанным способом. По измерениям лучевой скорости, полученным в 1998 г., для обеих половин профиля абсорбции в пределах ошибок были найдены одинаковые значения периода и амплитуды: Рь = 13.5 ± 2.5й, Аь = 7.0 ± 1.0 км/с
Phase
Рис. 4. Фазовая кривая переменности лучевой скорости для линии H,3 в 1998 г. (а) для синей;
(b) для красной половины профиля с параметрами P = 13.5 ± 2.5d, A = 7.0 ± 1.0 км/с; (с) для бисектора абсорбции с параметрами P = 13.5 ± 2.5d,
A = 2.3 ± 0.6 км/с. Начальная фаза JD 2451053. Для наглядного представления сдвига фаз на величину порядка 0.6 между переменностью лучевой скорости для синей и красной половин профиля абсорбции лучевые скорости за JD 2451053—1054 обозначены светлыми квадратиками. В правом нижнем углу вертикальный отрезок указывает ошибку определения амплитуды, а ошибки измерения лучевой скорости приведены в разделе 2.
и Pr = 13.5 ± 2.5d, Ar = 7.5 ± 0.8 км/с, соответственно, для синей и красной половин профиля абсорбции. Фазовые кривые представлены на Рис. 4а, b. Но наблюдается сдвиг фазы (на при-
Рис. 5. Фазовая кривая переменности лучевой скорости для бисектора абсорбции линии Н,3 в 1999 г. с параметрами Р = 20.5 ± 2.5^, А =10.0 ± 1.0 км/с. В правом нижнем углу вертикальный отрезок указывает ошибку определения амплитуды.
мерно 0.6Р) между колебаниями лучевой скорости для синей и красной половин профиля абсорбции. Это приводит к тому, что амплитуда колебания лучевой скорости, измеренной по бисектору абсорбции, значительно уменьшается. Таким образом, для параметров колебания, найденным по лучевой скорости бисектора абсорбции были получены следующие значения: Ры3 = 13.5 ± 2.5й, Аы3 = 2.3 ± 0.6 км/с и 7 = —15.5 ± 0.5 км/с. Фазовая кривая представлена на Рис. 4с.
По данным, полученным в 1999 г., для обеих половин профиля абсорбции на всех уровнях т значения периода, амплитуды и фаза колебания переменности лучевой скорости в пределах ошибок измерения оказались одинаковыми. Следовательно, переменность лучевой скорости бисектора абсорбции на всех уровнях т реально отражает переменность лучевой скорости линии в целом. Поэтому на Рис. 5 приведена фазовая кривая переменности лучевой скорости для бисектора абсорбции на уровне т = 0.55. Параметры колебания следующие: Ры3 = 20.5 ± 2.5й, АЫз = 10.0 ± 1.0 км/с и 7 = —12.0 ± 0.5 км/с.
В 2000 г. картина переменности лучевой скорости со временем различается для синей и красной половин профиля абсорбции. Для уровня т = 0.55 эти кривые приведены на Рис. 6. Как видно из Рис. 6, для обеих половин профиля абсорбции в течение 5 ночей (JD 1743—1747) в пределах ошибок измерений переменность лучевой скорости не наблюдалась. Поэтому, как и в случае линии ионов, при поиске возможных периодов эти значения лучевой скорости не были использованы. На Рис. 6 синусоиды были приведены на основе параметров переменности лучевой
скорости со временем, найденных гармоническим анализом. Для синей и красной половин профиля абсорбции, соответственно, были найдены следующие параметры переменности лучевой скорости: Рь = 14.0 ± 1.5й, Аь = 7.5 ± 0.5 км/с, 7 = —57.0 ± 0.5 км/с и Рг = 22.0 ± 2.0й, Аг = 11.5 ± 0.5 км/с, 7 = 28.0 ± 0.5 км/с. Фазовые кривые представлены на Рис. 7.
3.2. Переменность фотометрических параметров линий
В работе [7] линии ионов по их эквивалентной ширине были разделены на 3 группы: слабые 50 < < 200, средние 200 < < 500 и сильные > 500 тА. Расчет модели атмосферы [7]
показал, что интенсивные линии формируются в верхних слоях атмосферы, в отличие от слабых линий ионов. Между указанными группами линий ионов разницы в лучевой скорости не обнаружено. Исследование переменности интенсивности линий ионов со временем показало, что в спектре а С;^ среднеквадратическое и максимальное отклонение от среднего не превышает 3 и 16 % соответственно.
Для повышения точности измерения эквивалентных ширин и интенсивностей линий мы поступили следующим образом. Участки спектра, где были проведены измерения, были сопоставлены друг с другом и путем (медианного) усреднения были выбраны опорные точки (по координате длины волны) для проведения непрерывного спектра. При необходимости для линий с широкими крыльями (например, для Нв) был использован непрерывный спектр соседних порядков. Для линий ионов слабой и средней интенсивности были отобраны более “чистые” (неблендированные) линии. В итоге было отобрано для изучения 10, 10 и 5 линий слабой, средней и сильной интенсивности соответственно.
Наши исследования показали, что переменность со временем эквивалентных ширин и центральной интенсивности линий ионов для сильной группы составляет около 3.0%, средней < 3.0% и <1.5%, а для слабой <5.0% и < 1%. Для линии Нв средние значения эквивалентной ширины и интенсивности, найденные по всему материалу, составляют = 2.500 ± 0.046 А (переменность < 2%),
то = 0.35 ± 0.01 (переменность <3%), максимальные отклонения от среднего не превышают 8%.
4. ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ И ЗАКЛЮЧЕНИЕ
Итак, мы показали, что в спектре а С;^ в пределах ошибок измерений переменность со временем эквивалентных ширин и интенсивности линий ионов и Нв не наблюдается. Как было сказано
1720 1740 1760 1780
Л) 2450000+
Рис. 6. Переменность лучевой скорости со временем для синей (а) и красной (Ь) половин профиля абсорбции линии Н,3 за 2000 г. Кривые (синусоиды) проведены на основе параметров переменности лучевой скорости, найденных методом гармонического анализа. По оси абсцисс указаны юлианские даты наблюдений. В правом нижнем углу вертикальный отрезок указывает ошибку определения амплитуды.
в работах [7, 13], это свидетельствует о том, что наблюдаемая переменность лучевой скорости обусловлена движениями типа пульсаций с малыми амплитудами.
Наш методический подход, описанный в разделе 2, позволил выявить картину переменности лучевой скорости со временем. По этим картинам (Рис. 1 и 6) было обнаружено, что иногда в атмосфере а С;^ переменность лучевой скорости не наблюдается в пределах ошибок измерений. Эти постоянные значения лучевой скорости наблюдаются в период затухания предыдущего и возникновения следующего колебания. Наши исследования показали, что обнаружение разных периодов для разных сезонов наблюдений связано с этим фактом. При поиске периодов необходимо исключить постоянные значения лучевой скорости.
Было обнаружено, что для всех исследованных линий эти постоянные значения в пределах ошибок измерений равны 7-скорости изменений лучевой скорости со временем. Для линий ионов 7-скорость, по данным трех сезонов
(1998—2000 гг.) равна —4.5 ± 0.5 км/с, что соответствует лучевой скорости центра массы звезды в пространстве. Наши исследования показали, что асимметрия профиля линии Нв, которая в основном вызвана изменением синего крыла и скорость расширения соответствующих слоев различается в разные сезоны наблюдений. Для уровня т & 0.55, представленного на Рис. 4—7 7-скорость меняется в пределах от —12.0 до —15.0 км/с. На уровне FWHM (т ~ 0.675) она составляет — 16.0... —19.0 км/с. Например, на Рис. 6 для синей и красной половин профиля абсорбции постоянное значение лучевой скорости, соответственно, около —57 и 27 км/с, следовательно, для бисектора абсорбции оно около —15 км/с.
Таким образом, постоянная лучевая скорость приводит к сдвигу фазы как между колебаниями разных сезонов, так и между колебаниями отдельных сетов наблюдений данного сезона. Для всех исследованных линий продолжительность времени, когда наблюдаются постоянные значения лучевой скорости, по нашим данным, составляет
5.0 ± 1.0й.
На Рис. 1с видно, что с ^ 1740 по 1745 и, по-видимому, с ^ 1779 по 1784 наблюдаются постоянные значения лучевой скорости около —4.5 км/с. Следовательно, по данным 2000 г., продолжительность колебания (^ ~ 1745—1779) около 34 дней. По данным 1999 г. продолжительность колебания (^ ~ 1432—1467) около 35 дней (Рис. 1Ь). Если предположить, что в 1998 г., как и в других двух сезонах (1999—2000 гг.), колебание начинается и затухает при постоянных значениях лучевой скорости около —4.5 км/с, то продолжительность колебания будет (^ ~ 1050—1083) около 33 дней (Рис. 1а). Таким образом, продолжительность колебания, найденная по линиям ионов, составляет около 34.0 ± 1.0 дней.
Для периодов пульсаций в радиальной фундаментальной моде в работах [4] и [7] для а С;^ были получены величины Р0 & 14.0 и Р0 & 12.0й. В работах [1, 14] с учетом структуры расширяющейся атмосферы и звездного ветра была детально проанализирована атмосфера сверхгиганта а С;^. Путем сопоставления расчетного спектра с наблюдаемым были уточнены основные атмосферные и фундаментальные звездные параметры, в том числе МЪо1 и Те^. Но, если по последним двум параметрам (например, из работы [14], где совпадение расчетного и наблюдаемого спектра для фотосферных и ветровых линий очень хорошие, и МЪо1 = —8.49, а Teff = 8525К) по формуле из работы Лови и др. [15] определять величину Р0, она получится около 30 дней. Возможно это свидетельствует о том, что для уточненных параметров,
Phase
Рис. 7. Фазовая кривая переменности лучевой скорости для (а) синей и (Ь) красной половин профиля абсорбции линии Нв в 2000 г. с параметрами указанными в тексте. В правом нижнем углу вертикальный отрезок указывает ошибку определения амплитуды.
определенных современными методами, уравнение из [15] должно корректироваться.
Если выбирать величины Мъог и Те^ из литературы (ссылки в работах [1, 6, 14]), то по Лови и др. [15] значение Р0 будет в пределах от 11 до 15 дней. Поэтому будем считать, что найденные нами периоды 12.5 и 14.0 дней, в пределах ошибок,— это периоды радиальной пульсации в основной фундаментальной моде.
Наши исследования показали, что колебание слоев, где формируются линии ионов, начинается с расширения этих слоев (Рис. 1Ь, с). В начальной фазе расширения (^ 1433 на Рис. 1Ь и ^ 1746 на Рис. 1с) наблюдается незначительная разница между лучевыми скоростями слабых и сильных линий ионов. Например, в датах ^ 1433 и 1476 лучевая скорость по слабым и сильным линиям ионов, соответственно, составляла —8.0 ± 1.0, —6.2 ± 0.5 км/с и —6.5 ± 1.0, —4.8 ± 0.3 км/с. Это свидетельствует о том, что процессы (движения), которые приводят к возникновению пульсаций, происходят в нижележащих слоях фотосферы. Далее, например, начиная с ^ 1747 (Рис. 1с), лучевые скорости для всех групп линий ионов ста-
Рис. 8. Переменность лучевой скорости для слабой (темные квадратики), сильной (светлые квадратики) линий ионов в 1999 (а), и 2000 гг. (Ь), когда происходит затухание колебания. Темные треугольники — средняя скорость, найденная по всем линиям, светлые треугольники — лучевая скорость центра массы звезды. Вертикальные размеры символов соответствуют ошибкам измерения.
новились одинаковыми и показывали одинаковый ход и темп движения. Следующая разница в лучевых скоростях между слабыми и сильными группами линий была нами зафиксирована примерно через 31—32 дня, ^ 1465—1468 на Рис. 1Ь и ^ 1777—1779 на Рис. 1с. Переменность лучевой скорости в этих датах отдельно для сильных и слабых групп, а также для среднего значения по всем трем группам линий представлена на Рис. 8. Видно, что в течение указанных ночей разница лучевых скоростей между сильной и слабой группами линий уменьшается, а средняя лучевая скорость постепенно становится равной 7-скорости линий ионов. Таким образом, постепенно происходит затухание колебания слоев, где формируются линии ионов.
Для линий Неї параметры переменности лучевой скорости со временем с меньшей точностью измерения совпадают с таковым для линий ионов: Р = 12.5 ± 2.0й; А = 6.0 ± 1.5 км/с; 7-скорость: —5.5 ± 1.5 км/с. С учетом сдвига фазы данные за
1999 и 2000 гг. хорошо ложатся на фазовую кривую на Рис. 3. В атмосфере сверхгиганта линии Неї преимущественно фотосферного происхождения, что еще раз свидетельствует о достоверности полученных результатов для линии ионов.
В отличие от линий Неї и ионов, для линии Нв характер и параметры переменности лучевой скорости различаются для разных сезонов наблюдений. В общей сложности мы получили два периода и амплитуды: Р & 14й с А & 7 км/с и Р & 21й с А & 11 км/с, а также три характерных поведения переменности лучевой скорости для синей и
красной половин профиля абсорбции. В 1998 г. параметры переменности, т.е. период и амплитуда для обеих половин профиля абсорбции, одинаковые, но фазы колебаний различаются. В 1999 г. как параметры переменности, так и фаза колебания одинаковы для обеих половин профиля абсорбции. В 2000 г. параметры переменности различаются для синей и красной половин профиля абсорбции. Минимальная продолжительность этих колебаний по нашим данным около 40—45 дней.
Отметим, что наши исследования как для линий ионов, так и для линии Нв подтвердили результат, полученный в работе [7]. Возбуждение новых колебаний происходит после исчезновения (затухания) предыдущих. Полученные три характерные переменности лучевой скорости для синей и красной половин профиля абсорбции линии Нв позволяют нам судить об особенности колебания этих слоев атмосферы. Вначале возбуждается колебание (с определенными параметрами колебания) в поднимающихся слоях атмосферы, которое передается (отражается) в опускающиеся слои атмосферы. Через некоторое время (по нашим данным через 40—45 дней) происходит сдвиг фазы между колебаниями поднимающихся и опускающихся слоев атмосферы, который постепенно может привести к затуханию колебания. При этом диссипация энергии может являться источником возбуждения очередного колебания как в нижних (где формируется линия Нв), так и для вышележащих слоев атмосферы.
По нашим исследованиям, средняя скорость расширения слоев атмосферы, где формируется линия Нв, около —13.0 км/с и она меняется со временем с амплитудой около 3.0 км/с. Суперпозиция этого колебания на колебание с Р = 14.0й и А = 7.0 км/с вполне может привести к колебанию с большей амплитудой и с большим периодом (Р = 22.0й, А = 11.0 км/с).
В спектре а С;^ различие параметров переменности для синей и красной половин профиля абсорбции наблюдается и для линии На (статья III). Наши исследования также показали, что в спектре звезды HD 93521 (09.5 V) для всех исследованных линий — На, Нв и Не1, которые формируются в слоях атмосферы, где истекает звездный ветер, параметры переменности различаются для двух половин профиля абсорбции [16]. Все это свидетельствует о том, что, во-первых, у ветровых линий синяя и красная половины профиля абсорбции формируются в разных областях атмосферы. Во-вторых, влияние переменности звездного ветра различается для этих двух половин.
Отметим, что в 2000 г. по лучевым скоростям, измеренным для бисектора абсорбции линии Нв
были найдены два значимых периода со значениями около 14.0 и 22.0 дня. Обычно в таких случаях принято говорить, что одновременно возбуждаются две моды колебания. Но в нашем случае это означает, что параметры переменности лучевой скорости различаются для двух половин профиля абсорбции.
Наш методический подход позволил выявить новые интересные особенности в картине переменности линий, подтверждение и объяснение которых требует дальнейших исследований, а также модельных расчетов.
БЛАГОДАРНОСТИ
Автор благодарит В. Г. Клочкову, Е.Л.Ченцова и В. Е. Панчука за обсуждение результатов, В. Д. Бычкова за предоставление программ гармонического анализа. Работа частично поддержана грантом РФФИ (проекты № 08-02-00072).
СПИСОК ЛИТЕРАТУРЫ
1. J. P. Aufdenberg, P. H. Hauschildt, E. Baron, et al., Astrophys. J. 570, 344 (2002).
2. H. A. Abt, Astrophys. J. 126,138(1957).
3. G. F. Paddock, Lick. Obs. Bull. 17, 39(1935).
4. L. B. Lucy, Astrophys. J. 206, 499 (1976).
5. M. O. Inoue, Publ. Astronom. Soc. Japan 31, 11 (1979).
6. A. Kaufer, O. Sthal, B. Wolf, et al., Astronom. and Astrophys. 305,887(1996).
7. A. Kaufer, O. Sthal, B. Wolf, et al., Astronom. and Astrophys. 320, 273 (1997).
8. А. Х. Рзаев и др., Астрофизический бюллетень 62, 61 (2007).
9. А. Х. Рзаев, Астрофизический бюллетень 63, 25 (2008).
10. Ф. А. Мусаев, Письма в АЖ 19,776(1993).
11. А. Х. Рзаев, Н. О. Гасанов, Х. М. Микаилов и др., Циркуляр ШАО 95, 3 (1999).
12. Г. А. Галазутдинов, Препринт САО РАН №92, (Н. Архыз, 1992).
13. W. Unno, Y. Osaki, H. Ando, and H. Shibahashi, Nonradial Oscillations of stars, (University of Tokyo Press, 1979), p.30ff.
14. F. Schiller and N. Przybilla, Astronom. and Astrophys. 479, 849 (2008).
15. D. Lovy, A. Maeder, A. Noels and Gabriel M., Astronom. and Astrophys. 133, 307 (1984).
16. А. Х. Рзаев, Астрофизический бюллетень 63, 282 (2008).
NONSTATIONARITY OF a CYG ATMOSPHERE. IV: SOME PECULIARITIES OF LINE
PROFILES TIME VARIATION
A. Kh. Rzaev
We use 240 CCD spectra taken in 1998—2000 with the coude echelle spectrograph of the 2-m telescope of the National Academy of Sciences of Azerbaijan to study temporal radial velocity and line profile variations of the ion, Hel, and lines in the spectrum of the a Cyg supergiant. We demonstrate that these variations are caused by pulsation-type motions in the star’s atmosphere. Ion and Hel lines oscillate in the main fundamental mode with a period of about 12.0 ± 0.5d and an amplitude of 5.0 ± 0.5 km/s. These ion-line oscillations continue for about 35 days. Then the difference between the radial velocities of strong and weak ion lines results in a gradual decay of oscillations over a time interval of about 5.0 ± 1.0d. Thereafter the process repeats itself. For the H^ line we found two significant periods, two amplitudes, and three characteristic radial velocity variability behaviors for the blue and red halves of the absorption profile: with equal variability parameters (period P and amplitude A); with equal P and A, but with a phase shift between the radial velocity variations of the blue and red halves of the absorption profile; with different P and A for the two halves of the absorption profile. The star’s center of mass radial velocity as inferred from the 7-velocity is —4.5 ± 0.5 km/s. The average expansion velocity of the atmospheric layers, where the H^ line forms, amounts to about —16.5 ± 0.5 km/s and varies temporally with an amplitude of about 3.0 km/s.
Keywords: stars: variable and pequliar
