CC BY
16
УДК 524.38-327:520.84
ЗВЕЗДЫ С РАЗЛИЧАЮЩИМИСЯ ВЕЛИЧИНАМИ v sin i, ОПРЕДЕЛЕННЫМИ ПО ЛИНИЯМ Ca II 3933 И MgII 4481 A. VII. HD9531 (SB), HD31592 (SB2), HD 129174 (SB?)
©2018 Ю. Зверко12*, И. Илиев3, И. И. Романюк4, И. Статева3, Д. О. Кудрявцев4, Е. А. Семенко4
1Институт теоретической физики и астрофизики, Татранска Ломница, 05960 Словакия 2Масариков университет, Брно, 60177 Чешская Республика 3Институт астрономии, НАО Рожен, София, 1784 Болгария 4Специальная астрофизическая обсерватория РАН, Нижний Архыз, 369167 Россия Поступила в редакцию 4 апреля 2018 года; принята в печать 14 мая 2018 года
В данной работе серии мы анализируем три объекта из списка звезд с различающимися значениями v sin i: HD9531 и HD 31592, которые демонстрируют переменность лучевых скоростей, характерную для спектральной двойственности, а также HD 129174, которая является звездой класса пекулярности Mn с возможными магнитным полем. Мы подтверждаем переменность лучевой скорости, полученной для HD9531 по водородным линиям, а также по линии CaII на 3933 A . Профиль линии кальция также переменный, что, учитывая оцененное магнитное поле Be = —630 ± 1340 Гс, свидетельствует в пользу того, что кальций неравномерно распределен по поверхности звезды. Мы отождествили линии вторичного компонента в спектре HD 31592, показывая таким образом, что она является двойной звездой типа SB2 с компонентами классов B9.5V и A0V. Главный компонент вращается со скоростью v sin i = 50 кмс-1, а вторичный — быстрее: v sin i = 170 кмс-1. Только 60% линий Mn, идентифицированных в спектре HD 129174, удается промоделировать единственным значением содержания; остальные линии или сильнее, или слабее модели. Мы также обнаружили две линии Xe II на 5339.33 A и 5419.15 A и оценили их log gf.
Ключевые слова: звезды: двойные — звезды: вращение — звезды: индивидуальные: HD9531, HD 31592, HD129174
1. ВВЕДЕНИЕ
В данной серии работ мы исследуем спектры звезд, проекции скоростей вращения которых отличаются в зависимости от того, какие спектральные линия использовались для определения. Мы сравнили каталог Палмера и др. [1], основанный на линии ионизованного кальция на 3933 A, и каталог Вольффа и Престона [2] с данными
по линии ионизованного магния на 4481 A. В предыдущих шести работах мы разобрали 17 из 24 звезд выборки. Семь из девяти звезд с v sin i (CaII) ^ v sin i (MgII) [3] имеют малые значения проекции скоростей вращения вследствие наличия узкой межзвездной компоненты; одна из двух оставшихся звезд (HD 44783) оказалась Be-звездой, другая (HD 25152) окружена околозвездным веществом. Из восьми оставшихся объектов
*E-mail: zve@ta3.sk
HD 8837 является звездой с оболочкой, остальные семь — двойные [4—8].
В настоящей работе мы анализируем три звезды с возможными признаками вторичных компонентов в спектрах.
2. НАБЛЮДЕНИЯ
ПЗС-спектры с высоким отношением S/N получены в Болгарской национальной астрономической обсерватории Рожен (НАО) и в Специальной астрофизической обсерватории Российской академии наук (САО), также один спектр был извлечен из публичного архива ELODIE [9]. Спектры НАО были получены с помощью ПЗС камеры Photometries AT200 1024 х 1024 на третьей оптической камере куде-спектрографа установленного на 2-м телескопе системы Ричи—Кретьена с разрешением R = 22 000. Область Ca II находится в диапазоне от 3898 до 3967 A, а область Mg II — от
4447 до 4550 A. Для обработки спектров использовались стандартные процедуры IRAF1. Спектры САО были получены на 6-метровом телескопе БТА с помощью эшелле-спектрометра фокуса Нэсмита (НЭС), оборудованного ПЗС камерой 2048 х 2048 [10], с разрешением R = 43 000 в спектральной области 4200—6700 A. Обработка осуществлялась в пакете программ REDUCE [11].
Кроме того, мы анализировали фотографические спектры, полученные на куде-спектрографе 2-м телескопа Астрономического института Академии наук Чешской Республики в городе Ондржеев. Спектральный диапазон составлял Л3650—4900 A, снимки регистрировались на фотоэмульсии Kodak IIaO или IIaOb, величина обратной линейной дисперсии была равна 8.5 A мм-1. Фотографические спектры были оцифрованы в Астрономическом институте Академии наук Словакии в Татранской Ломнице с помощью микроденситометра и обработаны кодом SPEFO (Шкода [12]), который позволяет построить дисперсионную и калибровочную кривые, приведение к шкале интенсивностей и нормализацию континуума. Определение лучевой скорости было произведено методом кросс-корреляции (CCF) [13]. Исправление спектров было проведено таким образом, что точки максимальной интенсивности между взаимно пересекающимися крыльям бальмеровских линий выбирались к качестве точек континуума. Тот же метод исправления был применен к синтетическим спектрам. Такая процедура позволяет избежать отклонений, вызванных прогрессирующим затуханием интенсивности спектра в направлении увеличения номера линий бальмеровской серии, к которому метод CCF чувствителен. После применения этой процедуры значение лучевой скорости, полученное кросс-корреляцией синтетических и фотографических спектрограмм, соответствует системе бальме-ровских линий как доминантной черты спектра.
Код SYNSPEC [14, 15] использовался для расчета синтетических спектров и детальных профилей линий, а также для определения содержания химических элементов путем сравнения с наблюдаемыми линиями. ЛТР-модели атмосфер были получены путем интерполяции сетки моделей Кастел-ли и Куруца [16]. В дополнение к "списку линий gf0800", используемому в коде SYNSPEC, для получения данных по атомным линиям мы использовали базы VALD [17-20] и NIST [21], а также
'¡RAF распространяется Национальной обсерваторией
оптической астрономии (NOAO), под управлением Ассоциации университетов для исследований в астрономии
(AURA), в рамках совместного соглашения с Национальным научным фондом (NSF).
работу Вейса и Фюра [22]. Значения эффективных температур и поверхностных тяготений были получены с помощью кодов UVBYBETA Муна и Дворецки [23] и TEFFLOGG Смолли и Дворец-ки [24] по данным uvbyв-фотометрии, взятыми из базы данных SIMBAD. Поправки содержаний, где они указаны, относятся к солнечным [25].
Суммарные спектры двойной были смоделированы путем совмещения спектров индивидуальных компонентов с помощью кода SUMSPEC, предоставленного Р. Комжиком [26].
3. АНАЛИЗ 3.1. HD 9531
HD9531 (HR446, HIP321, KKAnd, B9IV, V = 5.91). Хоффлайт и Варрен [27] классифицируют объект как CP-звезду класса a2 CVn с аномалиями Si, обладающей переменностью с коротким периодом в 0.6684 дня и амплитудой 0.012 величины в полосе V, а также отмечают долговременную вариацию периода. Палмер и др. [1] оценили v sin i = 215 км с-1 по линии Ca II на 3933 A, а Вольфф и Престон [2], которые использовали линию MgII на 4481 A, получили sin i = 175 кмс-1. Ройер и др. [28] приводят значение v sin i = 176 км с-1.
Значения лучевой скорости, приведенные в литературе, меняются от —20 до +10 кмс-1. Хоффлайт и Варрен [27] приводят RV = —4 кмс-1. Абт [29] приводит пять значений от —20 до 8.9 кмс-1, наблюдаемых между 1918 и 1921 годами. Вольфф [30] измерил семь значений между —4.6 и +13.4 кмс-1 и отметил, что это SB звезда. Хрик [31] получил значения от —9.8 до +7.8 км с-1 путем измерения водородных линий в 21 спектре и утверждает, что лучевая скорость постоянна; однако, линии ScI, SiII, SrII, SII и PbI могут отражать период близкий к 0d7, найденный ранее Хильдебрандтом [32].
В настоящей работе мы используем спектр, полученный на ОЗСП БТА, охватывающий область MgII 4481 A, спектры НАО, покрывающие область Ca II 3933 A, а также один спектр ELODIE в диапазоне длины волн от 4200 A до Ha. Набор фотографических спектров был получен на 2-м телескопах обсерваторий Ондржеев и Рожен, три из них, полученные в обсерватории Ондржеев, были оцифрованы нами. Спектры, использованные в настоящей работе, приведены в Таблице 1.
Значения эффективной температуры, приведенные в литературе, варьируются в широком интервале от Teff = 10 400 K [33] до Teff = 14 000 K [31].
ЗВЕЗДЫ С РАЗЛИЧАЮЩИМИСЯ ВЕЛИЧИНАМИ V SIN I Таблица 1. Список спектров и измеренные лучевые скорости HD9531
Спектр HJD 2400000+ S/N Область RV, KM c_1 Спектр HJD 2400000+ S/N Область RV, км с-1
Ondrejov 42332.5174 42 H10-H/3 +5.3 ±4.5 NAO 56586.4050 96 3933 Á -0.5 + 4.1
44159.4342 34 Hll-H/3 -6.2 + 4.6 56587.4283 85 3933 Á -1.1 + 3.2
44819.5338 45 Hll-H/3 -11.3 + 4.8 56588.4166 70 3933 Á -3.4 + 3.3
ELODIE 53320.5132 243 4200—Ha +2.8 + 2.2 56613.3577 89 3933 Á -4.2 + 3.7
MSS 55075.4071 772 4481 Á +4.0 + 2.8 56640.2675 137 3933 A -2.1 + 3.7
NAO 56527.5612 70 3933 Á -8.4 + 3.6 56641.2326 116 3933 A -8.9 + 3.9
56528.5570 80 3933 Á -8.4 + 3.9 56642.2309 97 3933 A -13.7 + 4.2
Значения поверхностного тяготения лежат в узком интервале, а именно от log g = 4.0 [31] до logg = 4.19 [34].
Мы использовали индексы uvby/З из базы данных SIMBAD и получили Teff = 12 040 K и log g = 4.06.
Синтетический спектр, построенный для модели атмосферы с Teff = 12 040 K и log g = 4.06, нулевой микротурбулентностью и v sin i = 165 км с-1, хорошо описывает линии MgII на 4481 A, HeI на 4471 A и Ca II на 3933 A. При расчете спектра содержание магния было увеличено в 2.5 раза, а гелия — уменьшено в 0.9 раз.
В фотографических спектрах и спектре, взятом из архива ELODIE, доминируют бальмеровские линии, в то время как линии металлов размыты вследствие высокой скорости вращения; таким образом, лучевые скорости, полученные методом CCF, относятся к бальмеровским линиям.
В спектре ОЗСП с высоким отношением S/N присутствуют линии HeI 4471 A и MgII 4481 A, а также до 20 слабых линий с глубиной менее 0.02, которые могут быть использованы для определения лучевой скорости. Спектры в области линии MgII на 4481 A показаны на Рис. 1.
В области линии Ca II 3933 A отсутствуют различимые спектральные линии. Сама линия имеет относительную интенсивность примерно 0.06, а ее профиль искривлен на четырех из девяти спектров. Таким образом, использование метода CCF для оценки лучевой скорости не приведет к однозначному результату. Более того, низкая чувствительность ПЗС-матрицы в УФ области делает ниже отношение S/N, что проиллюстрировано на Рис. 2.
В пяти из девяти спектров наблюдаемые профили линий соответствуют синтетическим. Профили третьего, восьмого и девятого спектров на Рис. 2,
если считать сверху, искажены в левой части, а четвертого — в правой. Поэтому мы аппроксимируем каждый профиль графически, используя его неискаженную часть для определения лучевой скорости. После чего узкая деталь около центра Са II в каждом спектре приводилась к одной длине волны. Мы обнаружили, что в спектре И0 9531, полученном на ELODIE, присутствуют линии межзвездного натрия D1 и D2, таким образом, особенность в линии кальция также имеет межзвездную природу и может быть взята в качестве репера для измерения лучевой скорости по линии звезды. Лучевая скорость, измеренная по межзвездной линии, составила +0.55 км с-1.
Сравнение значений 3<г-интервала с величиной искажений показывает, что последние могут представлять собой реальные изменения. Наблюдается фотометрическая переменность с периодом 0.6684 дня, согласно Хильдебрандту [32]. Отметим, однако, что искажения не совпадают с фазами периода.
Искажения в профиле звездной линии Са II 3933 А могут указывать на связь с фотометрической переменностью, открытой Хильдебрандтом [32], и возникающей вследствие наличия магнитного поля, которое которое является причиной неравномерности в распределении химических элементов по поверхности звезды. Панов и Шонайх [35] исследовали HD9531 как магнитную звезду, а мы оценили величину продольного магнитного поля Ве = —630 ± 1340 О. Пятна повышенного содержания химических элементов проявляют себя в виде узких депрессий, гуляющих по профилю линии. В случае, показанном на Рис. 2, однако, эта гуляющая деталь больше похожа на небольшую эмиссию, которая может быть вызвана относительно большим поверхностным пятном с экстремально низким содержанием единожды ионизованного кальция.
Рис. 1. Область вокруг линии MgII в спектре HD 9531. Спектры сверху вниз: ELODIE, три фотографических и ОЗСП. На наблюдаемые спектры наложен синтетический (тонкая сплошная линия). Каждый наблюдаемый спектр имеет обозначения, слева направо: HJD, лучевая скорость и S/N, а также сдвинут по шкале длин волн на величину наблюдаемой лучевой скорости.
Wavelength, А
Рис. 2. Область вокруг линии Ca II в спектре HD9531. На каждом спектре слева отмечены HJD, справа приведены фазы, соответствующие периоду в 0.6684 дня согласно эфемериде [32]. Усы погрешности показывают уровень 3а соответствующий отношению S/N в континууме вокруг линии Ca II. Точки показывают положения межзвездной компоненты.
Мы делаем заключение, что значения приведенные в Таблице 1 указывают на вероятную переменность лучевой скорости HD9531. Это становится очевидным при рассмотрении линий водорода в длинных спектрах, несмотря на то, что ошибки фотографических спектрограмм относительно высоки. В случае с ультрафиолетовыми спектрами НАО, где ошибки в два раза меньше, чем у фотографических спектров, возможная переменность лучевой скорости подтверждается нашими оценками по неискаженным профилям линии Ca II
на 3933 A, в которых присутствует межзвездная компонента, имеющая постоянную лучевую скорость. Переменная лучевая скорость указывает на то, что звезда принадлежит спектрально двойной системе. Однако, мы не нашли следов вторичного компонента даже в спектре ОЗСП с высоким отношением S/N = 772 (3а ~ 0.004), что соответствует отношению интенсивностей компонентов примерно 260 и разности блеска около 6 величин. Это указывает, что невидимый второй компонент может быть звездой главной последовательности позднего класса G.
Периодическая фотометрическая переменность, а также наличие магнитного поля у HD 9531 также согласуются с моделью, в которой переменность профилей металлических линий объясняется существованием пятен с повышенным или пониженным содержанием химических элементов. Авторам не удалось найти доступных наблюдений временных рядов с высоким отношением S/N, охватывающих синюю область спектра с металлическими линиями, которые бы подходили для проведения исследования, направленного на подтверждение заключений Хрика [31]. Неискаженный профиль линии Ca II 3933 A моделируется нормальным содержанием кальция, поэтому отклонения четырех профилей линий, видимые на наших спектрах НАО, указывают скорее на пятно с пониженным содержанием ионизованного кальция. Следует отметить, что сильная вертикальная стратификация кальция, вызванная радиативной диффузией, ранее наблюдалась только в холодных CP звездах [36, 37]. В присутствии магнитного поля этот процесс может привести к появлению пятен химических элементов. Эффективная температура HD9531, полученная в настоящей работе, однако, выше, Teff = 12 040 K, хоть и значительно ниже значения Teff = 10 406 K, приведенного Макдо-нальдом и др. [33].
Разница между двумя значениями v sin i, полученными по двум разным линиям, Ca II 3933 A и MgII 4481 A, вызвана слабой интенсивностью самих линий. А именно, линия магния, имеющая
глубину примерно 0.1 с отношением S/N фотометрической эмульсии приблизительно 35—45, относительно хорошо различима, а линия кальция на наших фотографических спектрах имеет глубину только около 0.05—0.07, поэтому плохо оценивается не только ширина линии, даже обнаружить ее в спектре довольно сложно. Таким образом, разница в 40 км с-1 между этими двумя значениями очевидна.
3.2. HD 31592
HD 31592 (k Tau, HR 1590, HIP 23088, ADS 3547 CCDM J04581+2503A, A0V, V=5.79) является звездой в двойной системе. Палмер и др. [1] получили оценку v sin i = 170 км с-1 по линии Ca II на 3933 A, а Вольфф и Престон [2], которые использовали линию MgII на 4481 A, приводят sin i = 50 кмс-1. Доммангет и Нис [38] сообщают о наличии трех визуальных компонент, где компонента B (зв. вел. 13.3) и C (зв. вел. 9.9) удалены от компоненты A на 41.8 и 94.6 угловых секунд соответственно. Шмидтке и др. [39] обнаружили методом лунного покрытия, что компонент A сам является двойной, с компонентой на 0.22 величины слабее. Хоффлайт и Варрен [27] отмечают, что звезда является спектрально двойной.
Имеющиеся у нас спектры перечислены в Таблице 2.
В литературе может быть найдена альтернативная спектральная классификация звезды. Абт и Моррелл [40] приводят Am/B9.5/A0/A1. Хоффлайт и Варрен [27] определяют звезду как A0V, а Шмидтке и др. [39] — B9V.
Для определения эффективной температуры и поверхностного тяготения мы использовали индексы uvby/З, приведенные Кроуфордом и др. [41], и получили Teff = 10120 к и logg = 4.46. Существенно более низкая эффективная температура, Teff = 9132 K, была получена Макдональдом и др. [33], которые использовали метод моделей атмосфер для построения спектрального распределения энергии звезды. Хайнтце [42] приводит шкалу температур, полученную по индексам (U — — B). Используя два набора индексов UBV, приведенных в базе SIMBAD, и следуя трем формулам из вышеупомянутой работы (25a-c), мы вычислили Teff = 10 072 до 10 290 K при (U — B) = —0.04 и Teff = 9885 до 10 109 K при (U — B) = —0.01. Мы использовали модель атмосферы, полученную с помощью индексов uvby/З, и приближенную к средним значениям, полученным по индексу (U — B). Для наилучшего соответствия синтетического и наблюдаемого спектров содержания гелия и магния были уменьшены в 0.1 и 0.15 раз, а содержания
Таблица 2. Список спектров и измеренные лучевые скорости HD 31592
Спектр HJD S/N Область RV главная, RV вторичная,
2400000+ км c_1 км с-1
MSS 55075.5803 >500 MgII 4481 Á 44.5 + 0.4 0 + 4
NAO 56588.4457 102 Ca II 3933 Á 24.5 + 2.2 -0.7 + 4.6
56613.4136 86 16.7+1.3 25.6 + 3.7
56614.4771 108 13.2 + 1.6 34 + 4.6
56640.4471 135 -5.1 + 2.3 29 ::
56641.4409 133 -8.0 + 2.7 42 ::
Рис. 3. Область вокруг линии М^П в спектре HD 31592. Сверху вниз: синтетический спектр вторичного и первичного компонентов. Первичный и вторичный спектры сдвинуты по лучевой скорости на 44.5 км с-1 и 0 км с-1 соответственно. Суммарный спектр (внизу) наложен поверх наблюдаемого. Спектры сдвинуты по вертикали для наглядности.
кальция и железа были увеличены в 15 раз. Чтобы промоделировать линию кальция на 3933 A, понадобилось уменьшить его содержание в звезде в 0.01 раз.
Исходя из спектральной классификации MK [43], эффективная температура, полученная в настоящей работе, соответствует спектральному классу B9.5. Разница в спектральном классе второго компонента, который слабее на 0.22 величины, соответствует половине подтипа, следовательно, его спектр может быть определен как A0. Модель с Teff = 9720 Ки log g = 4.45 была использована для расчета спектров вторичного компонента. Наилучшее соответствие синтетического и наблюдаемого спектров вторичного компонента достигается при пониженном в 0.5 раз содержании кальция.
На Рис. 3 и 4 продемонстрирована часть спектра с ОЗСП. Асимметрично уширенные линии, как
и широкая депрессия между 4465 A и 4475 A, хорошо описываются после наложения двух спектров. Синтетические спектры главного и вторично-
го компонентов показаны сверху вниз, в положениях +44.5 км с-1 и 0 км с-1 соответственно. Мы использовали значения v sin i = 55 км с-1 для главного компонента и v sin i = 170 км с-1 — для вторичного. Третий профиль представляет собой наложение наблюдаемого и синтетического спектров, сложенные в соотношении /prim//sec = 0.551/0.449 для компенсации разности блеска. Результат хорошо согласуется с наблюдаемым спектром.
Часть спектра вокруг линии Ca II показана на Рис. 5. Пять наблюдаемых спектров, суммированные по соотношению интенсивностей и выровненные по индивидуальным значениям лучевой скорости, наложены на суммарный синтетический спектр и представлены на нижнем профиле. Синтетические спектры главного и вторичного компонентов в шкале лабораторных длин волн представлены сверху. Согласие между наблюдаемыми и сложенными синтетическими спектрами в целом удовлетворительно, за исключением области на длине волны 3930 A. На этой длине волны образуется линия Fe I, поэтому мы рассчитали спектр главного
Wavelength, A
Рис. 4. Часть спектра ИЭ 31592, полученного на ОЗСП. Главный (в середине) и вторичный (сверху) спектры сдвинуты по лучевым скоростям на 44.5 км с-1 и 0 км с-1 соответственно. Суммарный спектр (снизу) наложен на наблюдаемый. Спектры сдвинуты по вертикали для наглядности.
2456588.4457 2456613.4136 2456614.4471
3915 3920 3925 3930 3935 3940 3945 3950 3955
Wavelength, А
Рис. 5. Область вокруг линии Ca II в спектре HD 31592. Сверху вниз: синтетический спектр главного компонента ("горячий" — толстая сплошная линия, "холодный" — пунктирная линия). Спектры сдвинуты на соответствующее значение v sin i и уменьшены по интенсивности в соотношении /primary : /secondary = 0.551 : 0.449. Пять сложенных наблюдаемых спектров (толстая сплошная линия) снизу показаны вместе с суммарным синтетическим спектром (пунктирная линия) и со вторичным компонентом (тонкая сплошная линия); все спектры приведены к одной лучевой скорости.
компонента, используя более низкое значение эффективной температуры из работы [33], вследствие чего линия нейтрального железа должна стать более интенсивной. Этот холодный спектр наложен на более горячий вверху рисунка, однако разница между ними оказалась не существенной. Таким образом, источник абсорбционной детали на 3930 A остается неясным.
Синтетический спектр вторичного компонента, свернутый со значением v sin i = 170 кмс-1 и
уменьшенный соответственно в 0.45 раз, также показан на Рис. 5. Каждый из них сдвинут по лучевой скорости для достижения наилучшего визуального соответствия, которое, однако, не совпадает со значением, полученным методом CCF для двух нижних спектрах. В любом случае, неопределенность лучевой скорости вторичного компонента
возникает вследствие наличия особенности в линии кальция. Не только абсорбционная деталь на
3933 A, но и линия Fe II около 3936 A вносят вклад в наблюдаемый профиль.
Таким образом, мы показали, что заключения Палмера и др. [1], которые получили оценку
v sin i = 170 км с-1 по линии Ca II на 3933 A, и Вольффа и Престона [2], которые измеряли
положение линии MgII на 4481 A и получили sin i = 50 км с-1, согласуются с наблюдениями, так как в спектре присутствуют линии с двумя различными значениями v sin i. Линия магния в спектре главного компонента узкая и глубокая, а линия кальция в спектре главного компонента с уменьшенным содержанием слаба, над ней преобладает линия вторичного компонента.
3.3. HD129174
HD129174 (п1 Boo A, HR5475, HIP7176, ADS9338A, CCDMJ14407+ 1625A) — яркая, очень хорошо изученная CP звезда, классифицированная как B9IIIpHgMnSi [44]. Она является более ярким компонентом (V=4.94 величины [27]) в визуальной двойной системе, где более слабый компонент B (V=5.8 величины [27]) находится на расстоянии 5.6 угловых секунд.
Мы включили эту звезду в наблюдательную программу из-за различий между двумя определениями проекции скорости вращения: Палмер и др. [1] получили v sin i = 60 км с-1 по линии Ca II на 3933 A, а Вольфф и Престон получили v sin i = 16 км с-1 [2] по линии MgII на 4481 A. Список спектров приведен в Таблице 3.
3.3.1. Параметры атмосферы
Основываясь на индексах uvby/З из базы SIMBAD, мы оценили эффективную температуру Teff = 12 330 K и поверхностное тяготение logg = 3.45. В многочисленных исследованиях разные авторы приводят значения Teff от 11 000 до 14 400 K, используя фотометрию [45—47], спектро-фотометрию [48], профили водородных линий [49], модели атмосфер [50], степени ионизации Mg и Fe [51]. Таким образом, наши модельные параметры лежат около нижней границы температуры и тяготения по разным измерениям. Широкий диапазон оценок, ±1700 K от среднего значения, и относительно широкий диапазон поверхностных тяготений, от 3.45 (в настоящей работе) до 4.15 [46], вместе с оцененными значениями микротурбулентности между 3 км с-1 [52] и 7 км с-1 [49] и возможное наличие магнитного поля —75 Гс и +190 Гс, упомянутое Монтгомери [51], которое нетипично для химической пекулярности типа Mn, подчеркивают необычные свойства HD 129174.
3.3.2. Область линии Ca II на 3933 A
Мы рассчитали синтетический спектр, используя модель атмосферы с Teff = 12 330 K и log g = 3.45, а также микротурбулентностью £ = 1 км с-1. Кроме кальция, линий HeI (на 3926.5 A), Al 11 (3900.6 A), хрома (3905.6 и 3911.3 A), TiII (3900.6 и 3913.5 A), CII (3918.9 и 3920.6 A), Y (3950.3 A), присутствуют и другие линии MnII и FeII. Содержание марганца оказалось увеличенным в 300 раз, а иттрия — в примерно 2000; недостаток гелия составляет приблизительно 0.05, алюминия — примерно 0.05, а содержание магния уменьшено в 0.5 раз. Хром слегка избыточен, его содержание увеличено в 5 раз, а содержание титана — в 1.8 раз. Содержание железа также было уменьшено в 0.8 раз для получения наилучшего соответствия интенсивностей линий. Для того, чтобы промоделировать линию Ca II на 3933.6 A, нам пришлось увеличить содержание кальция в 1.7 раз. Наилучшее значение проекции скорости вращения составило v sin i = 18 км с-1.
По сравнению со спектром, полученным в оптическом диапазоне (например, спектры NES), ультрафиолетовые спектры вокруг линии Ca II на 3933 A имеют гораздо более низкие отношения S/N. Поэтому мы сложили пять индивидуальных спектров для получения одного с лучшим S/N. В этой процедуре спектры, которые предстояло сложить, сначала интерполировались эквидистантными шагами по лучевой скорости (здесь 0.079 км с-1), затем кросс-коррелировались для нахождения наилучшего сдвига, а уже потом суммировались. Таким образом, мы получили спектр с S/N œ 220. На Рис. 6 показана узкая область вокруг линии Ca II на 3933 A. Дополнительная абсорбция четко видна относительно синтетического профиля в левом крыле линии. Эта абсорбция может имитировать большую ширину линии на фотографических спектрах, что в свою очередь приводит к переоценке v sin i (например, у Палмера
и др. t1]).
Положение этой детали, по-видимому, распространяется от более длинных (верхний спектр) к более коротким волнам. В то же время верхний спектр, похоже, имеет еще одну слабую абсорбционную деталь у основания правого крыла линии, что похоже на гуляющую слабую и более широкую линию, которая является следствием присутствия некого компонента двойной системы. Это заключение, однако, не подтверждается наблюдаемыми глубинами спектральных линий, которые, в таком случае, должны быть в целом менее глубокими из-за дополнительного блеска от компаньона. (Похожая, но более слабая деталь также видна на пятом
ЗВЕЗДЫ С РАЗЛИЧАЮЩИМИСЯ ВЕЛИЧИНАМИ V SIN I Таблица 3. Список спектров и измеренные лучевые скорости HD 129174
Спектр HJD 2400000+ S/N Область RV, км с-1
NES 54964.2641 250 4228-5655 Á -1.48
NAO 56375.4372 128 Ca II 3933 Á -4.99 + 0.79
56410.3811 101 -5.11 + 0.83
56411.4136 130 -3.13 + 0.90
56412.4063 109 -0.92 + 0.68
56437.3815 78 -6.52 + 0.83
1.2 1
0.8
и 0.6
с
0
■Е 0.4 >
1 0.2 а>
* 0 -0.2 -0.4
- HJD 2456375.2641 -
2456410^3811 ^^
2456411.4136 ".....4. \
2456412.4063 '" \ /
2456437.3815 ""OX //"""" ...
/ f coadded of 5
/' -
3932.5 3933 3933.5 3934
Wavelength, А
3934.5
3935
Рис. 6. Узкая область вокруг линии Са II на 3933А в спектре HD 129174. Каждый отдельный спектр (толстая линия), а также суммарный спектр внизу, наложены на синтетические (тонкая линия). Еще одна абсорбционная деталь ясно видна в левом крыле наблюдаемой линии.
спектре, но у последнего гораздо хуже отношение S/N.)
3.3.3. Спектр NES
В спектре NES профили линий лучше всего воспроизводятся при значении проекции скорости вращения v sin i = 14 кмс-1. В этой области недостатки и/или избытки содержания отличаются
в сравнении с областью Ca II 3933 A. Детали линий марганца, например, лучше всего воспроизводятся при увеличении модельного содержания на 300. Тем не менее, из 160 наблюдаемых линий марганца в этой области, от сильных до очень слабых, только около 60% совпадают по интенсивности с модельными профилями, тогда как остальные линии оказываются сильнее или слабее в равных долях. Такие расхождения могут быть вызваны неточностями в силе осцилляторов, однако, данные
нашего списка линий не настолько отличаются от данных NIST. Такая же картина наблюдается и при параметрах Teff = 14 000 К, log g = 4.0, что близко к верхней границе температурных оценок других авторов.
Линии в этой области указывают на то, что содержание гелия понижено в 0.1 раз, а магний — уменьшено в 0.3 раз, за исключением линии на 4481 A, которая описывается при двукратном увеличении содержания. Избыток галлия составляет 7.9 х 104, а хрома — 3, тогда как азот и кислород демонстрируют нехватку 0.1 и 0.8 раза соответственно. Избыток фосфора составляет 8, скандия — 35, а циркония — 8.3.
Однако, мы нашли несколько слабых линий с центральными глубинами, достигающими примерно 0.05, которые остаются неидентифицированны-
5500
5505 5510 5515
\Л/ауе!епдШ, А
5520
Рис. 7. Часть спектра ИЭ 129174 около А5510 А. Самая сильная линия на 5510 А является блендой Б II, УII и Мп11 с отношением интенсивностей 5:1:4.
1.02 1
0.98 5 0.96
с аз
с 0.94
0.92 0.9 0.88 0.86
0.84
Мп II Мп II
4839 4840 4841 4842 4843 4844 4845 4846 4847 4848
\Л/ауе1епдШ, А
Рис. 8. Часть спектра ИЭ 129174 около линии ХеП А4844 А. Содержания марганца, ксенона и хрома увеличены в 450, 8000 и 3 раза соответственно.
4602 4603 4604 4605 4606 4607 4608 4609 4610 4611 \/Уауе1епдШ, А
Рис. 9. Часть спектра ИЭ 129174 около линии Хе II А4603 А. Содержания фосфора, ксенона, железа и иттрия увеличены в 7, 4500, 2, 2000 и 3 раза соответственно.
1.05 1
.£> 0.95
от
с
0J 0 9
с
о
> 0.85
(О
(|>
ог 0.8
0.75
0.7
Y II
Fell
Observed Computed
см
CM
cv a> см
Ю
a> X
5288 5290 5292 5294
Wavelength, A
5296
5298
Рис. 10. Часть спектра ИЭ 129174 около линии ХеП А5292 А. Содержания иттрия, железа, марганца и фосфора умножены на 8000, 0.6, 300 и 8 соответственно. Показанные линии марганца нельзя аппроксимировать однозначным значением содержания.
Таблица 4. Список линий, идентифицированных с помощью ШБТ
Длина волны, Идентификация log gf 1
А NIST Эмпирическое Другое
1 2 3 4
4251.1 4251.155 Gall 0.46 0.46 [49]
4254.0 4254.073 Gall -0.20
4255.7 4255.720 Gall 1.38 0.68 [49]
4263.1 4263.136 Gall -0.80
4603. 4603.03 ХеП -0.02 список линий §§Ю800
4844. 4844.33 ХеП 0.51 список линий §10800
5219.7 5219.658 Gall -0.20
5292. 5292.22 ХеП 0.38 список линий §Ю800
5338.2 5338.240 Ga II 0.00
5339.3 5339.33 ХеП 0.00
5360.4 5360.402 Ga II 0.20
5363.6 5363.585 Ga II -0.10
5416.3 5416.318 Gall 0.20
5419.2 5419.15 Хе II 0.10
ми. На Рис. 7 мы показываем часть спектра около А5510 А с пятью неотождествленными линиями.
Дворецки и др. [53] обнаружили линии ксенона у ряда СР звезд класса И§Мп. Они идентифицировали линии ХеП на 4603.03, 4844.33 и 5292.22 А и оценили переизбыток ксенона в 3.1—4.8 dex. Пользуясь этими данными, мы получили оценку содержания ксенона 4.5—8 х 103. Помимо трех данных линий мы также идентифицировали линии ХеП на
длинах волн 5339.33 и 5419.15 А; они приводятся в №БТ, но без величины сил осцилляторов. С другой стороны, линий на 4330.52 и 4464.19 А с интенсив-ностями, сравнимыми со значениями из №БТ, мы не обнаружили. На Рис. 8, 9 и 10 показаны спектры около линий Хе II на 4884.33 А, 4603.03 А и 5292 А, содержания которых подбирались по отдельности.
По данным №БТ мы также отождествили линии галлия, которые не были включены в наш список
линий. Эти линии Ga II и Xe II перечислены в Таблице 4, оценки значений log gf приведены для них в колонке (3).
Мы показали, что в коротковолновом крыле линии кальция Ca II 3933 A видна абсорбционная деталь, которая добавляет примерно 0.07 относительной интенсивности (Рис. 6). А это значит, что значение v sin i переоценено.
4. ЗАКЛЮЧЕНИЕ
Мы обнаружили свидетельства присутствия вторичного компонента в спектре HD 31592, известной спектрально двойной, и классифицировали объект как двойную систему типа SB2. Главный компонент — звезда B9.5V с v sin i = 50 км с-1, тогда как вторичный компонент — звезда A0V с v sin i = 170 км с-1. Этим обусловлено разногласие между значениями приведенными Палмером и др. [1], которые соответствуют более быстро вращающемуся вторичному компоненту, и Вольф-фом и Престоном [2], которые соответствуют более глубокой и узкой линии магния.
HD9531 и HD 129174 схожи в нескольких аспектах. Обе они — СР звезды, первая является переменной класса a2 CVn, вторая классифицирована как звезда класса Mn. Оценки эффективных температур, полученных разными методами, у обоих объектов имеют очень широкий диапазон. Также у всех из них обнаружено магнитное поле. Первая звезда демонстрирует фотометрическую переменность и изменения в профиле линии Ca II 3933 A. Вторая имеет неглубокую и широкую абсорбционную деталь в коротковолновом крыле линии Ca II 3933 A которая, по видимому, движется. Мы не исключаем вероятность, что HD 129174 является SB звездой схожей с HD 9531.
БЛАГОДАРНОСТИ
В данной работе использовалась базы SIMBAD Центра Астрономический Данных, Страсбург, Франция, NIST [21], а также база VALD для атомных данных Венского университета. Работа основана на спектральных наблюдениях на 6-м телескопе САО РАН и на 2-м телескопе Национальной Астрономической Обсерватории, Болгария, а также на данных публичного архива ELODIE. ИР, ДК и ЕС благодарны РНФ за частичную поддержку (грант РНФ 14-50-00043). ИИ и ИС благодарят за частичную поддержку Болгарский ННФ, гранты DN-08/1-2016 и DN-18/13-12.12.2017. ЮЗ благодарит Р. Комжика за поддержку программного обеспечения. Астрономический институт Словацкой академии наук предоставил компьютерное и инструментальное оснащение.
СПИСОК ЛИТЕРАТУРЫ
1. D. R. Palmer, E. N. Walker, D. H. P. Jones, and R. E. Wallis, Royal. Obs. Bull. 135, 385(1968).
2. S. C. Wolff and G. W. Preston, Astrophys. J. Suppl. 37,371 (1978).
3. J. Zverko, I. Iliev, I. Romanyuk, et al., Astrophysical Bulletin 68, 57(2013).
4. J. Zverko, J. Ziznovsky, I. Iliev, et al., Astrophysical Bulletin 66,325(2011).
5. J. Zverko, I. Iliev, I. Romanyuk, et al., Astrophysical Bulletin 67,57 (2012).
6. J. Zverko, I. Romanyuk, I. Iliev, ey al., Astrophysical Bulletin 68,442(2013).
7. J. Zverko, I. Romanyuk, I. Iliev, et al., Astrophysical Bulletin 71, 199(2016).
8. J. Zverko, I. Romanyuk, I. Iliev, et al., Astrophysical Bulletin 72, 17(2017).
9. J. Moultaka, S. A. Ilovaisky, P. Prugniel, and C. Soubiran, Publ. Astron. Soc. Pacific 116, 693 (2004).
10. V. Panchuk, V. Klochkova, M. Yushkin, and I. Najdenov, Journ. Opt. Technology 76, 87 (2009).
11. N. E. Piskunov and J. A. Valenti, Astron. and Astrophys. 385, 1095(2002).
12. P. Skoda, ASPConf. Ser. 101, 187(1996).
13. J. Zverko, J. Z iz ü ovsky, Z. Mikula s ek, and I.Kh. Iliev, Contrib. Astron. Obs. Skalnate Pleso 37, 49 (2007).
14. I. Hubeny, T. Lanz, and C. S. Jeffery, Newslett. Analysis Astron. Spectra, No. 20 (1994).
15. J. Krticka, PhD Thesis, (Masaryk University, Brno,1998) [in Czech language].
16. F. Castelli and R. L. Kurucz, ASP Conf. Ser. 288, A20 (2003).
17. N. E. Piskunov, F. Kupka, T. A. Ryabchikova, et al., Astron. and Astrophys. Suppl. 112,525(1995).
18. T. A. Ryabchikova, N. E. Piskunov, F. Kupka, and W. W. Weiss, Baltic Astronomy 6, 244 (1997).
19. F. Kupka, N. E. Piskunov, T. A. Ryabchikova, et al., Astron. and Astrophys. Suppl. 138, 119 (1999).
20. F. Kupka, T. A. Ryabchikova, N. E. Piskunov, et al., Baltic Astronomy 9, 590 (2000).
21. A. Kramida, Yu. Ralchenko, J. Reader, and NIST ASD Team (2018). NIST Atomic Spectra Database (ver. 5.5.2), [Online. Available: http://pgysics.nist.gov/asd [2018, February 17]. National Institute of Standards and Technology, Gaithersburg, MD.
22. W. L. Wiese and J. R. Fuhr, J. Phys. Chem. Ref. Data 38, 565 (2009).
23. T. Moon and M. M. Dworetsky, Monthly Notices Royal Astron. Soc. 217,782(1985).
24. B. Smalley and M. M. Dworetsky, Astron. and Astrophys. 293,446(1995).
25. N. Grevesse, M. Asplund, and A. J. Sauval, Space Sci. Rev., 130, 105(2007).
26. R. Komzik, private communication.
27. D. Hoffleit and J. H. Warren,Jr., The Bright Star Catalogue, 5th rev. ed. (Astron. Data Center, 1991).
28. F. Royer, J. Zorec, and A. E. Gomez, Astron. and Astrophys. 463, 671 (2007)
29. H. A. Abt, Astrophys. J. 19, 387 (1970).
30. S. C. Wolff, Astrophys. J. 222, 556 (1978).
31. L. Hric, Astrophys. and Space Sci. 169,177(1990).
32. G. Hildebrandt, Astronomische Nachrichten 296,277 (1975).
33. I. McDonald, A. A. Zijlstra, and M. L. Boyer, Monthly Notices Royal Astron. Soc. 427, 343 (2012).
34. P. North, Astron. and Astrophys. 334, 181 (1998).
35. K. Panov and W. Schoneich, Astronomische Nachrichten 297, 177(1976).
36. J. Babel, Astron. and Astrophys. 283, 189 (1994).
37. T. Ryabchikova, O. Kochukhov, and S. Bagnulo, Astron. and Astrophys. 480, 811 (2008).
38. J. Dommanget and O. Nys, Com. de l'Observ. Royal de Belgique 115, 1 (1994).
39. P.C. Schmidtke, J.R. Africano, and R. Quigley, Astron. J. 97,909(1989)
40. H. A. Abt and N. I. Morrell, Astrophys. J. Suppl. 99, 135(1995).
41. D. L. Crawford, J. V. Barnes, J. C. Golson, and D. P. Hube, Astron. J. 78, 738 (1973).
42. J. R. W. Heintze, IAU Symp. 54, 231 (1973).
43. A. N. Cox, Allen's Astrophysical Quantities, 4th ed. (Los Alamos, 1999), p. 381.
44. H. A. Abt and O. Cardona, Astrophys. J. 276, 266 (1984).
45. M. Zboril, Contrib. Astro. Obs. Skalnate Pleso, 16, 47(1987).
46. M. Takada, Publ. Astron. Soc. Japan 31, 473 (1979). 47 . C. R . Cowley and G . C . L. Aikman, Astrophys . J . 242,
684 (1980).
48. S. J. Adelman and D. M. Pyper, Astron. J. 84, 1603 (1979).
49. K. M. Strom, Astron. and Astrophys. 2, 182 (1969).
50. D. Mihalas and J. L. Henshaw, Astron. J. 144, 25 (1966).
51. J. Wm. Montgomery and L. H. Aller, Proc. Nat. Academy Sci. USA 63, 1039(1969).
52. L. Searle, W. T. Lungershausen, and E. L. W. Sargent, Astrophys. J. 145, 141 (1966).
53. M. M. Dworetsky, J. L. Persaud, and K. Patel, Monthly Notices Royal Astron. Soc. 385, 1523(2008).
Stars with Discrepant v sin i as Derived from the Ca II 3933 and MgII 4481 A Lines. VII. HD9531 (SB), HD31592 (SB2), HD 129174 (SB?)
J. Zverko, I. Iliev, 1.1. Romanyuk, I. Stateva, D. O. Kudryavtsev, and E. A. Semenko
In this paper of the series we analyze three stars listed among stars with discrepant v sin i: HD9531 and HD 31592, which also show radial velocity variations inherent to spectroscopic binaries, and HD 129174 which is an Mn-type star with a possible magnetic field. In HD 9531 we confirm the radial velocity derived from the hydrogen lines as well as from the Ca II line at 3933 A as variable. The profile of the calcium line also appears variable, and with the estimated magnetic induction Be = —630 ± 1340 G, this suggests that the abundance of calcium possibly varies over the surface of the star. We identified the lines of the secondary component in the spectrum of HD 31592 revealing thus it is an SB2 binary with B9.5V and A0V components. While the primary star rotates with v sin i = 50 kms-1, the secondary star is faster with v sin i = 170 kms-1. We find that only 60% of the Mn lines identified in the spectrum of HD 129174 can be fitted with a unique abundance value, whereas the remaining lines are stronger or fainter. We also identified two Xe II lines at 5339.33 A and 5419.15 A and estimated their log gf.
Keywords: stars: binary — stars: rotation — stars: individual: HD9531, HD31592, HD 129174
