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

Résumé technique : Retrait du premier plan dans la cartographie d'intensité HI sous des distorsions de faisceau dépendantes de la fréquence

Énoncé du problème
La cartographie d'intensité de l'hydrogène neutre (HI) est une sonde primaire pour la cosmologie de l'époque post-réionisation, offrant une vue tomographique de l'Univers. Cependant, l'extraction du faible signal HI est sévèrement entravée par les avant-plans astrophysiques (rayonnement synchrotron, rayonnement libre et sources ponctuelles) qui sont jusqu'à cinq ordres de grandeur plus brillants. Bien que la séparation aveugle de sources (BSS) et les méthodes d'ajustement de modèles aient montré des résultats prometteurs, elles supposent généralement un faisceau simplifié, indépendant de la fréquence, ou une résolution dégradée uniforme. Les instruments réalistes, tels que MeerKAT, présentent des distorsions de faisceau dépendantes de la fréquence, incluant des variations non lisses (oscillatoires) de la largeur à mi-hauteur (FWHM). Ces effets instrumentaux se couplent aux avant-plans pour créer des caractéristiques spectrales artificielles qui imitent le signal cosmologique, dégradant ainsi les performances des techniques de séparation standard.

Méthodologie
Les auteurs évaluent la performance de diverses techniques de retrait d'avant-plan en utilisant des données simulées de cartographie d'intensité HI couvrant la plage 900–1400 MHz, pertinente pour MeerKAT et SKA-Mid. Les simulations incluent des signaux HI cosmologiques réalistes (basés sur CRIME), des avant-plans galactiques (utilisant des modèles PySM3 pour le rayonnement synchrotron et libre), et des sources ponctuelles extragalactiques, tous convolués avec des faisceaux de télescope et contaminés par le bruit thermique.

Trois modèles de faisceau ont été testés pour augmenter progressivement le réalisme :

1. Gaussien dégradé : Une résolution uniforme obtenue en lissant tous les canaux à la pire résolution.
2. Gaussien évolutif : Un faisceau où la FWHM varie naturellement avec la fréquence.
3. Gaussien oscillant : Un modèle réaliste incorporant un terme oscillatoire dans l'évolution de la FWHM, imitant le comportement observé du faisceau de MeerKAT.

L'étude compare deux catégories de méthodes :

• Ajustement de modèle (Model-fitting) : Ajustement polynomial (ordre 6) et ajustement paramétrique (modélisation du rayonnement libre et du rayonnement synchrotron avec des indices spectraux fixes ou libres).
• Séparation aveugle de sources (BSS) : Analyse en composantes principales (PCA), Analyse de composantes indépendantes rapide (FastICA), Analyse de composantes morphologiques généralisée (GMCA) et la méthode proposée SDecGMCA.

Approche SDecGMCA
Le document introduit l'application de SDecGMCA (Spherical DecGMCA) à la cartographie d'intensité HI. Contrairement à la BSS standard qui suppose un mélange linéaire $X = AS$, SDecGMCA traite le faisceau dépendant de la fréquence en résolvant un problème conjoint de déconvolution et de séparation de sources dans le domaine des harmoniques sphériques. Le modèle est formulé comme $\hat{X}_{\ell m}^\nu = \hat{H}_\ell^\nu A_\nu \hat{S}_{\ell m} + \hat{N}_{\ell m}^\nu$, où $\hat{H}_\ell^\nu$ est la fonction de transfert du faisceau. L'algorithme estime de manière itérative les cartes de sources et la matrice de mélange tout en imposant la parcimonie dans le domaine des ondelettes (starlet) et la fidélité des données dans l'espace harmonique.

Résultats clés

• Faisceau uniforme : Sous le modèle « Gaussien dégradé », toutes les méthodes (ajustement et BSS) performent bien, récupérant les spectres de puissance angulaire et fréquentielle avec une précision supérieure à 10 % sur la plupart des échelles.
• Faisceau évolutif : La performance diverge. L'ajustement paramétrique échoue complètement. La GMCA se dégrade aux échelles intermédiaires en raison d'une rupture de la cohérence spatiale. La PCA, la FastICA et la SDecGMCA maintiennent leur précision, bien que les méthodes BSS peinent à démêler pleinement les distorsions dépendantes de l'échelle aux grandes échelles de fréquence ($k_\nu < 0,1$).
• Faisceau oscillant : C'est le scénario le plus difficile. L'ajustement polynomial et l'ajustement paramétrique échouent de manière catastrophique, produisant des spectres de puissance très biaisés. Parmi les méthodes BSS, la SDecGMCA est la seule technique capable de récupérer le signal HI avec une précision $\lesssim 5\,\%$ aux échelles angulaires intermédiaires ($10 < \ell < 200$). Crucialement, la SDecGMCA supprime un pic spectral parasite à $k_\nu \sim 0,3$ (attribué au couplage faisceau-avant-plan) qui persiste dans la PCA, la FastICA et la GMCA.
• Limites : Toutes les méthodes, y compris la SDecGMCA, présentent une instabilité et des pics parasites aux multipôles élevés ($\ell > 200$) en raison de l'instabilité intrinsèque de l'inversion du faisceau. De plus, un excès de grande échelle ($\ell < 10$) est observé dans le cas du faisceau oscillant, provenant probablement d'artefacts de prétraitement liés au masquage des sources HI brillantes.
• Stratégies de masquage : L'étude démontre que le masquage des régions galactiques brillantes (soit via le post-traitement, soit en simulant l'atténuation du plan galactique) améliore significativement la fidélité de la reconstruction. La SDecGMCA bénéficie le plus de cette stratégie, montrant une meilleure récupération aux faibles multipôles et une suppression accrue du pic parasite.
• Hyperparamètres : Le nombre de sources ($n_s$) est identifié comme un hyperparamètre critique. Une plage optimale de $n_s \approx 5-7$ est trouvée pour équilibrer le sous-ajustement et le sur-ajustement dans le cas du faisceau oscillant.

Signification et affirmations
L'article établit que les méthodes d'ajustement simples et les techniques BSS standard sont insuffisantes lorsque des distorsions de faisceau réalistes et dépendantes de la fréquence sont considérées. Il positionne la SDecGMCA comme une approche particulièrement prometteuse pour les futures campagnes de cartographie d'intensité à antenne unique car elle prend explicitement en compte la réponse de l'instrument pendant le processus de séparation.

Les auteurs affirment que la SDecGMCA fournit des reconstructions robustes en harmoniques sphériques du spectre de puissance HI, même sous des conditions de faisceau complexes et sans pré-masquage. Cependant, ils restent modestes quant à ses limites actuelles : la méthode nécessite toujours un réglage minutieux des hyperparamètres (spécifiquement $n_s$ et la régularisation), et l'instabilité intrinsèque de l'inversion du faisceau limite la reconstruction fiable à $\ell \lesssim 200$. L'article conclut que, bien que la SDecGMCA soit une étape importante, les travaux futurs doivent se concentrer sur l'intégration directe des stratégies de masquage dans le cadre de travail, le développement de diagnostics basés sur les données pour la sélection des hyperparamètres, et l'affinement des schémas de seuillage pour améliorer davantage les performances.