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

Technische Zusammenfassung: Einschränkung dunkler akustischer Oszillationen mit der UV-Leuchtkraftfunktion bei hoher Rotverschiebung

Problemstellung
Während das Standard-$\Lambda$CDM-Modell postuliert, dass Dunkle Materie (DM) kalt und kollisionsfrei ist, sagen verschiedene gut begründete Erweiterungen des dunklen Sektors „Dunkle Akustische Oszillationen“ (DAOs) voraus. Diese entstehen, wenn ein Bruchteil der DM in einem frühen Universum mit einer dunklen Strahlung (DR) interagiert, was zu einer Unterdrückung und gedämpften Oszillationen im Materieleistungsspektrum $P(k)$ führt. Bestehende Einschränkungen aus der kosmischen Mikrowellenhintergrundstrahlung (CMB) sind auf großen Skalen streng, begrenzen jedoch den Anteil der interagierenden DM auf $\lesssim 5\,\%$ für starke DAOs. Umgekehrt sind Einschränkungen aus dem Lyman-$\alpha$-Wald effektiv für schwache DAOs oder vollständig interagierende Modelle, stehen jedoch vor der Herausforderung, die Signatur starker DAOs von Warmen Dunkler Materie (WDM)-Effekten zu unterscheiden. Es besteht ein Bedarf an modellagnostischen Einschränkungen der kleinskaligen Struktur des Universums, insbesondere bei hohen Rotverschiebungen ($z \sim 3-9$), wo die nichtlineare Strukturentwicklung und der spezifische Abdruck von DAOs auf die Halo-Massenfunktion (HMF) untersucht werden können.

Methodik
Die Autoren entwickeln eine Pipeline, um DAO-Parameter mithilfe der hochrotverschobenen Ultravioletten Leuchtkraftfunktion (UVLF) von Galaxien einzuschränken. Die Methodologie besteht aus vier Komponenten:

1. Phänomenologisches Transferfunktionsmodell:
   Um die Rechenkosten für das Lösen von Boltzmann-Gleichungen für jeden Punkt eines Markov-Chain-Monte-Carlo-Scans (MCMC) zu vermeiden, führen die Autoren ein neues vierparametriges phänomenologisches Transferfunktionsmodell $T(k)$ relativ zu $\Lambda$CDM ein. Dieses Modell erfasst die Unterdrückung und die gedämpften Oszillationen von DAOs und erlaubt einen Teilanteil der interagierenden DM von unter eins ($f$). Die Parameter sind:

• $f$: Der Bruchteil der DM, der DAOs unterliegt.
• $k_{\text{peak}}$: Die Wellenzahl des ersten DAO-Peaks.
• $k_{\text{damp}}$: Die Dämpfungsskala der Oszillationen.
• $A$: Die Amplitude der Oszillationen.
  Das Modell wird gegen Leistungsspektren kalibriert, die mit dem modifizierten Boltzmann-Code CLASS-aDM (für Atomare Dunkle Materie) und CLASS (für Interagierende Dunkle Materie) generiert wurden, wobei die Parameter direkt auf physikalische Größen wie den dunklen Schallhorizont und die Diffusionsdämpfungsskala abbilden.

2. Erweiterte Press-Schechter-Formalismus (EPS) und N-Körper-Kalibrierung:
   Um die HMF aus dem linearen Leistungsspektrum vorherzusagen, verwenden die Autoren den erweiterten Press-Schechter-Formalismus. Da Standard-Fensterfunktionen (z. B. Top-Hat im Realraum) bei abgeschnittenen Leistungsspektren versagen, nutzen sie eine „Smooth-k“-Fensterfunktion. Entscheidend ist, dass sie eine Suite von 13 DM-nur N-Körper-Simulationen mit dem GIZMO-Code und MUSIC-Anfangsbedingungen durchführen. Diese Simulationen decken verschiedene Anteile interagierender DM ($f=0,1; 0,5; 1$) sowie DAO-Skalen ($k_{\text{DAO}} = 5, 20, 66, 190 \, h/\text{Mpc}$) ab. Die Simulationen werden verwendet, um die EPS-Parameter ($q, \beta, c$) zu kalibrieren, damit die semi-analytische HMF die Simulationsergebnisse über Rotverschiebungen von $z=4$ bis $12$ genau reproduziert.

3. Halo-Galaxie-Verbindung:
   Die Verbindung zwischen der HMF und der beobachtbaren UVLF wird mithilfe des GALLUMI-Codes modelliert. Diese Pipeline setzt voraus, dass die UV-Leuchtkraft durch Sternentstehungsraten (SFR) angetrieben wird, die im Durchschnitt einer Doppel-Potenzgesetz-Beziehung zur Halo-Akkretionsrate folgen. Das Modell enthält acht astrophysikalische Störparameter (Nuisance-Parameter), welche die Normalisierung, die Pivot-Skala und die Evolution der SFR-Halomasse-Beziehung sowie die Stochastizität in der Magnitude beschreiben. Diese Parameter werden marginalisiert, um die Unsicherheiten in der Hochrotverschiebungs-Galaxienbildung zu berücksichtigen.

4. Datenanalyse:
   Die Autoren kombinieren UVLF-Messungen des Hubble-Weltraumteleskops (HST, sowohl Blank- als auch Linseneffekte-Felder), des James-Webb-Weltraumteleskops (JWST), des Subaru-Teleskops und des Canada-France-Hawaii-Teleskops (CFHT), die Rotverschiebungen von $z=3$ bis $9$ abdecken. Sie führen Bayessche MCMC-Scans und frequentistische Profil-Likelihood-Analysen durch, wobei sie die $\Lambda$CDM-Parameter auf die Planck-2018-Best-Fit-Werte fixieren.

Wesentliche Beiträge

• Neues Phänomenologisches Modell: Einführung eines flexiblen Transferfunktionsmodells, das subdominante Anteile interagierender DM (wie in der Atomaren Dunklen Materie) aufnimmt und die Lücke zwischen spezifischen Teilchenphysikmodellen und allgemeinen Einschränkungen schließt.
• Simulationskalibrierung: Erste Kalibrierung des EPS-Formalismus für DAO-Modelle unter Verwendung eines dedizierten Satzes von N-Körper-Simulationen, was den semi-analytischen Ansatz zur Vorhersage der HMF in Anwesenheit gedämpfter Oszillationen validiert.
• Modellagnostische Einschränkungen: Ableitung von Einschränkungen für den DAO-Parameterraum ($f$ vs. $k_{\text{peak}}$), die unabhängig von spezifischer Teilchenphysik des dunklen Sektors sind und eine einfache Reinterpretation in verschiedenen Modellen (z. B. IDM, Atomare Dunkle Materie) ermöglichen.

Ergebnisse

• Einschränkungen für $f$ und $k_{\text{peak}}$: Die Analyse liefert eine starke Korrelation zwischen dem Anteil der interagierenden DM und der Skala der DAOs.
  • Wenn der Anteil der interagierenden DM signifikant ist ($f \gtrsim 0,07$), muss die Wellenzahl des ersten DAO-Peaks $k_{\text{peak}} \gtrsim 50 \, h/\text{Mpc}$ betragen, um mit den UVLF-Daten konsistent zu bleiben.
  • Umgekehrt, wenn $k_{\text{peak}}$ klein ist (was auf eine frühere Entkopplung oder stärkere Wechselwirkungen hindeutet), ist der Anteil der interagierenden DM auf $f \lesssim 0,07$ begrenzt.
• Vergleich mit anderen Sonden:
  • Für hohe Anteile interagierender DM ($f \sim 1$) verbessern die UVLF-Einschränkungen die durch die CMB abgeleiteten Grenzen um etwa eine Größenordnung in Bezug auf $k_{\text{peak}}$.
  • Für sehr kleine Skalen ($k_{\text{peak}} \lesssim 1 \, h/\text{Mpc}$) bleiben die CMB-Einschränkungen hinsichtlich der Häufigkeit interagierender DM ($f \lesssim 5\,\%$) strenger.
  • Die Einschränkungen für den DM-Baryon-Streuquerschnitt in IDM-Modellen sind um mehr als fünf Größenordnungen stärker als die CMB-Einschränkungen und konkurrenzfähig mit den Werten aus der Häufigkeit von Milchstraßen-Satellitengalaxien.
• Unbeschränkte Parameter: Die Analyse zeigt, dass keine signifikanten Einschränkungen für die DAO-Amplitude ($A$) oder das Verhältnis der Dämpfungsskala ($k_{\text{damp}}/k_{\text{peak}}$) gefunden werden, wenn über andere Parameter marginalisiert wird, da die UVLF-Sensitivität primär durch die anfängliche Unterdrückung des Leistungsspektrums und nicht durch die oszillatorischen Merkmale getrieben wird.

Bedeutung
Die Arbeit zeigt, dass Messungen der UV-Leuchtkraftfunktion bei hoher Rotverschiebung eine leistungsstarke, komplementäre Sonde zu CMB- und Lyman-$\alpha$-Wald-Daten darstellen, um die Physik der Dunklen Materie auf kleinen Skalen einzuschränken. Durch die erfolgreiche Kalibrierung des EPS-Formalismus gegen N-Körper-Simulationen und die Anwendung eines flexiblen Halo-Galaxie-Verbindungsmodells etablieren die Autoren robuste Grenzwerte für den Anteil der Dunklen Materie, der akustische Oszillationen unterliegen kann. Diese Ergebnisse verschieben den Entkopplungszeitpunkt der interagierenden Dunklen Materie zu früheren Zeiten (höheres $k_{\text{peak}}$) oder erfordern, dass der interagierende Anteil subdominant ist ($f < 7\,\%$) für spezifische Interaktionsskalen. Die Arbeit unterstreicht die Nützlichkeit von Beobachtungen bei hoher Rotverschiebung, um zwischen verschiedenen Paradigmen des dunklen Sektors zu unterscheiden, und legt den Grundstein für zukünftige Einschränkungen durch kommende Daten des Vera C. Rubin Observatory, des Nancy Grace Roman Space Telescope und des Square Kilometre Array.