The Impact of Galaxy Formation on Galaxy Biasing, and Implications for Primordial non-Gaussianity Constraints

Technische Samenvatting: De impact van de vorming van sterrenstelsels op de bias van sterrenstelsels, en implicaties voor beperkingen op primordiale niet-Gaussianiteit

Probleemstelling
De parameter $f_{\text{NL}}$ kwantificeert de lokale primordiale niet-Gaussianiteit (PNG) in de initiële energiefluctuaties van het Universum. Een detectie van een niet-nul waarde zou onderscheid maken tussen single-field en multi-field inflatiemodellen. Hoewel de volgende generatie sterrenstelsel-surveys streeft naar het beperken van $f_{\text{NL}}$ met statistische fouten van $\sigma_{f_{\text{NL}}} \sim 1$, domineren systematische fouten momenteel deze beperkingen. Specifiek is de meting van $f_{\text{NL}}$ via de clustering van sterrenstelsels degeneratief met de bias-parameters van sterrenstelsels: $b_1$ (respons op de massadichtheid) en $b_\phi$ (respons op primordiale potentiaalfluctuaties).

Huidige beperkingen vertrouwen vaak op de "universaliteitsrelatie", die een vaste schaling aanneemt tussen $b_\phi$ en $b_1$ (bijv. $b_\phi = 2\delta_c(b_1 - 1)$), afgeleid van donkere materie halo's. Echter, recent werk (bijv. Barreira et al. 2020) geeft aan dat deze relatie niet opgaat voor sterrenstelsels, aangezien $b_\phi$ afhankelijk is van de selectie van sterrenstelsels en de specifieke fysica van de vorming van sterrenstelsels. De complexiteit van het modelleren van de vorming van sterrenstelsels introduceert een aanzienlijke onzekerheid in $b_\phi$, wat de precisie van toekomstige $f_{\text{NL}}$-beperkingen bedreigt.

Methodologie
Om de onzekerheid in $b_\phi$ aan te pakken die voortvloeit uit de modellering van de vorming van sterrenstelsels, maken de auteurs gebruik van het Santa Cruz Semi-Analytische Model (SC-SAM), ingebed in aparte-universum N-body simulaties.

• Simulaties: Drie AREPO N-body simulaties werden uitgevoerd met een doosgrootte van $L = 205 h^{-1}$ cMpc en $N = 1280^3$ deeltjes. Deze simulaties delen identieke Gaussische initiële condities, maar verschillen in de amplitude van het primordiale vermogensspectrum ($\sigma_8$), waarbij de techniek van het aparte universum wordt gebruikt om de respons van de vorming van sterrenstelsels op langgolvige perturbaties efficiënt te meten.
• Modelvariaties: De auteurs voerden 56 unieke instanties van de SC-SAM uit. Deze variaties omvatten het modificeren van drie feedbackparameters:
  1. $\epsilon_{\text{SN}}$ (amplitude van de door sterren gedreven massaverplaatsing, gevarieerd via factor $A_{\text{SN1}}$).
  2. $\alpha_{\text{rh}}$ (exponent van de massaverplaatsingsrelatie, gevarieerd via factor $A_{\text{SN2}}$).
  3. $\kappa_{\text{radio}}$ (amplitude van de AGN radio-mode feedback, gevarieerd via factor $A_{\text{AGN}}$).
• Metingen:
  • $b_\phi$ wordt berekend via eindige differentiatie van de aantalconcentraties van sterrenstelsels tussen hoog-$\sigma_8$, fiducieel en laag-$\sigma_8$ simulaties.
  • $b_1$ wordt afgeleid van het kruis-vermogensspectrum van sterrenstelsels en materie.
  • Sterrenstelsels werden geselecteerd op basis van Stellaire Massa ($M_*$), Stervormingssnelheid (SFR) en specifieke SFR (sSFR) over roodverschuivingen $z=0, 1, 2$.
  • Een "realisme"-metriek ($\chi^2$) werd aan elk model toegewezen op basis van de overeenstemming met geobserveerde $z=0$ stellaire massafuncties (Bernardi et al. 2018).
  • De studie beoordeelde ook de impact van galaxy assembly bias door sterrenstelsels te verschuiven naar halo's met een vergelijkbare massa maar verschillende vormingsgeschiedenissen.

Belangrijkste Bijdragen

1. Systematische Exploratie van Feedback: In tegen tegenovercede studies die vaak een enkele hydrodynamische simulatie of een vaste SAM gebruikten, varieert dit werk systematisch de stellaire en AGN feedbackparameters om de resulterende landschappen van $b_\phi$ en $b_1$ in kaart te brengen.
2. Toepassing van de Apart-Universum SAM: De auteurs passen de aparte-universum techniek toe op een semi-analytisch model, wat een efficiënte meting van $b_\phi$ mogelijk maakt over een brede parameterruimte die computationeel onhaalbaar zou zijn met volledige hydrodynamische simulaties.
3. Analyse van Selectieafhankelijkheid: De paper vergelijkt expliciet hoe verschillende selectiecriteria voor sterrenstelsels ($M_*$, SFR, sSFR) reageren op variaties in de fysica van de vorming van sterrenstelsels, waarbij wordt geïdentificeerd welke selecties het meest robuust zijn tegen modelleringsonzekerheden.
4. Cross-Model Validatie: De resultaten worden vergeleken met eerdere bevindingen uit de IllustrisTNG hydrodynamische simulatie, de Munich SAM en de Galacticus SAM om te bepalen of de SC-SAM parameterruimte het gedrag van andere modellen van sterrenstelselvorming omvat.

Resultaten

• Afwijking van Universaliteit: De SC-SAM modellen bevestigen dat de universaliteitsrelatie ($b_\phi = 2\delta_c(b_1 - 1)$) niet accuraat is voor sterrenstelsels. De $b_\phi(b_1)$ relatie varieert aanzienlijk afhankelijk van de selectie van sterrenstelsels en de specifieke feedbackparameters.
• Stellaire Massa Selecties: $b_\phi$ neemt over het algemeen toe met de stellaire massa. De $b_\phi(b_1)$ relatie is echter steiler dan die van IllustrisTNG (Barreira et al. 2020), met name voor selecties met een lage massa. De relatie is gevoelig voor de stellaire feedbackparameters ($A_{\text{SN1}}, A_{\text{SN2}}$).
• SFR Selecties: De $b_\phi$ waarden voor op SFR geselecteerde sterrenstelsels zijn chaotisch en vaak dicht bij nul. De relatie is zeer gevoelig voor feedbackparameters omdat verschillende parameters de SFR–halo massa relatie op complexe wijze verschuiven.
• sSFR Selecties (Kernbevinding): Selecties gebaseerd op de specifieke stervormingssnelheid (sSFR) zijn opmerkelijk robuust tegen variaties in de modellering van de vorming van sterrenstelsels. Over alle 56 modellen volgen de $b_\phi$ vs. sSFR en de $b_\phi(b_1)$ relaties een strak, consistent pad. Dit wordt toegeschreven aan de zelfregulerende aard van de SC-SAM, die een consistente sSFR–halo massa relatie handhaaft, ongeacht de variaties in feedbackparameters.
• Assembly Bias: De studie vindt dat galaxy assembly bias de $b_1$ metingen beïnvloedt (het vermindert clustering voor stellaire massa selecties, en verhoogt het voor SFR/sSFR selecties wanneer de vormingsgeschiedenis wordt gehusseld). Echter, deze assembly bias signalen hebben een minimale impact op de gemeten $b_\phi(b_1)$ relaties.
• Vergelijking met Andere Modellen: De SC-SAM variaties omvatten de $b_\phi$ gedragingen die gemeten zijn in de Munich SAM en Galacticus. Hoewel er spanningen zijn met de IllustrisTNG resultaten (waarschijnlijk door baryonische effecten die afwezig zijn in SAM's of verschillen in de constructie van merger trees), dekken de SC-SAM variaties een breed bereik van $b_\phi$ waarden dat consistent is met andere modellen.

Betekenis en Claims
Het artikel claimt dat onzekerheden in de modellering van de vorming van sterrenstelsels een kritieke systematische fout vormen voor toekomstige $f_{\text{NL}}$ beperkingen. Door aan te tonen dat $b_\phi$ niet universeel is en sterk afhangt van selectie en modellering, betogen de auteurs dat het aannemen van een vaste universaliteitsrelatie het risico loopt op bevooroordeelde (biased) beperkingen.

De primaire betekenis van dit werk ligt in het identificeren van sSFR als een robuuste selectie van sterrenstelsels voor $f_{\text{NL}}$ studies. Omdat sSFR selecties consistente $b_\phi$ waarden opleveren over een brede reeks fysica van de vorming van sterrenstelsels, bieden ze een pad om systematische fouten in toekomstige surveys te verminderen. De auteurs concluderen dat hoewel geen enkel model de het universum perfect kan beschrijven, het selecteren van sterrenstelsels op basis van sSFR de gevoeligheid van $f_{\text{NL}}$ beperkingen voor de specifieke details van de vorming van sterrenstelsels minimaliseert. Bovendien biedt het werk een raamwerk voor het bouwen van informatieve priors voor $b_\phi$ op basis van observationele data, bewegend naar fysisch gemotiveerde beperkingen op de primordiale niet-Gaussianiteit.