Spectral Data-cube Cleaning for CCAT Deep Spectroscopic Survey. I. Effect of correlated noise and filtering on the power spectrum

満員のスタジアムの中で、たった一つの微かな囁き声を聞こうとしている場面を想像してみてください。それは、初期宇宙の地図を作ろうとしている天文学者たちが直面している課題そのものです。彼らは、数十億年前の「再電離の時代」として知られる時期に形成された古代銀河から放出される、特定の種類の「光」（電波）を探しています。これらの囁き声を捉えるために、彼らは超高感度な検出器を備えた巨大な望遠鏡を使用します。しかし、私たちの頭上にある大気は、水蒸気が絶えずパチパチと音を立て、うねる海のように騒がしく、宇宙からの信号をかき消してしまいます。この「ノイズ」はランダムではなく、時間の経過とともにゆっくりと変化する低周波の唸りであり、実際の信号から分離することを非常に困難にしています。これを解決するために、科学者たちは「輝度強度マッピング（Line-Intensity Mapping: LIM）」と呼ばれる手法を用います。遠方の個々の銀河を鮮明に撮影しようとする（それらはあまりに暗すぎて明確に見ることができません）代わりに、空の一画に存在するすべての銀河の総和としての輝きを測定するのです。これは、一人ひとりの声を聞こうとするのではなく、群衆の集団的なハミングを聞こうとするようなものです。目標は、宇宙の構造の3Dマップを作成し、宇宙が若かった頃に星やガスがどのように分布していたかを明らかにすることです。

本論文では、そのノイズの多いデータをクリーンアップするという極めて難しい問題に取り組んでいます。著者らは、Fred Young Submillimeter Telescope (FYST) に搭載される予定の次世代観測装置「EoR-Spec」による観測をシミュレーションしました。彼らは、「フィルター・アンド・ビン（Filter-and-Bin: F&B）」と呼ぶ特定のクリーニングツールのセットが、宇宙の信号を誤って捨ててしまうことなく、大気のノイズを効果的に取り除くことができるかどうかを確認したいと考えました。データを泥の混じった小川だと考えてください。このパイプラインは、泥（ノイズ）を取り除きつつ、金塊（信号）を保持するように設計された一連のふるいやフィルターです。チームは、このパイプラインをテストするために数千回のコンピュータ・シミュレーションを実行しました。その結果、F&Bパイプラインは低周波の大気の唸りを除去することに非常に効果的であり、ノイズを約10,000分の1に低減できることが分かりました。しかし、強力なフィルターには常に副作用があります。それは、データ内の最も大きく広がったパターンを平滑化してしまう傾向があることです。結果は、この手法が宇宙の「粒状」の詳細（個々の銀河がノイズに寄与する部分）を検出するには完璧である一方で、滑らかで大規模な構造を復元することには苦労することを示唆しています。しかし、特定の目的である「ショットノイズ」領域（ガスや星の混沌とした小規模な集積）を測定するという点においては、このパイプラインは素晴らしく機能しており、今後のサーベイによって初期宇宙の炭素と一酸化炭素の複合的な輝きを検出できる可能性を約束しています。

宇宙の囁きの物語

設定：騒がしい宇宙
あなたが微かなラジオ局を聞こうとしているが、ラジオが置かれている部屋の空気自体が静電気でブーンと鳴っている場面を想像してください。天文学において、この静電気は地球の大気、特に水蒸量から発生します。望遠鏡が空を見上げるとき、単に星を見ているだけでなく、大気の「呼吸」を見ています。これは変動し、$1/f$ ノイズとして知られる低周波の唸りを生み出します。このノイズは、天文学者が測定しようとしているゆっくりとした大規模な変化と似ているため、非常に厄介です。もしこれを掃除しなければ、あなたの宇宙マップは、大気の天候によるぼやけた塊になってしまいます。

ツール：フィルター・アンド・ビン・パイプライン
これを修正するために、著者らは「フィルター・アンド・ビン（F&B）」と呼ぶデジタル・クリーニング・クルーを設計しました。彼らは、仮想的な望遠鏡として機能する大規模なコンピュータ・シミュレーションを用いて、このクルーをテストしました。彼らは、探したい信号（炭素と一酸化炭素ガスの光）と、厄介な大気ノイズを含む架空の宇宙をシミュレーションに投入しました。そして、データを4つの明確なクリーニング・ステップに通しました：

1. 多項式デトレンド（Polynomial De-trending）： まず、データのゆっくりとした緩やかなドリフトを取り除きました。これは、グラグラするテーブルの脚を水平にするような作業です。
2. スキャン同期補正（Scan-Synchronous Correction）： 次に、望遠鏡が前後に動くことによって生じるパターンを修正しました。これは、首を振るたびに静電気の音が大きくなることに気づき、その特定のパターンを差し引くようなものです。
3. コモンモード減算（Common-Mode Subtraction）： 大気のノイズはすべての検出器に同時に影響を与えるため、彼らは「平均的なノイズ」を計算し、全員のデータから差し引きました。これは、合唱団の指揮者が、全員に対して背景の同じ音を止めるよう求めるようなものです。
4. PCAフィルタリング（PCA Filtering）： 最後に、主成分分析（PCA）と呼ばれるスマートな数学的手法を用いて、他のステップで見逃した残りの執拗なパターンを見つけ出し、除去しました。彼らは上位3つのノイズの「パターン」を取り除き、残ったのはほとんど無害なランダムな静電気だけでした。

結果：キャッチはあるがクリーンな信号
シミュレーションによれば、このクリーニング・クルーは大成功を収めました。最低周波数（最も遅く、最も厄介な唸り）において、パイプラインはノイズを約4桁（つまり10,000倍静かに！）減少させました。残ったデータは、ほぼ完全にランダムなホワイトノイズのように見え、これはクリーンな測定にまさに求められる状態です。

しかし、トレードオフが存在します。ふるいが砂を通すが小石は止めるように、このパイプラインは小さく詳細な信号は通しますが、非常に大きく滑らかな波はブロックしてしまいます。著者らは、小規模な詳細（波長数 $k$ が 0.5 Mpc$^{-1}$ より大きい場合）については、パイプラインが信号の80%以上を維持することを発見しました。しかし、最も大きく広がった構造（$k$ が 0.1 Mpc$^{-1}$ 未満の場合）については、パイプラインが信号を抑制しすぎてしまい、残るのは20%未満となります。

これが将来にとって何を意味するか
本論文は、このF&BパイプラインがCCATディープ・スペクトロスコピック・サーベイ（DSS）が計画している任務に適したツールであると結論付けています。このサーベイは、赤方偏移が 3.5 から 8.0 の間（宇宙がもっと若かった時期）の宇宙をマッピングすることを目的としています。シミュレーションにおいて、単一の望遠鏡モジュールが 50% の効率で動作する場合でも、小規模なスケールでの信号を検出できることが分かりました。もし2つのモジュールを 100% の効率で使用すれば、彼らが探しているすべてのスケールにおいて信号を検出できます。

著者らは、このパイプラインは「ショットノイズ」領域（ガスや星の混沌とした小規模な集積）には適していますが、宇宙の滑らかな大規模構造を測定するには最適ではないことを注意深く述べています。そのためには、将来さらに高度な手法が必要になるでしょう。しかし現時点では、このパイプラインは、嵐をフィルタリングして宇宙の物語がようやく聞こえるようにするための、確実で信頼できる方法を彼らに提供しているのです。