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

Technisches Resümee: JWST Advanced Deep Extragalactic Survey (JADES) Data Release 5: Wisp-Subtraktion mit dem Non-negative Matrix Factorization Algorithmus

Problemstellung
Wisps sind markante Streulicht-Artefakte in der Kurzwellen-Bildgebung (SW) des James Webb Space Telescope (JWST)/NIRCam, die durch die Streuung von Licht abseits der optischen Achse an der Stützstruktur des Sekundärspiegels verursacht werden. Diese Artefakte erscheinen in festen Regionen der Detektoren mit einer Morphologie, die sich zwischen den Belichtungen in ihrer Struktur und Intensität leicht, aber sichtbar verändert. Diese Variation erschwert die Subtraktion, da Single-Template-Ansätze (Ansätze mit einem einzelnen Referenzbild) oft signifikante Residuen hinterlassen. Bestehende Methoden stehen vor einem Zielkonflikt: Median-gestapelte Templates (z. B. STScI Pipeline) bieten eine hohe Sensitivität gegenüber diffusen Strukturen, lassen aber wenig Spielraum für morphologische Variationen, während dynamische, bel exposures-spezifische Methoden zwar beliebige Morphologien zulassen, aber unter geringem Signal-Rausch-Verhältnis (SNR) leiden und auf einer genauen Quellen-Subtraktion basieren müssen. Diese Residuen kontaminieren die Photometrie schwacher Quellen sowie die morphologische Analyse und verfälschen Messungen für Quellen, die jensefer 25 AB mag liegen.

Methodik
Die Autoren stellen eine neue Subtraktionsmethode vor, die den Non-negative Matrix Factorization (NMF) Algorithmus nutzt, um Wisps zu modellieren und zu entfernen. Der Ansatz konstruiert Multi-Komponenten-Wisp-Templates, die filterspezifisch und detektorspezifisch sind, unter Verwendung tiefer NIRCam-Beobachtungen aus dem JADES-Survey und anderen Programmen.

• Algorithmischer Rahmen: Das Problem wird als Matrixapproximation formuliert, bei der die Datentrix $D$ (repräsentiert durch $m$ Belichtungen und $n$ Pixel) in nicht-negative Template-Matrizen $W$ und Amplituden-Matrizen $A$ zerlegt wird, um die Verlustfunktion $L(A, W) = \frac{1}{2}\|D - AW\|^2$ zu minimieren. Die Norm ist gewichtet nach der inversen Rauschvarianz.
• Template-Konstruktion:
  • Daten aus Breitband- (F090W, F115W, F150W, F200W) und Mittelband-Filtern werden vorverarbeitet (Flat-Fielding, Hintergrundsubtraktion, 1/f-Rauschkorrektur, Ausreißer-Eliminierung).
  • Quellen und Artefakte werden maskiert, um Kontaminationen zu verhindern.
  • NMF wird iterativ angewendet. Die Autoren verwenden eine sequentielle Initialisierungsstrategie, die mit einem einzelnen Template (dem gemittelten Stack-Bild) beginnt und sukzessive weitere Komponenten hinzufügt (bis zu drei für A3/B4-Detektoren, eine für A4/B3) unter Verwendung der vorherigen Ergebnisse als Initialwerte.
  • Die Nicht-Negativitäts-Beschränkung ist entscheidend, da sie die physikalische Konsistenz mit der additiven Natur des Streulichts sicherstellt und Modell-Degenerationen reduziert.
  • Die Templates werden auf Vollauflösungsdaten verfeinert, und eine Bias-Korrektur wird angewendet, um den statistischen Bias zu berücksichtigen, der durch die Nicht-Negativitäts-Beschränkung in rauschdominierten Regionen eingeführt wird.
• Subtraktionsstrategie: Für einzelne Belichtungen werden die Wisp-Amplituden mittels Non-negative Least Squares (NNLS) innerhalb von Bereichen mit hohem SNR der Wisps bestimmt. Der Prozess beinhaltet eine iterative 1/f-Rauschkorrektur, um die Degenerierung zwischen den beiden Artefakten zu adressen. Unsicherheiten aus der Amplitudenanpassung sowie aus Template-Imperfektionen werden in die finalen Fehlerkarten propagiert.

Zentrale Beiträge

• NMF-basierte Templates: Die Arbeit stellt öffentlich zugängliche, Multi-Komponenten-Wisp-Templates bereit, die sowohl die dominante Wisp-Morphologie als auch die Variationen zwischen den Belichtungen erfassen.
• Verbesserte Subtraktionsleistung: Die Methode erreicht im Vergleich zu den STScI v4.0 Median-Templates und den rauschgewichteten PCA-Ansätzen einen niedrigeren residuellen Root Mean Square (RMS) in den Wisp-betroffenen Regionen.
• Photometrische Genauigkeit: Die Methode reduziert den photometrischen Bias und die Streuung für Quellen in Wisp-kontaminierten Regionen auf ein Niveau, das nahezu konsistent mit sauberen Detektorbereichen ist.
• Pipeline-Integration: Die Methode ist in die JADES-Datenreduktions-Pipeline integriert und so konzipiert, dass sie ohne die Konstruktion neuer Templates für jedes Feld auf andere Beobachtungsprogramme anwendbar ist.

Ergebnisse

• Visuelle und statistische Verbesserungen: Im F150W-Band (NIRCam B4 Detektor) reduziert die NMF-Methode den mittleren überschüssigen Residual-RMS in Wisp-Regionen auf etwa 50 % des Wertes der STScI-Templates. Die schärfsten Residual-Merkmale werden auf etwa 30 % des STScI-Überschuss-RMS reduziert.
• Photometrie: Tests der Apertur-Photometrie an Quellen, die sowohl in Wisp-betroffenen als auch in sauberen Regionen erscheinen, zeigen, dass die NMF-Methode einen geringeren Bias und eine reduzierte Streuung im Vergleich zu STScI-Templates liefert, wobei die Verteilungen konsistent mit unbeeinflussten Regionen sind.
• Mosaik-Qualität: Durch NMF-Templates verarbeitete Mosaike weisen sauberere Hintergründe und eine verbesserte Bildqualität auf, was die Fehlidentifikation von Quellen reduziert und die Messgenauigkeit verbessert, verglichen mit Prozessen, die mit STScI-Templates bearbeitet wurden, welche oft auch nach dem Stacking restliche Wisps aufweisen.
• Stabilität: Es wurde festgestellt, dass der Wisp-Fluss über verschiedene Felder (GOODS-S und GOODS-N) und Zeiten hinweg nahezu konstant ist, mit einer Standardabweichung von nur etwa 12 %, was die Verwendung universeller Templates für verschiedene Programme unterstützt.

Bedeutung und Behauptungen
Das Paper behauptet, den grundlegenden Zielkonflikt zwischen Sensitivität und Flexibilität bei der Wisp-Subtraktion gelöst zu haben, indem es die Sensitivität tiefer, gestapelter Daten mit der Flexibilität einer Multi-Komponenten-Nicht-Negativitäts-Template-Basis kombiniert. Die Autoren geben an, dass diese Methode eine präzisere Subtraktion von Wisps im Vergleich zu bestehenden Ansätzen ermöglicht und somit die präzise photometrische und morphologische Analyse schwacher Quellen direkt unterstützt. Die Arbeit betont, dass die NMF-Templates robust und öffentlich verfügbar sind sowie für zukünftige extragalaktische NIRCam-Surveys anwendbar sind, wodurch eine wesentliche Quelle der Kontamination für die Identifizierung von Hochrot-Galaxien und andere Studien schwacher Quellen gemindert wird. Die Autoren weisen auf Einschränkungen hinsichtlich „heller Wisps“ (seltene, distinkte Merkmale mit anderem Ursprung) und „Claws“ (sternabhängige Artefakte) hin, die durch diesen Template-basierten Ansatz nicht adressiert werden.