Half-wave-plate non idealities propagated to component separated CMB $B$-modes

1. Problématique et Contexte

La détection des modes B primordiaux dans la polarisation du Fond Diffus Cosmologique (CMB) est l'un des objectifs majeurs des expériences actuelles et futures (comme le Simons Observatory et LiteBIRD). Ces signaux, extrêmement faibles, sont la signature des ondes gravitationnelles primordiales générées lors de l'inflation cosmique.

Cependant, l'extraction de ce signal se heurte à plusieurs défis :

• Le bruit instrumental et les systématiques : Les détecteurs et l'optique introduisent des biais.
• Les avant-plans galactiques : L'émission thermique de la poussière et le rayonnement synchrotron masquent le signal cosmologique.
• Les modulateurs de polarisation : Pour supprimer le bruit basse fréquence (bruit $1/f$) et les effets systématiques différentiels, les expériences utilisent des lames demi-onde (HWP) en rotation continue.

Le problème central de l'article réside dans la nature non idéale de ces HWPs. Une HWP idéale ne fonctionne parfaitement qu'à une fréquence unique. Les HWPs réels, souvent constitués de piles de couches (ex: saphir) pour élargir la bande passante, présentent une réponse dépendante de la fréquence. Cette dépendance provoque :

1. Un mélange des paramètres de Stokes ($Q$ et $U$) qui varie avec la fréquence.
2. Une fuite d'intensité vers la polarisation ($I \to Q,U$).
3. Une fuite de polarisation circulaire vers linéaire ($V \to Q,U$).

Si ces effets ne sont pas correctement modélisés lors de la reconstruction des cartes et de la séparation des composantes, ils génèrent des résidus systématiques qui biaisent l'estimation du rapport tenseur-scalaire ($r$), faussant ainsi les contraintes cosmologiques.

2. Méthodologie

Les auteurs ont développé une chaîne de traitement de bout en bout (end-to-end) pour quantifier et atténuer ces effets.

A. Modélisation Instrumentale

Ils utilisent la formalisme de la matrice de Mueller pour décrire la propagation de la lumière polarisée à travers l'instrument.

• Modèle Non Idéal : Une HWP à trois couches (saphir avec revêtements anti-reflets) est modélisée en tenant compte de la dépendance fréquentielle, de l'angle d'incidence et des propriétés optiques des matériaux (indices de réfraction, épaisseurs).
• Données Simulées : Des données temporelles (TOD - Time Ordered Data) sont générées en intégrant la réponse de l'instrument sur toute la bande passante des détecteurs, incluant les effets de la HWP non idéale.

B. Stratégies d'Analyse Comparées

L'étude compare plusieurs approches pour traiter ces données :

1. Approche Standard (Cartes) :

   • Modèle HWP "Efficace" : On approxime la réponse de la HWP par un simple déphasage et un facteur d'amplitude constants par bande de fréquence.
   • Modèle HWP "Empilé" (Stack) : On moyenne la matrice de Mueller de la pile de couches sur la bande passante, en tenant compte du spectre de l'avant-plan (poussière) pour pondérer l'intégration.
   • Ces modèles sont utilisés dans l'étape de cartographie (map-making) pour reconstruire les cartes $Q$ et $U$ par fréquence, suivies d'une séparation de composantes paramétrique standard.

2. Approche Généralisée (Cartes) :

   • Cette méthode reconnaît qu'une matrice de pointage standard ne peut pas reconstruire exactement les cartes $Q$ et $U$ pures en présence de mélange fréquentiel.
   • Elle reconstruit d'abord des cartes de Stokes mixtes (combinaisons linéaires de $Q$ et $U$ pondérées par la réponse instrumentale).
   • Ces cartes mixtes sont ensuite utilisées comme entrée pour une séparation de composantes généralisée, où la matrice de mélange inclut explicitement la réponse instrumentale (matrice de Mueller) en plus des lois spectrales astrophysiques.

3. Approche dans le Domaine Temporel :

   • Extension de la séparation de composantes directement sur les données temporelles (TOD) plutôt que sur les cartes. Cela permet de traiter de manière cohérente le bruit corrélé et les effets dépendants du temps.

C. Outils

L'analyse est implémentée via FURAX, une bibliothèque open-source basée sur JAX, permettant des calculs différentiables, une compilation Just-In-Time (JIT) et une accélération GPU pour gérer la complexité des opérateurs linéaires inverses.

3. Résultats Clés

Les simulations, basées sur une configuration similaire au Simons Observatory (SO), montrent les résultats suivants :

• Impact des modèles simplifiés :

  • Ignorer la dépendance fréquentielle de la HWP entraîne une contamination significative des modes B et un biais important sur $r$.
  • L'utilisation d'un modèle "efficace" (déphasage simple) réduit le biais mais ne l'élimine pas totalement, laissant un résidu correspondant à $r \sim 10^{-2}$.
  • Le modèle "empilé" (moyenné sur la bande) améliore la situation mais reste insuffisant pour atteindre les objectifs de précision ($r \lesssim 10^{-3}$), avec un biais résiduel autour de $r \sim 7 \times 10^{-3}$.

• Efficacité de la séparation généralisée :

  • L'approche généralisée, qui intègre la réponse instrumentale directement dans la matrice de mélange lors de la séparation de composantes, supprime drastiquement les biais systématiques.
  • Elle permet de récupérer une valeur de $r$ compatible avec la valeur d'entrée (zéro dans les simulations de test) avec une incertitude réduite à $r \sim 7 \times 10^{-4}$.

• Validation dans le domaine temporel :

  • L'extension de la méthode au domaine temporel a été validée avec succès. Les résultats obtenus sur les TOD sont cohérents avec ceux obtenus sur les cartes (après rescaling du bruit), confirmant la robustesse de la méthodologie et sa capacité à gérer des modèles de bruit plus réalistes.

• Paramètres spectraux :

  • Les modèles simplifiés biaisent l'estimation de l'indice spectral de la poussière ($\beta_d$) de plusieurs pourcents (ex: 5.8% pour le modèle efficace).
  • La méthode généralisée réduit ce biais à environ 1.3%, et à des niveaux négligeables dans les cas de validation où le modèle correspond parfaitement à la simulation.

4. Contributions et Signification

Contributions principales :

1. Quantification précise : L'article démontre de manière rigoureuse que la dépendance fréquentielle des HWPs non idéales est une source majeure de biais systématique pour les mesures de modes B, même avec des modèles de correction de cartes avancés.
2. Nouveau cadre d'analyse : Introduction d'une méthode de séparation de composantes généralisée qui traite conjointement les effets instrumentaux et astrophysiques, brisant la séparation traditionnelle entre "cartographie" et "séparation".
3. Implémentation logicielle : Développement et utilisation de FURAX pour permettre des analyses complexes et différentiables, facilitant l'optimisation conjointe des paramètres instrumentaux et cosmologiques.

Signification pour la cosmologie :

• Objectifs de précision : Pour atteindre les objectifs des prochaines générations d'expériences (recherche de $r < 0.001$), il est impératif de dépasser les modèles de correction de cartes standards.
• Changement de paradigme : L'article plaide pour une modélisation unifiée où la réponse instrumentale n'est pas traitée comme une étape de prétraitement isolée, mais comme une partie intégrante du modèle d'ajustement global.
• Faisabilité : La démonstration de la faisabilité de la séparation de composantes dans le domaine temporel ouvre la voie à des analyses futures incluant des modèles de bruit atmosphérique et des systématiques temporels complexes, essentiels pour les observations au sol.

En conclusion, ce travail établit que la prise en compte explicite et complète des non-idéalités des HWPs via des cadres d'analyse unifiés est une condition sine qua non pour garantir l'absence de biais dans la détection des ondes gravitationnelles primordiales.