Dynamical Test of Cosmic Acceleration: $k$-nearest Neighbor Cross Correlation of Cosmic Microwave Background and Cosmic Infrared Background

技術要約：kNN相互相関による宇宙加速の動的テスト

問題と動機
宇宙加速の物理的な起源は、現代宇宙論における中心的な未解決問題である。幾何学的なプローブ（例：Ia型超新星、BAO）は距離・赤方偏移関係を制約するが、動的なプローブは、大規模構造の成長と重力ポテンシャルの進化を直接追跡することにより、独立した情報を提供する。積分サックス・ヴォルフェ効果（ISW）およびリース・サイマ効果（RS）—これらを総称してISWRSと呼ぶ—は、ダークエネルギーが支配的な宇宙において、宇宙マイクロ波背景放射（CMB）光子が時間変化する重力ポテンシャルを通過する際に発生する。ISWRS信号の検出は、誘起される温度異方性が一次的なCMBのゆらぎよりも一桁小さいため、困難を伴う。大規模構造トレーサー（例：SDSS、2MASS、WISE）を用いた相関による従来の検出では、通常、有意水準 $\gtrsim 4\sigma$ に達するために複数のサーベイをスタッキングする必要があった一方、空洞（voids）や銀河団のスタッキング解析では、より低い有意水準（$\lesssim 3\sigma$）しか得られていない。

手法
本研究では、$k$ 近傍累積分布関数（kNN-CDF）という新しい統計的手法を用いて、CMBと宇宙赤外背景（CIB）の間の相互相関を測定する。

• データセット:

  • CIBトレーサー: 著者らは、Y.-K. Chiang (2023, C23) によって再構成された、前景放射除去済みの 100 $\mu$m CIB マップを利用している。このマップは、$z \approx 2.5$ まで広がる約6億個の Wide-field Infrared Survey Explorer (WISE) 銀河から導出されており、実質的に銀河–CMB 相互相関トレーサーとして機能する。このマップは、分光参照サンプルとの相関を用いることで、銀河系内のダスト成分を銀河外成分から分離している。
  • CMBデータ: Planck Public Data Release 3 (PR3) の温度マップ、特に Sunyaev–Zel'dovich (SZ) 効果による汚染を最小限に抑えるための「SZフリー」（SMICA-nosz）製品を使用する。
  • 模擬データ: $\Lambda$CDM モデルに対する有意性を定量化するために、著者らは Millennium-XXL (MXXL) N体シミュレーションを用いて、模擬的な CIB および ISWRS マップを構築する。模擬 CIB は、シミュレーションのライトコーン・ハローに対してハロ赤外放射モデル (W18) を適用し、観測された CIB パワースペクトルおよび確率密度関数 (PDF) に一致するようにスケーリングとノイズ項を加えることで作成される。ISWRS マップは、視線方向に沿った重力ポテンシャルの時間微分の積分によって生成される。

• 統計的枠組み:

  • kNN-CDF 統計量 $\psi$ は、2つの離散的または連続的な場において、与えられた体積内に近傍点が存在する結合確率を測定する。ここでは、トップハット平滑化された CMB および CIB マップに適用される。
  • 解析には、固定されたパーセンタイル閾値（$T^*, I^* \in \{50\%, 70\%, 90\%\}$）と角平滑化スケール（$\theta \in \{0.1^\circ, \dots, 0.9^\circ\}$）を用い、結果として45要素のデータベクトルが得られる。
  • 共分散行列は、200個のランダムな CIB および CMB マップを用いて推定される。逆共分散行列を補正するために Hartlap 係数が適用される。
  • 系統誤差については、CIB 再構成誤差に基づくマスクの構築と、外れ値ピクセルの除外（基準カット：CMBは0%、CIBは2%）によって対処される。

主な貢献と結果

1. ISWRS の検出: 本研究は、100 $\mu$m CIB と CMB の間の正の相互相関を報告している。

   • 相関がないという帰無仮説を棄却すると、有意性は $p = 0.02$ ($\chi^2_{\text{null}} \approx 66.6$、自由度 45、これは $2.3\sigma$ に相当) となる。
   • ISWRS 効果を含む $\Lambda$CDM 模擬データに対する振幅を定量化すると、最良適合振幅は $A_{\text{kNN}} = 0.95 \pm 0.20$ である。これは $4.8\sigma$ の検出有意度に対応する。
   • 個々の角ビンにおける帰無仮説からの偏差は、$1.8\sigma$ から $4.1\sigma$ の範囲にある。

2. kNN-CDF の優位性:

   • 著者らは、同じデータとマスクを用いた標準的な 2 点角パワースペクトル解析 ($C_\ell$) と、自らの kNN-CDF 結果を比較している。
   • 2 点解析による検出有意度は、線形理論で $\approx 3.9\sigma$、シミュレーションベースで $\approx 4.0\sigma$ である。
   • kNN-CDF 法は、検出有意度を約 20% （$4.8\sigma$ 対 $4.0\sigma$）向上させる。著者らは、この利得は単に広いフーリエモードへのアクセスによるものではなく、kNN 統計量が非ガウス的な高次クラスタリング情報を捉える能力によるものであるとしている。

3. 系統誤差と頑健性:

   • 様々な CMB 成分分離パイプライン（SMICA, NILC, SEVEM, Commander, L-GMCA）を用いたテストにより、残留熱的 SZ 汚染を持つマップ（例：標準的な SMICA）は振幅が低くなることが示され、高密度領域における SZ 残留が信号を抑制することを裏付けている。「SZ フリー」マップが最も頑健な結果を提供している。
   • 銀河系ダストマップとの相互相関において、有意な汚染は見られなかった。

4. トモグラフィー解析:

   • ダークエネルギーの進化を調べるため、CIB マップを $z_{\text{cut}} \approx 0.7$ で 2 つのトモグラフィー・ビンに分割した。
   • 低赤方偏移ビンおよび高赤方偏移ビンの振幅は、それぞれ $A_{\text{low-}z} = 1.41 \pm 0.47$ ($3.0\sigma$) および $A_{\text{high-}z} = 1.66 \pm 0.81$ ($2.0\sigma$) である。
   • 不確実性が大きいため、これらをフルマップの結果や $\Lambda$CDM 予測と明確に区別するには至っていないが、結果は $\Lambda$CDM と矛盾しない。

意義と主張
本論文は、kNN-CDF 法が、従来の 2 点相関関数よりも強力な宇宙加速の動的テストを提供し、前景放射除去済みの単一の CIB トレーサーを用いて ISWRS 効果の $4.8\sigma$ 検出を達成したと主張している。この検出は $\Lambda$CDM モデルと一致している。著者らは、この有意性の向上は、手法が非ガウス的なクラスタリング情報に対して敏感であることに由来すると強調している。

本研究は、CIB のより精密なトモグラフィー再構成、特に SPHEREx のような将来のサーベイによるマルチバンドデータを用いることが、進化するダークエネルギーの制約をさらに強固にすることを示唆して結論づけている。CIB 再構成と kNN-CDF 統計量の組み合わせは、現在のおよび今後のサーベイの高い全天被覆率と膨大な銀河数を活用した、将来のダークエネルギー探査のための強力なフレームワークとして提示されている。