🔭 astrophysics

Galaxy assembly bias in cosmological hydrodynamical simulations: a comparison between SIMBA and IllustrisTNG

Deze studie vergelijkt de galaxy assembly bias in de SIMBA- en IllustrisTNG-simulaties en concludeert dat de bias voornamelijk wordt veroorzaakt door de kosmische dichtheid, met een opvallend divergerende evolutie tussen de twee modellen na z=2 waarbij SIMBA een afname en TNG een toename toont.

Oorspronkelijke auteurs: Hong-Gang Yang, Marcos Pellejero, John Peacock, Romeel Davé

Gepubliceerd 2026-03-11

📖 5 min leestijd🧠 Diepgaand

CC BY 4.0

Oorspronkelijke auteurs: Hong-Gang Yang, Marcos Pellejero, John Peacock, Romeel Davé

Oorspronkelijk artikel gelicentieerd onder CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Dit is een AI-gegenereerde uitleg van het onderstaande artikel. Het is niet geschreven of goedgekeurd door de auteurs. Raadpleeg het oorspronkelijke artikel voor technische nauwkeurigheid. Lees de volledige disclaimer

De Verborgen Invloed van de Buurman: Hoe Galaxieën Groeperen in het Heelal

Stel je het heelal voor als een gigantisch, donker bos. In dit bos zijn de bomen niet willekeurig verspreid; ze groeien in groepjes. De "bomen" zijn hier de sterrenstelsels (galaxieën), en de "bosjes" zijn onzichtbare wolken van donkere materie die we halo's noemen.

Vroeger dachten wetenschappers dat het alleen maar om de grootte van het bosje ging: een groot bosje (een zware halo) trok automatisch meer bomen (galaxieën) aan dan een klein bosje. Het leek een simpele regel: Hoe zwaarder de halo, hoe meer sterrenstelsels erin zitten.

Maar deze nieuwe studie, geschreven door Hong-Gang Yang en collega's, vertelt ons dat het iets ingewikkelder is. Het is alsof je merkt dat twee bomen van precies dezelfde grootte, in twee verschillende buurten, toch een heel ander aantal bomen hebben. Waarom? Omdat de omgeving er toe doet. Dit fenomeen noemen ze Galaxy Assembly Bias (GAB).

Hier is wat ze hebben ontdekt, vertaald naar alledaagse taal:

1. Twee Simulaties, Twee Verhalen

De onderzoekers keken naar twee van de meest geavanceerde computersimulaties van het heelal die we hebben: SIMBA en IllustrisTNG.

De Analogie: Stel je voor dat twee verschillende koks (SIMBA en TNG) proberen dezelfde soep te maken met dezelfde ingrediënten (de wetten van de natuurkunde), maar ze gebruiken verschillende recepten voor hoe ze de kruiden (de fysica van gas en sterren) toevoegen.
Het Resultaat: Tot een bepaald punt in de tijd (toen het heelal ongeveer half zo oud was als nu, zo'n 10 miljard jaar geleden), gedroegen de twee soepen zich hetzelfde. Maar daarna gingen ze uit elkaar.
- In IllustrisTNG bleef het effect van de "buurt" (de omgeving) sterk toenemen naarmate het heelal ouder werd.
- In SIMBA verdween dit effect juist weer; de sterrenstelsels leken zich weer alleen te laten leiden door de grootte van hun halo.

Dit betekent dat de manier waarop sterrenstelsels zich groeperen, sterk afhankelijk is van hoe we de "fysica" in onze modellen beschrijven.

2. De Dichte Buurt vs. De Tijden

Wat is dan die "omgeving"? De onderzoekers keken naar twee dingen:

De Dikte van de Buurt (Overdichtheid): Zit je in een drukke stad (veel materie) of op een afgelegen eiland (weinig materie)?
De Vorm van de Buurt (Tidal Anisotropy): Is de buurt rond je een ronde kring of een langgerekte strook?

De Grote Ontdekking:
Het bleek dat de dichtheid (hoe druk het is in de buurt) de belangrijkste drijver is.

Analogie: Stel je voor dat je in een drukke stad woont (een hoge dichtheid). Dan is de kans groter dat je een huis bouwt, zelfs als je grondstuk (de halo) niet heel groot is. In een dorpje (lage dichtheid) bouw je misschien geen huis op een even groot stuk grond.
De "vorm" van de buurt (de tijdelijke anisotropie) bleek bijna geen rol te spelen. Het maakt niet uit of de buurt rond of langwerpig is; het gaat puur om hoe druk het er is.

3. De "Kleine" Halos zijn de Sleutel

Een verrassend detail: dit effect gebeurt vooral bij de kleine halo's.

De Analogie: Grote, zware halo's zijn als enorme stadskernen. Die zijn altijd druk, ongeacht waar ze zitten. Maar kleine halo's zijn als kleine dorpjes. Een klein dorpje in een drukke regio (hoge dichtheid) heeft veel meer kans om een sterrenstelsel te "bevatten" dan een klein dorpje in een afgelegen regio.
De onderzoekers ontdekten dat als je alleen naar de grote halo's kijkt, er geen sprake is van dit "buurt-effect". Het effect zit hem volledig in de kleinere, lichtere halo's.

4. Waarom is dit belangrijk?

Vandaag de dag proberen wetenschappers het heelal te begrijpen door naar miljoenen sterrenstelsels te kijken (zoals met de nieuwe telescopen DESI en Euclid). Als ze denken dat alleen de grootte van de halo telt, maken ze een fout in hun berekeningen.

Het Gevolg: Als je dit "buurt-effect" negeert, kun je de manier waarop het heelal is opgebouwd verkeerd interpreteren. Het kan leiden tot een foutieve schatting van hoeveel donkere energie er is of hoe snel het heelal uitdijt.
De Oplossing: De auteurs hebben een nieuwe formule bedacht. In plaats van te zeggen: "In een halo van grootte X zitten Y sterrenstelsels", zeggen ze nu: "In een halo van grootte X zitten Y sterrenstelsels, maar als die halo in een drukke buurt zit, zijn het er misschien Z."

Samenvatting in één zin

Deze studie laat zien dat sterrenstelsels niet alleen worden bepaald door hun eigen gewicht, maar ook door hoe druk het er om hen heen is, en dat verschillende computermodellen hierover nog steeds heel verschillend denken – wat ons waarschuwt om voorzichtig te zijn met onze conclusies over het heelal.

Het is alsof je merkt dat je niet alleen bepaalt wie je vrienden zijn door je eigen karakter, maar ook door in welke wijk je woont. En in de kosmos maakt die wijk een groot verschil in hoe het universum eruitziet.

Probleemstelling

De ruimtelijke verdeling van sterrenstelsels is een cruciale observabele voor de kosmologie, omdat deze de onderliggende grote-schaalstructuur (LSS) van het universum traceert. In het kader van hiërarchische structuurvorming worden sterrenstelsels gevormd binnen donkere-materiehalo's. Traditionele modellen, zoals de Halo Occupation Distribution (HOD), gaan ervan uit dat de bezetting van een halo uitsluitend wordt bepaald door de massa van die halo.

Echter, dit veronderstelde onafhankelijkheid van de kosmische omgeving is in twijfel getrokken door het fenomeen van Galaxy Assembly Bias (GAB). GAB verwijst naar de afhankelijkheid van de clustering van sterrenstelsels van halo-eigenschappen die verder gaan dan alleen massa (zoals vormingstijd, concentratie, of de lokale omgeving). Hoewel GAB in simulaties is waargenomen, blijft de fysieke oorsprong ervan onduidelijk, en er is geen consensus over de omvang van het effect of welke secundaire eigenschappen (zoals dichtheid of getijdenkrachten) de belangrijkste drijvende factoren zijn. Dit is een kritiek probleem voor de precisie-kosmologie van de toekomst (bijv. DESI, Euclid), omdat het negeren van GAB systematische fouten in de interpretatie van data kan veroorzaken.

Methodologie

De auteurs onderzoeken de sterkte en de roodverschuivings-evolutie van GAB door twee state-of-the-art kosmologische hydrodynamische simulaties te vergelijken: SIMBA en IllustrisTNG. Hoewel beide simulaties de waargenomen sterrenmassa-functie goed reproduceren, gebruiken ze fundamenteel verschillende subgrid-prescripties voor baryonische fysica, met name voor feedbackprocessen (sterren en AGN).

De kernmethodologie omvat:

Shuffle-techniek: Om GAB te kwantificeren, worden sterrenstelsels willekeurig herverdeeld (geshuffled) naar halo's binnen dezelfde massa-bin. Dit verwijdert elke afhankelijkheid van secundaire halo-eigenschappen. De GAB-signaal wordt gedefinieerd als de verhouding tussen de clustering van de originele catalogus en de geshuffled catalogus: $[w_p(r_p)/w_p(r_p)_{shuffled} - 1]$ .
Omgevingskarakterisering: De auteurs analyseren de bijdrage van twee specifieke omgevingsparameters:
- Gesmoorde materie-overdichtheid ( $\delta_5$ ): Gemeten binnen een straal van 5 $h^{-1}$ Mpc.
- Getijdenanisotropie ( $\alpha_5$ ): Een maat voor de vorm van het getijdenveld.
Dubbel-shuffling: Om de bijdrage van deze parameters te isoleren, worden sterrenstelsels herverdeeld binnen bins van zowel halo-massa als $\delta_5$ (of $\alpha_5$ ). Als het GAB-signaal hierdoor verdwijnt, wordt de bias volledig toegeschreven aan die parameter.
HOD-analyse: De auteurs onderzoeken hoe het gemiddelde aantal centrale en satellietstelsels varieert met de halo-massa en de omgevingsdichtheid.
Massa-gesplitste analyse: Er wordt gekeken naar de bijdrage van lage-massa versus hoge-massa halo's aan het totale GAB-signaal.

Belangrijkste Resultaten

1. Divergerende evolutie tussen SIMBA en TNG

Bij hoge roodverschuiving ( $z=5$ ) is GAB in beide simulaties verwaarloosbaar.
Tegen $z=2$ neemt GAB toe tot ongeveer 5% in beide simulaties.
Na $z=2$ $z = 2$ vertonen de simulaties sterk verschillende gedragingen:
- SIMBA: Het GAB-signaal neemt af en is bij $z=0$ bijna nul.
- IllustrisTNG: Het GAB-signaal blijft toenemen en bereikt ongeveer 10% bij $z=0$ .
De auteurs suggereren dat dit verschil wordt gedreven door de verschillende implementaties van AGN-feedback. SIMBA gebruikt stabiele, bipolaire jets die gas verplaatsen, terwijl TNG willekeurige, niet-ontkoppelde jets gebruikt. Dit leidt tot verschillen in de baryonfractie en het circum-galactische medium.

2. Dominantie van Overdichtheid ( $\delta_5$ )

De dubbel-shuffle analyse toont aan dat de gesmoorde overdichtheid ( $\delta_5$ ) verantwoordelijk is voor het merendeel van het GAB-signaal in beide simulaties.
De bijdrage van getijdenanisotropie ( $\alpha_5$ ) is verwaarloosbaar in SIMBA en slechts klein (maar niet verwaarloosbaar) in TNG. Dit suggereert dat de lokale dichtheid de belangrijkste omgevingsfactor is voor GAB.

3. Rol van Lage-Massa Halo's

De HOD-analyse onthult dat de afhankelijkheid van de omgeving voornamelijk optreedt bij lage-massa halo's (rond de "turn-on" massa waar stelsels beginnen te vormen).
In hoge-massa halo's is de inhoud van stelsels vrijwel onafhankelijk van de omgeving.
Shuffling-tests bevestigen dat bijna het volledige GAB-signaal afkomstig is van lage-massa halo's ( $M_h < 10^{11.5} M_\odot$ ). Halo's in dichtere gebieden hebben een hogere kans om stelsels te herbergen, maar dit effect is sterker in TNG dan in SIMBA.

4. Invloed van Sterrenmassa-drempel

Een analyse met een hogere drempel voor sterrenmassa ( $>10^{10} M_\odot$ ) in TNG toont aan dat het GAB-signaal iets afneemt, maar dat de trend (toename naar lage $z$ ) en de dominantie van $\delta_5$ behouden blijven.

Bijdrage en Parametrisatie

De auteurs presenteren een nieuwe parametrisatie voor de HOD die rekening houdt met de omgevingsdichtheid. In plaats van alleen de massa te gebruiken, modelleren ze de helling van de HOD-curve ( $\mu$ en $\nu$ ) als lineair afhankelijk van de rangorde van de overdichtheid ( $\delta_5$ ). Dit stelt onderzoekers in staat om realistischere "mock" catalogi te genereren die de effecten van GAB incorporeren, wat essentieel is voor de analyse van toekomstige surveys.

Betekenis en Conclusie

De studie concludeert dat Galaxy Assembly Bias een significant effect heeft op de clustering van sterrenstelsels (tot ~10% in amplitude) en dat dit effect sterk afhankelijk is van de gebruikte baryonische fysica in simulaties.

Fysische Oorsprong: De belangrijkste drijvende kracht is de verhoogde kans dat lage-massa halo's in dichtere gebieden stelsels vormen. Dit effect wordt echter onderdrukt door AGN-feedback uit nabijgelegen massieve halo's, wat verklaart waarom SIMBA en TNG verschillende resultaten geven.
Observationele Impact: Het negeren van GAB kan leiden tot fouten in de schatting van kosmologische parameters. De auteurs waarschuwen dat de discrepantie tussen hun simulaties en sommige observationele claims (die een grotere rol voor getijdenanisotropie suggereren) verder onderzoek vereist.
Toekomst: De ontwikkelde parametrisatie biedt een praktische tool om GAB te modelleren in HOD-gebaseerde simulaties, waardoor de nauwkeurigheid van interpretaties van data van instrumenten zoals DESI en Euclid wordt verbeterd.