The quest for Magrathea planets. II. Orbital stability of exoplanets formed around double white dwarfs

Resumen Técnico: Estabilidad Orbital de Exoplanetas Formados alrededor de Dobles Enanas Blancas

Planteamiento del Problema
Aunque los modelos teóricos sugieren que los planetas de "segunda generación" pueden formarse a partir de discos circumbinarios creados por el material expulsado durante la fase de envolvente común de binarias de enanas blancas dobles (DWD, por sus siglas en inglés), la estabilidad dinámica a largo plazo de tales sistemas permanece sin probar. Estudios previos han demostrado el potencial de formación de planetas gigantes en estos entornos (Ledda et al. 2023), pero no está claro si estos sistemas pueden sobrevivir a las interacciones gravitatorias caóticas inherentes a las arquitecturas circumbinarias de múltiples planetas durante escalas de tiempo de millones de años. Además, la estabilidad de estos sistemas es crítica para evaluar su detectabilidad por la próxima misión de la LISA (Laser Interferometer Space Antenna), la cual tiene como objetivo detectar exoplanetas orbitando DWDs mediante modulaciones de frecuencia de ondas gravitacionales (GW).

Metodología
Los autores realizaron simulaciones de $N$-cuerpos de sistemas planetarios circumbinarios orbitando tres configuraciones específicas de DWD (DWD2, DWD3 y DWD4) extraídas de la población de DWD de LISA. El estudio empleó un integrador simpléctico híbrido (Chambers et al. 2002) adaptado específicamente para sistemas binarios para manejar tanto la evolución Kepleriana a largo plazo como los encuentros cercanos/colisiones.

Las simulaciones cubrieron cuatro poblaciones distintas basadas en condiciones iniciales:

• Población A: Sistemas no resonantes con planetas de menor masa ($M_p \leq 1.2 M_J$) alrededor de binarias de menor masa ($M_{bin} < 0.6 M_\odot$).
• Población B: Sistemas no resonantes con planetas de mayor masa ($M_p \geq 2.5 M_J$) alrededor de binarias más masivas ($M_{bin} \sim 1 M_\odot$).
• Población C: Contrapartes resonantes de la Población B, donde los planetas están inicialmente bloqueados en resonancias de movimiento medio 2:1.
• Población D: Sistemas de cinco planetas en una resonancia específica de "NMS" (una combinación de resonancias 4:3 y 5:3) alrededor de DWD*4.

Las condiciones iniciales incluyeron distribuciones de masas planetarias basadas en Ledda et al. (2023), semiejes mayores entre el límite de estabilidad crítica ($a_{crit}$) y 1 ua (o escalados para casos resonantes), y bajas excentricidades iniciales ($e < 0.01$) e inclinaciones ($1^\circ-3^\circ$). Las simulaciones se ejecutaron durante $3 \times 10^6$ a $10^7$ años.

Para analizar los resultados, los autores utilizaron varias métricas:

1. Déficit de Momento Angular Normalizado (NAMD): Para cuantificar la excitación dinámica (temperatura) del sistema.
2. Seguimiento de Eventos Catastróficos: Registro de expulsiones (dispersión más allá de 100 ua), colisiones y aproximaciones cercanas a la binaria.
3. Espaciamiento Orbital y Concentración de Masa: Análisis de los cambios en el espaciamiento orbital ($S_s$) y la fracción de masa del planeta más masivo ($S_m$).
4. Desplazamiento del Centro de Masa (CoM): Cálculo del desplazamiento del baricentro del sistema planetario para inferir qué planetas se perdieron.

Resultados Clave

• Estabilidad y Multiplicidad: Los sistemas planetarios que orbitan DWDs pueden permanecer estables durante los tiempos de estudio, pero la estabilidad depende fuertemente de la multiplicidad inicial. Los sistemas que comienzan con dos planetas son los más propensos a sobrevivir sin alteraciones. Por el contrario, los sistemas con tres, cuatro o cinco planetas experimentan frecuentemente inestabilidades catastróficas.
• Resultados de las Poblaciones:
  • Población A (Baja masa, no resonante): El 82.5% retuvo su conteo inicial de planetas; el 17.5% experimentó eventos catastróficos.
  • Población B (Alta masa, no resonante): Solo el 49% retuvo su multiplicidad inicial; el 51% experimentó pérdida de planetas.
  • Población C (Alta masa, resonante): El 65.6% retuvo su multiplicidad inicial. Aunque la resonancia mejoró la estabilidad para sistemas de baja multiplicidad, los sistemas resonantes de cuatro planetas fueron más propensos a la disrupción que sus contrapartes no resonantes.
  • Población D (Cinco planetas resonantes): El 100% de los sistemas experimentó eventos catastróficos. La mayoría (58%) perdió exactamente dos planetas y el 17.2% fue completamente disrumpido (todos los planetas perdidos).
• Evolución Dinámica: Los sistemas que sufrieron inestabilidad típicamente evolucionaron hacia estados "dinámicamente calientes" (NAMD > $1.3 \times 10^{-3}$), caracterizados por un mayor espaciamiento orbital y una mayor concentración de masa en el planeta sobreviviente más masivo. Se observó una fuerte anticorrelación entre la multiplicidad final y el NAMD: los sistemas que terminaron con menos planetas generalmente exhibieron una mayor excitación dinámica.
• Efectos de Resonancia: Si bien las resonancias 2:1 iniciales ayudaron a preservar los sistemas de baja multiplicidad (2-3 planetas), no evitaron la disrupción de las cadenas de mayor multiplicidad (4-5 planetas), las cuales casi universalmente rompieron la resonancia y perdieron planetas.
• Detectabilidad por LISA:
  • Sistemas de un Solo Planeta: Aproximadamente el 7% de los sistemas simulados terminaron como sistemas de un solo planeta. La mayoría de estos sobrevivientes poseen masas y periodos que entran en la ventana de detección de LISA, siempre que la ubicación en el cielo y la distancia sean favorables.
  • Sistemas de Multi-planetas: Los autores derivaron una fórmula (Ec. 15) para estimar la amplitud global del desplazamiento de frecuencia de onda gravitacional ($\Delta f$) inducido por múltiples planetas. Para la mayoría de los sistemas simulados, este desplazamiento está por debajo de la resolución de frecuencia de LISA ($\sim 7.92 \times 10^{-9}$ Hz para una misión de 4 años). Solo un puñado de sistemas resonantes (Población C) mostró desplazamientos potencialmente detectables dentro de una misión de 8 años.
  • Aceleración del CoM: El estudio señaló que las expulsiones planetarias pueden inducir aceleraciones medibles en el centro de masa de la binaria (hasta $\sim 0.9$ ua/año$^2$ en casos específicos), lo que podría dejar una impronta en la señal de GW, aunque aún no se ha definido un umbral de detección específico para LISA.

Significancia y Reivindicaciones
El artículo afirma que los sistemas planetarios de segunda generación alrededor de DWD son candidatos dinámicamente viables para la observación, pero sus arquitecturas finales están fuertemente esculpidas por las inestabilidades dinámicas tempranas. Los autores enfatizan que:

1. La Supervivencia es Común pero Modificada: Aunque la disrupción completa es rara (7.7% del total de los sistemas), la pérdida de planetas es común, lo que conduce a un aumento significativo (122%) en la población de sistemas de dos planetas derivados de progenitores de mayor multiplicidad.
2. NAMD como Diagnóstico: La métrica NAMD sirve como una herramienta crucial para distinguir entre sistemas que se formaron con baja multiplicidad y aquellos que fueron dinámicamente excitados y perdieron planetas. Los valores altos de NAMD en sistemas de dos planetas observados pueden indicar un historial de inestabilidad.
3. Potencial de LISA: El estudio respalda la viabilidad de detectar planetas "Magrathea" (planetas de segunda generación alrededor de DWD) con LISA, particularmente los sobrevivientes de un solo planeta. Sin embargo, la detección de sistemas de múltiples planetas mediante modulación de frecuencia es poco probable para la mayoría de las configuraciones debido a los pequeños desplazos de frecuencia, lo que requiere estudios de forma de onda específicos y, potencialmente, depender de las firmas de aceleración del CoM.

Los autores concluyen que, si bien sus simulaciones cubren la evolución dinámica temprana (hasta $\sim 10$ millones de años), los resultados proporcionan un marco fundamental para interpretar los futuros datos de LISA y comprender el destino de los planetas en sistemas binarios evolucionados. Declaran explícitamente que sus afirmaciones de detectabilidad son ilustrativas y que se requieren análisis de forma de onda dedicados para confirmar las perspectivas observacionales.