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

Resumen Técnico: Un Modelo Mejorado para el Efecto de la Absorción de Si III Correlacionada en el Espectro de Potencia de la Foresta de Lyman-$\alpha$ Unidimensional

Planteamiento del Problema
El espectro de potencia unidimensional de la foresta de Lyman-$\alpha$ (Ly$\alpha$) es una herramienta primordial para restringir los parámetros cosmológicos y los modelos de materia oscura. Sin embargo, su modelado preciso se complica por la contaminación de líneas de absorción metálicas, específicamente aquellas correlacionadas con la foresta de Ly$\alpha$. Los metales no correlacionados pueden estimarse y restarse utilizando bandas laterales (sidebands), pero los metales correlacionados como el silicio triplemente ionizado (Si III, $\lambda1207$) se encuentran mezclados dentro de la foresta.

Estudios previos, notablemente McDonald et al. (2006, M06), modelaron la contaminación de Si III asumiendo que el contraste de flujo transmitido de Si III es una versión reescalada linealmente del contraste de flujo de Ly$\alpha$. Este ansatz describe con éxito las características oscilatorias en el espectro de potencia inducidas por el desfase de velocidad entre Ly$\alpha$ y Si III en grandes escalas ($k \lesssim 0.02$ s km$^{-1}$). Sin embargo, este artículo identifica que el modelo M06 falla en escalas más pequeñas ($k > 0.02$ s km$^{-1}$). Específicamente, el enfoque de M06 no logra contabilizar:

1. Perfiles de Línea Distintos: Las líneas de Si III son térmicamente más estrechas que las líneas de Ly$\alpha$ debido a la mayor masa del silicio, lo que conduce a un aumento de la potencia en escalas pequeñas que el modelo de reescalado lineal no logra capturar.
2. Razones de Profundidad Óptica Variables: La razón entre las profundidades ópticas de Si III y Ly$\alpha$ coetáneas no es constante, sino que depende de la densidad del gas, la temperatura y la abundancia de silicio.
3. Evolución con el Desplazamiento al Rojo (Redshift): La amplitud de la contaminación de Si III puede no seguir la escala previamente asumida con la transmisión media de Ly$\alpha$, particularmente debido al endurecimiento del fondo UV metagaláctico en $z \lesssim 3.5$.

Metodología
Los autores utilizan la suite de simulaciones hidrodinámicas cosmológicas de alta resolución Sherwood-Relics para construir espectros de cuásares simulados (mock).

• Configuración de la Simulación: Las simulaciones utilizan un volumen de $40 h^{-1}$ cMpc con $2 \times 2048^3$ partículas. Emplean un resolvedor de termoquímica de no equilibrio y un fondo UV espacialmente uniforme.
• Modelado de Si III: A diferencia de simulaciones previas que pueden carecer de modelos de producción de metales, los autores asumen un modelo de abundancia de silicio fiduciario dependiente de la sobredensidad del gas ($\Delta$) y del desplazamiento al rojo ($z$), calibrado con las restricciones observacionales de Schaye et al. (2003) y Aguirre et al. (2004). La fracción de ionización de Si III ($x_{\text{Si III}}$) se calcula utilizando el código de fotoionización Cloudy combinado con el espectro del fondo UV.
• Construcción de Espectros: Los autores extraen 5,000 líneas de visión periódicas en intervalos de desplazamiento al rojo de $\Delta z = 0.1$. Computan las profundidades ópticas de Ly$\alpha$ y Si III, aplicando convoluciones de perfil de Voigt (con las líneas de Si III remuestreadas en una rejilla más fina para resolver sus anchos térmicos más estrechos).
• Análisis del Espectro de Potencia: Computan el espectro de potencia 1D ($P(k)$) para el flujo total, y para los componentes de solo Ly$\alpha$ y solo Si III. Comparan estos resultados directos de la simulación contra el modelo analítico de M06 y su nueva función de ajuste propuesta.

Contribuciones Clave y Desarrollo del Modelo
El artículo propone una nueva función de ajuste analítica (Ec. 18) para reemplazar el ansatz de M06 para modelar la contaminación de Si III. El nuevo modelo introduce cuatro parámetros libres para capturar dos efectos físicos críticos que fueron omitidos en el trabajo anterior:

1. Aumento Dependiente de la Escala: El modelo incluye términos exponenciales ($\exp(k/k_{\text{auto}})$ y $\exp(k/k_{\text{cross}})$) para dar cuenta de los perfiles de línea más estrechos de Si III y de la relación de profundidades ópticas dependiente de los píxeles.
2. Descorrelación: El modelo permite la descorrelación entre la absorción de Si III y Ly$\alpha$ en escalas pequeñas. Esto ocurre porque los absorbedores fuertes de Si III suelen corresponder a líneas de Ly$\alpha$ saturadas, rompiendo la relación lineal asumida en M06.

El espectro de potencia total propuesto viene dado por:
$$P_{\text{tot}}(k) = P_{\text{Ly}\alpha}(k) + P_{\text{Si III}}(k) + P_{\text{Ly}\alpha-\text{Si III}}(k)$$
Donde:

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

• Validación: El nuevo modelo fue validado contra las simulaciones Sherwood-Relics en el rango de desplazamientos al rojo $2.2 \le z \le 5.0$ y números de onda $k < 0.2$ s km$^{-1}$.
• Desempeño: En grandes escalas ($k < 0.02$ s km$^{-1}$), tanto el modelo M06 como el nuevo modelo ajustan bien los datos. Sin embargo, en escalas más pequeñas, el modelo M06 subestima significativamente la potencia, mientras que el nuevo modelo reproduce los datos de la simulación con una desviación fraccional de menos del 1% para el caso fiduciario.
• Sensibilidad de los Parámetros:
  • Los parámetros de escala ($k_{\text{auto}}$ y $k_{\text{cross}}$) son principalmente sensibles a la temperatura del gas ($T_0$). Un gas más caliente produce líneas más anchas y diferentes escalas de descorrelación.
  • Los parámetros de amplitud ($a_{\text{auto}}$ y $a_{\text{cross}}$) son principalmente sensibles a la abundancia de silicio ($[\text{Si/H}]$).
• Evolución con el Desplazamiento al Rojo: Contrario a la suposición de M06 de que la potencia de Si III escala inversamente con la transmisión media ($1 - \bar{F}$), las simulaciones muestran que la amplitud del espectro de potencia de Si III aumenta hacia desplazamientos al rojo menores ($z \lesssim 3.5$). Esto es impulsado por el endurecimiento del fondo UV, que reduce la fracción de ionización de Si III ($x_{\text{Si III}}$) a desplazamientos al rojo menores, pero es contrarrestado por la evolución del modelo de abundancia de silicio utilizado.

Significancia y Reivindicaciones
Los autores afirman que su modelo mejorado es una consideración esencial para las futuras mediciones de alta precisión del espectro de potencia de la foresta de Ly$\alpha$, particularmente aquellas que utilizan espectrógrafos de alta resolución (ej. Keck/HIRES, VLT/UVES) que sondean escalas más allá de $k \sim 0.02$ s km$^{-1}$.

Respecto a las restricciones cosmológicas, el artículo sostiene que adoptar un prior motivado físicamente sobre la abundancia de silicio significa que el nuevo modelo tiene poco efecto en las restricciones existentes de materia oscura tibia (WDM) derivadas de la foresta de Ly$\alpha$. Sin embargo, los autores enfatizan que ignorar el aumento de potencia en escalas pequeñas podría conducir a errores sistemáticos en futuros análisis a medida que las incertidumbres estadísticas disminuyan. El artículo no propone nuevas campañas observacionales, sino que proporciona el marco teórico necesario para interpretar con precisión los próximos datos de alta resolución.