Constraining Dark Acoustic Oscillations with the High-Redshift UV Luminosity Function

Résumé Technique : Contraindre les Oscillations Acoustiques Sombres avec la Fonction de Luminosité UV à haut décalage vers le rouge

Énoncé du Problème
Bien que le modèle standard $\Lambda$CDM postule que la matière noire (DM) est froide et sans collision, diverses extensions bien motivées du secteur sombre prédisent des « Oscillations Acoustiques Sombres » (DAO). Celles-ci surviennent lorsqu'une fraction de la DM interagit avec une espèce de rayonnement sombre (DR) dans l'univers primordial, entraînant une suppression et des oscillations amorties dans le spectre de puissance de la matière $P(k)$ aux petites échelles. Les contraintes existantes issues du Fond de Rayonnement Cosmique Micro-onde (CMB) sont strictes sur les grandes échelles mais limitent la fraction de DM interagissante à $\lesssim 5\%$ pour des DAO fortes. Inversement, les contraintes du forêt Lyman-$\alpha$ sont efficaces pour les DAO faibles ou les modèles entièrement interagissants, mais font face à des difficultés pour distinguer les signatures de DAO fortes des effets de la Matière Noire Tiède (WDM). Il existe un besoin de contraintes agnostiques vis-à-vis du modèle pour la structure à petite échelle de l'univers, particulièrement aux décalages vers le rouge élevés ($z \sim 3-9$), où l'évolution non linéaire des structures et l'empreinte spécifique des DAO sur la fonction de masse des halos (HMF) peuvent être sondées.

Méthodologie
Les auteurs développent un pipeline pour contraindre les paramètres des DAO en utilisant la fonction de luminosité UV (UVLF) de haute fréquence des galaxies. La méthodologie se compose de quatre composantes :

1. Modèle de Fonction de Transfert Phénoménologique :
   Pour éviter le coût computationnel de la résolution des équations de Boltzmann pour chaque point d'un balayage par chaîne de Markov Monte Carlo (MCMC), les auteurs introduisent un nouveau modèle de fonction de transfert phénoménologique à quatre paramètres $T(k)$ par rapport à $\Lambda$CDM. Ce modèle capture la suppression et les oscillations amorties des DAO et permet une fraction sub-unitaire de DM interagissante ($f$). Les paramètres sont :

• $f$ : la fraction de la DM subissant des DAO.
• $k_{\text{peak}}$ : le nombre d'onde du premier pic de DAO.
• $k_{\text{damp}}$ : l'échelle d'amortissement des oscillations.
• $A$ : l'amplitude des oscillations.
  Le modèle est calibré par rapport aux spectres de puissance générés par le code de Boltzmann modifié CLASS-aDM (pour l'Atomic Dark Matter) et CLASS (pour l'Interacting Dark Matter), démontant que les paramètres se rapportent directement à des quantités physiques telles que l'horizon sonore sombre et l'échelle de diffusion d'amortissement.

2. Formalisme de Press-Schechter Étendu (EPS) et Calibration N-corps :
   Pour prédire la HMF à partir du spectre de puissance linéaire, les auteurs utilisent le formalisme de Press-Schechter étendu. Reconnaissant que les fonctions de fenêtre standard (par exemple, le top-hat dans l'espace réel) échouent pour les spectres de puissance tronqués, ils utilisent une fonction de fenêtre « smooth-k ». Crucialement, ils effectuent une suite de 13 simulations N-corps de matière noire seule avec le code GIZMO utilisant des conditions initiales MUSIC. Ces simulations couvrent diverses fractions de DM interagissante ($f=0.1, 0.5, 1$) et diverses échelles de DAO ($k_{\text{DAO}} = 5, 20, 66, 190 \, h/\text{Mpc}$). Les simulations sont utilisées pour calibrer les paramètres EPS ($q, \beta, c$) afin de s'assurer que la HMF semi-analytique reproduit fidèlement les résultats des simulations à travers les redshifts $z=4$ à $12$.

3. Connexion Halo-Galaxie :
   La connexion entre la HMF et l'observable UVLF est modélisée à l'aide du code GALLUMI. Ce pipeline suppose que la luminosité UV est pilotée par les taux de formation stellaire (SFR) qui suivent, en moyenne, une relation en loi de puissance double avec le taux d'accrétion du halo. Le modèle inclut huit paramètres de nuisance astrophysiques (décrivant la normalisation, l'échelle de pivot et l'évolution de la relation SFR-masse de halo, ainsi que la stochasticité de la magnitude). Ces paramètres sont marginalisés pour tenir compte des incertitudes de la physique de la formation des galaxies à haut redshift.

4. Analyse des Données :
   Les auteurs combinent les mesures de l'UVLF provenant du télescope spatial Hubble (HST, champs vierges et lentillés), du télescope spatial James Webb (JWST), du télescope Subaru et du Canada-France-Hawaii Telescope (CFHT) couvrant les redshifts $z=3$ à $9$. Ils effectuent des balayages MCMC bayésiens et des analyses de vraisemblance de profil fréquentistes pour contraindre les paramètres des DAO, en fixant les paramètres $\Lambda$CDM aux valeurs de meilleur ajustement de Planck 2018.

Contributions Clés

• Nouveau Modèle Phénoménologique : Introduction d'un modèle de fonction de transfert flexible qui accommode des fractions de DM interagissante sub-dominantes, comblant le fossé entre les modèles de physique des particules spécifiques (comme l'Atomic DM) et les contraintes générales.
• Calibration par Simulation : Première calibration du formalisme EPS pour les modèles de DAO à l'aide d'un ensemble dédié de simulations N-corps, validant l'approche semi-analytique pour prédire la HMF en présence d'oscillations amorties.
• Contraintes Agnostiques du Modèle : Dérivation de contraintes sur l'espace des paramètres des DAO ($f$ vs. $k_{\text{peak}}$) indépendantes de la physique des particules spécifique du secteur sombre, permettant une réinterprétation facile dans divers modèles (par exemple, IDM, Atomic DM).

Résultats

• Contraintes sur $f$ et $k_{\text{peak}}$ : L'analyse produit une forte corrélation entre la fraction de DM interagissante et l'échelle des DAO.
  • Si la fraction de DM interagissante est significative ($f \gtrsim 0.07$), le nombre d'onde du premier pic de DAO doit être $k_{\text{peak}} \gtrsim 50 \, h/\text{Mpc}$ pour rester cohérent avec les données de l'UVLF.
  • Inversement, si $k_{\text{peak}}$ est petit (indiquant un découplage plus précoce ou des interactions plus fortes), la fraction de DM interagissante est contrainte à $f \lesssim 0.07$.
• Comparaison avec d'autres Sondes :
  • Pour des fractions élevées de DM interagissante ($f \sim 1$), les contraintes de l'UVLF améliorent les limites dérivées du CMB d'environ un ordre de grandeur en $k_{\text{peak}}$.
  • Pour des échelles très petites ($k_{\text{peak}} \lesssim 1 \, h/\text{Mpc}$), les contraintes du CMB restent plus strictes sur l'abondance de la DM interagissante ($f \lesssim 5\%$).
  • Les contraintes sur la section efficace de diffusion DM-baryon dans les modèles de Matière Noire Interagissante (IDM) sont trouvées être plus de cinq ordres de grandeur plus fortes que les contraintes du CMB et compétitives avec les limites dérivées des abondances de satellites de la Voie Lactée.
• Paramètres Non Contraints : L'analyse ne trouve aucune contrainte significative sur l'amplitude des DAO ($A$) ou sur le rapport de l'échelle d'amortissement ($k_{\text{damp}}/k_{\text{peak}}$) lors de la marginalisation sur les autres paramètres, car la sensibilité de l'UVLF est principalement pilotée par la suppression initiale du spectre de puissance plutôt que par les caractéristiques oscillatoires.

Signification
Ce papier démontre que les mesures de la fonction de luminosité UV à haut redshift constituent une sonde puissante et complémentaire aux données du CMB et de la forêt Lyman-$\alpha$ pour contraindre la physique de la matière noire à petite échelle. En calibrant avec succès le formalisme EPS par rapport aux simulations N-corps et en employant un modèle flexible de connexion halo-galaxie, les auteurs établissent des limites robustes sur la fraction de la matière noire pouvant subir des oscillations acoustiques. Ces résultats poussent l'époque de découplage de la matière noire interagissante vers des temps plus précoces (plus hauts $k_{\text{peak}}$) ou exigent que la fraction interagissante soit sub-dominante ($f < 7\%$) pour des échelles d'interaction spécifiques. Ce travail souligne l'utilité des observables à haut redshift pour distinguer entre différents paradigmes du secteur sombre et pose les bases pour les futures contraintes utilisant les données à venir de l'Observatoire Vera Rubin, du télescope spatial Nancy Grace Roman et de l'SKA.