Euclid preparation: TBD. Cosmic Dawn Survey: evolution of the galaxy stellar mass function across 0.2<z<6.5 measured over 10 square degrees

Titel en Context

Titel: Euclid voorbereiding: Evolutie van de sterrenmassa-functie van sterrenstelsels over 0,2 < z ≤ 6,5 gemeten over 10 vierkante graden.
Auteurs: Euclid Collaboration (o.a. L. Zalesky, J.R. Weaver, et al.)
Publicatie: Astronomy & Astrophysics (Pre-publicatie manuscript, 2025).
Doel: Dit werk dient als een "pre-launch" benchmark voor de Euclid-missie, gebruikmakend van de Cosmic Dawn Survey (DAWN) catalogi om de evolutie van de sterrenmassa-functie (Stellar Mass Function, SMF) te kwantificeren in een ongekend groot volume van het vroege heelal.

---

1. Het Probleem

De evolutie van de sterrenmassa-functie (SMF) is cruciaal voor het begrijpen van hoe sterrenstelsels groeien en hoe ze worden beïnvloed door processen zoals gasinstroom, fusies en feedback (van sterren en actieve galactische kernen, AGN). Echter, eerdere studies hebben te kampen met twee belangrijke beperkingen:

1. Cosmische Variatie: Diepe ruimtetelescopen (zoals HST en JWST) hebben een zeer klein gezichtsveld. Dit maakt het moeilijk om zeldzame, zeer massieve sterrenstelsels (boven de karakteristieke massa $M^*$) statistisch significant te tellen, vooral bij hoge roodverschuivingen ($z > 3$). De onzekerheid door "cosmische variatie" (verschillen in dichtheid tussen verschillende gebieden van het heelal) is hierdoor groot.
2. Compleetheid: Bestaande surveys missen vaak de diepte in het midden-infrarood (mid-IR) die nodig is om de sterrenmassa's van oude of stoffige sterrenstelsels nauwkeurig te meten, of ze missen de ruimtelijke afmeting om zeldzame objecten te vinden.

Er is een dringende behoefte aan een dataset die zowel diep is (voor nauwkeurige massametingen) als breed (om zeldzame objecten en hun omgeving te bestuderen) om de theorieën over sterrenvormingsefficiëntie en feedback in het vroege heelal te testen.

---

2. Methodologie

Data en Observaties:

• Dataset: De auteurs gebruiken de Cosmic Dawn Survey Pre-launch (DAWN PL) catalogi. Deze omvatten het Euclid Deep Field North (EDF-N) en Euclid Deep Field Fornax (EDF-F).
• Oppervlak: Een effectief oppervlak van 10,13 deg² met uniforme, diepe data.
• Bandbreedtes: Combinatie van diepe UV/optische data (van CFHT/MegaCam en Subaru/HSC) en diepe mid-infrarood data (Spitzer/IRAC [3,6µm] en [4,5µm]). De IRAC-data heeft een diepte van $m \sim 25$ mag (5$\sigma$), wat de grootste oppervlakte is die ooit tot deze diepte in het infrarood is gemapt.
• Redshift-bereik: $0,2 < z \le 6,5$, wat een tijdsbestek van ongeveer 10,2 miljard jaar bestrijkt.

Data-analyse en Selectie:

• Fotometrie: Gebruikmakend van The Farmer (gebaseerd op The Tractor) voor modelgebaseerde fotometrie in drukke velden.
• Fotometrische Roodverschuivingen (Photo-z): Berekenen met LePHARE (gekozen boven EAZY vanwege betere prestaties bij het vermijden van onfysische kleuren door de grote golflengte-gaten in de data).
• Selectiecriteria:
  • Bronnen geselecteerd op basis van een compositie van HSC $r+i+z$ beelden.
  • S/N > 3 in alle bands (met uitzondering van HSC $r$ voor $z > 3,5$).
  • Sterren verwijderd via SED-fitting en Gaia-gegevens.
  • Kwaliteitsfilters toegepast op de $\chi^2$ van de SED-aanpassing.
• Compleetheid: De sterrenmassa-limiet ($M_{lim}$) wordt empirisch bepaald met de methode van Pozzetti et al. (2010). Voor het totale monster is $M_{lim} \approx 4,31 \times 10^8 (1+z)^{1,74} M_\odot$.
• Classificatie: Stervormende vs. rustende (quiescent) sterrenstelsels worden gescheiden met de NUV-$r$, $r$-$J$ kleur-diagrammethode (Ilbert et al. 2013), geldig tot $z \approx 3$.

Statistiek en Onzekerheid:

• Volumecorrectie: Toepassing van de $1/V_{max}$ correctie (Schmidt 1968) om Malmquist-bias te corrigeren.
• Onzekerheidsanalyse: De totale onzekerheid ($\sigma_\Phi$) is de kwadratische som van:
  1. Poisson-onzekerheid (tellen van objecten).
  2. Cosmische variatie ($\sigma_{CV}$): Geschat met de UniverseMachine simulaties (Jespersen et al. 2024), wat een significant lagere onzekerheid oplevert dan eerdere werken.
  3. SED-fitting onzekerheid ($\sigma_{SED}$): Geschat via Monte Carlo simulaties van photo-z en massa.
• Eddington Bias: Correctie voor de bias die ontstaat doordat onzekerheid in massa objecten van lage massa naar hogere massa-bins "verspreidt". Dit wordt gemodelleerd met een convolutiekernel (Gaussisch + Lorentzian) en afgeleid uit de posterior-verdelingen van de massa.

Modellering:

• De intrinsieke SMF wordt afgeleid door de waargenomen SMF te fiten met een Schechter-functie (of een dubbele Schechter-functie voor $z \le 2$) die is geconvoled met de Eddington-bias kernel.
• MCMC-analyse (emcee) wordt gebruikt om de parameters ($\Phi^*$, $M^*$, $\alpha$) en hun onzekerheden te bepalen.

---

3. Belangrijkste Bijdragen

1. Ongekend Volume: De studie bestrijkt een volume van 0,62 Gpc³, wat een orde van grootte groter is dan eerdere werken (zoals COSMOS2020/W23) en een factor 20 groter dan Shuntov et al. (2024). Dit reduceert de cosmische variatie drastisch.
2. Pre-launch Benchmark: Het biedt de eerste robuuste meting van de SMF in dit volume zonder Euclid-data, dienend als referentiepunt voor toekomstige "post-launch" analyses.
3. Omgevingsafhankelijkheid: Voor het eerst wordt er bewijs gevonden voor een omgevingsafhankelijkheid van de SMF tot $z \sim 3,5$, gebaseerd op lokale dichtheidsmetingen.
4. Extreme Sterrenvormingsefficiëntie: De data dwingen tot de conclusie dat zeer massieve sterrenstelsels in het vroege heelal ($z > 3,5$) een geïntegreerde sterrenvormingsefficiëntie vereisen die veel hoger is dan de "universele piekefficiëntie" die vaak wordt aangenomen.

---

4. Resultaten

Evolutie van de SMF:

• De waargenomen SMF toont een snelle toename in het aantal sterrenstelsels met massa $M \sim 10^{10,5} M_\odot$ tussen $z \sim 6$ en $z \sim 3,25$ (factor > 100), gevolgd door een langzamere groei tot $z \sim 0,35$. Dit bevestigt het "downsizing"-effect.
• De resultaten zijn over het algemeen consistent met eerdere metingen (D17, W23), maar bieden nu vaste schattingen waar eerder alleen bovengrenzen waren, vooral voor de meest massieve systemen ($M > 10^{11} M_\odot$).

Stervormende vs. Rustende Populaties:

• Rustende sterrenstelsels: Een aanzienlijk deel van de meest massieve rustende sterrenstelsels ($M > 10^{11,7} M_\odot$) is al aanwezig bij $z \sim 2,5$. Hun dichtheid verandert weinig tot $z \sim 0,2$, wat suggereert dat ze weinig massa hebben opgebouwd via fusies na hun vorming.
• Beperkingen: Vanwege de optische selectie (ontbreken van NIR-data) zijn rustende sterrenstelsels bij $z > 3$ ondervertegenwoordigd; de steekproef is voornamelijk blauw en stervormend.

Sterrenvormingsefficiëntie en Feedback:

• Door de intrinsieke SMF te vergelijken met de donkere-materie-halo-massa-functie (HMF), concluderen de auteurs dat massieve sterrenstelsels bij $z > 3,5$ een geïntegreerde sterrenvormingsefficiëntie ($\epsilon$) hebben van 0,25 tot 0,5 (en soms hoger).
• Dit is veel hoger dan de gebruikelijke piekefficiëntie van $\sim 0,2$ bij lagere roodverschuivingen.
• Implicatie: Feedback-mechanismen (zoals AGN en sterrenwinden) die normaal gesproken de vorming van sterren in massieve halo's onderdrukken, zijn in het vroege heelal ($z > 3,5$) ineffectief. Dit staat in lijn met recente simulaties (FIRE-project) die aangeven dat bij hoge dichtheden gas sneller wordt omgezet in sterren dan dat het kan worden uitgestoten.

Omgevingsafhankelijkheid:

• Er is een correlatie gevonden tussen sterrenmassa en lokale dichtheid tot $z \sim 3,5$.
• Massieve sterrenstelsels zijn oververtegenwoordigd in dichte omgevingen, terwijl minder massieve sterrenstelsels vaker in lage-dichtheid omgevingen voorkomen.
• Bij $z \sim 1$ is er duidelijk bewijs voor "omgevingsdemping" (environmental quenching) van rustende sterrenstelsels in dichte gebieden.

Ultra-massieve Kandidaten:

• De studie identificeert enkele ultra-massieve kandidaten, waaronder een sterrenstelsel bij $z \approx 5,15$ met een geschatte massa van $10^{12,5} M_\odot$. Hoewel de onzekerheid groot is (door stof en ontbrekende NIR-data), suggereert de aanwezigheid van deze objecten dat extreme evolutiepaden mogelijk zijn in het vroege heelal.

---

5. Betekenis en Conclusie

Deze studie markeert een mijlpaal in de studie van de evolutie van sterrenstelsels door de combinatie van diepte en breedte van de DAWN PL data.

• Validatie van theorieën: De bevindingen van hoge sterrenvormingsefficiënties bij hoge $z$ ondersteunen het idee dat de "koele" feedback-mechanismen die het vroege heelal domineren, nog niet effectief genoeg zijn om de vorming van supermassieve sterrenstelsels te stoppen.
• Voorbereiding op Euclid: De resultaten dienen als een cruciale benchmark. De toekomstige Euclid-data (met diepe NIR-fotometrie en spectroscopie) zal de onzekerheden in de sterrenmassa's en roodverschuivingen verder verkleinen, de detectie van rustende sterrenstelsels bij $z > 3$ mogelijk maken en de statistische kracht vergroten voor het bestuderen van de "Cosmic Dawn".
• Toekomstperspectief: De combinatie van dit werk met toekomstige Euclid-observaties zal het mogelijk maken om de evolutie van de SMF door de periode van re-ionisatie ($z \sim 8$) te volgen en de rol van de lokale omgeving in de vorming van de zwaarste structuren in het heelal definitief te ontrafelen.

Kortom, dit werk levert de eerste robuuste, volumegedreven beperkingen op de overvloed van de zwaarste sterrenstelsels in het vroege heelal, wat fundamentele inzichten biedt in de fysica van sterrenvorming en feedback voordat de Euclid-missie officieel van start gaat.