Euclid Quick Data Release (Q1). The average far-infrared properties of Euclid-selected star-forming galaxies

Technische Samenvatting: De gemiddelde ver-infrarood eigenschappen van door Euclid-geselecteerde stervormende sterrenstelsels

Probleem en Motivatie
De eerste Quick Data Release (Q1) van de Euclid-missie biedt een ongekende optische en nabij-infrarood (NIR) dataset, die miljoenen sterrenstelsels bevat met robuuste fotometrische roodverschuivingen, stellaire massa's en stervormingssnelheden (SFR's) afgeleid van spectrale energieverdeling (SED) fitting. Echter, deze op optische/NIR gebaseerde SFR's rusten op aannames over stofabsorptie en de onderliggende stellaire populatie. Om de gemiddelde ver-infrarood (far-IR) eigenschappen van deze sterrenstelsels onafhankelijk te beperken en de consistentie van SFR-schattingen te verifiëren, onderzoekt deze studie de thermische emissie van stofkorrels. Hoewel eerdere stacking-analyses van far-ir data (bijv. in de UDS, COSMOS of GAMA velden) beperkt werden door kleine hemelgebieden of roodverschuivingsbereiken, maakt de Euclid Q1-catalogus het mogelijk om miljoenen sterrenstelsels te stacken over tientallen vierkante graden, wat een statistisch robuuste karakterisering van de hoofdsequentie (MS) van stervormende sterrenstelsels tot $z \approx 3$ mogelijk maakt.

Methodologie
De auteurs maakten gebruik van de Euclid Q1 merging-catalogus (specifiek de sample beschreven door Euclid Collaboration: Enia et al. 2025), waarbij stervormende sterrenstelsels werden geselecteerd binnen het roodverschuivingsbereik $0.2 < z < 3.0$ en het stellaire massabereik $8.3 < \log_{10}(M_*/M_\odot) < 11.5$. De sample werd gekruislinkt met bestaande far-ir beeldvorming van de Herschel Space Observatory (PACS bij 100/160 $\mu$m; SPIRE bij 250/350/500 $\mu$m) en het SCUBA-2 instrument (450/850 $\mu$m). De analyse besloeg een totaal overlappend gebied van 17,6 deg$^2$ voor SPIRE en 2,4 deg$^2$ voor SCUBA-2.

Om de grove hoekresolutie van far-ir instrumenten (10–30 boogseconden) vergeleken met de Euclid-resolutie (~0,2 boogseconden) aan te pakken, gebruikten de auteurs een stacking-analyse met het SimStack algoritme. Deze methode lost simultaan de gemiddelde fluxdichtheden op van meerdere bins gedefinieerd door stellaire massa en roodverschuiving:

• SimStack Implementatie: Het oplossen van een lineair systeem $y = S_1 X_1 + \dots + S_N X_N$ om de stacking-amplitudes ($S_b$) voor $N$ bins te bepalen, waarbij rekening wordt gehouden met galaxy clustering en bin-tot-bin correlaties zonder een Poisson-verdeling aan te nemen.
• Masking: Rigoureuze masking van regio's die gecontamineerd zijn door heldere sterren, randeffecten en Galactische bronnen (bijv. NGC 6543) om bias in de gestackte flux te voorkomen.
• SED Fitting: De gestackte fluxdichtheden werden gefit met een gemodificeerd zwart lichaam-model dat overgaat in een machtswet bij hoge frequenties. De stof-emissiviteitsindex ($\beta$) en de machtswet-exponent ($\alpha$) werden vastgezet op waarden afgeleid van lokale far-ir-geselecteerde sterrenstelsels ($\beta=1.96, \alpha=2.3$) vanwege beperkte fotometrische restricties.
• Parameter Afleiding: Uit de best-fit SED's leidden de auteurs de gemiddelde stoftemperaturen ($T_d$), stofmassa's ($M_d$) en far-ir-afgeleide SFR's af.

Belangrijkste Resultaten

1. SFR Consistentie: De far-ir-afgeleide SFR's vertonen een goede overeenkomst met de SFR's die worden voorspeld door de stervormende hoofdsequentie (MS) afgeleid van optische/NIR fotometrie, met een gewogen gemiddelde ratio van 1,0.
2. Stoftemperatuur Evolutie: De gemiddelde stoftemperatuur ($T_d$) neemt af met de kosmische tijd (toenemende roodverschuiving). De evolutie wordt goed beschreven door de functie $T_d(t) = T_2 + (T_1 - T_2)e^{-t/\tau}$, waarbij $T_1 = 79,7 \pm 7,4$ K (geëxtrapoleerd naar $t=0$), $T_2 = 23,2 \pm 0,1$ K (asymptotische waarde), en $\tau = 1,6 \pm 0,1$ Gyr.
   • Cruciaal is dat $T_d$ wordt gevonden als onafhankelijk van de stellaire massa over de bestudeerde bins.
   • De convergentie naar een niet-nul waarde ($T_2 \approx 23$ K) onder $z=1$ suggereert dat stof in MS-sterrenstelsels steeds meer wordt verhit door de aanwezige koelere, oudere stellaire populatie in plaats van uitsluitend door hete jonge sterren in stervormingsgebieden.
3. Stofmassa en Ratio's: Stofmassa's ($M_d$) en stof-tot-stellaire massaverhoudingen ($M_d/M_*$) nemen toe van $z=0.2$ naar $z \approx 1$, waarna ze een plateau bereiken of dalen. De stof-tot-stellaire massaverhouding neemt af met toenemende stellaire massa bij alle roodverschuivingen, wat wordt toegeschreven aan lagere specifieke SFR's (sSFR's) in massievere sterrenstelsels.
4. Kosmische Infrarood Achtergrond (CIB): De Euclid Q1-catalogus lost meer dan 60% van de CIB op bij 250, 350 en 500 $\mu$m. Bij 100 $\mu$m wordt het opgeloste fractie geschat op 30–80%, terwijl het bij 850 $\mu$m 30–40% is.
5. Vergelijking met Simulaties: Wanneer vergeleken met de MAMBO mock-catalogus (gebaseerd op de Millennium Simulatie en de Empirical Galaxy Generator), komen de gesimuleerde stofmassa's en SFR's redelijk goed overeen. De simulatie voorspelt echter een monotone afname van de stoftemperatuur met de tijd, en slaagt er niet in om het geobserveerde plateau bij lage roodverschuivingen ($z < 1$) te reproduceren.

Betekenis en Claims
De paper claimt dat deze studie de grootste stacking-analyse van dit type tot nu toe vertegenwoordigt, waarbij gebruik wordt gemaakt van de unieke diepte en het gebied van de Euclid Q1-release. De primaire betekenis ligt in de ontdekking dat de gemiddelde stoftemperatuur van stervormende sterrenstelsels convergeert naar een constante waarde bij lage roodverschuivingen, onafhankelijk van de stellaire massa. Deze bevinding daagt modellen uit waarbij de stoftemperatuur strikt proportioneel is aan de SFR of sSFR, en suggereert een transitie in het verwarmingsmechanisme van stof van jonge, massieve sterren naar de algemene interstellaire stralingsvelden van oudere stellaire populaties naarmate het universum ouder wordt.

Verder demonstreren de auteurs dat de correlatie tussen stoftemperatuur en SFR wordt gemedieerd door stellaire massa, aangezien $T_d$ zelf mass-onafhankelijk is. De resultaten bieden ook een benchmark voor toekomstige Euclid-datareleases, wat aangeeft dat diepere beeldvorming vereist zal zijn om de volledige CIB volledig op te lossen en de far-ir eigenschappen van sterrenstelsels bij hogere roodverschuivingen ($z > 1.5$) en lagere stellaire massa's te onderzoeken. De studie concludeert dat hoewel de huidige Q1-data de volledige CIB niet volledig oplost, het de meerderheid van de achtergrond bij SPIRE-golflengten account, wat het potentieel van Euclid valideert om de thermische geschiedenis van de stervormende populatie van het universum in kaart te brengen.