X-Shooting ULLYSES: Massive stars at low metallicity XIV. Properties of SMC late-O and B supergiants reveal the metallicity dependence of winds in the Magellanic Clouds

Resumen Técnico: Propiedades de supergigantes tardías-O y B en la Nube Menor de Magallanes (SMC) y la dependencia de la metalicidad en los vientos estelares

1. Problema y Contexto

Las estrellas masivas ($M > 8 M_\odot$) pierden masa a través de vientos impulsados por radiación, un proceso crucial que determina su evolución, su destino final (supernovas, agujeros negros o estrellas de neutrones) y su retroalimentación química y mecánica en el medio interestelar. La teoría de vientos impulsados por líneas predice que las propiedades de estos vientos dependen fuertemente de la metalicidad ($Z$) del entorno, ya que los metales (especialmente especies tipo hierro) proporcionan la opacidad necesaria para transferir el momento de la radiación al gas.

Sin embargo, existe una discrepancia persistente entre las predicciones teóricas y las derivaciones empíricas de las tasas de pérdida de masa en entornos de baja metalicidad. Además, se debate la existencia de un "salto de bi-estabilidad" (un aumento abrupto en la pérdida de masa y una disminución en la velocidad terminal del viento alrededor de 21-25 kK) en estrellas B. Este estudio busca cuantificar la dependencia de la metalicidad en los vientos de supergigantes azules en la Nube Menor de Magallanes (SMC, $Z \approx 0.2 Z_\odot$) y compararlos con estudios previos en la Nube Mayor de Magallanes (LMC, $Z \approx 0.5 Z_\odot$) para validar las recetas utilizadas en modelos evolutivos.

2. Metodología

El estudio se basa en el análisis espectroscópico de una muestra de 20 supergigantes tardías-O y B en la SMC.

• Datos Observacionales: Se utilizaron datos combinados del programa ULLYSES (HST, espectroscopía UV de alta resolución) y el legado XShootU (VLT/X-shooter, espectroscopía óptica de 3000-10200 Å). También se incluyeron datos de FUSE (UV lejano) y MIKE (Magellan, alta resolución óptica) para estrellas específicas.
• Modelado Atmosférico: Se empleó el código de transferencia radiativa CMFGEN (en equilibrio no-LTE y simetría esférica) para modelar las atmósferas estelares y los vientos. El modelo considera el "apantallamiento" de líneas de hierro y la inhomogeneidad del viento (clumping) bajo la aproximación de clumps ópticamente delgados.
• Procedimiento de Ajuste:
  • Se determinaron parámetros estelares ($T_{\text{eff}}$, $\log g$, abundancia de helio) mediante el equilibrio de ionización y el ensanchamiento de las líneas de Balmer.
  • Se ajustaron los parámetros del viento ($\dot{M}$, velocidad terminal $v_\infty$, parámetro de aceleración $\beta$, factor de llenado volumétrico $f_{\text{vol},\infty}$) utilizando líneas de emisión/absorción como H$\alpha$ y perfiles P Cygni en el UV (Si IV, C IV, etc.).
  • Se utilizó un método de inferencia bayesiana (similar a BONNSAI) para estimar las masas evolutivas y edades, comparándolas con las masas espectroscópicas derivadas de la gravedad superficial.
• Comparación: Los resultados se compararon directamente con un estudio previo de la LMC (Paper XIII) realizado con las mismas técnicas y códigos, permitiendo un análisis controlado de la variable metalicidad.

3. Contribuciones Clave

• Homogeneidad: Es uno de los primeros estudios que analiza una muestra amplia de supergigantes en la SMC utilizando un único código de modelado (CMFGEN) y un conjunto de datos UV/óptico coherente (ULLYSES/XShootU), eliminando la variabilidad introducida por el uso de diferentes códigos y conjuntos de datos en la literatura.
• Receta Empírica de Momento de Viento: Derivación de una nueva receta matemática que describe el momento del viento en función de la luminosidad y la metalicidad para la SMC.
• Análisis de la Bi-estabilidad: Una investigación exhaustiva sobre la existencia del salto de bi-estabilidad en la SMC, utilizando datos de alta calidad que cubren el rango de temperatura crítico (21-25 kK).
• Implicaciones Evolutivas: Evaluación crítica de la capacidad de los vientos de supergigantes azules para despojar las capas de hidrógeno y formar estrellas Wolf-Rayet (WR) en entornos de baja metalicidad.

4. Resultados Principales

• Dependencia de la Metalicidad: Se confirma una fuerte dependencia de la metalicidad en el momento del viento ($D_{\text{mom}}$). La relación derivada es:
  $$\log D_{\text{mom}} = (1.64 - 0.75 \log Z/Z_\odot) \log (L_{\text{bol}}/10^6 L_\odot) + 1.38 \log Z/Z_\odot + 29.17$$
  Esta dependencia se debe principalmente a las tasas de pérdida de masa ($\dot{M}$), que son significativamente más bajas en la SMC que en la LMC, mientras que la velocidad terminal ($v_\infty$) no muestra una dependencia clara con la metalicidad en este rango.

• Tasas de Pérdida de Masa ($\dot{M}$):

  • Las tasas derivadas son generalmente más bajas que las predichas por recetas numéricas recientes (como Vink & Sander 2021 o Björklund et al. 2023) para estrellas B frías ($T_{\text{eff}} < 22$ kK).
  • Las recetas numéricas sobreestiman $\dot{M}$ en el régimen de estrellas B, mientras que coinciden mejor con las estrellas O.
  • No se encontró evidencia de un aumento abrupto en $\dot{M}$ (salto de bi-estabilidad) entre 25 y 21 kK. El comportamiento de $\dot{M}$ es continuo o incluso decreciente con la temperatura.

• Velocidad Terminal ($v_\infty$): La relación entre la velocidad terminal y la velocidad de escape ($v_\infty/v_{\text{esc}}$) disminuye monótonamente con la temperatura, sin mostrar la caída característica esperada en el salto de bi-estabilidad.

• Discrepancia de Masas: Se encontró una discrepancia significativa (40% de la muestra) entre la masa evolutiva y la masa espectroscópica, siendo la masa evolutiva sistemáticamente mayor. Esto sugiere que los modelos actuales de pérdida de masa podrían no estar capturando completamente la física de la evolución de estas estrellas.

• Composición Química: Se observó un enriquecimiento significativo de nitrógeno en la mayoría de las estrellas, lo que no puede explicarse únicamente por la rotación inicial alta (dado que la mayoría son rotadores lentos). Esto apunta a la interacción binaria o a una mezcla rotacional más eficiente de lo previsto.

5. Significado y Conclusiones

El estudio tiene implicaciones profundas para la evolución estelar y la formación de estrellas Wolf-Rayet (WR):

1. Formación de Estrellas WR en la SMC: Las tasas de pérdida de masa actuales de las supergigantes azules en la SMC son demasiado bajas para despojar la envoltura de hidrógeno de la estrella y formar una estrella WR clásica a través de la evolución de una sola estrella dentro de su vida restante.
2. Necesidad de Mecanismos Alternativos: Para explicar la población observada de estrellas WR en la SMC (muchas de las cuales parecen ser de una sola estrella), es necesario recurrir a mecanismos adicionales, como:
   • Erupciones masivas tipo Variable Luminosa Azul (LBV) que provoquen una pérdida de masa episódica extrema.
   • Interacciones binarias (transferencia de masa, envoltura común) que eliminen la envoltura de hidrógeno.
3. Validación de Modelos: Los resultados desafían las recetas de pérdida de masa estándar utilizadas en los códigos de evolución estelar, sugiriendo que estas necesitan ser recalibradas para entornos de baja metalicidad, especialmente para estrellas B.

En resumen, este trabajo proporciona una base empírica sólida y homogénea para entender cómo la metalicidad moldea los vientos estelares, revelando que la física de la pérdida de masa en la SMC es más compleja y menos eficiente que lo predicho por los modelos teóricos actuales, lo que obliga a reconsiderar los caminos evolutivos hacia las estrellas Wolf-Rayet en galaxias de baja metalicidad.