[論文レビュー] Fast CMB Analyses via Correlation Functions

István Szapudi, S. Prunet+3 more|arXiv (Cornell University)|Oct 12, 2000

Reservoir Engineering and Simulation Methods参考文献 1被引用数 7

ひとこと要約

本稿では、2点相関関数 $\xi(\theta)$ を最初に計算し、その後ルジャンドル多項式による投影を用いて多極成分に分解することで、CMBの角度スペクトル $\mathcal{C}_\ell$ を高速かつ普遍的に推定する手法を提示する。この手法は $N^2$（可能性として $N\log N$）のスケーリングを示し、ノイズや天球幾何の対称性を仮定せず、メガピクセルマップにおいて不偏でほぼ最適な誤差バーを提供する。これは、将来のミッション（例：Planck や MAP）に適している。

ABSTRACT

We propose and implement a fast, universally applicable method for extracting the angular power spectrum C_l from CMB temperature maps by first estimating the correlation function ξ(θ). Our procedure recovers the C_l's using N^2 (but potentially N logN), operations, where N is the number of pixels. This is in contrast with standard maximum likelihood techniques which require N^3 operations. Our method makes no special assumptions about the map, unlike present fast techniques which rely on symmetries of the underlying noise matrix, sky coverage, scanning strategy, and geometry. This enables for the first time the analysis of megapixel maps without symmetries. The key element of our technique is the accurate multipole decomposition of ξ(θ). The C_l error bars and cross-correlations are found by a Monte-Carlo approach. We applied our technique to a large number of simulated maps with Boomerang sky coverage in 81000 pixels. We used a diagonal noise matrix, with approximately the same amplitude as Boomerang. These studies demonstrate that our technique provides an unbiased estimator of the C_l's. Even though our method is approximate, the error bars obtained are nearly optimal, and converged only after few tens of Monte-Carlo realizations. Our method is directly applicable for the non-diagonal noise matrix. This, and other generalizations, such as minimum variance weighting schemes, polarization, and higher order statistics are also discussed.

研究の動機と目的

• 大規模で非対称なマップから、ノイズや天球幾何の対称性を仮定せずに、CMBの角度スペクトル $\mathcal{C}_\ell$ を高速かつ汎用的に推定する手法の開発。
• 将来のミッション（例：MAP や Planck）向けに、特にメガピクセルCMB調査において、標準的な最尤法の $N^3$ スケーリングが非現実的になる問題を克服すること。
• システムティクス、ノイズ、foregroundのモンテカルロスタディを劇的に短縮された解析時間で可能にする。
• 非対角ノイズ行列やその他の複雑な観測条件（不規則な天球カバー、偏光など）に対しても一般化可能であるように、手法を拡張すること。
• 相関関数が、CMB宇宙論における従来のスペクトル推定法の実用的で効率的な代替手段であることを示すこと。

提案手法

• まず、CMB温度マップから高分解能のビンで2点相関関数 $\xi(\theta)$ を推定する。
• ルジャンドル多項式の根を中心としたガウス平滑化を適用し、$\xi(\theta)$ の多極成分への効率的分解を可能にする。
• $\mathcal{C}_\ell$ の係数は、ピクセル対の重み付き和により回復され、重みはルジャンドル多項式の投影から導出される。これは2次推定器を形成する。
• 誤差バーと相関行列は、ノイズ行列の直接逆行列を計算せずに、人工ノイズ実現を用いたモンテカルロ手法により計算される。
• 非対角ノイズに対しては、ノイズ行列の逆行列を近似するためにモンテカルロ実現を用いる。これにより、フル行列を保存できない場合でも応用可能となる。
• 反復的ピクセルペア重み付けと $k$-点相関関数を用いることで、最小分散重み付け、偏光、高次統計への一般化が可能となる。

実験結果

リサーチクエスチョン

• RQ1 対称性の仮定なしに、$N^2$ またはそれ以下のスケーリングで、ほぼ最適な $\mathcal{C}_\ell$ 誤差バーを達成できるか。
• RQ2 バルーン搭載実験のような非一様ノイズを持つ大規模で非対称なCMBマップにおいて、この手法はどの程度の性能を示すか。
• RQ3 大きなノイズ行列の逆行列を計算せずに、誤差バーを効率的に推定するためにモンテカルロ手法をどの程度活用できるか。
• RQ4 現実の観測状況において、非対角ノイズ行列や最小分散重み付けスキームに一般化可能か。
• RQ5 2点相関関数は、将来的な高分解能調査において、直接的な $\mathcal{C}_\ell$ 推定の実用的で効率的な代替手段として適しているか。

主な発見

• 81,000ピクセルのシミュレーテッドBoomerang類似マップと対角ノイズを用いた実験で、標準的な2次推定器と同等の不偏な $\mathcal{C}_\ell$ 係数推定が達成された。
• 誤差バーは、宇宙分散とノイズによって制約される理論的最小値から最大で10%以内に収まり、ほぼ最適性を示している。
• 解析時間は、標準的な $N^3$ 手法で数か月かかっていたのを、数時間にまで短縮した。これは大規模CMB解析において顕著な性能向上を示している。
• 基本形では $N^2$ スケーリングを示し、高度なアルゴリズムを用いることで $N\log N$ に最適化可能であり、将来のミッションの $N \sim 10^6 - 10^7$ マップに対しても実用的である。
• ノイズの非対称性、銀河平面カット、不規則な天球カバーに対して、この手法は、ノイズの対称性やマップ幾何の仮定を一切行わないため、頑健である。
• 短い解析時間のおかげで、システムティクスやforegroundの実用的モンテカルロスタディが可能となり、$N^3$ 手法では以前は非現実的だった。