Synthesis imaging with a lunar orbit array: I. global sky map and its systematics

技術要約：月周回アレイによる合成イメージング：I. 全天マップとその系統誤差

問題提起
30 MHz以下の周波数における地上ベースの電波天文学は、電離層による吸収と地球上の無線周波数干渉（RFI）によって深刻な制限を受けている。DSL（Discovering Sky at the Longest wavelength）ミッション、別名「紅盟（Hongmeng）」ミッションは、これらの制限を克服するために、月軌道を周回する線形衛星アレイを提案している。このようなアレイによって収集される可視度データは、スカイマップに対して線形に関連しているが、全天マップの再構成には、地上ベースの干渉計とは異なる特有のアルゴリズム上の課題が存在する。これらの課題には、月が空の一部を遮ることによる不完全な$uvw$カバレッジ、不良設定の逆問題（ill-posed inversion problems）を解くための正則化の必要性、および空の有限なピクセル化に起因する系統誤差が含まれる。具体的には、本論文は、ビーム行列の構築中に適切に処理されない場合に、ピクセルよりも小さいスケールのノイズ（サブピクセルノイズ）がいかにして重大なエイリアシング効果を誘発するかを調査している。

手法
著者らは、DSLアレイ構成における全天マップ再構成アルゴリズムを評価するために、数値シミュレーションを用いている。

• スカイモデル: 既存の全天モデル（GSM、SSMなど）における小スケール成分の不足に対処するため、著者らは、高解能な銀河拡散放射モデル（拡張されたULSAモデルに基づく）と、GLEAMサーベイのカウントから導出された点源の集団を組み合わせた、模擬フルスカイマップを構築した。このマップは$N_{side} = 2048$で生成され、シミュレーション入力用に$N_{side} = 1024$へとダウングレードされている。
• アレイ構成: シミュレーションでは、円形月軌道（高度300 km、傾斜角30°）上にある8基の娘衛星（daughter satellites）をモデル化しており、これらは1.3年の歳差運動周期を持つ。衛星は「ブリージング（呼吸）」戦略を採用しており、基線長を変化させるために、そのフォーメーションを周期的に圧縮および伸長させている。基線分布は100 mから100 kmの範囲に最適化されている。
• 信号モデリング: 可視度データは、熱雑音、有限の時間分解能、および有限の帯域幅効果を組み込んだ、空の強度に対する線形写像としてモデル化されている。月面反射は、基線長に対する遅延のため、コヒーレンスへの影響は無視できると推定されているが、ソースのフラックスをわずかに増加させる。
• 再構成アルゴリズム: 著者らは、チホノフ正則化を用いて線形逆問題を解いている。ビーム行列（$B$）を構築するための2つの手法を比較している：
  1. ピクセル中心法（Pixel-center method）: 各ピクセルの中心におけるビーム応答を評価する。
  2. ピクセル平均法（Pixel-averaging method）: より高解像度で応答を評価し、ターゲット解像度までダウンサンプリングして平均化することで、サブピクセル変動を平滑化する。
• 指標: 再構成の品質は、入力マップと再構成マップの相関係数（$\rho_l$）、球面調和関数空間における信号対雑音比（SNR）、および熱雑音と不完全なビーム効果の誤差共分散行列のトレースを用いて定量化される。

主な貢献と結果

1. ピクセル平均化によるエイリアシングの軽減: 本研究は、ナイーブにピクセル中心法を用いてビーム行列を構築すると、ナイキスト・サンプリング限界（$b > \lambda/\theta_p$）を超える基線が含まれた場合に、深刻なエイリアシング・アーティファクトが生じることを示している。これらのアーティファクトは、重大な歪みと相関係数の低下として現れる。著者らは、ピクセル平均法がこのエイリアシングを効果的に軽減し、マップの品質を低下させることなく、より長い基線の導入を可能にすることを示している。
2. 基線の選択: 解析の結果、$b < 2b_p$（ここで$b_p$はピクセルサイズに対応する基線）までの基線を含めることは、再構成に十分であることが明らかになった。より長い基線（$b > 2b_p$）を含めても、平均化プロセスによってその寄与がほぼゼロに減衰するため、最終的なマップへの改善は無視できる程度である。
3. 正則化パラメータのトレードオフ: 本論文は、正則化パラメータ（$\epsilon$）の選択について調査している。熱雑音を抑制すること（大きな$\epsilon$が有利）と、ビームの不完全性による誤差を最小限に抑えること（小さな$\epsilon$が有利）の間にはトレードオフが存在する。最適な$\epsilon$は、誤差曲線の交点で見出されるが、著者らはこれはケース依存であると注記している。
4. 極域の系統誤差: 再構成において、軌道傾斜角に起因する短い投影基線の欠如により、極域で系統的な負のバイアス（暗化）が生じる。この効果は、投影基線が波長単位で実質的に長くなる高周波数（例：10 MHz）においてより顕著になる。事前マップ（prior map）を用いることで、このパワー損失を緩和できるが、これは事前情報の正確さへの感度を導入することになる。
5. 周波数依存性: システム温度が高いために低周波数ほどノイズの多い画像が得られるという予想に反して、シミュレーションでは、3 MHzにおける明るいソースの信号対雑音比は10 MHzよりも実際には良好であることが示された。これは、強いソースの輝度温度が高いこと、および（波長が長いため）より多くの基線が再構成に含まれることに起因している。

意義と主張
本論文は、DSLミッションの合成イメージングのための堅牢なアルゴリズム的枠組みを提供することを主張している。その主な意義は、標準的な地上ベースの手法とは異なる、月周回アレイ・イメージングにおける決定的な系統誤差であるサブピクセル・エイリアシング問題を特定し、解決した点にある。ピクセル平均法の妥当性を検証することで、著者らは、独特な$uvw$カバレッジや軌道制約にもかかわらず、高品質な全天マップが再構成可能であることを示している。著者らは、結果が理想化されたモデルに基づいていることを明示しており、計器のキャリブレーション、衛星の位置・姿勢誤差、月面反射効果などの系統誤差は本研究では無視されており、後続の論文で扱われる予定であることを認めている。提示された結論は、衛星アレイが設計パラメータを満たしている限り、そのような系統誤差に対して堅牢であると主張されている。