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

Technische Samenvatting: JWST Advanced Deep Extragalactic Survey (JADES) Data Release 5: Wisp-subtractie met het Non-negative Matrix Factorization Algoritme

Probleemstelling
Wisps zijn prominente artefacten van verstrooid licht in de kortgolvige (SW) NIRCam-afbeeldingen van de James Webb Space Telescope (JWST), veroorzaakt door lichtverstrooiing van de secundaire spiegelsteun (strut) buiten de optische as. Deze artefacten verschijnen in vaste regio's op de detectoren met een morfologie die weliswaar licht varieert, maar zichtbaar verschilt tussen blootstellingen, met name in de relatieve sterkte van de substructuren. Deze variatie bemoeilijkt subtractie, aangezien single-template benaderingen vaak significante residuen achterlaten. Bestaande methoden kampen met een afweging: median-stacked templates (bijv. STScI pipeline) bieden een hoge gevoeligheid voor diffuse structuren maar missen flexibiliteit voor morfologische variaties, terwijl dynamische, op blootstelling gebaseerde methoden flexibel zijn voor willekeurige morfologieën maar lijden onder een lage signaal-ruisverhouding (SNR) en een afhankelijkheid van nauwkeurige bronsubtractie. Deze residuen vervuilen de fotometrie van zwakke bronnen en morfologische analyses, wat metingen voor bronnen zwakker dan 25 AB mag biasert.

Methodologie
De auteurs introduceren een nieuwe subtractiemethode die gebruikmaakt van het Non-negative Matrix Factorization (NMF) algoritme om wisps te modelleren en te verwijderen. De aanpak construeert multi-component, filter- en detectorspecifieke wisp-templates met behulp van diepe NIRCam-waarnemingen uit de JADES-survey en andere programma's.

• Algoritmisch Kader: Het probleem wordt geformuleerd als een matrixbenadering waarbij de datamatrix $D$ (die $m$ blootstellingen en $n$ pixels vertegenwoordigt) wordt gedecomposeerd in niet-negatieve template-matrices $W$ en amplitude-matrices $A$ om de verliesfunctie $L(A, W) = \frac{1}{2}\|D - AW\|^2$ te minimaliseren. De norm wordt gewogen door de inverse ruisvariantie.
• Template Constructie:
  • Data van breedbandige (F090W, F115W, F150W, F200W) en mediumbandige filters worden voorbewerkt (flat-fielding, achtergrondsubtractie, 1/f-ruiscorrectie en outlier-rejectie).
  • Bronnen en artefacten worden gemaskeerd om contaminatie te voorkomen.
  • NMF wordt iteratief toegepast. De auteurs maken gebruik van een sequentiële initialisatiestrategie, beginnend met een enkele template (de gemiddelde gestapelde afbeelding) en het geleidelijk toevoegen van componenten (tot drie voor A3/B4 detectoren, één voor A4/B3), waarbij de vorige resultaten als initiële waarden dienen.
  • De non-negativiteitsbeperking is cruciand, omdat deze zorgt voor fysieke consistentie met de additieve aard van verstrooid licht en modeldegeneratie vermindert.
  • Templates worden verfijnd op volledige resolutie data, en een bias-correctie wordt toegepast om rekening te houden met de statistische bias die wordt geïntroduceerd door de non-negativiteitsbeperking in ruisgedomineerde regio's.
• Subtractiestrategie: Voor individuele blootstellingen worden de wisp-amplitudes bepaald met behulp van Non-negative Least Squares (NNLS) binnen regio's met een hoge SNR voor de wisps. Het proces omvat iteratieve 1/f-ruiscorrectie om degeneratie tussen de twee artefacten aan te pakken. Onzekerheden van zowel de amplitude-fitting als de imperfecties van de templates worden gepropageerd naar de uiteindelijke foutenkaarten.

Belangrijkste Bijdragen

• NMF-gebaseerde Templates: Het artikel levert publiekelijk beschikbare, multi-component wisp-templates die zowel de dominante wisp-morfologie als de variaties tussen blootstellingen en omgevingen vastleggen.
• Verbeterde Subtractieprestaties: De methode bereikt een lagere residuele root mean square (RMS) in wisp-geaffecteerde regio's vergeleken met de STScI v4.0 median templates en ruis-gewogen PCA-benaderingen.
• Fotometrische Nauwkeurigheid: De methode vermindert de fotometrische bias en spreiding voor bronnen in wisp-gecontamineerde regio's tot niveaus die bijna consistent zijn met schone detectorgebieden.
• Pipeline Integratie: De methode is geïntegreerd in de JADES-datareductiepipeline en is ontworpen om gemakkelijk toepasbaar te zijn op andere observatieprogramma's zonder dat er voor elk nieuw veld nieuwe templates geconstrueerd hoeven te worden.

Resultaten

• Visuele en Statistische Verbeteringen: In de F150W-band (NIRCam B4 detector) vermindert de NMF-methode de gemiddelde excessieve residuele RMS in wisp-regio's tot ongeveer 50% van die verkregen met STScI-templates. De scherpste residuele kenmerken worden gereduceerd tot ongeveer 30% van de STScI excess RMS.
• Fotometrie: Tests met apertuurfotometrie op bronnen die zowel in wisp-geaffecteerde als in schone regio's voorkomen, laten zien dat de NMF-methode een kleinere bias en verminderde spreiding oplevert vergeleken met STScI-templates, met distributies die consistent zijn met ongeaffecteerde regio's.
• Mosaic Kwaliteit: Mosaics die verwerkt zijn met NMF-templates vertonen schonere achtergronden en verbeterde beeldkwaliteit, wat bronmisidentificatie vermindert en de meetnauwkeurigheid verbetert vergeleken met afbeeldingen die verwerkt zijn met STScI-templates, die vaak zelfs na het stapelen nog residuele wisps bevatten.
• Stabiliteit: De wisp-flux blijkt vrijwel constant te zijn over verschillende velden (GOODS-S en GOODS-N) en tijd, met een standaarddeviatie van slechts ~12%, wat het gebruik van universele templates over verschillende programma's ondersteunt.

Betekenis en Claims
Het artikel beweert de fundamentele afweging tussen gevoeligheid en flexibiliteit bij wisp-subtractie op te lossen door de gevoeligheid van diepe, gestapelde data te combineren met de flexibiliteit van een multi-component, non-negatieve template-basis. De auteurs stellen dat deze methode een nauwkeurigere subtractie van wisps levert dan bestaande benaderingen, wat directe ondersteuning biedt voor precieze fotometrische en morfologische analyse van zwakke bronnen. De auteurs benadrukken dat de NMF-templates robuust en publiek beschikbaar zijn en toepasbaar zijn op toekomstige NIRCam extragalactische surveys, waardoor een belangrijke bron van contaminatie voor de identificatie van hoog-rood verschoven sterrenstelsels en andere studies naar zwakke bronnen wordt gemitigeerd. De auteurs merken beperkingen op met betrekking tot "heldere wisps" (zeldzame, onderscheidende kenmerken met een andere oorsprong) en "claws" (ster-afhankelijke artefacten), die niet worden geadresseerd door deze template-gebaseerde methode.