Recalibration of the H$\alpha$ surface brightness-radius relation for planetary nebulae using Gaia DR3: new distances and the Milky Way oxygen radial gradient

Resumo Técnico: Recalibração da Relação Brilho Superficial Hα–Raio para Nebulosas Planetárias usando Gaia DR3

Declaração do Problema
A distribuição espacial dos elementos químicos no disco galáctico serve como uma restrição crítica para modelos de evolução de galáxias. No entanto, estudos que utilizam nebulosas planetárias (NPs) como traçadores têm sido historicamente limitados por grandes incertezas na determinação de distâncias. Embora escalas de distância estatísticas, como a relação brilho superficial Hα–raio estabelecida por Frew et al. (2016, doravante F16), tenham melhorado as estimativas de distância, elas permanecem sujeitas a incertezas que frequentemente excedem fatores de dois. O advento da missão Gaia, particularmente o Lançamento de Dados 3 (DR3), oferece paralaxes astrométricas precisas para centenas de estrelas centrais de nebulosas planetárias (ECNPs), proporcionando uma oportunidade para medir distâncias diretamente e recalibrar escalas estatísticas a fim de reduzir essas incertezas.

Metodologia
Os autores apresentam uma recalibração abrangente da relação brilho superficial Hα–raio usando paralaxes do Gaia DR3. A metodologia envolve várias etapas-chave:

1. Identificação e Seleção de Fontes: Um catálogo de 1.200 NPs com coordenadas confiáveis e diâmetros angulares foi compilado. Os autores consultaram o arquivo Gaia para fontes dentro de um raio de busca definido pela geometria da nebulosa. Para garantir a confiabilidade astrométrica, as fontes foram filtradas com base no Erro de Peso Unitário Renormalizado (RUWE < 1,4) e no número de períodos de visibilidade (> 8).
2. Correspondência de ECNP: Um método cor-distância, aprimorado a partir de trabalhos anteriores, foi empregado para identificar as ECNPs mais prováveis. Isso envolveu a aplicação de restrições de temperatura (Teff > 13.000 K) e relações cor-temperatura. As identificações resultantes foram classificadas em grupos de confiabilidade (A, B, C) e cruzadas com o catálogo de Chornay & Walton (2021). Fontes com valores de extinção visual que se desviavam significativamente dos mapas de poeira galáctica foram excluídas, resultando em uma amostra final de 436 NPs com associações de alta probabilidade.
3. Derivação de Distância:
   • Recalibração: Usando 135 NPs com incertezas fracionárias de paralaxe ($f = \sigma_\varpi/\varpi$) $\le 0,20$, os autores realizaram uma regressão de Mínimos Quadrados Ordinários (OLS) com simulações de Monte Carlo (MC) para recalibrar a relação $\log R_{pc}$ vs. $\log S_{H\alpha}$.
   • Distâncias Bayesianas: Para uma subamostra de 415 NPs com paralaxes do Gaia, distâncias bayesianas foram derivadas. Esta abordagem usa a distância estatística recalibrada como prior e a paralaxe do Gaia como verossimilhança. Para as 711 NPs restantes sem contrapartes do Gaia, a distribuição posterior foi adotada como distribuição prior.
4. Análise de Abundância Química: As abundâncias de oxigênio (O/H) foram compiladas a partir de Stanghellini & Haywood (2018) e Bucciarelli & Stanghellini (2023, doravante BS23), recalculadas usando um fator de correção de ionização (ICF) consistente. A amostra foi separada em populações de disco fino, disco grosso, bojo e halo com base em velocidades peculiares ($V_{pec}$) derivadas dos movimentos próprios e velocidades radiais do Gaia.
5. Ajuste de Gradiente: O gradiente radial de O/H para 231 NPs do disco foi analisado usando regressão linear por partes. Três modelos foram testados: um ajuste linear único, um ajuste segmentado com uma quebra e um ajuste segmentado com duas quebras. A seleção do modelo foi guiada pelo Critério de Informação de Akaike (AIC). Adicionalmente, um mapa de abundância 2D foi gerado usando interpolação de krigagem universal.

Principais Contribuições e Resultados

• Escala de Distância Recalibrada: A nova calibração produz uma relação mais plana em comparação com F16, resultando em distâncias estatísticas mais curtas. A relação adotada é $\log R_{pc} = (-0,257 \pm 0,008) \times \log S_{H\alpha} - (1,51 \pm 0,03)$. Esta escala demonstra métricas de ajuste melhoradas ($\langle k \rangle = 1,020$, $\langle \sigma \rangle = 0,142$) em comparação com a calibração original F16 aplicada aos mesmos dados.
• Comparações de Distância: As distâncias bayesianas derivadas são, em média, 16% menores do que as de Chornay & Walton (2021) e 28% menores do que as distâncias F16. O método bayesiano reduz significativamente a dispersão e o viés para NPs com incertezas fracionárias de paralaxe $f < 0,5$.
• Gradiente Radial de Oxigênio: A análise de 231 NPs do disco revela que um gradiente linear único não é o melhor ajuste. O modelo com uma quebra é estatisticamente preferido (menor AIC), indicando uma mudança na incluição perto do raio solar ($R \approx 7,25 \pm 1,00$ kpc).
  • Disco Interno ($R \lesssim 8$ kpc): O gradiente é plano ou ligeiramente positivo ($0,026 \pm 0,028$ dex/kpc).
  • Disco Externo ($R \gtrsim 8$ kpc): O gradiente torna-se mais íngreme e negativo ($-0,034 \pm 0,009$ dex/kpc).
  • Um modelo de duas quebras também foi testado, mas produziu resultados inconclusivos sobre a segunda quebra, embora tenha sugerido uma estrutura em degrau perto do Sol.
• Disco Fino vs. Disco Grosso: Ao separar a amostra por cinemática, as NPs do disco fino exibem um gradiente mais íngreme ($-0,035 \pm 0,012$ dex/kpc) em comparação com o disco grosso ($-0,005 \pm 0,017$ dex/kpc). Os autores sugerem que a aparente quebra na amostra combinada pode resultar da superposição dessas populações estelares distintas.
• Assimetria Azimutal: O mapa de krigagem 2D revela uma modesta assimetria azimutal, com abundâncias aumentadas perto da barra galáctica em longitudes positivas. O círculo solar parece atuar como uma região de transição, criando uma estrutura de abundância bimodal entre as regiões internas e externas.

Significado e Alegações
O artigo afirma que a recalibração da relação brilho superficial Hα–raio usando paralaxes do Gaia DR3 fornece um indicador de distância mais robusto para NPs galácticas, abordando incertezas de longa data que dificultaram o estudo da evolução química.

O principal significado reside na evidência de um gradiente radial de oxigênio não linear no disco da Via Láctea. A detecção de uma quebra perto do raio solar desafia a suposição de um gradiente linear constante. Os autores propõem vários mecanismos físicos para esta característica:

1. Efeitos Dinâmicos: A quebra pode refletir mudanças na eficiência de formação estelar impulsionadas pela barra galáctica ou pela ressonância de corotação dos braços espirais.
2. Populações Estelares: A quebra observada pode ser um artefato da superposição dos discos fino e grosso, que possuem diferentes gradientes de metalicidade radial.
3. Evolução Química: Os resultados fornecem novas restrições para modelos de evolução química galáctica, particularmente no que diz respeito ao impacto de estruturas não axissimétricas (barras e espirais) nos fluxos de gás e distribuições de metalicidade.

Os autores concluem que, embora os dados apoiem uma mudança na inclinação perto do raio solar, são necessárias investigações adicionais com amostras maiores e homogêneas e incertezas reduzidas para confirmar definitivamente a origem dessas quebras e refinar os modelos da evolução química da Via Láctea.