[논문 리뷰] 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명|arXiv (Cornell University)|2019. 02. 05.

Radio Astronomy Observations and Technology참고 문헌 30인용 수 9

한 줄 요약

rPICARD는 CASA 기반의 오픈소스 파이프라인으로, 글로벌 Schwab-Cotton 간섭계 피팅과 S/N 적응형 해상도 간격을 사용하여 VLBI 데이터의 종단 간보정 및 이미징을 수행한다. 이는 강력한 위상 간섭계 보정을 가능하게 하며, M87의 AGN 제트의 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.

연구 동기 및 목표

• 현대 소프트웨어 생태계에서 자동화, 확장성, 사용자 우호성에 뒷받침되는 VLBI 간섭계 보정 파이프라인의 부족을 해결한다.
• CASA의 프레임워크를 활용하여 고주파수, 저감도, 희박한 VLBI 어레이의 강력한 간섭계 보정을 가능하게 한다.
• 사용자 입력 최소화로 과학용으로 바로 사용 가능한 간섭계 보정 데이터를 제공하는 비상호작용형, 자가조정 파이프라인을 제공한다.
• MPI 기반의 CPU 확장성을 통해 향후 고대역폭 VLBI 데이터를 지원한다.
• 다중 파장 및 다중 장비 VLBI 실험(예: EHT 관측 포함)에 대한 일관되고 재현 가능한 데이터 감소를 촉진한다.

제안 방법

• 위상 간섭계 보정을 위해 글로벌 Schwab-Cotton 간섭계 피팅 알고리즘을 사용하여 이전의 AIPS 기반 워크플로우를 대체한다.
• 각 스캔에 대해 신호 대 잡음비(S/N) 분석을 수행하여 각 천체망원경에 최적의 해상도 간격을 결정함으로써, 저감도 영역에서의 탐지 성능을 향상시킨다.
• 표준 천체망원경 메타데이터를 사용하여 진폭 간섭계 보정을 수행하고, 고주파수에서 대기 투과율 감쇠를 해결하기 위한 강력한 알고리즘을 포함한다.
• 표준 CASA 작업(예: tclean)을 사용하여 CLEAN 이미징 및 자기 간섭계 보정을 수행함으로써 호환성과 재현 가능성을 확보한다.
• 대규모 대역폭 VLBI 데이터셋에 대한 효율적인 CPU 사용을 위해 MPI 기반 병렬 처리를 구현한다.
• 새로운 polsolve 작업을 통해 향후 확장 기능(예: 기저선 스택킹, 스펙트럼선 간섭계 보정, 고급 편광 간섭계 보정)을 지원한다.

실험 결과

연구 질문

• RQ1 CASA 기반 파이프라인이 상호작용 없이도 7 mm VLBA 데이터에 대해 신뢰할 수 있고 고해상도 간섭계 보정을 수행할 수 있는가?
• RQ2 S/N 적응형 해상도 간격 선택 방식이 저감도 VLBI 관측에서 간섭계 탐지 및 간섭계 보정 정확도를 어떻게 향상시키는가?
• RQ3 rPICARD가 실제 천체 관측 데이터에서 기존의 VLBI 소프트웨어(AIPS, Difmap 등)의 결과를 어느 정도 재현할 수 있는가?
• RQ4 rPICARD는 복잡한 대기 효과를 동반한 고주파수, 고다이나믹 레인지 VLBI 데이터의 간섭계 보정을 지원할 수 있는가?
• RQ5 rPICARD는 EHT와 같은 다중 관측소, 다중 주파수 실험에서 어떻게 성능을 발휘하는가?

주요 결과

• rPICARD는 기본 설정을 사용하여 M87 중심 라디오 소스의 7 mm VLBA 데이터를 성공적으로 간섭계 보정하고 과학용 결과를 제공하였다.
• 재구성된 제트 이미지에서는 복잡한 콘자구조와 이전 연구와 일치하는 가장자리 밝은 형태를 보였다.
• 약 10%의 밝기로 접근 제트의 반대편에 해당하는 잠재적 반대 제트가 탐지되었다.
• 반대 제트 탐지 결과는 작은 기울기 각도에서 핵 부근의 제트 속도가 약 0.5c로 제약됨을 시사한다.
• rPICARD의 결과는 AIPS 및 Difmap의 결과와 뛰어난 일치를 보이며, 그 신뢰성과 정확성을 확인하였다.
• rPICARD는 1.3 mm EHT 2017년 4월 관측에서 세 가지 독립적인 감소 경로 중 하나로 사용되었으며, 블랙홀 그림자 첫 번째 이미지 제작에 기여하였다.