Wolf-Rayet stellar evolution models with improved treatment of the atmosphere

Resumen Técnico: Modelos de Evolución Estelar de Wolf-Rayet con un Tratamiento Mejorado de la Atmósfera

Enunciado del Problema
Los modelos evolutivos de estrellas masivas han dependido históricamente de la aproximación de atmósfera gris de Eddington para definir las condiciones de contorno externas. Esta aproximación falla al capturar fenómenos físicos complejos inherentes a las atmósferas de estrellas masivas, como los efectos de no equilibrio termodinámico local (no-LTE), la transferencia radiativa no gris, los vientos estelares y el oscurecimiento por líneas. En consecuencia, los modelos estándar predicen temperaturas efectivas ($T_{\text{eff}}$) para las estrellas Wolf-Rayet (WR) que son significativamente más altas que las derivadas de análisis espectrales utilizando códigos de atmósfera detallados. Esta discrepancia dificulta la colocación precisa de estrellas masivas evolucionadas en el diagrama Hertzsprung-Russell (HR) y limita la comprensión cuantitativa de su papel en la evolución galáctica.

Metodología
Los autores desarrollaron un método para integrar atmósferas modelo de vanguardia directamente en el código de evolución estelar STAREVOL (versión 3.40) para estrellas masivas. La estrategia central consiste en eludir el acoplamiento en tiempo real, que consume mucho tiempo, entre los códigos de evolución y de transferencia radiativa, mediante el uso de una cuadrícula precalculada de atmósferas.

1. Construcción de la Cuadrícula de Atmósferas: Se calculó una cuadrícula de atmósferas modelo no-LTE, simétricas esféricamente, utilizando CMFGEN. Los parámetros de la cuadrícula incluyeron luminosidad ($L$), temperatura efectiva ($T_{\text{eff}}$), gravedad superficial ($g$), tasa de pérdida de masa ($\dot{M}$) y abundancias químicas superficiales.
2. Esquema de Interpolación: Para proporcionar condiciones de contorno durante la evolución, los autores emplearon un método de interpolación basado en la triangulación de Delaunay (utilizando el algoritmo delaunaysparse). La interpolación se realizó en un espacio de parámetros normalizado definido por $T(\tau_{\text{match}})$ (temperatura en una profundidad óptica de coincidencia $\tau_{\text{match}} = 25$), $\log(g)$ y $\log(\dot{M})$. La luminosidad no se utilizó como variable de interpolación, sino que se corrigió post-interpolación para coincidir con la luminosidad del modelo de evolución.
3. Acoplamiento a la Evolución:
   • La superficie numérica del modelo de evolución se estableció en $\tau_0 = 10$, lo suficientemente profundo para garantizar la validez de la aproximación de difusión.
   • En cada paso de tiempo, el código interpoló los perfiles de temperatura, densidad y radio desde la cuadrícula de CMFGEN basándose en los parámetros estelares actuales.
   • Estos perfiles interpolados sirvieron como condiciones de contorno para las ecuaciones de la estructura interna (específicamente temperatura y densidad en $\tau_0$).
   • El perfil de radio interpolado se utilizó para definir el radio fotosférico ($R_{\text{eff}}$) y calcular la temperatura efectiva ($T_{\text{eff}}$) y la gravedad superficial para la siguiente iteración.
4. Parámetros del Modelo: El método se aplicó a modelos no rotatorios con masas iniciales de 50 $M_\odot$ y 60 $M_\odot$ a metalicidad solar. Las prescripciones de pérdida de masa se adaptaron para diferentes fases (OB, LBV, WR), con una modificación específica para la fase WR para asegurar que pudieran calcularse estructuras de atmósfera hidrodinámicamente consistentes cerca del límite de Eddington.

Resultados Clave

• Reducción de la Temperatura Efectiva: La inclusión de atmósferas detalladas reduce significativamente las temperaturas efectivas predichas durante la fase WR. En comparación con los modelos grises clásicos de Eddington, los modelos basados en CMFGEN muestran reducciones de $T_{\text{eff}}$ que superan los 20 kK durante aproximadamente la mitad del tiempo de evolución WR.
• Expansión Fotosférica: La reducción en $T_{\text{eff}}$ se impulsa principalmente por una expansión del radio fotosférico (de hasta ~80% en casos específicos) causada por el viento estelar, una característica ausente en las atmósferas grises planas y estáticas.
• Estructura Interna Inalterada: A pesar de los cambios significativos en las propiedades superficiales, la estructura interna (perfiles de temperatura, densidad y radio a profundidades ópticas altas) y la evolución química de los modelos permanecen virtualmente idénticos a los cálculos clásicos. La retroalimentación desde la atmósfera hacia el interior profundo es despreciable.
• Acuerdo con las Observaciones: Las trayectorias revisadas muestran un acuerdo mucho mejor con las observaciones de estrellas WR de la Vía Láctea y de la Gran Nube de Magallanes (LMC) (tanto tipos WN como WC). Específicamente, los modelos reproducen con éxito la posición de las estrellas WNh y la transición a estrellas WN y WC libres de hidrógeno a temperaturas efectivas consistentes con los datos observacionales (por ejemplo, estrellas WC apareciendo alrededor de 85,000 K).
• Equivalencia en el Post-Procesamiento: Los autores demostraron que aplicar una corrección de post-procesamiento a los modelos clásicos, utilizando el mismo método de interpolación para recalcular $T_{\text{eff}}$ y $R_{\text{eff}}$ basándose en los parámetros internos del modelo clásico, produce trayectorias evolutivas casi indistinguibles de las calculadas con la inclusión directa de atmósferas. Esto confirma que el impacto principal recae en la definición de la superficie y no en la evolución interna.

Significado y Afirmaciones
El artículo afirma que la inclusión de modelos de atmósfera detallados es necesaria para reproducir correctamente la posición de las estrellas masivas evolucionadas en el diagrama HR, resolviendo la discrepancia de larga data entre la $T_{\text{eff}}$ teórica y el análisis espectral observacional. Sin embargo, los autores son modestos respecto al alcance de esta mejora:

• El método corrige la apariencia (propiedades superficiales) de los modelos, pero no altera la evolución interna o química.
• El estudio se limita a modelos no rotatorios de metalicidad solar de 50 y 60 $M_\odot$.
• Las tasas de pérdida de masa utilizadas se redujeron en un factor de 3 para garantizar la estabilidad numérica en los cálculos de atmósfera; por lo tanto, las reducciones de temperatura reportadas representan un límite inferior del efecto esperado para las fases WR libres de hidrógeno.
• Los autores reconocen limitaciones en el tratamiento de las zonas de convección subsuperficiales (por ejemplo, el pico de opacidad del hierro) y la física específica del mecanismo de lanzamiento del viento en WNh, señalando que se requiere más trabajo para lograr un tratamiento completamente autoconsistente a través de todas las fases evolutivas.

En conclusión, el estudio establece que, si bien las atmósferas detalladas son cruciales para definir las propiedades observables de las estrellas WR, la evolución interna de estas estrellas es robusta frente a la elección de las condiciones de contorno atmosféricas, validando el uso de correcciones de post-procesamiento como una alternativa eficiente al acoplamiento directo en muchos contextos.