derivación de la función de masa binaria para órbita excéntrica

Al investigar métodos para detectar exoplanetas, aprendí sobre la función de masa binaria (BMF) que podría aplicarse para obtener la velocidad radial y la masa. He derivado BMF para una órbita circular, pero obtengo una respuesta incorrecta cuando se trata de una órbita excéntrica.

METRO 2 3 METRO t o t 2 = 4 π 2 GRAMO PAG 2 a 1 3
de:
GRAMO METRO 1 METRO 2 a ( 1 mi ) + 1 2 METRO 1 v 1 2 + 1 2 METRO 2 v 2 2 = GRAMO METRO 1 METRO 2 a ( 1 + mi ) + 1 2 METRO 1 v 1 2 + 1 2 METRO 2 v 2 2 METRO 1 v 1 a 1 ( 1 mi ) = METRO 1 v 1 a 1 ( 1 + mi ) METRO 2 v 2 a 2 ( 1 mi ) = METRO 2 v 2 a 2 ( 1 + mi ) METRO 1 v 1 = METRO 2 v 2 a = a 1 METRO t o t METRO 2
Obtuve:
a 1 = 1 + mi 1 mi GRAMO METRO 2 3 METRO t o t 2 v 1 2
entonces:
METRO 2 3 METRO t o t 2 = PAG v 1 3 2 π GRAMO ( 1 mi 1 + mi ) 3 / 2 = PAG k 3 2 π GRAMO pecado 3 i ( 1 mi 1 + mi ) 3 / 2
Sin embargo, el resultado de Wikipedia es:
F = METRO 2 3 pecado 3 i METRO t o t 2 = PAG k 3 2 π GRAMO ( 1 mi 2 ) 3 / 2

Me pregunto qué está mal con mi cálculo...

Hola y bienvenido a Astronomy SE. ¿Cuál es exactamente tu pregunta aquí?
@PierrePaquette la pregunta se ha actualizado muy bien; Creo que fue publicado accidentalmente antes de que estuviera terminado.

Respuestas (1)

Traté de encontrar una derivación detallada en línea después de mi falla. Finalmente lo encontré en The Exoplanet Handbook de Michael Perryman :

En realidad, K, o la llamada semiamplitud de velocidad radial , no es la velocidad radial máxima (que entendí antes) sino algún tipo de amplitud de la fluctuación de la velocidad radial . Y debe ser considerado en el espacio 3D.
sea ​​z = posición de la estrella solo la línea de visión = r ( t ) pecado ( ω + θ ) pecado i bosquejo

entonces la velocidad radial v r = z ˙ = pecado i [ r ˙ pecado ( ω + θ ) + r θ ˙ porque ( ω + θ ) ] .
por otro lado, tenemos:

r = a 1 ( 1 mi 2 ) 1 + mi porque θ porque θ = porque mi mi 1 mi porque mi mi mi pecado mi = 2 π PAG t
dónde mi = mi ( t ) es la llamada anomalía excéntrica .
Después de un largo cálculo, obtenemos:
v r ( t ) = 2 π a 1 PAG pecado i 1 mi 2 [ mi porque ω + porque ( ω + θ ( t ) ) ]
y aquí aparece K (amplitud):
k = 2 π a 1 PAG pecado i 1 mi 2

usa esto para sustituir a 1 en la Tercera Ley de Kepler:
METRO 2 3 METRO t o t 2 = 4 π 2 GRAMO PAG 2 a 1 3
por fin hay:
METRO 2 3 pecado 3 i METRO t o t 2 = PAG k 3 2 π GRAMO ( 1 mi 2 ) 3 / 2