Научная статья на тему 'ПЯТЕННАЯ АКТИВНОСТЬ НОВОЙ ПЕРЕМЕННОЙ ТИПА WUMa GSC3599-2569'

ПЯТЕННАЯ АКТИВНОСТЬ НОВОЙ ПЕРЕМЕННОЙ ТИПА WUMa GSC3599-2569 Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

CC BY
271
37
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
звёздные пятна — двойные: затменные — звёзды: индивидуальные: GSC3599-2569 / starspots—binaries: eclipsing—stars: individual: GSC 3599-2569 АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ

Аннотация научной статьи по наукам о Земле и смежным экологическим наукам, автор научной работы — С. Ю. Горда, A. П. Ляпцев, А. М. Соболев

В окрестности молодой яркой массивной переменной V645 Cyg обнаружена новая затменная переменная звезда типа WUMa подкласса A. На интервале наблюдений в четыре года (2009–2013 гг.) зарегистрированы смещения кривой блеска вдоль оси звездных величин с амплитудой 0.m 08 и изменения глубин минимумов на величину до 0.m 055. Получены оценки параметров компонентов этой новой затменной системы. Заподозрено изменение потока излучения системы с периодом 940 суток. Сделано предположение, что возможной причиной таких изменений может быть пятенная активность компонентов. Найдены параметры пятенной структуры для случая максимума пятенной активности.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по наукам о Земле и смежным экологическим наукам , автор научной работы — С. Ю. Горда, A. П. Ляпцев, А. М. Соболев

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Spot Activity of the New WUMa-Type Variable GSC 3599-2569

A new A-subtype WUMa eclipsing binary is found in the vicinity of the young bright massive variable V645 Cyg. Over the four-year long observation interval (2009–2013), the star’s light curve shifted along the magnitude axis with an amplitude of 0.m 08, and the depths of its minima varied by up to 0.m 055. The parameters of the components of this new eclipsing system are estimated, and indications of long-term flux variation with a period of 940 days are found. It is suggested that such variations can be explained by the spot activity of the components. The parameters of the spot structure are determined for the case of maximum spot activity.

Текст научной работы на тему «ПЯТЕННАЯ АКТИВНОСТЬ НОВОЙ ПЕРЕМЕННОЙ ТИПА WUMa GSC3599-2569»

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ, 2015, том 70, № 1, с. 113-121

УДК 524.33/38:520.82

ПЯТЕННАЯ АКТИВНОСТЬ НОВОЙ ПЕРЕМЕННОЙ ТИПА WUMa

GSC 3599-2569

2015 С. Ю. Горда* 1*, A. П. Ляпцев2, А. М. Соболев1

1Коуровская астрономическая обсерватория, Уральский федеральный университет,

Екатеринбург, 620000 Россия

2Уральский федеральный университет, Екатеринбург, 620000 Россия Поступила в редакцию 5 октября 2014 года; принята в печать 24 ноября 2014 года

В окрестности молодой яркой массивной переменной V 645 Cyg обнаружена новая затменная переменная звезда типа WUMa подкласса A. На интервале наблюдений в четыре года (2009—2013 гг.) зарегистрированы смещения кривой блеска вдоль оси звездных величин с амплитудой 0m08 и изменения глубин минимумов на величину до 0 m055. Получены оценки параметров компонентов этой новой затменной системы. Заподозрено изменение потока излучения системы с периодом 940 суток. Сделано предположение, что возможной причиной таких изменений может быть пятенная активность компонентов. Найдены параметры пятенной структуры для случая максимума пятенной активности.

Ключевые слова: звёздные пятна — двойные: затменные — звёзды: индивидуальные: GSC 3599-2569

1. НАБЛЮДЕНИЯ И ПЕРВИЧНАЯ ОБРАБОТКА

Во время проведения четырехлетнего (2009— 2013 гг.) фотометрического мониторинга молодой переменной звезды неправильного типа V 645 Cyg [ 1 ] был получен большой ряд ПЗС-кадров окрестностей этой звезды. Наблюдения проводились на телескопе АЗТ-3 (D = 0.45 м, ANewton = 2.0 м) астрономической обсерватории Уральского федерального университета, оснащенного ПЗС-камерой Apogee Alta-U6 с чипом Kodak KAF-1001E, 1048 х 1048, 24 мкм. На кадрах регистрировалась область неба размером 40' х 40' с центром в районе V 645 Cyg. Первоначальная обработка ПЗС-кадров с целью измерения блеска звезд проводилась в пакете Muniwin.1

Переменность звезды GSC 3599-2569 была обнаружена случайно в процессе поиска звезд, обладающих переменностью неправильного типа и находящихся в ближайшей окрестности переменной V 645 Cyg.

В качестве звезд сравнения и контрольной в непосредственной близости от переменной были выбраны две рядом расположенные звезды, обозначенные нами как C и C1. Данные о звездах приведены в таблице 1. В результате фотометрической обработки было получено 1233 оценки

E-mail: Stanislav.Gorda@urfu.ru

1http://c-munipack.sourceforge.net

блеска переменной звезды в фильтре V и 915 — в R-фильтре.

Предварительный анализ данных показал, что изменение блеска звезды происходит периодически с характерной формой, присущей затменным переменным с коротким значением периода. За весь интервал наблюдений было получено 25 моментов минимумов, приведенных в таблице 2. Кривая блеска, полученная 8 октября 2010 г. (JD 2455478), содержала оба минимума, что позволило получить достаточно уверенную оценку периода P = 0d4 и вычислить значения E для каждого момента минимума, полученного в другие ночи. Точные величины периода, эпохи минимума и погрешностей были получены методом наименьших квадратов путем аппроксимации моментов минимумов линейной зависимостью от E и приведены ниже:

JD©(!min) = 2454979.1586 + 0.4029112 х E.

±0.0004 ±0.0000002

(1)

Кривые блеска переменной звезды в фильтрах V и R с фазами, вычисленными на основе коэффициентов полученного эфемеридного уравнения (1) и с использованием данных всего периода наблюдений приведены, на рис. 1. В верхней части рисунка для сравнения помещены кривые разностей блеска звезды сравнения и контрольной звезды (C — C 1)в фильтре V в тех же фазах.

8 АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 №1 2015

114

ГОРДА и др.

03

я

Е

о

0.4

0.5

0.6

0.7

0.8

0.9

1.0

1.1

1.2

1Ц*|ДМИИГ| И%»»<И1в#

R

-0.3

-0.2

-0.1

0.5

0.6

0.7

0.8

0.9

Е

О

О

О)

га

Е

ОС

<1

-0.50 -0.25 0.00 0.25 0.50 0.75 1.00

Phase

Рис. 1. Полные кривые блеска новой переменной в фильтрах V и R. В верхней части рисунка — кривая разностей блеска звезд сравнения и контрольной в фильтре V.

Таблица 1. Обозначения и координаты звезд

зис GSC Обозначение а(2000) 5(2000) my mR

281-202918 3599-2569 variable 21h39m03?99 +50°09,36,/8 13m26 12m39

281-203195 3599-2386 С, comparison 21h39m19?55 +50°07,39,/9 12m53 llm57

281-203149 - Cl, check 21h39m16?64 +50°07,19,/7 13m02 12m21

Можно видеть, что даже на внезатменных участках на фазах 0.25 и 0.75 вблизи максимума блеска кривые блеска показывают значительные изменения. Глубины обоих минимумов 5т = Ammin — Ammax сравнимы по величине. Разности глубин главного, более глубокого, и вторичного минимумов составляют всего 0m01 и 0m03 в фильтрах в V и R соответственно. Перечисленные выше характеристики кривых блеска, наряду с малой величиной периода, позволяют однозначно отнести данную переменную звезду к типу WUMa (KW по классификации Общего каталога переменных звезд).

Следует также отметить, что глубина главного минимума одинакова в обоих фильтрах 5т = 0m385. На участках кривых блеска, полученных в отдельные ночи на фазах вторичного затмения, хорошо заметно, что вторичный минимум имеет плоское дно. Это указывает на полное затмение в системе и означает, что угол наклона i должен быть значительным, более 80°.

На рис. 1 видно, что разброс точек на кривых

блеска переменной звезды заметно больше величины разброса точек на кривой разностей блеска звезды сравнения и контрольной. Вместе с тем на кривых блеска переменной звезды, полученных в отдельные ночи, наблюдаемый разброс точек оказался значительно меньше. Причиной сильного разброса точек на суммарных кривых блеска являлся сдвиг участков кривых блеска, полученных в разные ночи, друг относительно друга вдоль оси звездных величин. Определенный вклад в разброс точек на суммарных кривых блеска вносит также изменение глубин минимумов.

Амплитуда смещений кривых блеска, полученных в отдельные ночи, и величина изменения глубины минимумов больше в V-фильтре. Амплитуды смещений составляют 0m08 и 0m06 в фильтрах V и R соответственно. Глубины минимумов изменяются в меньшей степени. Амплитуды их изменений составляли 0m055 и 0m045 в фильтрах V и R. При этом разности блеска звезды сравнения и контрольной звезды, полученные в те же ночи, оставались практически неизменными в

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 №1 2015

ПЯТЕННАЯ АКТИВНОСТЬ НОВОЙ ПЕРЕМЕННОЙ ТИПА WUMa GSC3599-2569

115

Таблица 2. Моменты минимумов новой переменной

JD 2450000+ Е О-С JD 2450000+ Е О-С

4979.3582 0.5 -0.0018 5478.1637 1238.5 -0.0004

5035.3654 139.5 0.0007 5478.3643 1239.0 -0.0012

5040.4000 152.0 -0.0011 5533.1606 1375.0 -0.0008

5089.3565 273.5 0.0017 5533.3620 1375.5 -0.0009

5092.3787 281.0 0.0021 5543.2347 1400.0 0.0005

5153.2173 432.0 0.0011 5570.2282 1467.0 -0.0011

5176.1815 489.0 -0.0006 5582.1153 1496.5 0.0001

5223.1225 605.5 0.0012 6205.2168 3043.0 -0.0005

5399.3940 1043.0 -0.0009 6212.2675 3060.5 -0.0007

5413.2958 1077.5 0.0004 6274.1149 3214.0 -0.0002

5414.3027 1080.0 0.0001 6281.1660 3231.5 -0.0001

5441.2971 1147.0 -0.0006 6517.2741 3817.5 0.0021

5442.3058 1149.5 0.0008

пределах разброса точек. Величина среднеквадратичного уклонения разброса точек этих разностей, характеризующая точность наблюдений, не превышала ау = ±0m006. Следует отметить, что значение среднеквадратичного уклонения разброса точек на участках кривых блеска переменной звезды, полученных в отдельные ночи, было сравнимо по величине с этой оценкой.

Отмеченные выше особенности проиллюстрированы на рис. 2, где приведены кривые блеска переменной звезды, построенные для двух периодов наблюдений и показывающие максимально смещенные кривые блеска. Можно видеть, что величины смещений кривых блеска заметно превосходят величины случайного разброса точек на кривых блеска, полученных в отдельные ночи. В тоже время значения разностей блеска звезд сравнения и контрольной (C — C1), полученные в эти же даты, не показывают каких-либо смещений.

Известно, что временные вариации кривых блеска присущи переменным типа W UMa. Наиболее часто отмечаются изменения общего блеска системы, приводящие к смещению участков кривых блеска, полученных в разные ночи, относительно друг друга. Также наблюдаются изменения глубин минимумов и дифференциальные изменения величин обоих максимумов кривой блеска. Наиболее вероятной причиной таких изменений принято считать появление горячих и холодных пятен на

поверхности маломассивных компонентов спектральных классов F—G этих контактных тесных двойных систем.

Для определения величин смещений участков кривых блеска, полученных в разные ночи, с целью уменьшения величины разброса точек в суммарных кривых блеска были построены нормальные (средние) кривые блеска в каждом из фильтров, которые использовались как своеобразные нуль-пункты. Процесс усреднения значений фаз и разностей блеска не представлял особых затруднений, поскольку количество наблюдательных точек было достаточно большим, и они равномерно покрывали кривые блеска в обоих фильтрах. В среднем в одну нормальную точку усреднялось по 15 точек отдельных наблюдений, при этом интервал усреднения по фазе не превышал величины 0.012, что практически не приводило к какому-либо заметному искажению кривых блеска. В результате было получено 68 нормальных точек в фильтре R и 91 — в фильтре V.

Далее вычислялись средние величины смещений кривых блеска, полученных в отдельные ночи, согласно следующему соотношению:

5m = ^2(Am — Amnorm)/N, (2)

где Ami — разность блеска (переменная минус звезда сравнения) отдельной точки на кривой блеска, Am”0™ — разность блеска соответствующей точки нормальной кривой блеска, N — количество измерений за ночь.

Поскольку в основном наблюдался только сдвиг кривых блеска, полученных в отдельные ночи, логичным решением было уменьшить разброс точек на суммарных кривых блеска путем внесения полученных значений сдвигов в разности блеска с обратным знаком (своеобразное сжатие кривых блеска). В результате проведенной процедуры разброс точек на суммарных кривых блеска уменьшился в 2.5—3 раза. При этом форма кривых блеска не претерпела существенных изменений. Исправленные кривые блеска приведены на рис. 3. Именно эти кривые блеска были использованы для определения параметров данной тесной двойной системы.

2. РЕШЕНИЕ КРИВЫХ БЛЕСКА

Решение кривых блеска было проведено c использованием пакета PHOEBE,2 реализующего метод синтеза кривых блеска на основе кода Вилсона—Девиннея [2]. Поиск оптимальных значений искомых параметров в пакете производился методом Нелдера—Мида [3].

2http://phoebe-project.org/1.0/?q=node/78

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 № 1

2015

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

8

116

ГОРДА и др.

03

я

Е

О

О

оз

я

Е

ОС

<1

Рис. 2. Кривые блеска новой переменной и разностей звезд сравнения и контрольной (C — C1), полученные в течение двух периодов: черные кружки — 8 октября 2010 г., светлые кружки — октябрь—декабрь 2012 г.

оз

я

Е

<1

0.5

0.7

0.8

1.0

1.2

V

R

-0.1

0.0

0.1

оз

я

Е

О

6

0.5

0.6 03

я

Е

0.7

<1

0.8

0.9

-0.50 -0.25 0.00 0.25 0.50 0.75 1.00

Phase

Рис. 3. Скорректированные кривые блеска и разности O — C (зависимость O — C в фильтре V смещена на —01“ 1); кружки — наблюдательные точки, сплошные линии — синтетические кривые блеска.

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 №1 2015

ПЯТЕННАЯ АКТИВНОСТЬ НОВОЙ ПЕРЕМЕННОЙ ТИПА WUMa GSC3599-2569

117

Параметры системы находились в классе решений для контактных систем типа W UMa. Для получения в процессе решения реалистичных оценок масс и размеров компонентов, соответствующих маломассивным контактным системам, расстояние между центрами компонентов а принималось равным трем радиусам Солнца, что соответствует среднему значению для систем типа W UMa. Начальные значения поверхностных температур компонентов задавались в предположении спектральных классов F5—G0, наиболее характерных для звезд типа W UMa.

Определение элементов орбиты и физических параметров компонентов затменной системы проводилось с одновременным использованием данных, полученных в фильтрах V и R. Соответствующие полученному решению синтетические кривые блеска, а также разности O — C для обоих фильтров, приведены на рис. 3. Можно видеть, что приближение наблюдательных кривых блеска синтетическими кривыми достаточно хорошее. Зависимости O — C не показывают заметных систематических отклонений.

Полученные значения параметров системы приведены в таблице 3. Здесь i — угол наклона плоскости орбиты системы к картинной плоскости; т\ и r2 — средние радиусы компонентов, выраженные в долях величины a; k — отношение радиусов компонентов; gphot — отношение масс компонентов; Teff д и Teff,2 — эффективные температуры компонентов; ji/j2 — отношение поверхностных яркостей компонентов; f — степень контакта (см. ниже); Q — относительный гравитационный потенциал Копала на поверхности компонентов; Qli и Ql2 — относительные гравитационные потенциалы первой и второй точек Лагранжа; L1 и L2 — относительные светимости компонентов (L1 + L2 = 1).

Приведенные в таблице 3 эффективные температуры и их формальные погрешности соответствуют полученной в процессе решения модели системы. Реальные значения эффективных температур компонентов необходимо определить из спектральных или специально для этой цели проведенных фотометрических наблюдений. Значения температур компонентов приведены в таблице 3 исключительно для того, чтобы показать, что в результате решения кривых блеска более горячим получился первый компонент.

На основе полученных данных можно сделать следующие заключения. Конфигурация системы является контактной, значение степени контакта

f = 100 х (QLi — П)/(Пы — Ql2) = 24.7%.

Вторичный, менее глубокий, минимум синтетической кривой блеска имеет плоское дно (полное

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 № 1

Таблица 3. Параметры относительной орбиты новой переменной звезды

Параметр Значение Параметр Значение

i, deg 85.8 ± 0.5 hlh 1.17

Г1 0.54 ±0.05 /,% 24.7

г2 0.26 ±0.06 п 2.152

к 0.47 p(Ti) 2.181

#phot 0.18 ± 0.01 Р(Т2) 2.066

Teff.l, К 7029 ±10 Li 0.84

ДЯ,2, К 6896 ±10 L2 0.16

затмение) (см. рис. 3). Главный, более горячий компонент (Teff д > Teff 2), значительно массивнее вторичного (<?phot = M2/M1 0.2). Перечисленные

выше значения параметров этой W UMa-системы, особенно малое значение отношения масс компонентов (q < 0.4), позволяют отнести данную W UMa-систему к подклассу А. Для окончательного вывода необходимы спектральные наблюдения.

3. ХАРАКТЕР ИЗМЕНЕНИЯ ОБЩЕГО БЛЕСКА СИСТЕМЫ

Относительно небольшой временной интервал наблюдений данной звезды и отсутствие спектральных данных не позволяют однозначно определить причину изменения общего (незатменного) блеска системы, индикатором которого является смещение кривой блеска по оси звездных величин. Можно перечислить три основных механизма, которые могут быть причиной подобных вариаций:

• пятенная активность на поверхности компонентов;

• наличие в системе третьего тела, обладающего переменностью блеска;

• наличие временного газового диска вокруг системы.

Существование третьего тела, составляющего совместно с затменной тройную систему, пока не может быть подтверждено, т.к. интервал наблюдений этой затменной системы слишком мал для обнаружения явления изменения периода, а спектральные данные, по которым можно было бы обнаружить «лишний» компонент в спектральных линиях, отсутствуют. Тем не менее, гипотеза переменности блеска третьего компонента полностью не может быть отвергнута.

2015

118

ГОРДА и др.

Наличие временного газового диска вокруг тесной двойной также требует дополнительных доказательств его существования и объяснения механизмов его образования в контактной системе.

Таким образом, наиболее вероятным остается предположение о наличии на поверхности компонентов переменной во времени пятенной структуры, которая может быть причиной небольшого изменения потока излучения системы. Такая гипотеза рассматривалась, например, в статье Дюрасевича и др. [4] для объяснения малоамплитудных изменений уровня блеска затменной переменной типа W UMa U Peg.

Интересно также выяснить, имеется ли какая-либо временная закономерность в смещениях кривых блеска. По общему виду полученной зависимости 5m(JD) можно было предположить наличие периодических изменений. Объем имеющихся данных, полученных в течение 39 ночей, распределенных на интервале в четыре года, вполне достаточен для выявления периодической зависимости с характерным временем в 1—4 года.

Поскольку точная форма кривой 5m(JD) неизвестна, было решено представить ее обычной синусоидой (гармоническим колебанием) и, варьируя значения периода, вычислить величину суммы квадратов уклонений S отдельных смещений 5m(JD) от соответствующей синусоиды. За искомое значение периода принять то, которое соответствует минимальному значению суммы квадратов уклонений.

В результате аппроксимации было найдено следующее значение периода, соответствующее минимуму суммы квадратов отклонений: Pmin = 940 ± 130 сут. На рис. 4 представлена кривая 5m(JD) и аппроксимация синусоидой, построенной с этим значением вероятного периода. Видно, что полученная кривая гармонического колебания с найденным значением периода достаточно хорошо согласуется с временным трендом изменения общего блеска двойной системы. Однако из-за пробела в наблюдениях с января 2011 г. по октябрь 2012 г. делать утверждение о наличии периодической зависимости можно только с некоторой степенью вероятности. Данные об изменении общего блеска системы в течение этого интервала времени отсутствуют. К тому же все имеющиеся в наличии данные не покрывают и одного цикла периода.

4. АПРОБАЦИЯ МОДЕЛИ ПЯТЕННОЙ АКТИВНОСТИ

Участки кривых блеска, полученные в ночи увеличенного общего блеска затменной системы (максимальные отрицательные значения 5m), могут соответствовать как минимальному количеству

темных пятен на поверхности компонентов, так и наличию ярких областей. Выбрать один из двух вариантов, исходя из данных наблюдений, не представляется возможным. Иными словами, определить уровень блеска системы, соответствующий отсутствию пятен на поверхности компонентов, мы не можем. В целях определенности было принято первое предположение. К сожалению, построить серию карт распределения пятен на поверхности компонентов для прослеживания эволюции изменения пятенной структуры по нашим наблюдениям не представляется возможным, поскольку при получении наблюдательных данных продолжительность непрерывных наблюдений в отдельные ночи была небольшой. При организации данных наблюдений ставилась задача поиска молодых объектов в окрестности V 645 Cyg, блеск которых мало меняется в течение суток. По этой причине только небольшие участки кривых блеска новой затмен-ной переменной регистрировались в отдельные ночи. Интервалы между датами наблюдений составляли от десятков дней до нескольких месяцев, что представляло определенную трудность в объединении данных нескольких ночей и построения гладкой кривой блеска, соответствующей определенной конфигурации распределения пятен.

Следует отметить, что для построения карты пятенной структуры на поверхности компонентов затменной системы необходимо иметь данные по возможности на всех фотометрических фазах, т.е. иметь практически полную кривую блеска, включающую оба максимума и минимума блеска и — в нашем случае — полученную за достаточно короткий интервал времени.

Как отмечено выше (см. рис. 2), в наблюдениях имеются данные, полученные в течение двух достаточно коротких интервалов времени, на основе которых можно построить две полные кривые блеска, причем заметно смещенные одна относительно другой. Первая из этих кривых блеска была получена в течение одной ночи 2010 г. (JD 2455478), вторая сформирована из данных, полученных в течение шести ночей в интервале порядка двух месяцев 2012 г. (JD 2456212—Jd 2456281). Первая кривая блеска (2010 г.) была получена в период, когда общий блеск системы был ниже зарегистрированного среднего уровня (см. рис. 4). Мы предполагаем, что в этот период пятенная активность компонентов затменной системы могла быть высокой. Данные, соответствующие второй кривой блеска (2012 г.), относятся к периоду, когда зарегистрированный общий блеск затменной системы был максимальным, что соответствует правому максимуму кривой 5m(JD) на рис. 4. В этот период площадь поверхности компонентов, покрытая пятнами, должна быть минимальной, либо темные пятна должны отсутствовать совсем.

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 №1 2015

ПЯТЕННАЯ АКТИВНОСТЬ НОВОЙ ПЕРЕМЕННОЙ ТИПА WUMa GSC3599-2569

119

Рис. 4. Зависимость 3m(JD): кружки — фильтр V, квадраты — фильтр R, штриховая линия — аппроксимация синусоидой.

Исходя из того, что параметры тесной двойной системы на четырехлетнем интервале наших наблюдений должны оставаться практически неизменными, для аппроксимации этих кривых блеска были использованы синтетические кривые, вычисленные на основе параметров (таблица 3), определенных по скорректированной суммарной кривой блеска.

Поскольку, согласно нашему предположению, кривая блеска 2012 г. была получена в период максимума общего блеска системы, предположительно без наличия темных пятен на поверхности компонентов, либо при минимальном их количестве, аппроксимация ее синтетической кривой должна осуществляться простым смещением последней вдоль оси интенсивностей до максимально возможной степени совмещения. Иными словами, при выполнении процедуры аппроксимации варьированию подвергался только один параметр luminosity, ответственный за общий блеск системы. Полученное в результате аппроксимации новое значение параметра luminosity было принято за основной блеск системы.

При аппроксимации кривой блеска 2010 г. (JD 2455478) все параметры, в том числе и новое значение параметра luminosity, фиксировались. Аппроксимационное соответствие достигалось исключительно введением параметров модели пятенной структуры. Наилучший результат был достигнут в модели наличия обширного темного пятна в районе одного из полюсов главного, более горячего, компонента. Это обусловило снижение

общего блеска системы, существенно не повлияв на форму кривой блеска. Небольшая разница в высоте максимумов кривой блеска 2010 г. хорошо моделировалась введением небольшого темного пятна в экваториальной области главного компонента на долготе, соответствующей фазе

0.75. Все вычисления были проведены в пакете PHOEBE. Результаты аппроксимации приведены на рис. 5 и 6, а параметры модели пятенной структуры — в таблице 4. Здесь <р — широта центра пятна, отсчитываемая от полюса звезды, Л — долгота центра пятна, отсчитываемая от линии, соединяющей центры компонентов, против часовой стрелки, если смотреть со стороны полюса, для которого <р = 0, r — радиус пятна, Ts/Tf — отношение эффективных температур пятна и поверхности звезды.

Следует отметить, что полученная в данном случае простая модель пятенной структуры достаточно хорошо аппроксимирует наблюдательные данные, но не является единственно возможной. Так, для

Таблица 4. Параметры пятенной структуры поверхности главного компонента, полученные по наблюдательным данным за 8 октября 2010 г.

№ К, deg A, deg г, deg Ts/Tf

1 3 180 33 0.68

2 80 88 8 0.80

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 № 1

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

2015

120

ГОРДА и др.

Рис. 5. Синтетические (сплошные линии) и наблюдаемые кривые блеска, полученные 8 октября 2010 г. во время значимой пятенной активности (черные кружки) и в период октября—декабря 2012 г., когда пятна, предположительно, отсутствовали (светлые кружки).

Рис. 6. Зависимости O — C, обозначения те же, что и на рис. 5. Во избежание наложения графиков три верхние зависимости O — C приведены со смещениями —01“05, —01“ 15 и —0“20.

объяснения падения блеска, вероятно, наиболее реалистичным будет предположить наличие групп пятен в районах обоих полюсов более горячего компонента или на обоих компонентах сразу. Для моделирования разной высоты максимумов можно ввести горячее пятно, например, на вторичном компоненте и т.п. Изменение глубин минимумов можно моделировать наличием пятен на затмеваемых областях поверхностей компонентов. В любом случае мы получим теоретическую кривую, лежащую внутри коридора разброса точек кри-

вой блеска, полученной из наблюдений, либо в той или иной степени выходящей за его пределы. Наличие только фотометрической информации не позволяет сделать однозначный выбор конкретной модели пятенной структуры. Важно, что в рамках изменяемой во времени пятенной модели можно с достаточной степенью убедительности объяснить небольшие флуктуации общего блеска этой за-тменной системы. Для более детального исследования необходимы дальнейшие наблюдения данной звезды по специально разработанной программе.

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ том 70 №1 2015

ПЯТЕННАЯ АКТИВНОСТЬ НОВОЙ ПЕРЕМЕННОЙ ТИПА WUMa GSC3599-2569

121

5. ЗАКЛЮЧЕНИЕ

Открытая нами новая переменная звезда GSC 3599-2569 по виду кривой блеска и короткому периоду была классифицирована как затменная переменная типа WUMa. На основе значений параметров компонентов и системы в целом, полученных в процессе решения кривой блеска, было установлено, что эта контактная маломассивная затменная система типа W UMa принадлежит к подклассу A. Показано, что обнаруженные на интервале наблюдений в четыре года небольшие флуктуации общего блеска системы, вероятно, вызваны пятенной активностью компонентов и могут иметь периодический характер. Для объяснения причин флуктуаций общего блеска этой затмен-ной системы нельзя исключать из рассмотрения предположение о наличии третьего, переменного, компонента, а также образование временного газового диска вокруг системы. Для однозначного выбора того или иного предположения необходимо получение новых фотометрических и особенно спектральных данных об этой системе.

БЛАГОДАРНОСТИ

Авторы выражают благодарность В. П. Го-ранскому за конструктивные замечания. Работа выполнена при частичной финансовой поддержке Министерства образования и науки РФ (госзадание № 3.1781.2014/K; базовая часть госзадания, регистрационный номер 01201465056).

СПИСОК ЛИТЕРАТУРЫ

1. A. M. Sobolev, S. Yu. Gorda, and O. A. Davidova, IBVS, № 6061 (2013).

2. R. E. Wilson and E. J. Devinney, Astrophys. J. 166, 605(1971).

3. J .A. Nelder and R. Mead, Computer J., 7, 308 (1965).

4. G. Djurasevich, H. Rovithis-Livaniou, P Rovithis, et al. Astron. and Astrophys. 367, 840 (2001).

Spot Activity of the New W UMa-Type Variable GSC 3599-2569

S. Yu. Gorda, A. R. Lyaptsev, and A. M. Sobolev

A new A-subtype WUMa eclipsing binary is found in the vicinity of the young bright massive variable V 645 Cyg. Over the four-year long observation interval (2009—2013), the star’s light curve shifted along the magnitude axis with an amplitude of OP08, and the depths of its minima varied by up to 0m055. The parameters of the components of this new eclipsing system are estimated, and indications of long-term flux variation with a period of 940 days are found. It is suggested that such variations can be explained by the spot activity of the components. The parameters of the spot structure are determined for the case of maximum spot activity.

Keywords: starspots—binaries: eclipsing—stars: individual: GSC 3599-2569

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ

том 70

№ 1

2015

i Надоели баннеры? Вы всегда можете отключить рекламу.