Pulsar searches of Fermi-LAT gamma-ray sources with the MWA

기술 요약: MWA를 이용한 Fermi-LAT 감마선 천체 대상 펄서 탐색

문제점
밀리초 펄서(MSP)는 중력파 검출을 위한 펄서 타이밍 어레이를 포함하여 천체 물리학에서 매우 중요한 도구입니다. 알려진 은하계 MSP의 약 5분의 1은 Fermi Large Area Telescope(Fermi-LAT)가 식별한 미연관 감마선 천체를 대상으로 한 표적 라디오 탐색을 통해 발견되었습니다. 그러나 이러한 탐색은 역사적으로 분산(dispersion) 및 산란(scattering) 효과를 완화하기 위해 300 MHz 이상의 라디오 주파수에서 수행되어 왔습니다. 이는 가파른 라디오 스펙트럼을 가진 MSP가 탐지될 수 있는 300 MHz 미만의 파라미터 공간에 공백을 남깁니다. 저주파 서베이는 차가운 이온화된 성간 매질(IISM)로 인한 색수차 효과, 특히 큰 분산 지연과 채널 내 스미어링(intrachannel smearing)을 보정하는 데 드는 계산 비용 문제로 인해 상당한 어려움을 겪으며, 이는 단주기 펄서에 대한 민감도를 제한합니다. LOFAR(Low Frequency Array) 서베이는 세미 코히어런트 디스퍼션(semi-coherent dedispersion) 전략을 활용하여 135 MHz에서 세 개의 MSP를 성공적으로 발견했으나, MWA(Murchison Widefield Array)를 이용한 유사한 서베이는 수행되지 않았습니다.

방법론
저자들은 남반구 하늘의 MWA Rapid Two-metre(SMART) 서베이의 아카이브 데이터를 사용하여 미연관 Fermi-LAT 소스를 대상으로 하는 새로운 MWA 맞춤형 탐색 파이프라인을 개발했습니다. 방법론은 다음과 같습니다:

• 데이터 소스: 서베이는 154 MHz에서의 SMART 서베이로부터 얻은 20분간의 빔포밍 관측 데이터를 활용하였으며, 308개의 미연관 Fermi-LAT 소스(4FGL-DR4 카탈로그에서 적위 < +25°인 천체)를 대상으로 했습니다.
• 파이프라인 구조: 파이프라인은 계산 비용과 분산 스미어링 사이의 균형을 맞추기 위해 세미 코히어런트 디스퍼션 방식을 채택합니다.
  • 위상 코히어런트 디스퍼션(Phase-Coherent Dedispersion): GPU 가속 CDMT 소프트웨어를 사용하여 복소 전압(complex voltages)을 1.5에서 79.5 cm⁻³ pc 범위의 27개 시도 분산량(DM) 단계(3 cm⁻³ pc 간격)로 코히어런트하게 디스퍼션합니다.
  • 인코히어런트 디스퍼션(Incoherent Dedispersion): 이후 DEDISP 라이브러리를 사용하여 코히어런트 단계 사이의 간극을 메우기 위해 인코히어런트하게 디스퍼션하며, 결과적으로 총 27,000개의 DM 시도를 생성합니다.
  • 가속도 탐색(Acceleration Search): PRESTO의 ACCELSEARCH 루틴을 통해 바이너리 펄서를 탐지하기 위한 푸리에 도메인 가속도 탐색(FDAS)을 수행하며, $z$ 값이 -200에서 +200 사이인 궤도 가속도를 탐색합니다.
  • 후보군 선별(Candidate Sifting): 후보들은 신호 대 잡음비 임계값, 고조파 일관성 및 DM 유효성을 기준으로 필터링됩니다. 최종 후보들은 PREPFOLD를 사용하여 폴딩(folding) 및 최적화됩니다.
• 검증: 파이프라인은 SMART 데이터 내에서 알려진 다섯 개의 MSP(네 개는 단주기 궤도 바이너리)를 블라인드 탐지함으로써, 저주파에서도 스파이더 시스템(spider systems)을 탐지할 수 있는 능력을 확인했습니다.

주요 기여

1. 최대 규모의 저주파 서서베이: 본 연구는 현재까지 미연관 Fermi-LAT 소체를 대상으로 한 가장 큰 규모의 라디오 서베이이며, 300 MHz 미만에서 수행된 두 번째 서베이입니다.
2. 파이프라인 구현: 저자들은 MWA 전압 캡처 시스템(VCS) 데이터에 적합한 세미 코히어런트 디스퍼션 파이프라인을 성공적으로 적응 및 구현하였으며, 이를 통해 기존의 완전 인코히어런트 방식보다 더 깊은 수준의 MSP 탐색을 가능하게 했습니다.
3. 민감도 분석: 본 연구는 1.4 GHz에서의 28% 듀티 사이클을 가진 2 ms 주기 펄서에 대해 약 30220 mJy(1.4 GHz 기준 0.75.2 mJy에 해당)의 최소 탐지 가능 플럭스 밀도($S_{min}$)를 계산하여, 154 MHz에서의 MWA에 대한 상세한 민감도 추정치를 제공합니다.

결과

• 발견 수율: 308개의 소스에 대한 서베이에서 새로운 펄서가 식별되지 않았습니다. 저자들은 이 무효 결과(null result)를 불충분한 민감도 때문으로 돌렸습니다.
• 블라인드 탐지: 파이프라인은 타이드 어레이 빔(tied-array beam)의 사이드로브(sidelobe) 내에 있는 두 개의 알려진 비재생 펄서(PSR J1752−2806 및 PSR J2048−1616)를 성공적으로 탐지하여 파이프라인의 기능을 입증했습니다.
• 민감도 한계: 서베이는 알려진 라디오 밝은 감마선 펄서 인구(산란에 의한 DM 제한 하한선 아래)의 8%에서 44% 사이의 민감도를 가졌습니다. 반면, LOFAR 감마선 서베이는 인구의 86~95%에 대한 민감도를 가졌으며 세 개의 MSP를 발견했습니다. 저자들은 MWA 서베이가 LOFA의 수율과 일치했더라도 기대되는 발견 수가 1개 미만이었을 것이며, 이는 무효 결과와 일치한다고 언급했습니다.
• 블라인드 테스트: 파이프라인 검증을 위해 사용된 다섯 개의 테스트 MSP 중 네 개가 블라인드로 탐지되었습니다. 하나(PSR J0952−0607)는 낮은 신호 대 잡음비와 상당한 분율 스미어링(fractional smearing)으로 인해 탐지되지 않았습니다.

의의 및 향후 전망
본 논문은 파일럿 서베이가 새로운 발견을 이끌어내지는 못했지만, MWA를 이용한 표적 MSP 탐색을 위한 검증된 워크플로우를 구축했다고 결론짓습니다. 저자들은 무효 결과가 방법론의 실패라기보다는 서베이의 민감도 한계와 일치한다고 주장합니다.

이 연구의 의의는 향후 업그레이드를 위한 준비에 있습니다:

• MWA Phase III: 타일 수를 256개로 두 배 늘리는 예정된 Phase III 업그레이드는 순시 민감도를 2배로 높일 것으로 예상됩니다. 이는 동일한 적분 시간(20분) 동안 탐지 가능한 감마선 펄서의 기대 수율을 약 30% 증가시킬 것입니다.
• 실시간 빔포머(Real-Time Beamformer, RTB): RTB의 개발은 소스에 더 가깝게 중심을 둔 표적 관측을 가능하게 하여, 빔 오프셋(beam offset)으로 인한 민감도 손실을 줄이고 더 긴 적분 시간을 가능하게 할 것입니다.
• 적분 시간 연장: 적분 시간을 60분으로 늘리면 기대 수율을 약 60% 증가시킬 수 있지만, 이는 매우 짧은 바이너리 궤도에 대한 민감도를 낮출 수 있습니다.
• 광범위한 응용: 개발된 세미 코히어런트 파이프라인은 초신성 잔해, 구상 성단 및 가파른 스펙트럼을 가진 라디오 소스 탐색에도 적용 가능하며, 이는 향후 SKA-Low 서베이의 토대가 됩니다.

저자들은 산란 제한과 적분 시간이 효과적으로 관리된다면, 낮은 주파수에서 관측하는 것이 가파른 스펙트럼을 가진 펄서들에게 여전히 유리하다는 점을 강조합니다.