Cluster Infall for Mass Calibration in the Stage-IV Era

Titel: Cluster Infall voor Massakalibratie in het Stage-IV Tijdperk

Auteurs: Connor Sweeney en Eduardo Rozo (Universiteit van Arizona)
Datum: 25 maart 2026

---

1. Probleemstelling en Motivatie

Sterrenstelselclusters zijn de zwaarste gravitationeel gebonden structuren in het heelal en vormen krachtige probes voor de kosmologie, vooral via hun massa-abundantie als functie van massa. Het gebruik van deze gevoeligheid wordt echter bemoeilijkt door niet-lineaire processen binnen de clusters, zoals feedback van actieve galactische kernen en andere baryonische effecten. Deze effecten introduceren onzekerheden in de massabepaling, wat de nauwkeurigheid van kosmologische constraints beperkt.

Aan de andere kant van het spectrum zijn dynamische studies op grote stralen ($r \gtrsim 2 h^{-1}$ Mpc) relatief ongevoelig voor deze baryonische effecten, maar wel gevoelig voor het clusterpotentieel en de expansie van het heelal. Eerdere studies hebben de line-of-sight (LOS) snelheidsdispersie gebruikt als een "standaardliniaal", maar deze benadering maakt vaak onderscheid tussen verschillende populaties (bijv. "orbiting" vs. "infalling") en gebruikt voornamelijk dispersie in plaats van de volledige snelheidsverdeling.

Het doel van dit werk is om de nuttigheid van de volledige LOS-snelheidsverdeling van in vallende (infalling) sterrenstelsels op de buitenranden van clusters (beyond $r \ge 5 h^{-1}$ Mpc) te karakteriseren als een hulpmiddel voor het kalibreren van clustermassa's met een hoge precisie, robuust tegenover baryonische fysica.

2. Methodologie

De auteurs hanteren een benadering die gebaseerd is op numerieke simulaties en een geavanceerd statistisch model.

• Data en Simulaties:

  • Gebruik van de Multi-Dark Planck 2 (MDPL2) simulatie ($3840^3$ donkere materiedeeltjes).
  • Sterrenstelsels zijn gepopuleerd in de halos met het UniverseMachine-algoritme.
  • Selectie van een snapshot bij $z \approx 0.194$ (representatief voor SDSS en DESI).
  • Focus op halos met $M_{vir} \ge 8 \times 10^{13} h^{-1} M_\odot$ en sterrenstelsels met $M_* \ge 10^{10} h^{-1} M_\odot$.
  • Analysegebied: $5 h^{-1}$ Mpc tot $30 h^{-1}$ Mpc van het clustercentrum, om orbitale sterrenstelsels te vermijden.

• Model voor Snelheidsverdeling:
  De auteurs modelleren de gezamenlijke verdeling van radiale ($v_r$) en tangentiële ($v_t$) snelheden, $P(v_r, v_t | r, M)$, als een product van een gemarginaliseerde verdeling en een conditionele verdeling:
  $$P(v_r, v_t | r, M) = P(v_r | r, M) \cdot P(v_t | v_r, r, M)$$

  1. Radiale Snelheid ($v_r$):

     • De verdeling is niet symmetrisch (scheef en met lange staarten).
     • Gepatenteerd met de Johnson's SU (JSU) verdeling.
     • In plaats van vier vrije parameters, worden de parameters gekoppeld aan twee interpreteerbare parameters: de pieksnelheid ($v_{r,peak}$) en de variantie ($\sigma^2_{vr}$).
     • De afhankelijkheid van straal ($r$) en massa ($M$) wordt gemodelleerd via machtswetten en lineaire correcties.

  2. Tangentiële Snelheid ($v_t$):

     • De conditionele verdeling $P(v_t | v_r, r)$ wordt gemodelleerd als een Student's t-verdeling (met 5 vrijheidsgraden).
     • De schaalparameter $\sigma_{vt}$ hangt af van $v_r$, $r$ en $M$, en wordt geparametriseerd met een kubische functie van $v_r$.

• Projectie naar Line-of-Sight (LOS):

  • De 3D-snelheidsverdeling wordt geprojecteerd naar de LOS-snelheid ($v_{LOS}$) door integratie over de LOS-afstand, gewogen door het dichtheidsprofiel van de in vallende sterrenstelsels ($\rho_{inf}$).
  • Het dichtheidsprofiel wordt afgeleid uit het halo-galaxy correlatiefunctie-model van Salazar et al. [21], specifiek het in vallende component, om contaminatie door orbitale sterrenstelsels te voorkomen.

• Forecasting:

  • De auteurs voorspellen de nauwkeurigheid van massabepalingen met data van het Dark Energy Spectroscopic Instrument (DESI).
  • Ze gebruiken de Fisher-matrix-benadering om de onzekerheid op de herwonnen massa te berekenen, gebaseerd op het aantal waargenomen sterrenstelsel-halo paren.

3. Belangrijkste Bijdragen

1. Geavanceerd Snelheidsmodel: Een nieuwe, robuuste parametrisatie van de gezamenlijke snelheidsverdeling van in vallende materie op grote schaal, die de complexiteit van de verdeling (scheefheid, kurtosis) nauwkeurig vastlegt met een beperkt aantal parameters.
2. Massa-afhankelijkheid: Het expliciet modelleren van de massa-afhankelijkheid van de snelheidsverdelingen, wat essentieel is voor het kalibreren van clusters over een breed massabereik.
3. Robuustheid tegen Baryonische Effecten: Het bewijs dat het gebruik van de buitenste regio's ($r > 5 h^{-1}$ Mpc) van clusters, waar de dynamiek wordt gedomineerd door gravitatie en expansie, leidt tot een methode die minder gevoelig is voor de onzekerheden van baryonische fysica dan traditionele methoden.
4. DESI Forecast: Een kwantitatieve forecast die aantoont dat spectroscopische data van DESI clustermassa's kan kalibreren met een precisie die vergelijkbaar is met of beter is dan die van Stage-IV weak lensing surveys.

4. Resultaten

• Modelnauwkeurigheid: Het ontwikkelde model past de simulatiedata uitstekend aan. De voorspelde LOS-snelheidsverdelingen komen overeen met de simulatieresultaten binnen een nauwkeurigheid van ongeveer 5%.
• Massakalibratie:
  • Voor DESI (Bright Galaxy Sample en LRG) wordt een onzekerheid op de clustermassa voorspeld die sub-percent is (minder dan 1%) voor bepaalde massabins en rodeverschuivingen.
  • Bij hoge rodeverschuivingen ($z \approx 0.4 - 0.8$) presteert deze methode beter dan gestapelde weak lensing (WL) metingen.
  • Bij lagere massabins en hoge rodeverschuivingen wordt een nauwkeurigheid bereikt die aanzienlijk beter is dan die van WL.
• Vergelijking met Weak Lensing: Traditioneel wordt weak lensing gezien als de "gouden standaard" voor massametingen. Deze studie toont aan dat de kinematische methode via in vallende sterrenstelsels een evenwaardig of superieur alternatief biedt, vooral omdat deze minder afhankelijk is van systematische fouten door baryonische feedback.

5. Betekenis en Conclusie

Deze studie opent een nieuwe weg voor kosmologische analyses met sterrenstelselclusters. Door de kinematiek van in vallende sterrenstelsels op de buitenranden van clusters te benutten, kunnen onderzoekers clustermassa's kalibreren met een precisie die nodig is voor de volgende generatie kosmologische experimenten (Stage-IV).

De belangrijkste implicaties zijn:

• Onafhankelijkheid: Het biedt een onafhankelijke methode voor massakalibratie die complementair is aan weak lensing, wat cruciaal is voor het doorbreken van degeneraties in kosmologische parameters.
• Robuustheid: De methode is inherent robuuster tegenover de complexe baryonische fysica die binnenin clusters plaatsvindt.
• Toekomstperspectief: Met de komst van grote spectroscopische surveys zoals DESI, wordt deze techniek direct toepasbaar en belooft deze de precisie van cluster-abundantie metingen aanzienlijk te verbeteren, wat leidt tot nauwkeurigere beperkingen op de kosmologie en de aard van donkere energie.

De auteurs concluderen dat, hoewel de forecast perfect modelkalibratie veronderstelt en kosmologische afhankelijkheid negeert, de resultaten sterke motivatie bieden voor het integreren van "infall-regio" analyse in toekomstige cluster-kosmologische studies.