Forecast Cosmological Constraints with the 1D Wavelet Scattering Transform and the Lyman-$α$ forest

Énoncé du Problème
La forêt Lyman-$\alpha$ constitue une sonde critique pour les structures à petite échelle, offrant des perspectives sur la matière noire, la masse des neutrinos, la physique de l'inflation et l'histoire thermique du milieu intergalactique (IGM). Cependant, l'extraction d'informations cosmologiques à partir de ce régime non linéaire est complexe. Les analyses actuelles reposant sur le spectre de puissance du flux 1D deviennent dominées par les erreurs systématiques dans les grands relevés comme SDSS et DESI. De plus, les statistiques d'ordre supérieur traditionnelles, telles que le bispectre, souffrent d'une augmentation de la variance et d'une robustesse réduite face aux valeurs aberrantes dans des distributions hautement non gaussiennes. Il existe un besoin démontré pour de nouvelles statistiques de résumé capables d'extraire l'information non gaussienne du champ de densité non linéaire sans les exigences d'entraînement ou les problèmes d'interprétabilité associés aux approches d'apprentissage profond.

Méthodologie
Les auteurs proposent l'utilisation de la Transformée de Diffusion par Ondelettes (WST) 1D comme statistique de résumé pour la cosmologie de la forêt Lyman-$\alpha$. La WST est une technique de traitement du signal déterministe et mathématiquement robuste qui convolve le champ de flux d'entrée avec des ondelettes à valeurs complexes à différentes échelles, suivi d'opérations de module et d'un filtrage passe-bas. Ce processus sépare l'information à travers les échelles et extrait des statistiques d'ordre supérieur (jusqu'au second ordre dans cette étude) qui sont stables contre les petites déformations et le bruit.

Pour prévoir le pouvoir de contrainte de la WST, les auteurs ont utilisé :

• Simulations : Une suite de 10 simulations hydrodynamiques (utilisant MP-Gadget) contenant $2 \times 256^3$ particules dans une boîte de $40 \text{ Mpc}/h$. À partir de celles-ci, 32 000 spectres d'absorption Lyman-$\alpha$ ont été générés aux redshifts $z=2, 3, 4$.
• Paramètres : L'analyse s'est concentrée sur quatre paramètres cosmologiques et thermiques : l'indice spectral de la puissance primaire ($n_P$) et son amplitude ($A_P$), ainsi que les paramètres thermiques de l'IGM ($H_A$ et $H_S$).
• Modélisation du bruit : Un bruit observationnel réaliste a été simulé en ajoutant un bruit gaussien à chaque spectre, le rapport continuum-sur-bruit (CNR) étant tiré d'une distribution log-normale ajustée aux données SDSS DR16Q.
• Cadre statistique : Une matrice d'information de Fisher a été construite pour estimer les contraintes de paramètres. La matrice de covariance a été estimée à l'aide de 2 000 réalisations de bruit aléatoire. Les auteurs ont comparé l'information de Fisher dérivée des coefficients de la WST par rapport à celle du spectre de puissance du flux 1D standard.

Contributions Clés et Résultats
L'étude démontre que la WST capture une information cosmologique significative absente du spectre de puissance du flux.

• Contenu d'information : Les coefficients d'ondelettes du second ordre, qui sondent l'information non gaussienne, contribuent environ 50 % d'information supplémentaire à la matrice de Fisher par rapport aux seuls coefficients du premier ordre.
• Contraintes de paramètres : La WST fournit des contraintes nettement plus serrées sur les paramètres cosmologiques que le spectre de puissance. Plus précisément, les incertitudes postérieures estimées ($\sigma_X$) pour les paramètres $n_P$, $A_P$, $H_S$ et $H_A$ sont améliorées d'un facteur de 30 à 60 en utilisant les coefficients combinés de la WST.
• Réduction de dimensionnalité : Les auteurs estiment que la WST réduit l'espace des paramètres 4D d'un facteur de $10^6$ par rapport au spectre de puissance.
• Sensibilité : Les coefficients de la WST sont moins sensibles au flux moyen ($\tau_{\text{eff}}$) que le spectre de puissance (changeant de $\sim 2,5 \%$ pour un changement de flux de 10 %, contre $\sim 7 \%$ pour le spectre de puissance), suggérant une réduction potentielle de la dégénérescence entre $\tau_{\text{eff}}$ et l'amplitude primaire $A_P$.
• Dépendance au redshift : Le contenu d'information culmine à $z=3$ pour la plupart des paramètres, tandis que le paramètre thermique $H_A$ est le plus contraint à $z=2$.

Signification et Revendications
L'article affirme que la Transformée de Diffusion par Ondelettes 1D offre une approche nouvelle et puissante pour extraire l'information cosmologique de la forêt Lyman-$\alpha$, surpassant le spectre de puissance du flux traditionnel d'un ordre de grandeur ou plus en termes de pouvoir de contrainte. Les auteurs postulent que si ces améliorations se concrétisent dans les données à venir de l'instrument de spectroscopie de l'énergie sombre (DESI), la WST pourrait permettre :

• Des contraintes sur la course de l'inflation ($\alpha_s$) suffisantes pour tester les modèles inflationnaires courants prédisant $\alpha_s = -6 \times 10^{-4}$.
• Des contraintes sur la masse des neutrinos suffisamment améliorées pour permettre une détection à $5\sigma$ de la masse minimale des neutrinos.

Les auteurs restent modestes quant aux limites actuelles de leur prévision, notant que l'approche de la matrice de Fisher néglige les dégénérescences de paramètres et que la petite taille de la boîte de simulation conduit probablement à des prévisions pessimistes pour l'indice spectral $n_P$. Ils reconnaissent qu'une analyse de vraisemblance complète et un traitement robuste des erreurs systématiques (tels que l'ajustement du continuum) sont des étapes futures nécessaires avant d'appliquer ces méthodes aux données réelles de DESI. Cependant, ce travail établit que la WST préserve nettement plus d'informations que le spectre de puissance, ouvrant la possibilité d'améliorer les contraintes de paramètres de plusieurs ordres de grandeur.