The Back-in-time Void Finder: dynamical identification of cosmic voids through optimal transport reconstruction

技術要約：Back-in-time Void Finder (BitVF)

問題提起
宇宙空洞（コスミック・ボイド）は強力な宇宙論的プローブであるが、その同定は重大な観測的および手法的な課題によって妨げられている。従来のトポロジカルな手法や密度ベースのボイド検出器は、トレーサーのサンプリングが希薄であり、銀河バイアスが高まる低密度領域において、ショットノイズに対して非常に敏感である。さらに、赤方偏移空間における同定は、赤方偏移空間歪み（RSD）およびアルコック・パチニ（AP）効果によって損なわれる。これらはボイドの形状を歪ませ、視線方向（LOS）に沿って中心をシフトさせ、実空間のボイドと赤方偏移空間のボイドとの間の一対一の対応関係を崩してしまう。これらの限界により、動的に動機付けられ、サンプリングの変動に対して堅牢であり、かつ複雑なサーベイ幾何学やRSDを扱うことができるボイド検出アプローチが必要とされている。

手法
著者らは、現在観測されるトレーサー分布から再構成されたラグランジュ変位場における負のダイバージェンス（発散）領域としてボイドを同定する、動的なアルゴルズムであるBack-in-time Void Finder (BitVF) を導入する。

1. 最適輸送（OT）による再構成:

   • 本手法は、離散最適輸送アルゴリズムを用いて、観測されたトレーサー（オイラー的な位置）を、初期の均質なラグランジュ分布へと写像する。
   • 直接的な密度推定とは異なり、このアプローチは、現在の不均質な分布と一様な原始状態との間の最も効率的な写像を見出すために、グローバルな輸送コスト（モンジュ・コスト関数）を最小化する。
   • 実装には、ペアワイズ・スワッピングを行うPIZA法を改良した、効率的な反復的**クァーテット・スワップ・スキーム（quartet-swap scheme）**が用いられ、一致したトレーサーとランダムなペア間の距離の二乗和を最小化する。
   • この再構成は、事前の宇宙論的仮定を必要とせずにトレーサー・バイアスを自然に取り込む。バイアスは、変位場の線形スケーリングを通じて考慮される。

2. ダイバージェンス計算と平滑化:

   • 変位場 $\Psi$ は、時間を遡る形で再構成される。ボイドは、$\nabla_q \cdot \Psi < 0$ となる領域、すなわち、将来に向けて進化させた際に局所的な質量流出に対応する領域として定義される。
   • ダイバージェンス場は、ガウスの定理を用いて、離散グリッド（周期境界条件のボックスにはデカルト座標、ライトコーンには曲線適応座標）上で計算される。
   • ショットノイズを軽減するため、手法は複数のランダムな実現（通常10〜50回）にわたって平均化を行い、オプションとしてダイバージェンス場へのガウス平滑化を適用する。

3. ウォーターシェッド分割:

   • ボイドは、平滑化されたダイバージェンス場に適用されるウォーターシェッド・アルゴリズムを介して同定される。セルは、ダイバージェンス場の局所的な極小値へと反復的に割り当てられ、一貫した流出の連結した盆地（basin）を定義する。
   • ボイドの中心は、近傍のダイバージェンス値に基づく重み付き補正を用いて精緻化される。

4. 赤方偏移空間の補正:

   • 赤方偏移空間のデータに対しては、バイアス補正されたカイザー・フレームワークを用いてRSDをモデル化できる。仮定の宇宙論とトレーサー・バイアスを用いることで、再構成された変位場は実空間の位置の回復を可能にし、RSDによる系統的なシフトを効果的に緩和する。

主な貢献

• 斬新なアルゴリズム: BitVFは、密度やトポロジーの閾値ではなく、ラグランジュ変位を通じて動的にボイドを定義する、ステージIV銀河サーベイ（例：DESI、Euclid）向けに特別に設計された初のフレームワークである。
• ノイズへの堅牢性: 非局所的で一貫した大規模質量流に符号化された変位場に依存することで、BitVFは密度ベースの手法よりもショットノイスやトレーサーの希薄さに対して本質的に堅牢である。
• RSDの緩和: 本手法は、赤方偏移空間において実空間のボイド・サイズ関数を保持する内在的な能力を示しており、仮定の宇宙論を用いて実空間のカタログを再構成する経路を提供する。
• 実装: アルゴリズムは CosmoBolognaLib パッケージに実装されており、並列実行（OpenMP）を備え、数百万のトレーサーに対しても計算量 $O(N^{1.1})$ 程度のスケーラビリティを持つ。

結果
著者らは、高解像度のN体シミュレーション（Aletheia suite）および現実的なライトコーン・モック（Buzzard/DESI様）を用いてBitVFを検証した。

• REVOLVERとの比較: ボロノイ・テッセレーションに基づくトポロジカルなボイド検出器であるREVOLVERと比較して、BitVFは以下の特徴を持つボイドカタログを生成する：
  • より滑らかな密度プロファイル: BitVFのボイドは、トポロジカルな手法でよく見られる「トップハット」構造ではなく、一貫した質量流を反映した、内部と境界の間のより滑らかな遷移を示す。
  • 安定した存在量: REVOLVERが希薄なサンプルにおいて顕著な不安定性と偽の構造生成を示すのに対し、BitVFのボイド数およびサイズ分布は、大幅なトレーサーのサブサンプリング（粒子集団の2.5%まで）やショットノイズの下でも安定している。
  • サイズの差異: BitVFは、REVOLVARよりも系統的に小さな有効半径を持つボイドを同定する。これは、ダイバージェンスが符号を変える境界で膨張を停止し、補償的な過密シェルを除外するためである。
• 速度再構成: OT変位場から導出された再構成速度は、真のハロー速度とのピアソン相関係数 $R \approx 0.77$ を示し、動的な再構成の妥当性を検証している。
• ライトコーンにおける性能: DESI様モックカタログに適用された際、BitVFはRSDの系統的な影響をうまく緩和し、トポロジカルな手法よりも実空間のボイド統計を忠実に保持することに成功した。

意義と主張
本論文は、BitVFが既存のボイド検出器に対する、物理的に動機付けられた堅牢な代替手段を提供すると主張している。特に、トレーサー密度が変化し、RSDの影響が顕著となるステージIVサーベイの時代において、その有用性が高い。動的な質量流出にボイド同定の根拠を置くことで、本手法は以下を実現する：

1. 系統誤差への感受性の低減: トポロジカルなアプローチよりも、ショットノイズやトレーサーの希薄さを本質的にうまく扱う。
2. 解釈性の向上: 動的な定義は宇宙構造の物理的な進化と一致しており、偽のボイドの同定を減少させる可能性がある。
3. 汎用性: 本フレームワークは、理想的なシミュレーション・ボックスと、複雑なマスク付きライトコーン・サーベイデータの両方に適用可能である。

著者らは、BitVFがボイド宇宙論における重要な進歩であり、将来のサーベイで期待される大規模なデータセットに対して、精度、計算効率、および適用性のバランスが取れたツールを提供すると結論付けている。研究を促進するため、コードとツールは公開されている。