Optimal Transport Reconstruction of Biased Tracers in Primordial Non-Gaussian Fields

Technische Samenvatting: Optimal Transport Reconstructie van Vertekende Tracers in Primordiale Niet-Gaussiaanse Velden

Probleemstelling
Optimal Transport (OT) heeft zichzelf gevestigd als een efficiënt raamwerk voor het reconstrueren van het kosmische verplaatsingsveld, waarbij de initiële posities worden gekoppeld aan de huidige locaties (of vice versa). Voorgaand werk heeft aangetoond dat OT de verplaatsingen van vertekende (biased) tracers nauwkeurig herstelt onder de aanname van Gaussiaanse initiële condities. Een kritieke openstaande vraag is echter: blijft dit raamwerk geldig wanneer de initiële condities niet-Gaussiaans zijn? Dit is een actuele uitdaging, gezien de toenemende inspanningen om de primordiale niet-Gaussiaanseheid (PNG) te beperken met behulp van kosmologische surveys van de volgende generatie. Een specifieke moeilijkheid ontstaat omdat OT-implementaties vaak uitgaan van een uniforme initiële ruimtelijke verdeling. Ho'wel geldt dit voor het volledige donkere materieveld, maar niet voor vertekende deelverzamelingen (zoals halo's). In Gaussiaanse scenario's kan de onwaargenomen massa ("stof") effectief worden gemodelleerd; echter, in niet-Gaussiaanse scenario's vertoont het stof zelf specifieke schaalafhankelijke vertekkingen die rekening moeten worden gehouden met om een nauwkeurige reconstructie te garanderen.

Methodologie
De auteurs onderzoeken de toepasbaarheid van OT-reconstructie op "lokale" primordiale niet-Gaussiaanse velden, waarbij het gravitationele potentieel $\phi_{NG}$ wordt gemodelleerd als een machtsreeks in een Gaussiaans veld $\phi_G$ (getrunceerd in tweede orde, $f_{NL} \neq 0$). In dit model bevat de dichtheidscontrast $\delta_{NG}$ een term proportioneel aan het product van Gaussiaanse velden, wat leidt tot een kenmerkende schaalafhankelijke vertekking in de clustering van vertekende tracers op grote schalen (kleine $k$).

De studie maakt gebruik van de Quijote-PNG simulatieset, bestaande uit 20 periodieke boxen ($L = 1 h^{-1}$ Gpc) met $512^3$ deeltjes, die evolueren van $z=127$ naar $z=0$ onder een vlak $\Lambda$CDM-kosmologie. Vertekende tracers worden geïdentificeerd als halo's met een massa $> 20 m_p$.

De kern van de methodologische innovatie betreft het modelleren van het onwaargenomen "stof" (het massafractie $1-f_b$ dat niet in de halo-catalogus zit) om te voldoen aan de restrictie dat het volledige materieveld geen schaalafhankelijke vertekking mag vertonen. De auteurs leiden af dat als vertekende tracers een schaalafhankelijke vertekking vertonen die wordt gekarakteriseerd door parameters $b_\delta$ en $b_\phi$ (waarbij de vertekking schaalt als $b(k) \approx b_\delta + 2f_{NL}b_\phi \alpha(k)$), het stof een compenserende schaalafhankelijke vertekking met een tegengesteld teken moet bezitten.

Twee reconstructiebenaderingen worden getest:

1. OT-K: Gebruikt een Wiener-filter model om de stofverdeling te schatten op basis van waargenomen tracers. Dit vereist het schatten van de combinatie $2f_{NL}b_\phi$ (die geconstreerd kan worden door survey-data) en het massafractie $f_b$.
2. OT-W: Gebruikt de werkelijke stofdeeltjes uit de simulatie (een willekeurige steekproef) in combinatie met de massa-gewogen halo's om de reconstructie uit te voeren, wat dient als een bovengrens voor de prestaties.

De kwaliteit van de reconstructie wordt geëvalueerd door de gereconstrueerde posities te vergelijken met de werkelijke protohalo-posities en door de schaalafhankelijkheid van de divergentie van het verplaatsingsveld ($\theta_{OT}$) te analyseren, wat dient als een proxy voor het initiële dichtheidsveld.

Belangrijkste Bijdragen en Resultaten

• Karakterisering van Stof: Het artikel toont aan dat in lokale PNG-modellen het "stof" niet een passieve achtergrond is, maar zijn eigen karakteristieke schaalafhankelijke vertekking vertoont. Specifiek wordt de amplitude van de vertekking van het stof bepaald door de vereiste dat de som van het vertekende tracer- en stofveld een schaal-onafhankelijk totaalveld oplevert.
• Uitbreiding van OT naar Niet-Gaussiaanseheid: De auteurs tonen aan dat de OT-methodologie succesvol wordt uitgebreid naar niet-Gaussiaanse initiële condities. Door een realistisch model voor het stof te integreren dat rekening houdt met de $f_{NL}$-afhankelijke vertekking, blijft de reconstructie van vertekende tracers nauwkeurig.
• Nauwkeurigheid van de Reconstructie:
  • Tracer Posities: De OT-K reconstructie (met het Wiener-filter stofmodel) herstelt correct de schaalafhankelijke vertekking van de protohalo's. De OT-W reconstructie (met echt stof) vertoont een nog betere overeenkomst, wat aangeeft dat de methode robuust is maar gevoelig voor de getrouwheid van het stofmodel.
  • Volledig Veld Proxy: De divergentie van het OT-verplaatsingsveld ($\theta_{OT}$) vertoont bijna geen schaalafhankelijkheid bij kleine $k$, wat het gedrag van het werkelijke initiële veld weerspiegelt. Het OT-W model breidt deze schaal-onafhankelijkheid uit naar grotere $k$-waarden vergeleken met OT-K, waarbij een factor 2 verbetering wordt bereikt in het bereik waar $b(k) \approx 1$.
• Robuustheid ten opzichte van $f_{NL}$: De gemiddelde kwadratische waarde van het gereconstrueerde veld $\langle \theta_{OT}^2 \rangle$ blijft stabiel over verschillende waarden van $f_{NL}$ ($+100, 0, -100$), met variaties binnen de standaarddeviatie, wat suggereert dat de methode niet beïnvloed wordt door de aanwezigheid van niet-Gaussiaanseheid zelf.

Betekenis en Claims
Het artikel beweert dat Optimal Transport een nuttig raamwerk biedt voor het reconstrueren van kosmische structuren, zelfs wanneer de initiële condities niet-Gaussiaans zijn, mits het model voor het onwaargenomen "stof" voldoende realistisch is. De auteurs benadrukken dat het succes van de methode berust op het feit dat de kenmerkende schaalafhankelijke signatuur van PNG aanwezig is in zowel de vertekende tracers als het stof, maar met tegengestelde tekens, waardoor het volledige veld onvertekend kan blijven.

De auteurs stellen bescheiden dat zij niet direct een verbetering van de restricties op $f_{NL}$ uit dit werk alleen demonstreren. In plaats daarvan bieden zij een fundamenteel raamwerk dat toestaat dat reconstructie-gebaseerde analyses consistent kunnen worden toegepast in PNG-kosmologieën. Zij suggereren dat dit werk een noodzakelijke stap is naar het integreren van reconstructietechnieken in toekomstige hoog-precisie PNG-analyses, wat potentieel kan helpen om geobserveerde vertekende tracers dichter bij hun primordiale Lagrangiaanse tegenhangers te brengen en subtiele niet-lineaire of stochastische bijdragen te mitigeren die relevant zijn voor surveys van de volgende generatie. Het artikel concludeert dat verbeterde modellering van het stof (bijv. via iteratieve methoden of machine learning) verdere winst zou kunnen opleveren.