Explorando el Planeta Nueve con Python usando poliastro

This post was originally published on this site

No cuesta nada repetirlo: no se conoce la órbita del Planeta Nueve. Batygin y Brown en sus análisis lo que han hecho ha sido, parafraseando a Daniel Marín, estrechar el cerco en torno a él. Sin embargo, en el artículo trabajan con valores característicos y eso sí podemos aprovecharlo – hasta cierto punto, como en seguida veremos.

En el resumen del artículo tenemos esta frase:

We find that the observed orbital alignment can be maintained by a distant eccentric planet with mass $gtrsim 10 m_{bigoplus}$ whose orbit lies in approximately the same plane as those of the distant KBOs, but whose perihelion is $180^{circ}$ away from the perihelia of the minor bodies. [Énfasis mío]

Al final del artículo, en la figura 9, tenemos los datos de la órbita tentativa utilizada para las simulaciones:

  • Ascensión del nodo $Omega = 100^{circ}$
  • Argumento del perigeo $omega = 150^{circ}$
  • Inclinación $i = 30^{circ}$
  • Semieje mayor $a = 700~text{AU}$
  • Excentricidad $e = 0.6$

(La masa que tenga el planeta no nos afecta)

Con esto ya podríamos representar la forma de la órbita, pero nos falta un parámetro crucial… ¿dónde está el Planeta Nueve dentro de esta órbita? Si lo supiéramos no estaríamos hablando de especulaciones. Averiguarlo no es fácil, porque con este semieje mayor estamos hablando de unos períodos larguísimos. ¡Manos a la obra!

Comments are closed.