Pagina 2 di 5

Re: Il mio primo "diagramma H-R"!

MessaggioInviato: 13 mar 2018, 22:56
da Paolo
Ti ringrazio Lorenzo, le vostre preziose informazioni (tue ma anche degli altri amici fotometristi) hanno rappresentato un fortissimo stimolo per me. Senza di queste credo che non avrei avuto modo di fare un'esperienza del genere.

Ho letto adesso il tuo post ma avevo gia' preparato le due trasformazioni tenendo a mente i concetti contenuti nei tuoi fogli Excel, con alcune diversita' perche' il mio dato di base non e' una magnitudine strumentale ma direttamente la magnitudine stimata a valle della calibrazione (Astroart usa le magnitudini di un grande numero di stelle per trovare la funzione polinomiale Mag = k1 + k2 log (ADU)). Quindi, per ogni stella, ho calcolato la differenza tra la mag di catalogo (Apass B e V) e quella osservata (B e V).

Ho ordinato gli elenchi per luminosita' ed ho fatto diversi grafici considerando tutte le stelle (434), le 100 stelle piu' luminose dell'ammasso e infine le 50 stelle piu' luminose. L'ultima scelta mi sembra la migliore in quanto, ampliando il campione, c'e' troppa dispersione dei dati probabilmente per il modesto SNR delle stelle meno luminose (ricordo che ho usato il Mak 90 e 15 minuti di esposizione per filtro).

Trasformazione per il Baader V
Trasf_V.jpg

Trasformazione per l'Astronomik B
Trasf_B.jpg

Mi sembra che ci sia una bella differenza tra i due filtri!

Particolare l'intersezione della linea di tendenza sullo zero (Apass B - B osservata) intorno a B-V ~ 1.1 per il filtro V e B-V ~ 0.6 per il B. Me l'aspettavo nelle rispettive bande di colore (la calibrazione fotometrica del frame comunque non e' stata fatta con queste 50 stelle ma con tutte quelle che Astroart ha potuto riconoscere, dovrei provare a limitare lo SNR)...

Appena ho un attimo correggo le magnitudini B e V e aggiorno il diagramma colore-magnitudine con i nuovi dati. Non mi aspetto grandissime variazioni, non perche' le trasformazioni siano poco utili (al contrario!) ma perche' molta dispersione verra' comunque introdotta a causa dell'incertezza associata alle misure delle stelle deboli. Servirebbe senz'altro ripetere l'esperienza usando un'apertura maggiore e un tempo di integrazione piu' lungo.

Paolo

Re: Il mio primo "diagramma H-R"!

MessaggioInviato: 14 mar 2018, 0:47
da LFranco
Ciao Paolo,
il termine Tv = 0.0274 è piccolo e vicino a zero e questo significa che il filtro V (Baader) è conforme allo standard fotometrico. Il filtro B (Atronomik) invece mostra un termine Tb = 0.2108 e quindi si discosta dallo standard. Detto in altri termini ... si può usare il filtro V anche senza trasformazioni, mentre la stessa cosa non si può fare per il filtro B per il quale è consigliabile usare il coefficiente Tb.

Saluti
Lorenzo Franco

Re: Il mio primo "diagramma H-R"!

MessaggioInviato: 14 mar 2018, 17:43
da Paolo
Ok Lorenzo! Ho applicato la trasformazione e il grafico si e' modificato leggermente sugli estremi dei valori B-V (dove le misure B avevano la maggiore deviazione). Adesso il diagramma dovrebbe essere un po' piu' fedele. Nel centro rimane una certa dispersione dei dati. Sono convinto che molto dipende dalla scala dell'immagine, infatti questa non permette di isolare al meglio tante stelle che sul sensore sono troppo vicine fra loro. Usando un'ottica piu' grande (ed una focale maggiore) dovrebbe andare molto meglio.

M67_HR_schema_colore_trasf.png

Ho provato a stimare la distanza dell'ammasso. Sebbene sia bassa, l'estinzione interstellare non e' trascurabile, quindi ho corretto i grafici con il valore E(B-V) = 0.041 (https://arxiv.org/abs/0904.2907). Sarebbe stato interessante provare a ricavarlo con le misure B-V e V-R.

Il metodo si basa sulla comparazione del diagramma colore-magnitudine di M67 con un diagramma standard H-R in cui e' nota la relazione tra tipo spettrale/indice di colore e la magnitudine assoluta Mv. Un esempio nel sito sotto.

http://webhome.phy.duke.edu/~kolena/hou ... rties.html

HR_std.JPG


Approssimando, le magnitudini apparenti V dell'ammasso differiscono da quelle assolute del diagramma standard di un valore chiamato "modulo di distanza". Per ricavarlo occorre in pratica mettere i dati nello stesso grafico e spostare verticalmente (variando la magnitudine) il plot del diagramma standard fino a farlo coincidere con quello di M67. La regione dei due diagrammi da confrontare e' quella della sequenza principale (vedi tabella sopra).

Nel solito foglio Excel ho inserito una seconda serie contenente i valori B-V e relative Mv del diagramma H-R standard. Il delta magnitudine trovato per far sovrapporre i dati della sequenza e' di 9.6 (con i dati di M67 cosi' dispersi, specie nella parte bassa del grafico, il posizionamento e' approssimato).

M67_dist_9.6.PNG

Due simulazioni con valori del modulo di distanza (DM) differenti (solo per mostrare la traslazione).

M67_dist_8.5.PNG

M67_dist_10.5.PNG


Utilizziamo adesso la legge secondo cui l'intensita' della luce si attenua in modo inversamente proporzionale al quadrato della distanza (i relativi rapporti). Avendo a che fare con le magnitudini, la distanza (espressa in parsec) si calcola mediante la formula: d = 10^((V-Mv+5)/5). Per DM = 9.6 si ottiene d ~ 830 parsec. Su wikipedia e' riportata una distanza compresa tra 800-900 pc per l'ammasso M67.

Un'altra interessante esperienza, disponendo di dati osservativi migliori, e' la stima dell'eta' dell'ammasso. Occorre infatti un punto di turn-off piu' definito nel grafico in quanto viene usato per stimare la massa della stella che ha appena lasciato la sequenza principale (vedi relazione tra B-V e massa nella tabella standard sopra). In alternativa si possono usare le isocrone (https://it.wikipedia.org/wiki/Isocrona), identificando quella che meglio approssima il plot dei dati. In questo caso pero' occorrono dei calcoli piu' complessi che tengono conto dei modelli evolutivi stellari.

Paolo

Re: Il mio primo "diagramma H-R"!

MessaggioInviato: 14 mar 2018, 19:51
da LFranco
Bel lavoro Paolo,
questa si chiama astrofisica applicata.
Complimenti!

Saluti
Lorenzo Franco

Re: Il mio primo "diagramma H-R"!

MessaggioInviato: 14 mar 2018, 23:53
da Fab
...io penso che intorno ad un lavoro così, detto da uno come Paolo,
che ha sempre sostenuto che di fotometria non ci capiva nulla, c'è da rimanere
ammirati.

Grazie ;)

Re: Il mio primo "diagramma H-R"!

MessaggioInviato: 15 mar 2018, 12:25
da Paolo
Sono io che vi ringrazio! Il merito va a tutte le persone che condividono i propri lavori esponendoli in modo didattico e ai mezzi che oggi ci consentono di reperirli (il nostro ed altri forum, i siti web), oltre naturalmente alla strumentazione ed ai software di cui disponiamo. Siamo veramente fortunati, anni fa questo non si poteva fare e la curva di apprendimento sarebbe stata ben diversa...

In ogni caso la fotometria impiegata per questa esperienza e' realmente semplice. Parliamo di una fotometria differenziale, nulla a che vedere con quella AllSky dove devi tenere sotto controllo un notevole numero di parametri "ballerini" e per la quale e' mandatorio lavorare sotto un cielo perfetto.

Comunque... che emozione premere il pulsante "Inserisci grafico a dispersione" di Excel dopo aver selezionato un'ammasso di numeri apparentemente privi di senso e vedere apparire un diagramma H-R come quelli che nel corso degli anni ho trovato nei libri e nelle riviste di astronomia! Nonostante la fotometria semplificata e nonostante la mia scarsa esperienza.

Paolo

Re: Il mio primo "diagramma H-R"!

MessaggioInviato: 26 apr 2019, 19:15
da LFranco
Ciao a tutti,

approfittando di questo periodo festivo ho provato anch'io a costruire il diagramma H-R dell'ammasso aperto M67 partendo da osservazioni fatte qualche anno fa con i filtri fotometrici B,V,R. E' semplicemente meraviglioso vedere come di dispongono i vari punti delle magnitudini V rispetto all'indice di colore B-V. Il gomito (turn-off) mostra il punto dove le stelle di massa maggiore si sono evolute abbandonando la sequenza principale.
M67_HR_diagram_image.png

Nel grafico, oltre alla linea della sequenza principale, ho inserito anche le due tracce evolutive corrispondenti a stelle di 1 e 1.25 masse solari. In corrispondenza di quest'ultima vediamo il punto di turn-off che ci fornisce una indicazione sull'età dell'ammasso (circa 5-6 miliardi di anni). Infine il fit con la linea della sequenza principale ci fornisce una indicazione sul modulo di distanza (m-M) e quindi sulla distanza (circa 800-900 parsec).

Saluti
Lorenzo Franco

Re: Il mio primo "diagramma H-R"!

MessaggioInviato: 26 apr 2019, 20:39
da tonyven
Fantastico! Veramente molto interessante. Sono rimasto a bocca aperta. Non pensavo si potesse realizzare un H-R di un ammasso. Sono sicuro che ci sarebbe da scrivere un libro su questa cosa. Che argomento affascinante! Complimenti vivissimi!
Saluti Tonino