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$フォレストと相関している金属（例：三価のケイ素、Si III、$\lambda1207$）は、フォレスト内に混在している。

先行研究（特にMcDonald et al. 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$よりも狭くなっており、これがM06の線形スケーリングモデルが見落としている小スケールのパワーの増強をもたらしている。
2. 変動する光学的厚さの比： 同時期のSi IIIとLy$\alpha$の光学的厚さの比は一定ではなく、ガスの密度、温度、およびケイ素存在量に依存する。
3. 赤方偏移進化： Si IIIの汚染の振幅は、以前に想定されていた平均Ly$\alpha$透過率とのスケーリング（$1 - \bar{F}$）に従わない可能性がある。これは、$z \lesssim 3.5$におけるメタ銀河UV背景放射の硬化によるものである。

手法
著者らは、高解像度の宇宙論的流体シミュレーションであるSherwood-Relicsスイートを用いて、模擬クエーサースペクトルを構築している。

• シミュレーション設定： シミュレーションは、$40 h^{-1}$ cMpcのボックスを用い、$2 \times 2048^3$個の粒子を使用している。非平衡熱化学ソルバーと空間的に一様なUV背景を採用している。
• Si IIIモデリング： 金属生成モデルを欠いている可能性がある従来のシミュレーションとは異なり、著者らはSchaye et al. (2003) および Aguirre et al. (2004) の観測的制約に基づいて較正された、ガス過密度（$\Delta$）と赤方偏移（$z$）に依存する標準的なケイ素存在量モデルを仮定している。Si IIIの電離分率（$x_{\text{Si III}}$）は、UV背景スペクトルと組み合わせた光電離コード Cloudy を用いて計算されている。
• スペクトル構築： 著者らは、赤方偏移間隔 $\Delta z = 0.1$ で5,000本の周期的な視線方向を抽出している。Ly$\alpha$およびSi IIIの光学的厚さを計算し、Voigtプロファイルの畳み込みを適用している（Si IIIのラインは、そのより狭い熱的幅を解像するために、より細かいグリッド上に再サンプリングされている）。
• パワースペクトル解析： 全フラックス、Ly$\alpha$のみ、およびSi IIIのみの成分について、一次元パワースペクトル（$P(k)$）を計算している。これら直接的なシミュレーション結果を、M06の解析モデルおよび提案された新しいフィッティング関数と比較している。

主要な貢献とモデル開発
本論文は、Si III汚染をモデル化するためのM06の仮定に代わる、新しい解析的フィッティング関数（式18）を提案している。この新モデルは、先行研究で無視されていた2つの重要な物理的効果を捉えるために、4つの自由パラメータを導入している。

1. スケール依存の増強： モデルには、Si IIIのより狭いラインプロファイルとピクセル依存の光学的厚さの比を考慮するために、指数関数項（$\exp(k/k_{\text{auto}})$ および $\exp(k/k_{\text{cross}})$）が含まれている。
2. 脱相関（Decorrelation）： モデルは、小スケールにおける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}]$）に敏感である。
• 赤方偏移進化： Si IIIのパワーは平均透過率（$1 - \bar{F}$）に反比例するというM06の仮定に反して、シミュレーションはSi IIIのパワースペクトルの振幅が低い赤方偏移（$z \lesssim 3.5$）に向かって増加することを示している。これは、UV背景の硬化がSi IIIの電離分率（$x_{\text{Si III}}$）を減少させるものの、使用されたケイ素存在量モデルの進化によってそれが相殺されるためである。

意義と主張
著者らは、改良されたモデルが、将来の高精度なLyman-$\alpha$フォレスト・パワースペクトル測定（特に $k \sim 0.02$ s km$^{-1}$ 以上のスケールを探索する Keck/HIRES や VLT/UVES などの高解像度分光器を用いるもの）において不可欠な検討事項であると述べている。

宇宙論的制約に関しては、ケイ素存在量に関する物理的に動機付けられた事前分布を採用することで、新しいモデルは、Ly$\alpha$フォレストから導出される既存の温いダークマター（WDM）制約にはほとんど影響を与えないと主張している。しかし、著者らは、Si IIIによる小スケールのパワー増強を無視することは、統計的な不確かさが減少するにつれて、将来の解析において系統誤差を招く可能性があることを強調している。本論文は新しい観測キャンペーンを提案するものではなく、今後の高解像度データを正確に解釈するための必要な理論的枠組みを提供するものである。