Comparing cosmic shear nulling methods for Stage-IV surveys

Technische Samenvatting: Vergelijking van Cosmic Shear Nulling-methoden voor Stage-IV-surveys

Probleemstelling
Cosmic shear-metingen van Stage-IV-surveys (bijv. Euclid, Rubin Observatory) hebben als doel de fysica van donkere energie tot op procentniveau te beperken. Deze metingen worden echter verontreinigd door baryonfeedbackprocessen (bijv. AGN-feedback, supernova's), die het vermogensspectrum van materie onderdrukken op kleine schalen ($k \gtrsim 10 \, \mathrm{Mpc}^{-1}$). Huidige hydrodynamische simulaties voorspellen onderdrukkingsamplitudes die aanzienlijk variëren (10% tot >30%), en subgrid-fysica-prescripties blijven onzeker. Standaardanalyse mitigeert dit door hoekschalensnijden toe te passen (beperking van $\ell_{\max}$), maar deze aanpak is suboptimaal omdat een gegeven hoekschaal bijdragen integreert uit een breed scala aan fysieke 3D-schalen als gevolg van de brede radiale lenskern. Bijgevolg zijn agressieve schaalverwijderingen vereist om onbevooroordeelde kosmologische beperkingen te garanderen, wat leidt tot een aanzienlijk verlies aan statistische precisie. Dit artikel onderzoekt of "nulling"-transformaties toegepast op de datavektor de bijdragen van hoge-$k$ (kleine schaal) effectiever kunnen isoleren en verwijderen dan standaard schaalverwijderingen.

Methodologie
De auteurs voeren een Fisher-prognoseanalyse uit waarbij drie verschillende nulling-strategieën worden vergeleken met een standaard, niet-genuleerde cosmic shear-analyse. De analyse maakt gebruik van een mock datavektor gebaseerd op een Euclid-achtige surveyopstelling met zes tomografische bronbinnen, een dichtheid van 30 sterrenstelsels/boogminuut$^2$, en een hemeloppervlak van 15.000 deg$^2$. Baryonfeedback wordt gemodelleerd met HMCode2020 met een fiduciaire AGN-verwarmingsparameter ($\Theta_{\mathrm{AGN}} = 7.8$). De auteurs gaan uit van een $w_0w_a$CDM-kosmologie en stellen alle systematiek (bijv. fotometrische roodverschuivingen, shear-bias) vast op nul om de impact van baryonfeedback te isoleren.

De drie vergeleken methoden zijn:

1. LU Nulling (LUnul): Een nieuwe implementatie van een methode die oorspronkelijk werd voorgesteld door Huterer & White (2005). Het behandelt het Limber-integraal als een discrete matrixvermenigvuldiging. Door een LU-factorisatie van de kernmatrix uit te voeren, construeert de methode lineaire combinaties van de tomografische vermogensspectra. Dit schakelt effectief de gemeenschappelijke hoge-$k$-bijdragen uit (die vergelijkbaar zijn tussen binnen bij lage roodverschuiving) terwijl verschillen bij lagere $k$ behouden blijven. De methode staat een "zachte" snij toe waarbij specifieke rijen van de getransformeerde datavektor worden verworpen op basis van een tolerantieparameter ($tol=0.01$) voor hoge-$k$-bijdragen.
2. BNT-transformatie: De Bernardeau-Nishimichi-Taruya-transformatie werkt op het lensveld (shear-kaarten). Het gebruikt een lineaire, inverteerbare transformatie van de tomografische binnenindices om de lenskernen in roodverschuivingsruimte te lokaliseren. Dit creëert een schonere mapping tussen hoekmultipool $\ell$ en 3D-golfgetal $k$, waardoor kruiscorrelaties tussen niet-aangrenzende binnen worden onderdrukt.
3. Kruiscorrelatie: Aan de CMB-lensingliteratuur ontleend, gebruikt deze methode een voorgrondsterrenstelselclusteringsteekproef (laag-roodverschuiving lensbinnen) om het shear-signaal te decorreleren. Door het geprojecteerde dichtheidsveld (gewogen met het kruiscorrelatiespectrum) van het shear-veld af te trekken, wordt de bijdrage van laag-roodverschuiving, kleine-schaalstructuren onderdrukt. In deze studie worden de twee laagste roodverschuivingsbinnen herbestemd als lenzen om de drie hoogste roodverschuivingsbronbinnen te "schoonmaken".

Belangrijkste bijdragen

• Implementatie van LUnul: De auteurs bieden een publieke implementatie van de LU-factorisatiemethode voor cosmic shear, een techniek die sinds de oorspronkelijke voorstellen niet eerder in de literatuur is onderzocht.
• Toepassing van kruiscorrelatie op cosmic shear: Het artikel past de kruiscorrelatie-nullingtechniek, eerder gebruikt voor CMB-lensing, toe op de context van Stage-IV cosmic shear gecombineerd met fotometrische sterrenstelselclustering.
• Systematische vergelijking: De studie biedt een directe vergelijking van deze drie verschillende benaderingen (veldgebaseerd, statistiekgebaseerd en tracergebaseerd) binnen een verenigd Fisher-prognosekader, waarbij hun effectiviteit wordt beoordeeld in het verminderen van bias terwijl informatie behouden blijft.

Resultaten
De analyse evalueert de afweging tussen parameterbias (nauwkeurigheid) en foutmarges (precisie) over een reeks maximale golfgetallen ($k_{\max}$):

• Biasreductie: Alle drie de methoden verminderen effectief de bias op kosmologische parameters ($S_8, \Omega_m, w_0, w_a$) veroorzaakt door niet-gemodelleerde baryonfeedback in vergelijking met de standaardcase.
• LUnul-prestaties: Deze methode is het meest agressief in het onderdrukken van hoge-$k$-modi. Het bereikt de laagste bias (vermindering van $S_8$-bias tot $<0.5\sigma$ bij $k_{\max} = 0.3 \, h/\mathrm{Mpc}$), en presteert beter dan andere methoden qua nauwkeurigheid. Dit gaat echter ten koste van precisie; het beperkingsvermogen (Figure of Merit, FoM) is aanzienlijk lager dan de standaardcase bij vergelijkbare schaalverwijderingen omdat het meer informatie verworpt.
• BNT-prestaties: De BNT-transformatie behoudt meer informatie dan LUnul en heeft een rigoureuze theoretische onderbouwing. Echter, in deze specifieke opstelling met slechts zes tomografische binnen, toont het slechts een marginale verbetering ten opzichte van de standaardcase. De auteurs merken op dat BNT bekend staat als effectiever met een groter aantal binnen.
• Kruiscorrelatie-prestaties: Deze methode is effectief in het verminderen van absolute biases in $S_8$, maar vereist aanvullende informatie van sterrenstelselclustering. Hoewel het strakkere foutmarges produceert (door de opname van dichtheidsveldinformatie), blijft de relatieve bias (bias gedeeld door fout) hoog. Bovendien is het minder beperkend voor donkere-energieparameters ($w_0, w_a$) in vergelijking met andere methoden. De methode wordt beperkt door de overlap tussen lens- en bronsteekproeven in fotometrische surveys.
• Parametervolatiliteit: $S_8$ is consistent de meest bevooroordeelde parameter over alle methoden heen vanwege zijn gevoeligheid voor de amplitude van het vermogensspectrum.

Betekenis en claims
Het artikel concludeert dat nullingtechnieken de door baryonfeedback geïntroduceerde bias weliswaar aanzienlijk kunnen mitigeren, maar dit over het algemeen doen ten koste van statistische precisie. De auteurs beweren dat:

1. LUnul de meest robuuste methode is voor het verwijderen van bias en het doordringen tot kleinere schalen, hoewel het zorgvuldige numerieke afstemming vereist en precisie opoffert.
2. BNT theoretisch onderbouwd is en informatie behoudt, maar mogelijk meer tomografische binnen vereist (karakteristiek voor toekomstige Stage-IV-datareleases) om zijn volledige potentieel boven standaardverwijderingen te realiseren.
3. Kruiscorrelatie een haalbare weg biedt om bias te verminderen met bestaande 3$\times$2-pt-gegevens, maar wordt beperkt door de moeilijkheid om lens- en populaties perfect te scheiden in fotometrische surveys.

De auteurs benadrukken dat deze methoden een alternatief bieden voor complexe baryonfeedbackmodellering. Zij stellen dat Stage-IV-surveys waarschijnlijk toch sterke schaalverwijderingen zullen moeten toepassen of vertrouwen moeten hebben op nauwkeurige baryonmodellen, ongeacht de gekozen nulling, maar dat deze nullingstrategieën een manier bieden om datavektoren te construeren die immuner zijn voor kleine-schaal astrofysische onzekerheden. De resultaten bevestigen de gevestigde afweging tussen nauwkeurigheid (biasreductie) en precisie in cosmic shear-analyses.