[論文レビュー] MasQU: Finite Differences on Masked Irregular Stokes Q,U Grids

Jude Bowyer, Andrew Jaffe+1 more|arXiv (Cornell University)|Jan 3, 2011

Geophysics and Gravity Measurements参考文献 4被引用数 7

ひとこと要約

本稿では、特にHEALPixに適した球面グリッド上の数値微分の一般アルゴリズムを用いて、マスク処理された不規則なストークスQおよびUグリッドにおける有限差分を計算するソフトウェアパッケージであるMasQUを提示する。ハーモニック法と比較して、特に大スケールにおいて、EモードからBモードへの漏れを少なくとも1桁以上低減する。また、極特異性の問題に対処し、ガウス求積法を用いて高精度なハーモニック係数推定を可能にする。

ABSTRACT

The detection of B-mode polarization in the CMB is one of the most important outstanding tests of inflationary cosmology. One of the necessary steps for extracting polarization information in the CMB is reducing contamination from so-called "ambiguous modes" on a masked sky, which contain leakage from the larger E-mode signal. This can be achieved by utilising derivative operators on the real-space Stokes Q and U parameters. This paper presents an algorithm and a software package to perform this procedure on the nearly full sky, i.e., with projects such as the Planck Surveyor and future satellites in mind; in particular, the package can perform finite differences on masked, irregular grids and is applied to a semi-regular spherical pixellization, the HEALPix grid. The formalism reduces to the known finite-difference solutions in the case of a regular grid. We quantify full-sky improvements on the possible bounds on the CMB B-mode signal. We find that in the specific case of E and B-mode separation, there exists a "pole problem" in our formalism which produces signal contamination at very low multipoles l. Several solutions to the "pole problem" are presented; one proposed solution facilitates a calculation of a general Gaussian quadrature scheme, which finds application in calculating accurate harmonic coefficients on the HEALPix sphere. Nevertheless, on a masked sphere the software represents a considerable reduction in B-mode noise from limited sky coverage.

研究の動機と目的

• マスク処理されたスカイにおけるCMB偏光解析において、非局所的信号混合のためハーモニック法が失敗するというEモードからBモードへの漏れの問題に対処すること。
• 特にHEALPixピクセル化に適した、不規則でマスク処理された球面グリッド上の有限差分を計算する汎用アルゴリズムの開発。
• 球面有限差分法における極特異性問題を克服し、これにより低多重数Bモード信号が汚染されないようにすること。
• 導出されたガウス求積スキームを用いて、HEALPix球面上でのハーモニック係数の高精度な計算を可能にすること。
• マスク領域でのサンプリングを回避する実空間代替手法を提供し、原始的Bモード信号検出のための信号忠実度を向上させること。

提案手法

• 任意のグリッドタイプおよび微分次数に対して、固定ステンシルに依存しないアルゴリズム的アプローチにより有限差分重みを定式化する。
• カーテシアンおよびHEALPixグリッドにこの手法を適用し、θおよびφ座標の特異性が標準的差分法を破壊する極部に対して特別な処理を行う。
• スピン重み付き球面調和関数の摂動展開とWigner d関数の再帰関係を用いて数値的重みを導出する。
• 小d関数とジャコビ多項式の直交性を活用し、調和空間における微分作用素の畳み込み積分を評価する。
• 導出された形式から一般化されたガウス求積スキームを構築し、HEALPix球面上での高精度なハーモニック係数推定を実現する。
• EBEX型実験設定を含むノイズとスカイマスクを伴うモックデータを用いて、標準的なハーモニック技術との妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 不規則でマスク処理された球面グリッド（特にHEALPix）上で、極部を含めて有限差分をどのように高精度に計算できるか。
• RQ2 提案手法が、マスク処理されたスカイにおける従来のハーモニック分解と比較して、EモードからBモードへの漏れをどの程度低減するか。
• RQ3 HEALPix極帽子におけるピクセル配置が、有限差分導関数の精度に与える影響は何か。
• RQ4 有限差分形式から一般化されたガウス求積スキームを導出可能か。これにより球面上でのハーモニック係数推定がどの程度向上するか。
• RQ5 ノイズおよび部分的スカイカバレッジの下で、特に低多重数Bモード信号検出において、この手法の性能はいかにか。

主な発見

• MasQUアルゴリズムは、EモードからBモードへの漏れを少なくとも1桁以上低減し、特に大きなスケール（低多重数）でその改善効果が増大する。
• 極問題を効果的に緩和するため、一般化されたガウス求積スキームを導出し、高精度なハーモニック係数計算を可能にした。
• HEALPixグリッド上では、特にスカイマスクが存在する状況において、標準的なハーモニック技術よりも顕著に低いBモード汚染を達成した。
• この形式は、正則グリッド上では既知の有限差分解に還元され、一般的手法の整合性と正しさを確認した。
• マスク領域でのサンプリングを回避するため、完全スカイ解析が可能となり、ハーモニック分解に内在するグローバル信号混合効果を回避した。
• 球面上での微分作用素の畳み込み演算子は、多重数l, l′およびモードm, m′の間でパワーを混合するが、Wigner d関数の再帰関係を用いて解析的に定量化された。