The Manticore Project II: Bayesian digital twins of cosmic structure across the SDSS and BOSS volumes

技術要約：マンティコア・プロジェクト II：SDSSおよびBOSSボリュームにおけるベイズ的デジタルツインによる宇宙構造の再構成

問題と動機
宇宙の大規模構造（LSS）に関する従来の解析は、二点相関関数やパワースペクトルといった低次の要約統計量に依存している。これらの手法は効果的ではあるが、観測データを圧縮してしまい、進化した物質場に刻み込まれた高次の相関や位相情報を破棄してしまう。フィールドレベルの推論は、根本的に異なるアプローチを提供し、観測された銀河分布を生じさせた完全な三次元の密度場および速度場の回復を目指すものである。一方、ベイズ的なフィールドレベルの手法（例：BORG）は、近傍の宇宙（$\lesssim 200 \, h^{-1}\text{Mpc}$）には成功裏に適用されてきたが、計算上の制約や不均質な深層サーベイデータの処理の必要性から、宇宙論的に代表的な体積（$\sim \text{Gpc}^3$）へとこの高解像度かつ完全なベイズ再構成を拡張することは、依然として未解決の課題であった。

手法
本論文では、赤方偏移 $z \approx 0.7$ まで及ぶ共動体積 $(4 \, h^{-1}\text{Gpc})^3$、解像度 $\sim 4 \, h^{-1}\text{Mpc}$ の高解像度ベイズ・フィールドレベル推論である Manticore-Deep を提示する。この手法は、コンパニオン論文である「Manticore-Local」解析で確立されたフレームワークに基づいているが、いくつかの重要な拡張を導入している。

• マルチサーベイ統合： 単一の階層的ベイズ・フレームワーク（BORGアルゴリズム）内で、5つの銀河レッドシフトサーベイ（2M++, 6dFGS, 2dFGRS, SDSS, BOSS）を共同で制約する。
• タイル化推論戦略： $1024^3$ ボクセルの推論を扱いやすくするため、親ドメインを64個の立方体タイル（各辺 $1024 \, h^{-1}\text{Mpc}$）からなる $4 \times 4 \times 4$ グリッドに分割する。各タイルは $256^3$ グリッド上で独立して推論される。境界部のアーティファクトを軽減し連続性を確保するために、隣接するタイル間で銀河を複製するバッファ領域（$64 \, h^{-1}\text{Mpc}$）を採用し、外側のバッファ領域は純粋に $\Lambda$CDM プライアの下で進化させる。サンプリング後、バッファ領域は破棄され、信号を運ぶコア部分が結合されてグローバルな場を形成する。
• フォワードモデル： 本モデルは、原始的なガウス型ホワイトノイズ初期条件を、非線形COLAベースの重力ソルバー、赤方偏移空間歪みのモデリング、および局所的な非線形銀河バイアス（Sigmoid-Truncated Double Power-Law）を介して、現在の観測量へと写像する。尤度は、過分散を考慮するために一般化ポアソン分布を利用する。
• 事後再シミュレーション： 推論された初期条件は、高精度な $N$ 体再シミュレーション（SWIFTソルバー、$4096^3$ ダークマター粒子）を用いて進化させ、ハローカタログ、マージツリー、およびライトコーン出力を生成する。このステップは、推論の検証を行い、より高解像度の構造へのアクセスを提供する。
• 検証フレームワーク： 物理的な忠実度は、推論に使用されていない2つの独立したテンプレートフリーの事後予測可能量に対してテストされる：
  1. CMBレンシング相互相関： 再構成された物質密度を、Planck PR3 CMBレンシング収束マップと比較する。
  2. 運動学的スニヤエフ・ゼルドビッチ（kSZ）スタッキング： Planck 217 GHzマップ上での銀河クラスターの速度重み付きスタッキングにより、3D速度場をテストする。

主な結果

• 統計的一致性： 事後実現（posterior realisations）は $\Lambda$CDM と統計的に一致している。推論された初期条件は、タイル化戦略による方向性のアーティファクトを検出することなく、ガウス的かつ等方的なパワースペクトルを示す。進化した物質パワースペクトル、バイスペクトル、および $z=0$ でのハロー質量関数は、制御用の $\Lambda$CDM シミュレーションとすべての分解能スケールにおいて一致している。
• CMBレンシングの検出： 再構成された物質場と Planck CMB レンシングマップとの相互相関により、保守的な累積検出有意度 7.4 $\sigma$ が得られた。これは、銀河バイアスやテンプレートフィッティングを仮定することなく、再構成された3D密度場が観測された宇宙と同様にCMBを整合的にレンズしていることを裏付けている。
• kSZ検出： 64,750個の銀河クラスターを用いた速度重み付きスタッキングにより、中央値で 3.5 $\sigma$（事後アンサンブル全体で 2.9 から 4.3 \sigma} の範囲）の運動学的スニヤエフ・ゼルドビッチ検出が得られた。後退・接近によるクラスターのモデルに依存しない分割は、推論された速度が真のクラスター運動と統計的に一致していることを確認している。
• BOSS Great Wall ケーススタディ： 再構成によって、BOSS Great Wall が一貫した過密領域として回収された。シリンダー質量解析により、ランダムな視線に対する開口半径 $\sim 100 \, h^{-1}\text{Mpc}$ における $\sim 3\sigma$ の過密が明らかになり、このような巨大構造が観測データに基づいた条件下で $\Lambda$CDM の期待値と整合していることを示している。

意義と主張
本論文は、Manticore-Deep を、サーベイ深度における制約付き宇宙論的デジタルツインの新たなベンチマークとして確立する。その主要な意義は以下の通りである：

1. スケーラビリティ： 新しいタイル化戦略を通じて、高解像度の完全なベイズ・フィールドレベル推論を、局所的な体積からマルチギガパーセックスケールへと拡張可能であることを示した。これにより、物理的な忠実度と統計的一致性が維持されている。
2. フィールドレベルの検証： 投影された密度と3D速度場の両方を検証する、厳格なテンプレートフリーのフィールドレベル診断（CMBレンシング相互相関およびkSZスタッキング）を導入し、適用した。これらの指標は、宇宙論的な深さにおいて利用可能なオブジェクトレベルの検証が不足している現状に対し、再現可能なクロスメソッド・ベンチマークを提供する。
3. 物理的忠実度： 観測された大規模構造に対する因果的かつ物理学に基づいた説明を提供し、密度と速度の両方の情報を回収し、独立した観測プローブによって検証されている。

著者らは、これらの結果がタイル化推論アプローチの堅牢性を支持しており、将来の制約付きシミュレーションのための標準的なコミュニティ・ベンチマークとして、これらのフィールドレベルの指標を採用することを推奨すると結論づけている。