Beyond Linear Bias Expansions for AbacusSummit Halos at z = 8

Probleemstelling

Op hoge roodverschuivingen (zoals $z=8$) zijn massieve halos, die heldere sterrenstelsels kunnen herbergen, zeldzaam en vormen ze zich uitsluitend in de meest overdense gebieden van het universum. Hierdoor vertonen ze een sterke "bias" (voorkeur) ten opzichte van de onderliggende materie. Traditionele lineaire bias-modellen ($\delta_h = b_1 \delta$) zijn ontoereikend omdat de fluctuaties in de dichtheid van deze tracers groot zijn (orde 1 of groter), wat leidt tot sterk niet-Gaussische clustering.

De uitdaging ligt in het bepalen van welke termen in een bias-uitbreiding (bias expansion) noodzakelijk zijn om deze systemen correct te modelleren. Bestaande theorieën veronderstellen vaak dat termen met hogere orde in de materie-overdichtheid ($\delta^n$) verwaarloosbaar zijn vergeleken met lagere orde termen, tenzij de bias-coëfficiënten extreem groot zijn. Bij hoge $z$ wordt echter vermoed dat de verhouding tussen de bias-coëfficiënten ($b_2/b_1$, $b_3/b_2$, etc.) zo groot is dat hogere-orde termen in de materie-overdichtheid dominant worden, zelfs boven termen die gerelateerd zijn aan de vorm van het dichtheidsveld (zoals getijdenkrachten).

Methodologie

De auteurs analyseren halos uit de AbacusSummit suite van N-body-simulaties, specifiek gericht op roodverschuiving $z=8$ en ter verificatie bij $z=5$.

1. Data: Ze gebruiken 25 simulatieboxen met een zijde van $2 h^{-1}$ Gpc en een resolutie van $6912^3$ deeltjes. Ze selecteren halos met een massa $M \geq 1 \times 10^{11} h^{-1} M_\odot$.
2. Observabelen: De analyse is gebaseerd op het vermogensspectrum ($P(k)$) en het bispectrum ($B(k_1, k_2, k_3)$) van de halo-overdichtheid.
3. Bias Uitbreiding: Ze gebruiken een Euleriaanse bias-uitbreiding die lokale termen in de materie-overdichtheid ($\delta$) en getijdentermen ($K_{ij}$) omvat. De uitbreiding omvat:
   • Lineaire term: $b_1 \delta$
   • Kwantumtermen: $b_2 \delta^2$, $b_3 \delta^3$
   • Getijdenterm: $b_{K^2} K^2$
   • Niet-lineaire materie: $\delta^{(2)}$ (de tweede orde perturbatie van de materie).
4. Renormalisatie: Om afhankelijkheid van de schaalafsnijding ($\Lambda$) te elimineren, worden de bias-operatoren renormaliseerd (bijv. $[\delta^2] = \delta^2 - \sigma^2(\Lambda)$).
5. Fitting: Ze passen de bias-coëfficiënten aan op de simulatiedata door een $\chi^2$-fitting uit te voeren op het vermogensspectrum en het bispectrum. Ze nemen rekening met shot noise en stochasticiteit door extra parameters ($b_0, b'_0, b''_0$) toe te voegen.
6. Vergelijking: Ze testen verschillende scenario's, waaronder het variëren van de massadrempel en het vergelijken van $z=8$ met $z=5$ (waarbij de massadrempel wordt aangepast om dezelfde variantie in het haloveld te behouden).

Belangrijkste Bijdragen

• Validatie van Hogere Orde Termen: Het artikel biedt een kwantitatieve onderbouwing dat bij hoge roodverschuivingen termen zoals $\delta^2$ en $\delta^3$ cruciaal zijn en vaak belangrijker zijn dan getijdentermen, ondanks dat ze in de perturbatietheorie een hogere orde hebben.
• Niet-Gaussische Kwantificering: Het toont aan dat de niet-Gaussische aard van de clustering bij $z=8$ voornamelijk wordt gedreven door de lokale dichtheid (magnitude van overdichtheid) en niet door de vorm (getijdenkrachten) van het dichtheidsveld.
• Efficiënt Model: Het stelt een vereenvoudigd model voor dat alleen lineaire en kwadratische bias-termen (plus een zwakke detectie van $\delta^3$) nodig heeft, mits het onderliggende materieveld als Gaussisch wordt behandeld voor de basis, maar met expliciete niet-lineaire correcties in de bias-uitbreiding.

Resultaten

De analyse levert de volgende specifieke resultaten op voor halos bij $z=8$:

1. Dominantie van Dichtheidstermen:

   • De term $b_2$ (kwadratisch in $\delta$) wordt zeer sterk gedetecteerd (17$\sigma$).
   • De term $b_3$ (kubisch in $\delta$) wordt zwak gedetecteerd op het 1.3$\sigma$ niveau.
   • De getijdenbias ($b_{K^2}$) is niet significant gedetecteerd; de afwijking van nul is slechts $\sim 0.1\sigma$.
   • De niet-lineariteit van het materieveld zelf ($\delta^{(2)}$) wordt ook gedetecteerd op het 1.3$\sigma$ niveau.

2. Vergelijking Vermogensspectrum vs. Bispectrum:

   • De beperkingen uit het vermogensspectrum en het bispectrum zijn consistent met elkaar.
   • Het vermogensspectrum geeft strakkere beperkingen op $b_1$ en $b_3$, terwijl het bispectrum aanzienlijk strakkere beperkingen geeft op de getijdenbias ($b_{K^2}$).

3. Evolutie met Roodverschuiving ($z=8$ vs $z=5$):

   • Bij $z=5$ (met een aangepaste massadrempel voor gelijke variantie) zijn de bias-coëfficiënten kleiner, wat overeenkomt met de verwachte lineaire evolutie van het materieveld.
   • Bij $z=5$ wordt de getijdenbias wel gedetecteerd op het 1$\sigma$ niveau, wat aangeeft dat getijdenkrachten relatief belangrijker worden naarmate de roodverschuiving daalt, maar ze blijven ondergeschikt aan de dichtheidstermen.
   • De bias-coëfficiënten bij $z=5$ kunnen tot op 10% nauwkeurig worden voorspeld door de $z=8$-waarden lineair te evoluteren, wat suggereert dat het dezelfde populatie halos is die evolueert.

4. Modelprestaties:

   • Het model past de data uitstekend zonder systematische fouten (zoals zichtbaar in de residuals van het bispectrum).
   • Termen die lokaal zijn in de materie-overdichtheid zijn essentieel; het weglaten ervan leidt tot slechte fits.

Significantie en Conclusie

De studie concludeert dat de klassieke perturbatieve uitbreiding, waarbij $\delta$ als de enige verstoorparameter wordt gezien, niet ideaal is voor hoge roodverschuivingen. Omdat de verhoudingen tussen de bias-coëfficiënten ($b_n/b_1$) groot zijn, moeten termen met hogere orde in $\delta$ (zoals $\delta^2$ en $\delta^3$) worden meegenomen, zelfs als ze in de theorie "hoger orde" zijn.

• Praktische Implicatie: Voor toekomstige surveys van hoge-roodverschuiving-galaxieën (bijv. met JWST) is het voldoende om bias-uitbreidingen te gebruiken die bestaan uit alleen $b_1$ en $b_2$ (en eventueel $b_3$), gebaseerd op het lineaire materieveld. Getijdenbijdragen kunnen voor een efficiënte beschrijving worden genegeerd of als verwaarloosbaar worden beschouwd.
• Fysisch Inzicht: De vorming van vroege halos wordt primair bepaald door de magnitude van de lokale materie-overdichtheid (hoe dicht het gebied is) en niet door de vorm (getijdenkrachten) van het dichtheidsveld. Dit bevestigt het fysieke beeld dat halos bij hoge $z$ alleen in de uiterst extreme pieken van het dichtheidsveld ontstaan.

Dit werk biedt een robuust kader voor het interpreteren van clustering-data uit de kosmische dageraad en onderstreept de noodzaak van niet-lineaire bias-modellen voor nauwkeurige kosmologische analyses in de toekomst.