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

Resumo Técnico: Um Modelo Melhorado para o Efeito da Absorção de Si III Correlacionada no Espectro de Potência de um Dímenso do Floresta Lyman-$\alpha$

Definição do Problema
O espectro de potência unidimensional da floresta Lyman-$\alpha$ (Ly$\alpha$) é uma ferramenta primária para restringir parâmetros cosmológicos e modelos de matéria escura. No entanto, a modelagem precisa é complicada pela contaminação de linhas de absorção de metais, especificamente aquelas correlacionadas com a floresta Ly$\alpha$. Enquanto metais não correlacionados podem ser estimados e subtraídos usando bandas laterais (sidebands), metais correlacionados como o silício triplamente ionizado (Si III, $\lambda1207$) estão misturados dentro da floresta.

Estudos anteriores, notadamente McDonald et al. (2006, M06), modelaram a contaminação de Si III assumindo que o contraste de fluxo transmitido de Si III é uma versão escalonada linearmente do contraste de fluxo de Ly$\alpha$. Este ansatz descreve com sucesso as características oscilatórias no espectro de potência induzidas pelo deslocamento de velocidade entre Ly$\alpha$ e Si III em grandes escalas ($k \lesssim 0.02$ s km$^{-1}$). No entanto, este artigo identifica que o modelo M06 falha em escalas menores ($k > 0.02$ s km$^{-1}$). Especificamente, a abordagem M06 falha ao não contabilizar:

1. Perfis de Linha Distintos: As linhas de Si III são termicamente mais estreitas do que as linhas de Ly$\alpha$ devido à massa mais elevada do silício, levando a um aumento do poder em pequenas escalas que o modelo de escalonamento linear não captura.
2. Razões de Profundidade Óptica Variáveis: A razão das profundidades ópticas de Si III e Ly$\alpha$ coevas não é constante, mas depende da densidade do gás, temperatura e abundância de silício.
3. Evolução com o Redshift: A amplitude da contaminação de Si III pode não seguir a escala anteriormente assumida com a transmissão média de Ly$\alpha$, particularmente devido ao endurecimento do fundo UV metagaláctico em $z \lesssim 3.5$.

Metodologia
Os autores utilizam a suíte de simulações hidrodinâmicas cosmológicas de alta resolução Sherwood-Relics para construir espectros de quasares fictícios (mock).

• Configuração da Simulação: As simulações utilizam uma caixa de $40 h^{-1}$ cMpc com $2 \times 2048^3$ partículas. Elas empregam um resolvedor de termoquímica de não-equilíbrio e um fundo UV espacialmente uniforme.
• Modelagem de Si III: Diferente de simulações anteriores que podem carecer de modelos de produção de metais, os autores assumem um modelo de abundância de silício fiduciar dependente da sobredensidade do gás ($\Delta$) e do redshift ($z$), calibrado contra restrições observacionais de Schaye et al. (2003) e Aguirre et al. (2004). A fração de ionização de Si III ($x_{\text{Si III}}$) é calculada usando o código de foto-ionização Cloudy combinado com o espectro do fundo UV.
• Construção do Espectro: Os autores extraem 5.000 linhas de visada periódicas em intervalos de redshift de $\Delta z = 0.1$. Eles computam as profundidades ópticas de Ly$\alpha$ e Si III, aplicando convoluções de perfil de Voigt (com as linhas de Si III reamostradas em uma grade mais fina para resolver suas larguras térmicas mais estreitas).
• Análise do Espectro de Potência: Eles computam o espectro de potência 1D ($P(k)$) para o fluxo total, apenas Ly$\alpha$, e apenas Si III. Eles comparam esses resultados diretos da simulação contra o modelo analítico M06 e a nova função de ajuste proposta.

Contribuições Principais e Desenvolvimento do Modelo
O artigo propõe uma nova função de ajuste analítica (Eq. 18) para substituir o ansatz M06 para modelar a contaminação de Si III. O novo modelo introduz quatro parâmetros livres para capturar dois efeitos físicos críticos negligenciados no trabalho anterior:

1. Aumento Dependente da Escala: O modelo inclui termos exponenciais ($\exp(k/k_{\text{auto}})$ e $\exp(k/k_{\text{cross}})$) para contabilizar os perfis de linha mais estreitos do Si III e a razão de profundidades ópticas dependente do pixel.
2. Descorrelação: O modelo permite a descorrelação entre a absorção de Si III e Ly$\alpha$ em pequenas escalas. Isso ocorre porque absorvedores fortes de Si III frequentemente correspondem a linhas de Ly$\alpha$ saturadas, quebrando a relação linear assumida em M06.

O espectro de potência total proposto é dado por:
$$P_{\text{tot}}(k) = P_{\text{Ly}\alpha}(k) + P_{\text{Si III}}(k) + P_{\text{Ly}\alpha-\text{Si III}}(k)$$
Onde:

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

Resultados

• Validação: O novo modelo foi validado contra as simulações Sherwood-Relics em redshifts $2.2 \le z \le 5.0$ e números de onda $k < 0.2$ s km$^{-1}$.
• Desempenho: Em grandes escalas ($k < 0.02$ s km$^{-1}$), tanto o modelo M06 quanto o novo modelo ajustam bem os dados. No entanto, em escalas menores, o modelo M06 subestima significativamente o poder, enquanto o novo modelo reproduz os dados da simulação com um desvio fracional de menos de 1% para o caso fiduciar.
• Sensibilidade dos Parâmetros:
  • Os parâmetros de escala ($k_{\text{auto}}$ e $k_{\text{cross}}$) são primariamente sensíveis à temperatura do gás ($T_0$). Gases mais quentes levam a linhas mais largas e diferentes escalas de descorrelação.
  • Os parâmetros de amplitude ($a_{\text{auto}}$ e $a_{\text{cross}}$) são primariamente sensíveis à abundância de silício ($[\text{Si/H}]$).
• Evolução com o Redshift: Contrariamente à suposição M06 de que o poder de Si III escala inversamente com a transmissão média ($1 - \bar{F}$), as simulações mostram que a amplitude do espectro de potência de Si III aumenta em direção a redshifts mais baixos ($z \lesssim 3.5$). Isso é impulsionado pelo endurecimento do fundo UV, que reduz a fração de ionização de Si III ($x_{\text{Si III}}$) em redshifts mais baixos, mas é contrabalançado pela evolução do modelo de abundância de silício utilizado.

Significância e Alegações
Os autores afirmam que o seu modelo melhorado é uma consideração essencial para futuras medições de alta precisão do espectro de potência da floresta Ly$\alpha$, particularmente aquelas utilizando espectrógrafos de alta resolução (ex: Keck/HIRES, VLT/UVES) que sondam escalas além de $k \sim 0.02$ s km$^{-1}$.

Em relação às restrições cosmológicas, o artigo afirma que, ao adotar um prior fisicamente motivado sobre a abundância de silício, o novo modelo tem pouco efeito nas restrições existentes de matéria escura morna (WDM) derivadas da floresta Ly$\alpha$. No entanto, os autores enfatizam que ignorar o aumento de poder em pequenas escalas proveniente do Si III poderia levar a erros sistemáticos em análises futuras à medida que as incertezas estatísticas diminuírem. O artigo não propõe novas campanhas observacionais, mas fornece o arcabouço teórico necessário para interpretar adequadamente os próximos dados de alta resolução.