Joint Multiband Photometry with crowdsource

1. Probleemstelling

Grote infraroodopnames zoals de Wide-field Infrared Survey Explorer (WISE) staan voor aanzienlijke uitdagingen in drukke velden door bronvermenging, ruimtelijk variërende achtergronden en grote pixelschalen (2,75″/pixel). Bestaande fotometrische pipelines (zoals DAOPHOT, DOLPHOT) en de originele crowdsource-pipeline verwerken beeldbanden doorgaans onafhankelijk van elkaar. Deze onafhankelijke passing leidt tot:

• Inconsistente Bronposities: Bronnen worden gedetecteerd op licht verschillende coördinaten in verschillende banden.
• Flux-inconsistenties: Afgeleide kleuren en spectraalenergieverdelingen (SED's) lijden onder fouten door misalignement.
• Verminderde Gevoeligheid: Zwakke bronnen die detecteerbaar zijn in een diepere band, kunnen gemist worden in een ruizigere band als de bronlijst niet wordt gedeeld.
• Bias: Bij enkel-band passing kan de flux van een zwakke bron worden geabsorbeerd in het luchtachtergrondmodel, wat leidt tot onderschatte fluxen en vertekende kleuren.

Het artikel adresseert de behoefte aan een unifyend raamwerk dat kruisbandconsistentie afdwingt, terwijl het gebruik maakt van de computerefficiëntie van de crowdsource-pipeline.

2. Methodologie

De auteurs presenteren een multiband-extensie voor de crowdsource-pipeline, die gebruikmaakt van grote packages voor lineaire algebra met schaarse matrices om overlappende bronnen te passen. De kerninnovatie is een gelijktijdig, zelfconsistent model over meerdere beeldbanden.

Belangrijkste algoritmische componenten:

• Gecombineerde Detectie (Significantie-afbeeldingen):

  • In plaats van bronnen in elke band apart te detecteren, construeert het algoritme een gezamenlijke significantie-afbeelding ($S_{joint}$).
  • Dit combineert getallen voor gematchte filters en varianties uit alle banden, gewogen naar de aangenomen bron-SED.
  • $S_{joint} = N_{joint} / \sqrt{V_{joint}}$, waarbij $N_{joint}$ en $V_{joint}$ gewogen sommen zijn van individuele bijdragen per band.
  • Dit maakt detectie mogelijk van zwakke bronnen die marginaal zijn in individuele banden maar significant wanneer ze gecombineerd worden, wat vooral de diepte in ruizigere banden verbetert.

• Unifyend Lineair Kleinste-Kwadratenkader:

  • De fotometrische passing wordt geformuleerd als een lineair systeem $y \approx Ax$.
  • Data-vector ($y$): Gestapelde pixelwaarden uit alle banden.
  • Parameter-vector ($x$): Bevat gedeelde centroidparameters (posities) voor alle bronnen over alle banden, maar bandspecifieke fluxen en luchtparameters.
  • Ontwerpmatrix ($A$): Een blokgestructureerde matrix waarbij centroidkolommen worden gedeeld over banden, terwijl flux- en luchtkolommen bandspecifiek blijven.
  • Oplosser: Het systeem wordt iteratief opgelost met de LSQR-oplosser voor schaarse matrices.

• Iteratief Proces:

  1. Luchtschatting: Gladde achtergrondaftrekking.
  2. Borndetectie: Pieken geïdentificeerd in de gezamenlijke significantie-afbeelding.
  3. Lineaire passing: Gelijkmatige oplossing voor gedeelde posities en bandspecifieke fluxen.
  4. Verfijning: Iteratieve updates van posities en PSF-modellen.

3. Belangrijkste bijdragen

• Algoritmische Generalisatie: Succesvolle uitbreiding van de crowdsource-pipeline van enkel-band naar multiband passing zonder in te leveren op computerefficiëntie.
• Gedeelde Positie-Constraint: Afdwingen van één enkele, zelfconsistente set bronposities over alle banden, waardoor positiespreiding tussen filters wordt geëlimineerd.
• Optimale Detectiestrategie: Ontwikkeling van een statistisch consistente methode voor gezamenlijke brondetectie die zich aanpast aan verschillende PSF's en ruiseigenschappen van elke band.
• Implementatie: De code is geïmplementeerd in Python en publiek beschikbaar, gebouwd op de bestaande crowdsource-infrastructuur.

4. Resultaten en Validatie

De methode werd getest op unWISE samengevoegde tegels in het COSMOS-veld (spaarzaam tot gemiddelde dichtheid) en een Galactische Bulge-veld (hoge dichtheid), waarbij WISE W1 (3,4 $\mu$m) en W2 (4,6 $\mu$m) banden werden vergeleken met enkel-band passing en diepe Spitzer/IRAC-referentiedata.

• Flux-consistentie & Bias-reductie:

  • Multiband W2-fluxen zijn systematisch helderder dan enkel-band W2-fluxen aan het zwakke einde.
  • Dit corrigeert een bias waarbij enkel-band passing het licht van zwakke bronnen absorbeert in het luchtmodel. De multiband passing wijst dit licht correct toe aan bronnen, wat de fotometrie van zwakke bronnen verbetert met ongeveer 5% bij 17 mag.
  • Kleurstabiliteit ($W1-W2$) is aanzienlijk verbeterd, met bronnen die dichter bij de fysische hoofdreekslocus clusteren.

• Astrometrische Stabiliteit:

  • De gezamenlijke passing produceert één enkele centroid per bron, waardoor de positiespreiding van band tot band die bij onafhankelijke passing wordt waargenomen, wordt geëlimineerd.

• Bronherstel (Volledigheid):

  • De multiband-catalogus herstelde ongeveer 10.000 extra bronnen vergeleken met de W1 enkel-band catalogus.
  • Kruisverwijzing met Spitzer-data toonde aan dat ongeveer 80% van de nieuwe multiband-detecties (die geen enkel-band tegenhangers hadden) echte astrofysische objecten zijn.
  • Dieptewinst: De methode verbeterde de limietmagnitude in de ruizigere W2-band met ongeveer 1 magnitude vergeleken met enkel-band passing, terwijl een hoge volledigheid in W1 werd behouden.

• Betrouwbaarheid (Zuiverheid):

  • Multiband passing behield een hoge betrouwbaarheid (zuiverheid) zelfs bij zwakkere limieten, waardoor spookdetecties die diepe enkel-band zoektochten vaak teisteren, effectief werden onderdrukt.

5. Betekenis

Dit werk biedt een robuuste basis voor toekomstige multiband-catalogi voor drukke velden. Door positiemconsistentie af te dwingen en gebruik te maken van gecorreleerde informatie over banden heen, stelt het algoritme het volgende in staat:

• Meer accurate fotometrische roodverschuivingen en bronclassificatie.
• Verbeterde eigenbeweging-metingen door consistente astrometrie over epoches en banden heen.
• Maximalisatie van de wetenschappelijke opbrengst van bestaande opnamegegevens (zoals WISE/NEOWISE) door detectielimieten dieper te duwen zonder nieuwe waarnemingen te vereisen.
• Een flexibel raamwerk dat kan worden uitgebreid naar tijdsdomein-astronomie (meerdere epoches) en variërende brongewichten voor specifieke wetenschappelijke doelen (bijvoorbeeld hoge-roodverschuiving quasarren).

De auteurs plannen om deze pipeline toe te passen op de volledige diepte, multiband samengevoegde data van de gehele WISE-missie om de diepste mogelijke catalogus van infraroodbronnen tot nu toe te creëren.