Cosmological constraints from the DESI DR1 joint power spectrum and bispectrum analysis

Titel: Kosmologische beperkingen uit de gezamenlijke analyse van het vermogensspectrum en het bispectrum van DESI DR1

Auteurs: S. Novell-Masot et al. (DESI Collaboration)
Publicatie: Voorbereid voor indiening bij JCAP (arXiv:2603.19356v1)

---

1. Het Probleem en de Context

De kosmologie bevindt zich in een fase van steeds grotere precisie, waarbij het standaardmodel $\Lambda$CDM (Lambda-Koude Donkere Materie) tot nu toe uitstekend presteert. Er blijven echter uitdagingen bestaan, zoals de aard van donkere energie en donkere materie, en de aanhoudende spanningen tussen vroege- en late-universe metingen (bijvoorbeeld de $H_0$-spanning).

Traditionele analyses van grote roodverschuivingsonderzoeken (zoals DESI) focussen vaak op twee-punts statistieken (het vermogensspectrum, $P(k)$) en gebruiken gecomprimeerde observabelen zoals Baryon Acoustic Oscillations (BAO) en Redshift-Space Distortions (RSD). Hoewel deze methode robuust is, verwaarloost ze informatie in de "breedband-vorm" van het spectrum en kan het bepaalde parameterdegeneraties niet volledig doorbreken.

Het centrale probleem: Hoe kunnen we de volledige informatie uit de clustering van sterrenstelsels benutten, inclusief hogere-orde statistieken zoals het bispectrum (drie-punts correlatie), om de kosmologische beperkingen te verscherpen zonder afhankelijk te zijn van complexe, modelafhankelijke "full-shape" fittingen die zware aannames vereisen?

2. Methodologie

De auteurs gebruiken de ShapeFit-methode, een tussenweg tussen standaard BAO-RSD-compressie en volledige modellering van het vermogensspectrum.

• Data: De analyse maakt gebruik van DESI Data Release 1 (DR1).
  • Gebruikt: Luminous Red Galaxies (LRG) en Quasars (QSO) voor de gezamenlijke analyse van het vermogensspectrum en het bispectrum.
  • Aanvullend: Reconstructie-BAO-resultaten voor Bright Galaxy Survey (BGS), Emission Line Galaxies (ELG) en het Lyman-$\alpha$ bos.
• ShapeFit Compressie: In plaats van direct een kosmologisch model op de data te fitten, worden de data gecomprimeerd naar een set tussenparameters:
  • $\alpha_{iso}$ en $\alpha_{AP}$: Schalingsfactoren voor BAO (isotroop en anisotroop).
  • $f(z)$ en $\sigma_{s8}(z)$: Groeisnelheid van structuur en amplitude van fluctuaties.
  • $m(z) + n(z)$: Een gecombineerde parameter die de helling van het vermogensspectrum rond de materie-straling gelijkheidsschaal modelleert.
• Bispectrum Integratie: Voor het eerst worden de ShapeFit-constraints voor LRG en QSO afgeleid uit een gezamenlijke analyse van het vermogensspectrum (multipolen) en het bispectrum (monopool).
  • Het bispectrum helpt de degeneratie tussen $f$ en $\sigma_8$ te doorbreken, wat niet mogelijk is met alleen het vermogensspectrum.
  • De analyse gebruikt een maximum $k_{max}$ van $0.15 \, h \text{Mpc}^{-1}$ voor het vermogensspectrum en $0.12 \, h \text{Mpc}^{-1}$ voor het bispectrum.
• Inferentie: Een Markov Chain Monte Carlo (MCMC) algoritme wordt gebruikt om de posterior-verdeling van kosmologische parameters te schatten, gebaseerd op de gecomprimeerde parameters en hun covariantiematrix.

3. Belangrijkste Bijdragen

1. Eerste ShapeFit-resultaten met Bispectrum: Dit is het eerste werk dat ShapeFit-kosmologische inferentie toepast op het bispectrum van DESI DR1.
2. Verscherpte Beperkingen: Het tonen aan dat het toevoegen van het bispectrum de onzekerheid op de amplitude van fluctuaties ($\ln(A_s \times 10^{10})$) met ongeveer 20% verkleint ten opzichte van een analyse met alleen het vermogensspectrum.
3. Robuustheidstest: Het valideren dat de ShapeFit-methode, zelfs met de toevoeging van het bispectrum, resultaten oplevert die consistent zijn met officiële DESI-analyses en met "Full Modelling" benaderingen.
4. Uitgebreide Modellen: Toepassing van de methode op uitgebreide modellen, waaronder variabele donkere energie ($w_0w_a$) en de som van neutrino-massa's ($\sum m_\nu$), zonder sound-horizon-afhankelijkheid.

4. Resultaten

A. Baseline $\Lambda$CDM Model

• De resultaten zijn breed in overeenstemming met eerdere DESI DR1-analyses en Planck 2018.
• Materiedichtheid: $\Omega_m = 0.310 \pm 0.012$
• Hubble-constante: $H_0 = 68.92 \pm 0.97 \, \text{km s}^{-1}\text{Mpc}^{-1}$
• Amplitude: $\ln(A_s \times 10^{10}) = 3.125 \pm 0.079$ (met bispectrum).
• Het toevoegen van het bispectrum verschuift de posteriors lichtjes richting de Planck-resultaten en vermindert de foutmarges aanzienlijk.

B. Sound-Horizon Marginalisatie (Vrij van $r_d$)

• De auteurs testen een scenario waarbij de geluidshorizon ($r_d$) niet als een vaste waarde uit de vroege fysica wordt aangenomen, maar als een vrije parameter.
• Dit leidt tot een "sound-horizon-free" bepaling van $H_0$.
• Resultaat: $H_0 = 70.9 \pm 2.9 \, \text{km s}^{-1}\text{Mpc}^{-1}$.
• De onzekerheid neemt toe (van ~1.3% naar ~5%), maar het resultaat is consistent met zowel Planck als lokale afstandsladder-metingen, wat de robuustheid van de methode bevestigt.

C. Evoluerende Donkere Energie ($w_0w_a$)

• Gebruikmakend van de CPL-parametrizatie ($w(a) = w_0 + w_a(1-a)$).
• De resultaten tonen geen sterke afwijking van het kosmologische constant ($\Lambda$).
• Resultaat: $w_0 = -0.62 \pm 0.29$, $w_a = -1.34 \pm 0.82$.
• De waarden liggen binnen $2\sigma$ van $w_0 = -1, w_a = 0$. Er is geen significant bewijs voor dynamische donkere energie.

D. Neutrino-massa ($\sum m_\nu$)

• De analyse combineert DESI data met een gecomprimeerde Planck CMB-likelihood.
• Resultaat: De posterior is consistent met de normale hiërarchie ($\sum m_\nu \approx 0.06$ eV).
• De omgekeerde hiërarchie ($\sum m_\nu \gtrsim 0.1$ eV) wordt uitgesloten op een betrouwbaarheidsniveau van 95%.
• Bovenlimiet: $\sum m_\nu < 0.10$ eV (95% CL). Dit is iets minder streng dan de volledige modellering van DESI ($<0.071$ eV), wat wordt toegeschreven aan het niet combineren van full-shape en BAO-data voor dezelfde tracers vanwege ontbrekende covariantie-termen.

5. Significatie en Conclusie

Deze studie demonstreert dat de ShapeFit-methode, wanneer gecombineerd met het bispectrum, een krachtige en model-onafhankelijke tool is voor kosmologische inferentie.

• Efficiëntie: Het biedt bijna dezelfde statistische kracht als complexe "Full Modelling" analyses, maar is modularer en makkelijker aan te passen voor verschillende kosmologische modellen.
• Validatie: De resultaten bevestigen de robuustheid van de DESI DR1-data en de consistentie tussen verschillende analysemethoden.
• Toekomst: Hoewel er momenteel geen sterke aanwijzingen zijn voor nieuwe fysica buiten $\Lambda$CDM, zal de toenemende volume van DESI (met een factor ~10 in de toekomst) de statistische fouten verkleinen. Dit vereist verdere verfijning van de methodologie (bijv. betere modellering van vensterfuncties voor het bispectrum) om systematische fouten onder controle te houden.

Kortom, dit werk legt een fundament voor toekomstige analyses van hoge-precisie kosmologische data door het succesvol te integreren van hogere-orde statistieken in een gecomprimeerde, flexibele analyseframework.