Cosmic variance or galaxy bias? Disentangling finite-volume and galaxy formation effects in cosmological analysis

Technische Samenvatting: Ontrafelen van eindige-volume- en sterrenstelselformatie-effecten in kosmologische analyse

Probleemstelling
Huidige en toekomstige grootschalige galaxy redshift-surveys (bijv. DESI, Euclid) streven naar het bepalen van kosmologische parameters met ongekende precisie. Het bereiken van onbevooroordeelde resultaten vereist echter een strikte controle over theoretische systematische fouten. Een primaire uitdaging is het onderscheiden van twee sterk gekoppelde onzekerheidsbronnen: kosmische variantie, voortkomend uit de intrinsieke variantie van het initiële Gaussische willekeurige veld in een eindig volume, en galaxy bias, die de complexe niet-lineaire processen van sterrenstelselformatie beschrijft.

In standaard perturbatietheorieën worden bias-parameters geschat door observaties te fitten aan het ensemble-gemiddelde materie-vermogensspectrum. Deze benadering gaat impliciet ervan uit dat het survey-volume het ensemble-gemiddelde vertegenwoordigt, waardoor de specifieke kosmische variantie van de geobserveerde realisatie wordt genegeerd. Gevolglik worden de geschatte bias-parameters degeneraat, waarbij de werkelijke fysica van sterrenstelselformatie wordt samengevoegd met de statistische fluctuaties van het specifiek onderzochte kosmische volume. Deze degeneratie bemoeilijkt toepassingen zoals Baryon Acoustic Oscillation (BAO) reconstructie en de inferentie van kosmologische parameters, waar aannames over bias-waarden cruciaal zijn.

Methodologie
De auteurs ontwikkelen een nieuw formalisme gebaseerd op Euleriaanse Perturbatietheorie (EPT) om kosmische variantie van galaxy bias te ontrafelen.

1. Kosmische variantie herformuleren als bias: De auteurs stellen voor om het dichtheidsveld van een willekeurige eindige-volume-realisatie ($\delta'$) te behandelen als een "biased" tracer van het ensemble-gemiddelde dichtheidsveld ($\delta$). Door een nieuwe set bias-parameters ($b_{1,dm}, b_{2,dm}$) te introduceren, drukken zij het eindige volume-veld uit als:
   $$\delta'(x) = b_{0,dm} + b_{1,dm}\delta(x) + \frac{1}{2}b_{2,dm}\delta^2(x) + \dots$$
   In dit kader, naarmate het volume naar oneindig nadert, geldt $b_{1,dm} \to 1$ en $b_{2,dm} \to 0$. Deze formulering maakt het mogelijk om de "kosmische variantie-bias" te scheiden van de fysieke "galaxy bias" ($b_{1,g}, b_{2,g}$).

2. Toepassing op BAO-verschuiving: De auteurs passen dit formalisme toe op de niet-lineaire verschuiving van de BAO-schaal. Gebruikmakend van de analytische expressies afgeleid door [22], ontleden zij de totale BAO-verschuiving ($\Delta \alpha$) in twee afzonderlijke componenten:

   • Een component gedreven door de intrinsieke kosmische variantie van het donkere materie-veld (gecontroleerd door $b_{1,dm}, b_{2,dm}$).
   • Een component gedreven door galaxy biasing (gecontroleerd door $b_{1,g}, b_{2,g}$).
     Deze decompositie onthult dat kosmische variantie alleen al een niet-verwaarloosbare BAO-verschuiving kan induceren, met name in kleine volumes, wat degeneratief is met de verschuiving veroorzaakt door niet-lineaire gravitationele groei.

3. Afhandeling van de Integrale Constante: Het artikel behandelt de integrale constrictie (de bias die ontstaat door de gemiddelde dichtheid van een survey-volume op nul te dwingen). De auteurs demonstreren dat dit effect kan worden geabsorbeerd in de constante term van de bias-expansie ($b'_0$), waardoor het geen extra bias-parameters introduceert of de scheiding van kosmische variantie bemoeilijkt.

4. Numerieke Validatie: Om de analytische argumenten te valideren, maken de auteurs gebruik van de UNIT simulatie-suite. Zij analyseren vier verschillende sets van $N$-body simulaties ($2 \times 10^3$ deeltjes in $(1 h^{-1}\text{Gpc})^3$ volumes):

   • 100 standaard simulaties.
   • 100 fixed-amplitude simulaties (onderdrukking van amplitude-variantie).
   • 2 $\times$ 50 gepaarde simulaties (onderdrukking van fase-variantie).
   • 2 $\times$ 50 fixed-and-paired simulaties (onderdrukking van zowel amplitude als fase).
     Zij meten de BAO-piekpositie op een model-onafhankelijke wijze (door een Gaussische functie te fitten aan de residuele correlatiefunctie) om de reductie van de variantie te kwantificeren.

Belangrijkste Resultaten

• Analytische Ontrafeling: Het nieuwe formalisme slaagt erin de bijdragen van kosmische variantie en galaxy bias aan de niet-lineaire BAO-verschuiving te scheiden. De analyse bevestigt dat kosmische variantie alleen al een verschuiving kan induceren die vergelijkbaar is met die van galaxy bias, wat extra zorg vereist bij het gebruik van perturbatieve benaderingen op eindige volumes.
• Variantie-onderdrukking: Numerieke resultaten uit de UNIT-simulaties tonen aan dat het controleren van de begincondities de onzekerheid in de BAo-piekpositie aanzienlijk vermindert:
  • Fixing (controleren van amplitude) vermindert de onzekerheid met een factor $\sim 1.5 - 2.1$.
  • Pairing (controleren van fase) vermindert de onzekerheid met een factor $\sim 1.5 - 1.7$.
  • Gecombineerd (Fixed-and-Paired): De combinatie van beide technieken levert een variantie-onderdrukkingsfactor op van $\sim 3.4 - 5.4$ vergeleken met standaard ensembles.
• Convergentie: De BAO-piekpositie blijkt goed geconvergeerd te zijn op een sub-percentage niveau over alle simulatietypes heen, zonder dat er duidelijke systematische fouten worden geïntroduceerd door de fixing of pairing technieken ten opzichte van het standaard ensemble.

Betekenis en Voorgestelde Strategieën
Het artikel concludeert dat huidige methoden voor het schatten van bias-parameters onbedoeld de degeneratieve effecten van kosmische variantie bevatten. Om dit aan te pakken, stellen de auteurs twee praktische strategieën voor voor toekomstige kosmologische analyses:

1. Optimale Multi-Tracer Combinatie: In plaats van enkel de statistische precisie te maximaliseren, zouden toekomstige analyses complementaire tracers (bijv. sterrenstelsels en kosmische voids) kunnen combineren om specifiek de kosmische variantie te elimineren. Dit zou het herdefiniëren van pre-reconstructie monsters kunnen inhouden die complementaire dichtheidsregimes volgen, of het optimaliseren van gewichten in de gecombineerde correlatiefunctie om het effect van gepaarde realisaties na te bootsen.
2. Variantie-onderdrukkende Reconstructie: De auteurs suggereren het toepassen van variantie-onderdrukkende technieken (fixing en pairing) op zowel de geobserveerde data (via reconstructiemethoden zoals backward inversion of Bayesian forward reconstruction) als op de mock-catalogi die gebruikt worden voor covariantie-estimatie. Dit zou de analyse van een "variantie-onderdrukt" versie van het geobserveerde Universum en zijn gepaarde anti-Universum mogelijk maken. Zij merken op dat dit een zorgvuldige implementatie vereist om te voorkomen dat de BAO-verschuiving wordt beïnvloed door de gehanteerde kosmologie tijdens het reconstructieproces.

De auteurs benadrukken dat hoewel deze strategieën veelbelovend zijn, ze verdere tests vereisen op gedetailleerde mock-catalogi (bijv. DESI-achtig) en integratie met state-of-the-art analystechnieken zoals BAO-reconstructie en Effective Field Theory.