Collapse versus Disruption: The Fate of Compact Stellar Systems in Ultralight Dark Matter Halos

Technische Samenvatting: Collaps versus Disruptie in Ultralichte Donkere Materia Halo's

Probleemstelling
Modellen voor ultralicht donkere materie (ULDM) voorspellen dat de interferentie van ULDM-golven stochastische dichtheidsgranulatie genereert, wat leidt tot tijdsafhankelijke fluctuaties in het gravitationele potentiaal. Deze fluctuaties induceren dynamische verhitting in ingebedde stelsels, waardoor deze expanderen. Hoewel dit mechanisme is gebruikt om strikte ondergrenzen vast te stellen voor de deeltjesmassa van ULDM ($m_a$) met behulp van kleine systemen zoals nucleaire sterrenclusters en ultra-zwakke dwergstelsels (UFD's), blijft het evolutionaire gedrag van systemen die aanzienlijk kleiner zijn dan de de Broglie-golflengte ($\lambda_{dB}$) van ULDM slecht begrepen. Specifiek is het onduidelijk hoe de wisselwerking tussen ULDM-geïnduceerde verhitting en interne stellaire relaxatie (gedreven door twee-lichamen gravitatieverstrooiing) het lot van compacte stellaire systemen beïft. Eerdere studies suggereerden dat verhitting ineffectief zou zijn voor zeer kleine systemen, terwijl anderen opmerkten dat compacte systemen een relaxatietijd hebben die korter is dan de leeftijd van het universum, wat potentieel de kerncollaps (core collapse) dominant maakt.

Methodologie
De auteurs maken gebruik van een tweestaps numeriek simulatie-framework om de uiteenlopende ruimtelijke en temporele schalen te resolveren:

1. Fase 1 (Halo Evolutie): De ULDM-halo wordt gemodelleerd als een geheel, evoluerend binnen de getijdenvelden van de Melkweg (MW) en de Grote Magelhaense Wolk (LMC). Het stellaire systeem wordt behandeld als een puntmassa ($M_\star = 10^3 M_\odot$) om het traject en de tijdsafhankelijke getijden-tensor geïnduceerd door de ULDM-onregelmatigheden vast te leggen. De simulatie hanteert een realistische assemblagegeschiedenis waarbij de halo ongeveer 3 Gyr buiten de MW evolueert voordat deze op een traject invalt dat vergelijkbaar is met dat van de UFD Segue 1.
2. Fase 2 (Interne Evolutie): Met behulp van de directe $N$-body code PeTar wordt de interne evolutie van het stellaire systeem gesimuleerd onder de tijdsafhankelijke getijden-tensor afgeleid uit Fase 1. Het stellaire systeem wordt geïnitialiseerd met een Plummer-profiel in benaderde dynamische evenwicht.
   • Parameters: Zes representatieve systemen (S1–S6) worden gesimuleerd met variërende totale stellaire massa's ($M_\star$), deeltjesaantallen ($N$) en halve-massa-stralen ($R_h$).
   • Vergelijking: Elk systeem wordt zowel geëvolueerd met ULDM-geïnduceerde verhitting als in isolatie (zonder ULDM) om het verhittings-effect te isoleren.
   • Fysica: De ULDM-versnelling wordt benaderd via de getijden-tensor die werkt op de stellaire deeltjes relatief aan het zwaartepunt van het systeem. De simulaties draaien gedurende 13 Gyr.

Belangrijkste Bijdragen en Resultaten
De studie onthult dat de evolutie van compacte stellaire systemen in ULDM-halo's wordt beheerst door een competitie tussen interne stellaire relaxatie en ULDM-geïnduceerde verhitting, wat leidt tot drie verschillende evolutionaire regimes:

1. Relaxatie-gedomineerd Regime (Kerncollaps):

   • In voldoende compacte systemen (bijv. S1, S3, S4) drijft interne twee-lichamen relaxatie de kerncollaps aan.
   • Contraintuïtieve Bevinding: ULDM-verhitting voorkomt de collaps niet, maar versnelt deze juist. Door de buitenste sterren weg te strippen, vermindert ULDM-verhitting de gebonden stellaire massa en de halve-massa-straal. Deze reductie verkort de relaxatietijd ($t_{rh}$) en verlaagt de snelheidsdispersie, waardoor de efficiëntie van twee-lichamen relaxatie wordt vergroot. Bijgevolg blijven deze systemen gebonden en worden ze zelfs meer geconcentreerd dan ze in isolatie zouden zijn.

2. Verhitting-gedomineerd Regime (Disruptie):

   • In meer uitgebreide systemen (bijv. S2, S5) is de initiële stellaire dichtheid laag genoeg zodat de relaxatietijd lang is. Hier domineert ULDM-verhitting, waarbij energie wordt geïnjecteerd die het systeem naar een dispergerende staat drijft.
   • Voorbij een kritische drempel wordt het systeem volledig gedisrumpeerd (disruptie).

3. De "Segue 1" Analogie:

   • De auteurs identificeren een regime nabij disruptie waarvan de resten lijken op geobserveerde UFD's zoals Segue 1.
   • Simulatie S5, die rond $t \approx 10$ Gyr de disruptie nadert, produceert stellaire oppervlakte-dichtheid en snelheidsdispersie-profielen die breed overeenkomen met de observationele data voor Segue 1.
   • Een systeem met meer realistische stellaire massa's (S6) reproduceert eveneens Segue 1-achtige eigenschappen op een jongere leeftijd ($\sim 6,5$ Gyr). Dit suggereert dat voor $m_a = 10^{-22}$ eV, Segue 1 geïnterpreteerd kan worden als een systeem dat zich momenteel nabij disruptie bevindt.

Fasediagram en Kwantificering
Om deze bevindingen te generaliseren, introduceren de auteurs een dimensieloze parameter $\alpha$, gedefinieerd als de ratio van de halve-massa relaxatietijd ($t_{rh}$) tot de karakteristieke ULDM disruptie-tijdschaal ($t_{dis}$):
$$\alpha \equiv \frac{t_{rh}}{t_{dis}}$$

• $\alpha < 0,05$: Relaxatie domineert, wat leidt tot kerncollaps.
• $\alpha > 0,05$: Verhitting domineert, wat leidt tot disruptie.
• $t_{rh} > 4,7$ Gyr (en $t_{dis} > 85$ Gyr): Beide effecten zijn verwaarloosbaar over de leeftijd van het universum; systemen blijven stabiel.

Met behulp van uitgebreide simulaties construeren de auteurs een tweedimensionaal fasediagram in het $(N, M_\star R_h^{-3})$ vlak. Dit diagram scheidt de parameterruimte duidelijk in regio's van kerncollaps, disruptie en verwaarloosbare evolutie, en toont aan dat het evolutionaire lot wordt bepaald door de initiële systeemparameters ($M_\star, N, R_h$) en de deeltjesmassa van de ULDM.

Betekenis en Claims
Het artikel beweert dat de dynamica van compacte stellaire systemen in ULDM-halo's "rijk en niet-triviaal" is, wat de simplistische visie uitdaagt dat ULDM-verhitting altijd tot expansie leidt of dat kleine systemen immuun zijn voor verhittingseffecten.

• Robuuste Constraints: De auteurs betogen dat precieze modellering van de wisselwerking tussen relaxatie en verhitting essentieel is voor het afleiden van robuuste beperkingen op de deeltjesmassa van ULDM. Het negeren van de versnelling van kerncollaps of de specifieke disruptiedrempels zou kunnen leiden tot onjuiste grenzen voor de deeltjesmassa van ULDM.
• Observationele Interpretatie: De resultaten bieden een plausibel evolutionair scenario voor UFD's zoals Segue 1, wat suggereert dat zij mogelijk systemen zijn in een transiënte staat nabij disruptie in plaats van statische evenwichtsobjecten.
• Beperkingen: De auteurs geven expliciet aan dat hun analyse wordt uitgevoerd bij een vaste ULDM deeltjesmassa en halo-configuratie. Zij merken op dat het integreren van aanvullende fysica, zoals een realistische stellaire initiële massafunctie (IMF) en stellaire evolutie, het beeld verder zou compliceren en noodzakelijk is voor het vestigen van robuustere verbindingen met observaties.