Measurement of Substructure from the Kinematics of the GD-1 Stellar Stream

Technische Samenvatting: Meting van Substructuur vanuit de Kinematica van de GD-1 Sterrenstroom

Probleemstelling
In het Cold Dark Matter (CDM) paradigma wordt voorspeld dat sterrenstelsels een enorme populatie donkere materie-subhalo's herbergen, variërend van galactische schalen tot aardmassa-schalen. Subhalo's onder de drempelwaarde voor stelselvorming ($\approx 10^8 M_\odot$) worden echter verwacht vrij te zijn van baryonische materie, waardoor ze onzichtbaar zijn voor directe observatie. Bijgevolg vereist het testen van het CDM-paradigma op kleine schalen indirecte gravitationele tracers. Sterrenstromen, die kinematisch koude structuren zijn gevormd door de getijdenverstoring van bolstaande sterrenhopen, zijn gevoelig voor verstoringen door laagmassieve subhalo's. Hoewel eerdere studies gebruik hebben gemaakt van de morfologische gaten en uitsteeksels (spurs) in stromen zoals GD-1 om subhalo-interacties af te leiden, richt dit werk zich op de noodzaak om de intrinsieke snelheidsdispersie van de stroom te beperken. Een gemeten snelheidsdispersie die significant hoger is dan die voorspeld door gladde Melkweg-potentialen, zou wijzen op verhitting gedreven door donkere materie-substructuur.

Methodologie
De auteurs analyseren de GD-1 sterrenstroom met behulp van een combinatie van observationele data en perturbatietheorie-modellering:

1. Dataverzameling en Meting van Snelheidsdispersie:

   • De studie maakt gebruik van 160 unieke lidsterren geïdentificeerd uit de datagedreven GD-1 catalogus (Starkman et al. 2025) met een lidmaatschapskans $>75\%$.
   • Radiale snelheden ($v_r$) zijn samengesteld uit vier spectrografische surveys: DESI (early data release), SDSS DR9, LAMOST DR8 en MMT Hectochelle.
   • Om de contaminatie door dubbelsterren te beperken, gebruiken de auteurs herhaalde observaties van 34 sterren met verschillende instrumenten.
   • Een hiërarchisch Bayesiaans model wordt toegepast om observationele fouten, door binaire systemen veroorzaakte snelheidsschommelingen ($\sigma_{binary}$) en de intrinsieke snelheidsdispersie van de stroom ($\sigma_{\Delta v_r}$) te scheiden. De stroom wordt in drie segmenten langs zijn lengte in segmenten verdeeld ($\phi_1 \in [-60, -40], [-40, -20], [-20, 0]$ graden) om de lokale dispersie te meten.

2. Dynamische Modellering:

   • De auteurs maken gebruik van het streamsculptor framework, dat perturbatietheorie gebruikt om het cumulatieve effect van honderden subhalo-impacten op de stroom te modelleren.
   • Onverstoorde stroommodellen worden gegenereerd met een particle-spray benadering voor progenitor-massa's van $2 \times 10^4 M_\odot$ en $10^5 M_\odot$ en dynamische leeftijden van 3,5, 5 en 8 Gyr, binnen een statisch Melkweg-potentiaal (MilkyWayPotential2022).
   • Subhalo's worden gemodelleerd als Hernquist-profielen. De studie verkent twee scenario's om degeneraties tussen de subhalo-getalsdichtheid en interne structuur te doorbreken:
     • Model I: Gaat uit van een Warm Dark Matter (WDM) scenario met een onderdrukking in de subhalo-massa-functie (SHMF) onder een half-mode massa ($M_{hm}$) en een constante massa-grootte relatie ($r_s/r_{s,cdm}$).
     • Model II: Gaat uit van een standaard CDM SHMF ($M_{hm}=0$) maar staat toe dat de massa-grootte relatie varieert als een gebroken machtswet met de massa, waarbij wordt getest op massa-afhankelijke concentratie.

3. Inference (Afleiding):

   • De auteurs gebruiken Approximate Bayesian Computation (ABC) om gesimuleerde stroomrealisaties te vergelijken met de geobserveerde snelheidsdispersie.
   • Simulaties worden gegenereerd voor diverse subhalo-populatieparameters (normalisatie, $M_{hm}$, concentratie). Monsters worden geaccepteerd als de Euclidische afstand tussen de gesimuleerde en geobserveerde snelheidsdispersie (en stroombreedte-beperkingen) binnen een tolerantie $\epsilon$ valt.

Belangrijkste Bijdragen en Resultaten

• Intrinsieke Snelheidsdispersie: De studie meet de intrinsieke radiale snelheidsdispersie van GD-1 op ongeveer $2\text{--}5$ km s$^{-1}$ over de gehele lengte. Het gebied met de hoogste dispersie ($\phi_1 \in [-40, -20]$ graden) vertoont een waarde van $4,8^{+1,2}_{-0,7}$ km s$^{-1}$. Dit vertegenwoordigt een $4\sigma$ afwijking van onverstoorde stroommodellen gevormd in een glad potentiaal, die dispersies van slechts $\approx 0,7\text{--}1,1$ km s$^{-1}$ voorspellen.
• Beperkingen op de Subhalo-populatie: De geobserveerde dispersie kan worden verklaard door een populatie van donkere materie-subhalo's. De afgeleide fractie van de massa van de Melkweg in subhalo's is $f_{sub} = 0,05^{+0,08}_{-0,03}$ (68% betrouwbaarheid), wat consistent is met CDM-verwachtingen ($\approx 7,5\%$).
• Compactheid van Subhalo's: Een primaire bevinding is de voorkeur voor subhalo's die compacter zijn dan de standaard CDM-voorspellingen.
  • In Model I, als de subhalo-massa fractie $10\%$ is, wordt de schaalradius beperkt tot $r_s/r_{s,cdm} = 0,3^{+0,3}_{-0,2}$ (68% betrouwbaarheid).
  • In Model II is er een voorkeur voor steeds compactere subhalo's bij lage massa's ($M < 10^7 M_\odot$). Voor subhalo's met $M_{200} = 10^6 M_\odot$, wordt de concentratie $c_{-2}$ beperkt tot $>190$ bij 68% betrouwbaarheid, wat aanzienlijk hoger is dan de CDM-voorspelling van $\approx 43$.
• Alternatieve Scenario's: De data kunnen worden verklaard door ofwel talrijke impacts met laagmassieve subhalo's, ofwel door een enkele impact met een zeer compacte subhalo ($M \gtrsim 10^8 M_\odot$). De auteurs merken op dat een enkele impact met een compacte subhalo voldoende is om de geobserveerde dispersie te generen, wat consistent is met eerdere interpretaties van het GD-1 "spur"-kenmerk.

Betekenis en Claims
Het artikel beweert een nieuwe, onafhankelijke beperking te bieden op de donkere materie-substructuur met behulp van de kinematica van sterrenstromen, wat complementair is aan eerdere beperkingen afgeleid van de dichtheidsmorfologie van stromen. De resultaten suggereren een mogelijke afwijking van de standaard CDM massa-grootte relatie bij lage subhalo-massa's, waar subhalo's naar verwachting puur door donkere materie worden gedomineerd. Met name de voorkeur voor hogere concentraties (of kleinere schaalradii) dan de CDM-voorspellingen komt overeen met alternatieve donkere materie-theorieën zoals Self-Interacting Dark Matter (SIDM) of Atomic Dark Matter, die hogere centrale dichtheden voorspellen.

De auteurs benadrukken dat hun beperkingen empirisch en agnostisch zijn ten opzien van specifieke deeltjesfysica-modellen, en dienen als een benchmark om verschillende donkere materie-theorieën te testen. Zij erkennen dat er degeneraties bestaan tussen het aantal subhalo-impacten en hun concentratie (de "rate-concentration degeneracy"), die in toekomstig werk gezamenlijk gemodelleerd kunnen worden door zowel de stroomdichtheid als de kinematica. De studie concludeert dat hoewel de subhalo-massa fractie consistent is met CDM, de interne structuur van deze subhalo's mogelijk fysica vereist die verder gaat dan het standaard CDM-paradigma, met name bij lage massa's.