Beyond the Power Spectrum: A New Framework for Non-Stationary Fields with Applications to Light-Cone Effects in Line Intensity Mapping

1. Le Problème : L'Effet Cône de Lumière et les Limites du Spectre de Puissance

Les relevés cosmologiques modernes (comme DESI, Euclid, ou les expériences de cartographie d'intensité de raies comme HERA) couvrent des volumes extrêmement vastes, atteignant des échelles de plusieurs gigaparsecs. Une hypothèse fondamentale de la cosmologie standard est l'invariance par translation (stationnarité) des champs cosmologiques, ce qui permet d'utiliser le spectre de puissance comme statistique résumée optimale.

Cependant, sur de tels volumes, l'observation le long de la ligne de visée (LoS) englobe une période significative de l'histoire cosmique. Cela introduit l'effet cône de lumière : les propriétés statistiques du champ observé évoluent le long de la ligne de visée (par exemple, durant l'Ère de la Réionisation, le gaz passe d'un état neutre à un état ionisé).

• Conséquence : Le champ n'est plus stationnaire le long de la LoS.
• Limitation du spectre de puissance : Le spectre de puissance suppose que les modes de Fourier sont non corrélés (la matrice de covariance est diagonale). En présence de l'effet cône de lumière, les modes de Fourier le long de la LoS deviennent corrélés. Le spectre de puissance perd ces informations de covariance non diagonales, rendant l'estimation des paramètres astrophysiques sous-optimale.

2. Méthodologie : Une Nouvelle Base Non-Fourier

Les auteurs proposent un cadre mathématique pour capturer ces corrélations en développant une nouvelle base de décomposition inspirée de la décomposition en valeurs propres (eigen decomposition).

A. Construction de la Matrice de Covariance

Au lieu de supposer l'invariance par translation, les auteurs calculent la matrice de covariance le long de la ligne de visée ($Cov_{LS}$) en utilisant des simulations cosmologiques (réalisations de la température de brillance du signal 21 cm).

• L'invariance par translation est préservée dans les directions angulaires (perpendiculaires à la LoS), où l'on continue d'utiliser une base harmonique (sphérique ou de Fourier).
• Le long de la LoS, la covariance est calculée en moyennant sur de nombreuses lignes de visée.

B. Décomposition en Valeurs Propres et Nouvelle Base

La matrice de covariance $Cov_{LS}$ est diagonalisée numériquement.

• Résultat : Les vecteurs propres obtenus forment une nouvelle base optimale qui diagonalise la covariance du champ évolutif.
• Caractéristiques : Contrairement aux modes de Fourier purs (sinusoïdes constantes), ces nouveaux vecteurs propres sont des sinusoïdes modulées dont l'amplitude varie en fonction du redshift (z). Ils capturent naturellement l'évolution physique du signal (par exemple, la transition absorption/émission du signal 21 cm).

C. La Statistique Résumée Hybride ($M(k_\perp, i)$)

Pour éviter les problèmes de convergence numérique et le surajustement (overfitting) liés à l'estimation de covariance sur un nombre fini de lignes de visée, les auteurs proposent une base hybride :

1. On conserve les premiers vecteurs propres (ceux avec les plus grandes valeurs propres) qui convergent rapidement et capturent la majeure partie de la variance physique.
2. Le reste de la base est complété par des modes de Fourier analytiques standards pour les modes d'ordre supérieur (dominés par le bruit).
3. Un processus d'orthogonalisation (Gram-Schmidt) assure l'orthogonalité stricte de cette base hybride.

La nouvelle statistique résumée, notée $M(k_\perp, i)$, est définie comme le carré du module des coefficients d'expansion de la carte dans cette nouvelle base, binnés sur les anneaux de l'espace de Fourier transversal ($k_\perp$).

3. Contributions Clés

1. Formalisme Mathématique : Dérivation rigoureuse d'une base optimale pour les champs non stationnaires le long de la ligne de visée, généralisable à tout champ cosmologique évolutif.
2. Statistique Résumée Optimisée : Introduction de $M(k_\perp, i)$, une statistique qui préserve les informations de covariance entre les modes de Fourier le long de la LoS, contrairement au spectre de puissance 2D traditionnel $P(k_\perp, k_\parallel)$.
3. Stratégie Hybride : Développement d'une méthode pratique combinant décomposition en valeurs propres et modes de Fourier pour rendre l'estimation robuste face aux limitations des données réelles (nombre fini de lignes de visée).
4. Application à HERA : Mise en œuvre concrète pour une expérience de cartographie d'intensité de raie (Line Intensity Mapping - LIM) réaliste, incluant le bruit thermique et les stratégies d'évitement des avant-plans (foregrounds).

4. Résultats

Les auteurs ont réalisé une prévision de Fisher (Fisher forecast) pour comparer leur nouvelle méthode avec deux approches Fourier classiques :

• Spectre de Puissance 2D Complet (Full Lightcone) : Transformée de Fourier sur tout le volume cône de lumière.
• Spectre de Puissance 2D par Tranches (Slicing) : Division du volume en sous-boîtes pour supposer la stationnarité locale.

Constats principaux :

• Gain d'Information : La méthode hybride $M(k_\perp, i)$ fournit des contraintes sur les paramètres astrophysiques de l'Ère de la Réionisation (EoR) environ 30 % plus serrées que les meilleures méthodes Fourier existantes.
• Réduction des Dégénérescences : L'amélioration provient principalement d'une réduction des dégénérescences entre les paramètres (les termes hors-diagonale de la matrice de Fisher sont réduits), permettant une inversion plus précise pour obtenir les incertitudes sur les paramètres.
• Robustesse : La base reste qualitativement similaire même avec des variations de paramètres astrophysiques (ex: variation de 10 % de l'efficacité de fuite des photons ionisants), ce qui rend la méthode applicable en pratique (on peut estimer les paramètres grossièrement avec un spectre de puissance classique, puis affiner avec la nouvelle base).

5. Signification et Perspectives

Ce travail marque une avancée significative pour l'analyse des données cosmologiques de nouvelle génération :

• Optimisation de la Compression de Données : Il démontre qu'il est possible de compresser l'information cosmique de manière plus efficace que le spectre de puissance sans passer par des méthodes de corrélation à deux points brutes (qui sont trop coûteuses en calcul, $O(N^2)$). La nouvelle statistique reste compacte ($O(N)$).
• Exploitation des Grands Volumes : Avec des relevés couvrant des volumes gigantesques, l'effet cône de lumière ne peut plus être ignoré. Cette méthode permet d'exploiter pleinement l'information contenue dans l'évolution le long de la ligne de visée.
• Applicabilité Future : Bien que testé sur l'Ère de la Réionisation (21 cm), le cadre est généralisable à d'autres époques (Aube Cosmique, Âges Sombres) et à d'autres sondages (galaxies, CMB) où l'évolution cosmique le long de la ligne de visée est significative.

En conclusion, les auteurs proposent un outil statistique robuste qui dépasse les limitations du spectre de puissance traditionnel, offrant une voie prometteuse pour extraire la richesse informationnelle maximale des futurs relevés cosmologiques à grand volume.