Probing gravity with non-linear clustering in redshift space

Probleemstelling
De waargenomen versnelde expansie van het universum wordt doorgaans toegeschreven aan donkere energie of een kosmologische constante ($\Lambda$), maar alternatieve theorieën stellen modificaties van de Algemene Relativiteitstheorie (GR) voor op kosmologische schalen voor. Een primaire uitdaging bij het testen van deze theorieën is het doorbreken van de degeneratie tussen modellen van gemodificeerde zwaartekracht (MG) en modellen van donkere energie, met name omdat afschermingsmechanismen in MG-theorieën (zoals het chameleon-mechanisme in $f(R)$-zwaartekracht) afwijkingen in gebieden met een hoge dichtheid, zoals ons Zonnestelsel, onderdrukken. De $E_G$-schatter werd voorgesteld als een schaalonafhankelijke probe voor zwaartekracht op grote schalen, gedefinieerd als de verhouding tussen het krommingspotentiaal (gerelateerd aan zwakke lenswerking) en het snelheidsveld (gerelateerd aan redshift-ruimte distorties, RSD). Huidige analyses gaan echter vaak uit van een schaalonafhankelijke groeisnelheid en vertrouwen op lineaire theorie, wat onvoldoende kan zijn voor het detecteren van subtiele afwijkingen in MG-modellen waar de groeisnelheid schaalafhankelijk wordt. Bovoven houdt ook de vraag open of de bekwaamheid van $E_G$ om levensvatbare MG-modellen (zoals het Hu & Sawicki $f(R)$-model met $|f_{R0}| = 10^{-5}$, aangeduid als F5) te onderscheiden van GR in de aanwezigheid van niet-lineaire clustering.

Methodologie
De auteurs voeren een end-to-end kosmologische analyse uit met hoogwaardige gesimuleerde galaxy-mocks, afgeleid van twee identieke twin N-body simulaties: één gebaseerd op GR en de andere op het F5 Hu & Sawicki $f(R)$-model. Beide simulaties delen identieke begincondities en kosmologische parameters ($\Omega_m=0.3089$, $\Omega_\Lambda=0.6911$, etc.), maar verschillen in de wet van de zwaartekracht.

1. Datageneratie: Galaxy-mocks worden gegenereerd met een combinatie van Halo Occupation Distribution (HOD) en Sub-Halo Abundance Matching (SHAM) toegepast op de simulaties. De studie maakt gebruik van een full-sky lightcone output met spectroscopische redshifts, waarbij zes verschillende galaxy-samples (full, central, red, blue, bright en faint) over drie redshift-bins ($z=0.35, 0.55, 0.8$) worden geselecteerd.
2. Constructie van de Schatter: De auteurs hanteren de verfijnde $E_G$-schatter beschreven in Ref. [106], die hoekvermogensspectra ($C_\ell$) in configuratieruimte gebruikt. De schatter is geformuleerd als:
   $$\hat{E}_G^\ell(z_{eff}) \approx \Gamma(z_{eff}) \frac{C_{\kappa g^*}^\ell}{\beta C_{gg}^\ell}$$
   waarbij $C_{gg}^\ell$ het auto-vermogensspectrum van de galaxies is, $C_{\kappa g^*}^\ell$ het gewogen convergentie-galaxy cross-vermogensspectrum is, en $\beta = f/b_1$ wordt afgeleid van RSD-analyses.
3. Schatting van Groeisnelheid en Bias:
   • Lineaire Bias ($b_1$): Wordt onafhankelijk geschat uit de hoek-autocorrelatie van de galaxies met behulp van de code pyCCL, wat een Gaussische prior biedt voor daaropvolgende RSD-fits.
   • Groeisnelheid ($f$): Wordt geschat door te fitten op de monopole en quadrupole van de 2-punts correlatiefunctie in configuratieruimte. De auteurs gebruiken de COMET-VDG emulator, die de Velocity Difference Generating function model implementeert om niet-lineaire RSD-effecten (inclus_n de Finger-of-God effect) nauwkeuriger te vangen dan standaard Effective Field Theory (EFT) benaderingen op kleine schalen.
4. MG Boosting: Omdat COMET niet van nature MG-kosmologieën ondersteunt, passen de auteurs de power spectrum boost toe van de e-mantis emulator (specifiek voor $f(R)$) op het matter power spectrum voordat de theoretische multipolen worden gegenereerd.
5. Statistische Analyse: Een Markov Chain Monte Carlo (MCMC) benadering met behulp van MultiNest wordt gebruikt om de modelparameters te fitten. De covariantie-matrix wordt geschat via 100 Jackknife resampling regio's. De uiteindelijke $E_G$-distributie wordt berekend met een ratio-distributie van de observables om rekening te houden met niet-Gaussianiteit.

Belangrijkste Bijdragen

• End-to-End Pipeline: De studie presenteert een uitgebreide pipeline voor het berekenen van $E_G$ met realistische galaxy-mocks in gemodificeerde zwaartekracht, waarbij gebruik wordt gemaakt van state-of-the-art niet-lineaire modellering (COMET-VDG) en schaalafhankelijke groeisnelheid-schatting.
• Schaalafhankelijke Analyse: In tegen tegenstelling tot eerdere werken die vaak schaalonafhankelijkheid veronderstellen, test dit artikel de $E_G$-schatter expliciet over verschillende schaalbereiken (klein versus groot) om de specifieke signaturen van $f(R)$-zwaartekracht te onderzoeken.
• Null Test Voorstel: De auteurs stellen een "null test" van zwaartekracht voor op basis van de 2-punts clustering statistieken in de redshift-ruimte. Deze test vergelijkt direct de verhouding van de datavectoren (F5 versus GR), waarbij de complexiteit van de $E_G$-schatter wordt omzeild om de intrinsieke discriminatieve kracht van kleine-schaal clustering te beoordelen.
• Robuustheidstests: De analyse bevat uitgebreide robuustheidstests met betrekking tot de selectie van de galaxy-sample (kleur, magnitude, centraliteit), de keuze van de covariantie-matrix (Jackknife versus theoretisch versus diagonaal) en de aannames over de priors van niet-lineaire bias parameters.

Resultaten

• Discriminatieve Kracht van $E_G$: Zelfs voor een ideale, all-sky, van systematische fouten vrije survey, kan de $E_G$-schatter niet duidelijk onderscheid maken tussen de F5- en GR-modellen. De foutmarges zijn te groot en de voorspelde verschillen tussen de modellen (ongeveer 5%) zijn vergelijkbaar met de statistische onzekerheden.
• Schaalafhankelijkheid: De geschatte $E_G$-waarden vertonen een trend die tegengesteld is aan de theoretische voorspellingen wanneer ze worden gesplitst naar schaal (overvoorspeld op kleine schalen, ondervoorspeld op grote schalen), wat waarschijnlijk te wijten is aan modelleringsbeperkingen en projectie-effecten in het niet-lineaire regime.
• Bias in Groeisnelheid: De groeisnelheid $f$ wordt met enige bias gerecupereerd op kleine schalen, met name voor het F5-model, waar deze te laag wordt ingeschat. Deze bias wordt echter gedeeltelijk gecompenseerd door de clustering ratio term in de $E_G$ berekening.
• Succes van de Null Test: In contrast met de $E_G$-schatter, vertoont de voorgestelde null test gebaseerd op de ratio van de 2-punts correlatiefunctie multipolen (monopole en quadrupole) aanzienlijke discriminatieve kracht. De datavectoren van F5 en GR verschillen op een significantieniveau van $\sim 3\sigma$ tot $4\sigma$ over de redshift-bins, wat aangeeft dat de informatie aanwezig is in de data, maar wordt vertroebeld door de specifieke constructie van de $E_G$-schatter en huidige modelleringsbeperkingen.
• Modelleringsbeperkingen: De hoge $\chi^2_\nu$ waarden die worden waargenomen in de fits, vooral op kleine schalen, suggereren dat huidige perturbatieve modellen (VDG/EFT) moeite hebben met het accuraat beschrijven van clustering op schalen onder $\sim 10-20$ Mpc/$h$, waar de sterkste afwijkingen van GR worden verwacht.

Betekenis en Claims
Het artikel claimt bescheiden dat hoewel de $E_G$-schatter een theoretisch robuuste, bias-onafhankelijke test van zwaartekracht is, de praktische bruikbaarheid ervan voor het onderscheiden tussen levensvatbare MG-modellen (zoals F5) en GR beperkt wordt door de huidige capaciteiten van niet-lineaire modellering en statistische ruis, zelfs in geïdealiseerde scenario's. De auteurs benadrukken dat het falen om modellen te onderscheiden via $E_G$ niet noodzakelijkerwijs een falen van de data is, maar eerder een beperking van de constructie van de schatter, die de schaalafhankelijke signalen waar MG-effecten het sterkst zijn, middelt.

De studie concludeert dat toekomstige vooruitgang vereist:

1. Verbeterde Modellering: Het ontwikkelen van MG-consistente RSD-templates en emulators (bijv. BACCO) die clustering op kleinere schalen ($<10$ Mpc/$h$) accuraat kunnen modelleren.
2. Alternatieve Probes: Het benutten van de volledige informatie van de redshift-ruimte clustering (zoals aangetoond door de null test) of het combineren van meerdere observables (bijv. 3x2pt analyses) om degeneraties tussen nuisance parameters (zoals galaxy bias) en kosmologische parameters te doorbreken.
3. Realistische Voorspellingen: De gepresenteerde resultaten zijn bovengrenzen onder ideale omstandigheden; echte surveys zullen te maken krijgen met extra systematische effecten (baryonische fysica, intrinsic alignments, fotometrische fouten) die de discriminatieve kracht van $E_G$ verder zullen verslechteren.

Uiteindelijk dient dit werk als een benchmark voor de beperkingen van huidige $E_G$-methodologieën en motiveert het de ontwikkeling van meer geavanceerde, niet-lineaire modelleringskaders om het potentieel van de volgende generatie galaxy-surveys (zoals Euclid, DESI en LSST) volledig te benutten in het testen van zwaartekracht.