Galaxy formation in modified gravity -- II. galaxy halo connection and assembly bias

Technische Samenvatting: Galactievorming in gemodificeerde zwaartekracht – II. de verbinding tussen sterrenstelsel en halo en assembly bias

Probleemstelling
Moderne kosmologische surveys (bijv. DESI, Euclid) beogen kosmologische parameters te beperken en theorieën over zwaartekracht te testen door de grootschalige structuur (LSS) van het Universum in kaart te brengen. De interpretatie van deze gegevens wordt echter beperkt door onzekerheden in de "galaxy-halo connection"—de mapping tussen geobserveerde sterrenstelsels en de onderliggende donkere materie-halo's. Standaardanalyses vertrouwen vaak op het Halo Occupation Distribution (HOD) model, dat ervan uitgaat dat de bezetting van sterrenstelsels uitsluitend afhangt van de halo-massa. Deze aanname negeert "assembly bias" (AB), waarbij de vorming van sterrenstelsels ook afhangt van secundaire halo-eigenschappen zoals de vormingsgeschiedenis, omgeving en concentratie. Hoewel AB bekend staat om het beïnvloeden van clustering in het standaard $\Lambda$CDM-model, blijft het gedrag ervan in Modified Gravity (MG) theorieën, met name $f(R)$-zwaartekracht, slecht begrepen. Het complexe samenspel tussen MG-screeningmechanismen (bijv. chameleon screening) en de fysica van sterrenstelselvorming zou significante biases kunnen introduceren in de schatting van kosmologische parameters als dit niet correct wordt meegewogen.

Methodologie
De auteurs maken gebruik van geavanceerde hydrodynamische simulaties uitgevoerd met de MG-AREPO code, die de volledige baryonische fysica van het IllustrisTNG-model incorporeert maar aangepast voor $f(R)$-zwaartekracht. De simulatiesuite bevat een reeks Hu-Sawicki $f(R)$-modellen met variërende sterkte van de vijfde kracht, aangeduid door de achtergrond-scalarveld-amplitude vandaag, $|f_{R0}| = 10^{-x}$ (waarbij $x \in \{4.0, 4.5, 5.0, 5.5, 6.0\}$), naast een standaard $\Lambda$CDM-controlegroep ($f_{R0}=0$).

Uit deze simulaties genereren de auteurs mock-catalogi voor twee populaties sterrenstelsels die relevant zijn voor Stage-IV surveys:

1. Luminous Red Galaxies (LRGs): Geselecteerd op basis van stellaire massa-rangschikking.
2. Emission Line Galaxies (ELGs): Geselecteerd op basis van de specifieke stervormingssnelheid (sSFR) en een stellaire massagrens ($M_\star > 2.5 \times 10^{10} M_\odot$) om resolutie-effecten te mitigeren.

De studie verloopt in drie fasen:

1. HOD Fitting: De auteurs passen standaard HOD-modellen toe (Zheng et al. 2005 voor LRGs; Alam et al. 2020 voor ELGs) op de simulatiedata om te bepalen hoe HOD-parameters evolueren met roodverschuiving en de sterkte van de zwaartekracht.
2. Secundaire Afhankelijkheidsanalyse: De sterrenstelsel-catalogi worden opgesplitst in kwartielen op basis van secundaire halo-eigenschappen: omgevingsdichtheid (gedefinieerd door overdensiteit in een sferische annulus) en halo-concentratie. HOD-parameters worden onafhankelijk gefit voor elk kwartiel om de afhankelijkheid van de bezetting van sterrenstelsels van deze eigenschappen te kwantificeren.
3. Kwantificering van Assembly Bias: De auteurs genereren mock-sterrenstelsel-catalogi met drie HOD-voorschriften: (i) de basis massale HOD, (ii) een omgeving-afhankelijke HOD, en (iii) een concentratie-afhankelijke HOD. Ze vergelijken de grootschalige correlatiefunctie ($\xi_{LS}$) van deze mocks met de "ground truth" die direct uit de hydrodynamische simulaties is gemeten om de resterende assembly bias te kwantificeren.

Belangrijkste Bijdragen en Resultaten

• Evolutionaire Sequentie van MG-effecten: De studie identificeert een duidelijke evolutionaire sequentie in hoe MG de galaxy-halo verbinding beïnvloedt. Naarmate de vijfde kracht de halo's in de loop van de tijd "ontscreened", bewegen modellen zich door regimes van initiële ontscreening, inhaalbeweging (catch-up) en quenching. Deze sequentie verklaart het niet-monotone gedrag van HOD-parameters (bijv. $M_c$, $p_{max}$, $M_1$) over verschillende zwaartekrachtmodellen en roodverschuivingen heen. Bijvoorbeeld, in zwakkere MG-modellen vindt de transitie van gescreende naar ongescreende toestanden later plaats, wat leidt tot andere HOD-gedragingen vergeleken met sterkere modellen bij dezelfde roodverschuiving.
• Verzadigingseffect: De auteurs bevestigen een "verzadigingseffect" waarbij de impact van MG op de vorming en clustering van sterrenstelsels een limiet bereikt. Sterkere MG-modellen (bijv. F4.0, F4.5) halen uiteindelijk een punt waarop ze wat betreft HOD-parameters "bij komen" met elkaar, wat suggereert dat het universum voor een vastgesteld selectiecriterium voor sterrenstelsels niet oneindig meer sterrenstelsels kan produceren; MG versnelt primair de vormingstijdlijn in plaats van een oneindig overschot aan sterrenstelsels te creëren.
• Secundaire Afhankelijkheden:
  • Omgeving: Er is een sterke afhankelijkheid van HOD-parameters van de omgevingsdichtheid in zowel $\Lambda$CDM als MG-modellen. In MG wordt dit verergerd door het chameleon-screeningmechanisme, dat ervoor zorgt dat ontscreening eerst optreedt in omgevingen met een lage dichtheid.
  • Concentratie: Hoewel HOD-parameters een significante afhankelijkheid vertonen van de halo-concentratie, vertaalt deze afhankelijkheid zich niet naar een significante vermindering van de assembly bias wanneer ze worden gebruikt om mocks te genereren.
• Magnitude van Assembly Bias:
  • Het negeren van assembly bias (het gebruik van een basis massale HOD) onderschat de clusteringsterkte van sterrenstelsels met 10–20% in $\Lambda$CDM bij $z \lesssim 1$.
  • In MG-modellen is het AB-effect niet-monotoon en kan het aanzienlijk verschillen van $\Lambda$CDM. Bijvoorbeeld, in het F5.5-model is het AB-effect sterker dan in GR, terwijl het in F4.0 en F4.5 juist zwakker kan zijn.
  • Cruciaal is dat de richting en magnitude van het AB-effect in MG kunnen verschillen van $\Lambda$CDM, wat potentieel kan leiden tot een foutieve identificatie van het zwaartekrachtmodel als standaard analyse-pipelines worden gebruikt.
• Mitigatiestrategie: Het uitbreiden van het basis HOD-model met omgevingsdichtheid als secundaire variabele vermindert het assembly bias-effect tot 2–3% voor $z \lesssim 0.5$ voor alle modellen (inclusief $\Lambda$CDM en diverse $f(R)$-scenario's). In tegenstelling hiertoe leidt het uitbreiden van het model met halo-concentratie niet tot een significante vermindering van het AB-effect.

Betekenis en Claims
Het artikel claimt dat het basis HOD-model onvoldoende is voor precisie-kosmologische tests met betrekking tot gemodificeerde zwaartekracht, met name voor Emission Line Galaxies (ELGs). De auteurs stellen dat het complexe samenspel tussen MG-screening en de vorming van sterrenstelsels assembly bias-signaturen creëert die kwalitatief en kwantitatief verschillen van de standaard zwaartekracht.

De primaire betekenis ligt in het aantonen dat het opnemen van omgevingsdichtheid in het HOD-model een goed gemotiveerd en effectief startpunt biedt om systematische fouten door assembly bias te minimaliseren. Dit is bijzonder relevant voor aankomende surveys zoals DESI en Euclid, die schalen onderzoeken waar assembly bias sterk aanwezig is. De auteurs waarschuwen dat het niet rekening houden met deze effecten de kosmologische parameter-inferentie kan beïnvloeden of kan leiden tot onjuiste conclusies over de aard van de zwaartekracht. Ze concluderen dat, hoewel het omgeving-afhankelijke HOD de AB reduceert tot een residu-niveau ($\lesssim 3\%$ bij lage $z$), verder werk met grotere, hogere-resolutie simulaties nodig is om deze systematische fouten volledig te elimineren en meer geavanceerde HOD-uitbreidingen te testen.