An improved model for the effect of correlated Si-III absorption on the one-dimensional Lyman-$α$ forest power spectrum

Technische Samenvatting: Een verbeterd model voor het effect van gecorreleerde Si III-absorptie op het eendimensionale Lyman-$\alpha$-bos power spectrum

Probleemstelling
Het eendimensionale Lyman-$\alpha$ (Ly$\alpha$)-bos power spectrum is een primair instrument voor het beperken van kosmologische parameters en modellen van donkere materie. De nauwkeurige modellering wordt echter bemoeilijkt door contaminatie van metaalabsorptielijnen, specifiek die gecorreleerd zijn met de Ly$\alpha$-bos. Terwijl ongecorreleerde metalen kunnen worden geschat en afgetrokken met behulp van sidebands, zijn gecorreleerde metalen zoals triplelijk geïoniseerd silicium (Si III, $\lambda1207$) vermengd binnen de bos.

Eerdere studies, met name McDonald et al. (2006, M06), modelleerden de Si III-contaminatie door aan te nemen dat het geabsorbeerde Si III-fluxcontrast een lineair geschaalde versie is van het Ly$\alpha$-fluxcontrast. Deze ansatz beschrijft succesvol de oscillerende kenmerken in het power spectrum veroorzaakt door de snelheidsverschuiving tussen Ly$\alpha$ en Si III op grote schalen ($k \lesssim 0.02$ s km$^{-1}$). Echter, dit artikel identificeert dat het M06-model bezwijkt bij kleinere schalen ($k > 0.02$ s km$^{-1}$). Specifiek faalt de M06-aanpak om rekening te houden met:

1. Distinctieve lijnprofielen: Si III-lijnen zijn thermisch smaller dan Ly$\alpha$-lijnen vanwege de hogere massa van silicium, wat leidt tot een versterking van de kleine-schaal power die het lineaire rescalingsmodel mist.
2. Variabele optische diepteverhoudingen: De verhouding tussen coëvale Si III- en Ly$\alpha$-optische dieptes is niet constant, maar hangt af van gasdichtheid, temperatuur en de abundantie van silicium.
3. Redshift-evolutie: De amplitude van de Si III-contaminatie volgt mogelijk niet de eerder aangenomen schaling met de gemiddelde Ly$\alpha$-transmissie, met name door het verharden van de metagalactische UV-achtergrond bij $z \lesssim 3.5$.

Methodologie
De auteurs maken gebruik van de Sherwood-Relics suite van hoog-resolutie kosmologische hydrodynamische simulaties om mock quasar-spectra te construeren.

• Simulatie-opzet: De simulaties gebruiken een $40 h^{-1}$ cMpc box met $2 \times 2048^3$ deeltjes. Ze maken gebruik van een niet-evenwicht thermo-chemische solver en een ruimtelijk uniforme UV-achtergrond.
• Si III Modellering: In tegen tegenstelling tot eerdere simulaties die mogelijk geen modellen voor metaalproductie bevatten, nemen de auteurs een fiduciaal silicium-abundantiemodel aan dat afhankelijk is van gasoverdichtheid ($\Delta$) en redshift ($z$), gekalibreerd tegen observationele beperkingen van Schaye et al. (2003) en Aguirre et al. (2004). De Si III ionisatiefractie ($x_{\text{Si III}}$) wordt berekend met behulp van de foto-ionisatiecode Cloudy in combinatie met het UV-achtergrondspectrum.
• Constructie van Spectra: De auteurs extraheren 5.000 periodieke gezichtslijnen (sightlines) bij redshift-intervallen van $\Delta z = 0.1$. Ze berekenen de Ly$\alpha$- en Si III-optische dieptes, waarbij Voigt-profiel convoluties worden toegepast (waarbij Si III-lijnen op een fijner raster worden hersampled om hun smallere thermische breedtes te resolveren).
• Power Spectrum Analyse: Ze berekenen het 1D power spectrum ($P(k)$) voor de totale flux, de Ly$\alpha$-alleen en de Si III-alleen componenten. Ze vergelijken deze directe simulatie-resultaten met het analytische M06-model en hun voorgestelde nieuwe passende functie.

Belangrijkste Bijdragen en Modelontwikkeling
Het artikel stelt een nieuwe analytische passende functie (Eq. 18) voor om de M06-ansatz te vervangen voor het modelleren van Si III-contaminatie. Het nieuwe model introduceert vier vrije parameters om twee cruciale fysische effecten te vangen die in eerdere werken werden verwaarloosd:

1. Schaalafhankelijke Versterking: Het model bevat exponentiële termen ($\exp(k/k_{\text{auto}}$ en $\exp(k/k_{\text{cross}})$) om de smallere Si III-lijnprofielen en de pixel-afhankelijke verhouding van optische dieptes te verklaren.
2. Decorrelatie: Het model staat decorrelatie toe tussen Si III- en Ly$\alpha$-absorptie op kleine schalen. Dit gebeurt omdat sterke Si III-absorbers vaak overeenkomen met verzadigde Ly$\alpha$-lijnen, wat de lineaire relatie die in M06 wordt aangenomen, doorbreekt.

Het voorgestelde totale power spectrum wordt gegeven door:
$$P_{\text{tot}}(k) = P_{\text{Ly}\alpha}(k) + P_{\text{Si III}}(k) + P_{\text{Ly}\alpha-\text{Si III}}(k)$$
Waarbij:

• $P_{\text{Si III}}(k) = a_{\text{auto}} \exp(k/k_{\text{auto}}) P_{\text{Ly}\alpha}(k)$
• $P_{\text{Ly}\alpha-\text{Si III}}(k) = a_{\text{cross}} \cos(\Delta v k) \exp(k/k_{\text{cross}}) P_{\text{Ly}\alpha}(k)$

Resultaten

• Validatie: Het nieuwe model werd gevalideerd tegen de Sherwood-Relics simulaties over redshifts $2.2 \le z \le 5.0$ en wavenumbers $k < 0.2$ s km$^{-1}$.
• Prestaties: Op grote schalen ($k < 0.02$ s km$^{-1}$) passen zowel het M06-model als het nieuwe model de data goed. Echter, op kleinere schalen onderschat het M06-model de power aanzienlijk, terwijl het nieuwe model de simulatiegegevens reproduceert met een fractionele afwijking van minder dan 1% voor het fiduciaal geval.
• Parametergevoeligheid:
  • De schaalparameters ($k_{\text{auto}}$ en $k_{\text{cross}}$) zijn primair gevoelig voor de gastemperatuur ($T_0$). Heter gas leidt tot bredere lijnen en andere decorrelatieschalen.
  • De amplitude parameters ($a_{\text{auto}}$ en $a_{\text{cross}}$) zijn primair gevoelig voor de silicium-abundantie ($[\text{Si/H}]$).
• Redshift-evolutie: In tegen tegenstelling tot de M06-aanname dat de Si III-power schaalt met de inverse van de gemiddelde transmissie ($1 - \bar{F}$), laten de simulaties zien dat de Si III-power spectrum amplitude toeneemt naar lagere redshifts ($z \lesssim 3.5$). Dit wordt gedreven door het verharden van de UV-achtergrond, wat de Si III ionisatiefractie ($x_{\text{Si III}}$) bij lagere redshifts vermindert, maar dit wordt gecompenseerd door de evolutie van het gebruikte silicium-abundantiemodel.

Betekenis en Claims
De auteurs stellen dat hun verbeterde model een essentiële overweging is voor toekomstige, hoog-precieze Ly$\alpha$-bos power spectrum metingen, met name die gebruikmaken van hoog-resolutie spectrografen (bijv. Keck/HIRES, VLT/UVES) die schalen voorbij $k \sim 0.02$ s km$^{-1}$ onderzoeken.

Wat betreft kosmologische restricties, claimen de auteurs dat het adopteren van een fysisch gemotiveerde prior op de silicium-abundantie betekent dat het nieuwe model weinig effect heeft op bestaande warm donkere materie (WDM) restricties afgeleid van de Ly$\alpha$-bos. De auteurs benadrukken echter dat het negeren van de kleine-schaal power-versterking van Si III tot systematische fouten zou kunnen leiden in toekomstige analyses naarmate de statistische onzekerheden afnemen. Het artikel stelt geen nieuwe observationele campagnes voor, maar biedt het noodzakelijke theoretische kader om toekomstige hoog-resolutie data accuraat te interpreteren.