CC BY
28
УДК 524.38-852/36
ХИМИЧЕСКИЙ СОСТАВ АТМОСФЕРЫ И ЭВОЛЮЦИОННЫЙ СТАТУС СПЕКТРАЛЬНОЙ И СПЕКЛ-ИНТЕРФЕРОМЕТРИЧЕСКОЙ ДВОЙНОЙ
СИСТЕМЫ HD 10009
© 2009 Ю. Ю. Балега1,2*, В. В. Леушин1,2**, М. К. Кузнецов2
1Специальная астрофизическая обсерватория, Нижний Архыз, 369167 Россия 2Южный федеральный университет, г. Ростов-на-Дону
Поступила в редакцию 14 октября 2008 г.; принята в печать 20 октября 2009 г.
Система ИЭ 10009 является спектральной и спекл-интерферометрической двойной с практически одинаковыми компонентами солнечного типа. Она исследовалась спекл-интерферометрически на БТА и спектрофотометрически на 2-м телескопе Терскольского филиала ИНАСАН в спектральной области 3500—9600 А. Проведен детальный анализ атмосферы главного компонента (компонент 1) системы, получены содержания некоторых элементов в его атмосфере. Оценены светимость и температура каждого из компонентов системы: Ь1= 2.9 Ь&, Ь2 = 1.0 Ь&; 1=6017 К, 2=5930 К.
Содержание железа равно [Ре]=—0.27±0.05. Полученные результаты позволяют оценить эволюционный статус компонентов системы. Ее возраст составляет 7.9 миллиардов лет, причем компонент 1 уже близок к области красных гигантов, а компонент 2 все еще находится на стадии выгорания водорода на главной последовательности.
Ключевые слова: звезды — переменные и пекулярные
1. ВВЕДЕНИЕ
Звезда HD 10009 (HR 466, HIP 7580, BD—10 343, LTT 885) — двойная система с большим собственным движением (р,а = 0.2494"/год, Us = 0.111"/год)[1] на удалении примерно 40 пк от Солнца. Орбитальный период системы около 29 лет. Ее интегральный спектр соответствует спектральному классу F7V, а спектры отдельных компонент разделить сложно из-за малой разницы скоростей. Первым двойственность HD 10009 обнаружил Койпер в 1934 г., наблюдая ее визуально на 36-дюймовом телескопе. Им же дана оценка разности блеска компонент, Am = 0.8. По имени первооткрывателя система получила наименование Kui 7. Параметры орбитального движения по визуальным наблюдениям вычислили Хайнц [2] и Бэз [3]. Обе орбиты оказались неверными, что связано с трудностями визуальных измерений тесной пары — на протяжении значительной части орбиты расстояние между компонентами остается меньше 0.1", что часто делает визуальные измерения ошибочными. Хайнц [4] в дальнейшем включил в исходные данные результаты спекл-интерферометрических наблюдений HD 10009
E-mail: balega@sao.ru
E-mail: leushin@sci.lebedev.ru
на 4-м телескопе обсерватории Китт Пик. Полученная орбита имела высокий эксцентриситет е = 0.78, но заниженную величину большой полуоси (а = 0.182"), что давало неверный динамический параллакс и суммарную массу системы (пауи = 0.013", 3.6М©).
Комплексное решение орбиты, учитывающее по методу наименьших квадратов как визуальные, так и спектроскопические измерения, было впервые получено Токовининым [5]. Однако из-за малого числа доступных измерений лучевых скоростей, сосредоточенных вблизи периастра, точность элементов орбиты оказалась невысокой. По орбитальным параметрам суммарная масса системы оценена как 2.7бМ© при отношении масс д = 0.82. Из корреляционных кривых скоростей компонент и суммарного показателя цвета В — V = 0.53 Токовининым найдена разность блеска Ат = 1.23, совпадающая с нашей поздней оценкой по данным спекл-интерферометрии на БТА. Впервые была дана и достаточно реальная оценка расстояния до системы порядка 45 пк, что всего на 20 процентов больше современных оценок.
Мартен и Миньяр [6] определили отношение масс в системе по астрометрическим измерениям ШррагсоБ, д = 1.37, то есть вторичный компонент, по их данным, оказался массивнее
Параметр Tokovinin [5] Hartkopf et al. [9] Pourbaix [10] Balega et al. [11]
Параллакс (х 10~3 ") 22.2 25.73 27 27
а (") 0.292 0.318 0.324
е 0.748 0.6816 0.798 0.79
Р (лет) 28.83 23.7 28.8 29
SM/M© 2.7 2.8 2.16±0.32 2.35±0.35
Дт 1.23 1.20
Sp (1) F5V F5V F5
Sp (2) G2V G G3
главного. Такой результат получен ими по явно заниженной оценке разности блеска между компонентами, Am = 0.41 ± 0.15. Более точный анализ наблюдений пары со спутника Hipparcos выполнен Содерйелмом [7], который получил Am = 1.03 и q = 0.93 ± 0.08. При суммарной массе EM = 2.18 ± 0.26M® массы компонентов оказались равны 1.13M® и 1.05M®, а параллакс п = (25.73 ± 0.93) х 10_3". Наиболее достоверная орбита пары опубликована Пурбо [8], который вслед за Токовининым использовал подбор параметров методом наименьших квадратов по имевшимся астрометрическим и спектральным данным: P = 28.8 ± 0.8 лет, е = 0.798 ± 0.007, a = (0.324 ± 0.005)". Если использовать параллакс Hipparcos, пнгр = 0.02615", получим суммарную массу EM = 2.29Mq . Имеются и более современные спекл-измерения пары, однако все они получены на малых телескопах и не ведут к уточнению параметров. В Табл.1 приведены некоторые данные по параметрам системы.
Все орбиты HD 10009, вычисленные по спекл-данным, подтверждают её высокий эксцентриситет (е = 0.8). Изучение таких двойных звезд имеет большое значение для теорий формирования и эволюции кратных систем. Особенно важно получить информацию о химическом составе, эволюционном статусе и физических параметрах компонентов. Однако наблюдение спектров каждого компонента высокоэксцентричной пары возможно только в эпоху прохождения периастра, когда составной спектр системы разделяется. Для HD 10009 следующее такое событие будет только в 2018 году. Нами выполнены спектроскопические исследования звезды в период, когда спектры компонентов полностью перекрывались. Ниже мы приводим результаты определения параметров атмосфер ком-
понент и предлагаем один из возможных сценариев формирования высокоэксцентричной системы.
2. НАЬЛЮДЕНИЯ И ОЬРАЬОТКА
Спектры HD і0009 получены нами 28—29 августа 2005 г. и 25—26 августа 2006 г. на эшель-ном спектрографе в фокусе кудэ 2-м телескопа Международного центра астрономических и медико-биологических исследований на пике Тер-скол [12] с использованием ПЗС-камеры с матрицей 1242 x 1152 элемента фирмы Write (Великобритания). В каждую ночь было получено по три спектра объекта в диапазоне 3500 — 9600 A с разрешением R = 45000 и отношением сигнал/шум в спектрах около 200. Обработка выполнялась с помощью пакета программ DECH20, созданного в Специальной астрофизической обсерватории РАН [13]. Во время наблюдений звезда находилась далеко от периастра, и спектральные линии не раздваивались.
3. ОЬЩИЕ ПАРАМЕТРЫ КОМПОНЕНТ И ХАРАКТЕРИСТИКИ ИХ АТМОСФЕР
3.1. Эффективные температуры
Спектральные линии звезд системы в полученных нами спектрах не разделяются, что затрудняет определение характеристик каждого из компонентов. Из Табл. 1 видно, что параметры объектов должны мало отличаться друг от друга. Поэтому можно считать, что их спектры довольно похожи.
Фотометрические измерения дают возможность оценить среднюю температуру системы, что при некоторых условиях позволяет оценить температуру каждого из компонентов. В Табл. 2 приведены фотометрические индексы и соответствующие им
Таблица 2. Цвета, температуры и [Fe] системы HD 10009
V 6.23 B-V +0.53 [15] U-B 0.00 [15] /3 2.620 [16] Ъ-у 0.331 [16] mi 0.160 [16] С1 0.387 [16]
Те//(при [Fe] = -0.13) Те//(при [Fe] = -0.27) 6050±130 6000±130 6060±140 6090±140 6115±110 6090±110
величины температур, полученные по зависимостям из работы [14]. Оценки температур получены для двух значений [Fe], первое ([Fe]=—0.13) найдено по индексам mi, c1 и b — y [17], второе ([Fe]=—0.27) определено в настоящей работе по результатам анализа спектров методом моделей атмосфер. Средняя эффективная температура системы при нашем значении [Fe] получается равной 6060±130 К.
Спектральные классы компонент системы, полученные разными авторами и приведенные в Табл. 1, показывают, что температуры компонент должны быть близки друг к другу (около 6000 К), причем первый компонент, вероятно, горячее второго.
Используя видимую звездную величину V=6m.23, известное расстояние п = 0".027 (берутся результаты последних работ [10, 11]) и разницу блеска компонентов Am = 1.20, получаем светимости L1 = 2.90 Lq и L2 = 1.00L©. Помещаем компоненты на общую эволюционную изохрону и оцениваем температуры каждого из них: Te// 1 = 6017 K и Teff2 = 5900 K. Одновременно здесь же оцениваются и ускорения силы тяжести на поверхностях компонентов: log^1=4.04, log g2=4.40.
Поскольку полученные температуры попадают в область величин, где ошибки определяются ошибками фотометрических индексов, мы рассчитали модели атмосфер с температурами от 5800 до 6250 К и logg от 4.00 до 4.50. Химический состав для моделей брался с дефицитом некоторых элементов до —0.3 dex (Табл. 3). Расчеты проводились по программе LLmodels [18] с учетом всех линий поглощения. В Табл. 4 приведены параметры наборов моделей для расчета синтетических спектров, которые сравниваются с наблюдаемыми на Рис. 1.
Модели использовались для расчетов профилей водородных линий в спектрах каждого из компонентов по программе SintVA [19]. Затем профили для первого и второго компонентов суммировались и сравнивались с наблюдаемыми (Рис. 1). Эти сопоставления показывают, что наблюдаемые профили водородных линий в крыльях совпадают с теоретическими для моделей с параметрами: Teff=6017 K, log g = 4.04 для компонента 1 и Te//= 5900 K, log g = 4.40 для компонента 2.
Изменение же температуры на ± 200 K приводит к заметным расхождениям между расчетом и наблюдениями.
Учитывая точность оценок по фотометрическим индексам, можно считать, что атмосферы компонент удовлетворительно описываются параметрами: Te//=6017±100 K, log g = 4.04±0.10 для компонента 1 и Te// = 5900±100K, log g = 4.40±0.10 для компонента 2 при химическом составе атмосферы, указанном в Табл. 3.
3.2. Содержание железа, титана, кремния и микротурбулентная скорость в атмосфере главного компонента HD 10009
Для определения содержания железа в атмосферах компонентов системы мы измерили эквивалентные ширины 50 линий FeI и 43 линии FeII в спектре HD 10009. В связи с тем, что параметры обоих компонентов очень близки, интенсивности их линий, а, следовательно, и эквивалентные ширины должны быть практически одинаковыми. В момент получения наших спектров компоненты системы имели равные лучевые скорости, кроме того они имеют практически равные скорости вращения, что позволяет определить эквивалентную ширину в спектре каждого из них, используя равенство Wa1 = WXobs[n/(n + 1) + WX2 l(WXi(n + 1))]-1.
Здесь Wa1 и Wa2 — эквивалентные ширины линий поглощения в спектрах первого и второго компонента, WA0bs — наблюдаемая эквивалентная ширина линии в спектре системы, n = L1/L2 — отношение светимостей. Различия в температурах и радиусах компонент учитывается через теоретически рассчитанные отношения эквивалентных ширин компонент. В нашем случае для Am=1.2 и n = 3. Отношение WA2/WA1 определялось нами из расчетов эквивалентных ширин для моделей атмосфер с параметрами Te//1 = 6017 K, log g1 = 4.04 и Te//2 = 5900 K, log g2 = 4.40 с химическим составом, представленным в Табл. 3. Вычисленные по указанной формуле эквивалентные ширины вместе с параметрами линий [20] и определенным содержанием элемента могут быть получены по адресу: leushin@sci.lebedev.ru. Расчеты проводились для ЛТР с использованием программы KONTUR [21]. По результатам определения
6540 6550 6560 6570
Wavelength, А
6580
4840 4850 4860 4870
Wavelength, А
4880
r
r
Рис. 1. Сопоставление наблюдаемых (точки) профилей линий Иа и Ир синтетическим спектрам, рассчитанным для моделей с параметрами из Табл. 4. Верхний теоретический профиль соответствует первому набору параметров, средний, практически совпадающий с наблюдениями, — второму, нижний — третьему набору параметров из Табл. 4.
Таблица 3. Содержание элементов в атмосферах моделей и Солнца
Элемент С N О Ne Na Mg Al Si S Ca Sc Ті V Cr Mn Fe
log N(El)* 8.43 7.68 8.56 7.75 6.29 7.24 6.13 7.21 6.84 6.02 3.06 4.65 3.96 5.63 5.05 7.25
log N(El)e 8.39 7.78 8.66 7.84 6.17 7.53 6.37 7.51 7.14 6.31 3.05 4.90 4.00 5.64 5.39 7.45
log N (Fe) для различных значений микротурбу-лентной скорости Vturb строились линейные регрессии log N (Fe) = log N (Fe)0 + k х WA, которые позволяют выбрать реальную величину Vturb в атмосфере исследуемой звезды. Соответствующие регрессии показывают, что микротурбулент-
ная скорость в атмосфере первого компонента ИЭ 10009 близка к 2.00±0.10 км/с.
Та же процедура, что и для линий железа, была проведена для линий Тії и ТІІІ.
Средние содержания железа и титана, полученные по линиям Реї, Pe.II, Тії и ТІІІ для атмосферы первого компонента ИЭ 10009, равны:
Таблица 4. Параметры моделей атмосфер компонентов HD 1OOO9
Компонент Teff Л logs' Li/L2 Vt, км/ с V sin i, км/с
1 5850 3.95 2.9 2.0 10.0
2 5800 4.00 2.0 10.0
1 6017 4.04 2.9 2.0 10.0
2 5900 4.40 2.0 10.0
1 6250 4.18 2.9 2.0 10.0
2 5850 4.00 2.0 10.0
log N(FeI) = 7.18 ± 0.02,
log N (Fell) = 7.22 ± 0.02,
log N (Til) = 4.70 ± 0.04,
log N (Till) = 4.64 ± 0.04.
Среднеквадратичные ошибки среднего содержания характеризует внутреннюю точность полученных величин и показывают для каждого элемента совпадение оценок содержаний по разным ионам в пределах ошибок.
Для определения содержания кремния мы измерили 18 линий в спектре HD10009, и те же линии были измерены для солнечного спектра. Наши эквивалентные ширины линий Солнца практически совпадают с опубликованными ранее [22]. При анализе линий кремния были использованы величины сил осцилляторов из работы [23]. По данным из [20] содержание кремния, полученное по разным линиям, как для атмосферы Солнца, так и для HD10009 отличается между собой на ± 1 dex. В то же время с использованием данных из [23] различия не превышают нескольких сотых dex.
Средние значения содержания кремния:
log N(Si)HDi0009 = 7.33 ± 0.01,
log N(Si)© = 7.57 ± 0.01.
Полученное обилие для Солнца совпадает с опубликованными ранее работами [24, 25] и свидетельствует о том, что в рамках использованных моделей, различия наших определений химического состава атмосферы HD 10009 по сравнению с солнечными реальны.
В Табл. 5 приведена разность содержания железа и титана, определенного по линиям нейтральных и ионизованных атомов Ae(El) = log N(ElI) —
для моделей с разными температурами. Ионный баланс, полученный для ионов железа и титана, подтверждает правильность выбора параметров модели атмосферы главного компонента: Teff 1 = 6017 K, log gi = 4.04.
Таблица 5. Изменение Де = log N (ElI) — log N (ElII) с Tef f модели атмосферы
Teff, к Ae(Fe) Дє(Тг)
5850 0.01 0.13
6017 0.04 0.06
6150 0.33 0.27
3.3. О величинах эквивалентных ширин в спектрах компонентов системы
Мы видим, что для главного компонента ИЭ 10009 наблюдается дефицит содержания металлов, что несколько необычно для звезд в окрестности Солнца. Возникает вопрос, является ли этот эффект реальным или он вызван неправильным учетом двойственности. Если фундаментальные параметры компонент определены достаточно уверенно, то оценки содержания элементов могут быть ошибочными из-за неправильного определения эквивалентных ширин линий в спектрах компонентов по наблюдаемому спектру. Мы считали, что эквивалентные ширины линий компонентов почти одинаковы (теоретически рассчитанные значения для всех линий мало
отличаются от единицы) и почти равны наблюдаемым в суммарном спектре. Это справедливо, если суммируется два практически одинаковых спектра. Если же спектры не одинаковы, то наблюдаемые эквивалентные ширины могут существенно отличаться от величин эквивалентных ширин линий в спектрах каждого из компонентов.
Низкое содержание элементов в таком случае 1оможеГЩ)ть обусловлено ослаблением линий в суммарном спектре. В том случае, если центры линий не смещены относительно друг друга из-за относительных скоростей компонентов и ушире-ния, вызванные вращением, одинаковы, различия эквивалентных ширин могут быть связаны
Элемент log N(£^)Hdiooo9 \og N(El)Q Элемент log N(£^)Hdiooo9 \og N(El)Q
С 8.33±0.10 8.39±0.04 V 3.96±0.15 4.00±0.03
N 7.84±0.10 7.83±0.11 Cr 5.40±0.10 5.65±0.05
О 8.71 ±0.15 8.60±0.05 Mn 5.14±0.07 5.50±0.03
Na 6.03±0.15 6.30±0.03 Fe 7.20±0.05 7.47±0.03
Si 7.33±0.05 7.54±0.02 Co 4.73±0.11 4.91 ±0.03
S 6.99±0.15 7.19±0.04 Ni 5.90±0.09 6.22±0.03
Са 6.27±0.10 6.34±0.03 Ba 1.79±0.15 2.18±0.03
Ti 4.63±0.05 4.92±0.03
только с различием глубин линий и определяются формулой из раздела 3.2. В предельном случае отсутствия измеримых линий в спектре второго компонента (что может осуществиться при очень быстром вращении компонента и расширении линий до полного исчезновения) линии в суммарном спектре становятся в ^ХоЬв/^А1 = и/(и + 1) раз меньше линий первого компонента. По оценкам разных авторов разность блеска компонентов Ат = 2.51с^(£1 /Ь2) варьируется от 0.75т до 1.2т. Это дает нам отношение и = Ь1/Ь2 : 2 ^ и ^ 3. Таким образом, в случае Ат = 0.75 и полном отсутствии линий поглощения в спектре второго компонента, эквивалентные ширины главного компонента будут в 1.5 раза больше наблюдаемых, при этом содержания элементов будут близки к солнечным. Однако из предыдущих наблюдений, проводившихся ближе к периастру, когда линии компонентов разделяются [5], следует, что высказанное выше предположение об отсутствии в спектре линий второго компонента неверно, и величина п ближе к 3, что и использовалось нами при расчетах. Таким образом, можно заключить, что дефицит металлов в ИЭ 10009 — установленный факт.
3.4. Содержание С, N О и некоторых других
элементов в атмосфере первого компонента НО 10009
Поскольку интенсивности линий С, N и О в спектре ИЭ 10009 очень малы и практически все блендируются линиями других элементов, прежде всего линиями Ре1, определение содержания этих элементов по эквивалентным ширинам линий затруднено и подвержено значительным ошибкам. Поэтому определение содержания в атмосфере первого компонента ИЭ 10009 проводилось по
расчетам синтетических спектров в окрестностях линий CI, NI и OI.
Параметры линий выбирались из списка VALD [20]. Механизмы уширения спектральных линий учитывались в соответствии с методикой, описанной в работе [26]. Сопоставления наблюдаемых и вычисленных спектров для участков с линиями CI представлены на Рис. 2, для линий NI — на Рис.3 и для линий OI — на Рис. 4. При расчетах выбранных участков синтетических спектров учитывались все линии, попадающие в этот интервал и имеющие глубины больше
0.01. Среднее содержание углерода, по нашим определениям, составляет:
log N(C)hdiooo9 = 8.33 ± 0.10; азота:
log N(N)HD10009 = 7.84 ± 0.10. Оно оценивалось всего по двум слабым (на уровне шума) линиям и может иметь ошибку существенно больше указанной.
Для кислорода содержание равно:
log N(0)hdiooo9 = 8.71 ± 0.15 .
Содержание кислорода определялось с учетом поправок за отклонения от ЛТР [27]. При определении среднего содержания кислорода использовались эти исправленные значения.
Средние величины содержания элементов по всем определениям приведены в Табл. 6, где сопоставлены содержания элементов, определенные нами относительно солнечных данных [25]. Ошибки, указанные в Табл. 6, соответствуют изменениям содержаний при ошибке моделей в ДT = ±100 и Дк^д = ±0.10.
Wavelength, А
Рис. 2. Сопоставление наблюдаемых (точки) участков спектра с линиями СІ с рассчитанными синтетическими спектрами. Синтетические спектры рассчитаны при меняющихся в два раза величинах содержаний соответствующих элементов.
4. ЭВОЛЮЦИОННЫЙ СТАТУС СИСТЕМЫ ИЭ 10009
Выше упоминалось, что наблюдаемое отношение светимостей компонент (Ат = 1.20) и их общая светимость, оцениваемая по видимому блеску (V = 6т .23) и известному расстоянию (п = 0".027) дает светимости компонентов Ь1 = 2.90Ь& и Ь2 = 1.00£©. Полученные параметры компонентов лучше всего совмещаются с эволюционной изохроной [28] с 1с^ £ = 9.90 (£ — время эволюции в годах) для звезд с уменьшенным содержанием металлов (г = 0.008). По изохронам мы можем уточнить температуры и массы компонентов: Teff 1 = 6017 К, Teff2 = 5900 К и М1 = (1.00 ± 0.15)М© и М2 = (0.9 ± 0.14) М©. Величина содержания металлов г = 0.008 соответствует уменьшенному примерно в 2 раза обилию тяжелых элементов по сравнению с Солнцем. Оценка [Ре]=-0.27±0.05 для ИЭ 10009 показывает, что химический состав этой системы в
пределах ошибок соответствует выбранному г. Сумма масс компонентов в пределах ошибок совпадает с полученной в работах [10, 11] и кажется соответствующей действительности.
Соответствующие изохроны и положения компонент представлены на Рис. 5. Другие комбинации светимостей и температур хуже согласуются друг с другом и наблюдениями и не позволяют поместить оба компонента на одну изохрону. В то же время естественное предположение об общности происхождения компонентов системы требует их одинакового возраста. Поэтому можно считать, что полученные здесь характеристики соответствуют действительности. Общая светимость системы равна 3.9Ь©, (Ь1 = 2.9Ь©, Ь2 = 1.0Ь©), а массы компонентов практически равны солнечной массе. Возраст системы приближается к 8 миллиардам лет и первый компонент находится в стадии окончания выгорания водорода в ядре, а второй на-
Wavelength, A
Рис. 3. То же, что и на Рис. 2, для линий NI.
ходится в начальной стадии эволюции на главной последовательности.
5. ЗАКЛЮЧЕНИЕ
При расстоянии до системы 40 пк обнаруженный дефицит металлов [Ре]=—0.27±0.05 несколько необычен для нашей области Галактики, хотя он не противоречит современным представлениям о химическом составе звезд в окрестностях Солнца [29]. ИШ0009 является системой с большим собственным движением [1]. Исходя из этого, можно предположить, что она образовалась 8 миллиардов лет назад в сферической подсистеме Галактики в области, слабо обогащенной тяжелыми элементами, и за это время достигла окрестностей Солнца. Наиболее точно определенные содержания элементов (С, Si, Т1, Ре и N1) указывают на
слабое увеличения дефицита (по сравнению с Солнечным содержанием) с возрастанием атомного номера элемента. Небольшой избыток кислорода не меняет этой общей тенденции.
Компоненты системы имеют практически солнечные массы, но по возрасту гораздо старше Солнца. В то же время такие динамические характеристики компонентов системы, как микро-турбулентные скорости в атмосферах и скорости вращения, имеют существенно большие значения, чем у Солнца. Наблюдаемая скорость вращения первого компонента сейчас около 10 км/с (близка к этой величине и скорость вращения второго компонента [5]), что при наблюдаемом радиусе первого компонента в 1.6 Rq (радиус определяется через массу 1.0 Mq и log д = 4.04 компонента) соответствует периоду вращения в 8.1 дня. То есть его вращение вокруг своей оси в настоящий момент
6153
6155
6157
6159
6161
6163
Wavelength, А
Рис. 4. То же, что и на Рис. 2, для линий OI.
bg Ге1Т
Рис. 5. Положение компонентов системы HD 10009 на изохронах эволюционных треков.
существенно быстрее, чем движение по орбите. Для звезд солнечных масс нельзя получить по химическому составу свидетельства о перемешивании на главной последовательности и сделать заключение о скорости вращения компонентов системы вокруг своих осей в начальный период эволюции, как это делается для звезд больших масс [30]. Разумно предположить, что скорость вращения компонентов в момент рождения определялась условием устойчивости Vrot = (GM/R)1/2, что для звезды типа Солнца составляет 437 км/с, при этом период вращения равен 0й.116. Наблюдаемые скорости вращения (V sini) F—G-звезд иногда бывают большими (HD 10156, F5, V sin i = 440 км/с; HD 125519, G5, V sin i = 300 км/с [31]), что близко к скорости, определяемой пределом устойчивости. Скорости меньшие, но того же порядка, более многочисленны [31, 32], а с учетом угла наклона часть из них может быть близка к предельной. Причем наибольшие значения скоростей наблюдаются именно при достижении звездой начальной главной последовательности [33, 34].
Если, в соответствии со сделанными нами в предыдущей работе [35] предположениями, считать, что в начальный период система состояла из компонентов со скоростями вращения, обусловленными пределом устойчивости, то их вращательная энергия была равна Erot= 1.58 х 1049 эрг. Если, к тому же, в начальный период орбита системы была близка к круговой с радиусом, равным нынешнему расстоянию в периастре r = а(1 — е) = 370 х 106 км, то ее орбитальный период был около 1000й, что на много больше периода вращения компонентов. И, таким образом, здесь мы можем говорить о приливных выступах [35], которые из-за более быстрого осевого вращения компонентов, по сравнению с их вращением по орбите, приводят к возможности возникновения ускорения движения компонентов по орбите, причем энергия этого ускорения может черпаться из энергии вращения компонент.
Следовательно, первоначальная круговая орбита может перейти в орбиту (а = 1.77х109 км, е = 0.79, T = 10610й), которая наблюдается сейчас. При этом орбитальная энергия системы Eorb = —G х М1М2/2а увеличивается на
1.33 х 1047 эрг, что составляет малую долю от первоначальной вращательной энергии компонент. Здесь мы вновь ссылаемся на работы [36, 37], которые показывают, что предложенный механизм увеличения периода и большой полуоси орбиты может работать в случае двойных звездных систем, аналогично тому, что мы наблюдаем на примере увеличения орбитального периода и полуоси орбиты вращения Луны вокруг Земли.
СПИСОК ЛИТЕРАТУРЫ
1. The Hipparcos and Tycho Catalogues, SP-1200 (ESA, 1997).
2. W. D. Heintz, Astronom. J. 73,565,(1968).
3. P Baize, Astronom. and Astrophys. Suppl. Ser. 60, 333(1985).
4. W. D. Heintz, Astronom. and Astrophys. Suppl. Ser. 72,543(1988).
5. A. Токовинин, Письма в АЖ 19,3(1993).
6. C. Martin and F. Mignard, Astronom. and Astrophys. 330,585(1998).
7. S. Soederhjelm, Astronom. and Astrophys. 341,121 (1999).
8. D. Pourbaix, Astronom. and Astrophys. Suppl. Ser. 131,377(1998).
9. W. I. Hartkopf, B. D. Mason, H. A. McAlister, et al., Astronom. J. 111,370(1996).
10. D. Pourbaix, Astronom. and Astrophys. Suppl. Ser. 145,215(2000).
11. I. I. Balega, Yu. Yu. Balega, K.-H. Hofmann, et al., Astronom. and Astrophys. 385, 87 (2002).
12. Ф. А. Мусаев, Г. А. Галазутдинов, А. В. Сергеев и др., Кинемат. и физика неб. тел., 15, 282 (1999).
13. G. A. Galazutdinov, ftp.sao.ru.
14. A. Alonsj, S. Arribas, and C. Martinez-Roger, Astronom. and Astrophys. 313, 873 (1996).
15. H. L. Johnson, R. I. Mitchell, B. Iriarte, et al. Comm. Lun. Plan. Lab. 4, 99, (1966).
16. P. B. Stetson, Astronom. J. 102,589(1991).
17. B. Nordstrom, M. Mayor, J. Andersen, et al., Astronom. and Astrophys. 418, 989 (2004).
18. D. Shulyak, V. Tsymbal, T. Ryabchikova, et al., Astronom. and Astrophys. 428, 993 (2004).
19. В. В. Цимбал, частное сообщение, (2002).
20. F Kupka, N. E. Piskunov, T. A. Ryabchikova, et al., Astronom. and Astrophys. Suppl. Ser. 138, 119
(1999).
21. В. В. Леушин, Г. П. Топильская, Астрофизика 22, 121 (1985).
22. C. E. Moore, M. G. J. Minnaert, and J. Houtgast, The Solar Spectrum 2935 A to 8770 A, NBS Monograph 61 (Washington, 1966).
23. J. Luck, Astrophys. J. Suppl. 75,759(1991).
24. C. A. Prieto, M. Asplund, R. J. C. Lopez, et al., Astronom. J. 576, 544 (2002).
25. E. Lodders, Astronom. J. 591,1220(2003).
26. Ю. Ю. Балега, В. В. Леушин и Г. Вайгельт, Астрон. ж. 82,247(2005).
27. M. Asplund, N. Grevesse, A. J. Sauval, et al., Astronom. and Astrophys. 417, 751 (2004).
28. L. Girardi, A. Bressan, G. Bertelli, and C. Chiosi, Astronom. and Astrophys. Suppl. Ser. 141, 371
(2000).
29. В. А. Марсаков, Т. А. Боркова, Письма в АЖ 32, 419 (2006).
30. Ю. Ю. Балега, В. В. Леушин, М. К. Кузнецов и др., Астрон. ж. 85, 1 (2006).
31. R. Glebocki, P. Gnacinski, and A. Stawikowski, Acta Astron. 50, 509 (2000).
32. J. R. Stauffer, L. W. Hartmann, C. F Prosser, et al., Astrophys. J. 479,776(1997).
33. J. Bouvier, R. Wichmann, K. Grankin, et al., Astronom. and Astrophys. 3І8, 495 (i997).
34. R. Neuhauser, G. Torres, M. F. Sterzik, and
S. Randich, Astronom. and Astrophys. 32Б, 647 (i997).
35. Yu. Yu. Balega, V. V. Leushin, and M. K. Kuznetsov, Astron. Astrophys. Transact. 2б, i39(2007).
36. O. Struve, B. Lynds, and H. Pillans, Elementary Astronomy (New York, Oxford University Press, 1959).
37. J. Chapront and M. Chapront-Touze, Celest. Mech. 66,31 (1997).
EVOLUTIONARY STATUS AND CHEMICAL COMPOSITION OF THE ATMOSPHERE OF THE SPECTROSCOPIC AND SPECKLE INTERFEROMETRIC BINARY HD 10009
Yu.Yu. Balega, V.V. Leushin, M.K. Kuznetsov
HD 10009 is a spectroscopic and speckle interferometric binary with almost identical solar-type components. It was studied via speckle interferometry using the 6-m telescope of the S pecial Astrophysical Observatory of the Russian Academy of Sciences and had its spectrophotometry performed in the 3500—9600 A wavelength interval with the 2-m telescope of the Terskol Branch of the Institute of Astronomy of the Russian Academy of Sciences. A detailed analysis of the atmosphere of the primary component (component 1) of the binary yielded the abundances of some of the elements. The luminosities and temperatures of the components are found to be L1= 2.9 Le, L2= 1.0 L& and Teffi=6017 K, Teff 2 =5930 K, respectively The iron abundance is [Fe]=—0.27±0.05. Our results make it possible to assess the evolutionary status of the system. The binary is 7.9 Gyr old and the primary component is close to become a red giant, whereas the secondary component is still in the hydrogen-burning stage near the main sequence.
Key words: stars: variable and pequliar
