[論文レビュー] rPICARD: A CASA-based Calibration Pipeline for VLBI Data. Calibration and imaging of 7 mm VLBA observations of the AGN jet in M87

Michaël Janssen, C. Goddi+9 more|arXiv (Cornell University)|Feb 5, 2019

Radio Astronomy Observations and Technology参考文献 30被引用数 9

ひとこと要約

rPICARD は、S/N に応じた解の時間間隔を採用するグローバルな Schwab-Cotton 位相揃えフィッティングを用いた、CASAbased のオープンソースでエンドツーエンドの校正および画像化パイプラインであり、7 mm VLBA データの校正に成功し、複雑で縁が明るいジェット構造と、10% の輝度を持つ可能性のある反ジェットを明らかにした。これは、小さな傾斜角に対してジェット速度が約 0.5c であることを制限するものである。

ABSTRACT

(Abridged) The CASA software suite, can now reduce very long baseline interferometry (VLBI) data with the recent addition of a fringe fitter. Here, we present the Radboud PIpeline for the Calibration of high Angular Resolution Data (rPICARD), which is an open-source VLBI calibration and imaging pipeline built on top of the CASA framework. The pipeline is capable of reducing data from different VLBI arrays. It can be run non-interactively after only a few non-default input parameters are set and delivers high-quality calibrated data. CPU scalability based on a message-passing interface (MPI) implementation ensures that large bandwidth data from future arrays can be processed within reasonable computing times. Phase calibration is done with a Schwab-Cotton fringe fit algorithm. For the calibration of residual atmospheric effects, optimal solution intervals are determined based on the signal-to-noise ratio (S/N) of the data for each scan. Different solution intervals can be set for different antennas in the same scan to increase the number of detections in the low S/N regime. These novel techniques allow rPICARD to calibrate data from different arrays, including high-frequency and low-sensitivity arrays. The amplitude calibration is based on standard telescope metadata, and a robust algorithm can solve for atmospheric opacity attenuation in the high-frequency regime. Standard CASA tasks are used for CLEAN imaging and self-calibration. In this work we demonstrate the capabilities of rPICARD by calibrating and imaging 7 mm VLBA data of the central radio source in the M87 galaxy. The reconstructed jet image reveals a complex collimation profile and edge-brightened structure. A potential counter-jet is detected that has 10 % of the brightness of the approaching jet. This constrains jet speeds close to the radio core to about half the speed of light for small inclination angles.

研究の動機と目的

• 現代のソフトウェアエコシステムにおいて、自動的でスケーラブルかつユーザーフrndly な VLBI キャリブレーションパイプラインの不足に対処すること。
• CASA のフレームワークを活用して、高周波数、低感度、スパarsな VLBI アレイの校正を堅牢に行えるようにすること。
• 最小限のユーザ入力で科学的利用に適した校正済みデータを提供する、非インタラクティブで自己調整可能なパイプラインを提供すること。
• MPI を用いた CPU スケーラビリティにより、将来的な高帯域幅 VLBI データをサポートすること。
• EHT 観測を含む、多波長・多機器 VLBI 実験における一貫性があり再現可能なデータ還元を促進すること。

提案手法

• 位相校正にグローバルな Schwab-Cotton 位相揃えフィッティングアルゴリズムを採用し、レガシーデータの AIPS 依存ワークフローに代わるものである。
• 各スキャンごとの信号対雑音比（S/N）分析を用いて、各電波望遠鏡ごとに最適な解の時間間隔を決定し、低 S/N 境界でも検出性能を向上させる。
• 標準的な望遠鏡メタデータを用いて振幅校正を実施し、高周波数における大気透過率の減衰を効果的に補正するアルゴリズムを組み込む。
• CLEAN 画像化と自己校正には標準的な CASA タスク（例：tclean）を採用し、互換性と再現性を確保する。
• 大規模帯域幅 VLBI データセット向けに、MPI を用いた並列化を実装し、CPU 使用量を効率的にスケーリングする。
• 新しい polsolve タスクを介して、将来の拡張として基線スタッキング、スペクトルライン校正、高度なストークス偏光校正をサポートする。

実験結果

リサーチクエスチョン

• RQ1 CASA ベースのパイプラインは、インタラクティブ操作なしで 7 mm VLBA データの信頼性の高い高精度な校正を達成できるか？
• RQ2 S/N に応じた解の時間間隔選択は、低感度 VLBI 観測における位相揃え検出と校正精度をどのように向上させるか？
• RQ3 rPICARD は、実宇宙データにおいて、AIPS や Difmap といった既存の VLBI ソフトウェアの結果をどの程度再現できるか？
• RQ4 rPICARD は、複雑な大気効果を伴う高周波数・高ダイナミックレンジの VLBI データの校正をサポートできるか？
• RQ5 rPICARD は、EHT のような複数観測所・多周波数実験において、どの程度の性能を示すか？

主な発見

• rPICARD はデフォルト設定で M87 の中心電波源の 7 mm VLBA データを校正し、科学的利用に適した結果を出力した。
• 再構成されたジェット画像は、複雑なコラimationプロファイルと、先行研究と整合する縁が明るい形状を示した。
• 反ジェットが、アプローチングジェットの約 10% の輝度温度で検出された。
• 反ジェットの検出結果から、小さな傾斜角に対して、核付近のジェット速度が光速の約 0.5 倍であることが制限された。
• rPICARD の結果は AIPS や Difmap の結果と優れた一致を示し、その信頼性と正確性を確認した。
• rPICARD は 2017 年 4 月の 1.3 mm EHT 観測において、3 つの独立した還元経路の一つとして使用され、ブラックホールの影の初画像に貢献した。