La stella 80 della costellazione dell’Orsa Maggiore è di tipo A5V e si trova a 81.7 A.L. dalla nostra stella. La sua massa così come il raggio è praticamente il doppio di quello del Sole. La temperatura effettiva è pari a circa 8’500 K e forma un sistema binario con una stella nana rossa di classe spettrale M3. Inoltre recentemente si è dimostrato che Alcor e la sua compagna formano un sistema legato con Mizar (ζUMa).

L’immagine in figura rappresenta lo spettro di 80UMa ripreso il giorno 03 maggio 2012 alle ore 23.07 (TMEC) da Briosco (MB) con un telescopio Newton SkyWatcher 200mm f/4 + reticolo di diffrazione StarAnalyser 100 + camera Magzero MZ-5m. Una prima analisi dello spettro di assorbimento mostra in modo marcato la presenza della serie di Balmer ed in particolare le linee dalla Hα alla Hη. Di seguito è riportato l’istogramma dellla regione di interesse. L’offset ottenuto tramite fit gaussiano è 16.4161 pixel.

Spettro di assorbimento non calibrato (Visual Spec)

A questo punto si è proceduto alla realizzazione dello spettro MAX-assorbimento al fine di semplificare le future operazioni di fit. Prendiamo in particolare in considerazione la regione dello spettro tra la linea Hβ ed il vicino UV (circa 400 nm) mostrata nella figura sottostante:

Spettro di 80UMa nella regione compresa tra 380 e 520 pixel

Effettuando un fit gaussiano sui picchi otteniamo:

• 389.1 pixel = 3830.2 Å – linea Hη dell’HI (3835.0 Å)
• 394.6 pixel = 3886.5 Å – linea Hζ dell’HI (3888.6 Å)
• 402.7 pixel = 3968.9 Å – linea Hε dell’HI (3969.7 Å)
• 415.6 pixel = 4102.1 Å – linea Hδ dell’HI (4101.3 Å)
• intorno a 426 pixel – assorbimento anomalo (CaI 4227 Å?)
• 439.0 pixel = 4342.0 Å – linea Hγ dell’HI (4340.0 Å)
• intorno a 450 pixel – assorbimento anomalo (HeI 4471 Å )
• 457.9 pixel = 4535.9 Å – linea da identificare
• 467.4 pixel = 4632.9 Å – linea del FeII (4630 Å)
• 489.9 pixel = 4864.2 – linea Hβ dell’HI (4860.8 Å)

Consideriamo ora la seconda parte dello spettro dalla linea Hβ dell’Idrogeno a 800 nm, mostrato nella figura sottostante:

Spettro di 80UMa nella regione compresa tra 500 e 800 pixel

Effettuando un fit gaussiano sui picchi otteniamo:

• 512.2 pixel = 5092.5 Å – linea da identificare
• 521.9 pixel = 5191.8 Å – linea del FeII (5197.6 Å)
• 529.1 pixel = 5265.4 Å – linea da identificare
• 535.7 pixel = 5333.1 Å – linea da identificare
• 550.8 pixel = 5488.7 Å – linea da identificare (NiI 5476.9 Å?)
• 560.3 pixel = 5585.6 Å – linea da identificare (OI 5577 Å o OV 5572 – 5598 Å da nebulose, FeI 5572.8 Å)
• 581.8 pixel = 5806.2 Å – linea da identificare (HgI, banda tellurica O3, …?)
• 590.3 pixel = 5893.8 Å – linea dell’NaI (5890 Å – 5896 Å)
• intorno ai 615 pixel – assorbimento anomalo
• 625.7 pixel = 6256.2 Å – linea da identificare
• 655.0 pixel = 6557.6 – linea Hα dell’HI (6562.1 Å)
• 686.3 pixel = 6878.5 Å – banda tellurica O2 (6884 Å)
• 692.0 pixel = 6936.7 Å – linea da identificare
• intorno ai 705 pixel – assorbimento anomalo
• 707.6 pixel = 7096.3 Å – linea da identificare
• 716.2 pixel = 7184.5 Å – banda tellurica H2O (7160 – 7400 Å)
• 723.3 pixel = 7257.6 Å – banda tellurica H2O (7160 – 7400 Å)
• 726.5 pixel = 7291.1 Å – banda tellurica H2O (7160 – 7400 Å)
• 758.5 pixel = 7619.0 Å – banda tellurica O2 (7621 Å)
• 772.0 pixel = 7757.5 Å – linea dell’OI (7771 Å)
• 782.2 pixel = 7862.5 Å – linea da identificare
• 790.2 pixel = 7944.2 Å – linea da identificare

Consideriamo infine la terza parte dello spettro da 800 nm a 1000 nm, mostrato nella figura sottostante:

Spettro di 80UMa nella regione compresa tra 500 e 1000 pixel

Effettuando un fit gaussiano sui picchi otteniamo:

• 817.2 pixel = 8220.5 Å – linea limite di Paschen dell’HI (8204 Å)
• da 820 a 850 pixel struttura di assorbimento complessa
• 860.8 pixel = 8668.7 – linea del CaII (8662 Å)
• 867.0 pixel = 8731.5 – linea da identificare
• da 880 a 1000 struttura di assorbimento complessa

Il processo di calibrazione di uno spettro elettromagnetico è il primo passo per riuscire a dare una descrizione quantitativa di quanto acquisito per esempio con un reticolo di diffrazione come lo StarAnalyser 100. Per fare questo si è deciso di acquisire lo spettro di 9 stelle differenti di classe A. La scelta di questa classe spettrale deriva dal fatto che a quella temperatura sono ben visibili le linee della serie di Balmer dell’Idrogeno neutro HI, in grado con la banda tellurica dell’Ossigeno O2 a 7620 Å, di coprire l’intero spettro del visibile.

A questo punto la procedura seguita è la seguente:

• acquisizione dello spettro elettromagnetico della stella (file AVI) modificando tempi di esposizione e messa a fuoco al fine di ridurre eventuali errori sistematici;
• elaborazione del filmato AVI con Registax 6 al fine di ottenere l’immagine somma (file TIFF)
• rotazione dello spettro in modo da portarlo in posizione orizzontale. Selezione dello stesso ed espansione dell’immagine a 200 pixel di altezza;
• conversione del file TIFF in file PIC con IRIS;
• importazione del file PIC in Visual Spec, creazione dell’istogramma luminosità vs pixel in formato DAT;
• elaborazione di quest’ultimo con il software di analisi ROOT.

Per avere maggiori informazioni sui singoli passi riportati qui sopra fate riferimento agli altri post dedicati su ASTROtrezzi.

Obiettivo della calibrazione è trovare la relazione tra pixel e lunghezza d’onda λ della radiazione elettromagnetica. Fortunatamente per i reticoli di diffrazione questa relazione è in prima approssimazione di tipo lineare. Lunghezza d’onda zero corrisponde all’immagine dell’ordine zero della sorgente, ovvero nel nostro caso la stella in esame.

Quindi l’ideale sarebbe avere nel sistema di coordinate luminosità vs pixel come lunghezza d’onda zero il valore di pixel zero che in termini pratici si traduce nel porre la stella nell’origine del sistema di coordinate. Visual Spec purtroppo misura la distanza in pixel partendo dal lato sinistro dell’immagine che importiamo. Dato che la stella non è mai sul bordo, allora questa assumerà una posizione in pixel diversa da zero (offset). Per risolvere questo problema abbiamo fittato con una distribuzione gaussiana il picco di luminosità corrispondente all’immagine della stella in esame e shiftato lo spettro del valore medio di tale distribuzione (che dovrebbe coincidere con la posizione del centro della stella misurato in pixel). Dato che ogni stella in esame si troverà rispetto al lato sinistro dell’immagine acquisita ad una distanza diversa allora anche il relativo offset sarà differente. Di seguito riportiamo i valori degli offset per il campione di stelle utilizzato:

• 80 UMa (A5): 16.4 pixel
• ε UMa (A0): 17.2 pixel
• α Gem (A2): 19.4 pixel
• β Leo (A3): 24.3 pixel
• δ UMa (A3): 17.9 pixel
• β UMa (A1): 21.9 pixel
• ζ UMa (A0): 12.6 pixel
• γ UMa (A0): 29.8 pixel
• δ Leo (A3): 19.6 pixel

A questo punto tutti gli spettri hanno lo stesso offset e se volete potete sovrapporli. Da una rapida occhiata noterete che tutti i picchi della serie di Balmer stanno praticamente nella stessa posizione. Questo non è del tutto vero in quanto potrebbero esserci degli errori. Per prenderli correttamente in considerazione abbiamo pertanto determinato la posizione in pixel dei picchi Hα, Hβ, Hγ, Hδ, Hε, Hζ, Hη ed O2(7620 Å) per ogni stella in esame. I valori ottenuti nei singoli set di dati ed il valore medio con relativa incertezza statistica sono riportati in tabella:

Valore in pixel dei vari picchi di assorbimento ottenuti per mezzo di un fit gaussiano.

Ora non ci resta che riportare in un grafico il valore teorico delle lunghezze d’onda per la serie di Balmer (calcolabili analiticamente) e per la banda tellurica O2 approssimata ad un picco a λ = 7620 Å, in funzione del rispettivo valore misurato in pixel.

Un fit lineare del tipo y = mx permetterà di ottenere il valore sperimentale di m (dispersione del reticolo). Questo purtroppo non ha fornito risultati soddisfacenti portando ad errori sulla determinazione delle lunghezze d’onda spesso superiori a 30 Å. Si è quindi proceduto effettuando fit con polinomi di ordine sempre maggiore ottenendo al sesto grado un errore inferiore ai 10 Å. Questo ha il vantaggio di permettere misure precise di lunghezze d’onda, ma essendo un polinomio la relazione tra lunghezza d’onda e pixel perdeva di significato fisico per λ > 800 nm.

Dato che gli spettri acquisiti mostrano strutture anche per λ > 800 nm allora abbiamo abbandonato questa strada. Il risultato finale mostrato in figura è stato invece ottenuto utilizzando come curva di fit y = mx + q.

Determinazione della curva di calibrazione (fit lineare)

Il valore di q non nullo (comunque inferiore alla risoluzione dello strumento utilizzato) è da ricercarsi probabilmente in un offset intrinseco legato all’errore nella determinazione del centro della stella.

La curva di calibrazione che utilizzeremo per l’analisi degli spettri eseguiti con la strumentazione Newton SkyWather Wide Photo 200 mm f/4 + StarAnalyser 100 + camera Magzero MZ-5m è la seguente:

lunghezza d’onda = 10.2574387597*canale + 7.0294099901

Il massimo spred ottenuto tra i dati misurati e quelli teorici utilizzando questa curva di calibrazione è pari a 15.5 Å e sarà utilizzato come errore sulla misure effettuate. Un veloce test effettuato sulla stella αCMa ripresa con la medesima strumentazione mesi prima delle misure effettuate per la calibrazione mostra un accordo per la serie di Balmer e banda O2(7620 Å) con un errore inferiore ai 13.5 Å.