Analytical method for computing the covariance matrix of cosmic shear two-point correlation function

技術的サマリー：宇宙せん断の2点相関関数の共分散行列を計算するための解析的手法

問題提起
現在のおよび将来の広域弱重力レンズ調査を用いた宇宙論的解析において、宇宙せん断統計の共分散行列の正確な推定は不可欠である。サンプル共分散は多数のシミュレーション実現から測定可能であるが、このアプローチは、全天シミュレーションから少数の独立した実現しか抽出できない10,000平方度を超える天域を持つ第4段階（Stage-IV）の調査においては、計算量的に実行不可能となる。したがって、共分散の解析的モデリングが必要となる。しかし、調査マスクはモード結合効果をもたらすため、特に実空間の2点相関関数（2PCF）において、ガウス共分散項の計算を複雑にする。既存の解析的手法は、複雑な調査幾何学が存在する状況下で、共分散構造を正確に再現することにしばしば苦慮している。

手法
著者らは、有限な調査幾何学（調査ウィンドウ効果）を考慮した、宇宙せん断2PCF $\xi(\theta)$ のガウス共分散行列を計算するための解析的手法を調査する。本研究では、HSC 第3年（HSC-Y3）調査マスクを用いて生成された1,000個のガウス実現の宇宙せん断場をベンチマークとして利用し、解析的予測をテストする。

解析的枠組みは以下の通り進行する：

1. 調和空間から実空間へのマッピング：本研究は、ウィンドウ付きパワースペクトル $\tilde{C}_\ell$ と測定された相関関数 $\xi_i$ との関係を確立する。観測された相関関数は、ウィンドウの影響を受けたパワースペクトル $\tilde{C}_\ell$ から導出され、調査領域内の銀河ペアの数（$N_{pair}$）と平均数密度（$n$）に依存する因子によって補正された不偏推定量として扱われる。
2. 共分散の射影：実空間の2PCF の共分散 $\text{Cov}(\xi_i, \xi_j)$ は、調和空間の共分散 $\text{Cov}(\tilde{C}_\ell, \tilde{C}_{\ell'})$ の二重ルジャンドル変換として形式的に表現される。
3. 解析的推定量：調和空間の共分散を推定するための3つの具体的なアプローチが評価される：
   • 改良狭核近似（iNKA）：モード結合行列が $\ell \simeq \ell'$ 付近で鋭く局在化していると仮定し、$\text{Cov}(\tilde{C}_\ell, \tilde{C}_{\ell'})$ の高速計算を可能にする。
   • $f_{sky}$ 近似：共分散が対角成分のみを持ち、天域割合によってスケーリングされると仮定する標準的な近似。
   • 重み付き4乗カウント法：Schneider ら（2002）の形式に従い、$\xi \times \xi$ で重み付けされた銀河ペアの4乗カウントを用いて、構成空間で直接2PCF 共分散を計算する。

主要な結果
本研究は、1,000個のガウスシミュレーションから直接測定されたサンプル共分散と解析的予測を比較する：

• 枠組みの検証：ウィンドウ畳み込みパワースペクトルと相関関数の間の関係（式2.6）が検証され、疑似-$C_\ell$ 枠組みが調査ウィンドウ効果を正しく記述し、調和空間と実空間の間の正確なマッピングを提供することが確認された。
• iNKA の限界：iNKA は調和空間共分散 $\text{Cov}(\tilde{C}_\ell, \tilde{C}_{\ell'})$ の対角成分を10%以内で正確に再現するが、非対角成分を正確にモデル化することはできない。これらの調和空間共分散がルジャンドル変換を介して実空間へ射影されると、非対角項の残差誤差が顕著に伝播する。その結果、iNKA に基づく $\xi(\theta)$ の共分散は、2PCF 共分散行列の対角成分であっても、シミュレーションから測定された共分散を再現しない。
• 4乗カウント法の優位性：対照的に、重み付き4乗カウント法はシミュレーション結果と著しく良い一致を示す。この手法は構成空間におけるペアベースの相関を明示的にモデル化し、iNKA の対角優位性の仮定に依存することなく、共分散構造を成功裡に再現する。
• $f_{sky}$ 近似：$f_{sky}$ 近似もシミュレーションと比較して不一致を示すが、テストされた解析的推定量の中では4乗カウント法が最も正確であることが強調される。

意義と主張
本論文は、調査ウィンドウ効果が存在する場合、実空間の弱重力レンズ統計の信頼性の高い共分散推定を得るためには、調和空間共分散の非対角構造の正確なモデリングが不可欠であると結論づける。著者らは、ルジャンドル変換中に非対角誤差が増幅されるため、単に iNKA を実空間統計に拡張するだけでは不十分であることを実証する。

この研究の主な貢献は、2PCF 共分散の文脈における iNKA の特定の失敗様式を特定し、第4段階（Stage-IV）調査のためのより頑健な解析的代替手段として重み付き4乗カウント法を検証した点にある。これらの結果は、精密な宇宙論的パラメータ推定を必要とする将来の広域調査が、調和空間における単純な対角近似を超えてモード結合を捉える解析的手法を採用しなければならないことを示唆している。