[論文レビュー] CMB polarization TE power spectrum estimation with non-circular beam

Fidy A. Ramamonjisoa, Subharthi Ray+2 more|arXiv (Cornell University)|Sep 18, 2013

Radio Astronomy Observations and Technology参考文献 3被引用数 3

ひとこと要約

本論文は、擬似-Cℓ推定器を用いた高速な半解析的フレームワークを提示し、非円形ビームを伴うCMB TEパワー スペクトルを計算する。ビーム非対称性は摂動的調和展開によりモデル化される。𝓁_max = 500までで、1秒未満のバイアス行列計算が達成され、Planck LFI 30 GHzでは最大20%、WMAPバンドでは12–16%のビーム由来の不確実性が明らかになった。これは高解像度偏光実験における重要な系統誤差を示している。

ABSTRACT

Precise measurements of the Cosmic Microwave Background (CMB) anisotropy have been one of the foremost concerns in modern cosmology as it provides valuable information on the cosmology of the universe. However, an accurate estimation of the CMB power spectrum faces many challenges as the CMB experiments sensitivity increases. Furthermore, for polarization experiments the precision of measurement is complicated by the fact that the polarisation signal is very faint compared to the measured total intensity, and could be impossible to detect in the presence of high level of systematics. One of the most important source of errors in CMB polarization experiment is the beam asymmetry. For large data set the estimation of the CMB polarization power spectrum with standard optimal Maximum Likelihood (ML) is prohibitive for high resolution CMB experiments due to the enormous required computation time. In this paper, we present a semi-analytical framework using the pseudo-$C_{l}$ estimator to compute the power spectrum TE of the temperature anisotropy and the E-component of the polarization radiation field using non-circular beams. We adopt a model of beams obtained from a perturbative expansion of the beam around a circular (axisymmetric) one in harmonic space, and compute the resulting bias matrix which relates the true power spectrum with the observed one using an efficient algorithm for rapid computation. We show that for a multipole up to $l_{max}=500$, the bias matrix can be computed in less than one second with a single CPU processor at 2.53 GHz. We find that the uncertainties induced by the beam asymmetry in the polarization power spectrum at the peak of the bias matrix for WMAP and Planck experiments can be as large as a few 10 to 20% (upper limit for Planck LFI 30 GHz)

研究の動機と目的

• 高解像度実験における非円形（非対称）ビームによって引き起こされるCMB偏光パワー スペクトル推定の系統誤差を解消すること。
• 大規模データセットでは実行が困難な最適な最尤推定法の計算効率の良い代替手法を開発すること。
• ビーム楕円度がTEパワー スペクトルの不確実性に与える影響を、特に高多重度において定量化すること。
• 近い将来のCMB実験（例：Planckや将来の調査）におけるビーム由来バイアスの正確な補正を可能にすること。

提案手法

• 楕円度パラメータを用いて、軸対称ビームを基準として、調和空間における非円形ビームの摂動的展開によるモデル化。
• ビーム非対称性によって生じる真のパワー スペクトルと観測スペクトルとの関係を記述するバイアス行列を、高速アルゴリズムを用いて導出。
• 計算コストを削減しつつ精度を保持するため、バンド制限付きアプローチ（|𝓁 − 𝓁′| ≤ 20）を用いて多重度間の結合を計算。
• 非円形幾何形状を正確に捉えるために、Tモードではm = 2,4、Eモードではm = 4,6の高次ビーム補正を組み込む。
• 事前計算されたWigner-d関数と効率的なClebsch-Gordan係数評価を活用し、行列計算を高速化。
• 結果の妥当性をWMAPおよびPlanckのビームパラメータと比較し、𝓁_max = 3000への外挿も実施。

実験結果

リサーチクエスチョン

• RQ1 非円形ビームが高解像度実験におけるCMB TEパワー スペクトル推定に及ぼすバイアスはいかなるものか？
• RQ2 高多重度まで非円形ビームのバイアス行列を計算する上で、計算の実行可能性はどの程度か？
• RQ3 ビーム楕円度とビーム幅は、TEスペクトルにおけるバイアスピークの振幅と多重度位置にどのように影響するか？
• RQ4 時間領域データにおける有効ビーム楕円度は、スカイメップにおける真の系統誤差をどの程度過大評価するか？
• RQ5 計算コストの高い最尤推定器に代わる、高速で半解析的なフレームワークが、ビーム系統誤差補正に利用可能か？

主な発見

• 2.53 GHz CPU 1個を用いて、𝓁_max = 500 までで、非円形ビームのバイアス行列を1秒未満で計算可能であり、迅速な系統誤差評価が可能になった。
• Planck LFI 30 GHzでは、TEパワー スペクトルにおけるビーム由来の不確実性が最大20%に達し、有効楕円度ε ≈ 1.36であった。
• WMAPでは、Qバンド（ε ≈ 1.15）で12%、Vバンド（ε ≈ 1.09）で16%の不確実性が示され、高解像度において顕著な系統誤差が確認された。
• バイアスピークはビーム幅の逆数に相当する多重度（𝓁_peak ≈ 1/σ）に現れ、ビーム楕円度が大きいほど振幅が増加する。
• バイアス行列は強く非対角ではないが、多重度間の結合が顕著（𝓁σ ≳ 1）であり、対角成分から離れるに従い急速に減衰する。
• スカイメップにおける有効ビーム楕円度は、時間領域ビームモデルに比べて抑制されていることから、シミュレーションによるビーム効果は実際の系統誤差を過大評価している可能性がある。