1 Programm des heutigen Tages 1 Erluterungen zum

  • Slides: 31
Download presentation
1

1

Programm des heutigen Tages 1. Erläuterungen zum Vd. S-FG-Anfängerprojekte durch Dr. Thomas Hunger 2.

Programm des heutigen Tages 1. Erläuterungen zum Vd. S-FG-Anfängerprojekte durch Dr. Thomas Hunger 2. Themenwünsche des heutigen Kurses 3. Feed back der Teilnehmer zum Erst-Kurs 4. Eingegangene Themenvorschläge mit Betonung auf mehr Praxis (WINDOWSProgramme!) Interpretation des Kontinuums in einem Sternspektrum Was kann man daraus lesen? Wiederholung: Reduktion, Normierung, Kalibration 5. Fragen: wie gewinne ich eine Spektralaufnahme? wie wird nachgeführt? wie ermittle ich die Belichtungszeit? welches ist die optimale Spaltbreite? welches Gitter muss ich verwenden? wie muss der Spektrograph am Teleskop montiert sein? 6. Anwendungsbeispiele zu den Programmen VSpec und IRIS (Hugo Kalbermatten) 7. Fragen und Gedankenaustausch der DADOS-Freunde 8. Demo des Echelle-Spektrographen von Tobias Feger 2

Reduktion von Sternspektren Bei der Bearbeitung und Auswertung von Sternspektren ist zu beachten, dass

Reduktion von Sternspektren Bei der Bearbeitung und Auswertung von Sternspektren ist zu beachten, dass jeder Detektor (ob CCD oder Film) und jeder Spektrograph (ob Gitter, Prisma usw. ) seine eigenen Charakteristiken besitzt, die im Verlauf der Datenreduktion berücksichtigt und korrigiert werden müssen 3

α Cygni Eigenes Spektrum Wahres Spektrum 4

α Cygni Eigenes Spektrum Wahres Spektrum 4

Das Sternspektrum erscheint auf dem CCD als schmaler Strich. Dieser Spektralfaden enthält leider nicht

Das Sternspektrum erscheint auf dem CCD als schmaler Strich. Dieser Spektralfaden enthält leider nicht nur das Signal vom Stern. Das CCD selber hinterlässt seine Spuren im Rohbild, auch das Teleskop und vor allem der Spektrograph beeinflussen den Intensitätsverlauf des Spektrums. 5

Eigenes Spektrum Instrum. Funktion = wahres Spektrum : Wahres Spektrum = Instrum. Funktion /

Eigenes Spektrum Instrum. Funktion = wahres Spektrum : Wahres Spektrum = Instrum. Funktion / Eigenes Spektrum 6

Das Rohbild, so wie es vom CCD kommt, ist geprägt von: 1. Dunkelstrom 2.

Das Rohbild, so wie es vom CCD kommt, ist geprägt von: 1. Dunkelstrom 2. Rauschen 3. CCD-Flat-Field (lokale Empfindlichkeitsunterschiede) 4. Vignettierungen 5. spektralen Empfindlichkeitsfunktion 6. Streulicht im Spektrographen 7. wellenlängenabhängiger Effizienz des Spektrographen 8. Licht vom Himmelshintergrund 9. Sternsignal 7

Dunkelstrom … kann ermittelt werden aus einem Mittelwert aus mehreren (4 -6) Dunkelstromaufnahmen, die

Dunkelstrom … kann ermittelt werden aus einem Mittelwert aus mehreren (4 -6) Dunkelstromaufnahmen, die mit der gleichen Belichtungszeit und der gleichen CCD-Temperatur wie das Sternspektrum aufgenommen wurden. Rauschen Jede CCD-Aufnahme enthält zusätzlich zum erwünschten Signal auch noch Rauschen. Dieses Rauschen kann etwas zurückgedrängt werden, wenn man Einzelbilder mittelt. Bei der Addition von N Aufnahmen steigt die Signalstärke um den Faktor N, das Rauschen aber nur um √N. Im gemittelten Bild wird also das Rauschen vermindert. CCD-Flatfield Mit dem CCD-Flatfield werden die lokalen Empfindlichkeitsunterschiede einzelner Pixel des CCD herausgerechnet. Dazu fertigen wir ein Flatfield des CCD an. Der Chip wird dabei ohne jede Optik gleichmäßig beleuchtet. Das CCD-Flatfield muss wieder ein sehr rauscharmer Mittelwert von mehreren Einzelaufnahmen sein, damit es wirklich die unterschiedlichen Pixelempfindlichkeiten repräsentiert. Wenn man einmal so ein rauscharmes CCDFlatfield erzeugt hat, kann man es über einen längeren Zeitraum verwenden. 8

Der Himmelshintergrund Die Linien des Nachthimmels verlaufen nicht nur im Bereich neben dem Sternspektrum,

Der Himmelshintergrund Die Linien des Nachthimmels verlaufen nicht nur im Bereich neben dem Sternspektrum, sondern auch über dieses hinweg. Ohne Korrektur würden somit dort Linien vorgetäuscht werden, die gar nicht vom Stern herrühren. 9

Bevor wir also das eindimensionale Spektrum erzeugen, werden Streulicht und Himmelshintergrund entfernt. Gute Resultate

Bevor wir also das eindimensionale Spektrum erzeugen, werden Streulicht und Himmelshintergrund entfernt. Gute Resultate werden schon mit einer einfachen Methode erzielt. In den Flächen oberhalb und unterhalb des Spektrums bestimmt man den Mittelwert der Pixelintensitäten und zieht diesen von allen Pixeln des Spektrums ab. Der Median ist dazu noch besser geeignet als das gewöhnliche arithmetische Mittel, weil beim Median einzelne Pixel mit stark abweichenden Werten (z. B. Cosmics) schwächer gewichtet werden 10

Um nun eindimensionales Spektrum zu erzeugen, könnte man einfach die Zeile mit der höchsten

Um nun eindimensionales Spektrum zu erzeugen, könnte man einfach die Zeile mit der höchsten Intensität extrahieren. Besser ist jedoch, über mehrere Pixelzeilen zu mitteln oder dies einfach aufzuaddieren. Damit unterdrückt man wirksam das Rauschen, weil Informationen aus einem größeren Bereich des Spektrums verwendet werden. Das Ergebnis einer solchen Extraktion ist in der Abbildung gezeigt. 11

Spektrum Himmel Subtraktion: Spektrum - Himmel 12

Spektrum Himmel Subtraktion: Spektrum - Himmel 12

Wenn man nun nach Abzug des Himmelshintergrundes eindimensionales Spektrum erhalten hat, muss der Kontinuumsverlauf

Wenn man nun nach Abzug des Himmelshintergrundes eindimensionales Spektrum erhalten hat, muss der Kontinuumsverlauf gefittet werden. Achtung: die Absorptionslinien mancher Sterne haben sehr breite Flügel, besonders die bei Sternen mit hohen Schwerebeschleunigungen (z. B. Weiße Zwerge), man darf daher nicht zu sehr in der Nähe dieser Spektrallinien den Verlauf des Kontinuums festmachen. Gewöhnlich bestimmt man per Mausklick interaktiv den Verlauf des Kontinuums im Spektrum. Das ist in linienarmen Spektren sicher 13 leichter, als im linienreichen Spektren

Spektrumscan Pseudokontinuum Normierung = Spektrumscan/Pseudokontinuum 14

Spektrumscan Pseudokontinuum Normierung = Spektrumscan/Pseudokontinuum 14

3750 4000 4250 4500 4750 5000 5250 5500 Wellenlänge [Ǻ] Kontinuumnormiertes Spektrum alle Intensitäten

3750 4000 4250 4500 4750 5000 5250 5500 Wellenlänge [Ǻ] Kontinuumnormiertes Spektrum alle Intensitäten des Kontinuums (außerhalb der Absorptionslinien) liegen auf oder in der Nähe des Wertes 1. 15

Praktisches Beispiel einer Kontinuumsnormierung Verwendete WINDOWS-Programme: 1. GIOTTO zur Bildung des Summenspektrums 2. MAXIM-DL

Praktisches Beispiel einer Kontinuumsnormierung Verwendete WINDOWS-Programme: 1. GIOTTO zur Bildung des Summenspektrums 2. MAXIM-DL zur Spektrum-Scannung 3. MK 32 zur Normierung des Pseudokontinuums 4. EXCEL zur Endbearbeitung 16

Kalibration in Wellenlänge Ein linienreiches Spektrum ist ohne Hilfsmittel kaum in der Wellenlänge zu

Kalibration in Wellenlänge Ein linienreiches Spektrum ist ohne Hilfsmittel kaum in der Wellenlänge zu kalibrieren. Wir benutzen statt dessen das einfacher aufgebaute Neon. Spektrum einer Glimmlampe und versuchen die Wellenlängen der abgebildeten Peaks in einem Graphen darzustellen. Methode 1: Kalibration mit EXCEL Methode 2: Kalibration mit VSpec 17

Frage: Was kann man im Kontinuum eines Sternspektrums lesen? Ein Stern verhält sich annähernd

Frage: Was kann man im Kontinuum eines Sternspektrums lesen? Ein Stern verhält sich annähernd wie ein Schwarzer Körper, das heißt, er emittiert ebensoviel Strahlung wie er absorbiert, er befindet sich im thermischen Gleichgewicht Die Verteilung der abgestrahlten Intensität I der Strahlung eines Schwarzen Körpers auf die Wellenlängen hängt nach dem Planckschen Strahlungsgesetz nur von der Temperatur ab: 18

Planck-Kurven für verschiedene Temperaturen. der Planckkonstanten h = 6. 625 x 10 -34 Js

Planck-Kurven für verschiedene Temperaturen. der Planckkonstanten h = 6. 625 x 10 -34 Js und der Lichtgeschwindigkeit c = 2, 9979 x 108 ms-1 19

Das Plancksche Strahlungsgesetz liefert also den Zusammenhang der abgestrahlten Strahlungsleistung I (λ, T) als

Das Plancksche Strahlungsgesetz liefert also den Zusammenhang der abgestrahlten Strahlungsleistung I (λ, T) als Funktion der Temperatur T und der Wellenlänge λ. Dabei liegt das Maximum der Ausstrahlung nach dem Wienschen Verschiebungsgesetz bei der Wellenlänge λmax = 2, 9 x 10 -3 [m. K] / T Das Stefan-Boltzmann-Gesetz beschreibt die gesamte Intensität der Strahlung eines schwarzen Körpers. I (T) = σT 4. mit Stefan-Boltzmann-Konstante σ = 5, 670× 10 -8 [W/m 2 K 4] 20

Fragen: 1. 2. 3. 4. 5. 6. 7. wie muss der Spektrograph am Teleskop

Fragen: 1. 2. 3. 4. 5. 6. 7. wie muss der Spektrograph am Teleskop montiert sein wie wird nachgeführt welches ist die optimale Spaltbreite welche Kamera kann ich verwenden welches Gitter muss ich verwenden wie ermittle ich die Belichtungszeit wie gewinne ich eine Spektralaufnahme 21

Wie muss der Spektrograph am Teleskop montiert sein ? 22

Wie muss der Spektrograph am Teleskop montiert sein ? 22

Wie wird nachgeführt ? (Forumdiskussion) 23

Wie wird nachgeführt ? (Forumdiskussion) 23

Welches ist die optimale Spaltbreite ? Mit der Spaltbreite beeinflusst man in erster Linie

Welches ist die optimale Spaltbreite ? Mit der Spaltbreite beeinflusst man in erster Linie die Auflösung, d. h. je schmaler der Spalt ist, desto weniger wird das Sternenlicht über die Wellenlänge "verschmiert" (z. B. durch Nachführfehler oder schlechtes Seeing) Bei der Festlegung der Spaltbreite sollte man berücksichtigen, dass mit enger werdendem Spalt immer weniger Licht in den Spektrographen gelangt, was natürlich bei gleicher Belichtungszeit zu einem schlechteren S/N führt. 24

Welches ist die optimale Spaltbreite ? Ein Spalt der größer als das Sternenbild ist,

Welches ist die optimale Spaltbreite ? Ein Spalt der größer als das Sternenbild ist, verliert seine wesentliche Funktion, nämlich eine konstante hohe Auflösung zu garantieren. Ein guter Kompromiss ist eine Spaltbreite etwas kleiner als der Durchmesser des stellaren Bildpunktes (Beugungsscheibchens) im Teleskopfokus. Der lineare Durchmesser Ø des Beugungsscheibchens, das bei einer Richtungsszintilation (Seeing) σ in der Brennebene des Teleskops entworfen wird, ist Ø = f * tan σ (σ im Bogenmass, f = Teleskopbrennweite) Beispiel am C 14 Teleskop: f = 4000 mm, σ = 2'' (in rad = 9, 7*10 -6 mm) Ø = 9, 7*10 -6 mm * 4000 mm = 39 µ 25

Welche Kamera kann ich verwenden ? Digitale Spiegelreflexkamera Canon EOS 20 D – Ein

Welche Kamera kann ich verwenden ? Digitale Spiegelreflexkamera Canon EOS 20 D – Ein Detektor für die Sternspektroskopie? (v. Bernd Hanisch, Frankfurt/Oder) erschienen in FG-Mitteilungsblatt SPEKTRUM 33/2007 26

Welches Gitter muss ich verwenden ? β α φ Der (Blaze)Winkel φ bestimmt die

Welches Gitter muss ich verwenden ? β α φ Der (Blaze)Winkel φ bestimmt die Wellenlänge λm des Energiemaximums der Strahlung λm = 2 sinφ / k * n k = Ordnung; n = Anzahl der Gitterlinien L/mm Beispiel: In der 1. Ordnung mit einem 600 L/mm Gitter und einem Blazewinkel φ von 22° liegt λm bei 630 nm für φ = 15° liegt λm bei 420 nm 27

Welches Gitter muss ich verwenden ? β α φ In der gewünschten Ordnung k

Welches Gitter muss ich verwenden ? β α φ In der gewünschten Ordnung k und der Wellenlänge λ gilt die Gittergleichung sinα + sinβ = kλ / g (g=Gitterkonstante) Die Winkeldispersion des Gitters ist unabhängig von der Wellenlänge dβ/dλ = k /g cosβ 28

Welches Gitter muss ich verwenden ? f Ob jek tiv Die Lineardispersion d in

Welches Gitter muss ich verwenden ? f Ob jek tiv Die Lineardispersion d in der Fokalebene (des CCD) ist gegeben durch: d = f k / g cosβ (f=Brennweite der Abbildungsoptik, g=Gitterkonstante, k=Ordnung) 29

Wie ermittle ich die Belichtungszeit ? Die erforderliche Belichtungszeit für eine Spektralaufnahme ist klar

Wie ermittle ich die Belichtungszeit ? Die erforderliche Belichtungszeit für eine Spektralaufnahme ist klar abhängig von der gesamten Durchlässigkeitschrakteristik von Teleskop & Spektrograph, den Öffnungsverhältnissen und der Detektorempfindlichkeit. Die Einflüsse dieser Faktoren können zwar näherungsweise rechnerisch bestimmt werden, jedoch scheint es für den Amateur eher angebracht, empirisch vorzugehen. 30

Wie gewinne ich eine Spektralaufnahme ? 31

Wie gewinne ich eine Spektralaufnahme ? 31