CC BY
19
АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ, 2018, том 73, № 4, с. 498-504
УДК 524.31.01:520.84
СВЕРХБЫСТРАЯ ПЕРЕМЕННОСТЬ ПРОФИЛЕЙ ЛИНИЙ В СПЕКТРАХ OBA-ЗВЕЗД: B1-ЗВЕЗДА р Leo
©2018 А. Ф. Холтыгин1,2*, А. А. Батраков1, С. Н. Фабрика3, А. Ф. Валеев3, И. М. Туманова1, О. А. Циопа4
1Санкт-Петербургский государственный университет, Санкт-Петербург, 198504 Россия 2Институт астрономии РАН, Москва, 119017 Россия 3Специальная астрофизическая обсерватория РАН, Нижний Архыз, 369167 Россия 4Главная (Пулковская) астрономическая обсерватория РАН, Санкт-Петербург, 196140 Россия Поступила в редакцию 26 апреля 2018 года; принята в печать 3 августа 2018 года
Настоящая работа является продолжением исследований сверхбыстрой переменности профилей линий в спектрах звезд ранних спектральных классов. Мы исследовали переменность профилей линий в спектре OBA-звезд с многорежимным фокальным редуктором SCORPIO на 6-м телескопе BTA. Обнаружены короткопериодические регулярные вариации линий H и He в спектрах звезды B11 звезды р Leo с периодами от 2 до 90 минут и нерегулярные вариации профилей линий на промежутке времени менее 1 минуты. Обсуждается возможное происхождение быстрых спектральных вариаций.
Ключевые слова: звезды: ранние типы — звезды: индивидуальные: р Leo
1. ВВЕДЕНИЕ
Профили линий в спектрах OBA-звезд переменны на временных шкалах от дней до часов [1—3]. Вариации профилей на минутных шкалах до последнего времени не были исследованы. В 2014 г. Хубриг и др. [4] обнаружили изменения профилей линий Si II и Fe II в спектрах A0-сверхгиганта HD 92207 на промежутках 1—2 минуты.
Эта работа стала триггером наших исследований сверхбыстрой переменности на минутных шкалах в спектрах звезд ранних классов. Чтобы проверить, широко ли распространены короткопериодические спектральные вариации среди OBA-звезд, мы проанализировали вариации профилей линий в спектрах ряда OBA-звезд, полученных с высоким (минута и меньше) временным разрешением при использовании редуктора светосилы 6-м телескопа БТА SCORPIO [5]. В работе [6] представлены результаты анализа сверхбыстрой переменности в спектре звезды HD 93521 (O9.5III) по наблюдениям на БТА в 2015 г. Были обнаружены регулярные вариации с периодами 4—5 и 32—36 минут.
В настоящей работе представлены результаты поиска быстрых вариаций в спектрах медленно вращающегося [7] сверхгиганта р Leo (BIa). Звезда
E-mail: afkholtygin@gmail.com
р Leo уже наблюдалась нами на БТА в январе-феврале 2004 г. и на 1.8-м телескопе Бохинсан-ской оптической астрономической обсерватории (Южная Корея) [8]. Получено 47 спектров звезды с высоким временным разрешением (4-10 мин), отношением S/N от 300 до 700 и спектральным разрешением от 45 000 до 60 000. Обнаружена переменность профилей линий H I, He I, Si II, Si III и N II. Найдены регулярные компоненты переменности профилей линий с периодами от 3h8 до 1d8. Сделан вывод о том, что переменность профилей связана с их вращательной модуляцией и нерадиальными фотосферными пульсациями (Non-Radial Pulsations, NRP).
Звезда р Leo включена в программу исследования магнитных полей ярких звезд северного неба [9, 10]. В 2005 г. были выполнены наблюдения р Leo на БТА с использованием спектрографов NES и ОЗСП с анализатором круговой поляризации [8]. Было обнаружено умеренное магнитное поле, которое может быть описано в модели вращающегося магнитного диполя с напряженность на полюсах 240 ± 50 Гс и углом наклона оси диполя к оси вращения 59° ± 30°.
В настоящей статье анализируются наблюдения р Leo на БТА с фокальным редуктором SCORPIO со сверхвысоким временным разрешением 6 секунд.
В разделе 2 даны основные сведения о звезде и описаны наблюдения и обработка спектров. Вариации профилей линий проанализированы в разделе S. В разделе 4 обсуждаются результаты анализа переменности профилей. Выводы даны в разделе 5.
2. ОСНОВНЫЕ СВЕДЕНИЯ ОБ ОБЪЕКТЕ, НАБЛЮДЕНИЯ и ОБРАБОТКА СПЕКТРОВ
Сверхгигант p Leo (HD 91316) является относительно медленно вращающейся звездой спектрального класса B1 Iab. Эффективная температура звезды Teff определяется неуверенно: Morel et al. [11] на основе информации о спектральном подклассе приводят значение Teff = 20 260 K, в статье Gies et al. [12] по величине бальмеровского скачка и из сравнения профилей водородных линий с модельными указывают значение 24 200 K. В результате моделирования спектра звезды с учетом вклада в профили ветра звезды (Crowther et al. [13]) получено значение Teff = 22 000 K.
Звезда p Leo на диаграмме ГР находится в области переменных звезд типа ß Cep ранних подклассов спектрального класса B (см., например, [14]). Параметры звезды представлены в таблице 1: Teg — эффективная температура звезды, M — масса звезды, M — темп потери массы звездой, L — болометрическая светимость, V^ — терминальная скорость звездного ветра, V sin i — скорость вращения звезды, Prot — период вращения. В фигурных скобках приведено содержание элементов в логарифмической шкале.
Наблюдения на БТА проводились 19/20 января 2015 г. со светосильным спектрографом SCORPIO. Спектральное разрешение составляло 2000 (размер щели 0'!5), диапазон длин волн 4050-5830 A, S/N и 2000. Был получен 1271 спектр звезды. При экспозиции одного спектра от 1 до 2 секунд полное время с учетом считывания CCD составляло 6 секунд. Обработка спектральных данных была выполнена с помощью пакета MIDAS.
На рис. 1 показан нормированный на континуум участок спектра звезды p Leo в области ЛЛ 4074-5795, усредненный по 1271 наблюдению на SCORPIO.
3. ВАРИАЦИИ ПРОФИЛЕЙ ЛИНИЙ
На рис. 2 показаны вариации профилей линий H7 и Hß в спектре p Leo. Видно, что переменность профилей линий прослеживается не только для самих линий, но и для бленд линий на профиле линии H7. Профили линии Hß, полученные через
Таблица 1. Параметры p Leo
Параметр Значение Ссылка
Teff, К 24200 [12]
lg g 3.09 [12]
M/MQ TI [И]
R/Rq 37.4 [13]
Ko 1110 [13, 15]
-lg M/Mq -6.20 [И]
lg l/lq 5.18 [И]
V sin i, км с ^ 49 [16]
Prot, сут. 26.8 [7]
{С} 7.5 [13]
{N} 8.3 [13]
{0} 8.4 [13]
1332 и 1338 секунд после начала наблюдений (выделены штрих-пунктирными линиями) проходят на 2—4% выше (в единицах потока в континууме), чем остальные профили линий. Такое же превышение потока обнаружено для профилей линии HeI 4471. Для профилей линий Ш и H@ превышение потока на временах 1332 и 1338 секунд также присутствует, но амплитуда его меньше, чем для линии H@.
Наличие подобных профилей в спектре р Leo, существенно отклоняющихся от среднего, может быть свидетельством их нерегулярных вариаций на секундной шкале времени. Такие вариации могут быть связаны с высвечиванием горячего газа в неоднородностях звездного ветра после прохождения по ним ударной волны [17].
3.1. Регулярные вариации профилей линий
Для поиска регулярных компонентов вариаций профилей определим разностные профили линий. Пусть в результате проведенных наблюдений получено N спектров исследуемого объекта. Обозначим через Fi(X), i = 1,...,N — нормированный на континуум поток в i-м спектре звезды на длине волны Л. Пусть Fi{Л) — средний по всем наблюдениям поток на длине волны Л. Тогда разностный профиль линии
<Ц\)=Ъ(\)-Щ\). (1)
При анализе разностных профилей удобнее использовать вместо длины волны доплеровские смещения V от лабораторной длины волны Ло линии V = с(Л/Л0 — 1), где c — скорость света.
4200
4400
4600
4800
5000
5200
5400
5600
Wavelength, А
Рис. 1. Средний спектр р Leo
120
4328
4336
4344
4352
80
Ф
Е
40
0.06 0.05 0.04 0.03 0.02 0.01 0
-0.01 -0.02 -0.03 -0.04
-600
Velocity, km s
600 -i
4848 4856 4864 4872
o
Wavelength , A
Рис. 2. Профили линий H7 и He (сверху вниз) 1271 спектра р Leo. Отмечены бленды линии H7. Толстая светлая пунктирная линия — средний профиль линии. Штрих-пунктирными линиями отмечены профили линии H,3, полученные через 1332 и 1338 с после начала наблюдений.
При использовании спектров существенно разного качества следует при расчете среднего и разностных профилей линий использовать разные веса профилей д^, пропорциональные квадрату от-
Рис. 3. Динамический спектр вариаций профиля линии Hß.
ношения S/N в области континуума вблизи линии. Однако, так как во всех анализируемых нами профилях отношения S/N близки, можно положить 9i = 1.
На рис. 3 представлен динамический спектр вариаций профилей линии Hß в спектре р Leo. Видны регулярные изменения профилей линий со временем. Присутствуют как вариации на шкале времени порядка 20—30 минут, так и быстрые вариации на минутной шкале.
Для поиска регулярных компонентов вариации профилей линий в спектре р Leo нами был выполнен их Фурье-анализ методом CLEAN [18].
Фурье-спектры вариаций разностных профилей линий HeI 4471 и Hß (периодограммы) представ-
0
0
25
20
- 15
- 10
0.00
-800 -400 0 400 800 Velocity, km s-1
0.00
35
30
25
20
15
10
-800 -400 0 400 800 Velocity, km s -1
Рис. 4. Фурье-спектр спектр вариаций профиля линий He I 4471 (а)иН^ (b).
лены на рис. 4. В Фурье-спектрах видны регулярные компоненты с частотами = 0.031 мин-1 (Р ~ 33 мин) и = 0.014 мин-1 (Р ~ 88 мин). Такие же пики обнаружены нами в Фурье-спектрах других линий бальмеровской серии водорода и линий Не I.
Кроме того, в Фурье-спектрах вариаций профилей линий Н и Не виден сильный пик на частоте щ = 0.0045 мин-1. Соответствующий период Ро & 220 минут (3 .7) больше полного времени наблюдения 128 минут, поэтому для подтверждения реальности этого компонента вариации профилей нужны дополнительные наблюдения.
Период Р0 близок к периоду Р10 = 3 ^8 ± 0 . 1, обнаруженному при анализе вариаций профилей
линий в спектре р Leo по наблюдениям на БТА и 1.8-м телескопе Бохинсанской оптической астрономической обсерватории в 2004 г. [19].
Этот факт может свидетельствовать в пользу предположения о том, что обнаруженные нами частоты (периоды) регулярных вариаций профилей линий р Leo являются продолжением частот (периодов) вариаций, найденных по наблюдениям 2004 г., в область коротких периодов.
3.1.1. Оконное Фурье-преобразование
Проведенный нами анализ временных рядов вариаций профилей различных линий в спектре р Leo показал, что, хотя в динамических спектрах вариаций профилей линий отчетливо видны быстрые изменения профилей на шкалах 1—3 минуты, при построении периодограммы для всех спектров звезды в Фурье-спектре вариаций профилей линии обнаруживаются только регулярные компоненты с периодами больше 30 минут.
В то же время при анализе отрезков этих временных рядов в Фурье-спектрах находятся корот-копериодические компоненты. По этой причине мы использовали для поиска короткопериодических компонент вариаций профилей оконное преобразование Фурье (см., например, [20]):
<х
F(to,v) = I f(to)W(to — r)e-i2nvTdr, (2)
где Ш(¿) — оконная функция. В дальнейшем будем использовать прямоугольное окно переменной ширины: ш(¿) = ш0 Ат(¿) = 1 при г е [¿о, ¿0 + ДТ] и Ш(г) = 0 для всех остальных значений г.
При таком выборе оконной функции использование оконного преобразования Фурье эквивалентно применению Фурье-преобразования к отрезку ряда со значениями в промежутке [г0,г0 + ДТ]. Получаемый в результате применения оконного преобразования Фурье-спектр зависит от трех параметров: частоты V, г0 — начала окна и ДТ — ширины окна для каждого из значений V — доплеровского смещения от центральной длины волны линии.
Для полного профиля линии профиля мы имеем функцию четырех переменных:
Т (¿0^) = Т (¿0, ДТ,и^). (3)
Такую функцию весьма сложно анализировать, поэтому мы просуммировали значения указанной функции для всех доплеровских смещений V в пределах профиля линии:
м
Т (¿0, ДТ, V) = ^ Т (¿0, ДТ,^), (4) к=1
5
0
5
0
где суммирование выполнено для всех M значений V внутри профиля соответствующей линии.
При фиксированном значении ДТ (ширины окна) функция
F(t0,AT,u) =FAT(t0,v)
становится функцией только двух параметров, to и v. Анализ результатов применения оконного преобразования Фурье к вариациям профилей линий в спектре р Leo показал, что оптимальное значение ДТ = 20 минут.
На рис. 5 представлен результат применения оконного Фурье-преобразования к вариациям профилей бальмеровских линий в спектре р Leo для частот в интервале [0,1.2] мин-1.
В нижней части оконных Фурье-периодограмм видны регулярные компоненты вариаций профилей, обнаруженные при анализе 1271 спектра р Leo. В то же время, кроме этих компонент, обнаруживаются и новые. Наиболее ярко выражен компонент с переменной частотой v = 0.4 мин-1 (период P ~ 2.5 минут) при t0 = 0 и v = 0.5 мин-1 при t0 = 100 минут (P ~ 2 минуты). В области t0 & 70 минут данный компонент раздваивается. Нижняя ветвь примерно соответствует частотам вариаций при t0 = 0.
Кроме того, на периодограммах виден менее выраженный компонент с переменной частотой v = 0.2 мин-1 (P ~ 5 минут) при t0 = 0 и v = 0.1 мин-1 (P ~ 10 минут) при t0 = 110 минут. Возможно также присутствие компонента c v ~ 0.75 мин-1 c периодом P ~ 1.3 минут. Однако этот компонент прослеживается только до значения t0 < 10 минут.
На рис. 5 видны также слабые детали в области v > 0.8 мин-1, однако для подтверждения их реальности требуется дополнительное исследование. Такие же, как описано выше, регулярные компоненты обнаружены нами в оконных Фурье-спектрах линий He I, Fe II, N II, О II и Si II.
4. РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ
Обнаруженные по наблюдениям профилей спектров р Leo регулярные вариации с периодами P = 4h — 6h связаны с нерадиальными пульсациями звезды в квадрупольной (l = 2) и более высоких модах пульсаций [19]. Обнаружение изменений профилей линий в области коротких периодов P = 30—70 минут может свидетельствовать о присутствии высоких мод нерадиальных пульсаций с l = 6-12.
Необычной особенностью звезды р Leo является высокая скорость макротурбулентных движений
1-0
0 20 40 60 80 100 Frame initial time, min
0 20 40 60 80 100 Frame initial time, min
20 40 60 80 100 Frame initial time, min
Рис. 5. Оконный Фурье-спектр вариаций профиля линии HJ, H7 и Hß (сверху вниз) для значения AT = 20 минут в зависимости от положения начала окна t0 (initial time).
в фотосфере Vm = 72 км с 1 [7] при том, что скорость вращения V sin i ~ 50 км с-1 (см. таблицу 1). Такое высокое значение макротурбулентной скорости может быть объяснено присутствием крупномасштабных тангенциальных движений [21 ]. Обнаруженные нами регулярные вариации профилей с периодами P = 30-70 минут могут быть связаны с такими движениями.
Использование оконного Фурье-преобразования показало наличие в вариациях профилей линий в спектре р Leo короткопериодических регулярных компонент с переменной частотой. Такие компоненты ранее не были детектированы и природа их не вполне ясна.
Присутствие компонент вариаций профилей линий с короткими периодами в спектре H 93521 было нами ранее интерпретировано как свидетельство существования высоких мод NRP с l = 20-60 [6]. Однако присутствие высоких мод NRP не объясняет переменность частот регулярных компонент вариаций профилей.
Возможно, однако, что наличие таких компонент в регулярных вариациях профилей линий может быть объяснено при предположении, что высокие моды нерадиальных пульсаций нестабильны и могут как затухать, так и генерироваться на коротких шкалах времени порядка 10-100 минут. В случае компоненты с переменной частотой v = 0.4-0.5 мин-1 мы, возможно, видим затухание компоненты v = 0.4 мин-1 и ее новое появление при t0 & 70 минут и генерацию компоненты v = 0.5 мин-1 при t0 > 20 минут.
Особый интерес представляет собой присутствие в спектре р Leo профилей линий, существенно отклоняющихся от среднего. Появление подобных профилей в спектре говорит о возможном вкладе слабой эмиссионной компоненты в полный профиль линии и может быть свидетельством нерегулярных вариаций профилей на секундной шкале времени. Такие вариации могут быть связаны с высвечиванием горячего газа в неоднородностях звездного ветра после прохождения по ним ударной волны.
Используя функцию высвечивания горячей плазмы в области температур Т = 106-107 K [22], можно оценить время высвечивания т для звездного ветра р Leo c n ~ 1011 -1012 см-3. Величина т составляет 10-100 секунд, что близко к характерному времени нерегулярных вариаций профилей в спектре р Leo. Отметим также, что нерегулярные вариации профилей линий в спектре р Leo могут быть связаны с нестабильностью конфигурации вещества звездного ветра, захваченного в области магнитного экватора магнитным полем р Leo (см., например, [23]).
5. ЗАКЛЮЧЕНИЕ
В настоящей работе исследованы вариации
профилей линий в спектре звезды р Leo с высоким
временным разрешением. Обнаружены регулярные
компоненты вариаций профилей с периодами
примерно 30 и 90 минут. Возможно присутствие
быстрых нерегулярных вариаций на шкале времени менее 1 минуты.
Обнаружены высокочастотные (v > 0.1 мин-1) компоненты вариаций профилей с переменной частотой. Мы предполагаем, что наличие таких компонент может быть связано с нестабильностью высоких мод нерадиальных пульсаций, их генерацией
и затуханием на коротких временных шкалах в области 1 — 10 минут.
Отметим, что вариации профилей линий в спектре звезды могут быть обусловлены спорадическим изменением размера турбулентного диска звезды вследствие атмосферных флуктуаций. При этом вызванные их влиянием вариации профилей вряд ли будут регулярными, однако вклад таких флук-туаций в нерегулярные вариации профилей нельзя полностью исключить. Для оценки этого влияния предполагается выполнение одновременных наблюдений спектра звезды сравнения, расположенной на щели вблизи от исследуемой звезды.
БЛАГОДАРНОСТИ
Настоящее исследование выполнено за счет гранта Российского научного фонда (проект № 1812-00423).
СПИСОК ЛИТЕРАТУРЫ
1. L. Kaper, H. F. Henrichs, A. W. Fullerton, et al., Astron. and Astrophys. 327, 281 (1997).
2. A. F. Kholtygin, D. N. Monin, A. E. Surkov, and S. N. Fabrika, Astronomy Letters 29, 175 (2003).
3. V. V. Dushin, A. F. Kholtygin, G. A. Chuntonov, and D. O. Kudryavtsev, Astrophysical Bulletin 68, 184 (2013).
4. S. Hubrig, M. Scholler, and A. F. Kholtygin, Monthly Notices Royal Astron. Soc. 440, 1779 (2014).
5. V. L. Afanasiev and A. V. Moiseev, Astronomy Letters 31, 194(2005).
6. A. F. Kholtygin, S. Hubrig, V. V. Dushin, et al., ASP Conf. Ser. 510, p. 299 (2017).
7. C. Aerts, D. M. Bowman, S. Simon-Diaz, et al., Monthly Notices Royal Astron. Soc. 476, 1234 (2018).
8. A. F. Kholtygin, G. A. Chountonov, S. N. Fabrika, et al., in Proc. Conf on Physics of Magnetic Stars, Ed. by I. I. Romanyuk, D. O. Kudryavtsev, O. M. Neizvestnaya, and V. M. Shapoval (Spec. Astrophys. Obs., Nizhny Arkhyz, 2007), 262.
9. D. N. Monin, S. N. Fabrika, and G. G. Valyavin, Astron. and Astrophys. 396, 131 (2002).
10. A. F. Kholtygin, S. N. Fabrika, N. A. Drake, et al., Astronomy Letters 36, 370 (2010).
11. T. Morel, S. V. Marchenko, A. K. Pati, et al., Monthly Notices Royal Astron. Soc. 351, 552 (2004).
12. D. R. Gies and D. L. Lambert, Astrophys. J. 387, 673 (1992).
13. P. A. Crowther, D. J. Lennon, and N. R. Walborn, Astron. and Astrophys. 446, 279 (2006).
14. A. A. Pamyatnykh, Acta. Astron. 49, 119 (1999).
15. I. D. Howarth, K. W. Siebert, G. A. J. Hussain, and R. K. Prinja, Monthly Notices Royal Astron. Soc. 284,265(1997).
16. S. Simon-Diaz and A. Herrero, Astron. and Astrophys. 562, A135 (2014).
17. F. A. Kholtygin, C. J. Brown, P. J. Cassinelli, et al., Astron. Astrophys. Transactions 22, 499 (2003).
18. D. H. Roberts, J. Lehar, and J. W. Dreher, Astron. J. 93,968(1987).
19. A. F. Kholtygin, S. N. Fabrika, T. E. Burlakova, et al., Astronomy Reports 51, 920 (2007).
20. И. Добеши, Десять лекций по вейвлетам, (НИЦ «Регулярная и хаотическая динамика», Ижевск, 2001).
21. C. Aerts, J. Puls, M. Godart, and M.-A. Dupret, Astron. and Astrophys. 508, 409 (2009).
22. E. Landi and M. Landini, Astron. and Astrophys. 347, 401 (1999).
23. A. Ud-Doula, S. P. Owocki, and R. H. D. Townsend, Monthly Notices Royal Astron. Soc. 392, 1022 (2009).
Super-Fast Line-Profile Variability in the Spectra of OBA-Stars: B1-Star p Leo
A. F. Kholtygin, A. A. Batrakov, S. N. Fabrika, A. F. Valeev, I. M. Tumanova, and O. A. Tsiopa
This study continues our investigation of the super-fast variability of line profiles in the spectra of early-type stars. We have investigated the line-profile variability in spectra of OBA-stars with the SCORPIO multi-mode focal reducer mounted at the 6-m SAO RAS telescope. Regular short-period variations of the H and He lines in the spectra of the B1I star p Leo were detected with periods ranging from 2 to 90 minutes, as well as irregular line-profile variations on time intervals less than 1 minute. A possible origin of rapid spectral variations is discussed.
Keywords: stars:early-type—stars: individual: pLeo
