Detailed Investigation of a W UMa Contact Binary with an Ultralow Mass Ratio and a Third-Body as a Potential Merger Candidate

Resumen Técnico: Investigación Detallada de una Binaria de Contacto W UMa con una Relación de Masas Ultrabaja y un Tercer Cuerpo como Candidata Potencial a Fusión

Planteamiento del Problema
El límite inferior de la relación de masas ($q = M_2/M_1$) en los sistemas binarios de contacto W UMa sigue siendo un tema de incertidumbre, con evidencia observacional que sugiere la existencia de sistemas por debajo de las predicciones de estabilidad teórica. La estabilidad dinámica de tales sistemas de relación de masas extremadamente baja ($q < 0.1$) es objeto de debate, ya que son potenciales progenitores de fusiones estelares. Estudios previos de sistemas como VSX J082700.8+462850 y TYC 4002-2628-1 han dependido en gran medida del modelado fotométrico, que a menudo carece de confirmación espectroscópica del componente secundario debido a su intrínseca debilidad. Un estudio reciente realizado por Li et al. (2024) analizó el objetivo TYC 3801-1529-1, reportando una relación de masas de $q \approx 0.0356$. Sin embargo, los autores de este estudio señalaron discrepancias en el ajuste de la curva de luz y buscaron reevaluar los parámetros del sistema, las variaciones del periodo orbital y el estado de estabilidad utilizando datos de TESS y técnicas de modelado avanzado.

Metodología
Los autores realizaron un reanálisis exhaustivo de TYC 3801-1529-1 (TIC 445065429) utilizando los siguientes enfoques:

1. Modelado Fotométrico: Los datos de la curva de luz de TESS de los sectores 20, 47, 60 y 74 fueron procesados utilizando el Quick-Look Pipeline (QLP). El análisis empleó el código PHOEBE 2.4.9 para la validación inicial y la aplicación BSN (un acelerador MCMC basado en Windows) para la optimización final. La aplicación BSN, configurada con 28 caminantes (walkers) y 1200 iteraciones, se utilizó para derivar las incertidumbres y distribuciones de los parámetros, muestreando cinco parámetros principales: temperaturas efectivas ($T_1, T_2$), relación de masas ($q$), factor de llenado (fillout factor, $f$) e inclinación orbital ($i$).
2. Estimación de Parámetros Absolutos: Utilizando la relación de masas derivada y asumiendo la masa primaria consistente con Li et al. (2024), los autores calcularon la masa secundaria, el semieje mayor (vía la tercera ley de Kepler), los radios estelares y las luminosidades.
3. Análisis de la Variación del Periodo Orbital (O-C): Los autores compilaron 408 tiempos de mínimo de TESS, uno de ZTF y 17 de observaciones terrestres (incluyendo nuevos datos del Observatorio de San Pedro Mártir), cubriendo un lapso de seis años. Aplicaron el método de Observado menos Calculado (O-C) para investigar los cambios de periodo.
4. Análisis de Tercer Cuerpo y de Aliasing: Para distinguir entre cambios de periodo intrínsecos e influencias externas, los autores modelaron los residuos O-C mediante una función sinusoidal para probar el efecto de tiempo de viaje de la luz de un tercer cuerpo. También realizaron un análisis de periodograma de Lomb-Scargle y pruebas de preblanqueamiento (prewhitening) para descartar efectos de aliasing causados por patrones de muestreo estacionales.
5. Análisis de Estabilidad y de Derivadas: El estudio calculó la relación entre el momento angular de rotación y el orbital ($J_{spin}/J_{orb}$) utilizando radios de giro ($k^2$) actualizados y apropiados para las estructuras estelares específicas. Además, se aplicó el método de la derivada de Kouzuma (2023, 2025) a la curva de luz para estimar de forma independiente la relación de masas.

Resultados Clave

• Relación de Masas Ultrabaja: El reanálisis arrojó una relación de masas de $q = 0.024^{+0.002}_{-0.001}$, significativamente menor que el $0.0356$ reportado por Li et al. (2024). Esta representa la relación de masas determinada experimentalmente más baja para una binaria de contacto hasta la fecha.
• Propiedades de los Componentes: Se estima que el componente secundario tiene una masa de $M_2 \approx 0.050 M_\odot$, situándolo dentro del régimen de enanas marrones ($0.013–0.08 M_\odot$). La masa primaria es $M_1 \approx 2.096 M_\odot$. El sistema se clasifica como una binaria de contacto de subtipo A con un factor de llenado medio ($f \approx 0.315$).
• Presencia de un Tercer Cuerpo: El diagrama O-C exhibe una clara variación cíclica. Tras eliminar una tendencia lineal, los residuos se ajustaron mejor a un modelo sinusoidal con un periodo de $P_{mod} \approx 8.75$ años. La función de masa derivada y las simulaciones de Monte Carlo sugieren un tercer cuerpo con una masa de $M_3 \approx 0.49 M_\odot$ (probablemente una enana de tipo M de baja masa) que contribuye con una luz insignificante ($l_3 \approx 0.001$).
• Variaciones de Periodo: Una vez eliminada la contribución del tercer cuerpo, las variaciones intrínsecas del periodo de la binaria pudieron describirse mediante una tendencia lineal con pendiente negativa o una tendencia cuadrática con curvatura descendente. Las métricas estadísticas (AIC/BIC) favorecieron el modelo lineal más simple, aunque los autores señalan que un ajuste cuadrático implicaría un periodo orbital decreciente.
• Inestabilidad Dinámica: La relación calculada entre el momento angular de rotación y el orbital es $J_{spin}/J_{orb} \approx 1.18$. Este valor excede significativamente el umbral de estabilidad de Darwin ($\sim 0.33$), lo que indica que el sistema es dinámicamente inestable.
• Validación de la Relación de Masas: La relación de masas derivada mediante el método de la derivada de Kouzuma ($q \approx 0.031$) resultó ser consistente con los resultados de MCMC y con Li et al. (2024), reforzando la conclusión de que el sistema posee una relación de masas extremadamente baja.

Significancia y Reivindicaciones
El artículo sostiene que TYC 3801-1529-1 representa un candidato notable para una fusión binaria debido a su configuración dinámicamente inestable y su relación de masas ultrabaja. El estudio destaca los desafíos en la detección de segundas componentes tenues en tales sistemas, que a menudo conducen a degeneraciones en las soluciones fotométricas. Al demostrar que las variaciones de periodo del sistema son probablemente impulsadas por un compañero terciario en lugar de la actividad magnética (el mecanismo de Applegate fue descartado debido a la insuficiencia de la variación del momento cuadrupolar), los autores proporcionan una imagen más clara de la arquitectura del sistema.

Los autores enfatizan que, si bien el sistema es un fuerte candidato a fusión, la dirección definitiva del cambio del periodo orbital (creciente vs. decreciente) sigue siendo ambigua sin más observaciones a largo plazo. El estudio subraya la necesidad de un seguimiento espectroscópico de alta resolución para confirmar la relación de masas y la naturaleza de la secundaria, que actualmente se infiere que es una enana marrón. Los hallazgos contribuyen al debate continuo sobre la relación de masas mínima para las binarias de contacto y las vías evolutivas que conducen a las fusiones estelares.