MESA Isochrones and Stellar Tracks. II. Models with alpha-enhanced chemical composition

Resumen Técnico: MIST II - Modelos con Composición Química Enriquecida en $\alpha$

1. Problema y Motivación

Los modelos estelares son fundamentales para la astrofísica moderna, desde la interpretación de galaxias de alto corrimiento al rojo hasta la determinación de masas de exoplanetas. Sin embargo, la mayoría de los modelos anteriores asumen una composición química "escalada al solar", donde las abundancias relativas de los elementos se mantienen fijas respecto a la mezcla solar, variando solo la metalicidad total ($Z$) o el hierro ([Fe/H]).

Este enfoque es insuficiente porque:

• La nucleosíntesis en diferentes entornos (supernovas de colapso del núcleo, explosiones termonucleares de enanas blancas, quemado de capas en estrellas AGB) produce patrones de abundancia distintos.
• La desviación más común observada es en los elementos de captura $\alpha$ (O, Ne, Mg, Si, S, Ar, Ca, Ti) en relación con los elementos del pico de hierro.
• La mayoría de las estrellas en la Vía Láctea con [Fe/H] < -1, casi todos los cúmulos globulares y muchas galaxias masivas exhiben patrones de abundancia no escalados al solar (enriquecidos en $\alpha$).
• Existían pocas mallas de modelos que cubrieran el rango completo de masas (0.1 - 100 $M_\odot$), metalicidades y fases evolutivas variando simultáneamente [Fe/H] y [α/Fe].

2. Metodología

Los autores actualizan y expanden la base de datos MIST (MESA Isochrones and Stellar Tracks) utilizando el código de evolución estelar MESA (versión 11701).

• Química y Abundancias:

  • Se incluyen variaciones en [α/Fe] de -0.2 a +0.6 (en pasos de 0.2 dex) para un rango de [Fe/H] de -3 a +0.5.
  • Se adopta un enfoque donde cada elemento $\alpha$ se enriquece en la misma cantidad relativa a su abundancia solar escalada, igual a [α/Fe].
  • Se utiliza la mezcla solar Grevesse & Sauval (1998) (GS98) en lugar de la de Asplund et al. (2009) usada en el trabajo anterior, debido a nuevas tendencias en el análisis de abundancias solares y la disponibilidad de opacidades para composiciones GS98 enriquecidas.
  • Se define una relación de enriquecimiento de helio $\Delta Y/\Delta Z \approx 1.32$, basada en valores primordiales y la calibración solar inicial.

• Física Estelar y Microfísica:

  • Opacidades y Ecuación de Estado (EOS): Se utilizan tablas OPAL y se complementan con tablas FreeEOS para metalicidades altas ($Z \ge 0.04$). Se han generado tablas de opacidad adicionales para cada valor de [α/Fe] para capturar correctamente los efectos en la estructura estelar.
  • Condiciones de Frontera Atmosférica: Se emplea un método híbrido. Para estrellas masivas ($M \ge 8 M_\odot$) se usa una fotosfera simple; para estrellas menos masivas, se integra la relación $T-\tau$ (Temperatura-Profundidad Óptica) verificando la convección local. Se han creado tablas de condiciones de frontera específicas para cada composición química única.
  • Rotación: Se incluyen modelos con y sin rotación. Se incorpora transporte de momento angular debido a la mezcla Tayler-Spruit (prescripción de Fuller et al. 2019) para frenar los núcleos compactos en la rama de gigantes rojas (RGB), mejorando el acuerdo con datos sísmicos.
  • Mezcla y Sobrepaso (Overshoot): Se recalibra el parámetro de sobrepaso convectivo bajo la zona convectiva superficial ($f_{ov,env}$) utilizando datos del "bulto de la rama de gigantes rojas" (RGBB). La nueva calibración sigue la tendencia observada: $f_{ov,env} = 0.016 - 0.027 \times [Fe/H]$.
  • Espectros Sintéticos: Se generan espectros sintéticos utilizando ATLAS12 y SYNTHE con listas de líneas ajustadas empíricamente, incluyendo bases de datos como ExoMol para TiO, para convertir los modelos físicos al plano observacional.

3. Contribuciones Clave

• Expansión de la Malla de Modelos: Se ha creado una malla de 74 composiciones químicas únicas, cubriendo masas de 0.1 a 300 $M_\odot$ y todas las fases evolutivas (desde la pre-secuencia principal hasta el enfriamiento de enanas blancas o la quema central de carbono).
• Consistencia Física: Las variaciones en [α/Fe] se implementan de manera autoconsistente desde el interior estelar (opacidades, EOS, reacciones nucleares) hasta los espectros sintéticos y las correcciones bolométricas.
• Nuevas Columnas de Datos: Se han añadido datos sobre oscurecimiento gravitacional (para modelos rotantes), coordenadas de la zona de convección superficial y tiempos de recambio convectivo.
• Reproducibilidad: Se han liberado todos los archivos de entrada, directorios de trabajo y la instalación modificada de MESA necesaria para reproducir los modelos.

4. Resultados

• Impacto de [α/Fe]: La variación en [α/Fe] tiene efectos significativos en la evolución estelar, especialmente a metalicidades bajas. Aumenta la vida útil de las estrellas (principalmente debido al aumento de oxígeno) y modifica la opacidad en las capas superficiales.
• Incapacidad de la Escalación Simple: Se demuestra que simplemente escalar la metalicidad total ($Z$) de modelos con abundancias solares no puede replicar la morfología de las isocronas enriquecidas en $\alpha$, especialmente en edades avanzadas y en el diagrama Color-Magnitud (CMD).
• Validación Observacional:
  • Las nuevas isocronas reproducen correctamente la tendencia observada en la magnitud del bulto de la rama de gigantes rojas (RGBB) en función de la metalicidad, corrigiendo una sobreestimación de ~0.5 mag presente en la versión anterior (Paper 1) para cúmulos pobres en metales.
  • Se muestra un buen acuerdo con los cúmulos globulares galácticos (M92, M3, 47 Tuc, NGC 6397, NGC 6752) en diagramas CMD de HST (óptico e infrarrojo cercano), tanto en la rama de gigantes como en la secuencia principal baja.
  • Se observa que el enriquecimiento en $\alpha$ afecta fuertemente la formación de moléculas (H2O, TiO) en el infrarrojo y la absorción UV, alterando los colores estelares de manera distinta a la variación de [Fe/H].

5. Significancia

Este trabajo proporciona la herramienta más completa y actualizada disponible para estudiar poblaciones estelares con composiciones químicas no solares. Al ofrecer una malla que varía independientemente [Fe/H] y [α/Fe], permite:

• Interpretar con mayor precisión la edad y metalicidad de poblaciones estelares antiguas (cúmulos globulares, galaxias elípticas masivas).
• Mejorar la caracterización de exoplanetas alrededor de estrellas con composiciones químicas diversas.
• Entender la evolución química del universo, desde las primeras estrellas hasta las galaxias masivas actuales.
• Los datos están disponibles públicamente en el sitio web del proyecto MIST y en Zenodo, facilitando su uso en la comunidad astrofísica.