The Density Profile of Dynamical Halos

Technische Samenvatting: Het Dichtheidsprofiel van Dynamische Halo's

Probleemstelling
De fundamentele eigenschappen van donkere materie-halo's, met name hun dichtheidsprofielen, worden traditioneel gekenmerkt door analytische fits (bijv. NFW, Einasto) op de totale massaverdeling binnen een referentieradius zoals $R_{200m}$. Deze benadering combineert echter twee verschillende dynamische populaties: deeltjes die gravitationeel gebonden zijn en binnen de halo in een baan ronddraaien, en deeltjes die wel binnenkomen (infall) maar nog niet gebonden zijn. Recente voorstellen suggereren om een "dynamische halo" strikt te definiëren als de verzameling deeltjes die in een baan ronddraaien (orbiting). Hoewel deze definitie het mengen van interne structuur met voortdurende accretie voorkomt, waren de structurele eigenschappen (gemiddelde en spreiding) van deze orbitale profielen nog niet volledig gekarakteriseerd. Specifiek was het onduidelijk welke structurele vrijheidsgraden aanwezig blijven bij een vaste orbitale massa, wat de omvang van de intrinsieke spreiding is, en of correlaties met de vormingsgeschiedenis genuanteerde eigenschappen zijn van de orbitale component of artefacten van het mengen van populaties.

Methodologie
De auteurs maakten gebruik van een Cold Dark Matter (CDM) simulatie ($1024^3$ deeltjes in een $1 \, h^{-1}$ Gpc box) uitgevoerd met GADGET-2. Ze genereerden een catalogus van dynamische halo's met behulp van de Rockstar halo-finder als startpunt, en pasten vervolgens het algoritme van García et al. toe om deeltjes te labelen als ofwel "orbiting" (in een baan ronddraaiend) of "infalling" (binnenkomend) op basis van hun accretietijd ($a_{acc}$) en tijdgemiddelde radiale snelheden.

• Definitie Orbiting: Deeltjes met vroege accretietijden en een nabij nul liggende tijdgemiddelde radiale snelheid worden gelabeld als "orbiting"; deeltjes met late accretietijden en negatieve radiale snelheden als "infalling".
• Massabinning: Halo's werden ingedeeld in negen nauwe massabins ($\log M_{orb} \in [13.40, 15.00]$).
• Parametrisch Model: De orbitale dichtheidsprofielen werden gefit met behulp van het model van Salazar et al., gedefinieerd als:
  $$\rho_{orb}(x|M_{orb}) = A \left(\frac{x}{\epsilon}\right)^{-\alpha(x|M_{orb})} \exp\left(-\frac{x^2}{2}\right)$$
  waarbij $x \equiv r/r_h$, $r_h$ de halo-radius is (exponentiële afschermingsschaal), en de helling $\alpha(x)$ overgaat van 0 naar een asymptotische waarde $\alpha_\infty$. De parameter $\epsilon$ werd vastgesteld op 0,037.
• Fit-strategie: Individuele halo's werden gefit om de best passende $r_h$ en $\alpha_\infty$ te bepalen. De auteurs onderzochten vervolgens de correlatie tussen deze parameters en de afhankelijkheid van de accretiegeschiedenis, gekarakteriseerd door een relatieve vormings-tijdvariabele $a_{RF}$.

Belangrijkste Bijdragen en Resultaten

1. Enkele Effectieve Vrijheidsgraad:
   De studie toont aan dat bij een vaste orbitale massa de twee vrije parameters van het profiel ($r_h$ en $\alpha_\infty$) sterk gecorreleerd zijn. Meer uitgebreide halo's (grotere $r_h$) bezitten steilere binnenste hellingen (grotere $\alpha_\infty$). Deze correlatie maakt het mogelijk om het profiel van een individuele dynamische halo te karakteriseren met een enkele vrije parameter: de halo-radius $r_h$. De relatie wordt gegeven door:
   $$\Delta\alpha_\infty = 0.05 + \ln R$$
   waarbij $R = r_h / r_{h,st}(M_{orb})$ en $\Delta\alpha_\infty$ de afwijking is van de gestapelde (stacked) helling.

2. Spreiding in Halo-radius:

   • Bij een vaste orbitale massa is de distributie van halo-radii ongeveer log-normaal. De intrinsieke spreiding in $\ln r_h$ is $\sigma_{\ln r_h | M_{orb}} \approx 0.16$ (16%).
   • De auteurs onderzochten de afhankelijkheid van deze spreiding van de accretiegeschiedenis. Ze definieerden een relatieve vormingstijd $a_{RF}$ op basis van het 60e percentiel van de deeltjes-accretietijdverdeling.
   • Er bestaat een correlatie tussen $r_h$ en $a_{RF}$: laat gevormde halo's zijn compacter (kleinere $r_h$).
   • Rekening houden met deze correlatie vermindert de spreiding in de halo-radius bij een vaste massa en vormingstijd tot $\sigma_{\ln r_h | M_{orb}, a_{RF}} \approx 0.11$ (11%). De auteurs merken op dat deze reductie marginaal is, wat suggereert dat de relatieve vormingstijd alleen onvoldoende is om de spreiding volledig te verklaren.

3. Vergelijking met Accretiesnelheid:
   De auteurs vergeleken de correlatie tussen de halo-radius en de relatieve vormingstijd met de correlatie tussen de halo-radius en de accretiesnelheid per eenheid massa. Ze vonden dat de trends kwalitatief vergelijkbaar zijn, waarbij de correlatie met de accretiesnelheid iets strakkere beperkingen (constraints) vertoont maar meer uitgesproken niet-Gaussische staarten heeft. Ze waarschuwen dat ze de accretiesnelheden voor dynamische halo's in deze specifieke analyse niet zelfconsistent konden schatten.

4. Modelprestaties:
   De parametrische fits aan de orbitale profielen leveren residuen op die vergelijkbaar zijn met NFW-fits toegepast op totale massaprofielen (mediaan residu $\approx -0,02$, spreiding $\sim 0,15$). Dit valideert het gebruik van het Salazar et al.-model voor de orbitale component.

Significantie en Claims
Het artikel claimt de eerste karakterisering te bieden van de afhankelijkheid van de accretiegeschiedenis voor de orbitale profielen van dynamisch gedefinieerde halo's. Door de orbitale component te isoleren, stellen de auteurs dat de gemeten correlaties (bijv. tussen structuur en vormingstijd) werkelijke eigenschappen zijn van het gebonden systeem en geen definitionele artefacten veroorzaakt door het mengen van binnenkomend materiaal.

De primaire significantie ligt in het vermogen om de massa en de vormingstijd van een halo "af te lezen" uit de structurele eigenschappen ($r_h$ en $\alpha_\infty$). In tegenstelling tot standaard "spherical overdensity" halo's, waarbij de radius definitief verbonden is aan de omsloten massa, leggen dynamische halo's bloot hoe de omvang en structuur van een halo voortkomen uit orbitale dynamica. Deze kalibratie dient als een noodzakelijke baseline voor forward modeling, wat leidt tot nauwkeurigere dynamische decomposities van halo-dichtheidsprofielen en een betere interpretatie van de uit de lensing- en clusteringdata afgeleide massaverdelingen. De auteurs benadrukken dat hoewel hun model de empirische fits voor totale materieverdelingen niet vervangt, het essentieel is voor het begrijpen van de onderliggende dynamische structuur.