DESI DR2 Reference Mocks: Clustering results from UCHUU ELGs and QSOs

Technische Samenvatting: DESI DR2 Referentiemocks van Uchuu ELGs en QSOs

Probleem en Motivatie
Het Dark Energy Spectroscopic Instrument (DESI) heeft als doel de materieverdeling van het universum in kaart te brengen met behulp van meer dan 50 miljoen spectra van sterrenstelsels en quasars om de expansiegeschiedenis en de groei van structuren te beperken. Om deze data volledig te kunnen benutten, met name voor de aanstaande DESI Data Release 2 (DR2) kosmologische analyses, zijn realistische mock-catalogi vereist. Deze mocks moeten de selectiefuncties en statistische kenmerken van DESI emuleren om pipelines te testen, systematische onzekerheden te beoordelen en theoretische voorspellingen te valideren.

Een specifieke uitdaging die in dit werk wordt aangepakt, is de modellering van emissielijnstelsels (ELGs) en quasars (QSOs). In tegenstelling tot Luminous Red Galaxies (LRGs), die naar verwachting compleet zijn in termen van halo-massa-occupatie, zijn ELGs en QSOs incomplete tracers in bekende parameterruimtes (zoals luminositeit of stellaire massa). Daarom zijn traditionele Subhalo Abundance Matching (SHAM)-technieken, die vertrouwen op een één-op-één koppeling tussen halo-eigenschappen en observeerbare grootheden van sterrenstelsels, ontoereikend. De auteurs streven ernaar om referentiemocks te construeren die de clusteringstatistieken van DESI DR1 en DR2 ELGs en QSOs reproduceren over roodverschuivingen van $0.8 < z < 2.1$ (en uitlopend tot $z=3.5$ voor QSOs) om de constructie van gesimuleerde lichtkegels te verbeteren en het begrip van de sterrenstelsel-halo-verbinding te vergroten.

Methodologie
De studie maakt gebruik van de Uchuu N-body simulatie, een hoog-resolutie, groot-volume ($2 h^{-1}$Gpc box) simulatie die $12.800^3$ donkere materiedeeltjes bevat, uitgevoerd binnen een vlakke $\Lambda$CDM-kosmologie gebaseerd op Planck 15-parameters. De simulatie biedt snapshots van $z=14$ tot $z=0$, waarbij donkere materie-halo's en subhalo's zijn geïdentificeerd met de Rockstar halo-finder.

Om deze halo's te populen met tracers, gebruiken de auteurs een gemodificeerd Subhalo Abundance Matching (SHAM) kader:

1. QSO Modellering:

   • QSOs worden aan halo's toegewezen op basis van de piek circulaire snelheid ($V_{peak}$), gedefinieerd als de maximale $V_{max}$ gemeten over alle roodverschuivingsoutputs.
   • Het model gaat uit van een Gaussische verdeling van $V_{peak}$ voor zowel centrale als satelliet-QSOs, gekenmerkt door een gemiddelde ($V_{mean}$) en een standaarddeviatie ($\sigma_V$).
   • Een vrije parameter, de satellietfractie ($f_{sat}$), controleert de relatieve normalisatie van de centrale en satellietverdelingen.
   • Om de grote onzekerheden in roodverschuiving die inherent zijn aan QSO-spectroscopie (door brede emissielijnen) te compenseren, worden Gaussische roodverschuivingsfouten ($\sigma \approx 500$ km/s) toegepast op de mock-catalogi.
   • Parameters worden afgestemd door de $\chi^2$-statistiek te minimaliseren tussen de mock-monopool van de twee-punt correlatiefunctie (2PCF) en de DESI-data over schalen van $3$ tot $70 h^{-1}$Mpc.

2. ELG Modellering:

   • In het besef dat stervorming in satellietstelsels kan worden onderdrukt door ramdruk (ram pressure stripping) als ze te snel bewegen ten opzichte van hun centrale halo, introduceren de auteurs specifieke selectiecriteria voor ELG-satellieten.
   • Satellieten worden alleen geselecteerd als ze voldoen aan: (i) $V_{max}/V_{peak} > 0.9$ (om majeure fusies te vermijden), en (ii) lage relatieve snelheden ten opzichte van de centrale halo ($V_{xyz,pec} - V_{central} < 300$ km/s) en een lage centrale $V_{peak}$ ($< 300$ km/s).
   • Net als bij de QSOs gebruikt het model een Gaussische $V_{peak}$-verdeling met parameters $V_{mean}$, $\sigma_V$ en $f_{sat}$, maar de selectiecriteria worden toegepast vóór de abundance matching.

3. Lichtkegel Constructie:

   • Mock-lichtkegels worden geconstrueerd door de Uchuu-simulatieboxen te repliceren en ze in sferische schillen te snijden om het DESI-surveyvolume te dekken.
   • De mocks worden gemaskeerd om overeen te komen met de DESI DR2-voetafdruk en gewogen om de geobserveerde comoving getalsdichtheid $n(z)$ te reproduceren in roodverschuivingsbins van $\Delta z = 0.01$.
   • De best passende parameters worden onafhankelijk bepaald voor verschillende roodverschuivingsbins en vervolgens geparametriseerd als functies van de roodverschuiving (bijv. lineaire evolutie voor de QSO $f_{sat}$ en kwadratische evolutie voor de ELG $V_{mean}$).

Belangrijkste Bijdragen

• Referentiemocks: Het artikel presenteert de eerste set referentiemocks voor DESI DR2 die specifiek zijn afgestemd op het reproduceren van de clustering van ELGs en QSOs met behulp van de Uchuu-simulatie.
• Gemodificeerde SHAM Implementatie: De auteurs hebben de SHAM-techniek succesvol aangepast voor incomplete tracers. Voor ELGs introduceren zij een op snelheid gebaseerd selectiecriterium om de omgevingsafhankelijke quenching van stervorming in satellieten te verklaren, een modificatie die niet aanwezig is in standaard SHAM.
• Roodverschuivingsevolutie: De studie identificeert en parametriseert de evolutie van de SHAM-parameters met de roodverschuiving. Specifiek neemt de satellietfractie ($f_{sat}$) voor de QSOs lineair toe met de roodverschuiving, terwijl de gemiddelde snelheid ($V_{mean}$) voor de ELGs een kwadratische trend volgt. Dit maakt een efficiëntere parametrisatie mogelijk dan het fitten van onafhankelijke parameters in elke bin.
• Validatie: De mocks worden rigoureus gevalideerd tegen DESI DR1 en DR2 data, waarbij de getalsdichtheden en clusteringstatistieken over schalen van $1$ tot $100 h^{-1}$Mpc worden gereproduceerd.

Resultaten

• Clusteringstatistieken: De gemodificeerde SHAM-mocks vertonen een uitstekende overeenstemming met de DESI DR2-data voor zowel de monopool als de kwadrupool van de roodverschuivingsruimte twee-punt correlatiefunctie (2PCF) en de geprojecteerde correlatiefunctie. De overeenstemming houdt stand over het fittingbereik ($3$ tot $70 h^{-1}$Mpc) en strekt zich uit tot de BAO-schaal ($\sim 100 h^{-1}$Mpc).
• Vermogensspectrum (Power Spectrum): De mocks reproduceren de vermogensspectrum-monopool en -kwadrupool redelijk goed, met kleine afwijkingen bij hoge golfgetallen ($k$) voor ELGs, wat te verwachten is aangezien het vermogensspectrum niet in het proces van het fitten is gebruikt.
• Halo Occupation Distribution (HOD):
  • ELGs: De HOD bestaat uit een Gaussische component voor de centrale stelsels (piekend bij $M_{halo} \approx 10^{12} h^{-1}M_{\odot}$ met een bezettingsgraad van 0.095) en een machtswet-component voor de satellieten (dominant bij $M_{halo} > 10^{12.4} h^{-1}M_{\odot}$).
  • QSOs: De HOD vertoont vergelijkbaar een Gaussische centrale component (piekend bij $M_{halo} \approx 10^{12.4} h^{-1}M_{\odot}$ met een bezettingsgraad van 0.319) en een machtswet-component voor de satellieten (dominant bij $M_{halo} > 10^{12.8} h^{-1}M_{\odot}$).
• Grootschalige Bias: De grootschalige bias ($b$) gemeten uit de mocks komt overeen met de bias afgeleid van de DESI DR2-data in het bereik van $10$ tot $80 h^{-1}$Mpc. Beide tracers vertonen een monotone toename van de bias met de roodverschuiving, wat consistent is met het bezetten van meer massieve halo's op eerdere kosmische tijdstippen.

Betekenis en Claims
Het artikel beweert dat deze resultaten een cruciale brug vormen tussen theorie en observatie voor DESI-analyses. Door de geobserveerde getalsdichtheden en clusteringpatronen van ELGs en QSOs te reproduceren, dienen deze mocks als essentiële instrumenten voor:

1. Schatting van de Covariantiematrix: Ze kunnen worden gebruikt om covariantiematrices te genereren voor toekomstige DESI-kosmologische analyses.
2. Testen van Systematiek: Ze maken het mogelijk om kosmologische pipelines te testen en de systematische onzekerheden die inherent zijn aan de survey te beoordelen.
3. Fysiek Inzicht: De succesvolle toepassing van gemodificeerde SHAM-technieken biedt een verbeterd begrip van de sterrenstelsel-halo-verbinding, met name met betrekking tot de omgevingsafhankelijkheid van stervorming in satellietstelsels (voor ELGs) en de evolutie van QSO-populaties.

De auteurs benadrukken dat hoewel de modellen zijn afgestemd op de DESI-data, ze zijn gebouwd op de Planck-kosmologie en de Uchuu-simulatie, wat ervoor zorgt dat ze fysiek consistent zijn met de huidige kosmologische standaarden. Het werk claimt geen nieuwe fysica te ontdekken, maar biedt de noodzakelijke infrastructuur (realistische mocks) om een nauwkeurige kosmologische inferentie uit de DESI-survey mogelijk te maken.