Correlation Calibration: A Hybrid Calibration Technique for Radio Interferometric Arrays

技術要約：「Correlation Calibration: A Hybrid Calibration Technique for Radio Interferometric Arrays」

問題提起
ラジオ干渉計データの較正（キャリブレーション）は、特に要素が固定されており校正源を追跡できないドリフトスキャン・アレイ（例：HERA、CHIME、CHORD、HIRAX）において、根本的な課題である。従来の較正手法は、正確な電波天のモデルと観測装置の応答を必要とするスカイベースの手法、あるいはアンテナ配列の規則性に依存する冗長較正のいずれかに依存している。両方のアプローチには特有の限界がある。スカイベースの手法は、空または観測装置のモデリング誤差の影響を受けやすく、一方で冗長較称は、物理的な非冗長性（アンテナ配置の誤差やディッシュの個体差）によって制限され、絶対フラックススケールやフェーズセンター（位相中心）を制約できない。21cm宇宙論の文脈において、系統誤差が宇宙論的信号を圧倒する可能性があるため、モデリング誤差に対して堅牢で、かつ大規模なアレイに対して計算効率の高い較正戦略が切実に求められている。

手法
本論文では、計算効率の高い共分散ベースのフレームワーク内で、スカイベース較正と冗長較正の両方の強みを活用したハイブリッド較正技術であるCorrCal（Correlation Calibration）を紹介する。

• 定式化: データとモデル期待値の残差を最小化する（$\chi^2$を最小化する）従来の較正とは異なり、CorrCalは、観測された共分散に対してデータの共分散モデルを適合させることを目的としたパラメータ推定問題として較正を扱う。データは、熱雑音、点源放射、および拡散放射の3つの独立したガウス成分の混合としてモデル化される。
• 共分散モデリング: モデル共分散 $C$ は、以下の3つの項の和として構成される：
  1. 雑音 (Noise, $N$): 対角成分の熱雑音分散。
  2. 点源 (Point Source, $S$): 少数の明るい点源からの寄与を含む疎なソース行列 $\Sigma$ を用いてモデル化される。
  3. 拡散放射 (Diffuse Emission, $D$): 冗長グループ内のベースライン間共分散の固有値分解から導出された疎な拡散行列 $\Delta$ を用いてモデル化される。この項は、アレイの冗長性と拡散天の期待パワースペクトルをエンコードしている。
• 疎な表現 (Sparse Representation): 本研究の主要な革新は、完全な $2N \times 2N$ のベースライン間共分散行列（ここで $N$ はベースライン数）が、雑音行列、ゲイン行列、ソース行列、および拡散行列を用いて疎に表現できるという認識である。これにより、密な共分散行列を逆行列化することなく、Woodburyの恒等式と行列行列式補題（Matrix Determinant Lemma）を用いて、尤度とその勾配を効率的に計算することが可能となる。
• 最適化: 較正解は、測定されたデータがモデル共分散を与えるという負の対数尤度を最小化することによって得られる。アルゴリズムは、共役勾配法ソルバーを使用して、アンテナごとの複素ゲインを算出する。

主な貢献

1. ハイブリッド・フレームワーク: CorrCalは、冗長較正（拡散共分散項を介して）とスカイベース較正（点源項を介して）を単一の尤度関数へと統合する。これにより、外部の空の情報を利用しながら、内部的なデータの整合性を強制することが可能になる。
2. 計算効率: 共分散構造の疎性を活用することで、著者らはアルゴリズムが効率的にスケールすることを示している。その実装では、1,000要素のアレイに対して、単一のラップトップのコアを用いて1秒未満で尤度と勾配を評価できる。
3. モデリング誤差への堅牢性: 著者らは、広範なエンドツーエンドの検証テストを通じて、CorrCalが空のモデル（例：不正確なソースのフラックス密度、欠落したソース）やアレイの非冗長性（例：アンテナ配置の誤差）に対して鈍感であることを示している。
4. 正規化の改善: 著者らは、アンテナあたりの冗長グループの数を考慮していなかったために生じていたゲイン解の曲率のアーティファクトを修正するために、尤度を適切に正規化することで、オリジナルのCorrCalアルゴリズム（Sievers 2017）を改良した。

結果
本論文は、一連のシミュレーションを通じてCorrCalを検証している：

• 空のモデリング誤差: 較正ソースの数とソースのフラックス密度の精度を変化させたテストにより、較正誤差はバイアスがなく、熱雑音の閾値を十分に下回っていることが示された。解の質は、含まれるソースの数にほとんど依存しない（少数の明るいソースが含まれていればよい）。
• 非冗長性: アンテナ位置にガウス乱数摂動を加えた（現実世界の配置誤差を模した）シミュレーションでは、CorrCalは標準的な冗長較正を大幅に上回る性能を示した。モデルに正確なアンテナ位置を使用した場合、較正誤差は熱雑音に匹敵するサブパーセントレベルまで減少した。
• 非ガウス的な空: 現実的な空のモデル（GLEAMカタログ + Global Sky Model）を用いたテストにより、拡散前景放射における非ガウス性が較正の質に大きな影響を与えないことが確認され、拡散成分に対するガウスランダム場仮定の妥当性が示された。
• 位相勾配: 研究では、不完全な空のモデル（少ないソース）を使用すると、ゲイン解に空間的な位相勾配が導入され、推論される空の位置が実質的にシフトすることが特定された。この勾配は、明るいソースが天頂付近にあるときは制約されるが、それらがない場合には（冗長較正と同様に）制約されず、退化（デジェネレート）する挙動を示す。

意義と主張
本論文は、CorrCalが現在の、そして次世代のドリフトスキャン型ラジオ干渉計を較正するための、堅牢かつ計算効率の高いツールを提供すると主張している。共分散を可変のオブジェクトとして扱い、疎性を活用することで、CorrCalは完璧な空のモデルや完全に冗長なアレイを必要とせずに、高精度な較正を実現する。著者らは、この手法が系統誤差の制御が極めて重要となる21cm宇宙論実験において特に有用になると考えている。著者らは、CorrCalが較正誤差を大幅に低減するものの、不完全な空のモデリングによって誘発される小さな振幅オフセットや位相スロープを補正するために、実用的なアプリケーションでは補完的な絶対較正ステップが必要になる可能性があると述べている。著者らは、CorrCalが絶対較正の必要性を排除すると主張しているのではなく、むしろ、後続のステップの負担を最小限に抑える、非常に精度の高い相対較正を提供するものであると主張している。