Testing halo models for constraining astrophysical feedback with multi-probe modeling: I. 3D Power spectra and mass fractions

Probleemstelling

De komende Stage-IV astronomische surveys (zoals Euclid, LSST en DESI) zullen de verdeling van materie in het heelal met ongekende precisie meten. Om hieruit betrouwbare kosmologische parameters af te leiden, zijn theoretische modellen nodig die nauwkeurig zijn tot op het procentniveau. Een grote bron van onzekerheid in deze modellen is de astrofysische feedback geassocieerd met galaxievorming en -evolutie (zoals supernova's en feedback van actieve galactische kernen, AGN). Deze processen beïnvloeden de thermodynamica en verdeling van baryonen (gas en sterren) op schalen van enkele megaparsec, maar de onderliggende fysica is nog niet volledig begrepen.

Directe observaties van baryonen via verschillende methoden (thermisch en kinematisch Sunyaev-Zel'dovich-effect, Röntgenstraling, en dispersiemaat van snelle radiobursts) bieden een kans om deze feedback te beperken. Echter, er ontbreekt een unificerend analytisch kader dat deze verschillende waarnemingen (probes) simultaan en consistent kan modelleren. Bestaande modellen maken vaak verschillende aannames over de verdeling van materie en baryonen, en het is onduidelijk of ze flexibel genoeg zijn om alle relevante observables tegelijkertijd nauwkeurig te beschrijven zonder systematische fouten.

Methodologie

De auteurs testen drie bestaande halo-modellen uit de literatuur tegen elkaar en tegen hydrodynamische simulaties:

1. Mead20+: Een update van het model van Mead et al. (2020), dat feedback beïnvloedt via modificaties van de halo-concentratie en eenvoudige voorschriften voor gasverdeling.
2. Schneider19+: Een uitbreiding van de "baryonification"-methode van Schneider et al. (2019), die feedback modelleert door de dichtheidsprofielen van gas en donkere materie analytisch te vervormen, inclusief thermodynamische eigenschappen.
3. Arico24+: Een model van Aricò & Angulo (2024) dat thermodynamische grootheden (temperatuur, druk) integreert in een baryonificatie-approach, met een flexibele profielbeschrijving.

Simulatie:
De modellen worden getest tegen de Magneticum hydrodynamische simulatie (Box2b/hr), een onafhankelijke dataset die niet is gebruikt voor de oorspronkelijke kalibratie van deze modellen.

Analyse-aanpak:
De auteurs passen de modellen aan op de 3D auto- en kruis-spectra van verschillende velden:

• $P_{mm}$: Matter-matter (donkere materie).
• $P_{mp}$: Matter-druk (elektronen druk).
• $P_{ee}$: Elektronendichtheid-elektronendichtheid.

Ze analyseren drie scenario's:

1. Alleen de spectra van materie en druk ($R_{mm}$ en $R_{mp}$).
2. De spectra plus de gemeten massafracties van gebonden gas ($f_{BG}$) en sterren ($f_{star}$).
3. De spectra plus de elektronendichtheid ($R_{ee}$).

Ze vergelijken een minimale set (6 vrije parameters) met een uitgebreide set (11-15 parameters) om de flexibiliteit van de modellen te testen. De prestaties worden geëvalueerd op basis van de RMS-fout in de spectra en de nauwkeurigheid waarmee de modellen de onderliggende halo-eigenschappen (thermodynamische profielen en massafracties) kunnen herleiden.

Belangrijkste Bijdragen

• Unificerend Testkader: Dit is een van de eerste studies die halo-modellen systematisch test door ze simultaan te fitten op multiple probes (matter, druk, en dichtheid) in plaats van alleen op de materie-spectrum.
• Flexibiliteitsanalyse: De auteurs introduceren een methode om de "voorspellingsruimte" van modellen te vergelijken door de basisvectoren van parameters te projecteren. Dit toont aan in hoeverre het parameterbereik van het ene model de variaties van een ander model kan nabootsen.
• Validatie van Inferentie: Het artikel benadrukt dat het niet alleen gaat om het passen van spectra, maar om de fysische interpretatie: kunnen de modellen de werkelijke thermodynamische profielen (druk, temperatuur, dichtheid) van halos correct afleiden uit de geaggregeerde spectra?

Resultaten

1. Spectrale Aanduiding:

   • Alle drie de modellen kunnen de 3D power spectra reproduceren met een nauwkeurigheid van sub-percent tot enkele procenten, afhankelijk van de tracer-combinatie en het aantal vrije parameters.
   • Met een uitgebreide parameter-set (ca. 15 parameters) bereiken alle modellen een nauwkeurigheid van <1% voor alle drie de spectra ($R_{mm}, R_{mp}, R_{ee}$).
   • De Schneider19+ en Arico24+ modellen tonen een grotere flexibiliteit om de variaties van het Mead20+ model te reproduceren dan andersom.

2. Herstel van Halo-eigenschappen:

   • Hoewel de spectra goed worden gefit, varieert de nauwkeurigheid waarmee de modellen de onderliggende fysica (gasfracties en thermodynamische profielen) herleiden aanzienlijk.
   • Schneider19+ presteert het beste overall, met name bij het herleiden van thermodynamische profielen voor cluster-grootte halos (RMS fout < 10% wanneer massafracties worden meegenomen).
   • Mead20+ en Arico24+ tonen grotere afwijkingen (RMS fouten van 20-80%) in de afgeleide profielen, vooral voor groep-grootte halos. Dit suggereert dat hun aannames (zoals constante temperatuur voor uitgeworpen gas of een vaste relatie tussen gas- en donkere materie-concentratie) fysisch minder robuust zijn.

3. Degeneratie en Parameters:

   • Parameters zijn sterk gedegeerd wanneer slechts één spectrum wordt gebruikt. Het combineren van verschillende probes (matter, druk, dichtheid) breekt deze degeneraties effectief en beperkt de parameteronzekerheid.
   • Het toevoegen van massafracties ($f_{BG}, f_{star}$) als extra constraints verbetert de nauwkeurigheid van de afgeleide profielen aanzienlijk, vooral voor de Schneider19+ en Arico24+ modellen.

4. Overfitting Risico:

   • Het verhogen van het aantal vrije parameters verbetert niet altijd de herstelbaarheid van de fysische profielen. Soms leidt meer flexibiliteit tot "compenserend gedrag" waarbij parameters verschuiven om de spectra te passen, maar de onderliggende profielen onnauwkeurig blijven.

Significantie en Conclusie

De studie concludeert dat bestaande halo-modellen voldoende flexibel zijn om multi-probe data (matter, tSZ, kSZ, FRB) binnen een paar procent nauwkeurig te modelleren. Dit is een cruciale stap voor toekomstige Stage-IV surveys.

Echter, de fysische interpretatie blijft een uitdaging. Hoewel de modellen de geaggregeerde statistieken (spectra) goed kunnen reproduceren, zijn ze niet altijd in staat om de onderliggende thermodynamische profielen van individuele halo's nauwkeurig te voorspellen. Dit wijst op beperkingen in de huidige fysische aannames (zoals de behandeling van uitgeworpen gas en niet-thermische druk).

Voor betrouwbare kosmologische inferentie en een beter begrip van astrofysische feedback is het noodzakelijk:

• Meer multi-probe data te combineren om degeneraties te breken.
• De modellen te testen op een bredere reeks hydrodynamische simulaties.
• De focus te verleggen van puur empirische fitting naar modellen die gebaseerd zijn op eerste principes (first principles) van feedback-processen.

Dit werk vormt de basis voor vervolgstudies die deze modellen zullen testen in de "observable space" (met mock lightcones) om de impact van waarnemingssystematiek en projectie-effecten te kwantificeren.