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

Technische Zusammenfassung: Die durchschnittlichen ferninfraroten Eigenschaften von Euclid-selektierten sternbildenden Galaxien

Problem und Motivation
Die erste Quick Data Release (Q1) der Euclid-Mission stellt einen beispiellosen optischen und nahen Infrarot-Datensatz (NIR) bereit, der Millionen von Galaxien mit robusten photometrischen Rotverschiebungen, Sternenmassen und Sternentstehungsraten (SFRs) enthält, die mittels SED-Fitting (Spektrale Energieverteilung) abgeleitet wurden. Diese aus optischen/NIR-Daten abgeleiteten SFRs beruhen jedoch auf Annahmen über die Staubextinktion und die zugrunde liegende Sternpopulation. Um die durchschnittlichen ferninfraroten (far-IR) Eigenschaften dieser Galaxien unabhängig zu bestimmen und die Konsistenz der SFR-Schätzungen zu verifizieren, untersucht diese Studie die thermische Emission von Staubkörnern. Während frühere Stacking-Analysen von Far-IR-Daten (z. B. in den Feldern UDS, COSMOS oder GAMA) durch kleine Himmelsflächen oder Rotverschiebungsbereiche begrenzt waren, ermöglicht der Euclid Q1-Katalog das Stacking von Millionen von Galaxien über Zehn Quadratgrad hinweg, was eine statistisch robuste Charakterisierung der Hauptsequenz (MS) von sternbildenden Galaxien bis zu $z \approx 3$ ermöglicht.

Methodik
Die Autoren verwendeten den Euclid Q1 Merging-Katalog (speziell die von Euclid Collaboration: Enia et al. 2025 beschriebene Stichprobe) und wählten sternbildende Galaxien im Rotverschiebungsbereich von $0,2 < z < 3,0$ und im Sternenmassenbereich von $8,3 < \log_{10}(M_*/M_\odot) < 11,5$ aus. Die Stichprobe wurde mit bestehenden Far-IR-Bilddaten des Herschel Weltraumobservatoriums (PACS bei 100/160 $\mu$m; SPIRE bei 250/350/500 $\mu$m) und dem SCUBA-2-Instrument (450/850 $\mu$m) abgeglichen. Die Analyse deckte eine gesamte überlappende Fläche von 17,6 deg² für SPIRE und 2,4 deg² für SCUBA-2 ab.

Um die grobe Winkelauflösung von Far-IR-Instrumenten (10–30 Bogensekunden) im Vergleich zur Euclid-Auflösung (~0,2 Bogensekunden) zu adressieren, verwendeten die Autoren eine Stacking-Analyse unter Verwendung des SimStack-Algorithmus. Diese Methode löst gleichzeitig die durchschnittlichen Flussdichten mehrerer Bins, die durch Sternmasse und Rotverschiebung definiert sind. Zu den wesentlichen methodischen Merkmalen gehörten:

• SimStack-Implementierung: Lösung eines linearen Systems $y = S_1 X_1 + \dots + S_N X_N$, um die Stacking-Amplituden ($S_b$) für $N$ Bins zu bestimmen, wobei Galaxien-Clustering und Bin-zu-Bin-Korrelationen berücksichtigt werden, ohne eine Poisson-Verteilung vorauszusetzen.
• Maskierung: Rigorose Maskierung von Regionen, die durch helle Sterne, Randeffekte und galaktische Quellen (z. B. NGC 6543) kontaminiert sind, um eine Verzerrung des gestapelten Flusses zu verhindern.
• SED-Fitting: Die gestapelten Flussdichten wurden mit einem modifizierten Schwarzkörpermodell gefittet, das bei hohen Frequenzen in ein Potenzgesetz übergeht. Der Staubemissionsindex ($\beta$) und die Potenzgesetz-Steigung ($\alpha$) wurden auf Werte fixiert, die aus lokalen far-IR-selektierten Galaxien abgeleitet wurden ($\beta=1,96, \alpha=2,3$), aufgrund begrenzter photometrischer Constraints.
• Parameterableitung: Aus den Best-Fit-SEDs leiteten die Autoren mittlere Staubtemperaturen ($T_d$), Staubmassen ($M_d$) und far-IR-abgeleitete SFRs ab.

Wichtigste Ergebnisse

1. SFR-Konsistenz: Die far-IR-abgeleiteten SFRs zeigen eine gute Übereinstimmung mit den SFRs, die durch die optische/NIR-Photometrie der sternbildenden Hauptsequenz (MS) vorhergesagt werden, mit einem gewichteten Mittelwertverhältnis von 1,0.
2. Entwicklung der Staubtemperatur: Die durchschnittliche Staubtemperatur ($T_d$) nimmt mit der kosmischen Zeit (steigende Rotverschiebung) ab. Die Entwicklung wird gut durch die Funktion $T_d(t) = T_2 + (T_1 - T_2)e^{-t/\tau}$ beschrieben, wobei $T_1 = 79,7 \pm 7,4$ K (extrapoliert auf $t=0$), $T_2 = 23,2 \pm 0,1$ K ( asymptotischer Wert) und $\tau = 1,6 \pm 0,1$ Gyr ist.
   • Entscheidend ist, dass $T_d$ über die untersuchten Bins hinweg unabhängig von der Sternmasse ist.
   • Die Konvergenz gegen einen von Null verschiedenen Wert ($T_2 \approx 23$ K) unterhalb von $z=1$ deutet darauf hin, dass der Staub in MS-Galaxien zunehmend durch die vorhandene kühlere, ältere Sternpopulation erwärmt wird und nicht ausschließlich durch heiße junge Sterne in Sternentstehungsregionen.
3. Staubmasse und Verhältnisse: Staubmassen ($M_d$) und Staub-zu-Sternmasse-Verhältnisse ($M_d/M_*$) steigen von $z=0,2$ bis $z \approx 1$ an, wonach sie ein Plateau erreichen oder abfallen. Das Staub-zu-Sternmasse-Verhältnis sinkt mit zunehmender Sternmasse bei allen Rotverschiebungen, was auf niedrigere spezifische SFRs (sSFRs) in massereicheren Galaxien zurückzuführen ist.
4. Kosmischer Infrarot-Hintergrund (CIB): Der Euclid Q1-Katalog löst mehr als 60 % des CIB bei 250, 350 und 500 $\mu$m auf. Bei 100 $\mu$m wird der gelöste Anteil auf 30–80 % geschätzt, während er bei 850 $\mu$m 30–40 % beträgt.
5. Vergleich mit Simulationen: Im Vergleich zum MAMBO-Mock-Katalog (basierend auf der Millennium-Simulation und dem Empirical Galaxy Generator) stimmen die simulierten Staubmassen und SFRs vern reasonably gut überein. Die Simulation sagt jedoch eine monotone Abnahme der Staubtemperatur mit der Zeit voraus und kann das beobachtete Plateau bei niedrigen Rotverschiebungen ($z < 1$) nicht reproduzieren.

Bedeutung und Ansprüche
Das Paper behauptet, dass diese Studie die größte Stacking-Analyse ihrer Art darstellt, die den einzigartigen Tiefen- und Flächengehalt des Euclid Q1-Releases nutzt. Die primäre Bedeutung liegt in der Entdeckung, dass die durchschnittliche Staubtemperatur von sternbildenden Galaxien bei niedrigen Rotverschiebungen gegen einen konstanten Wert konvergiert, unabhängig von der Sternmasse. Dieser Befund stellt Modelle in Frage, in denen die Staubtemperatur strikt proportional zu SFR oder sSFR ist, und deutet auf einen Übergang im Heizmechanismus des Staubs von jungen, massereichen Sternen zu dem allgemeinen interstellaren Strahlungsfeld älterer Sternpopulationen hin, während das Universum altert.

Darüber hinaus zeigen die Autoren, dass die Korrelation zwischen Staubtemperatur und SFR durch die Sternmasse vermittelt wird, da $T_d$ selbst massenunabhängig ist. Die Ergebnisse liefern zudem einen Referenzwert für zukünftige Euclid-Datensätze und deuten darauf hin, dass tiefere Bildgebungen erforderlich sein werden, um den CIB vollständig aufzulösen und die far-IR-Eigenschaften von Galaxien bei höheren Rotverschiebungen ($z > 1,5$) und geringeren Sternmassen zu untersuchen. Die Studie kommt zu dem Schluss, dass die aktuelle Q1-Daten zwar nicht den gesamten CIB auflösen, aber den Großteil des Hintergrunds bei SPIRE-Wellenlängen abdecken, was das Potenzial von Euclid bestätigt, die thermische Geschichte der sternbildenden Population des Universums abzubilden.