A New PSF Deconvolution Algorithm: Simultaneous Spatial Resolution Enhancement and Point Source Removal for Morphological Analysis of AGN Host Galaxies

Resumo Técnico: Um Novo Algoritmo de Deconvolução de PSF para Análise de Galáxias Hospedeiras de Núcleos Galácticos Ativos

Enunciação do Problema
As galáxias hospedeiras de Núcleos Galácticos Ativos (NGA) são frequentemente difíceis de estudar morfologicamente porque a fonte pontual central brilhante (o NGA) contamina a luz da galáxia hospedeira. Embora telescópios espaciais de alta resolução, como o Hubble e o JWST, possam resolver essas estruturas, seus campos de visão estreitos limitam estudos estatísticos a pequenas amostras ($\lesssim 10^4$ objetos). Levantamentos baseados no solo (por exemplo, Subaru/HSC, Rubin, Euclid, Roman) oferecem grandes amostras, mas sofrem com resolução espacial inferior e significativo borrão da Função de Espalhamento de Ponto (PSF). Métodos tradicionais, como a decomposição de imagens 2D (por exemplo, GALFIT), são sensíveis à subtração excessiva da PSF e a subestruturas complexas, e lutam para separar a galáxia hospedeira do NGA sem introduzir incertezas em parâmetros estruturais, como o índice de Sérsic. Algoritmos existentes de deconvolução de PSF frequentemente falham em produzir simultaneamente um componente estendido separado, um componente de fonte pontual separado e uma imagem reconstruída combinada, todos necessários para uma análise morfológica e fotométrica robusta de hospedeiros de NGA fracos e de alto desvio para o vermelho.

Metodologia
Os autores propõem um novo algoritmo de deconvolução de PSF projetado para aprimorar simultaneamente a resolução espacial e remover a fonte pontual central. O método modela a imagem intrínseca $x$ como a soma de dois componentes distintos: um componente estendido $I_{sm}$ (a galáxia hospedeira) e um componente de fonte pontual $I_{sp}$ (o NGA). O modelo de observação assume que o fluxo observado segue uma distribuição Gaussiana baseada na convolução desses componentes com a PSF.

Para reconstruir a imagem intrínseca, o algoritmo minimiza uma função objetivo de log-posterior negativo sujeita a três restrições específicas:

1. Restrição de Suavidade ($I_{sm}$): Uma regularização tipo Tikhonov aplicada ao componente estendido para impor suavidade espacial, modelando a estrutura da galáxia hospedeira.
2. Restrição de Esparsidade ($I_{sp}$): Uma regularização da norma $L_1$ (LASSO) aplicada ao componente de fonte pontual para localizar o fluxo do NGA em um pequeno número de pixels.
3. Restrição de Equilíbrio de Fonte Pontual: Uma restrição inovadora baseada no produto pixel a pixel $I_{sm} \times I_{sp}$. Este termo penaliza a presença simultânea de valores de fluxo elevados tanto no componente estendido quanto no componente de fonte pontual na mesma localização de pixel. Isso previne o "vazamento" do fluxo da fonte pontual para a galáxia hospedeira (subtração excessiva) ou a supressão do fluxo central da galáxia hospedeira (subtração insuficiente).

A otimização é resolvida usando o Algoritmo Iterativo de Encolhimento-Thresholding (ISTA) com um esquema de retrocesso para a constante de Lipschitz, garantindo convergência estável e monótona. O algoritmo produz três imagens: a galáxia hospedeira separada ($I_{sm}$), o NGA separado ($I_{sp}$) e a reconstrução combinada ($I_{sm} + I_{sp}$).

Principais Contribuições

• Algoritmo Inovador: Desenvolvimento de um método de deconvolução de PSF com três saídas, especificamente adaptado para imagens de NGA, separando explicitamente a galáxia hospedeira e o NGA enquanto reconstrói a imagem total.
• Restrição de Equilíbrio de Fonte Pontual: Introdução de um novo termo de regularização ($I_{sm} \times I_{sp}$) que desacopla efetivamente os componentes estendido e de fonte pontual, abordando a interdependência de hiperparâmetros que aflige métodos anteriores.
• Aprimoramento de Resolução do Solo para o Espaço: Demonstração de que imagens baseadas no solo (Subaru/HSC) podem ser deconvolvidas para alcançar resolução espacial comparável às imagens do Telescópio Espacial Hubble (HST), removendo simultaneamente a contaminação do NGA.

Resultados
O algoritmo foi testado tanto em imagens artificiais de NGA quanto em dados reais de NGA do campo COSMOS do Subaru/HSC (objetos J0959 em $z=0.34$ e J1000 em $z=1.38$), utilizando imagens do HST/ACS como verdade terrestre.

• Análise Qualitativa: As imagens reconstruídas da galáxia hospedeira ($I_{sm}$) mostram resolução espacial significativamente aprimorada, revelando estruturas finas, como braços espirais e subestruturas agrupadas, com alto contraste. A fonte pontual central é efetivamente removida sem subtração excessiva ou insuficiente significativa. As imagens reconstruídas combinadas ($I_{sm} + I_{sp}$) assemelham-se estreitamente às imagens de referência do Hubble.
• Análise Quantitativa:
  • Convergência: A distância Hellinger entre as iterações diminuiu rapidamente, estabilizando após ~100 iterações.
  • Precisão: O Erro Quadrático Médio (RMSE) entre a imagem reconstruída (convolvida com a PSF do HSC) e a observação original do HSC foi extremamente baixo ($\sim 10^{-6}$ a $10^{-5}$). O RMSE entre a imagem reconstruída e a verdade terrestre do Hubble também foi pequeno ($\sim 10^{-2}$).
  • Aprimoramento de Resolução: A Largura Total a Meia Altura (FWHM) das estruturas reconstruídas da galáxia hospedeira melhorou por um fator de $\sim 2-3$ em comparação com as imagens originais do HSC, atingindo níveis comparáveis ao Hubble (por exemplo, de $1.06''$ para $0.57''$ para J0959).
  • Eficácia da Restrição: A inclusão da restrição de equilíbrio de fonte pontual ($\alpha_{bl}$) mostrou-se crítica para NGAs reais, reduzindo ainda mais os valores de FWHM e melhorando a remoção da fonte pontual em comparação com execuções sem essa restrição.

Significado
O artigo afirma que este algoritmo permite estudos morfológicos estatísticos de galáxias hospedeiras de NGA distantes usando dados de levantamentos de campo amplo de instalações atuais e futuras, incluindo o Observatório Vera C. Rubin, o Telescópio Espacial Euclid e o Telescópio Espacial Roman. Ao remover efetivamente a contaminação do NGA e aprimorar a resolução, o método permite a aplicação de índices morfológicos não paramétricos (por exemplo, concentração, assimetria, coeficiente de Gini) a grandes amostras ($>10^4$ objetos) sem os vieses anteriormente causados pelo fluxo do NGA. Espera-se que essa capacidade resolva a ambiguidade atual sobre a relação entre a atividade do NGA e fusões/interações de galáxias. Além disso, o método reduz a dependência da classificação visual e pode ser estendido para estudos de relações entre buracos negros supermassivos e bojos, bem como para a extração de outras estruturas, como aglomerados estelares gigantes.