A Coordinate System for Dynamical Instabilities in Hierarchical Systems in REBOUND

Resumen Técnico: Un Sistema de Coordenadas para Inestabilidades Dinámicas en Sistemas Jerárquicos en REBOUND

Enunciado del Problema
El problema de los $N$ cuerpos sigue siendo un desafío central en la mecánica clásica, particularmente para simular sistemas planetarios a escalas de tiempo de gigaaños. Aunque los integradores híbridos Wisdom-Holman (WH) se han convertido en el estándar de oro para la integración a largo plazo al dividir el hamiltoniano en un término kepleriano dominante y un término de interacción perturbativo, su eficacia depende en gran medida de la elección del sistema de coordenadas.

• Las coordenadas de Jacobi ofrecen alta precisión al resolver el problema de dos cuerpos exactamente, pero son inflexibles; fallan cuando el orden orbital cambia debido a inestabilidades dinámicas.
• Las Coordenadas Helicéntricas Democráticas (DHC) permiten cruces orbitales y son el estándar para integradores híbridos (por ejemplo, MERCURY, SyMBA, TRACE), pero desempeñan un papel deficiente en sistemas que carecen de una masa central dominante única, como sistemas binarios anchos donde un sistema planetario orbita uno de los componentes de una estrella binaria.

En sistemas binarios anchos, las DHC no logran mantener la jerarquía necesaria para una integración eficiente, mientras que las coordenadas de Jacobi no pueden manejar los encuentros cercanos y las inestabilidades orbitales comunes en estos entornos. Existe la necesidad de un sistema de coordenadas que soporte encuentros cercanos entre cualquier par de cuerpos (planetas, estrella anfitriona y estrella compañera) mientras retiene los beneficios de velocidad de los métodos WH.

Metodología
Los autores implementan un sistema de coordenadas derivado por J. E. Chambers et al. (2002), denominado "Coordenadas de Binario Ancho" (WB), en el integrador híbrido TRACE dentro del paquete de código abierto REBOUND.

1. Transformación de Coordenadas:
   El sistema WB define nuevas coordenadas canónicas ($\vec{Q}$) y momentos ($\vec{P}$) que separan el movimiento baricéntrico de la estrella central, el movimiento relativo de la estrella compañera y el movimiento relativo de los planetas con respecto a la estrella anfitriona. Esto permite dividir el hamiltoniano en tres componentes:

   • $H_{Jump}$: Movimiento baricéntrico de la estrella central.
   • $H_{Interaction}$: Interacciones planeta-planeta y planeta-compañero.
   • $H_{Kepler}$: Movimiento kepleriano de los planetas y el compañero alrededor de la estrella anfitriona.
     En ausencia de un compañero, estas coordenadas se reducen a las DHC estándar.

2. Esquema de Integración:
   Los autores adaptan el esquema de conmutación reversible en el tiempo de TRACE a las coordenadas WB. El integrador procede en subpasos (estilo saltamontes) avanzando $H_{Interaction}$, $H_{Jump}$ y $H_{Kepler}$.

   • Lógica de Conmutación: Cuando ocurren encuentros cercanos (planeta-planeta, planeta-estrella o planeta-compañero), el integrador desplaza dinámicamente términos desde $H_{Interaction}$ hacia $H_{Kepler}$ para asegurar que el término kepleriano permanezca dominante.
   • Manejo de Encuentros:
     • Si un planeta se acerca a la estrella anfitriona o a la estrella compañera dentro de un radio crítico (valor predeterminado: 0.1 radios de Hill para el compañero), el sistema revierte a coordenadas inerciales y utiliza un integrador tradicional (Bulirsch-Stoer o IAS15) para resolver las ecuaciones completas de movimiento.
     • Si los planetas se acercan entre sí, se mueven términos de interacción específicos a $H_{Kepler}$ y se resuelven utilizando un integrador Bulirsch-Stoer, mientras el resto del sistema continúa con el solucionador WH.

3. Pruebas de Referencia:
   El rendimiento de TRACE en coordenadas WB fue evaluado frente a:

   • IAS15: Un integrador de alta precisión con paso de tiempo adaptativo (sirviendo como la verdad fundamental).
   • TRACE en DHC: La implementación estándar para comparación.
   • WHFast en coordenadas de Jacobi: Utilizado para comparación en regímenes estables.

Resultados Clave

• Análisis del Espacio de Parámetros: Una encuesta de 10,000 simulaciones que variaban la masa del compañero ($1 M_J$ a $3 M_\odot$) y el semieje mayor (50–1000 UA) reveló que las coordenadas WB superan significativamente a las DHC para compañeros más grandes y más distantes. Los autores proporcionan una guía de ley de potencia (Ecuación 11) para determinar cuándo son necesarias las coordenadas WB:
  $$\log_{10}\left(\frac{M}{M_\odot}\right) \approx 0.9 \log_{10}\left(\frac{a}{1 \text{ UA}}\right) - 4.67$$
  Las coordenadas WB desempeñan de manera competitiva con WHFast en coordenadas de Jacobi en todo el espacio de parámetros explorado.

• Dispersión de Planetas: En simulaciones que involucraban tres planetas de masa joviana y un compañero de 0.5 $M_\odot$, TRACE en DHC produjo resultados estadísticamente incorrectos (sobreestimando semiejes mayores e inclinaciones, subestimando excentricidades) en comparación con IAS15. TRACE en coordenadas WB recuperó resultados estadísticamente similares a IAS15 mientras lograba una aceleración de 9×.

• Aproximaciones Estelares: En pruebas que simulaban aproximaciones estelares sobre un planeta similar a Neptuno, TRACE en DHC no logró reproducir las predicciones analíticas para los cambios de energía, excentricidad e inclinación. TRACE en coordenadas WB coincidió tanto con las predicciones analíticas como con los resultados de IAS15, ofreciendo una aceleración de 2× sobre IAS15.

• Oscilaciones ZLK: Para los ciclos de von Zeipel-Lidov-Kozai (ZLK), los integradores de paso de tiempo fijo a menudo luchan con las excentricidades extremas alcanzadas durante las fases de alta excentricidad.

  • En DHC, el espacio de parámetros válido para que TRACE coincidiera con la precisión de IAS15 era extremadamente limitado.
  • En coordenadas WB, TRACE coincidió con la precisión de IAS15 para inclinaciones binarias de hasta $70^\circ$ ($e_{max} \sim 0.9$).
  • En un caso de prueba específico que involucraba un planeta interior similar a la Tierra (forzando un paso de tiempo global pequeño), TRACE en coordenadas WB fue 20× más rápido que IAS15 manteniendo la precisión, mientras que DHC no logró reproducir el período de oscilación correcto.

Significado y Directrices
El artículo afirma que la implementación de coordenadas WB en TRACE aborda un régimen específico donde los métodos híbridos tradicionales fallan: sistemas binarios anchos propensos a inestabilidades dinámicas y encuentros cercanos. Los autores enfatizan que para grandes conjuntos de sistemas caóticos en estrellas binarias, este método puede ser el único enfoque de integración viable que equilibra la precisión con la viabilidad computacional.

Los autores proporcionan las siguientes directrices para seleccionar integradores en REBOUND:

1. IAS15: Utilizar cuando el tiempo de cómputo no es una restricción o se requieren trayectorias exactas.
2. WHFast: Utilizar cuando el tiempo de cómputo es una preocupación y no se esperan encuentros cercanos.
3. TRACE: Utilizar cuando el tiempo de cómputo es una preocupación, se esperan encuentros cercanos y trayectorias cualitativamente correctas son suficientes.
   • Elección de Coordenadas: Para sistemas alrededor de estrellas similares al Sol, utilizar la Ecuación 11 para decidir entre DHC y WB. WB es fuertemente preferido para compañeros estelares reales.

El trabajo demuestra que TRACE en coordenadas WB ofrece precisión suficiente para estudiar inestabilidades dinámicas, formación planetaria y aproximaciones estelares en sistemas binarios sin los costos computacionales prohibitivos asociados con integradores adaptativos de alta precisión.