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

kirx · 07 октября 2005, 20:19:43

Всем привет! :D
Пресказываю то, что писал 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 · 07 октября 2005, 20:52:45

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

kirx · 09 октября 2005, 11:49:40

Посчитал моменты главных минимумов на Октябрь. Время Московское летнее (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 · 09 октября 2005, 17:30:56

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

ctac · 10 октября 2005, 11:18:00

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

То, что в каждых сутках есть 25-тый час, это знает каждый аспирант ;D

kirx · 16 октября 2005, 15:32:41

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

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

Starcat · 16 октября 2005, 15:50:19

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

ctac · 16 октября 2005, 18:59:20

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

kirx · 16 октября 2005, 19:15:57

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

Starcat · 16 октября 2005, 21:20:11

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

ctac · 16 октября 2005, 22:20:05

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

kirx · 17 октября 2005, 01:22:12

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

kirx · 17 октября 2005, 13:58:46

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

ctac · 17 октября 2005, 15:59:40

Да, блин...

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

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

kirx · 19 октября 2005, 11:36:48

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

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

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

kirx · 27 октября 2005, 22:25:45

По наиболее правдоподобным (на данный момент) элементам:
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 · 28 октября 2005, 08:39:28

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

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

ctac · 28 октября 2005, 12:59:30

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

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

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

kirx · 13 ноября 2005, 22:10:30

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

Красные точки - кривая блеска одной из звёздочек в огромном 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 · 14 ноября 2005, 09:38:01

Пора ставить на поток работу с твоей прогой Кирилл ;)

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

kirx · 14 ноября 2005, 15:01:05

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

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

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

maxim · 14 ноября 2005, 15:38:21

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

ctac · 14 ноября 2005, 19:05:02

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

kirx · 14 ноября 2005, 20:11:16

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

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

Olegus · 15 ноября 2005, 15:18:07

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

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

ctac · 15 ноября 2005, 17:04:14

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

Olegus · 15 ноября 2005, 20:26:00

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

kirx · 20 ноября 2005, 19:52:55

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

ctac · 14 января 2006, 18:53:38

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

kirx · 10 ноября 2006, 00:17:16

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

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 · 10 ноября 2006, 00:29:53

YPA! :D

ctac · 08 января 2007, 06:02:36

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

kirx · 12 февраля 2007, 15:50:01

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 · 25 февраля 2007, 14:44:48

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

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

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

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

Балануца Павел · 25 февраля 2007, 15:19:55

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

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

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

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