Impact of subhalo dynamical friction heating on the formation of the first structures in the universe

Titel: Impact van dynamische wrijvingsverwarming door subhalo's op de vorming van de eerste structuren in het universum

Auteurs: Zhenyu Wu et al.
Publicatie: MNRAS (voorbereid voor publicatie in 2026)

1. Het Probleem

De thermische geschiedenis van baryonische materie in donkere-materiehopen (dark matter halos) is de bepalende factor voor sterformatie. In het vroege universum (bij hoge roodverschuivingen, $z \sim 15$) vormen zich de eerste sterren (Populatie III) in minihalos. Voor deze formatie moet gas kunnen afkoelen, voornamelijk via moleculaire waterstof ($H_2$).
Echter, verschillende feedbackmechanismen kunnen deze afkoeling onderdrukken en sterformatie vertragen of stoppen, wat de vorming van zware zwarte gaten (Direct Collapse Black Holes - DCBH) kan bevorderen. Hoewel processen zoals AGN-feedback en supernova's goed bestudeerd zijn in massieve halos, is de bijdrage van dynamische wrijving (Dynamical Friction - DF) van orbitale subhalo's aan gasverwarming in het vroege universum minder goed gekwantificeerd. De vraag is of deze verwarming voldoende sterk kan zijn om $H_2$-koeling te compenseren en zo de vorming van Pop III-sterren te onderdrukken.

2. Methodologie

De auteurs ontwikkelen een analytisch model om de verwarmingsrate door dynamische wrijving te berekenen en valideren dit met data uit de TNG50-simulatie (een hydrodynamische cosmologische simulatie).

• Analytisch Model:
  • Ze baseren zich op het model van Ostriker (1999) voor dynamische wrijving in een gassenmedium. De verwarmingsrate wordt berekend door de som van de wrijvingskrachten van alle subhalo's in een gastheerhalo.
  • De formule voor de verwarmingsrate ($\dot{E}_{DF}$) integreert de subhalo-massafunctie (SHMF) over de massa's van de subhalo's, waarbij rekening wordt gehouden met de gasdichtheid, de geluidssnelheid en een correctiefactor voor het Mach-getal.
• Subhalo Massa Functie (SHMF):
  • In plaats van standaard fittingen te gebruiken, deriveerden de auteurs de SHMF direct uit de TNG50-12 simulatie voor een breed scala aan roodverschuivingen ($z=15$ tot $z=0$).
  • Ze passen een fittingformule toe (een machtsfunctie met exponentiële afname) en analyseren de evolutie van de parameters ($\alpha, A, \beta, \omega$) als functie van de roodverschuiving.
  • Ze bestuderen de streuung (scatter) in het aantal subhalo's tussen verschillende halos, wat cruciaal is omdat de verwarming sterk afhankelijk is van het specifieke aantal subhalo's. Ze nemen aan dat deze fluctuaties Poisson-gedistribueerd zijn.
• Correcties:
  • Ze evalueren niet-lineariteitseffecten (wanneer de gravitatie van de subhalo sterk is ten opzichte van de gasdruk) en Mach-getal-correcties. Ze concluderen dat niet-lineariteit voor de meeste subhalo's een verwaarloosbaar effect heeft op de statistiek.
• Vergelijking:
  • De berekende verwarmingsrates worden vergeleken met koelingsrates (berekend met de Grackle-chemie- en koelingsbibliotheek) voor zowel massieve halos bij lage $z$ als minihalos bij hoge $z$.
  • Voor het vroege universum wordt de kritische $H_2$-fractie bepaald waarbij de koeltijd gelijk is aan de leeftijd van het universum, zowel met als zonder DF-verwarming.

3. Belangrijkste Bijdragen

1. Een nieuw model voor DF-verwarming: Een analytisch kader dat de verwarming door subhalo's kwantificeert op basis van een empirisch afgeleide SHMF uit TNG50, geldig van $z=15$ tot $z=0$.
2. Kwantificering van de variabiliteit: Het inzicht dat de verwarmingsrate voor een vaste gastheerhalo-massa kan variëren met twee orde van grootte door de natuurlijke spreiding in het aantal subhalo's (Poisson-fluctuaties).
3. Resolutie-eisen: De vaststelling dat om DF-verwarming correct te simuleren, subhalo's moeten worden opgelost met een massa-ratio $\psi = m_{sub}/M_{host} \gtrsim 0.05$.
4. Impact op de eerste sterren: Het aantonen dat DF-verwarming een significante rol speelt in de vorming van de eerste structuren, specifiek door de drempel voor Pop III-sterformatie te verhogen.

4. Resultaten

• Lage roodverschuiving ($z \lesssim 2$) en Massieve Halos:

  • In massieve halos ($M > 10^{13} M_\odot$) is DF-verwarming een belangrijk mechanisme voor het "quenchen" (onderdrukken) van sterformatie.
  • De verwarming is vergelijkbaar met of groter dan de koelingsrate, wat consistent is met eerdere studies over gravitationele verwarming.
  • Er is echter een grote spreiding: in halos met weinig subhalo's is DF-verwarming onvoldoende om de centrale sterformatie volledig te onderdrukken.

• Hoge roodverschuiving ($z \sim 15$) en Minihalos:

  • In metal-vrije minihalos ($M \sim 10^6 - 10^8 M_\odot$) kan de DF-verwarmingsrate vergelijkbaar zijn met de koelingsrate van moleculaire waterstof ($H_2$).
  • Dit verwarmingsmechanisme kan de afkoeling en fragmentatie van gas onderdrukken.
  • Kritieke $H_2$-fractie: Door DF-verwarming neemt de kritieke $H_2$-fractie ($f_{crit}$) die nodig is voor Pop III-sterformatie toe met tot één orde van grootte. Dit maakt de vorming van Pop III-sterren aanzienlijk moeilijker.

• Directe Inzinking Zwarte Gaten (DCBH):

  • DF-verwarming werkt synergistisch met het Lyman-Werner (LW) stralingsachtergrond. Zelfs als de gemiddelde $H_2$-fractie in de voorouder-minihalos laag is ($\sim 10^{-5}$), kan DF-verwarming de resterende $H_2$-koeling onderdrukken.
  • Dit bevordert de vorming van DCBH-zaden in atoomkoelende halos, zelfs in omstandigheden waar Pop III-sterformatie anders mogelijk zou zijn geweest.
  • De resultaten suggereren dat DCBH-vorming mogelijk is in een breder bereik van omstandigheden dan alleen bij zeer sterke LW-straling, mits er voldoende subhalo's aanwezig zijn om dynamische wrijving te genereren.

5. Betekenis en Conclusie

Dit onderzoek identificeert dynamische wrijving door subhalo's als een cruciaal, maar vaak verwaarloosd mechanisme in de evolutie van het vroege universum.

• Voor Pop III-sterren: Het mechanisme verhoogt de drempel voor hun vorming, wat betekent dat ze minder vaak zullen ontstaan dan modellen zonder DF-verwarming voorspellen.
• Voor DCBH's: Het biedt een natuurlijke, continue route voor de vorming van zware zwarte gatenzaden in het ΛCDM-model, zonder dat er noodzakelijkerwijs extreme pieken in externe straling (LW) nodig zijn.
• Voor Simulaties: De studie benadrukt dat hydrodynamische simulaties voldoende massa-resolutie moeten hebben om subhalo's met een massa-ratio van $\psi \gtrsim 0.05$ op te lossen. Zonder deze resolutie wordt een belangrijk verwarmingsmechanisme gemist, wat leidt tot onnauwkeurige voorspellingen over de vorming van de eerste sterren en zwarte gaten.

Samenvattend suggereert dit werk dat dynamische wrijving een fundamentele rol speelt in het reguleren van de thermische geschiedenis van het vroege universum en de overgang van Pop III-sterformatie naar de vorming van zware zwarte gaten.