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

Sintesi Tecnica: Un Modello Migliorato per l'Effetto dell'Assorbimento Correlato di Si III sullo Spettro di Potenza Unidimensionale della Foresta Lyman-$\alpha$

Definizione del Problema
Lo spettro di potenza unidimensionale della foresta Lyman-$\alpha$ (Ly$\alpha$) è uno strumento primario per vincolare i parametri cosmologici e i modelli di materia oscura. Tuttavia, una modellazione accurata è complicata dalla contaminazione da linee di assorbimento metalliche, specificamente quelle correlate con la foresta Ly$\alpha$. Mentre i metalli non correlati possono essere stimati e sottratti utilizzando i sideband, i metalli correlati come il silicio triplamente ionizzato (Si III, $\lambda1207$) sono mescolati all'interno della foresta.

Studi precedenti, in particolare McDonald et al. (2006, M06), hanno modellato la contaminazione da Si III assumendo che il contrasto del flusso trasmesso di Si III sia una versione linearmente riscalata del contrasto del flusso Ly$\alpha$. Questo ansatz descrive con successo le caratteristiche oscillanti nello spettro di potenza indotte dallo spostamento di velocità tra Ly$\alpha$ e Si III su grandi scale ($k \lesssim 0.02$ s km$^{-1}$). Tuttavia, questo articolo identifica che il modello M06 fallisce a scale più piccole ($k > 0.02$ s km$^{-1}$). Nello specifico, l'approccio M06 non riesce a tenere conto di:

1. Profili di Linea Distinti: Le linee di Si III sono termicamente più strette delle linee Ly$\alpha$ a causa della massa maggiore del silicio, portando a un aumento della potenza su piccola scala che il modello di riscalamento lineare non coglie.
2. Rapporti di Opacità Variabili: Il rapporto tra le opacità coeve di Si III e Ly$\alpha$ non è costante, ma dipende dalla densità del gas, dalla temperatura e dall'abbondanza di silicio.
3. Evoluzione con il Redshift: L'ampiezza della contaminazione da Si III potrebbe non seguire la scalatura precedentemente assunta con la trasmissione media di Ly$\alpha$, in particolare a causa dell'indurimento del fondo UV metagalattico per $z \lesssim 3.5$.

Metodologia
Gli autori utilizzano la suite di simulazioni idrodinamiche cosmologiche ad alta risoluzione Sherwood-Relics per costruire spettri di quasar fittizi.

• Configurazione della Simulazione: Le simulazioni utilizzano una scatola da $40 h^{-1}$ cMpc con $2 \times 2048^3$ particelle. Utilizzano un risolutore termochimico fuori equilibrio e un fondo UV spazialmente uniforme.
• Modellazione del Si III: A differenza di simulazioni precedenti che potrebbero mancare di modelli di produzione di metalli, gli autori assumono un modello di abbondanza di silicio fiduciale dipendente dalla sovradensità del gas ($\Delta$) e dal redshift ($z$), calibrato rispetto ai vincoli osservativi di Schaye et al. (2003) e Aguirre et al. (2004). La frazione di ionizzazione del Si III ($x_{\text{Si III}}$) è calcolata utilizzando il codice di foto-ionizzazione Cloudy combinato con lo spettro del fondo UV.
• Costruzione dello Spettro: Gli autori estraggono 5.000 linee di vista periodiche a intervalli di redshift $\Delta z = 0.1$. Calcolano le opacità Ly$\alpha$ e Si III, applicando convoluzioni di profilo Voigt (con le linee di Si III ricampionate su una griglia più fine per risolvere le loro larghezze termiche più strette).
• Analisi dello Spettro di Potenza: Calcolano lo spettro di potenza 1D ($P(k)$) per il flusso totale, per la componente solo Ly$\alpha$ e per quella solo Si III. Confrontano questi risultati diretti delle simulazioni contro il modello analitico M06 e la loro nuova funzione di fitting proposta.

Contributi Chiave e Sviluppo del Modello
L'articolo propone una nuova funzione di fitting analitica (Eq. 18) per sostituire l'ansatz M06 per la modellazione della contaminazione da Si III. Il nuovo modello introduce quattro parametri liberi per catturare due effetti fisici critici trascurati nel lavoro precedente:

1. Incremento Dipendente dalla Scala: Il modello include termini esponenziali ($\exp(k/k_{\text{auto}})$ e $\exp(k/k_{\text{cross}})$) per tenere conto dei profili di linea più stretti del Si III e del rapporto dipendente dai pixel delle opacità.
2. Decorrelazione: Il modello permette la decorrelazione tra l'assorbimento di Si III e Ly$\alpha$ su piccole scale. Ciò accade perché i forti assorbitori di Si III corrispondono spesso a linee Ly$\alpha$ saturate, rompendo la relazione lineare assunta in M06.

Lo spettro di potenza totale proposto è dato da:
$$P_{\text{tot}}(k) = P_{\text{Ly}\alpha}(k) + P_{\text{Si III}}(k) + P_{\text{Ly}\alpha-\text{Si III}}(k)$$
Dove:

• $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)$

Risultati

• Validazione: Il nuovo modello è stato validato contro le simulazioni Sherwood-Relics per $2.2 \le z \le 5.0$ e numeri d'onda $k < 0.2$ s km$^{-1}$.
• Performance: Alle grandi scale ($k < 0.02$ s km$^{-1}$), sia il modello M06 che il nuovo modello si adattano bene ai dati. Tuttavia, alle scale più piccole, il modello M06 sottostima significativamente la potenza, mentre il nuovo modello riproduce i dati della simulazione con una deviazione frazionaria inferiore all'1% per il caso fiduciale.
• Sensibilità dei Parametri:
  • I parametri di scala ($k_{\text{auto}}$ e $k_{\text{cross}}$) sono principalmente sensibili alla temperatura del gas ($T_0$). Gas più caldi portano a linee più larghe e diverse scale di decorrelazione.
  • I parametri di ampiezza ($a_{\text{auto}}$ e $a_{\text{cross}}$) sono principalmente sensibili all'abbondanza di silicio ($[\text{Si/H}]$).
• Evoluzione con il Redshift: Contrariamente all'assunzione M06 che la potenza del Si III scala inversamente con la trasmissione media ($1 - \bar{F}$), le simulazioni mostrano che l'ampiezza dello spettro di potenza del Si III aumenta verso redshift più bassi ($z \lesssim 3.5$). Ciò è guidato dall'indurimento del fondo UV, che riduce la frazione di ionizzazione del Si III ($x_{\text{Si III}}$) a redshift inferiori, ma è controbilanciato dall'evoluzione del modello di abbondanza di silicio utilizzato.

Significato e Rivendicazioni
Gli autori affermano che il loro modello migliorato è una considerazione essenziale per le future misurazioni di alta precisione dello spettro di potenza della foresta Ly$\alpha$, in particolare quelle che utilizzano spettrografi ad alta risoluzione (es. Keck/HIRES, VLT/UVES) che sondano scale oltre $k \sim 0.02$ s km$^{-1}$.

Per quanto riguarda i vincoli cosmologici, l'articolo afferma che l'adozione di un prior fisicamente motivato sull'abbondanza di silicio significa che il nuovo modello ha pochi effetti sui vincoli esistenti sulla materia oscura calda (WDM) derivati dalla foresta Ly$\alpha$. Tuttavia, gli autori sottolineano che ignorare l'incremento di potenza su piccola scala da parte del Si III potrebbe portare a errori sistematici nelle future analisi al diminuire le incertezze statistiche. Il lavoro non propone nuove campagne osservative, ma fornisce il quadro teorico necessario per interpretare correttamente i prossimi dati ad alta risoluzione.