Multivariate Statistical Analysis of Low Mass Ratio Contact Binaries: Definition, Dynamical Stability, and Parameter Relationships

Resumen Técnico: Análisis Estadístico Multivariante de Binarias de Contacto de Baja Relación de Masas

Planteamiento del Problema
Las binarias de contacto, específicamente los sistemas de tipo W Ursae Majoris (W UMa), son estrellas binarias cercanas que comparten una envoltura convectiva común. Un parámetro crítico que gobierna su evolución y estabilidad es la relación de masas ($q = M_2/M_1$). Mientras que los modelos teóricos sugieren que estos sistemas evolucionan hacia relaciones de masa más bajas, alcanzando potencialmente una inestabilidad de Darwin por marea que conduce a una fusión, el límite empírico preciso que separa los sistemas de "baja relación de masa" de los sistemas estables de mayor relación sigue siendo objeto de debate. Los puntos de corte teóricos varían ampliamente ($q \sim 0.04$ a $0.10$), y las restricciones observacionales suelen estar limitadas por muestras pequeñas, efectos de selección y la falta de un análisis uniforme en conjuntos de datos homogéneos. Además, la suposición de un radio de giro al cuadrado ($k^2$) constante para todos los componentes puede ser inexacta para sistemas de baja relación de masa donde las estrellas secundarias exhiben estructuras internas distintas (por ejemplo, totalmente convectivas).

Metodología
El estudio emplea un enfoque estadístico de múltiples etapas utilizando una muestra de 818 binarias de contacto y un subconjunto dedicado de 115 sistemas de baja relación de masa:

1. Definición del Rango de Baja Relación de Masas:

   • Los autores analizaron 818 binarias de contacto utilizando parámetros independientes: periodo orbital ($P$), relación de masa ($q$), factor de llenado ($f$) y temperaturas de los componentes ($T_{1,2}$).
   • Se adoptó inicialmente un umbral de $q = 0.25$ para comparar las distribuciones. Se utilizaron pruebas estadísticas (Mann-Whitney U) y análisis de árboles de decisión para evaluar la importancia de las características.
   • Para definir un límite empírico robusto, se aplicó una Estimación de Densidad de Kernel (KDE) tridimensional al espacio de parámetros conjunto de $\log P$, $f$ y $q$. Este método identificó caídas naturales de densidad en la distribución en lugar de imponer puntos de corte arbitrarios.

2. Análisis de Estabilidad Dinámica:

   • Se realizó un análisis de Monte Carlo (MC) sobre 501 sistemas para evaluar la relación entre el momento angular de rotación y el orbital ($J_{spin}/J_{orb}$).
   • A diferencia de estudios previos que asumen un $k^2$ fijo de $0.06$, este análisis muestreó radios de giro al cuadrado ($k^2_1$ y $k^2_2$) a partir de prioridades teóricas basadas en modelos de estructura estelar (por ejemplo, $k^2_2 = 0.18$ para secundarias de baja masa totalmente convectivas).
   • Las realizaciones de MC incorporaron incertidumbres observacionales y prioridades teóricas para calcular $J_{spin}/J_{orb}$ y evaluar la proximidad al límite de inestabilidad de Darwin ($0.333$).

3. Relaciones de Parámetros Empíricos:

   • Se compiló una muestra dedicada de 115 sistemas de baja relación de masa ($q < 0.27$) de la literatura.
   • Los parámetros absolutos se recalcularon utilizando paralajes de Gaia DR3. Los sistemas fueron filtrados basados en la extinción ($A_V \le 0.4$), el periodo orbital ($P < 0.6$ días) y la consistencia de los cálculos del semieje mayor.
   • Se realizaron regresiones lineales para derivar relaciones empíricas entre el periodo orbital y los parámetros absolutos (masa, radio, luminosidad, gravedad superficial, magnitud bolométrica), así como relaciones Masa-Luminosidad ($M-L$) y Masa-Radio ($M-R$).

Resultados Clave

• Umbral Empírico: El análisis de KDE reveló una caída significativa en la densidad del sistema en $q \approx 0.27$. Este valor se propone como el límite superior empírico para las binarias de contacto de baja relación de masa. Por debajo de este umbral, los sistemas exhiben factores de llenado significativamente más altos y periodos orbitales más cortos en comparación con los sistemas con $q \ge 0.25$.
• Importancia de las Características: El factor de llenado fue identificado como el discriminador más significativo entre los sistemas de baja y alta relación de masa (39.6% de importancia), seguido por el periodo orbital (18.7%).
• Estabilidad Rotacional: El análisis de MC demostró que, mientras el radio de giro al cuadrado de la primaria ($k^2_1$) permanece casi constante ($\approx 0.06$), el de la secundaria ($k^2_2$) y la relación $J_{spin}/J_{orb}$ disminuyen ligeramente a medida que la relación de masa aumenta. Esto confirma que la estructura interna de la estrella secundaria juega un papel crucial en la estabilidad rotacional, particularmente para relaciones de masa muy bajas donde la secundaria es totalmente convectiva.
• Relaciones Empíricas: El estudio proporciona un conjunto de ecuaciones empíricas lineales que relacionan el periodo orbital con los parámetros estelares para sistemas de baja relación de masa. Estas relaciones permiten la estimación de parámetros absolutos (masa, radio, luminosidad) sin requerir espectroscopia de alta resolución de ambos componentes.
• Estado Evolutivo: Al ser graficadas en diagramas de Masa-Luminosidad y Masa-Radio, las estrellas secundarias se encuentran por encima de la Secuencia Principal de Edad Terminal (TAMS), mientras que las estrellas primarias generalmente se sitúan cerca de la Secuencia Principal de Edad Cero (ZAMS), aunque algunas primarias también parecen evolucionadas. Los autores señan que las trayectorias evolutivas estándar de estrellas individuales deben interpretarse con cautela debido a las complejas interacciones en binarias de contacto.

Significación y Reivindicaciones
El artículo afirma proporcionar una definición empírica y estadísticamente robusta para las binarias de contacto de baja relación de masa ($q \lesssim 0.27$), yendo más allá de los puntos de corte teóricos arbitrarios. Al establecer este límite, el estudio ofrece un marco más claro para investigar la formación, evolución y destino final (fusión vs. estabilidad) de estos sistemas.

Los resúmenes estadísticos derivados de $k^2$ y $J_{spin}/J_{orb}$ sirven como referencias para futuros modelos, corrigiendo la suposición común de constantes estructurales constantes. Finalmente, el conjunto de datos compilado y las nuevas relaciones de parámetros empíricos proporcionan herramientas prácticas para estimar las propiedades físicas de las binarias de baja relación de masa, facilitando la investigación observacional y teórica adicional sobre sus vías evolutivas y su conexión con fenómenos como las estrellas de tipo azul errante (blue stragglers) y las estrellas de tipo FK Comae. Los autores enfatizan que estos hallazgos se basan en muestras observacionales actualizadas y un modelado consistente, abordando las brechas previas en el análisis uniforme.