Software SILICUPS
Jakub Kolář
ÚTFA, MUNI, Brno
Obsah
Základní funkce
Rozplétání světelných křivek - teorie
Rozplétání - praktické úlohy
2 / 44
SImple LIght CUrve Processing System
Pavel Cagaš
https://www.tcmt.org/software.html
Ukazuje časové řady, fázové křivky, poznámky
Fenomenologické modely fázových křivek (Mikulášek 2015)
Rozplétání světelných křivek
Okamžiky minim, O − C diagramy (pouze vykreslení)
Export fázových křivek, modely, rezidua
3 / 44
Nové hvězdné pole
4 / 44
Přidání objektu
5 / 44
Efemeriy, uložení pole
6 / 44
Přidání dat
formát čas mag chyba
desetinná tečka, bez dalších
věcí
.txt and .dat: obojí funkční
7 / 44
Časová série
8 / 44
Fázová křivka
9 / 44
Fitování fázové křivky
10 / 44
Parametry fitu
Amplitudy a délky zákrytů
(nemusí nutně představovat
reálné hodnoty, jen parametr ve
funkci)
Tvar zákrytů
Siny a cosiny
11 / 44
Parametry fitu
12 / 44
Je vhodné vložit přibližné počáteční hodnoty
13 / 44
Modelování
14 / 44
Dostaneme parametry fitu
15 / 44
Uložení parametrů
16 / 44
Určení času minima
17 / 44
Určení času minima
18 / 44
Určení času minima
19 / 44
Přidání do seznamu minim
20 / 44
Minima
21 / 44
Minima
22 / 44
Seznam minim
23 / 44
Vykreslení O − C diagramu
24 / 44
Rozplétání světelných křivek
Můžeme mít dva nebo více zdrojů změn
Dvojzákrytová soustava
Obecný postup je stejný pro různé typy signálů
Začátek: oba signály dohromady, data mohou vypadat divně
Iterativní postup, pár kroků obvykle stačí
25 / 44
Rozplétání
Máme dohromady páry A + B
Nemůžeme modelovat, zavádějící výpočty
Vezmeme si jednu z těch dvojhvězd, většinou pár A (v principu
na tom nezáleží, možná vzít tu, která se bude lépe modelovat, je
výraznější)
Zkusíme ručně stanovit parametry, je to první krok, nebude to
dokonalé
Měli bychom v tuto chvíli mít prvotní model páru A
26 / 44
Rozplétání
Máme model páru A
Máme také naše data, A + B
Odečteme od dat model páru A
Měl by zůstat pár B
27 / 44
Rozplétání
Teď máme samostatný pár B
Spočítáme jeho model
Návrat k původním datům s oběma páry
Odečtení páru B, zůstává pár A
28 / 44
Rozplétání
Dostali jsme oddělená data
Nový model páru A (předtím byl jen manuální odhad bez
výpočtu)
Návrat k původním datům, odečteme nový model A, dostáváme
nové B
Opakováním se dostaneme k dobrému výsledku, většinou stačí
pár iterací, záleží na světelných křivkách
29 / 44
Rozplétání v SILICUPSu
Čtyři nové objekty (převážně
kvůli praktickým důvodům)
Obě dvojhvězdy, perioda
podle A
Obě dvojhvězdy, perioda
podle B
Jen A
Jen B
30 / 44
S1 Leo (TIC 414969157)
Pár A: 4.6305080, 2458873.9590
Pár B: 6.9302830, 2458880.0734
Poměr period 1.4967
Dvě oddělené dvojhvězdy
Vložení dat s oběma páry
31 / 44
Oba páry, fázováno podle A
32 / 44
Oba páry, fázováno podle B
33 / 44
Oba páry s periodou A
Nemůžeme dát fitovat
kvůli kontaminaci
párem B
Dalo by nám to divné
tvary (křivka je
posunuta pro lepší
pohled)
Pro kontaktní
dvojhvězdy by byly
deformace mnohem
větší
34 / 44
Ruční nastavení parametrů
Modelová křivka by měla být trochu nad střední hodnotou hvězdné velikosti
35 / 44
Uložení modelu, odečtení od dat
36 / 44
Který pár nyní máme? Pár B
37 / 44
Modelování páru B
Snažte se nastavit nějaké vstupní hodnoty, někdy je potřeba po výpočtu zkontrolovat amplitudy
38 / 44
Uložení fitu, jít zpátky k oběma párům podle B
a načíst tento fit
39 / 44
Změnit vertikální posun a odečíst pár B z dat
Jiné parametry se nemohou měnit
40 / 44
Po odečtu máme A
41 / 44
Modelování A, uložení parametrů
42 / 44
Konec první iterace
Máme rozpletené páry
Zatím jen první krok
Máme nový model páru A, můžeme ho nahrát do společných dat
a odečíst
Potom získáme nový pár B, snažíme se ho nafitovat a dále použít
Opakování procesu, dokud nejsme spokojeni
43 / 44
Můžete vyzkoušet ASASSN-V J000432.60+632605.0
Pár A: 1.0979870, 2458772.295, pár B: 0.994421, 2458772.360
44 / 44