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

Resumo Técnico: Modelos de Evolução Estelar de Wolf-Rayet com Tratamento Aprimorado da Atmosfera

Declaração do Problema
Os modelos evolutivos de estrelas massivas historicamente dependeram da aproximação de atmosfera cinzenta de Eddington para definir as condições de contorno externas. Esta aproximação falha em capturar fenômenos físicos complexos inerentes às atmosferas de estrelas massivas, como efeitos de não equilíbrio termodinâmico local (não-LTE), transferência radiativa não cinzenta, ventos estelares e obscurecimento por linhas (line-blanketing). Consequentemente, os modelos padrão preveem temperaturas efetivas ($T_{\text{eff}}$) para estrelas Wolf-Rayet (WR) que são significativamente mais altas do que aquelas derivadas de análises espectrais utilizando códigos de atmosfera detalhados. Esta discrepância impede a colocação precisa de estrelas massivas evoluídas no diagrama de Hertzsprung-Russell (HR) e limita a compreensão quantitativa de seu papel na evolução galáctica.

Metodologia
Os autores desenvolveram um método para integrar atmosferas de modelo de última geração diretamente no código de evolução estelar STAREVOL (versão 3.40) para estrelas massivas. A estratégia central envolve contornar o acoplamento em tempo real, que consome tempo, entre códigos de evolução e de transferência radiativa, utilizando uma grade pré-calculada de atmosferas.

1. Construção da Grade de Atmosferas: Uma grade de atmosferas de modelo esféricas e simétricas em não-LTE foi calculada usando o CMFGEN. Os parâmetros da grade incluíram luminosidade ($L$), temperatura efetiva ($T_{\text{eff}}$), gravidade superficial ($g$), taxa de perda de massa ($\dot{M}$) e abundâncias químicas superficiais.
2. Esquema de Interpolação: Para fornecer condições de contorno durante a evolução, os autores empregaram um método de interpolação baseado em triangulação de Delaunay (utilizando o algoritmo delaunaysparse). A interpolação foi realizada em um espaço de parâmetros normalizado definido por $T(\tau_{\text{match}})$ (temperatura em uma profundidade óptica de correspondência $\tau_{\text{match}} = 25$), $\log(g)$ e $\log(\dot{M})$. A luminosidade não foi utilizada como variável de interpolação, mas foi corrigida pós-interpolação para corresponder à luminosidade do modelo de evolução.
3. Acoplamento à Evolução:
   • A superfície numérica do modelo de evolução foi definida em $\tau_0 = 10$, suficientemente profunda para garantir a validade da aproximação de difusão.
   • Em cada passo de tempo, o código interpolou os perfis de temperatura, densidade e raio da grade do CMFGEN com base nos parâmetros estelares atuais.
   • Esses perfis interpolados serviram como condições de contorno para as equações da estrutura interna (especificamente temperatura e densidade em $\tau_0$).
   • O perfil de raio interpolado foi utilizado para definir o raio fotosférico ($R_{\text{eff}}$) e calcular a temperatura efetiva ($T_{\text{eff}}$) e a gravidade superficial para a próxima iteração.
4. Parâmetros do Modelo: O método foi aplicado a modelos não rotativos com massas iniciais de 50 $M_\odot$ e 60 $M_\odot$ em metalicidade solar. As prescrições de perda de massa foram adaptadas para diferentes fases (OB, LBV, WR), com uma modificação específica para a fase WR para garantir que estruturas de atmosfera consistentes hidrodinamicamente pudessem ser calculadas próximas ao limite de Eddington.

Principais Resultados

• Redução da Temperatura Efetiva: A inclusão de atmosferas detalhadas reduz significativamente as temperaturas efetivas previstas durante a fase WR. Comparados aos modelos cinzentos clássicos de Eddington, os modelos baseados no CMFGEN mostram reduções de $T_{\text{eff}}$ superiores a 20 kK para aproximadamente metade do tempo de evolução WR.
• Expansão Fotosférica: A redução em $T_{\text{eff}}$ é impulsionada principalmente por uma expansão do raio fotosférico (de até ~80% em casos específicos) causada pelo vento estelar, uma característica ausente em atmosferas cinzentas planas e estáticas.
• Estrutura Interna Inalterada: Apesar das mudanças significativas nas propriedades superficiais, a estrutura interna (perfis de temperatura, densidade e raio em altas profundidades ópticas) e a evolução química dos modelos permanecem virtualmente idênticas aos cálculos clássicos. O feedback da atmosfera para o interior profundo é negligenciável.
• Concordância com Observações: As trajetórias revisadas mostram uma concordância muito melhor com as observações de estrelas WR da Via Láctea e da Grande Nuvem de Magalhães (LMC) (ambos os tipos WN e WC). Especificamente, os modelos reproduzem com sucesso a posição das estrelas WNh e a transição para estrelas WN e WC livres de hidrogênio em temperaturas efetivas consistentes com dados observacionais (por exemplo, estrelas WC aparecendo em torno de 85.000 K).
• Equivalência Pós-Processamento: Os autores demonstraram que aplicar uma correção de pós-processamento aos modelos clássicos — utilizando o mesmo método de interpolação para recalcular $T_{\text{eff}}$ e $R_{\text{eff}}$ com base nos parâmetros internos do modelo clássico — produz trajetórias evolutivas quase indistinguíveis daquelas calculadas com a inclusão direta da atmosfera. Isso confirma que o impacto principal está na definição da superfície e não na evolução interna.

Significado e Alegações
O artigo afirma que a inclusão de modelos de atmosfera detalhados é necessária para reproduzir corretamente a posição de estrelas massivas evoluídas no diagrama HR, resolvendo a discrepância de longa data entre a $T_{\text{eff}}$ teórica e a análise espectral observacional. No entanto, os autores são modestos quanto ao escopo desta melhoria:

• O método corrige a aparência (propriedades superficiais) dos modelos, mas não altera a evolução interna ou química.
• O estudo limita-se a modelos não rotativos de metalicidade solar de 50 e 60 $M_\odot$.
• As taxas de perda de massa utilizadas foram reduzidas por um fator de 3 para garantir estabilidade numérica nos cálculos de atmosfera; assim, as reduções de temperatura relatadas representam um limite inferior do efeito esperado para fases WR livres de hidrogênio.
• Os autores reconhecem limitações quanto ao tratamento de zonas de convecção subsuperficiais (por exemplo, o pico de opacidade do ferro) e a física específica do mecanismo de lançamento de ventos WNh, notando que trabalhos adicionais são necessários para alcançar um tratamento totalmente autoconsistente em todas as fases evolutivas.

Em conclusão, o estudo estabelece que, embora atmosferas detalhadas sejam cruciais para definir as propriedades observáveis das estrelas WR, a evolução interna dessas estrelas é robusta em relação à escolha das condições de contorno atmosféricas, validando o uso de correções de pós-processamento como uma alternativa eficiente ao acoplamento direto em muitos contextos.