Charting circumstellar chemistry of carbon-rich asymptotic giant branch stars. II. Abundances and spatial distributions of CS

Resumo Técnico: Mapeamento da química circumestelar de estrelas do ramo gigante assintótico ricas em carbono II. Abundâncias e distribuições espaciais de CS

Problema e Motivação
As envoltórias circumestelares (CSEs) de estrelas do ramo gigante assintótico (AGB) ricas em carbono são locais complexos de formação de poeira e química molecular. Embora a química dessas envoltórias seja bem estudada para a estrela de carbono arquetípica IRC +10 216, nossa compreensão de outras estrelas de carbono permanece limitada, frequentemente dependendo de observações de antena única (SD) que carecem de resolução espacial. Estudos anteriores de estrelas de carbono distintas de IRC +10 216 assumiram amplamente o equilíbrio termodinâmico local (LTE) para derivar abundâncias, o que fornece apenas estimativas de primeira ordem. Além disso, modelos existentes de transferência radiativa (RT) não-LTE foram tipicamente restringidos ou por espectros de SD (fornecendo informações de excitação) ou por mapas interferométricos de algumas linhas (fornecendo informações espaciais), mas raramente por ambos simultaneamente. Este artigo aborda a necessidade de restringir processos químicos circumstelares em uma amostra de estrelas de carbono para determinar se IRC +10 216 é um arquétipo representativo e para refinar as estimativas de abundância para o monossulfeto de carbono (CS), uma molécula precursora chave.

Metodologia
O estudo foca em cinco estrelas AGB ricas em carbono: IRAS 15194−5115, IRAS 15082−4808, IRAS 07454−7112, AFGL 3068 e IRC +10 216. Os autores empregaram uma abordagem de modelagem em múltiplas etapas:

1. Propriedades Estelares e de Poeira: Distribuições de Energia Espectral (SEDs) foram modeladas usando o código Monte Carlo MCMax. Espectros estelares de entrada foram derivados de atmosferas modelo CO-MARCS (ou corpos negros, quando necessário), e propriedades de poeira foram restringidas usando fluxos fotométricos do Gaia, 2MASS, IRAS, Akari e outros arquivos. Isso resultou em luminosidade estelar ($L_\star$), temperatura ($T_\star$) e taxas de perda de massa de poeira ($\dot{M}_d$).
2. Taxas de Perda de Massa Gasosa: Emissão de linha CO foi modelada usando um código Monte Carlo 1D não-LTE (MCP). Ao ajustar linhas de CO de arquivos de SD e do Atacama Compact Array (ACA), os autores derivaram taxas de perda de massa gasosa ($\dot{M}_g$) e perfis de temperatura cinética para as envoltórias.
3. Modelagem de Transferência Radiativa de CS: O núcleo do estudo envolveu modelagem detalhada de RT não-LTE de emissão de linhas de CS, $^{13}$CS e C$^{34}$S. Os autores utilizaram um código de iteração lambda acelerada (ALI) 1D. Crucialmente, os modelos foram restringidos simultaneamente por:
   • Informação Espacial: Intensidades integradas extraídas de mapas do ALMA usando aberturas progressivamente crescentes para a linha CS $J=2-1$ (e transições $J$ superiores onde disponíveis).
   • Informação de Excitação: Intensidades de linhas de observações de SD (APEX, IRAM 30m, Yebes 40m, OSO e Herschel/HIFI) cobrindo uma ampla gama de transições rotacionais ($J=1-0$ a $J=22-21$).
   • Incertezas Correlacionadas: A minimização do $\chi^2$ levou em conta a correlação entre espectros do ALMA extraídos do mesmo cubo de dados.

O perfil de abundância de CS foi modelado como uma distribuição gaussiana com pico central, $f(r) = f_0 \exp[-(r/R_e)^2]$, onde $f_0$ é a abundância fracional de pico e $R_e$ é o raio de decaimento exponencial. Para isotopólogos, a extensão espacial ($R_e$) foi fixada ao valor de melhor ajuste do isotopólogo principal, permitindo que apenas a abundância de pico variasse.

Principais Resultados

• Perfis de Abundância de CS: O estudo reproduziu com sucesso a emissão de linhas de CS em várias transições e aberturas para todas as cinco fontes. As abundâncias fracionais de pico derivadas ($f_0$) variam de $1 \times 10^{-6}$ a $4 \times 10^{-6}$. Os raios de decaimento exponencial ($R_e$) variam de $1,8 \times 10^{16}$ cm a $6,8 \times 10^{16}$ cm.
• Razões Isotópicas: Razões confiáveis de $^{12}$C/$^{13}$C e $^{32}$S/$^{34}$S foram derivadas da modelagem de isotopólogos. As razões de $^{12}$C/$^{13}$C variam de $\sim7$ (para IRAS 15194−5115) a $\sim65$ (para AFGL 3068). As razões de $^{32}$S/$^{34}$S são geralmente consistentes com o valor solar ($\sim22$).
• Comparação com Trabalhos Anteriores: Os resultados refinam estimativas de abundância baseadas em SD anteriores. Especificamente, para IRAS 07454−7112, a incerteza relativa no raio de decaimento exponencial do CS foi melhorada por um fator de $\sim2,5$ em comparação com trabalhos anteriores. As abundâncias de CS derivadas são significativamente menores (por um fator de 2,5–6,5) do que aquelas calculadas em uma análise LTE anterior (UK24) das mesmas fontes, destacando a necessidade de modelagem não-LTE.
• Comparação com Modelos Químicos: As extensões radiais observadas dos perfis de abundância de CS concordam bem com previsões de modelos cinéticos químicos avançados (usando o banco de dados UMIST Rate22) que incorporam fotodissociação. Os dados suportam o quadro de fotodissociação usado nesses modelos químicos, embora distinguir entre taxas específicas de fotodissociação (por exemplo, Pattillo et al. 2018 vs. Hrodmarsson & van Dishoeck 2023) permaneça difícil dentro das incertezas atuais.
• Física de Excitação: A modelagem confirma a importância do bombeamento IR (excitação por fótons IR da poeira) nas CSEs. Negligenciar níveis vibracionalmente excitados leva a uma subestimação das intensidades de linha de $\sim50-70\%$ e resultaria em abundâncias e raios superestimados.

Significado e Alegações
Os autores afirmam que este trabalho representa um avanço necessário na modelagem de emissão de linhas circumstelares ao combinar estritamente restrições espaciais e de excitação. O estudo demonstra que:

1. IRC +10 216 é Representativo: Além de variações impulsionadas pela densidade, os perfis de abundância de CS derivados para IRC +10 216 não mostram diferenças significativas das outras estrelas de carbono na amostra, apoiando seu status como arquétipo para a física e química de CSE de estrelas de carbono.
2. Robustez Metodológica: O uso simultâneo de dados interferométricos e de SD fornece restrições robustas sobre perfis de abundância que nenhum método poderia alcançar sozinho, particularmente na quebra de degenerescências entre abundância de pico e extensão espacial.
3. Validação Química: A concordância entre extensões radiais derivadas de RT e previsões de modelos químicos corrobora os quadros de fotodissociação usados na modelagem astroquímica para envoltórias ricas em carbono.

O artigo conclui que, embora a origem heterogênea dos dados introduza incertezas inerentes de calibração, a abordagem de modelagem combinada produz as estimativas mais precisas de abundâncias e distribuições espaciais de CS até a data para esta amostra, estabelecendo um padrão para futuros estudos de outras moléculas precursoras em CSEs de AGB.