Astronet forum

Heureka!

Vi har tidigare haft uppe RR Lyrae-stjärnan KR Lyr till diskussion här på Astronet. Fasdiagrammet ser kaotiskt ut. Men efter att ha gått igenom de publicerade originalobsarna så testade jag med en helt annan period än den som ges i GCVS. Och då blir plötsligt SVO-fasdiagrammet för stjärnan snyggt och prydligt, se nedan. Den korrekta perioden är förmodligen ca 0.655 dygn.

mvh
Hans Bengtsson

Vilken fin upptäckt, Hans! Jag föreslår att du skickar in en ändring till VSX.

/Gustav

Bra jobbat Hans! Detta får bli ett tidsserieprojekt till hösten!!!

Mvh,
Tomas W

Om vi sätter Hoffmeisters första observerade maximum som utgångsmaximum och ger perioden 0.654998 dygn så får vi följande (O-C)-värden för KR Lyr:

E = 0000 O = 30209.534 C = 30209.534 O-C = 0.000
E = 0388 O = 30463.574 C = 30463.673 O-C = -0.099
E = 0441 O = 30498.397 C = 30498.388 O-C = 0.009
E = 0516 O = 30547.474 C = 30547.513 O-C = -0.039
E = 0519 O = 30549.425 C = 30549.478 O-C = -0.053
E = 0878 O = 30784.582 C = 30784.622 O-C = -0.040
E = 0968 O = 30843.468 C = 30843.572 O-C = -0.104
E = 0974 O = 30847.447 C = 30847.502 O-C = -0.055
E = 1006 O = 30868.444 C = 30868.462 O-C = -0.018
E = 1154 O = 30965.355 C = 30965.402 O-C = -0.047
E = 1513 O = 31200.552 C = 31200.546 O-C = 0.006
E = 1519 O = 31204.474 C = 31204.476 O-C = -0.002
E = 1895 O = 31450.687 C = 31450.755 O-C = -0.068
E = 2055 O = 31555.547 C = 31555.555 O-C = -0.008
E = 2600 O = 31912.543 C = 31912.529 O-C = 0.014
E = 2867 O = 32087.411 C = 32087.413 O-C = -0.002
E = 2902 O = 32110.338 C = 32110.338 O-C = 0.000

Tills vidare tror jag detta är de element vi kan sätta för stjärnan. Fasdiagrammet i SVO blir då enligt nedanstående.

Edit:
Men jag ser ju också av det nya fasdiagrammet att vi kan behöva justera utgångsmaximets tid en smula, så att vi nu hamnar precis i fas 0.0 vid max.

mvh
Hans Bengtsson

Spännande! Vad är E, O & C värden?

Håller med, bra jobbat! Det var en stor skillnad på perioden du har upptäckt jämfört med vad som anges i GCVS (0.3957714).

mvh
Robin

Mycket välmotiverad fråga!

När man ska förutse tidpunkterna (julianskt datum) för en periodisk variabel stjärnas maxima så använder man ljusväxlingselement som brukar ha formen:

Max = M0 + P * E

M0 (M-noll) är tidpunkten för ett utgångsmaximum (initialmaximum).
P är periodlängden i dygn.
E är ett valfritt heltal, helt enkelt maximets "nummer i ordningen" sedan utgångsmaximet M0.

Vi behöver också känna till förkortningarna O och C:
O är tidpunkten (JD) för observerat maximum.
C är tidpunkten för beräknat (calculated) maximum enligt ljusväxlingselementen.

För varje maximum kan man således beräkna ett värde (O-C), som visar ifall det observerade maximet kommer för tidigt eller för sent i förhållande till ljusväxlingselementen. Om man har många bestämningar av (O-C)-värden kan man rita en (O-C)-kurva som visar hur perioden för en variabel förändras under årens lopp. Vi har en verklig specialist på (O-C)-studier i gruppen, nämligen Thomas Karlsson!

I det nu aktuella fallet, KR Lyr, visade det sig nu att upptäckaren Hoffmeister har tagit fel på periodlängden, han har helt enkelt valt fel värde. Det händer att fler än en periodlängd är matematiskt möjlig om man har ett begränsat obs-material. SVO-materialet visade ett kaotiskt fasdiagram (när man gör ett fasdiagram så hänför man alla observationer till en enda ljusväxlingscykel), så antingen måste det ha varit stora mätfel på flera SVO-obsar eller också var perioden en helt annan än den katalogiserade. Det senare visade sig vara fallet, jag hittade en andra möjlig period när jag kollade Hoffmeisters originaldata - och den visade sig ge ett utmärkt fasdiagram!

Som Tomas W antyder, så är det allra bästa om vi kan bestämma ett färskt och fräscht maximum med hjälp av tidsseriemätningar. Då får vi ett utgångsmaximum som är välbestämt, och perioden kan sedan beräknas t.ex. med hjälp av minstakvadratmetoden där vi utnyttjar både SVO-data och Hoffmeisters data.

Ganska kul att hitta något som ingen annan har känt till under alla de decennier som gått sedan stjärnan upptäcktes. Ytterligare en spännande bifångst i projektet med de 50 bortglömda mirorna!

Jag tror förresten knappast att Thomas K kommer att kunna motstå frestelsen att titta närmare på dessa data för KR Lyr ...

mvh
Hans Bengtsson

Fin upptäckt Hans! Var fick du tag i uppgifterna om Hoffmeisters maxima?

Det verkar som perioden han anger är 3/5 (0.6) av den nya, det är ett ganska vanligt fel om man har lite data och glest mellan punkterna i förhållande till perioden att man får andra multipler av perioden att stämma. Vi hade ju ett sånt problem med V2331 Cyg i början där olika perioder var möjliga.

Jag körde in våra obsar i Parenso och fick perioden 0.6550, snarlikt den du fick. Jag hittade också ett maximum till i http://www.konkoly.hu/cgi-bin/IBVSpdf?5918 på 54943.351. Från våra obsar kan man få ett maximum till 55989.63. Med dessa två nya maxima kanske du kan förfina perioden ytterligare.

EDIT: 54943.531 ska det vara ovan.

Intressant! Lyra är ju fotograferingsbar kring midnatt i dagarna. Det kanske skulle gå att få till en eller flera tidsserier innan ljuset tar över. Men det lär väl bli en kamp mot kalendern och molnen...

Mvh,
Tomas W

Det lustiga var hur oerhört snabbt och enkelt det gick att hitta den korrekta perioden. Jag började med att hämta hem Hoffmeisters maxima från GEOS ...

http://dbrr.ast.obs-mip.fr/

... och klistrade in dem i Excel. Därefter valde jag ut tre olika tidsintervall som fanns mellan Hoffmeisters maxima, och lät Excel kolla vilken period man fick om man dividerade intervallen med olika heltal. Satte upp dessa data i tre kolumner bredvid varandra och såg då att det fanns två olika periodlängder som förekom i alla tre kolumnerna: dels Hoffmeisters värde, dels ett värde kring 0.655 dygn. Lade in 0.655 dygn i SVO - och på en sekund förstod jag att jag hamnat rätt, när en urtypisk RRab-kurva uppenbarade sig på skärmen!

Sedan var det bara att kolla så alla Hoffmeisters 17 maxima stämde med den nya perioden, och det gjorde de inom felmarginalerna. Att perioden visade sig ligga så väldigt nära ett till synes rätt approximativt värde som 0.655 dygn var bara en slump, det är fullt möjligt att den kan vara 0.655000 dygn även om man får 0.654998 när man kollar intervallet mellan Hoffmeisters första och sista bestämning.

mvh
Hans Bengtsson

Jag utökade O-C diagrammet med IBVS- och vårt maximum

E=37761 O=54943.531
E=39358 O=55989.62

Jag gjorde O-C diagrammet i Excel och infogade en trendlinje i diagrammet och valde att visa trendlinjens ekvation. Bara man får E rätt får man från koefficienterna i denna ekvation vad perioden och startepoken ska justeras med för att trendlinjen ska bli vågrät och punkterna centrerade kring O-C = 0.

Genom detta fick jag periodelemenent Max = 2430209.4803 + 0.6550161 * E

Jag lägger in dessa värden i SVO så får vi lättare att se hur kommande obsar passar in på kurvan.

Utomordentligt. Precis vad jag hoppades på, proffsigt.

mvh
Hans Bengtsson

Jag har nu meddelat Sebastian Otero, VSX, de nya elementen.

mvh
Hans Bengtsson

Jag lyckades fånga en tidsserie av KR Lyr under natten. Start ca 2330 svensk tid och slut ca 0300. Exponeringar om tre minuter, jag tror att det blev 64 stycken, för att få hyfsat SNR. Jag tycker själv att mätningarna är lite fladdriga, det kan bero på att seeingen inte var helt hundra. Jag ska göra ett nytt försök så snart det går.

Mvh,
Tomas W

Jisses vilken fin mätserie av KR Lyr du nu levererat från natten 23-24 mars! Plus några andra mycket intressanta serier. Suveränt jobbat.

M v h
Hans Bengtsson