Kick Velocities and Mass Function of Free-Floating Planets from Dynamical Ejection in Hierarchical Three-Body Systems

Resumen Técnico: Velocidades de Patada y Función de Masa de Planetas de Flotación Libre a partir de la Eyección Dinámica en Sistemas de Tres Cuerpos Jerárquicos

Planteamiento del Problema
Los planetas de flotación libre (FFP, por sus siglas en inglés), o planetas errantes, son objetos subestelares no ligados a ninguna estrella anfitrona. Sus distribuciones de velocidad y masa sirven como fósiles cinemáticos de su historia de formación y eyección. Si bien el microlente gravitacional ha detectado una población de FFP que abarca desde masas terrestas hasta masas jovianas, los mecanismos dinámicos que gobiernan sus velocidades de eyección y la función de masa resultante aún no se han caracterizado completamente. Específicamente, existe la necesidad de mapear los resultados de la eyección de sistemas de tres cuerpos jerárquicos (una estrella anfitrona, un perturbador gigante masivo y un planeta interno más ligero) con las propiedades observables de la población galáctica de FFP, incluyendo sus velocidades asintóticas ($v_\infty$) y las escalas de tiempo de eyección ($t_{ejec}$).

Metodología
Los autores presentan un conjunto de simulaciones directas de $N$-cuerpos utilizando el paquete de integración rebound con el integrador IAS15 (paso de tiempo adaptativo de orden 15) para asegurar una alta precisión durante los encuentros cercanos. El estudio se centra en el problema de los tres cuerpos restringidos, tratando al planeta interno como una partícula de prueba respecto a una estrella anfitrona de masa solar ($M_\star = 1 M_\odot$) y un perturbador gigante masivo.

La rejilla de simulación está controlada y es sistemática:

• Parámetros Fijos: Masa de la anfitrona ($1 M_\odot$), semieje mayor del planeta gigante ($a_J = 5.2$ UA, un análogo de Júpiter) e inclinaciones iniciales ($0^\circ$–$3^\circ$).
• Parámetros Variados:
  • Masa del planeta eyectado ($m_p$): $10^{-5}$ a $10^{-1} M_{Jupiter}$.
  • Masa del planeta gigante ($M_J$): $1$a$75 M_{Jupiter}$.
  • Semieje mayor del planeta eyectado ($a_p$): $0.7$a$1.5 a_J$.
  • Excentricidades de ambos planetas ($e_p$ y $e_J$): $0.0$a$0.9$.
• Realizaciones: Cada configuración se simula con $N=10,000$ realizaciones independientes utilizando semillas aleatorias compartidas para las anomalías medias con el fin de aislar los efectos de los parámetros.
• Métricas: El estudio mide la velocidad de eyección asintótica ($v_\infty$), la velocidad de patada de radio finito ($v_{kick}$) en el momento de la confirmación de la eyección, y la escala de tiempo de eyección ($t_{ejec}$). La eyección se define por la positividad de la energía específica, una distancia de la estrella que excede $10 \times$ el semieje mayor máximo del sistema, y una velocidad radial hacia afuera que persiste durante diez periodos orbitales.

Marco Analítico
Los resultados numéricos se interpretan a través de un marco semianalítico que combina:

1. Presupuesto de Energía: El incremento de energía requerido para la eyección se deriva de la diferencia entre la energía orbital ligada y la energía cinética asintótica.
2. Honda Gravitatoria (Slingshot): La velocidad de patada máxima está limitada por el mecanismo de honda en el sistema de referencia del gigante, dependiendo de la velocidad de encuentro $u$ y la velocidad orbital del gigante.
3. Dispersión de la Escala de Hill: La velocidad media de eyección está gobernada por encuentros típicos dentro de la esfera de Hill del gigante, escalando como $v_H \propto M_J^{1/3}$.
4. Parámetro de Tisserand: Utilizado como un diagnóstico conservado para relacionar los elementos orbitales con las velocidades de encuentro.
5. Evolución Secular: Las escalas de tiempo para la inestabilidad están vinculadas al bombeo secular y al solapamiento de resonancias de movimiento medio (MMR).

Resultos Clave

1. Límite de Partícula de Prueba: La masa del planeta eyectado ($m_p$) tiene un impacto insignificante en las estadísticas de eyección sobre cuatro órdenes de magnitud ($10^{-5}$ a $10^{-1} M_{Jupiter}$). La fracción de eyección varía en menos de un 2.3%, y las velocidades medias permanecen estables. Esto confirma que para $m_p/M_J \lesssim 10^{-2}$, la dinámica está dominada por el perturbador gigante, y el cuerpo eyectado se comporta como una partícula de prueba.
2. Escalamiento de la Masa del Planeta Gigante: La masa del planeta gigante ($M_J$) establece la escala y el tiempo de eyección.
   • La velocidad media de eyección escala como $\langle v_\infty \rangle \propto M_J^{0.42}$.
   • La escala de tiempo de eyección escala como $\langle t_{ejec} \rangle \propto M_J^{-1.09}$, lo cual es consistente con la predicción secular $t_{sec} \propto M_J^{-1}$.
   • La fracción de eyección es no monotónica, alcanzando su máximo cerca de $10 M_{Jupiter}$ y disminuyendo en masas muy bajas y muy altas debido a las limitaciones de la escala de tiempo secular y la captura en la esfera de Hill, respectivamente.
3. Dependencia del Semieje Mayor: Para configuraciones dentro de la zona de solapamiento de resonancia ($|a_p - a_J|/a_J \lesssim 0.40$ para un perturbador de $10 M_{Jupiter}$), los tiempos de eyección son constantes ($\sim 1.9 \times 10^4$ años), gobernados por la escala de tiempo de solapamiento caótico rápido en lugar de la evolución secular. Las velocidades medias alcanzan su máximo cerca de $a_p = 0.9 a_J$.
4. Efectos de la Excentricidad (Respuesta de Dos Velocidades):
   • Excentricidad del Planeta Interno ($e_p$): Afecta principalmente a la cola de alta velocidad. Aumentar $e_p$ de 0 a 0.9 eleva la velocidad de patada máxima por un factor de $\sim 3$ (de $\sim 7.8$ a $\sim 23$ km/s) pero solo aumenta la velocidad media por un factor de $\sim 1.7$.
   • Excentricidad del Planeta Gigante ($e_J$): Tiene un efecto sustancialmente más fuerte en la cola. Aumentar $e_J$ de 0 a 0.9 eleva la velocidad de patada máxima por un factor de $\sim 10$ (de $\sim 7.8$ a $\sim 80$ km/s). Esto es impulsado por el crecimiento no lineal de la velocidad del pericentro del gigante ($v_{peri} \propto \sqrt{(1+e_J)/(1-e_J)}$), lo que permite encuentros de honda de alta energía poco frecuentes.
5. Envolvente Superior Analítica: Los autores derivan una fórmula analítica para la velocidad máxima de eyección basada en el mecanismo de honda y el parámetro de Tisserand. Esta envolvente predice con precisión el techo de las distribuciones de velocidad simuladas, particularmente para $e_p \gtrsim 0.5$.

Significancia e Implicaciones
El artículo afirma que estos resultados proporcionan un vínculo directo entre la arquitectura de los sistemas planetarios y las propiedades observables de la población galáctica de FFP:

• Función de Masa: Dado que la fracción de eyección es independiente de la masa del planeta, la función de masa actual de los FFP traza directamente la tasa de ocurrencia planetaria subyacente ($dN_{FFP}/dm_p \approx \bar{f}_{eject} \cdot dN_{planet}/dm_p$). Esto implica que si la tasa de ocurrencia aumenta hacia masas menores (como sugiere la data de Kepler), la población de FFP de masa sub-terrestre será numéricamente dominante.
• Firmas de Microlente: La distribución transformada del tiempo de escala de microlente ($dN/d\log t_E$) será plana en $\log t_E$ si la tasa de ocurrencia es plana, o estará aumentada en escalas de tiempo cortas si la tasa de ocurrencia aumenta hacia masas menores. Esto ofrece una predicción contrastable para el Roman Wide Field Survey y Euclid.
• Cinemática: La mayoría de los planetas eyectados tienen velocidades asintóticas de solo unos pocos km/s, lo que los hace cinemáticamente indistinguibles de su población estelar de origen (cambiando la dispersión de velocidad 3D en un $\sim 1\%$). Sin embargo, los sistemas con perturbadores gigantes altamente excéntricos ($e_J \gtrsim 0.5$) pueden producir una cola de alta velocidad ($v_\infty \sim 80$ km/s), creando potencialmente valores atípicos de movimiento propio detectables por Gaia o Roman.
• Escalas de Tiempo de Eyección: La mayoría de las eyecciones de sistemas masivos o excéntricos ocurren rápidamente ($10^3$–$10^4$ años), a menudo durante la fase de pre-secuencia principal. Por consiguiente, la distribución espacial de la población general de FFP es efectivamente indistinguible del campo, con solo los objetos jóvenes recientemente eyectados reteniendo memoria de su entorno de nacimiento.

Limitaciones
Los autores señen que la restricción a exactamente tres cuerpos es la principal limitación, ya que los sistemas reales pueden involucrar múltiples planetas que conducen a inestabilidades en cascada. Además, las simulaciones no incluyen discos protoplanetarios, pérdida de masa estelar, ni un rango completo de inclinaciones (específicamente el mecanismo de Kozai-Lidov), lo que podría alterar las estadísticas de eyección. Se identifica como necesario un trabajo futuro para incorporar estos factores en un modelo completo de síntesis de población.