Pierwsze fotometryczne koty za płoty SN2023ixf

[WojtasSP320]

Przyszła piękna pogoda a wraz z nią białe noce, które na moich szerokościach zaczęły się dobre kilka dni temu. Jest to zatem dobry czas na przerzucenie się na inne eksperymenty i gałęzie astrohobby.
Szczęśliwie się złożyło, że akurat teraz pewnej gwieździe w galaktyce M101 zachciało się wybuchnąć i w ten oto sposób możemy podziwiać supernową SN2023ixf.
Każdy już porobił swoje pamiątkowe fotografie, ale ja chciałem wypróbować coś jeszcze - chciałem podejść do tematu trochę bardziej naukowo.

Zatem jak każdy poważny obecnie naukowiec zrobiłem szybki risercz na pierwszej stronie googla (nie bądźmy szurami, nie będziemy zaglądać na drugą stronę wyników) i znalazłem ten oto szanowny artykuł 

Problem z nim (jak i z większością tego typu opracowań) był taki, że większość linków do zewnętrznych zasobów już nie działa. Ale ok udało mi się co nieco odkopać a co nieco znaleźć w inny sposób więc lecimy z tematem.
Dodam tylko, że jestem w tej dziedzinie zupełnym amatorem i nie znam się na tym ani trochę (robiłem metodą - odpal tutorial, rób dokładnie tak jak mówią i nie myśl za dużo, więc jeśli pojawiły się jakiekolwiek błędy to proszę o uwagi i krytykę).

Setup:

GSO 150/750mm

ZWO ASI1600MM Pro
ZWO B Filter

Akwizycja:
2023-05-22 23:11:00 00:35:00 (zbierane naprzemiennie L-R-G-B w klatkach co minutę, łącznie 20 minut kanału B )
Gain 0 Offset 21 Temp. -20°C

Parametry dobrałem tak, żeby jasność w kanale L supernowej miała wartość ok 60k (dla kanału B będzie to mniej)

Stack:
APP - kalibracja darkami, flatami i biasem > stack (nic więcej)

Fotometria:
Znalazłem gwiazdę z katalogu TYCHO-2 w pobliżu M101 w Stellarium:


Odnalazłem gdzieś w Internecie (http://simbad.u-strasbg.fr/simbad/sim-id?Ident=TYC+3852-521-1) jej dane fotometryczne (B 11.77):

 

Następnie po ciężkich bojach, żeby odpalić SalsaJ (w odpowiedniej wersji) udało mi się załadować stack kanału B i zaznaczyć obie gwiazdy (referencyjną TYC-3852-521-1 i SN2023ixf):


Trochę dziwnie (moim zdaniem, ale ja się nie znam) wygląda to zaznaczenie, ale było ustawione domyślnie w programie jako auto (FWHM). Dane jakie dostałem dla obu gwiazd to:


Potem podstawienie do wzoru Pogsona:


I wyliczenie jasności w Wolframie (no bo przecież nie na piechotę)


Wynik: 10.27mag (zły logarytm) 11.12mag

Zgadza się? Nie zgadza się? Popełniłem jakiś błąd? Jak dla mnie przynajmniej wygląda legitnie (a czy wszystko jest ok, to nie wiem).

Jako bonus dwa stacki:

2023-05-21 kanał L (110x60s):

 

2023-05-22 LRGB (50x60s + 3x20x60s):


A tu dla porównania z przed 2 miesięcy

 

Edytowane 24 Maja 2023 przez WojtasSP320
poprawka wyniku

[isset]

6 godzin temu, WojtasSP320 napisał(a):

Szczęśliwie się złożyło, że akurat teraz pewnej gwieździe w galaktyce M101 zachciało się wybuchnąć

 

Nie teraz, tylko 21 milionów lat temu 😄

[WojtasSP320]

Teraz, 21 milionów lat temu... nie kłóćmy się o takie drobne szczegóły  

A jeśli już chcesz być tak bardzo dokładny to musisz mi jeszcze napisać, czy stoisz, idziesz pieszo czy jedziesz rowerem
https://en.wikipedia.org/wiki/Rietdijk–Putnam_argument

I klasyk na koniec posta:

[Michal_IAUZ]

Hi,

 

pomiary wyglądają w miarę ok, rzeczywiście zaznaczenia wyglądają trochę dziwne, być może warto ręcznie zwiększyć apertury.

Co najważniejsze wynik nawet się jako tako zgadza, pewne różnice mogą wynikać z tego że znalazłeś pomiary w filtrze B systemu Johnsona, ty pewnie używasz B "kolorowego" który ma trochę inny zakres.

 

Co do oprogramowania to polecam przejście na AstroimageJ, zarówno do redukcji jak i fotometrii. Szczególnie przydatne przy większej ilości zdjęć, kiedy chcesz stworzyć krzywą blasku.
https://www.astro.louisville.edu/software/astroimagej/
https://britastro.org/wp-content/uploads/sites/Tutorial_AIJ_Photometry_v2.06.pdf
https://www.astro.louisville.edu/software/astroimagej/guide/AstroImageJ_User_Guide.pdf

 

Jeżeli chcesz zrobić pomiar dla jednego obrazka to możesz to zrobić nawet zliczając pixele ręcznie np w ds9 - https://sites.google.com/cfa.harvard.edu/saoimageds9

Możesz użyć regionów i zobaczyć sumę z liczeń z danego obszaru.

Łatwo możesz znaleźć też janości innych gwiazd porównania używając tego okienka w simbad

Możesz wybrać promień szukania (liczy od środka obrazka), następnie po kliknięciu logo SIMBAD dostaniesz oznaczenia wszystkich obiektów z bazy.

Edytowane 24 Maja 2023 przez Michal_IAUZ

[WojtasSP320]

1 godzinę temu, Michal_IAUZ napisał(a):

Co najważniejsze wynik nawet się jako tako zgadza, pewne różnice mogą wynikać z tego że znalazłeś pomiary w filtrze B systemu Johnsona, ty pewnie używasz B "kolorowego" który ma trochę inny zakres.

Tak właśnie było, ale jak się nie ma co się lubi...
Zastanawiam się czy dałoby się uwzględnić różnice między filtrem fotometrycznym B a moim zwykłym B znając charakterystyki obu filtrów (mam zmierzony B przez @dobrychemik).
Śliczne dzięki za obszerny post z sugestiami, spróbuję zmierzyć to jeszcze raz proponowanym sposobem.

[Piotrek Guzik]

@WojtasSP320

We wzorze Pogsona jest logarytm dziesiętny, a Ty w Wolframie użyłeś logarytmu naturalnego, dlatego dostałeś błędny wynik. Jeśli wstawisz tam logarytm dziesiętny, to wyjdzie Ci około 11.12 mag, czyli z dokładnością do kilku setnych części wielkości gwiazdowej tyle samo, co wyniki innych pomiarów zgłoszonych w tym czasie do AAVSO. Te kilka setnych magnitudo może wynikać właśnie z różnicy między filtrem "zwykłym" a fotometrycznym. Żeby tę różnicę uwzględnić, musiałbyś znać zarówno charakterystyki filtrów, jak i widmo obiektu, którego jasność chcesz zmierzyć.

[beny]

Fajne opracowanie  A w spajkach z luminancji jestem zakochany  Chyba czas założyć newtona na montaż

[WojtasSP320]

32 minuty temu, Piotrek Guzik napisał(a):

We wzorze Pogsona jest logarytm dziesiętny, a Ty w Wolframie użyłeś logarytmu naturalnego, dlatego dostałeś błędny wynik.

Ha! Właśnie takich błędów bałem się najbardziej. Rzeczywiście w Wolframie logarytm dziesiętny to log10(x) (jestem przyzwyczajony, że naturalny to ln(x) a dziesiętny to log(x)).
I tak jak mówisz po zmianie https://www.wolframalpha.com/input?i=-2.51+*+log10(517662.72%2F285205)%2B11.77 wychdzi 11.12mag czyli podobnie jak u innych.

A co do spajków to wydaje mi się, że trafił mi się wyjątkowo dobry seeing (chociaż warunki z których robiłem były paskudne, środek miasta obok nagrzanego bloku) bo prowadzenie było cały czas poniżej 1".
Dlatego są takie ostre i wąskie.