The host halo masses of AGNs and quasars at $z \sim 3-7$ with TNG-Cluster, FLAMINGO and other cosmological galaxy simulations

Met behulp van grootschalige kosmologische simulaties onthult deze studie dat hoewel de AGN-luminositeit over het algemeen toeneemt met de massa van de gast-halo tot een drempelwaarde, de relatie sterk niet-lineair is met aanzienlijke spreiding, wat aangeeft dat de meest lumineuze quasars bij $z \sim 3-7$ doorgaans zich in halo's met een intermediaire massa ($10^{12-12.5}$ M$_{\odot}$) bevinden in plaats van in de meest massieve, een bevinding die goed overeenkomt met observationele schattingen.

De kern

De Grote Visie: Het Zoeken naar de "Sweet Spot" voor Kosmische Supersterren

Stel je het universum voor als een gigantische, groeiende stad. In deze stad zijn de "donkere materie halo's" de wijken, en de "Active Galactic Nuclei" (AGN's) of "Quasars" zijn de enorme, verblindend heldere vuurtorens in het centrum van de belangrijkste gebouwen. Deze vuurtorens worden aangedreven door supermassieve zwarte gaten die gas verslinden.

Al een lange tijd proberen astronomen een simpele vraag te beantwoorden: Hoe groot moet een wijk (halo) zijn om de helderste vuurtoren (quasar) te huisvesten?

De meeste waarnemingen suggereren dat deze superheldere quasars in wijken wonen die ongeveer even groot zijn, ongeacht hoe ver je terugkijkt in de tijd (of hoe ver weg je kijkt). Ze lijken een "Goldilocks"-zone te verkiezen: niet te klein, niet te groot.

Dit paper gebruikt enorme computersimulaties om dat idee te testen. De onderzoekers bouwden digitale universa (met behulp van modellen zoals TNG-Cluster en FLAMINGO) om te zien of hun virtuele quasars zich gedragen als de echte.

De Belangrijkste Ontdekking: Het is Geen Rechte Lijn

De onderzoekers ontdekten dat de relatie tussen de grootte van de wijk en de helderheid van de vuurtoren geen rechte lijn is. Het is meer als een heuvel met een piek.

1. De Klim: Naarmate de wijk groter wordt (massiever), wordt de vuurtoren over het algemeen helderder. Dit is logisch; grotere wijken hebben meer gas om het zwarte gat te voeden.
2. De Piek: Deze trend gaat alleen omhoog tot een bepaald punt (ongeveer $10^{12}$ keer de massa van onze Zon).
3. De Daling: Zodra de wijk te massief wordt, wordt de vuurtoren zelfs minder helder of blijft deze gelijk.

De Analogie: Denk aan een feestje.

• In een klein huisje (kleine halo) zijn er maar een paar mensen, en is de muziek zacht.
• In een middelgroot huis (de sweet spot) heb je het perfecte publiek, en is het feest luid en energiek.
• In een enorm stadion (een enorme halo) wordt het feestje juist stiller. Waarom? Omdat de "gastheer" (het zwarte gat) zo gulzig en krachtig wordt dat hij de gasten letterlijk de kamer uit blaast (dit wordt AGN-feedback genoemd). Hij duwt het gas weg, waardoor hij zichzelf uithongert en het feestje niet luider kan worden.

De "Chaos"-factor: Enorme Spreiding

Een van de meest verrassende bevindingen is hoe rommelig de gegevens zijn. Het paper vergelijkt twee manieren om naar de gegevens te kijken:

• Scenario A: Als je een specifieke wijkgrootte kiest, hoe helder is de vuurtoren dan?
  • Resultaat: Wild onvoorspelbaar. De helderheid kan variëren met een factor 1.000 tot 10.000 (3 tot 4 "decades" verschil). Een huis van dezelfde grootte kan een klein nachtlampje hebben, terwijl een ander een verblindende spotlight heeft.
• Scenario B: Als je een specifieke helderheid kiest (bijv. een superheldere quasar), hoe groot is dan de wijk?
  • Resultaat: Veel voorspelbaarder. Als je een superheldere vuurtoren ziet, bevindt deze zich bijna zeker in een wijk van een specifieke grootteklasse.

De Analogie: Stel je voor dat je probeert te raden hoeveel een persoon weegt door alleen naar de lengte te kijken.

• Als je een specifieke lengte kiest (bijv. 1,80 meter), kunnen mensen ergens tussen de 65 kg en 120 kg wegen. Dat is een enorme spreiding.
• Maar als je een specifiek gewicht kixt (bijv. 90 kg), is die persoon bijna zeker rond de 1,80 meter lang. Het gewicht voorspelt de lengte veel beter dan de lengte het gewicht voorspelt.

Het paper concludeert dat de grootte van de wijk een slechte voorspeller is van hoe helder het zwarte gat op dit moment is. De helderheid van het zwarte gat is chaotisch en hangt af van veel andere dingen (zoals hoeveel gas er momenteel beschikbaar is), niet alleen van de grootte van de wijk.

De "Sweet Spot" Bevestigd

Ondanks de chaos bevestigden de simulaties het observationele idee: Quasars verkiezen een specifieke wijkgrootte.

• Of het universum nu 3 miljard jaar oud is of 7 miljard jaar oud, de helderste quasars wonen meestal in halo's met een massa van ongeveer $10^{12}$ tot $10^{12,5}$ zonnemassa's.
• Ze wonen zelden in de meest massieve halo's (de supersteden).
• Ze wonen zelen zelden in de kleinste halo's (de kleine dorpjes).

Dit suggereert dat er een "sweet spot" is voor de groei van zwarte gaten. Als de wijk te klein is, kan het zwarte gat niet genoeg eten krijgen. Als de wijk te groot is, stopt de eigen feedback van het zwarte gat (het wegblazen van gas) het proces van groei.

Waarom Dit Belangrijk Is

Het paper betoogt dat we niet simpelweg naar de grootte van de wijk van een sterrenstelsel kunnen kijken om te raden hoe helder het zwarte gat is. De connectie is te los en te chaotisch. Echter, als we een superheldere quasar zien, kunnen we vrij zeker zijn over de grootte van de woning.

De simulaties laten ook zien dat deze "sweet spot" gedurende de kosmische geschiedenis ongeveer dezelfde grootte behoudt, ook al zet het universum uit en verandert het. Dit impliceert dat de regels die het voeden van zwarte gaten en het afsluiten ervan reguleren, consistent zijn over miljarden jaren heen.

Samenvatting in één zin

Het paper stelt vast dat hoewel supermassieve zwarte gaten in een specifieke "Goldilocks"-grootte van kosmische wijk leven, de helderheid van het zwarte gat zo chaotisch is dat weten hoe groot de wijk is, je heel weinig zegt over hoe helder het zwarte gat op een gegeven moment is.

Technische samenvatting

Probleemstelling
Observationele studies en clusteringanalyses suggereren al lang dat quasars een nauw bereik van donkere materie-halo-massa's ($\sim 10^{12} - 10^{13} \, M_\odot$) bewonen over de kosmische tijd heen ($z \lesssim 7$). Deze gevolgtrekkingen zijn echter vaak indirect en rusten op flux-gelimiteerde steekproeven en aannames over de monotoniciteit van de relatie tussen de AGN-luminositeit en de gast-halo-massa. Er bestaat een kritische kloof in het begrip van de intrinsieke spreiding en de precieze aard van deze relatie vóór "Cosmic Noon" ($z \gtrsim 3$), met name met betrekking tot de vraag of de meest lumineuze AGNs zich in de meest massieve halo's bevinden en hoe feedbackmechanismen deze verbinding reguleren. Eerdere simulaties waren vaak beperkt door volume (niet in staat om zeldzame, massieve halo's te bemonsteren) of resolutie (niet in staat om kleinschalige gasdynamica die AGN-luminositeiten aanstuurt, op te lossen).

Methodologie
Dit werk maakt gebruik van een suite van kosmologische hydrodynamische simulaties om observationele beelden direct te vergelijken met zelfconsistente sterrenstelselvormingsmodellen. De studie richt zich op AGNs met bolometrische luminositeiten $L_{\text{AGN}}^{\text{bol}} \geq 10^{42} \, \text{erg s}^{-1}$ bij roodverschuivingen $z = 3 - 7$.

De primaire analyse steunt op twee grote volume-simulatiesuite:

1. TNG300 + TNG-Cluster: Onderdeel van het IllustrisTNG-project, biedt een hoge resolutie (baryonische doelmassa $\sim 1.1 \times 10^7 \, M_\odot$) over een volume van $300^3 \, \text{cMpc}^3$ en een zoom-in extensie gericht op massieve halo's tot $10^{15} \, M_\odot$.
2. FLAMINGO: Omvat volumes van $1 - 3 \, \text{Gpc}^3$ per zijde (runs L1_m8 en L2p8_m9), wat zorgt voor een superieure statistische bemonstering van het hoog-massa-einde van de halo-massa-functie, zij het met een lagere resolutie dan TNG.

Voor benchmarking bevatten de auteurs kleinere-volume ($\sim 100 \, \text{cMpc}^3$) simulaties uit eerdere generaties: Illustris, EAGLE, TNG100 en Simba.

AGN-bolometrische luminositeiten zijn geen directe simulatie-outputs, maar worden in de post-processing berekend vanuit de instantane SMBH-accretiesnelheden met behulp van een bimodaal conversiemodel (Churazov et al. 2005) dat onderscheid maakt tussen radiatief efficiënte en inefficiënte accretieregimes. Een uniforme radiatieve efficiëntie van $\epsilon_r = 0.1$ wordt in alle modellen gehanteerd om verschillen in de onderliggende accretiefysica te isoleren.

Belangrijkste Bijdragen en Resultaten

• Niet-lineaire en niet-monotone Relaties: In alle modellen is de mediane relatie tussen de AGN-bolometrische luminositeit en de gast-halo-massa sterk niet-lineair. Hoewel meer massieve halo's gemiddeld meer lumineuze AGNs huisvesten, houdt deze trend alleen stand tot een karakteristieke halo-massa van $M_{200,\text{crit}} \sim 10^{12} \, M_\odot$.

  • In TNG300+TNG-Cluster vlakt de mediane relatie af en vertoont zelfs een omslag (dip) bij $M_{200,\text{crit}} \gtrsim 10^{12} \, M_\odot$ voor $z < 5-6$.
  • In FLAMINGO vlakt de relatie af aan het hoog-massa-einde, maar vertoont niet dezelfde uitgesproken omslag.
  • Deze quenching bij hoge massa wordt toegeschreven aan AGN-feedback (kinetisch in TNG, thermisch in FLAMINGO), die gasaccretie naar SMBHs in massieve halo's onderdrukt.

• Grote Intrinsieke Spreiding: Het artikel kwantificeert een significante object-tot-object variatie. De spreiding in de AGN-bolometrische luminositeit bij een vaste gast-halo-massa is veel groter (tot 3 dex voor de 5e–95e percentielen) dan de spreiding in de gast-halo-massa bij een vaste AGN-luminositeit ($\lesssim 1$ dex). Dit impliceert een zwakke koppeling tussen halo-groei en instantane SMBH-accretie, wat suggereert dat halo-massa alleen een slechte voorspeller is van de instantane AGN-luminositeit.

• Quasar Gast-halo-massa's:

  • TNG300+TNG-Cluster: Quasars ($L_{\text{AGN}}^{\text{bol}} \sim 10^{45} - 10^{47} \, \text{erg s}^{-1}$) verblijven doorgaans in halo's van $10^{12} - 10^{12.5} \, M_\odot$ bij $z=3-6$.
  • FLAMINGO: Quasars strekken zich uit tot mediane massa's van $\sim 10^{12.8} \, M_\odot$ bij $z \sim 3-4$.
  • Cruciaal is dat de meest lumineuze AGNs niet de meest massieve halo's huisvesten die op elk gegeven moment beschikbaar zijn tot $z \approx 7$. In plaats daarvan verblijven ze in progressief zeldzamere halo's bij eerdere tijdperken, maar vermijden ze typisch de zeer hoog-massa systemen ($M_{200,\text{crit}} \gtrsim 10^{13} \, M_\odot$) vanwege feedback-quenching.

• Vergelijking met Observaties: De gesimuleerde gast-halo-massa's voor quasars komen in grote lijnen overeen met observationele schattingen afgeleid van clustering en cross-correlatiemetingen, die over het algemeen karakteristieke massa's van $10^{12} - 10^{13} , M_\odot impliceren. De grote spreiding voorspeld door de simulaties helpt de brede observationele beperkingen te verzoenen.

Betekenis en Claims
Het artikel claimt dat moderne kosmologische simulaties, ondanks verschillen in resolutie, volume en sub-grid fysica (seeding, accretie, feedback), convergeren naar een consistente visie: de AGN-luminositeit–halo-massa relatie wordt gekenmerkt door een "sweet spot" van een karakteristieke halo-massa ($\sim 10^{12} \, M_\odot$) waar quasars het meest prevalent zijn, geflankeerd door quenching in meer massieve halo's en onderdrukking in minder massieve halo's.

De auteurs benadrukken dat de grote intrinsieke spreiding in AGN-luminositeit bij een vaste halo-massa de aannames van veel semi-empirische halo-bezettingmodellen uitdaagt, die vaak een nauwere relatie veronderstellen. De resultaten suggereren dat hoewel de typische quasar zich in een specifieke massabericht bevindt, individuele systemen een enorme diversiteit vertonen, en dat de meest massieve halo's niet noodzakelijkerwijs de gastheren zijn van de meest lumineuze AGNs op elk gegeven moment. Dit werk biedt een kwantitatieve baseline voor het interpreteren van hoog-roodverschuiving quasar-observaties en benadrukt de rol van AGN-feedback bij het vormgeven van de demografie van de meest lumineuze zwarte gaten vóór Cosmic Noon.