La stella α della costellazione dei Gemelli è di tipo A1V e si trova a circa 51 A.L. dalla nostra stella. La sua massa così come il suo raggio è praticamente il doppio di quello del Sole. La temperatura effettiva è pari a 10’286 K e ruota su se stessa con una velocità di 5.2 km/s. Castore in realtà forma un sistema doppio con Castore B di classe spettrale A2Vm. Castore A è a sua volta un sistema doppio formato da due stelle orbitanti a soli 3 milioni di km di distanza con un periodo di 9.21 giorni. Castore Ab è una stella nana di classe spettrale incerta.

L’immagine in figura rappresenta lo spettro di αGem ripresa il giorno 03 maggio 2012 alle ore 22.17 (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 è 19.3740 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 αGem nella regione compresa tra 380 e 520 pixel

Effettuando un fit gaussiano sui picchi otteniamo:

• intorno ai 395 pixel – assorbimento anomalo (linea Hζ dell’HI 3888.6 Å)
• 403.5 pixel = 3946.8 Å – linea Hε dell’HI (3969.7 Å) difficoltà nel fit.
• 417.1 pixel = 4087.1 Å – linea Hδ dell’HI (4101.3 Å)
• 441.1 pixel = 4333.3 Å – linea Hγ dell’HI (4340.0 Å)
• intorno ai 470 pixel – assorbimento anomalo (HeI 4471 Å )
• 492.2 pixel = 4856.9 – 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 αGem nella regione compresa tra 500 e 800 pixel

Effettuando un fit gaussiano sui picchi otteniamo:

• 526.0 pixel = 5203.8 Å – linea del FeII (5197.6 Å)
• 552.2 pixel = 5472.0 Å – linea da identificare (NiI 5476.9 Å?)
• 584.3 pixel = 5801.3 Å – linea da identificare (HgI, banda tellurica O3, …?)
• 594.8 pixel = 5909.4 – linea dell’NaI (5890 Å – 5896 Å)
• 629.7 pixel = 6267.8 – assorbimento anomalo
• 658.5 pixel = 6562.8 Å – linea Hα dell’HI (6562.1 Å)
• intorno ai 690 pixel – assorbimento anomalo (banda tellurica O2 6884 Å?)
• 727.3 pixel = 7268.7 – banda tellurica H2O (7160 – 7400 Å)
• 761.6 pixel = 7620.1 Å – banda tellurica O2 (7621 Å)

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

Spettro di αGem nella regione compresa tra 500 e 1000 pixel

Effettuando un fit gaussiano sui picchi otteniamo:

• 818.5 pixel = 8204.1 Å – linea limite di Paschen dell’HI (8204 Å)
• tra i 845 ed i 850 pixel – struttura di assorbimento complessa
• 864.3 pixel = 8673.4 – linea del CaII (8662 Å)
• intorno ad 894 pixel – assorbimento anomalo
• intorno ai 930 pixel – assorbimento anomalo
• 968.2 pixel = 9739.7 Å – linea da identificare difficoltà nel fit

Riportiamo di seguito la seconda campagna di misure effettuate da Briosco (MB) il 03/12/2012 alle 22.05 (TMEC) . Telescopio Newton SkyWatcher 200mm f/4 + reticolo di diffrazione StarAnalyser 100 + camera Magzero MZ-5m. In figura è mostrato lo spettro di assorbimento e successivamente le linee identificate e misurate. Non sono stati considerati assorbimenti anomali. Offset misurato 19.5704 pixel, massimo dello spettro a 5119.6 Å. Scarica il file di testo della misura.

Spettro di αGem tra 300 e 1200 pixel

Linee identificate:

• 3894.7 Å
• 3976.9 Å
• 4107.5 Å
• 4346.1 Å
• 4867.0 Å
• 5206.4 Å
• 5498.8 Å
• 5553.8 Å
• 5601.3 Å
• 5792.3 Å
• 5900.8 Å
• 6261.0 Å
• 6551.6 Å
• 7287.5 Å
• 7604.9 Å

L’immagine visuale dello spettro elettromagnetico è riportata qui sotto.

Spettro di αGem ripreso il 03/12/2012

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 Å.