Bloodhound Unleashed: Particle-based Substructure Tracking for Cosmological Simulations

Titel: Bloodhound Unleashed: Deeltjesgebaseerde Substructuur-tracking voor Kosmologische Simulaties

Auteurs: Hyunsu Kong, Michael Boylan-Kolchin, en James S. Bullock
Publicatie: MNRAS (2026)

---

1. Het Probleem

Moderne studies over galaxievorming en de aard van donkere materie zijn afhankelijk van numerieke simulaties. Een cruciale uitdaging in dit veld is het nauwkeurig identificeren en volgen van zelfgebonden structuren (subhalo's) van sterren en donkere materie terwijl ze in de zware gravitatieput van een gastheerhalo (zoals de Melkweg) bewegen.

De standaardmethodiek, die vaak de combinatie ROCKSTAR (een halo-vinder per snapshot) en consistent-trees (een mergertree-bouwer) gebruikt, ondervindt ernstige beperkingen in dichte omgevingen:

• Vroegtijdig verlies van subhalo's: Subhalo's die diep in de gastheerhalo vallen, worden vaak "verloren" door de tracker, zelfs als ze nog steeds een aanzienlijke massa en een zelfgebonden kern hebben. Dit gebeurt vaak na peri-centrische passages (dichtste benadering van het centrum).
• Gebroken links (Broken-link trees): Wanneer een subhalo tijdelijk uit de catalogus verdwijnt en later weer verschijnt, wordt deze soms niet gekoppeld aan zijn oorspronkelijke geschiedenis. Dit leidt tot de schijnbare vorming van nieuwe subhalo's binnen de gastheerhalo, wat in strijd is met het $\Lambda$CDM-model.
• Onnauwkeurige statistieken: Deze fouten leiden tot een onderschatting van het aantal overlevende subhalo's, vooral in de binnenste regio's (< 50 kpc) van de Melkweg. Dit beïnvloedt conclusies over de "missing satellites"-problematiek en de eigenschappen van donkere materie.

2. Methodologie: Het Bloodhound-algoritme

Bloodhound is geen nieuwe halo-vinder, maar een post-processing tool die bestaande pipelines (zoals ROCKSTAR + consistent-trees) aanvult. In plaats van te proberen subhalo's per snapshot te identificeren en hun geschiedenis te "stikken" (stitching), volgt Bloodhound de deeltjes direct.

Kernprincipes:

1. Deeltjesselectie bij instorting (Infall):

   • Het algoritme identificeert een subhalo op het moment vlak voordat deze de viriale straal ($R_{vir}$) van de gastheerhalo binnendringt.
   • Het selecteert de 2% meest gebonden deeltjes van die subhalo op dat moment (met een minimum van 50 en maximum van 500 deeltjes).
   • Deze specifieke set deeltjes wordt vervolgens via hun unieke ID's door de tijd gevolgd tot $z=0$.

2. Positiebepaling:

   • Het zwaartepunt (Center of Mass) wordt berekend door iteratief een krimpende bol te gebruiken rondom de meest gebonden deeltjes. Dit vermijdt de verwarring met losse deeltjes in de dichte achtergrond van de gastheer.

3. Scheiding van gebonden en gestripte deeltjes:

   • Om de eigenschappen van de subhalo te berekenen, moet onderscheid worden gemaakt tussen de kern en gestripte deeltjes (tidal tails).
   • Bloodhound analyseert het dichtheidsprofiel $\rho(r)$. Het zoekt naar het punt waar de helling van het differentieel massaprofiel ($\rho(r) \times r^2$) verandert van negatief naar positief. Dit punt ($R_{trunc}$) markeert de grens waar gestripte, ongebonden materie begint.

4. Disruptie-criteria (Wanneer is een subhalo "vernietigd"?):

   • In plaats van alleen te kijken naar het aantal deeltjes, gebruikt Bloodhound de fysische concentratieparameter $c_V$ (de verhouding tussen de gemiddelde dichtheid binnen $R_{max}$ en de kritische dichtheid van het universum).
   • Een subhalo wordt als verstoord beschouwd als:
     1. $c_V$ abrupt daalt tot < 30% van de vorige snapshot.
     2. $c_V$ daalt met een factor van minstens 5 ten opzichte van de instortingswaarde.
     3. Deze lage waarde gedurende 5 opeenvolgende snapshots ($\sim$500 Myr) aanhoudt (om tijdelijke getijdenchokken uit te sluiten).

3. Data en Simulaties

De auteurs gebruiken de Phat ELVIS-suite van simulaties:

• Type: Cosmologische "zoom-in" simulaties van Melkweg-massa halo's ($M_v \approx 10^{12} M_\odot$).
• Vergelijking: Er zijn twee sets:
  1. DMO (Dark Matter Only): Alleen donkere materie.
  2. Disc: Dezelfde simulaties, maar met een toegevoegd potentieel voor een centrale sterrenschijf (om de verstoring door de Melkweg na te bootsen).
• Selectie: Er wordt gekeken naar subhalo's die in de standaardmethode als "verstoord" worden gemarkeerd, met een instortingsroodverschuiving $z_{infall} \le 3$ en een maximale rotatiesnelheid bij instorting $10 \le V_{infall} \lesssim 80$ km/s.

4. Belangrijkste Resultaten

• Verlengde Tracking-tijd:

  • Bloodhound volgt subhalo's gemiddeld 3 tot 4 miljard jaar (Gyr) langer dan de standaardmethode.
  • Voor subhalo's met zeer kleine pericentra (dicht bij het centrum) is dit verschil nog groter.
  • In de Disc-simulaties overleeft ongeveer 25% van de subhalo's die door de standaardmethode als vernietigd worden gemarkeerd, tot $z=0$ volgens Bloodhound. In DMO-simulaties is dit bijna 50%.

• Correctie van "Broken-link" Subhalo's:

  • Bloodhound lost het probleem van subhalo's die lijken te ontstaan binnen de gastheerhalo op. Het koppelt deze correct terug aan hun oorspronkelijke, zwaardere voorouders die eerder werden "verloren".
  • Dit resulteert in een significant hogere geschatte massa voor deze objecten. Subhalo's die in de standaardmethode een lage $V_{infall}$ (< 10 km/s) lijken te hebben, blijken vaak oorspronkelijk veel massiever te zijn geweest ($V_{infall} > 17$ km/s), wat relevant is voor het vormen van sterrenstelsels.

• Toename van Overlevende Subhalo's in het Centrum:

  • Binnen 50 kpc van het centrum van de Melkweg-halo identificeert Bloodhound 30-70% meer overlevende subhalo's dan de standaardmethode.
  • Voor zware subhalo's ($V_{peak} > 20$ km/s) vindt de standaardmethode in de Disc-simulaties geen enkele overlevende subhalo binnen 50 kpc, terwijl Bloodhound er gemiddeld 3 per halo vindt.
  • Binnen 300 kpc is de toename ongeveer 25% voor subhalo's met $V_{peak} \ge 10$ km/s.

• Disruptie bij Lagere Massa's:

  • De standaardmethode "verliest" subhalo's vaak bij een $V_{max}$ van ongeveer 10-20 km/s.
  • Bloodhound volgt ze tot een mediane $V_{max}$ van slechts 6 km/s (een factor 2 lager in snelheid, wat overeenkomt met een massareductie van factor 6-65).
  • Dit betekent dat subhalo's veel langer "overleven" in de simulaties dan eerder werd gedacht.

5. Betekenis en Conclusie

• Betrouwbaarheid van Donkere Materie-studies: De resultaten benadrukken dat de schijnbare afwezigheid van subhalo's in de binnenste regio's van de Melkweg deels een artefact is van onvoldoende tracking-methoden, en niet noodzakelijk een fysiek fenomeen.
• Invloed op Satellietgalaxie-tellingen: Het feit dat Bloodhound veel meer overlevende subhalo's vindt (vooral die met $V_{peak} > 17$ km/s, de drempel voor atomaire koeling), heeft grote implicaties voor de voorspelling van het aantal ultra-faint dwergsterrenstelsels rond de Melkweg.
• Rol van "Orphan" Galaxieën: De studie ondersteunt het idee dat "orphan" galaxieën (waarvan de donkere materie subhalo onder de detectiedrempel is gestript, maar die wel nog sterren bevatten) een belangrijke rol spelen in de totale populatie van satellietgalaxieën.
• Toekomstige Toepassingen: Bloodhound biedt een robuuste basis voor het interpreteren van toekomstige data van observatoria zoals LSST (Vera Rubin Observatory), waarbij nauwkeurige voorspellingen over substructuur essentieel zijn om de aard van donkere materie te testen.

Samenvattend: Bloodhound is een krachtige, deeltjesgebaseerde tool die de beperkingen van traditionele snapshot-gebaseerde methoden overwint in dichte omgevingen. Het levert een veel completer en nauwkeuriger beeld op van de substructuur van kosmische halo's, wat cruciaal is voor het begrijpen van de vorming van de Melkweg en de eigenschappen van donkere materie.