JWST Advanced Deep Extragalactic Survey (JADES) Data Release 5: Wisp Subtraction with the Non-negative Matrix Factorization Algorithm

Resumo Técnico: Lançamento de Dados 5 do JADES (JWST Advanced Deep Extragalactic Survey): Subtração de Wisps com o Algoritmo de Fatoração de Matriz Não Negativa

Definição do Problema
Os wisps são artefatos proeminentes de luz espalhada nas imagens de comprimento de onda curto (SW) do NIRCam do Telescópio Espacial James Webb (JWST), causados pelo espalhamento de luz fora do eixo no suporte do espelho secundário. Esses artefatos aparecem em regiões fixas dos detectores com uma morfologia que varia ligeiramente, mas visivelmente, entre as exposições, particularmente nas intensidades relativas das subestruturas. Essa variação complica a subtração, pois abordagens de template único (como o pipeline do STScI) deixam resíduos significativos. Os métodos existentes enfrentam um compromisso (trade-off): templates de empilhamento de mediana (ex: pipeline do STScI) oferecem alta sensibilidade a estruturas difusas, mas carecem de flexibilidade para variações morfológicas, enquanto métodos dinâmicos por nível de exposição permitem morfologias arbitrárias, mas sofrem de baixa relação sinal-ruído (SNR) e dependência de uma subtração de fontes precisa. Esses resíduos contaminam a fotometria de fontes tênues e a análise morfológica, enviesando medições para fontes mais fracas que 25 mag AB.

Metodologia
Os autores introduzem um novo método de subtração utilizando o algoritmo de Fatoração de Matriz Não Negativa (NMF) para modelar e remover os wisps. A abordagem constrói templates de wisps de múltiplos componentes, específicos por filtro e detector, utilizando observações profundas do NIRCam da pesquisa JADES e de outros programas.

• Estrutura Algorítmica: O problema é formulado como uma aproximação de matriz onde a matriz de dados $D$ (representando $m$ exposições e $n$ pixels) é decomposta em matrizes de template não negativas $W$ e matrizes de amplitude $A$ para minimizar a função de perda $L(A, W) = \frac{1}{2}\|D - AW\|^2$. A norma é ponderada pela variância inversa do ruído.
• Construção de Templates:
  • Dados de filtros de banda larga (F090W, F150W, F150W, F200W) e banda média são pré-processados (correção de flat-fielding, subtração de fundo, correção de ruído 1/f e rejeição de outliers).
  • Fontes e artefatos são mascarados para evitar contaminação.
  • O NMF é aplicado iterativamente. Os autores utilizam uma estratégia de inicialização sequencial, começando com um único template (a imagem de média empilhada) e adicionando gradualmente componentes (até três para os detectores A3/B4, um para A4/B3), usando os resultados anteriores como valores iniciais.
  • A restrição de não negatividade é crítica, garantindo consistência física com a natureza aditiva da luz espalhada e reduzindo a degenerescência do modelo.
  • Os templates são refinados em dados de resolução total, e uma correção de viés é aplicada para considerar o viés estatístico introduzido pela restrição de não negatividade em regiões dominadas pelo ruído.
• Estratégia de Subtração: Para exposições individuais, as amplitudes dos wisps são determinadas usando Mínimos Quadrados Não Negativos (NNLS) dentro de regiões de alto SNR de wisps. O processo inclui correção iterativa de ruído 1/f para tratar a degenerescência entre os dois artefatos. As incertezas provenientes tanto do ajuste de amplitude quanto das imperfeições do template são propagadas para os mapas de erro finais.

Principais Contribuições

• Templates Baseados em NMF: O artigo fornece templates de wisps de múltiplos componentes, liberados publicamente, que capturam tanto a morfologia dominante do wisp quanto as variações de exposição a exposição.
• Melhoria no Desempenho de Subtração: O método alcança um RMS residual inferior em regiões afetadas por wisps em comparação com os templates de mediana v4.0 do STScI e abordagens de PCA ponderadas pelo ruído.
• Precisão Fotométrica: O método reduz o viés fotométrico e a dispersão para fontes em regiões contaminadas por wisps para níveis quase consistentes com áreas de detectores limpas.
• Integração ao Pipeline: O método está integrado ao pipeline de redução de dados do JADES e foi projetado para ser facilmente aplicável a outros programas de observação sem a necessidade de nova construção de templates para cada campo.

Resultados

• Melhorias Visuais e Estatísticas: Na banda F150W (detector NIRCam B4), o método NMF reduz o RMS residual médio excedente em regiões de wisp para aproximadamente 50% do obtido com os templates do STScI. As características residuais mais nítidas são reduzidas para cerca de 30% do RMS residual excedente do STScI.
• Fotometria: Testes de fotometria de abertura em fontes que aparecem tanto em regiões afetadas por wisps quanto em regiões limpas mostram que o método NMF produz menor viés e redução de dispersão em comparação com os templates do STScI, com distribuições consistentes com regiões não afetadas.
• Qualidade do Mosaico: Mosaicos processados com templates NMF exibem fundos mais limpos e qualidade de imagem aprimorada, reduzendo a identificação errônea de fontes e melhorando a precisão das medições em comparação com aqueles processados com templates do STScI, que frequentemente retêm wisps residuais mesmo após o empilhamento.
• Estabilidade: A fluência (flux) do wisp foi encontrada como quase constante entre campos (GOODS-S e GOODS-N) e no tempo, com um desvio padrão de apenas ~12%, apoiando o uso de templates universais em diferentes programas.

Significância e Alegações
O artigo afirma resolver o compromisso fundamental entre sensibilidade e flexibilidade na subtração de wisps, combinando a sensibilidade de dados profundos e empilhados com a flexibilidade de uma base de template não negativa de múltiplos componentes. Os autores afirmam que este método entrega uma subtração de wisps mais precisa do que as abordagens existentes, apoiando diretamente a análise fotométrica e morfológica precisa de fontes tênues. O trabalho enfatiza que os templates NMF são robustos, estão disponíveis publicamente e são aplicáveis a futuros levantamentos extragalácticos do NIRCam, mitigando uma importante fonte de contaminação para a identificação de galáxias de alto redshift e outros estudos de fontes tênues. Os autores observam limitações em relação aos "bright wisps" (características distintas e raras com origens diferentes) e "claws" (artefatos dependentes de estrelas), que não são abordados por este método baseado em templates.