Predicted incidence of Jupiter-like planets around white dwarfs

Resumen Técnico: Incidencia Predicha de Planetas tipo Júpiter alrededor de Enanas Blancas

Planteamiento del Problema
Aunque se han confirmado miles de exoplanetas alrededor de estrellas de la secuencia principal, y los modelos teóricos sugieren que los planetas con órbitas más allá de unas pocas unidades astronómicas deberían sobrevivir a la evolución de su estrella anfitje hacia una enana blanca (WD), la evidencia observacional sigue siendo escasa. A pesar de las investigaciones de imagen de alto contraste y las búsquedas de tránsitos, solo se ha detectado un puñado de compañeros subestelares en enanas blancas. Actualmente, no está claro si esta escasez refleja limitaciones observacionales o las consecuencias reales de la evolución estelar (por ejemplo, la absorción durante las fases gigantes, interacciones de marea o pérdida de masa). Los estudios previos de síntesis de población se han visto limitados por supuestos genéricos sobre las poblaciones planetarias y una falta de dependencia en parámetros específicos de las enanas blancas, tales como la edad y la masa. Este estudio tiene como objetivo llenar este vacío mediante la realización de una síntesis de población detallada para predecir la fracción y las propiedades de las enanas blancas que albergan compañeros subestelares.

Metodología
Los autores emplearon un enfoque de síntesis de población de múltiples etapas, combinando códigos de evolución estelar con algoritmos de dinámica orbital:

1. Generación de Progenitores: Se generó una población representativa de 243,617 potenciales progenitores de enanas blancas utilizando la función de masa inicial de Kroupa (0.8–4 $M_{\odot}$). Se asignaron edades estelares asumiendo una tasa de formación estelar constante para un disco delgado de 10 Gyr, mientras que las metalicidades se derivaron de una relación edad-metalicidad global (Carrillo et al. 2023).
2. Asignación de Compañeros: Se asignaron compañeros subestelares (0.8–30 $M_J$) a los progenitores basándose en tasas de ocurrencia derivadas de sondeos de velocidad radial de estrellas gigantes (Wolthoff et al. 2022). La tasa de ocurrencia se modeló como una función de la masa estelar y la metalicidad, con periodos orbitales siguiendo una distribución log-normal y excentricidades basadas en estadísticas recientes de exoplanetas (Kane & Wittenmyer 2024).
3. Evolución Estelar y Orbital:
   • Evolución Estelar: La mayoría de la muestra fue evolucionada utilizando el código Single Star Evolution (SSE) para determinar la población actual de enanas blancas. Un subconjunto de 1,000 sistemas que inicialmente albergaban compañeros fue re-simulado utilizando los Modules for Experiments in Stellar Astrophysics (MESA) para capturar la física detallada, particularmente los pulsos térmicos y la pérdida de masa en la Rama Gigante Asintótica (AGB).
   • Dinámica Orbital: Se utilizó el algoritmo FATES (Final exoplAnet orbiTal dEstination around evolved Stars) para integrar la evolución orbital de los compañeros. Este código contabiliza la expansión orbital debido a la pérdida isotrópica de masa estelar y el decaimiento orbital debido a las mareas estelares.
   • Pruebas de Incertidumbre: El estudio probó la sensibilidad a varias prescripciones físicas, incluyendo dos modelos diferentes de escala de tiempo de rotación convectiva para las fuerzas de marea (Rasio et al. 1996 [R96] vs. Villaver & Livio 2009 [VL09]), diferentes eficiencias de pérdida de masa en la AGB (Bloecker 1995 con $\eta=0.02$ y $\eta=0.7$), y relaciones de edad-metalicidad alternativas (Fernández-Alvar et al. 2025).

Resultados Clave
Las simulaciones arrojan las siguientes predicciones cuantitativas para la población actual de enanas blancas en la Vía Láctea:

• Tasa de Ocurrencia: La fracción predicha de enanas blancas que albergan compañeros subestelares es baja, estimada en menos del 3 ± 1.5 %. Esta cifra es robusta ante las incertidumbres en las distribuciones iniciales, las mareas estelares y la pérdida de masa en la AGB.
  • Bajo el supuesto de fuerzas de marea débiles (VL09), la tasa es de ~1.3%.
  • Bajo fuerzas de marea fuertes (R96), la tasa cae a ~0.6%.
  • Si se ignoran las fuerzas de marea por completo (solo pérdida de masa), la tasa aumenta a 3.2%, lo que representa un límite superior estricto.
• Naturaleza del Compañero: Se predice que aproximadamente el 95% de los compañeros supervivientes son planetas gigantes gaseosos (masas < 13 $M_J$). La tasa de ocurrencia general para planetas tipo Júpiter específicamente es inferior al 2.9 ± 1.4 %.
• Propiedades Orbitales: Se predice que los semiejes mayores de los compañeros supervivientes oscilan entre 3 y 24 au, con una mediana de 11 au. Las órbitas finales son generalmente más anchas que las iniciales debido a la pérdida de masa, aunque las fuerzas de marea pueden circularizar las órbitas excéntricas.
• Características de la Estrella Anfitje: La tasa de ocurrencia alcanza su máximo para enanas blancas de baja masa (0.53–0.66 $M_{\odot}$) y sistemas relativamente jóvenes (1–6 Gyr), donde la tasa puede superar el 3%.
• Dependencia de la Metalicidad y el Entorno Local: Las predicciones son altamente sensibles a la relación edad-metalicidad asumida. Mientras que la predicción global permanece por debajo del 3%, la aplicación de una relación edad-metalicidad local (Fernández-Alvar et al. 2025) sugiere que la fracción de enanas blancas con compañeros cerca del Sol podría alcanzar el < 8 %.

Significancia y Conclusiones
Los autores concluyen que, si las distribuciones de compañeros derivadas de los sondeos de velocidad radial de estrellas gigantes de masa intermedia reflejan fielmente la población real, menos del 3 ± 1.5 % de las enanas blancas albergan compañeros subestelares. Esto sugiere que tales compañeros son intrínsecamente raros.

El artículo enfatiza que esta baja fracción probablemente representa un límite superior de los compañeros detectables. Un número significativo de compañeros supervivientes puede permanecer indetectable con las instalaciones actuales debido a sus amplias separaciones y la debilidad de los sistemas. El estudio destaca que la supervivencia de los compañeros depende críticamente de la fuerza de las mareas estelares y la eficiencia de la pérdida de masa durante la fase AGB. Además, los resultados subrayan que la relación edad-metalicidad asumida influye significativamente en las predicciones, particularmente para muestras locales.

Los autores afirman que sus hallazgos proporcionan una base teórica necesaria para interpretar los sondeos observacionales actuales y futuros (por ejemplo, con JWST y el ELT). Señalan que, si bien los escenarios de sistemas multiplanetarios y de supervivencia de la envoltura común podrían aumentar estas fracciones, el modelo actual de compañero único proporciona una estimación de primer orden que indica que la escasez de compañeros de enanas blancas detectados es probablemente un fenómeno astrofísico real y no meramente un sesgo observacional.