Cepheids with giant companions III. Evolutionary modeling of nine binary double Cepheids from the Milky Way and Magellanic Clouds

Resumen Técnico: Cefeidas con Compañeras Gigantes III. Modelado Evolutivo de Nueve Cefeidas Binarias Dobles

Problema y Contexto
Las cefeidas binarias dobles (BIND), sistemas que contienen dos variables cefeidas, ofrecen un laboratorio único para probar las teorías de evolución y pulsación estelar. Aunque actualmente se conocen diez de estos sistemas, solo uno (LMC-CEP-1718) ha tenido sus parámetros determinados completamente mediante el modelado de binarias eclipsantes. Existe una brecha observacional significativa respecto a los parámetros físicos (masa, radio, luminosidad) de las cefeidas en la Nube Pequeña de Magallanes (SMC), y la relación masa-luminosidad (ML) permanece insuficientemente restringida, particularmente en el extremo de menor masa (las predicciones teóricas abarcan de 3 a 11 $M_\odot$, mientras que las observaciones se limitan a 3.6–5 $M_\odot$). Además, muchos sistemas BIND exhiben diferencias inusualmente grandes en los períodos de pulsación entre sus componentes, lo que desafía los modelos evolutivos coetáneos estándar y sugiere posibles interacciones binarias, como la transferencia de masa o las fusiones.

Metodología
Los autores extienden el método q-PED (introducido en Espinoza-Arancibia & Pilecki 2025) para analizar nueve nuevas cefeidas BIND identificadas en la Vía Láctea (MW), la Gran Nube de Magallanes (LMC) y la Nube Pequeña de Magallana (SMC). Dado que las razones de masa dinámica directa ($q_s$) no están disponibles para la mayoría de estos sistemas, el método utiliza los períodos de pulsación de ambos componentes como aproximaciones para la razón de masa.

El enfoque integra:

1. Modelos Evolutivos: Una densa cuadrícula de trayectorias computadas utilizando MESA (versión r22.1.1) a través de metalicidades consistentes con la MW ($Z=0.009-0.011$), LMC ($Z=0.006-0.009$) y SMC ($Z=0.001-0.003$).
2. Modelos de Pulsación: Modelos de pulsación estelar radial lineal (RSP) generados mediante la funcionalidad RSP de MESA, restringidos al régimen lineal. Los modelos se seleccionan dentro de los límites empíricos de la franja de inestabilidad (IS) para cada galaxia, coincidiendo con los períodos observados dentro de un margen del 5%.
3. Restricciones: El método emplea fotometría multibanda y módulos de distancia conocidos (Pietrzyński et al. 2019 para LMC; Graczyk et al. 2020 para SMC) para restringir las magnitudes del sistema.
4. Análisis de Configuración: Los autores consideran todas las configuraciones evolutivas consistentes: primer cruce (1C), bucle azul (BL) y escenarios mixtos. Para los sistemas con razones de masa espectroscópicas ($q_s$), estos valores se utilizan para discriminar entre soluciones viables.

Resultados Clave

• Estimación de Parámetros: El estudio proporciona nuevas estimaciones de masa, radio, temperatura, luminosidad y edad para 18 cefeidas con parámetros físicos previamente desconocidos. Esto extiende el rango de masa de las cefeidas caracterizadas hasta los 2.3 $M_\odot$ en la SMC.
• Validación de Referencia: Para el sistema eclipsante LMC-CEP-1718, el método q-PED reproduce los parámetros fundamentales con un error de pocos puntos porcentuales. Sin embargo, el modelo predice una razón de masa ligeramente invertida en comparación con las observaciones espectroscópicas ($q_s = 0.84 \pm 0.04$), lo que sugiere que la primaria pudo haber experimentado una pérdida de masa impulsada por la pulsación, revirtiendo la razón de masa.
• Configuraciones Evolutivas:
  • Sistemas No Interactuantes: Sistemas como LMC-CEP-0571 y SMC-CEP-1526 (donde $q_s \approx 1.00$) se describen mejor como cefeidas BL coetáneas que evolucionan de forma independiente.
  • Sistemas Interactuantes: La razón de masa espectroscópica para GD-CEP-291 ($q_s = 0.84 \pm 0.04$) identifica de forma única una configuración 1C + BL, lo que implica que la primaria es producto de una interacción binaria (probablemente una fusión).
  • Razones de Masa: El estudio establece una razón de masa teórica mínima para las cefeidas BIND no interactuantes basada en isócronas y cruces de la IS. Para la metalicidad de la SMC ($Z=0.002$), este mínimo cae a $q \approx 0.936$, mientras que para LMC y MW es de $\approx 0.95$. Las razones observadas significativamente por debajo de estos umbrales (por ejemplo, $q < 0.85$) sugieren fuertemente interacciones binarias pasadas.
• Relaciones:
  • Se deriva una nueva relación Período–Masa–Radio (PMR): $\log P = -1.675 - 0.653 \log M + 1.741 \log R$. Esta relación produce residuales de período de menos del 3% para la muestra.
  • Se ajusta una nueva relación Masa–Luminosidad (ML) para cefeidas BL ($\log M = 2.922 \log P + 1.285$), la cual separa exitosamente a las estrellas en el segundo y tercer cruces de la IS. También se proporciona una relación ML preliminar para cefeidas 1C, la cual se sitúa significativamente por debajo de la relación BL.

Significancia y Conclusiones
Este trabajo representa las primeras estimaciones de masa para las cefeidas en la SMC, expandiendo significativamente el espacio de parámetros de las masas conocidas de las cefeidas. Los autores concluyen que las interacciones binarias en la evolución pasada de las cefeidas pueden ser comunes, afectando potencialmente hasta al 40% de los sistemas estudiados (dos casos confirmados, más dos más si se prefieren las soluciones BL específicas). La identificación de GD-CEP-291 como un sistema 1C+BL respalda la hipótesis de que las cefeidas de primer cruce en binarias pueden originarse frecuentemente a partir de fusiones.

El artículo enfatiza que, si bien el método q-PED proporciona parámetros preliminares robustos, la historia evolutiva definitiva de estos sistemas depende de la mejora de las razones de masa espectroscópicas de las campañas observacionales en curso. El estudio subraya la necesidad de combinar el modelado teórico con restricciones empíricas para resolver la degeneración entre los diferentes cruces evolutivos e historias de interacción.