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

技术摘要：一种改进的用于建模相关 Si III 吸收对一维 Lyman-$\alpha$ 森林功率谱影响的模型

问题陈述
一维 Lyman-$\alpha$ (Ly$\alpha$) 森林功率谱是约束宇宙学参数和暗物质模型的主要工具。然而，金属吸收线的污染使得准确建模变得复杂，特别是那些与 Ly$\alpha$ 森林相关的金属线。虽然不相关的金属可以通过侧带（sidebands）进行估计和减除，但像三价硅（Si III, $\lambda1207$）这样与 Ly$\alpha$ 森林混杂在一起的相关金属则难以处理。

先前的研究（特别是 McDonald 等人，2006，M06）通过假设 Si III 透射通量对比度是 Ly$\alpha$ 通量对比度的线性缩放版本来对 Si III 污染进行建模。这种假设成功描述了在大尺度（$k \lesssim 0.02$ s km$^{-1}$）上由 Ly$\alpha$ 与 Si III 之间的速度偏移所引起的功率谱振荡特征。然而，本文指出 M06 模型在小尺度（$k > 0.02$ s km$^{-1}$）上失效。具体而言，M06 方法未能考虑到以下几点：

1. 不同的谱线轮廓： 由于硅的质量较高，Si III 谱线在热力学上比 Ly$\alpha$ 谱线更窄，这导致了小尺度功率的增强，而线性缩放模型忽略了这一点。
2. 可变的光学深度比： 共时的 Si III 与 Ly$\alpha$ 光学深度之比并非定值，而是取决于气体密度、温度和硅丰度。
3. 红移演化： Si III 污染的振幅可能并不遵循先前假设的随平均 Ly$\alpha$ 透射率变化的比例关系，特别是由于在 $z \lesssim 3.5$ 时元银河系紫外背景（UVB）的硬化。

方法论
作者利用 Sherwood-Relics 系列高分辨率宇宙学流体动力学模拟来构建模拟类类星体光谱。

• 模拟设置： 模拟使用 $40 h^{-1}$ cMpc 的箱体，包含 $2 \times 2048^3$ 个粒子。他们采用了非平衡热化学求解器和空间均匀的紫外背景。
• Si III 建模： 与可能缺乏金属产生模型的先前模拟不同，作者假设了一个依赖于气体过密度（$\Delta$）和红移（$z$）的标称硅丰度模型，该模型根据 Schaye 等人 (2003) 和 Aguirre 等人 (2004) 的观测约束进行了校准。Si III 电离分数（$x_{\text{Si III}}$）是使用光电离代码 Cloudy 结合紫外背景光谱计算得出的。
• 光谱构建： 作者提取了在红移间隔 $\Delta z = 0.1$ 下的 5,000 条周期性视线。他们计算了 Ly$\alpha$ 和 Si III 的光学深度，并应用了 Voigt 轮廓卷积（将 Si III 谱线重采样到更细的网格上，以解析其更窄的热宽度）。
• 功率谱分析： 他们计算了总通量、仅 Ly$\alpha$ 成分以及仅 Si III 成分的一维功率谱（$P(k)$）。他们将这些直接模拟结果与解析的 M06 模型以及他们提出的新拟合函数进行了比较。

主要贡献与模型开发
论文提出了一种新的解析拟合函数（公式 18）来取代用于建模 Si III 污染的 M06 假设。新模型引入了四个自由参数，以捕捉前人工作中忽略的两个关键物理效应：

1. 尺度依赖性增强： 模型包含了指数项（$\exp(k/k_{\text{auto}})$ 和 $\exp(k/k_{\text{cross}})$），以解释更窄的 Si III 谱线轮廓和随像素变化的光学深度比。
2. 去相关性： 模型允许 Si III 与 Ly$\alpha$ 吸收在小尺度上的去相关。这是因为强 Si III 吸收体通常对应于饱和的 Ly$\alpha$ 线，从而打破了 M06 中假设的线性关系。

所提出的总功率谱为：
$$P_{\text{tot}}(k) = P_{\text{Ly}\alpha}(k) + P_{\text{Si III}}(k) + P_{\text{Ly}\alpha-\text{Si III}}(k)$$
其中：

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

结果

• 验证： 新模型在 $2.2 \le z \le 5.0$ 的红移范围和 $k < 0.2$ s km$^{-1}$ 的波数范围内，通过 Sherwood-Relics 模拟得到了验证。
• 性能： 在大尺度（$k < 0.02$ s km$^{-1}$）上，M06 模型和新模型都能很好地拟合数据。然而，在小尺度上，M06 模型显著低估了功率，而新模型在标称情况下与模拟数据的偏差小于 1%。
• 参数敏感性：
  • 尺度参数（$k_{\text{auto}}$ 和 $k_{\text{cross}}$）主要受气体温度（$T_0$）的影响。更热的气体会导致更宽的谱线和不同的去相关尺度。
  • 幅度参数（$a_{\text{auto}}$ 和 $a_{\text{cross}}$）主要受硅丰度（$[\text{Si/H}]$）的影响。
• 红移演化： 与 M06 认为 Si III 功率随平均透射率（$1 - \bar{F}$）反向变化的假设相反，模拟显示 Si III 功率谱振幅在向低红移（$z \lesssim 3.5$）移动时趋于增加。这是由于紫外背景的硬化降低了 Si III 的电离分数（$x_{\text{Si III}}$），但被所使用的硅丰度模型的演化所抵消。

意义与主张
作者指出，对于未来高精度的 Ly$\alpha$ 森林功率谱测量（特别是使用高分辨率光谱仪，如 Keck/HIRES, VLT/UVES，这些仪器探测的尺度超过 $k \sim 0.02$ s km$^{-1}$），改进后的模型是一个必不可少的考虑因素。

关于宇宙学约束，论文声称，通过采用基于物理动机的硅丰度先验，新模型对现有的由 Ly$\alpha$ 森林推导出的温暗物质（WDM）约束几乎没有影响。然而，作者强调，如果忽略来自 Si III 的小尺度功率增强，随着统计不确定性的减小，可能会导致未来的分析出现系统误差。本文并未提出新的观测活动，而是为准确解释即将到来的高分辨率数据提供了必要的理论框架。