EMU: Cross-correlating EMU Pilot Survey 1 with Dark Energy Survey to validate the radio galaxy bias and redshift distribution

Technische Samenvatting: EMU Cross-correlatie met DES voor Validatie van Radio-stelsel Bias en de Redshift-distributie

Probleemstelling
Radio-continuüm surveys bieden een krachtige methode voor het in kaart brengen van de grootschalige structuur (LSS) van het Universum vanwege hun vermogen om grote delen van de hemel te dekken en hun ongevoeligheid voor stofextinctie. Echter, een kritieke beperking voor kosmologische analyses met deze surveys is het gebrek aan directe redshift-informatie voor radiobronnen. Het dominante emissiemechanisme, synchrotronstraling, produceert een kenmerkloos machtswet-spectrum, waardoor schatting van de redshift extreem uitdagend is. Zonder een nauwkeurige redshift-distributie ($dN/dz$) en een robuust model voor galaxy bias, kunnen systematische fouten doorwerken in kosmologische parameterbeperkingen, zoals $\sigma_8$. Hoewel eerdere studies vertrouwden op cross-matching met optische surveys (wat fotometrische redshifts oplevert voor slechts een fractie van de bronnen) of simulaties (bijv. TRECS, SKADS), is er behoefte aan een methode die het volledige survey-oppervlak kan gebruiken om statistisch de redshift-distributie af te leiden zonder afhankelijk te zijn van individuele bronidentificatie.

Methodologie
Deze studie maakt gebruik van een cross-correlatietechniek tussen de tweedimensionale radio-continuümkaart van de Evolutionary Map of the Universe Pilot Survey 1 (EMU-PS1) en de fotometrische optische galaxy-sample van de Dark Energy Survey (DES). De analyse verloopt als volgt:

1. Data-voorbereiding:

   • EMU-PS1: Maakt gebruik van een catalogus gegenereerd door de PyBDSF source finder, bevattende $\sim$224.100 bronnen. De analyse past flux-dichtheid-snijpunten toe bij 180 $\mu$Jy en 400 $\mu$Jy om conservatieve samples te creëren. Weight maps worden geconstrueerd om rekening te houden met observationele ruis en variaties in bronenaantallen.
   • DES: Gebruikt de magnitude-limited (MagLim) galaxy sample, verdeeld in zes tomografische redshift bins ($0.20 \le z < 1.05$). Fotometrische redshifts worden afgeleid met het Directional Neighbourhood Fitting (DNF) algoritme.
   • Planck: Bevat de Planck CMB lensing convergence map ($\kappa$) om een onafhankelijke validatie van de radio bias en de redshift-distributie te bieden.

2. Power Spectrum Schatting:

   • Angular power spectra ($C_\ell$) worden berekend met een pseudo-$C_\ell$ estimator (NaMaster) om incomplete hemeldekking op te vangen.
   • De analyse meet het auto-power spectrum van EMU-PS1 ($C^{EMU}_\ell$), de auto-power spectra van DES bins ($C^{DES}_\ell$), de cross-power spectra tussen EMU-PS1 en DES bins ($C^{EMU \times DES}_\ell$), en het cross-power spectrum tussen EMU-PS1 en Planck lensing ($C^{EMU \times \kappa}_\ell$).
   • Een covariantie-matrix wordt geschat met behulp van 1.000 mock realisaties gegenereerd via GLASS, waarbij rekening wordt gehouden met de survey footprints en Poisson-ruis.

3. Modellering en Inferentie:

   • Redshift-distributie ($n(z)$): Twee modellen worden getest. Eerst de redshift-distributie van de Tiered Radio Extragalactic Continuum Simulation (TRECS), die bijdragen bevat van Active Galactic Nuclei (AGN's) en Star-Forming Galaxies (SFG's). Ten tweede een parametrisch model (geadapteerd van Nakoneczny et al. 2024) met drie vrije parameters ($z_0$, $r$, $\gamma$), ontworpen om de piek en de staart van de distributie fysiek te motiveren.
   • Galaxy Bias: Lineaire, schaalonafhankelijke bias-modellen worden gehanteerd. Voor de parametrische $n(z)$ wordt een constant amplitude bias-model gebruikt voor SFG's en een polynomiale evolutie voor AGN's (of een gecombineerde effectieve bias).
   • Parameter-schatting: Een maximum likelihood framework (gebruikmakend van PyMC) wordt ingezet om de bias-parameters en de redshift-distributie parameters te beperken door de gezamenlijke data-vector van de cross-power spectra te fitten.

4. Validatie:

   • De analyse-pipeline wordt gevalideerd met mock catalogi waar de ware $n(z)$ en bias bekend zijn. Dit bevestigt het vermogen van de methode om de locatie van de piek van de redshift-distributie te herstellen.
   • De impact van residuele incompletheid in de radio-catalogus wordt beoordeeld, wat aantoont dat hoewel auto-power spectra gevoelig zijn voor incompletheid, de cross-power spectra die voor inferentie worden gebruikt, robuust blijven.

Belangrijkste Resultaten

• Detectie Significantie: De cross-correlatie tussen EMU-PS1 en de zes DES tomografische bins wordt gedetecteerd met significante waarden variërend van $3.5\sigma$ tot $5.5\sigma$ in het multipool-bereik $50 < \ell < 250$.
• DES Bias Validatie: De galaxy bias die voor DES bins is geschat met enkel clustering-metingen, is consistent met eerdere DES-resultaten (bijv. Abbott et al. 2022) binnen $1\sigma$, wat de analyse-pipeline valideert.
• Redshift-distributie Herstel:
  • TRECS Model: De gesimuleerde $n(z)$ van TRECS, wanneer gecombineerd met een mixture model voor galaxy bias, past goed op de cross-power spectra, het EMU-PS1 auto-power spectrum en de cross-correlatie met CMB lensing. De totale $\chi^2$ waarden zijn 7.78 (180 $\mu$Jy) en 3.15 (400 $\mu$Jy) voor 20 vrijheidsgraden.
  • Parametrisch Model: Het parametrische model (Eq. 11) biedt ook een goede fit aan de cross-power spectra zonder dat voorkennis van de AGN/SFG-mixture vereist is. De best-fit distributie piekt bij $z \sim 0.85$. Deze piek is hoger dan wat voorspeld wordt door TRECS en LOFAR deep fields, maar is consistent met de COSMOS 3GHz en MIGHTEE surveys.
• Bias Evolutie: De herstelde galaxy bias-evolutie voor EMU-PS1 is consistent met zowel de TRECS-voorspellingen als onafhankelijke schattingen uit de literatuur, zelfs wanneer het parametrische model wordt gebruikt dat de verschillende bronpopulaties niet expliciet scheidt.
• Systematiek: De studie merkt een discrepantie op waarbij het parametrische model faalt om het EMU-PS1 auto-power spectrum te fitten bij de 180 $\mu$Jy snede, wat wijst op residuele systematiek in de radio-data (mogelijk incompletheid) die de cross-correlatie resultaten niet beïnvloedt.

Betekenis en Claims
De auteurs beweren dat dit werk aantoont dat het haalbaar is om de redshift-distributies en de galaxy bias van grootschalige radio-continuüm surveys direct te schatten via cross-correlaties met optische surveys, zonder afhankelijk te zijn van individuele bron-redshifts of enkel simulaties.

• Validatie van Methode: De cross-correlatietechniek slaagt erin de vorm van de redshift-distributie te herstellen en het bias-model te valideren, wat een complementaire benadering biedt voor directe optische matching.
• Kosmologische Bruikbaarheid: De resultaten valideren het gebruik van het parametrische $n(z)$-model en de TRECS-simulatie voor toekomstige kosmologische analyses. Dit is significant voor aankomende grootschalige surveys zoals de volledige EMU-survey en de Square Kilometre Array Observatory (SKAO), waarbij directe redshift-metingen voor alle bronnen onmogelijk zullen zijn.
• Robuustheid: De methodologie blijkt robuust tegen interne systematiek en residuele incompletheid in de radio-catalogus bij het gebruik van cross-power spectra, hoewel de auteurs waarschuwen dat het kleine survey-oppervlak van de pilot ($\sim$270 deg$^2$) de precisie van de beperkingen limiteert in vergelijking met wat bereikbaar zal zijn met de hoofd-survey.

Het artikel concludeert dat cross-correlatie met optische surveys een levensvatbare en noodzakelijke tool is voor het precisie-kosmologie tijdperk van de radioastronomie, wat de mitigatie van redshift-distributie onzekerheden mogelijk maakt die momenteel de wetenschappelijke opbrengst van radio-continuüm data beperken.