Intrinsic galaxy alignments in CAMELS simulations and the significant impact of baryon model

Technische Samenvatting: Intrinsieke Galactische Uitlijningen in CAMELS-simulaties en de Significante Impact van het Baryonische Model

Probleemstelling
Intrinsieke uitlijningen (IA) van sterrenstelsels dienen zowel als een sonde voor de vorming van sterrenstelsels als een belangrijke bron van systematische fouten in zwakke gravitationele lensanalyses. De huidige kosmologische inferentie uit fotometrische surveys wordt beperkt door onzekerheden in IA-modellering, met name met betrekking tot hoe baryonische feedbackprocessen (supernovae en actieve galactische kernen - AGN) en kosmologische parameters de uitlijningssignalen beïnvloeden. Hoewel er theoretische kaders bestaan (bijv. tidal alignment voor vroege-type sterrenstelsels en tidal torquing voor late-type), blijft de specifieke afhankelijkheid van IA van een breed scala aan astrofysische en kosmologische parameters slecht geconstraint, vooral bij het onderscheiden van verschillende populaties sterrenstelsels (stervormend versus rustend/quiescent).

Methodologie
De auteurs maken gebruik van de CAMELS (Cosmology and Astrophysics with MachinE Learning Simulations) suite, specifiek de Latin-Hypercube (LH) set van 1000 hydrodynamische simulaties met behulp van de IllustrisTNG-code. Deze simulaties variëren zes sleutelparameters: kosmologische parameters ($\Omega_m$, $\sigma_8$) en astrofysische feedbackparameters ($A_{SN1}$, $A_{SN2}$ voor supernovae; $A_{AGN1}$, $A_{AGN2}$ voor AGN).

• Dataconstructie: Sterrenstelselvormen worden gemeten op $z=0$ van subhalo's geïdentificeerd door het SUBFIND-algoritme, geselecteerd op basis van stellaire massa ($M_{2r} > 10^8 M_\odot$) en aantal deeltjes ($N_s > 10$). Ellipticiteiten worden berekend uit de tweede-momententensor van de stellaire dichtheid in drie Cartesiaanse projecties, wat resulteert in 3.000 catalogi.
• Statistische Maten: De studie meet de geprojecteerde cross-correlatiefunctie tussen de ellipticiteit van sterrenstelsels en het onderliggende donkere materiedichtheidsveld, aangeduid als $w_{m+}(r_\perp)$. Een nulltest wordt uitgevoerd met behulp van de cross-component correlatie $w_{m\times}$.
• Subsamples: Sterrenstelsels worden gesplitst in "star-forming" (stervormende) en "quiescent" (rustende) populaties op basis van een specifieke sterformingssnelheid (sSFR) drempelwaarde van $10^{-10.5} \text{ yr}^{-1}$.
• Oriëntatie-alleen Analyse: Om het effect van oriëntatie te isoleren van de magnitude van de ellipticiteit, construeren de auteurs "oriëntatie-alleen" catalogi waarbij alle ellipticiteitsmagnitudes genormaliseerd worden naar de gemiddelde waarde $\langle|\varepsilon|\rangle$, waarbij een overeenkomstige correlatiefunctie $\tilde{w}_{m+}$ wordt berekend.
• Parameterafhankelijkheid: Uitlijningsamplitudes worden gedefinieerd als gewogen gemiddelden van $w_{m+}$ over transversale scheidingsbins. Pearson correlatiecoëfficiënten en p-waarden worden berekend om lineaire afhankelijkheden tussen deze amplitudes en de simulatieparameters te testen. De ratio $w_{m+}/w_{mm}$ wordt ook geanalyseerd om consistentie met lineaire uitlijningsmodellen (NLA/TATT) te testen.

Belangrijkste Resultaten

1. Detectie en Populatieverschillen: Er is een significant intrinsiek uitlijningssignaal gedetecteerd in de CAMELS-simulaties. Rustende sterrenstelsels vertonen uitlijningsamplitudes die ongeveer een orde van grootte groter zijn dan die van stervormende sterrenstelsels.
2. Parameterafhankelijkheden:
   • Kosmologie: De uitlijningsamplitude $w_{m+}$ vertoont statistisch significante correlaties met $\Omega_m$ en $\sigma_8$. Echter, wanneer het signaal genormaliseerd wordt door de materiedichtheid auto-correlatie ($w_{mm}$), verdwijnt de afhankelijkheid van $\sigma_8$, wat consistent is met lineaire uitlijningsmodelvoorspellingen.
   • Feedback: Significante correlaties worden gevonden tussen de uitlijningsamplitudes en supernova-feedbackparameters ($A_{SN1}$ en $A_{SN2}$). Er wordt geen significante afhankelijkheid van AGN-feedbackparameters waargenomen, wat de auteurs toeschrijven aan het beperkte simulatievolume ($(25 h^{-1} \text{ Mpc})^3$) dat de range van grootschalige modi beperkt.
3. Stervormende vs. Rustende Mechanismen:
   • De correlatie tussen de uitlijningsamplitude en $\sigma_8$ is onbeduidend voor stervormende sterrenstelsels maar sterk voor rustende en volledige samples. De auteurs suggereren dat de afhankelijkheid in de volledige sample wordt gedreven door het toenemende aandeel sterk uitgelijnde rustende sterrenstelsels bij hogere $\sigma_8$-waarden.
   • Supernova-feedback ($A_{SN1}$) vertoont een sterkere correlatie met de uitlijning in stervormende sterrenstelsels dan in rustende sterrenstelsels. Omgekeerd reguleert $A_{SN1}$ in rustende sterrenstelsels primair de specifieke sterformingssnelheden in plaats van de uitlijning.
4. Ellipticiteit vs. Oriëntatie: Het normaliseren van de ellipticiteitsmagnitudes (oriëntatie-alleen analyse) behoudt de gevoeligheid voor supernova-feedbackparameters. Dit geeft aan dat supernova-feedback niet alleen de magnitude van de sterrenstelsel-ellipticiteit beïnvloedt, maar ook de oriëntatie van de vormen van sterrenstelsels ten opzichte van de donkere materieverdeling.
5. Ratio-analyse: De ratio $w_{m+}/w_{mm}$ verwijdert de $\sigma_8$-afhankelijkheid maar behoudt significante correlaties met $\Omega_m$ en supernova-feedbackparameters. De auteurs merken een verrassende negatieve correlatie op met $\Omega_m$ in deze ratio op, wat de eenvoudige NLA/TATT-voorspellingen tegenspreekt en verder onderzoek rechtvaardigt.

Betekenis en Claims
Het artikel beweert de eerste directe meting te bieden van de afhankelijkheid van de intrinsieke uitlijning van sterrenstelsels van kosmologische en astrofysische parameters over een breed, realistisch bereik met behulp van een grote suite van hydrodynamische simulaties. De resultaten bieden belangrijke inzichten voor:

• Modellering: Het aantonen dat baryonische feedback (specifiek supernovae) de uitlijningssignalen aanzienlijk beïnvloedt, wat suggereert dat modellen deze processen moeten meenemen om biases in de schatting van kosmologische parameters te voorkomen.
• Sterrenstelsel Evolutie: Het benadrukken van verschillende uitlijningmechanismen die opereren binnen verschillende typen sterrenstelsels, waarbij supernova-feedback de oriëntatieveranderingen in stervormende sterrenstelsels drijft terwijl het de sterformingssnelheden in rustende sterrenstelsels reguleert.
• Toekomstige Surveys: Het bieden van noodzakelijke beperkingen voor de validatie van pipelines en simulatie-gebaseerde inferentie voor volgende generatie zwakke lensing-surveys, waarbij wordt benadrukt dat mismatches tussen aangenomen baryonische feedbackmodellen en IA-voorschriften significante biases kunnen introduceren.

De auteurs blijven bescheiden over de AGN-feedbackresultaten en stellen expliciet dat het gebrek aan gedetecteerde impact waarschijnlijk te wijten is aan het kleine simulatievolume, en zij beweren niet het specifieke mechanisme achter de onverwachte $\Omega_m$-afhankelijkheid te hebben opgelost.