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

Resumen Técnico: Un Nuevo Algoritmo de Deconvolución de la PSF para el Análisis de Galaxias Anfitrionas de Núcleos Galácticos Activos

Enunciado del Problema
Las galaxias anfitrionas de Núcleos Galácticos Activos (NGA) son a menudo difíciles de estudiar morfológicamente porque la fuente puntual central brillante (el NGA) contamina la luz de la galaxia anfitriona. Aunque telescopios espaciales de alta resolución como Hubble y JWST pueden resolver estas estructuras, sus campos de visión estrechos limitan los estudios estadísticos a muestras pequeñas ($\lesssim 10^4$ objetos). Las encuestas desde tierra (por ejemplo, Subaru/HSC, Rubin, Euclid, Roman) ofrecen muestras grandes pero sufren de menor resolución espacial y un significativo desenfoque por la Función de Dispersión del Punto (PSF). Los métodos tradicionales, como la descomposición de imágenes 2D (por ejemplo, GALFIT), son sensibles a la sobre-restación de la PSF y a subestructuras complejas, y luchan por separar la galaxia anfitriona del NGA sin introducir incertidumbres en parámetros estructurales como el índice de Sérsic. Los algoritmos existentes de deconvolución de PSF a menudo fallan en generar simultáneamente un componente extendido separado, un componente de fuente puntual separado y una imagen reconstruida combinada, los cuales son todos necesarios para un análisis morfológico y fotométrico robusto de anfitriones de NGA tenues y de alto desplazamiento al rojo.

Metodología
Los autores proponen un nuevo algoritmo de deconvolución de PSF diseñado para mejorar simultáneamente la resolución espacial y eliminar la fuente puntual central. El método modela la imagen intrínseca $x$ como la suma de dos componentes distintos: un componente extendido $I_{sm}$ (la galaxia anfitriona) y un componente de fuente puntual $I_{sp}$ (el NGA). El modelo de observación asume que el flujo observado sigue una distribución gaussiana basada en la convolución de estos componentes con la PSF.

Para reconstruir la imagen intrínseca, el algoritmo minimiza una función objetivo de log-posterior negativo sujeta a tres restricciones específicas:

1. Restricción de Suavidad ($I_{sm}$): Una regularización tipo Tikhonov aplicada al componente extendido para imponer suavidad espacial, modelando la estructura de la galaxia anfitriona.
2. Restricción de Esparsidad ($I_{sp}$): Una regularización con norma $L_1$ (LASSO) aplicada al componente de fuente puntual para localizar el flujo del NGA en un pequeño número de píxeles.
3. Restricción de Equilibrio de Fuente Puntual: Una restricción novedosa basada en el producto píxel a píxel $I_{sm} \times I_{sp}$. Este término penaliza la presencia simultánea de valores de flujo altos en ambos componentes, extendido y de fuente puntual, en la misma ubicación de píxel. Esto previene la "fuga" de flujo de fuente puntual hacia la galaxia anfitriona (sobre-restación) o la supresión del flujo central de la galaxia anfitriona (sub-restación).

La optimización se resuelve utilizando el Algoritmo Iterativo de Umbralización por Encogimiento (ISTA) con un esquema de retroceso para la constante de Lipschitz para garantizar una convergencia estable y monótona. El algoritmo genera tres imágenes: la galaxia anfitriona separada ($I_{sm}$), el NGA separado ($I_{sp}$) y la reconstrucción combinada ($I_{sm} + I_{sp}$).

Contribuciones Clave

• Algoritmo Novel: Desarrollo de un método de deconvolución de PSF de tres salidas específicamente adaptado para imágenes de NGA, separando explícitamente la galaxia anfitriona y el NGA mientras reconstruye la imagen total.
• Restricción de Equilibrio de Fuente Puntual: Introducción de un nuevo término de regularización ($I_{sm} \times I_{sp}$) que desacopla efectivamente los componentes extendido y de fuente puntual, abordando la interdependencia de hiperparámetros que afecta a los métodos anteriores.
• Mejora de Resolución de Tierra a Espacio: Demostración de que las imágenes desde tierra (Subaru/HSC) pueden ser deconvolucionadas para lograr una resolución espacial comparable a las imágenes del Telescopio Espacial Hubble (HST) mientras se elimina la contaminación del NGA.

Resultados
El algoritmo fue probado tanto en imágenes artificiales de NGA como en datos reales de NGA del campo COSMOS de Subaru/HSC (objetos J0959 en $z=0.34$ y J1000 en $z=1.38$), utilizando imágenes de HST/ACS como verdad fundamental.

• Análisis Cualitativo: Las imágenes reconstruidas de la galaxia anfitriona ($I_{sm}$) muestran una resolución espacial significativamente mejorada, revelando estructuras finas como brazos espirales y subestructuras agrupadas con alto contraste. La fuente puntual central se elimina efectivamente sin sobre-restación o sub-restación significativa. Las imágenes reconstruidas combinadas ($I_{sm} + I_{sp}$) se asemejan estrechamente a las imágenes de referencia de Hubble.
• Análisis Cuantitativo:
  • Convergencia: La distancia de Hellinger entre iteraciones disminuyó rápidamente, estabilizándose después de ~100 iteraciones.
  • Precisión: El Error Cuadrático Medio (RMSE) entre la imagen reconstruida (convolucionada con la PSF de HSC) y la observación original de HSC fue extremadamente bajo ($\sim 10^{-6}$ a $10^{-5}$). El RMSE entre la imagen reconstruida y la verdad fundamental de Hubble también fue pequeño ($\sim 10^{-2}$).
  • Mejora de Resolución: El Ancho a Media Altura (FWHM) de las estructuras reconstruidas de la galaxia anfitriona mejoró en un factor de $\sim 2-3$ en comparación con las imágenes originales de HSC, alcanzando niveles comparables a Hubble (por ejemplo, de $1.06''$ a $0.57''$ para J0959).
  • Eficacia de la Restricción: Se demostró que la inclusión de la restricción de equilibrio de fuente puntual ($\alpha_{bl}$) fue crítica para NGA reales, reduciendo aún más los valores de FWHM y mejorando la eliminación de fuentes puntuales en comparación con las ejecuciones sin esta restricción.

Significado
El artículo afirma que este algoritmo permite estudios morfológicos estadísticos de galaxias anfitrionas de NGA distantes utilizando datos de encuestas de gran campo de instalaciones actuales y futuras, incluido el Observatorio Vera C. Rubin, el Telescopio Espacial Euclid y el Telescopio Espacial Roman. Al eliminar efectivamente la contaminación del NGA y mejorar la resolución, el método permite la aplicación de índices morfológicos no paramétricos (por ejemplo, concentración, asimetría, coeficiente de Gini) a muestras grandes ($>10^4$ objetos) sin los sesgos causados previamente por el flujo del NGA. Se espera que esta capacidad resuelva la ambigüedad actual con respecto a la relación entre la actividad del NGA y las fusiones/interacciones de galaxias. Además, el método reduce la dependencia de la clasificación visual y puede extenderse a estudios de las relaciones entre agujeros negros supermasivos y bochas, así como a la extracción de otras estructuras como cúmulos estelares gigantes.