Generation of Lognormal Synthetic Lyman-$α$ Forest Spectra for $P_{1D}$ Analysis

技術要約：P1D解析のための対数正規合成ライマンα森林スペクトルの生成

問題提起
ライマンα（Lyα）森林の一次元フラックスパワースペクトル（$P_{1D}$）は、宇宙再電離期における銀河間物質（IGM）の小スケール構造を調べるための重要なプローブであり、宇宙論パラメータ（例：ニュートリノ質量、ダークマターの性質）や天体物理学的プロセス（例：IGMの熱的履歴）に対して高い感度を持つ。しかし、これらのパラメータを抽出するには、観測機器の影響、連続成分（コンティニュアム）フィッティングの不確実性、および天体物理学的な汚染物質を含む系統誤差に対する厳密な検証が必要である。流体計算シミュレーションは高精度なIGMの実現を提供できるが、DESI（Dark Energy Spectroscopic Instrument）のような大量のデータを扱うサーベイにおいて、堅牢な誤差推定やパイプラインの検証に必要な大規模なスペクトルアンサンブルを生成するには、計算コストが極めて高い。従来の対数正規（lognormal）によるモック生成手法は、より高速ではあるものの、多くの場合、基礎となるガウス型パワースペクトルの静的な解析的形式に依存しており、現代のサーベイが持つ広範なパラメータ空間における観測された$P_{1D}$のスケール依存の形状や赤方偏移進化を正確に再現するための柔軟性に欠けていた。

手法
著者らは、$P_{1D}$解析に特化した合成1次元Lyα森林スペクトルを生成するための、改良された対数正規モックフレームワークを提示する。本手法はMcDonald et al. (2006) および Karaçaylı et al. (2020) による定式化に基づいているが、忠実度を高めるために2つの主要な革新を導入している。

1. 赤方偏移依存のガウス相関反転： 入力されるガウス型パワースペクトル（$P_G$）に固定された解析的形式を利用していた従来のアプローチとは異なり、本手法では構成空間における最適なガウス相関関数（$\xi_G$）を直接解く。フルシーケンスの対数正規変換を適用した後のフラックス相関関数（$\xi_F$）がターゲットとする$P_{1D}$と一致するように$\xi_G$を反復的に調整することで、一意で赤方偏移依存のガウス型パワースペクトル$P_G(k, z)$を導出する。この反転は、変換後の統計量とターゲットモデル（具体的にはDESI Early Data Releaseの測定値）の間の残差を最小化する非線形最小二乗ソルバーを用いて、各点ごとに実行される。
2. 解析的な平均フラックス・フィッティング： 本フレームワークは、平均透過フラックスの進化を個別に最適化する。著者らは有効光学的厚さモデル（$\tau(z) = \tau_0(1+z)^\gamma$）を採用し、DESI Data Release 1 (DR1) で観測された平均フラックスの進化を再現するために、4つの変換パラメータ（$\tau_0, \tau_1, \nu, \sigma^2_F$）をフィッティングする。これは、対数正規変換から導かれる解析的な平均フラックスとターゲットとなる観測データとの間の重み付き最小二乗コスト関数を最小化することによって達成される。

生成パイプラインは、1次元のランダムなガウス速度グリッドを作成し、導出された赤方偏移依存の$P_G(k, z)$を適用し、赤方偏移因子$a(z)$で場を進化させ、非ガウス性を近似するために二乗対数正規変換を適用した後、光学的厚さと透過フラックスへと変換する工程を含む。

主な結果
著者らは、DESI EDRの赤方偏移範囲（$2.0 \le z \le 3.8$）にわたって1,000個の独立した視線方向の実現例を生成することで、フレームワークの検証を行った。結果は以下の通りである。

• 平均フラックスの精度： モックは平均透過フラックスの進化をサブパーセントの精度で再現している。DESI EDRの全赤方偏移範囲において、平方根平均二乗（RMS）相対誤差は0.002であり、赤方偏移境界における最大偏差は0.5%であった。
• パワースペクトルの回収： 本手法は、$P_{1D}$の形状と振幅を高精度に再現する。DESI EDRの$k$範囲（$0.00075 \le k \lesssim 0.035 \text{ s km}^{-1}$）において、平均RMS相対誤差は0.02である。これは、同範囲において0.10の誤差を生じた従来のメソッド（Karaçaylı et al. 2020）と比較して、大幅な改善を示している。
• スケーラビリティ： 本フレームワークは、広範なスケール（$10^{-5} \le k \lesssim 0.1 \text{ s km}^{-1}$）においても堅牢な性能を維持しており、この広い範囲でのRMS相対誤差は0.11にとどまっている。これは、以前の手法による誤差18.15と比較して劇的な改善である。本手法におけるこの広い範囲での最大偏差は40%未満であり、従来のアプローチの300%と比較して極めて低い。
• 収束性： 収束解析により、1,000個の実現例があればRMS相対誤差が漸近的に一定値に近づくことが示されており、検証テストのための統計的安定性が確保されている。

意義と主張
本論文は、この改良された対数正規フレームワークが、大規模な合成スペクトルを生成するための、計算効率が高く解析的に扱いやすい、流体計算シミュレーションに代わる選択肢を提供すると主張している。平均フラックスの振幅と詳細なフラックス相関のモデリングを分離し、入力となるガウス型パワースペクトルを直接解くことで、本手法はターゲットとなる観測量に対する優れた制御を実現している。

著者らは、これらのモックは第一原理からの物理をモデル化していないため、完全な宇宙論的推論を目的としたものではないことを強調しているが、以下の用途には非常に適している。

• 解析パイプラインの検証および系統誤差（ノイズ、分解能、連続成分フィッティングなど）の定量化。
• FFTおよび尤度ベースのエスティメータを用いた、現在のおよび将来のDESIデータにおける$P_{1D}$測定のサポート。
• quickquasarsなどのツールと組み合わせた、天体物理学的な汚染物質（Damped Ly$\alpha$系、広幅吸収線など）の制御された研究。

本研究は、エミュレータやシミュレーションの出力が対数正規モック生成の直接的な入力として機能するエンドツーエンド解析の基礎を確立しており、理論モデルと観測的推論の架け橋となるものである。著者らは、分散への依存関係をより適切に扱うために、平均フラックスと$P_{1D}$を共同でフィッティングするなどの今後の改良により、本フレームワークの有用性をさらに高めることができると述べている。