Тема: Определение моментов минимумов новой затменной переменной.

[kirx]

Всем привет! 
Пресказываю то, что писал ctac'у.
В позапрошлом году по пластинкам ГАИШевской фототеки была открыта затменка типа Алголя. По пластинкам определили моменты двух минимумов, ещё по одному минимуму видели автоматические обзоры ASAS и NSVS. Всё хорошо, только ASAS'овский минимум не ложился на прямую на O-C диаграмме. Мы с моим научруком Антипиным прежде чем публиковать решили посмотреть. В прошлом году посмотрели её в Крыму (самая правая точка на картинке). Легче не стало - теперь NSVS'овский минимум на прямую не ложится. Эррорбары на картинке правдоподобные. Подозрение на изменение периода(!). А в этом году посмотреть не удалось: период ок 4.4d, за всё время наблюдений на ночное время пришлось только 2 минимума и (конечно) обе ночи была сплошная облачность.
Вырисовывается задача: прописать несколько минимумов этой затменки. Если ничего необычного не увидим - определяем более точно период и публикуем простое сообщение об открытии в IBVS (на подобии этих:http://www.konkoly.hu/cgi-bin/IBVS?5700). Ну а если увидим изменение периода - тогда будет совсем интересно. Для полноты картины интересно поискать вторичный минимум на который есть только слабый намёк.
Звёздочка в Андромеде, сейчас смотреть самое то. Фонарь - V(max)=12.5.
Времени понадобится несколько ночей, нужны серии снимков минимум по 4 часа (лучше больше).

На картинке O-C диаграмма. Красная прямая соответствует наиболее вероятному периоду.

[ctac]

 Думаю, что реализация данного наблюдательного проекта возможна с помощью MEADE 8" и ST-2000XM. Матрица прибудет из САО в начале наступающей недели, так что через неделю можно и начать работу
 Надеюсь на скорую встречу в обсерватории , Кирилл.

[kirx]

Посчитал моменты главных минимумов на Октябрь. Время Московское летнее (MSD).
2453653.345920 10/10/2005 0:18:07
2453657.748490 10/14/2005 9:57:50
2453662.151060 10/18/2005 19:37:32
2453666.553630 10/23/2005 5:17:14
2453670.956200 10/27/2005 14:56:56
Понятно, секундной точностью тут и не пахнет, просто удобно было копировать.
А в эти моменты интересно поискать вторичный минимум (предполагая, что он от главного сдвинут на 0.5 по фазе)
2453655.547205 10/12/2005 5:07:59
2453659.949775 10/16/2005 14:47:41
2453664.352345 10/21/2005 0:27:23
2453668.754915 10/25/2005 10:07:05
2453673.157485 10/29/2005 19:46:47
Вообще, хочу приехать на обсерваторию в первый же удобный момент, посмотреть как там у вас!

[ctac]

 Если не секрет, какова продолжительность главного минимума и глубина падения? Насколько быстро происходит падение?

[ctac]

 Только что вернулись с kirx из обсерватории. Там очень потрясно  ! Поставили новые компы и теперь можно прямо во время наблюдений обрабатывать материал
 Вот на скорую руку, что получили в MaximDL  - большое Спасибо Тимуру Крячко, который предоставил эту чудо программу!!!
 Тут не вся кривая и не со всеми точками, первой части нет - т.к. там сбита настройка времени, потом подрисуем
 
 То, что в каждых сутках есть 25-тый час, это знает каждый аспирант 

 

[kirx]

Обработал наблюдения. Впервые хорошо прописан минимум блеска. Время от второго до третьего контакта составляет 1 час 32 минуты. Поэтому на О-С диаграмме представленной на рисунке всем моментам минимумов можно приписать такую неопределённость (1.536h=0.064d) кроме последних двух, полученных по ПЗС наблюдениям когда прописывалась кривая блеска с хорошим разрешением по времени. ASAS-овский минимум (самая верхняя точка) выбивается за эту ошибку, чему можно придумать правдоподобное объяснение. Проблема в том, что минимумы, наблюдавшиеся в 2004 году в Крыму и сейчас (последние две точки), определённые с высокой точностью тоже между собой не очень стыкуются. Так что как обычно, понятно что пока ничего не понятно.
Дальше ещё веселее. В третью ночь затменная наблюдалась сразу на двух телескопах MEADE 8" и 14". И эти наблюдения плохо стыкуются. Этим ставится под серьёзное сомнение  остоверность обнаружения вторичного минимума во вторую ночь наблюдений: его амплитуда не превышает вероятной систематической ошибки. Очевидно, надо уделять большее внимание  орьбе с росой и качеству используемых плоских полей как наиболее вероятным источникам систематики. С дальнейшим планом исследований определюсь в ближайшие дни после консультаций со специалистами.


Кривая блеска за 3 дня наблюдений. Красные точки - 14" MEADE, зелёные - 8" MEADE.

[Starcat]

А гелиоцентрическая поправка вводится?  Может из-за этого несовпадение...

[ctac]

 Кирилл, а ты учитывал,  что мы симали на STL-6303E с фильтрами и только одну серию, как раз колибровачную без, а на ST-2000XM все снимали без фильтров! К тому же у этих чипов разные спектральные кривые чувствительности

[kirx]

Эх, образованые все стали, приятно поговорить!
2 Starcat, гелиоцентрическая поправка даёт минут 8, ну максимум 16, если смотреть в сторону Солнца. Крымский минимум с последним не сходятся почти на 0.04d = 57 минут(!).
2 ctac Меня беспокоит не то что наблюдения с STL-6303E систематичеси отличаются от ST-2000XM, а то что блеск звезды в максимуме вечером 11го измеренный с STL-6303E отличается на 0.2m от измеренного блеска в том же максимуме на той же матрице 12го. Или 12го мы снимали БЕЗ ФИЛЬТРА?!? Сейчас проверю...

[Starcat]

Хи!   
Доступа ко всем данным не имею, звиняйте нас - необразованых...
Насчет нестыковки по зведной величине - либо действительно с фильтрами и чувствительностями матриц пурга вышла, либо методы обработки разные (если апертурная фотометрия, то  одинаковый размер апертуры для разнофокусных изображений вполне может давать сдвиг на 0.1-0.2 звездных).

[ctac]

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

[kirx]

Да, это было бы слишком просто: 12го мы снимали с тем же фильтром что и предыдущие ночи (по крайней мере так следует из названий файла и записей в шапках).
Возможно действительно дело в фокусе, хотя я постарался взять апертуру с запасом 12pix, чтобы влезали и чуть поехавшие звёзды (типа сигнал сильный, шумов много нагрести не должны).
2 Starcat, ну 0.04 следует из картинки: расстояние по y между двумя последними точками, но картинку я сделал халтурную, без подписей и вообще не выразительную. Каюсь!

[kirx]

Вы будете смеяться... Долго, жестоко и совершенно справедливо.
Неделю назад ASAS открыл эту затменку по своим данным. (не прошло и 3х лет)
http://archive.princeton.edu/~asas/ там она называется 003933+2730.5
Публикацию сделаем, но это уже будет не открытие а исследование...

[ctac]

 Да, блин...

 Как всегда опаздали совсем чуть-чуть...

 Будем трудиться, может еще что найдем

[kirx]

Построил новую O-C. Здоровые неопределённости соответствуют длительности фазы полного затмения: когда есть только одна точка в минимуме - мы не можем точно сказать где середина минимума. Такой проблемы нет если мы прописываем всё затмение. Таких точных точек у нас две и они не ложатся на период, определённый по совокупности всех данных. Расхождение составляет 52 минуты. Варианта два: либо мы чётко видим изменение периода, либо на главном компьютере Ка-Дара часы идут с ошибкой ровно в час. Кто там будет в ближайшее время plz. проверьте! (невероятно, но проверить надо).

По оси Х отложена эпоха (Юлианская дата делёная на период), по Y - разница в сутках между наблюдавшимся и вычисленным моментами минимумов.


Вообще было бы здорово посмотреть ещё хоть один минимум. Типа один случай - это случай, 2 - это совпадение, а вот 3 - уже закономерность. Пока у нас вроде 2 совершенно достоверных случая.

[kirx]

По наиболее правдоподобным (на данный момент) элементам:
JDmin=2453653.40414+4.4027699*E
Вычислил моменты главных минимумов на ближайший месяц:
2453670.956200 10/27/2005 14:56:56
2453675.418048 11/01/2005 1:01:59
2453679.820830 11/05/2005 10:42:00
2453684.223612 11/09/2005 20:22:00
2453688.626394 11/14/2005 6:02:00
2453693.029175 11/18/2005 15:42:01
2453697.431957 11/23/2005 1:22:01
2453701.834739 11/27/2005 11:02:01
2453706.237520 12/01/2005 20:42:02

Если будет возможность, а меня не будет в обсерватории, plz смотрите!
Надо стараться прописать нисходящую и(хотябы или) восходящую
ветви затмения для определения момента минимума как точки пересечения прямых,
продолжающих эти ветви.
Координаты звёздочки: 00h39m33s +27 30 33"
обозначение: TYC1744-2329-1
На прилагаемой карте она отмечена как var. Звёзды сравнения и чек хорошо бы
тоже уместить в кадре, что не должно составить труда. При наблюдениях на мы
использовали 90сек. экспозиции.

[kirx]

А вторичный минимум там всётаки есть!
Глубина его в системе unfiltered ccd составляет всего 0.05m, но виден он чётко. Эксперимент показал, что с этим оборудованием можно надеяться поймать изменения блеска вплоть до 0.025m.

на графике зелёные крестики comp.-check, синие кружки - var-comp.+const.

[ctac]

 Очень интресно! Значит можем ловить и экзопланеты

 Рад, что все-таки вторичный максимум словили .

 Постараемся первого ноября словить главный минимум. А то следующий раз может быть только 23 (9-го мы  приедем из Пулково).

[kirx]

Похоже есть шанс отыграться за опубликованную ASAS'ом раньше нас затменку.

Красные точки - кривая блеска одной из звёздочек в огромном vixen'овском поле вокруг затменки. Зелёные - близкая на кадре звезда примерно такого же блеска. Для ASAS'а эта звезда близка к пределу и точности там нет никакой.
Кричать "эврика" пока рано, но проверить надо!
Это звезда GSC 01741-00747 = USNO B1.0 1164-0009591
00 41 34.26 +26 26 14.7
Требуется серия снимков длительностью в несколько часов.
Карты: http://saistud.sai.msu.ru/~kirx/img/vixen2-1.png - полный Vixen'овский кадр,  переменная отмечена маленьким зелёным кружком.
http://saistud.sai.msu.ru/~kirx/img/vixen2-2.png

По самой затменке сейчас сел руками переобрабатывать данные, но момент минимума определённый по (убогой) sextractor'ной фотометрии уже позволяет не сомневаться что ошибки нет, мы чётко видим изменение периода!

[ctac]

 Пора ставить на поток работу с твоей прогой Кирилл
 
 Вопрос, как расположены в поле кадра новая пересенная и звезда сравнения?
 Сколько всего звезд в этом поле Veixen'a было обработанно прогой?
 Новая переменная влезает в поле 14" или нет? а 8"+0.63Reductor+ST-2000XM?

[kirx]

Да, этого кандидата поймали моей прогой. Первая карта - это всё виксеновское поле, она на самом краю, так что там градуса 1.5 - ни во что не влезет,    надо отдельно наблюдать. На чём - пока не очень важно, главное убедиться что это не глюк (самый край поля зрения!).

А на поток ставить пора!!!!

Кстати, новая О-С любимой затменки:

[maxim]

 
  У Vixen на самом краю присуствует слабая кома и виньентирование, эти абберации существенно могут повлиять на результат фотометрирования. Все поле составляет 2°х1°20', а нормальное поле в пределах 1°30'x1°.

[ctac]

 Кирилл, объясни пожалуйста физический смысл "O-C", а то не все читающие могут понимать твою радость/горе по поводу именно такого распределения точек на графике.
 

[kirx]

2 maxim в том то и дело что край поля. Звезда сравнения (которая зелёненькая) находится всего в нескольких минутах от кандидата и такой систематики не показывает, но всякое бывает...

Про физический смысл O-C (точнее математический). Это метод анализа периодических процессов построения графика "Observed-Calculated": зависимости разности между вычесленным и отнаблюдённым моментом чего-нибудь от времени. Например, если по оси X мы откладываем номер периода 1,2,3....10000.... а по оси Y разность между вычисленным и наблюдаемым моментом минимума затменки в сутках и у нас получается прямая, не параллельная оси X, значит наш период неправельный, тангенс угла наклона прямой говорит на сколько (в сутках) нам нужно поправить наш период.
А вот если это вообще не прямая - значит период меняется.

[Olegus]

Меня беспокоит не то что наблюдения с STL-6303E систематичеси отличаются от ST-2000XM, а то что блеск звезды в максимуме вечером 11го измеренный с STL-6303E отличается на 0.2m от измеренного блеска в том же максимуме на той же матрице 12го.

Возник вопрос. Насколько я знаю  ST-2000XM (в отличии от STL-6303E) является ABG камерой со всей вытекающей нелинейностью....
Позволяет ли antiblooming камера делать более менее точную фотометрию?

[ctac]

 Эта нелинейность ABG мешает только для тех объектов у которых более 40 000 отсчетов - ведь эта фьюча начинает действовать когда переполнение ячейки. А так только в чувствительности теряем, но это не много.

[Olegus]

Станислав,
Спасибо за ответ. Теперь понял в чём дело.

[kirx]

Могу добавить, что за всю свою квазипрофессиональную карьеру в Ка-Даре я кажется только второй раз работал с НЕ-ABG камерами! Все профи, которых я видел, как-то предпочитают снимать на антиблюминг. Поэтому обычно на исследуемой звезде не допускают больше 20000 отсчётов в самом ярком пикселе.
Почему так - точно не знаю. Возможно дело в том что, как мы со ctac'ом узнали на собственном опыте , не ABG матрица даже после считывания некоторое время сохраняет "левые" отсчёты в передержанных на предыдущем кадре пикселях.

[ctac]

 Кирилл, как там история со статьей в "Переменные звезды"? Она скоро найдет свой счастливый конец?  У нас будет еще работа по этому объекту, а то я буду перемещать данную тему в раздел "Завершенные работы" Хотя все можно вернуть

[kirx]

Ура, товарищи!
Год с момента проведения основных наблюдений по этой программе в Ка-Даре отмечен выходом долгожданной статьи.


Peremennye Zvezdy (Variable Stars) 26, No. 7, 2006

Received 14 June; accepted 20 October.

Period Changes in the Algol-type Eclipsing Binary System TYC 1744 2329 1
K.V. Sokolovsky1,2, S.V. Antipin1,3, S.A. Korotkiy4

   1. Sternberg Astronomical Institute, 13, University Ave., Moscow 119992, Russia
   2. Astro-Space Center, Lebedev Physical Institute, Russian Academy of Sciences, Profsoyuznaya str., 84/32, Moscow 117997, Russia
   3. Institute of Astronomy, Russian Academy of Sciences, 48, Pyatnitskaya str., Moscow 119017, Russia
   4. Ka-Dar Public Observatory, Barybino, Domodedovo District, Moscow Region, Russia

Abstract: We present the discovery and the times of primary minima for the eclipsing binary TYC 1744 2329 1 = ASAS 003933+2730.5. Our CCD observations compared to times of minima derived from the archival photographic plates and from the NSVS and ASAS-3 data clearly show that the orbital period of the system is variable.


Полный текст лежит по адресу http://astronet.ru/db/varstars/msg/1217495 - в формате HTML, и тут http://images.astronet.ru/pubd/2006/11/09/0001217495/pz-26-007.pdf в PDF(ИМХО, смотрится посимпатичнее).


Есть и перспективы дальнейшего исследования этой звёздочки. Хотелось бы ответить на вопросы: является ли наблюдаемое изменение периода затмений действительно изменением орбитального периода двойной системы вследствие обмена массой между компонентами, или это вращение линии апсид (вызванное присутствием в системе третьего тела, несферичностью звёзд, эффектами ОТО или комбинацией этих факторов). Или наблюдаемое изменение периода - ни что иное как "лайт-тайм эффект", вызванный вращением затменной системы по орбите вокруг третьего тела (эффект Допплера)?
Для ответа на эти вопросы нужно установить следующие наблюдательные факты: является ли периодическим изменение периода (сорри за корявость)? Ясно, что если это так, то большой период составляет десятки лет. Хорошо бы наблюдать, скажем, по одному главному минимуму в год чтобы послеживать за этими изменениями. Кроме того, надо уверенно зарегистрировать и определить момент вторичного минимума (по данным NSVS он чётко виден!) и последить за ним несколько лет. Дело в том, что при движении линии апсид положение вторичного минимума на фазовой кривой блеска меняется (относительно главного), что позволяет отличить апсидальное движение от "лайт-тайм" эффекта.

Так что наблюдательная задача большая. Не на один год. Приглашаем всех любопытных и с телескопами принять посильное участие в её решении!   

[ctac]

 YPA! 

[ctac]

 Кирилл! Кинь эфемериды до февраля. Можем отснять

[kirx]

2 ctac  Сорри, давненько я не появлялся на форуме... 
Посчитал эфимериды на ближайшие 100 дней. Пусть под рукой лежат. 
Время везде всемирное.

Главный минимум:

Код: [Выделить]

Times of primary minima

from Tue Feb 13 00:00:00 2007 to Thu May 24 00:00:00 2007 (UT)

JD=2453675.43190+4.402620*E

JD                              Time

2454146.51224 Thu Feb 15 00:17:37 2007
2454150.91486 Mon Feb 19 09:57:23 2007
2454155.31748 Fri Feb 23 19:37:10 2007
2454159.72010 Wed Feb 28 05:16:56 2007
2454164.12272 Sun Mar 4 14:56:43 2007
2454168.52534 Fri Mar 9 00:36:29 2007
2454172.92796 Tue Mar 13 10:16:15 2007
2454177.33058 Sat Mar 17 19:56:02 2007
2454181.73320 Thu Mar 22 05:35:48 2007
2454186.13582 Mon Mar 26 15:15:34 2007
2454190.53844 Sat Mar 31 00:55:21 2007
2454194.94106 Wed Apr 4 10:35:07 2007
2454199.34368 Sun Apr 8 20:14:53 2007
2454203.74630 Fri Apr 13 05:54:40 2007
2454208.14892 Tue Apr 17 15:34:26 2007
2454212.55154 Sun Apr 22 01:14:13 2007
2454216.95416 Thu Apr 26 10:53:59 2007
2454221.35678 Mon Apr 30 20:33:45 2007
2454225.75940 Sat May 5 06:13:32 2007
2454230.16202 Wed May 9 15:53:18 2007
2454234.56464 Mon May 14 01:33:04 2007
2454238.96726 Fri May 18 11:12:51 2007
2454243.36988 Tue May 22 20:52:37 2007

Предположительные моменты вторичного минимума:

Код: [Выделить]

Times of  secondary minima
from Tue Feb 13 00:00:00 2007 to Thu May 24 00:00:00 2007 (UT)

JD=2453677.63321+4.402620*E

JD                              Time

2454148.71355 Sat Feb 17 05:07:30 2007
2454153.11617 Wed Feb 21 14:47:17 2007
2454157.51879 Mon Feb 26 00:27:03 2007
2454161.92141 Fri Mar 2 10:06:49 2007
2454166.32403 Tue Mar 6 19:46:36 2007
2454170.72665 Sun Mar 11 05:26:22 2007
2454175.12927 Thu Mar 15 15:06:08 2007
2454179.53189 Tue Mar 20 00:45:55 2007
2454183.93451 Sat Mar 24 10:25:41 2007
2454188.33713 Wed Mar 28 20:05:28 2007
2454192.73975 Mon Apr 2 05:45:14 2007
2454197.14237 Fri Apr 6 15:25:00 2007
2454201.54499 Wed Apr 11 01:04:47 2007
2454205.94761 Sun Apr 15 10:44:33 2007
2454210.35023 Thu Apr 19 20:24:19 2007
2454214.75285 Tue Apr 24 06:04:06 2007
2454219.15547 Sat Apr 28 15:43:52 2007
2454223.55809 Thu May 3 01:23:38 2007
2454227.96071 Mon May 7 11:03:25 2007
2454232.36333 Fri May 11 20:43:11 2007
2454236.76595 Wed May 16 06:22:58 2007
2454241.16857 Sun May 20 16:02:44 2007


[kirx]

В ночь с 24 на 24 февраля 2007, в нечеловеческих условиях мы всё-таки сделали это - отнаблюдали ещё один минимум нашей затменки! 
Нечеловечность условий включала в себя: аццкий дубак (на самом деле было всего-то -16, но от этого легче не становилось), ветер (который не только делал вполне умеренный дубак аццким, но ещё и заметно раскачивал монтировку!), высота объекта в момент явления 10 градусов над горизонтом (что плохо не только само по себе, но ещё и заставило вместо 350мм Мида использовать 100мм Виксен, который конечно хорош, но мал), перебитый USB-кабель и ещё много чего было...  Весь день 24го я проспал, а ночью (было уже 25е) сел за обработку. Она тоже проходила нелегко, по началу казалось, что из этих данных нельзя получить вообще никакую фотометрию... Но подобрав методом проб и ошибок оптимальную апертуру и складывая последние (сделанные на совсем уж маленькой высоте) кадры по два или три для измерения я наконец получил кривую блеска:

На рисунке показаны прямые (и их ошибки в 1 сигма) которые я вписал в ветвь спада и в плоскую часть чтобы определить момент наступления полной фазы затмения. Используя наблюдения, полученные в Ка-Даре 30 октября 2005, когда минимум блеска этой звезды был прописан намного лучше я нашёл, что момент наступления полной фазы отстоит от центра затмения на 0.05083d. Зная это, и вводя гелиоцентрическую поправку (которая составила почти ровно 8 минут) удалось определит момент минимума блеска:
HJD2454155.31587 +-0.01 (ИМХО, посчитанная формально ошибка завышена).
На О-С диаграмме можно поставить новую точку:

Из рисунка видно, что изменение периода (по сравнению с полученным совместно по фотографическим и NSVS данным) подтверждается (через все точки никак нельзя провести одну прямую в пределах ошибок). Наблюдавшаяся в последние годы тенденция - сохраняется. Пока нельзя сделать вывод о том, является ли наблюдаемое изменение периода вековым или периодическим. Для этого нужны ещё определения моментов наступления главного минимума с максимально достижимой точностью. Наблюдения хорошо бы проводить не реже чем 1-2 раза в год в течение нескольких (возможно многих) лет. Кроме того, интересно поймать вторичный минимум описать его и последить за его поведением. На существование  вторичного минимума указывают даные NSVS.

В общем, весело понаблюдали! 

[Балануца Павел]

Нечеловечность условий включала в себя: аццкий дубак (на самом деле было всего-то -16, но от этого легче не становилось), ветер (который не только делал вполне умеренный дубак аццким, но ещё и заметно раскачивал монтировку!)

  Понимаю. Сам в этот день ходил наблюдать на поле.

В общем, весело понаблюдали! 

И  главное с пользой!!! И НЛО увидели!!!