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

Resumen Técnico: Liberación de Datos 5 del JADES (Advanced Deep Extragalactic Survey) del JWST: Sustracción de Wisps mediante el Algoritmo de Factorización de Matrices No Negativas

Planteamiento del Problema
Los wisps son artefactos prominentes de luz dispersa en las imágenes de longitud de onda corta (SW) de la cámara NIRCam del Telescopio Espacial James Webb (JWST), causados por la dispersión de luz fuera de eje en el soporte del espejo secundario. Estos artefactos aparecen en regiones fijas de los detectores con una morfología que varía ligeramente, pero visiblemente, entre exposiciones, particularmente en las intensidades relativas de las subestructuras. Esta variación complica la sustracción, ya que los enfoques de plantilla única (como el proceso de STScI) suelen dejar residuos significativos. Los métodos existentes enfrentan un compromiso: las plantillas de apilado de la mediana (por ejemplo, el pipeline de STScI) ofrecen una alta sensibilidad a las estructuras difusas pero carecen de flexibilidad para las variaciones morfológicas, mientras que los métodos dinámicos a nivel de exposición permiten morfologías arbitrarias pero sufren de bajas relaciones señal-ruido (SNR) y dependencia de una sustracción de fuentes precisa. Estos residuos contaminan la fotometría de fuentes tenues y el análisis morfológico, sesgando las mediciones para fuentes más tenues que 25 AB mag.

Metodología
Los autores introducen un nuevo método de sustracción utilizando el algoritmo de Factorización de Matrices No Negativas (NMF) para modelar y eliminar los wisps. El enfoque construye plantillas de wisps multicomponente, específicas de cada filtro y detector, utilizando observaciones profundas de NIRCam de la encuesta JADES y otros programas.

• Marco Algorítmico: El problema se formula como una aproximación de matrices donde la matriz de datos $D$ (que representa $m$ exposiciones y $n$ píxeles) se descompone en matrices de plantillas no negativas $W$ y matrices de amplitud $A$ para minimizar la función de pérdida $L(A, W) = \frac{1}{2}\|D - AW\|^2$. La norma está ponderada por la varianza inversa del ruido.
• Construcción de Plantillas:
  • Los datos de filtros de banda ancha (F090W, F115W, F150W, F200W) y de banda media han sido preprocesados (corrección de campo plano, sustracción de fondo, corrección de ruido 1/f y rechazo de valores atípicos).
  • Las fuentes y los artefactos se enmascaran para evitar la contaminación.
  • Se aplica NMF de forma iterativa. Los autores emplean una estrategia de inicialización secuencial, comenzando con una sola plantilla (la imagen de la media apilada) y añadiendo gradualmente componentes (hasta tres para los detectores A3/B4, uno para A4/B3) utilizando los resultados previos como valores iniciales.
  • La restricción de no negatividad es crítica, asegurando la consistencia física con la naturaleza aditiva de la luz dispersa y reduciendo la degeneración del modelo.
  • Las plantillas se refinan en datos de resolución completa, y se aplica una corrección de sesgo para dar cuenta del sesgo estadístico introducido por la restricción de no negatividad en regiones dominadas por el ruido.
• Estrategia de Sustracción: Para exposiciones individuales, las amplitudes de los wisps se determinan mediante Mínimos Cuadrados No Negativos (NNLS) dentro de regiones de wisps de alta SNR. El proceso incluye una corrección iterativa de ruido 1/f para abordar la degeneración entre los dos artefactos. Las incertidumbres, tanto de la sustracción de amplitud como de las imperfecciones de la plantilla, se propagan a los mapas de error finales.

Contribuciones Clave

• Plantillas basadas en NMF: El artículo proporciona plantillas de wisps multicomponente, de acceso público, que capturan tanto la morfología dominante del wisp como las variaciones de exposición a exposición.
• Mejora en el Desempeño de la Sustracción: El método logra un RMS residual más bajo en las regiones afectadas por wisps en comparación con las plantillas de la mediana de STScI v4.0 y los enfoques de PCA ponderados por ruido.
• Precisión Fotométrica: El método reduce el sesgo fotométrico y la dispersión para fuentes en regiones contaminadas por wisps a niveles casi consistentes con las áreas de detectores limpias.
• Integración en el Pipeline: El método está integrado en el pipeline de reducción de datos de JADES y está diseñado para ser aplicable fácilmente a otros programas de observación sin requerir la construcción de nuevas plantillas para cada campo.

Resultados

• Mejoras Visuales y Estadísticas: En la banda F150W (detector NIRCam B4), el método NMF reduce el RMS del exceso de residuo medio en las regiones de wisps a aproximadamente el 50% de lo obtenido con las plantillas de STScI. Las características residuales más nítidas se reducen a aproximadamente el 30% del exceso de RMS de STScI.
• Fotometría: Las pruebas de fotometría de apertura en fuentes que aparecen tanto en regiones afectadas por wisps como en regiones limpias muestran que el método NMF produce un sesgo menor y una dispersión reducida en comparación con las plantillas de STScI, con distribuciones consistentes con las regiones no afectadas.
• Calidad de los Mosaicos: Los mosaicos procesados con plantillas NMF exhiben fondos más limpios y una calidad de imagen mejorada, reduciendo la identificación errónea de fuentes y mejorando la precisión de las mediciones en comparación con aquellos procesados con plantillas de STScI, que a menudo retienen wisps residuales incluso después del apilado.
• Estabilidad: Se encontró que el flujo de los wisps es casi constante a través de los campos (GOODS-S y GOODS-N) y del tiempo, con una desviación estándar de solo ~12%, lo que respalda el uso de plantillas universales en diferentes programas.

Significancia y Reivindicaciones
El artículo afirma resolver el compromiso fundamental entre sensibilidad y flexibilidad en la sustracción de wisps al combinar la sensibilidad de los datos profundos apilados con la flexibilidad de una base de plantillas no negativas de múltiples componentes. Los autores declaran que este método ofrece una sustracción de wisps más precisa que los enfoques existentes, apoyando directamente el análisis fotométrico y morfológico preciso de fuentes tenues. El trabajo enfatiza que las plantillas NMF son robustas, están disponibles públicamente y son aplicables a futuros estudios extragalácticos de NIRCam, mitigando así una fuente importante de contaminación para la identificación de galaxias de alto desplazamiento al rojo y otros estudios de fuentes tenues. Los autores señalan limitaciones respecto a los "bright wisps" (características raras y distintas con orígenes diferentes) y "claws" (artefactos dependientes de la estrella), que no son abordados por este método basado en plantillas.