Bayesian polarization calibration and imaging in very long baseline interferometry

宇宙が、絶えず地球に降り注いでいる巨大で目に見えない電波の海であると想像してみてください。この海の中に何が隠れているのかを見るために、天文学者は普通のカメラは使いません。その代わりに、世界中に散らばった電波望遠鏡を使って、巨大な「浮遊する網」を作り上げます。この手法は、超長基線電波干渉法（VLBI）と呼ばれます。これらの望遠鏡を連結させることで、地球サイズの仮想的な望遠鏡を作り出し、数千キロ先から見たコインが収まるほど微細な詳細を見ることができるのです。しかし、これらの網を通して宇宙を覗き見ることは困難です。私たちの望遠鏡の「レンズ」である電波皿自体が不完全なのです。そこには小さな「漏れ」があり、ある種類の偏波（電波の振動方向のこと）からの信号が別の種類へと漏れ出してしまいます。これは、耳が少し詰まった状態でオーケストラを聴こうとしているようなものです。ドラムの音がバイオリンのセクションに混じり込んでしまうのです。もし、これらの漏れを修正しなければ、構築される宇宙の画像はぼやけ、歪んでしまいます。これは、ブラックホールや銀河の周囲にある磁場を研究する際に特に困難です。なぜなら、私たちが探している信号は極めて微弱で、メインの電波の輝きの1000分の1以下の強さであることが多いからです。

本論文は、これらの漏れを修正し、より鮮明な画像を作り上げるための、よりスマートな新しい方法を紹介しています。従来の、パズルのピースがどこに合うかを推測して当てはめていくような古い方法ではなく、著者たちは「ベイズ的」なアプローチを採用しています。これは、単に推測するのではなく、あらゆる可能な組み合わせを考慮し、各シナリオの確率を重み付けし、最も可能性の高い画像を見つけ出すと同時に、個々のピースに対する確信度も提示してくれる、非常に賢い探偵のようなものです。著者たちは、この新しい探偵の手法を、クエーサーと呼ばれる3C273と、ブラザーと呼ばれるOJ287という2つの有名な天体の実データを用いてテストしました。その結果、彼らの手法は従来の方法よりも細かい詳細を見ることができ、以前は見ることができなかったこれらの天体の磁場の複雑な構造を明らかにできることが分かりました。決定的なのは、この新しい手法は単に画像を与えるだけでなく、すべてのピクセルに対して「信頼スコア」を与えることで、その画像をどの程度信頼できるかを正確に教えてくれる点です。

旧来の手法 vs 新しい探偵

長い間、天文学者はラジオデータを画像に変換するために「CLEAN」と呼ばれる手法を使用してきました。これは、数枚のぼやけた写真に基づいて雲の絵を描こうとしているようなものです。CLEAN法は、ぼやけた写真の中から最も明るい点を見つけ、そこに完璧な点を描き、その点を写真から差し引いて残ったものを見るという作業を行います。これを何度も繰り返します。それは、玉ねぎの皮を一層ずつ剥いていく作業に似ています。問題は、この方法には多くの人間の推測が伴うことです。コンピュータに対して、どこを見ればよいか、どのように点を差し引くべきかを指示しなければなりません。もし間違いを犯したり、あるいは雲が単なるいくつかの点ではなく複雑に渦巻く形状であった場合、最終的な画像には奇妙なアーティファクトが生じたり、詳細が欠落したりすることがあります。さらに、旧来の方法は、望遠鏡の「漏れ」（偏波エラー）を静的で単純なものとして扱っていますが、現代の望遠鏡から得られる複雑で高速なデータにおいては、それは事実ではありません。

論文の著者たちは、この古い方法が最適ではないと主張しています。それは熟練したユーザーによる手動の操作を多く必要とし、結果の不確実性を全く示してくれないからです。それは、降水確率が10%なのか90%なのかを教えずに、「雨が降ります」とだけ告げる天気予報のようなものです。

ベイズ的解決策：確率論的な探偵

本論文で提示されている新しい手法は、「ベイズ推論」と呼ばれるものを使用しています。玉ねぎを剥く比喩を別のものに変えてみましょう。暗い部屋のレイアウトを把握するために、壁に向かってボールを投げ、その反響を聞いているところを想像してください。旧来の方法では、ボールを数回投げ、壁の位置を推測して地図を描きます。しかし、新しい方法は、聞こえてきた反響を生み出す可能性のある「あらゆる」部屋のレイアウトを検討します。そして、部屋のあらゆる地点に対して「確率マップ」を計算します。

論文の言葉を使えば、彼らは「事後分布」を探求しています。これは、単に一つの「最良の」答えを選ぶのではなく、あり得るすべての答えの範囲（画像、望遠鏡の漏れ、および信号ゲイン）を一度に探索していることを意味します。彼らは、resolve というソフトウェアツールを使用してこれを行います。このツールは「変分推論」という技術を使用しており、これは、一つ一つの可能性をすべてチェックすることなく、確率マップの最も可能性の高い形状をコンピュータが迅速に見つけ出すための数学的なショートカットです。

得られた知見：より鮮明な画像と誠実な不確実性

チームは、この新しい探偵を2つの実際の宇宙ターゲットでテストしました。

第一に、15 GHzで観測されたクエーサー（非常に明るい活動銀河）である 3C273 を調査しました。彼らは新しい画像と従来のCLEAN画像を比較しました。新しい手法は、より鮮明で詳細な画像を作成しました。例えば、クエーサーから噴き出している物質の「ジェット」は、新しい画像ではより細く、明確に見えました。旧来の手法は、空が単純な点の集まりであると仮定していたため、これらの微細な詳細をぼやけさせてしまう傾向がありました。しかし、新しい手法は、空が連続的で複雑な景観であることを理解していました。また、観測に使用された8つの異なる周波数帯それぞれについて、「漏れ」（D項）を計算することにも成功しました。これにより、漏れが各帯域で異なっていることを示しましたが、これは旧来の手法ではしばら見落とされがちな詳細でした。

第二に、86 GHzで観測されたブラザー（ジェットが我々の方を向いている銀河）である OJ287 を調査しました。これは、信号がより弱く、大気がより多くの干渉を引き起こすため、さらに困難な課題です。ここで、彼らは自らの結果を、旧来の手法だけでなく、ehtim と呼ばれるもう一つの現代的な手法とも比較しました。その結果、彼らの新しい手法は、旧来の手法では見ることができなかった、湾曲したジェットのような複雑な構造を再構成できることが分かりました。また、新しい手法は「標準偏差」マップ、つまり不確実性のマップを提供しました。これは、データの質が悪いために画像がぼやけている場所と、鮮明な場所を正確に示すものです。これは科学において、何を知っているかと同じくらい、何を知らないかを知ることも重要であるため、非常に大きな利点となります。

ゲームのルール

著者たちは、新しい手法において物理学のルールに従うよう細心の注意を払いました。彼らは「偏波制約」というルールを課しました。これは、偏光した光の総量は全光量を超えてはならないというルールです。旧来の手法では、計算が複雑になり、偏光した光が全光量よりも強くなるという、現実にはあり得ない「非物理的」な画像が生成されてしまうことがありました。新しい手法は、このルールを数学の中に直接組み込むことで、画像が物理的に現実的であることを保証しています。

また、彼らの手法は「不均一な」データ、つまりサイズや感度が異なる様々な望遠鏡からのデータを扱うことができることも分かりました。旧来の手法はこれに苦戦することが多く、天文学者が多くの悪いデータを捨て去る（「フラグを立てる」）必要がありました。しかし、新しい手法は、データを適切に重み付けすることができ、データセットの一部を手動で削除することなく、より多くの情報を保持してより良い画像を作成することができました。

結論

この論文は、天文学におけるあらゆる問題を解決したと主張しているわけではありません。彼らがテストした特定のデータセット（3C273およびOJ287）において、彼らのベイズ的手法が高品質で詳細な画像を作成し、不確実性の正直な尺度を提供することを示しています。これは、手動の推測に頼る方法から、完全に自動化された確率論的なアプローチへと移行することで、宇宙の最も極端な環境における磁場や構造をより鮮明に見ることができるようになることを示唆しています。著者たちは、このパイプラインが将来のさらに強力な電波アレイで使用できる準備ができていると結論付けており、私たちが単に宇宙を見るだけでなく、自分たちが何を見ているかについてどれほど確信を持っているかを理解できる、高忠実度の偏光イメージングの新時代を約束しています。