Constraining Dark Acoustic Oscillations with the High-Redshift UV Luminosity Function

技術要約：高赤方偏移UV光度関数によるダーク音響振動の制約

問題提起
標準的な$\Lambda$CDMモデルは、ダークマター（DM）が冷たく衝突のないものであると仮定しているが、ダークセクターの様々な動機付けられた拡張モデルは「ダーク音響振動（Dark Acoustic Oscillations: DAOs）」を予測している。これらは、宇宙初期において一部のダークマターがダーク放射（DR）種族と相互作用することで発生し、物質パワースペクトル $P(k)$ の小スケールにおける抑制と減衰振動をもたらす。既存の宇宙マイクロ波背景放射（CMB）による制約は、大スケールでは厳格であるが、強いDAOを持つ場合の相互作用するダークマターの割合については $\lesssim 5\%$ に制限されている。一方、Lyman-$\alpha$ フォレストによる制約は、弱いDAOや完全に相互作用するモデルに対しては有効であるが、強いDAOのシグネチャーと温いダークマター（WDM）の効果を区別することに課題がある。したがって、非線形構造形成と、それがハロー質量関数（HMF）に与える具体的な刻印を調査できる、高赤方偏移（$z \sim 3-9$）における小スケールの構造に対するモデルに依存しない制約が必要とされている。

手法
著者らは、高赤方偏移の銀河の紫外線光度関数（UVLF）を用いてDAOパラメータを制約するパイプラインを開発した。その手法は、以下の4つのコンポーネッツで構成される：

1. 現象論的転送関数モデル：
   マルコフ連鎖モンテカルロ（MCMC）スキャンにおけるボルツマン方程式の計算コストを回避するため、著者らは$\Lambda$CDMに対する新しい4パラメータの現象論的転送関数 $T(k)$ を導入した。このモデルは、DAOによる抑制と減衰振動を捉え、1未満の相互作用するダークマターの割合（$f$）を許容する。パラメータは以下の通りである：

   • $f$: 振動を起こすダークマターの割合
   • $k_{\text{peak}}$: 第1次DAOピークの波数
   • $k_{\text{damp}}$: 振動の減衰スケール
   • $A$: 振動の振幅
     このモデルは、修正ボルツマンコード CLASS-aDM（原子ダークマター用）および CLASS（相互作用ダークマター用）によって生成されたパワースペクトルに対して較正されており、これらのパラメータがダーク音響地平線や拡散減衰スケールといった物理量に直接対応することを示している。

2. 拡張プレス・シェクター（EPS）定式化とN体シミュレーションによる較正：
   線形パワースペクトルからHMFを予測するために、著者らは拡張プレス・シェクター定式化を用いる。標準的なウィンドウ関数（例：実空間トップハット）は、切断されたパワースペクトルに対して機能しないことを認識し、「smooth-k」ウィンドウ関数を利用している。極めて重要な点として、著者らはMUSIC初期条件を用いたGIZMOコードによる13組のダークマターのみのN体シミュレーションを実施した。これらのシミュレーションは、様々な相互作用するダークマターの割合（$f=0.1, 0.5, 1$）およびDAOスケール（$k_{\text{DAO}} = 5, 20, 66, 190 \, h/\text{Mpc}$）をカバーしている。これらのシミュレーションは、赤方偏移 $z=4$ から $12$において半解析的なHMFがシミュレーション結果を正確に再現するように、EPSパラメータ（$q, \beta, c$）を較正するために使用される。

3. ハロー・銀河接続：
   HMFと観測可能なUVLFとの関係は、GALLUMIコードを用いてモデル化されている。このパイプラインは、UV光度が星形成率（SFR）によって駆動され、そのSFRは平均的にハロー降着率に対して二重べき乗則に従うと仮定している。モデルには、SFRとハロー質量の関係（正規化、ピボットスケール、進化）およびマグニチュードの確率性（stochasticity）を記述する8つの天体物理学的ノイズパラメータが含まれている。これらのパラメータは、高赤方偏移の銀河形成物理学に伴う不確実性を考慮するために周辺化される。

4. データ解析：
   著者らは、赤方偏移 $z=3$ から $9$をカバーする、ハッブル宇宙望遠鏡（HST、ブランクフィールドおよびレンズ場の両方）、ジェイムズ・ウェッブ宇宙望遠鏡（JWST）、スバル望遠鏡、およびカナダ・フランス・ハワイ望遠鏡（CFHT）からのUVLF測定値を組み合わせている。彼らは、$\Lambda$CDMパラメータをPlanck 2018の最良適合値に固定した上で、ベイズMCMCスキャンおよび頻度論的プロファイル尤度解析を実行している。

主な貢献

• 新しい現象論的モデル： 特定の粒子物理学モデル（原子ダークマターなど）と一般的な制約との間のギャップを埋める、サブドミナントな相互作用ダークマターの割合も受け入れ可能な、柔軟な転送関数モデルの導入。
• シミュレーションによる較正： 減衰振動が存在する状況下でのHMFを予測するための、専用のN体シミュレーションを用いたEPS定式化の初の較正。これにより、半解析的なアプローチの妥当性を検証した。
• モデルに依存しない制約： 特定のダークセクターの粒子物理学に依存しない、DAOパラメータ空間（$f$ 対 $k_{\text{peak}}$）の制約の導出。これにより、様々なモデル（IDM、原子ダークマターなど）への容易な再解釈が可能となる。

結果

• $f$ と $k_{\text{peak}}$ に対する制約： 解析により、相互作用するダークマターの割合とDAOのスケールの間に強い相関があることが示された。
  • もし相互作用するダークマターの割合が有意な場合（$f \gtrsim 0.07$）、UVLFのデータと整合するためには、第1次DAOピークの波数は $k_{\text{peak}} \gtrsim 50 \, h/\text{Mpc}$ でなければならない。
  • 逆に、$k_{\text{peak}}$ が小さい場合（デカップリングがより早い、あるいは相互作用がより強いことを示す）、相互作用するダークマターの割合は $f \lesssim 0.07$ に制約される。
• 他のプローブとの比較：
  • 高い割合の相互作用するダークマター（$f \sim 1$）の場合、UVLFによる制約は、CMBから導かれる境界よりも $k_{\text{peak}}$ に関して約1桁向上している。
  • 非常に小さなスケール（$k_{\text{peak}} \lesssim 1 \, h/\text{Mpc}$）においては、CMBによる制約が、相互作用するダークマターの存在量（$f \lesssim 5\%$）に対して依然としてより厳しい。
  • 相互作用ダークマター（IDM）モデルにおけるダークマターとバリオンの散乱断面積に対する制約は、CMBによる制約よりも5桁以上強く、天の川銀河の衛星銀河の存在量から得られる境界とも競合するレベルである。
• 制約されないパラメータ： 他のパラメータを周辺化した際、DAOの振幅（$A$）または減衰スケールの比（$k_{\text{damp}}/k_{\text{peak}}$）については、有意な制約は見られなかった。これは、UVLFの感度が、振動の特徴そのものよりも、むしろパワースペクトルの初期の抑制によって主に駆動されているためである。

意義
本論文は、高赤方偏移のUV光度関数測定が、小スケールのダークマター物理学を制約するための、CMBやLyman-$\alpha$ フォレストに対して強力かつ補完的なプローブであることを示している。専用のN体シミュレーションを用いてEPS定式化を校正し、柔軟なハロー・銀河接続モデルを採用することで、著者らは、音響振動を起こすダークマターの割合に対する堅牢な限界を確立した。これらの結果は、相互作用するダークマターのデカップリング時期をより早い時間（より高い $k_{\text{peak}}$）へと押し上げるか、あるいは特定の相互作用スケールに対して、相互作用する割合をサブドミナント（$f < 7\%$）にすることを要求している。本研究は、異なるダークセクターのパラダイムを区別するための高赤方偏移観測量の有用性を強調しており、ベラ・ルービン天文台、ナンシー・グレース・ローマン宇宙望遠鏡、およびスクエア・キロメートル・アレイによる今後のデータを用いた制約のための基礎を築くものである。