Impact of stellar rotation on type II supernova progenitor masses from pre-explosion imaging

1. Problema e Contexto

As estimativas das massas iniciais dos progenitores de supernovas do tipo II (SN II) são tradicionalmente derivadas de imagens pré-explosão, convertendo a luminosidade observada do progenitor em massa inicial utilizando modelos de evolução estelar. No entanto, a maioria desses modelos assume estrelas não rotativas ou uma rotação fixa (geralmente uma fração específica da velocidade crítica, como $\Omega_{ini}/\Omega_{crit} = 0.4$).

O problema central abordado neste trabalho é que as estrelas massivas apresentam uma ampla distribuição observada de velocidades rotacionais iniciais na sequência principal zero (ZAMS). A rotação estelar afeta significativamente a evolução estelar através de:

• Mistura química interna: Instabilidades rotacionais transportam novo combustível para o núcleo, estendendo as fases de queima nuclear e aumentando a massa do núcleo de hélio.
• Perda de massa: A rotação pode aumentar as taxas de perda de massa por ventos estelares devido a efeitos centrífugos.

Esses fatores alteram a relação entre a massa inicial e a luminosidade final. Se a rotação não for considerada corretamente (ou seja, se assumir-se uma rotação fixa ou nula quando a realidade é uma distribuição), as massas iniciais inferidas podem estar enviesadas. O objetivo do estudo é quantificar o impacto de utilizar a distribuição observada de velocidades rotacionais iniciais nas estimativas de massa dos progenitores de SN II, comparando-a com os modelos tradicionais não rotativos.

2. Metodologia

Os autores empregaram uma abordagem baseada em simulações de Monte Carlo (MC) combinando dados observacionais com novos grids de modelos estelares:

• Amostra de Dados: Foi compilada uma amostra de 21 progenitores de SN II detectados em imagens pré-explosão com estimativas diretas de luminosidade (excluindo limites superiores e casos com incertezas excessivas ou inconsistências, como SN 2003gd e SN 2020jfo).
• Modelagem Estelar (MESA):
  • Utilizou-se o código MESA (Modules for Experiments in Stellar Astrophysics) para calcular modelos de evolução de estrelas massivas de única massa com metalicidade solar ($Z = 0.0154$).
  • Os modelos cobrem massas iniciais de $9$a$30 M_\odot$ (o limite superior varia com a velocidade de rotação).
  • Rotação: Foram incluídos efeitos de mistura química induzida por rotação (instabilidades de cisalhamento, circulação de Eddington-Sweet, etc.) e transporte de momento angular (dínamo de Spruit-Tayler).
  • As velocidades rotacionais iniciais ($v_{ini}$) variaram de $0$a$600$ km/s.
• Distribuição de Velocidades: Em vez de assumir uma rotação fixa, os autores extraíram as velocidades iniciais de uma distribuição observada de estrelas O e B galácticas (Holgado et al. 2022), utilizando o algoritmo de deconvolução de Lucy (1974) para recuperar a distribuição intrínseca a partir das velocidades projetadas ($v \sin i$).
• Inferência de Massa (Método DBR):
  • Para cada progenitor, realizou-se um processo de Monte Carlo (10.000 iterações).
  • Em cada iteração, amostrou-se uma luminosidade (baseada na observação e sua incerteza) e uma velocidade rotacional inicial (baseada na distribuição observada).
  • Utilizando interpolação multidimensional nos grids de modelos rotativos, inferiu-se a massa inicial correspondente.
  • Este método foi denominado DBR (Distribution-Based Rotating).
  • Os resultados foram comparados diretamente com aqueles obtidos usando apenas modelos não rotativos.

3. Contribuições Principais

• Integração de Distribuições Observacionais: Pela primeira vez, uma análise sistemática de massas de progenitores de SN II incorporou explicitamente a distribuição completa de velocidades rotacionais observadas, em vez de um valor fixo ou zero.
• Análise Estatística Rigorosa: O estudo não apenas compara massas individuais, mas constrói funções de densidade de probabilidade (PDF), funções de distribuição cumulativa (CDF) e determina limites superiores de massa para a população inteira.
• Avaliação do "Problema das RSG": O trabalho reavalia o limite superior de massa dos progenitores de SN II (o chamado "Problema da Supergigante Vermelha" ou RSG problem), que sugere uma ausência de progenitores acima de $\sim 16-18 M_\odot$ em contraste com modelos teóricos.
• Exploração de Canais Binários: Inclui uma análise preliminar sobre o impacto de progenitores que são "ganhadores de massa" em sistemas binários, que podem atingir rotação crítica.

4. Resultados Chave

• Deslocamento para Massas Menores: Em todos os casos analisados, as distribuições de massa inferidas pelos modelos DBR (rotativos com distribuição observada) estão ligeiramente deslocadas para massas iniciais mais baixas em comparação com os modelos não rotativos. Isso ocorre porque a rotação aumenta a luminosidade final para uma dada massa inicial; portanto, para uma luminosidade observada fixa, a massa inicial necessária é menor.
• Magnitude do Efeito: O deslocamento é modesto. As diferenças típicas entre as medianas das massas inferidas variam de $0,2$a$1,1 M_\odot$. Essas diferenças são menores que as incertezas observacionais e de modelagem típicas.
• Limites Superiores de Massa:
  • Usando modelos não rotativos: $M_{up} = 20,8^{+2,8}_{-1,3} M_\odot$.
  • Usando modelos DBR: $M_{up} = 20,4^{+2,3}_{-1,9} M_\odot$.
  • A diferença entre os dois métodos é menor que a incerteza estatística.
• Problema das RSG: O limite superior de massa de $\sim 20,4 M_\odot$ (com limite superior de confiança de 95% em $28,5 M_\odot$) mantém a significância estatística do "Problema das RSG" em torno de $2\sigma$. Ou seja, ainda há uma tensão entre a ausência observada de progenitores muito massivos e a previsão de que estrelas acima de $25 M_\odot$ deveriam explodir como SN II, mas a rotação não resolve essa tensão por si só.
• Progenitores de Baixa Luminosidade: O estudo confirma que alguns progenitores (como SN 2008bk, 2009md) têm luminosidades inferidas abaixo do mínimo previsto pelos modelos de evolução padrão, sugerindo lacunas na modelagem das fases finais de estrelas de baixa massa ($7-10 M_\odot$) ou a necessidade de considerar canais como explosões de captura de elétrons (ECSN).

5. Significado e Conclusões

O estudo conclui que, embora a rotação estelar tenha um impacto físico significativo na evolução das estrelas massivas (alterando luminosidades e estruturas internas), sua inclusão explícita utilizando a distribuição observada de velocidades rotacionais tem um impacto modesto nas estimativas de massa dos progenitores de SN II dentro das incertezas atuais.

• Implicação Prática: Os modelos não rotativos continuam sendo suficientes para estimar as massas dos progenitores de SN II com precisão comparável aos modelos rotativos complexos, dado o nível atual de incertezas observacionais.
• Validação: A consistência entre os métodos reforça a robustez das estimativas anteriores de massa, mas destaca que a rotação deve ser considerada em estudos de alta precisão futura ou para entender a dispersão nas propriedades dos progenitores.
• Futuro: O trabalho aponta que a metalicidade e os canais de interação binária (especialmente ganhadores de massa) são fatores críticos que podem alterar mais significativamente as estimativas de massa do que a rotação isolada de estrelas de única massa, sugerindo direções para pesquisas futuras.

Em resumo, a rotação desloca as estimativas de massa para valores ligeiramente menores, mas não altera fundamentalmente as conclusões sobre a distribuição de massas dos progenitores de SN II ou a natureza do "Problema das RSG" com os dados atuais.