Foreground removal in HI 21 cm intensity mapping under frequency-dependent beam distortions

Technische Samenvatting: Voorgrondverwijdering in HI 21 cm intensiteitsmapping onder frequentieafhankelijke beam-distorties

Probleemstelling
Neutrale waterstof (HI) intensiteitsmapping is een primaire sonde voor kosmologie in het post-reïonisatie tijdperk, wat een tomografisch beeld van het Universum biedt. De extractie van het zwakke HI-signaal wordt echter ernstig bemoeilijkt door astrofysische voorgronden (synchrotron, free-free en puntbronnen) die tot vijf grootteordes helderder zijn. Hoewel Blind Source Separation (BSS) en model-fitting methoden veelbelovend zijn gebleken, gaan zij doorgaans uit van een vereenvoudigde, frequentie-onafhankelijke beam of een uniforme gedegradeerde resolutie. Realistische instrumenten, zoals MeerKAT, vertonen frequentieafhankelijke beam-distorties, inclusief niet-gladde (oscillerende) variaties in de Full Width at Half Maximum (FWHM). Deze instrumentele effecten koppelen met voorgronden om kunstmatige spectrale kenmerken te creëren die het kosmologische signaal nabootsen, wat de prestaties van standaard scheidingstechnieken verslechtert.

Methodologie
De auteurs evalueren de prestaties van diverse technieken voor voorgrondverwijdering met behulp van gesimuleerde HI intensiteitsmapping-data die het bereik van 900–1400 MHz beslaat, relevant voor MeerKAT en SKA-Mid. De simulaties bevatten realistische kosmologische HI-signalen (gebaseerd op CRIME), galactische voorgronden (gebruikmakend van PySM3-modellen voor synchrotron en free-free), en extragalactische puntbronnen, alles geconvolueerd met telescoopbeams en gecontamineerd met thermische ruis.

Drie beam-modellen werden getest om de realiteitsgraad stapsgewijs te verhogen:

1. Gaussian Degraded: Een uniforme resolutie verkregen door alle kanalen te middelen naar de slechtst denkbare resolutie.
2. Gaussian Evolving: Een beam waarbij de FWHM natuurlijk varieert met de frequentie.
3. Gaussian Oscillating: Een realistisch model dat een oscillerende term in de FWHM-evolutie incorporeert, wat het waargenomen gedrag van de MeerKAT-beam nabootst.

De studie vergelijkt twee categorieën methoden:

• Model-fitting: Polynoom-fitting (6e-orde) en parametrische fitting (het modelleren van free-free en synchrotron met vaste of vrije spectrale indices).
• Blind Source Separation (BSS): Principal Component Analysis (PCA), Fast Independent Component Analysis (FastICA), Generalized Morphological Component Analysis (GMCA), en de voorgestelde SDecGMCA.

SDecGMCA Benadering
Het artikel introduceert de toepassing van SDecGMCA (Spherical DecGMCA) op HI intensiteitsmapping. In tegen tegenstelling tot standaard BSS, dat uitgaat van een lineaire mengeling $X = AS$, adresseert SDecGMCA de frequentieafhankelijke beam door een gezamenlijk deconvolutie- en bron-scheidingsprobleem op te lossen in het sferische harmonische domein. Het model wordt geformuleerd als $\hat{X}_{\ell m}^\nu = \hat{H}_\ell^\nu A_\nu \hat{S}_{\ell m} + \hat{N}_{\ell m}^\nu$, waarbij $\hat{H}_\ell^\nu$ de beam transfer functie is. Het algoritme schat iteratief de bronkaarten en de mengmatrix, terwijl het schaarsheid afdwingt in het wavelet (starlet) domein en datatrouw handhaaft in het harmonische domein.

Belangrijkste Resultaten

• Uniform Beam: Onder het "Gaussian degraded" model presteren alle methoden (fitting en BSS) goed, waarbij ze hoek- en frequentiekrachtspectra met een nauwkeurigheid van beter dan 10% herstellen over de meeste schalen.
• Evolving Beam: De prestaties divergeren. Parametrische fitting faalt volledig. GMCA degradeert bij intermediaire schalen door verstoorde ruimtelijke coherentie. PCA, FastICA en SDecGMCA behouden hun nauwkeurigheid, hoewel BSS-methoden moeite hebben met het volledig ontwarren van schaalafhankelijke distorties bij grote frequentieschalen ($k_\nu < 0.1$).
• Oscillating Beam: Dit is het meest uitdagende scenario. Zowel polynoom- als parametrische fitting falen catastrofaal en produceren zeer bevooroordeelde krachtspectra. Onder de BSS-methoden is SDecGMCA de enige techniek die in staat is het HI-signaal te herstellen met een nauwkeurigheid van $\lesssim 5\%$ bij intermediaire hoekschalen ($10 < \ell < 200$). Cruciaal is dat SDecGMCA een spurieuze spectrale piek bij $k_\nu \sim 0.3$ (toegeschreven aan beam-voorgrondkoppeling) onderdrukt die persistent blijft in PCA, FastICA en GMCA.
• Beperkingen: Alle methoden, inclus inclusief SDecGMCA, vertonen instabiliteit en spurieuze pieken bij hoge multipolen ($\ell > 200$) vanwege de intrinsieke instabiliteit van beam-inversie. Daarnaast wordt een excess bij grote schalen ($\ell < 10$) waargenomen in het oscillerende beam-geval, wat waarschijnlijk voortkomt uit pre-processing artefacten gerelateerd aan het maskeren van heldere HI-bronnen.
• Masking Strategieën: De studie toont aan dat het maskeren van heldere galactische regio's (ofwel via post-processing of door het simuleren van een verzwakking van het galactische vlak) de reconstructiefidelity aanzienlijk verbetert. SDecGMCA profiteert het meest van deze strategie, waarbij een verbeterde reconstructie bij lage multipolen en verdere onderdrukking van de spurieuze piek wordt getoond.
• Hyperparameters: Het aantal bronnen ($n_s$) wordt geïdentificeerd als een kritieke hyperparameter. Een optimaal bereik van $n_s \approx 5-7$ wordt gevonden om een balans te vinden tussen underfitting en overfitting in het geval van de oscillerende beam.

Betekenis en Claims
Het artikel stelt vast dat eenvoudige fitting-methoden en standaard BSS-technieken ontoereikend zijn wanneer realistische, frequentieafhankelijke beam-distorties in overweging worden genomen. Het positioneert SDecGMCA als een bijzonder veelbelovende aanpak voor toekomstige single-dish intensiteitsmapping-surveys omdat het expliciet rekening houdt met de respons van het instrument tijdens het scheidingsproces.

De auteurs beweren dat SDecGMCA robuuste sferische harmonische reconstructies van het HI-krachtspectrum biedt, zelfs onder complexe beam-condities en zonder voorafgaande masking. Ze blijven echter bescheiden over de huidige beperkingen van de methode: de methode vereist nog steeds zorgvuldige afstemming van hyperparameters (met name $n_s$ en regularisatie), en de intrinsieke instabiliteit van beam-inversie beperkt de betrouwbare reconstructie tot $\ell \lesssim 200$. Het artikel concludeert dat hoewel SDecGMCA een belangrijke stap voorwaarts is, toekomstig werk zich moet richten op het integreren van masking-strategieën direct in het framework, het ontwikkelen van data-gedreven diagnostiek voor hyperparameter-selectie, en het verfijnen van thresholding-schema's om de prestaties verder te verbeteren.