Forecast Cosmological Constraints with the 1D Wavelet Scattering Transform and the Lyman-$α$ forest

Probleemstelling
De Lyman-$\alpha$-boson (Lyman-$\alpha$ forest) dient als een cruciale sonde voor kleinschalige structuren en biedt inzichten in donkere materie, neutrino-massa, inflatoire fysica en de thermische geschiedenis van het intergalactische medium (IGM). Het extraheren van kosmologische informatie uit dit niet-lineaire regime is echter uitdagend. Huidige analyses die vertrouwen op het 1D-flux-vermogensspectrum worden systematisch gedomineerd in grote surveys zoals SDSS en DESI. Bovendien lijden traditionele hogere-orde statistieken, zoals de bispectrum, aan een verhoogde variantie en een verminderde robuustheid tegen uitschieters in sterk niet-Gaussische verdelingen. Er is een aangetoonde behoefte aan nieuwe samenvattende statistieken die in staat zijn om niet-Gaussische informatie te extraheren uit het niet-lineaire dichtheidsveld zonder de trainingsvereisten of interpreteerbaarheidsproblemen die geassocieerd worden met deep learning-benaderingen.

Methodologie
De auteurs stellen het gebruik van de 1D Wavelet Scattering Transform (WST) voor als een samenvattende statistiek voor de Lyman-$\alpha$-boson kosmologie. De WST is een deterministische, wiskundig robuuste signaalverwerkingstechniek die het invoefluxveld convolveert met complex-waardige wavelets op verschillende schalen, gevolgd door modulus-operaties en laagdoorlaatfiltering. Dit proces scheidt informatie over schalen en extraheert hogere-orde statistieken (tot de tweede orde in deze studie) die stabiel zijn tegen kleine deformaties en ruis.

Om de voorspellende kracht van de WST te voorspellen, maakten de auteurs gebruik van:

• Simulaties: Een reeks hydrodynamische simulaties (met behulp van MP-Gadget) bevattende $2 \times 256^3$ deeltjes in een box van $40 \text{ Mpc}/h$. Hiervan werden 32.000 Lyman-$\alpha$-absorptiespectra gegenereerd bij roodverschuivingen $z=2, 3, 4$.
• Parameters: De analyse richtte zich op vier kosmologische en thermische parameters: de spectrale index ($n_P$) en amplitude ($A_P$) van het primordiale vermogensspectrum, en IGM-thermische parameters ($H_A$ en $H_S$).
• Ruismodellering: Realistische observationele ruis werd gesimuleerd door Gaussische ruis toe te voegen aan elk spectrum, waarbij de continuüm-ruisverhouding (CNR) werd getrokken uit een log-normale verdeling die gefit is aan SDSS DR16Q-data.
• Statistisch Kader: Een Fisher-informatie-matrix werd geconstrueerd om parameterbeperkingen te schatten. De covariantie-matrix werd geschat met behulp van 2.000 willekeurige ruisrealisaties. De auteurs vergeleken de uit de WST-coëfficiënten afgeleide Fisher-informatie met die van het standaard 1D-flux-vermogensspectrum.

Belangrijkste Bijdragen en Resultaten
De studie toont aan dat de WST significante kosmologische informatie vastlegt die niet aanwezig is in het 1D-flux-vermogensspectrum.

• Informatiegehalte: De tweede-orde wavelet-coëfficiënten, die niet-Gaussische informatie onderzoeken, dragen ongeveer 50% extra informatie bij aan de Fisher-matrix vergeleken met enkel eerste-orde coëfficiënten.
• Parameterbeperkingen: De WST levert substantieel nauwere beperkingen op kosmologische parameters dan het vermogensspectrum. Specifiek worden de geschatte posterieure onzekerheden ($\sigma_X$) voor parameters $n_P$, $A_P$, $H_S$ en $H_A$ met een factor 30 tot 60 verbeterd wanneer de gecombineerde WST-coëfficiënten worden gebruikt.
• Dimensiereductie: De auteurs schatten dat de WST de 4D-parameterruimte met een factor $10^6$ verkleint vergeleken met het vermogensspectrum.
• Gevoeligheid: De WST-coëfficiënten zijn minder gevoelig voor de gemiddelde flux ($\tau_{\text{eff}}$) dan het vermogensspectrum (verandering van $\sim 2,5\%$ voor een 10% fluxverandering, versus $\sim 7\%$ voor het vermogensspectrum), wat wijst op een potentiële reductie van de degeneratie tussen $\tau_{\text{eff}}$ en de primordiale amplitude $A_P$.
• Roodverschuivingsafhankelijkheid: Het informatiegehalte piekt bij $z=3$ voor de meeste parameters, terwijl de thermische parameter $H_A$ het meest wordt beperkt bij $z=2$.

Betekenis en Claims
Het artikel beweert dat de 1D Wavelet Scattering Transform een nieuwe en krachtige benadering biedt voor het extraheren van kosmologische informatie uit de Lyman-$\alpha$-boson, waarbij deze de traditionele 1D-flux-vermogensspectrum met een orde van grootte of meer overtreft in termen van beperkingskracht. De auteurs stellen dat, indien deze verbeteringen gerealiseerd worden in de komende data van het Dark Energy Spectroscopic Instrument (DESI), de WST in staat kan zijn om:

• Beperkingen op de inflatoire running ($\alpha_s$) te leveren die voldoende zijn om veelvoorkomende inflatoire modellen te testen die $\alpha_s = -6 \times 10^{-4}$ voorspellen.
• De neutrino-massa beperkingen zodanig te verbeteren dat een $5\sigma$-detectie van de minimale neutrino-massa kan worden bereikt.

De auteurs blijven bescheiden over de huidige beperkingen van hun voorspelling, waarbij zij opmerken dat de Fisher-matrix-benadering parameterdegeneraties verwaarloost en dat de kleine simulatieboxgrootte waarschijnlijk leidt tot pessimistische voorspellingen voor de spectrale index $n_P$. Zij erkennen dat een volledige likelihood-analyse en een robuuste behandeling van systematische fouten (zoals continuüm-fitting) noodzakelijke toekomstige stappen zijn voordat deze methoden op echte DESI-data kunnen worden toegepast. Desalniettemin vestigt dit werk de basis voor het feit dat de WST aanzienlijk meer informatie behoudt dan het vermogensspectrum, wat de mogelijkheid opent om parameterbeperkingen met meerdere ordes van grootte te verbeteren.