Ich frage mich, wie man die Exzentrizität eines Exoplaneten anhand seines Radialgeschwindigkeits-gegen-Phasen-Diagramms berechnet. Um meine Frage zu klären, nehme ich als Beispiel einen Exoplaneten namens WASP-14b 2 ( http://exoplanets.org/detail/WASP-14_b ).
Ein Diagramm der Radialgeschwindigkeit des Sterns über der Phase wird in der oberen linken Ecke angezeigt. Ich frage mich, wie ich möglicherweise die Exzentrizität des Exoplaneten mit diesem Diagramm (oder einigen anderen Werten aus den Originalmessungen) berechnen könnte. Ich habe ein paar Möglichkeiten gefunden, die Exzentrizität zu berechnen:
Dies verwendet den Exzentrizitätsvektor, der mit dieser Formel berechnet wird:
Das Problem dabei ist, dass diese Formel den spezifischen Drehimpulsvektor und den Positionsvektor benötigt, die ich aufgrund der Messungen nicht kenne. Es gibt jedoch eine andere Möglichkeit, die Exzentrizität zu berechnen:
Wo ist der Radius der Apoapsis und der Radius der Periose. Diese Werte sind nicht nur anhand der Messungen bekannt, aber ich glaube, es sollte möglich sein, sie zu berechnen, indem man das Integral der Sinusfunktion (Radialgeschwindigkeit vs. Phase) nimmt. Dies würde mir die Position des Sterns zu jedem beliebigen Zeitpunkt geben. Das Problem ist, dass ich die genauen Punkte, die in der Grafik angezeigt werden, nirgendwo finden kann, geschweige denn eine Sinusfunktion, die dazu passen würde.
Wenn ich ein Integral der Funktion bekomme, muss ich noch eines für den Planeten selbst erstellen, da dies die Bewegung des Sterns beschreibt. Ich kann die Masse des Planeten mit der folgenden Formel berechnen:
was mir die Entfernung zwischen dem Stern und dem Planeten gibt. Als nächstes kann ich die Geschwindigkeit des Planeten berechnen mit:
Und danach kann ich die Masse des Planeten mit dieser Formel berechnen:
Aber hier taucht ein weiteres Problem auf, wie es in einem Wikipedia-Artikel über Doppler-Spektroskopie heißt: „Beobachtungen eines echten Sterns würden ein ähnliches Diagramm erzeugen, obwohl Exzentrizität in der Umlaufbahn die Kurve verzerren und die Berechnungen unten erschweren wird.“
Wo finde ich die korrigierten Berechnungen und wie kann ich möglicherweise die Exzentrizität dieses Planeten nur mit diesen Werten berechnen ( und die Handlung, von der ich die genauen Punkte nicht finden kann)?
Weitere Quellen: http://adsabs.harvard.edu/abs/2009MNRAS.392.1532J
Es gibt eine Reihe von Optionen, wenn Sie eine Standardlösung zum Anpassen von RV-Kurven wünschen. Die vielleicht beste kostenlose ist Systemic Console .
Es ist jedoch nicht allzu schwer, etwas Grundlegendes selbst zu tun.
Definiere zunächst einige Begriffe:
ist die wahre Anomalie - der Winkel zwischen dem Perizentrum und der Position des Körpers um seine Umlaufbahn, gemessen vom Massenschwerpunkt der Ellipse.
ist die exzentrische Anomalie und wird durch die Gleichung definiert
Die mittlere Anomalie wird von gegeben
Das sagt uns die "Keplersche Gleichung".
Schließlich ist die Radialgeschwindigkeit gegeben durch
OK, also das Problem ist, dass die Radialgeschwindigkeit nicht explizit davon abhängt , sondern auf . Was Sie also tun, ist Folgendes:
Wählen Sie Werte für , , , Und ; Dies sind Ihre "freien Parameter, die die Umlaufbahn beschreiben. Je näher Sie Ihrer anfänglichen Schätzung kommen, desto besser.
Sie verwenden diese Parameter, um vorherzusagen, wie hoch die Radialgeschwindigkeiten zum Zeitpunkt der Beobachtung Ihrer RV-Datenpunkte sein würden. Sie tun dies, indem Sie rechnen unter Verwendung der obigen Gleichungen. Beginne mit der zweiten Gleichung und rechne . Dann müssen Sie die dritte Gleichung lösen, um zu erhalten . Dies ist transzendent, also müssen Sie eine Newton-Raphson-Methode oder etwas Ähnliches verwenden, um die Lösung zu finden. Sobald du hast dann verwenden Sie die erste Gleichung, um zu finden . Verwenden Sie dann die vierte Gleichung zur Berechnung zu jeder Ihrer Datenpunktzeiten.
Berechnen Sie ein Chi-Quadrat (oder eine ähnliche Gütezahl) aus dem Vergleich der vorhergesagten und gemessenen Werte von .
Iterieren Sie die Werte der freien Parameter und gehen Sie zurück zu Schritt 2. Fahren Sie fort, bis Ihre Anpassung konvergiert.
ProfRob