Periode und Amplitude der Radialgeschwindigkeitsvariation des Cepheiden Polaris
(v. Roland Bücke, Hamburg)
1. Einführung
Polaris (α Ursa Minor) gehört nicht nur aufgrund seiner exponierten Position nahe des Himmelspols seit über 100 Jahren zu den photometrisch und spektroskopisch intensiv beobachteten Sternen. Als Cepheid und Mitglied eines Mehrfachsystems mit einer spektroskopischen Komponente zieht er durch seine ungewöhnlichen Eigenschaften nach wie vor die Aufmerksamkeit von Fachastronomen auf sich. Eine Vielzahl von Publikationen bis in die jüngste Vergangenheit belegt dies überzeugend.
Polaris ist der hellste und zugleich uns am nächsten stehende Cepheid. Die Variabilität seiner Helligkeit wurde schon in der Mitte des 19. Jahrhunderts registriert (Seidel 1852, Schmidt 1856). 1911 wies Hertzsprung eine photometrische Periodizität von 4 Tagen nach und bereits 1913 erfolgte durch Pannekoek die Identifikation als Cepheid. Die Helligkeitsvariation nahm seither stetig ab und hat gegenwärtig mit ΔVvis ≈ 0.05 mag die kleinste bekannte Amplitude eines Cepheiden. Die Radialgeschwindigkeitsvariation von Polaris beobachtete als erster Champbell (1899). Diese ist ebenfalls von 5 – 6 km/s in der ersten Hälfte des 20. Jahrhunderts auf weniger als 2 km/s in den 1990er Jahren zusammen geschrumpft. Ferro (1983) und Dinshaw et. al. (1989) hielten sogar das gänzliche Verlassen der Instabilitätsphase mit der Einstellung der Pulsation für möglich. In den 1990er Jahren registrierte man ein Verharren bei ca. 1.5 km/s und neueste Beobachtungen (Engle et.al. 2004) belegen die gegenläufige Tendenz einer langsamen Zunahme der photometrischen Schwingungsamplitude, die sicher mit einer Vergrößerung der Radialgeschwindigkeitvariation korreliert.
Die in diesem Beitrag vorliegende Auswertung von Radialgeschwindigkeitsmessungen im Zeitraum von April bis Oktober 2005 bestätigt diese Tendenz. Eine Anpassung an eine Sinusfunktion ergab folgende Ephemeriden: Radialgeschwindigkeitsminimum JDRVmin = 2453484.82 ±0.04, Periode P = 3.9738 ±0.0009 Tage, Amplitude (Peak-to-Peak) A = 1.87 ±0.06 km/s.
2. Beobachtung
In der Zeit vom 21. April bis 05. Februar 2006 wurden 29 Radialgeschwindigkeiten (RV) gemessen (Tab. 1). Im Spektralbereich von 6050 bis 6600 Ǻ (Abb. 1) wurden hierfür die Linienzentren von 6 Fe I-Linien (6065.494, 6191.571, 6230.736, 6393.612, 6400.009 und 6411.658 Ǻ), 3 Ca I-Linien ( 6102.723, 6122.217 und 6162.173 Ǻ) sowie von der Hα-Linie bei 6562.808 Ǻ durch die Anpassung einer Gaußfunktion bestimmt.
Jede RV ist durch Mittelwertbildung aus mehreren Einzelmessungen (Spalte 2, Tab. 1) errechnet worden. Die Wellenlängenkalibrierung erfolgte mit dem Spektrum einer Neon-Glimmlampe. Zur Reduzierung des Fehlers durch das thermisch bedingte Abdriften des Spektrums während der gesamten Beobachtungsdauer wurde jede Einzelaufnahme mit Neonspektren kalibriert, die jeweils vor und nach dieser aufgenommen wurden. Die Unsicherheit der RV-Werte (Spalte 4, Tab. 1) resultiert dann aus der Standardabweichung der n Einzelmessungen.
Zur Beobachtung stand ein Newton-Reflektor mit 20 cm Öffnung zur Verfügung. Die Ankopplung des Gitterspektrographen mit einer Auflösung von 1.9 Ǻ an das Teleskop erfolgte mit einem 200 µm Lichtleiter (Bücke 2005).
| JD | n | RV | σ/√n | JD | n | RV | σ/√n |
| 2453000+ | (-) | (km/s) | (km/s) | 2453000+ | (-) | (km/s) | (km/s) |
| 482.41 | 9 | -19.002 | ±0.240 | 634.36 | 11 | -20.234 | ±0.226 |
| 484.38 | 9 | -17.537 | ±0.265 | 636.40 | 5 | -18.340 | ±0.262 |
| 485.37 | 8 | -17.092 | ±0.296 | 674.37 | 10 | -19.776 | ±0.271 |
| 503.40 | 11 | -18.482 | ±0.231 | 675.28 | 9 | -18.730 | ±0.357 |
| 518.45 | 18 | -19.246 | ±0.165 | 706.35 | 10 | -18.440 | ±0.160 |
| 519.42 | 11 | -19.148 | ±0.260 | 717.36 | 7 | -19.964 | ±0.500 |
| 612.47 | 8 | -18.501 | ±0.195 | 750.29 | 7 | -19.250 | ±0.230 |
| 613.37 | 10 | -19.756 | ±0.155 | 750.38 | 11 | -18.293 | ±0.216 |
| 613.46 | 8 | -20.094 | ±0.242 | 750.44 | 11 | -18.182 | ±0.175 |
| 614.37 | 10 | -19.969 | ±0.234 | 751.36 | 9 | -18.426 | ±0.271 |
| 617.36 | 11 | -20.011 | ±0.220 | 752.29 | 14 | -19.320 | ±0.548 |
| 619.39 | 9 | -18.825 | ±0.287 | 758.35 | 10 | -18.417 | ±0.276 |
| 627.34 | 14 | -18.038 | ±0.211 | 771.46 | 17 | -18.159 | ±0.139 |
| 630.34 | 8 | -20.222 | ±0.130 | 772.27 | 21 | -18.960 | ±0.131 |
| 631.33 | 9 | -18.837 | ±0.177 | ||||
| Tabelle 1: Messwerte | |||||||

Abb. 1: Spektrum von Polaris
3. Datenanalyse
An die Messwerte (Abb. 2) wurde eine Sinusfunktion nach der Methode der kleinsten Fehlerquadrate angefittet (Abb. 3). Diese Anpassung ergab folgende Ephemeriden für die RV-Variation: RV-Minimum JDRVmin = 2453484.80 ±0.067, Periode P = 3.9745 ±0.0013 Tage, Amplitude (Peak-to-Peak) A = 1.90 ±0.09 km/s. Eine lineare Regression an den Residuen der RV-Werte zeigt deutlich die Tendenz zur stetigen Abnahme der RV im Beobachtungszeitraum (Abb. 4).

Abb. 2: Radialgeschwindigkeitsmesswerte von Polaris

Abb. 3: Anpassung einer Sinusfunktion an die Messwerte

Abb. 4: Residuen der Sinusanpassung an die RV-Variation mit Regressionsgerade
Die RV-Variation im Beobachtungszeitraum resultiert also aus einer Überlagerung einer sinusförmigen Pulsation mit einer Periode P ~ 3,97 Tage und einem langfristigem Trend. Dieser resultiert mit hoher Wahrscheinlichkeit aus einem systematischen Fehler, bedingt durch die Dopplerverschiebung aufgrund der Bewegung Erde um die Sonne.
Diese stetige Veränderung der RV wurde unter Anwendung eines Polynoms dritter Ordnung (siehe Abb. 4) aus den Messwerten eliminiert. Eine erneute Sinusanpassung ergibt die vom langfristigen Trend beseitigten Parameter der Pulsation (Abb. 5):
RV-Minimum JDRVmin = 2453484.82 ±0.04, Periode P = 3.9738 ±0.0009 Tage, Amplitude (Peak-to-Peak) A = 1.87 ±0.06 km/s [1].

Abb. 5: Anpassung einer Sinusfunktion an die Messwerte nach Eliminierung des langfristigen Trends
4. Ergebnisse
4.1 Cepheid Polaris
Die gemessene Periode von P = 3.9738 ±0.0009 Tage weicht nur wenig von der Periode ab, die von Hatzes & Cochran (2000) für 1993 mit 3.9726769 ±0.00011 Tagen angegeben werden. Dinshaw at. al. (1989) fanden eine weitere Periode von 45.3 Tagen und Hatzes & Cochran (2000) eine Periode von 40.2 Tagen. Diese Perioden konnten durch die vorliegenden Messungen nicht nachgewiesen werden. Hierzu ist die Unsicherheit der Messungen zu groß bzw. die Anzahl der Messwerte zu klein.
Aus der Amplitude der Schwingungen lässt sich am sichersten eine Aussage über die Veränderungen der RV-Variation seit den zuletzt veröffentlichten Werten von 1997 (Kamper & Fernie 1998) ableiten. Die gemessenen Werte seit Beginn der Beobachtungen sind in Abb. 6 dargestellt. Die sprunghaft gestiegene Genauigkeit von RV-Messungen seit Mitte der 1980er Jahre ist beachtenswert. Der gemessene Wert von 1.87 ±0.06 km/s passt sich widerspruchsfrei in die zeitliche Entwicklung ein und zeigt eine Tendenz zu einer Vergrößerung der Schwingungsamplitude.

Abb. 6: RV-Amplitude von Polaris in der Zeit von 1958 bis 2005
4.2 Spektroskopischer Doppelstern Polaris
Die Periode der spektroskopischen Komponente des Mehrfachsystems Polaris beträgt nach Dinshaw et. al. (1989) 29.9 ±0.1 Jahre. Die einzige Datenquelle, aus deren sich über mehr als eine Periode hinweg der RV-Verlauf lückenlos konstruieren lässt, sind die RV-Werte, die aus den 1180 Platten des Lick Observatoriums am 36-Zoll Refraktor von 1896 - 1958 ermittelt wurden (Roemer 1964)(Abb. 7).

Abb 7: RV-Periode des spektroskopischen Doppelsterns Polaris
Die absolute mittlere RV dieser Arbeit passt sich in Anbetracht dessen, dass keine Kalibrierung an einem Standardstern vorgenommen wurde, gut in die Periode des spektroskopischen Doppelsterns ein.
5. Schlussbemerkung
Mit den in diesem Beitrag vorgestellten Ergebnissen zeichnet sich neben dem systematischen Monitoring von Be-Sternen ein weiteres wissenschaftlich orientiertes Betätigungsfeld für Amateurspektroskopiker ab. Es zeigt sich, dass selbst mit kleinen Teleskopen mit den heutigen Aufnahme- und Analysetechniken eine Präzision bei der Radialgeschwindigkeitsmessung erreicht werden kann, die noch vor wenigen Jahren nur der professionellen Astronomie mit Großgeräten möglich war.
Herrn Dr. Thomas Eversberg danke ich für die kritische Durchsicht des Manuskripts. Seine konstruktiven Vorschläge waren mir eine wertvolle Hilfe. Herrn Helmut Jahns von der Fachgruppe Computerastronomie des VdS danke ich für die Programmierung eines Windowsprogramms zur zuverlässigen Berechnung der Korrekturwerte zur Reduktion auf das heliozentrische Koordinatensystem.
Referenzen:
Bücke, R. (2005), VdS-J., 16, 81.
Champbell, W. W. (1899), Astrophys. J. 10, 180.
Dinshaw, N., Matthews, J. M., Walker, G. A. H., Hill, G. M. (1989), Astrophys. J. 98, 2249.
Engle, S. G., Guinan, E. F., Koch, R. H. (2004), American Astron. Soc. Meeting 204.
Ferro, A. A. (1983), Astrophys. J. 274, 755.
Hatzes, A. P., Cochran, W. D. (2000), Astron. J. 120, 979.
Hertzsprung, R. (1911), Astron. Nachr. 189, 89.
Kamper, K. W. & Fernie, J. D. (1998), Astrophys. J. 116, 936.
Kamper, K. W., Evans, N. R., Lyons, R. W. (1984), J. R. Astron. Soc. Can. 78, 173
Pannekoek, A. (1913), Königl. Akad. v. Wetenschappen, Amsterdam, 15, 1192.
Roemer, E. (1964), Astrophys. J. 141, 1415.
Schmidt, J. F. (1856), Astron. Nachr. 46, 293.
Seidel, L. (1852), Abh. d. II. Kl. d. königl. Akad. d. Wiss., München, 6, 564.
[1] Die angegebenen Unsicherheiten sind die mittleren Fehler des Mittelwertes der Residuen in RV-Richtung für die Amplitude, in der Zeitachse für das RV-Minimum und in der Zeitachse dividiert durch die Zahl der Perioden im Beobachtungszeitraum für die Unsicherheit der Periode.