미래 우주γ선 펄서 타이밍 배열: 은하계 MSP 탐사와 중력파 탐지의 새로운 지평

초록

본 논문은 현재 Fermi‑LAT이 관측한 약 100개의 밀리초 펄서(MSP)를 기반으로, 차세대 γ선 관측소가 탐지할 수 있는 MSP 수와 그에 따른 나노헤르츠(nHz) 중력파(GW) 감도 를 시뮬레이션한다. 0.1–5 GeV 에 최적화된 에너지 대역을 갖는 기구는 10³–10⁴개의 새로운 γ‑ray MSP를 발견하고, 현재 라디오 PTA와 동등하거나 그 이상으로 GW 배경을 탐지할 수 있음을 보여준다. 또한, 은하 중심 버즈(GCE)와 연관된 버즈 MSP 집단을 구분할 수 있는 가능성도 제시한다.

상세 분석

이 연구는 세 가지 MSP 인구 모델(S1, S2, S3)을 구축하고, 각각에 대해 스핀다운 파워와 자기장에 기반한 γ‑ray 스펙트럼을 고정밀로 예측하는 파이프라인을 개발하였다. S1은 순수 디스크 MSP만을 포함하고, S2와 S3는 디스크에 추가해 버즈(중심) MSP를 삽입한다. 버즈 MSP는 두 가지 형성 시나리오(백색왜성의 붕괴에 의한 AIC와 전통적인 디스크와 동일한 특성을 갖는 경우)로 나뉘며, 각각의 공간 분포는 GCE를 설명하는 r⁻² 프로파일을 따른다.

시뮬레이션 결과, γ‑ray 스펙트럼은 일반적인 파워‑로우와 지수 절단을 갖는 형태이며, 스핀다운 파워 ˙E가 클수록 높은 γ‑ray 효율을 보인다. 특히 0.1–0.3 GeV에서의 각도 해상도가 향상될수록 배경 대비 신호(S/N)가 크게 증가한다는 점을 강조한다. 이는 기존 Fermi‑LAT이 100 MeV 이하에서 감도 제한을 겪는 이유와 일치한다.

다음으로, 다양한 가상의 γ‑ray 관측소(예: 대형 위성형 전자‑포톤 탐지기, 컴프톤/MeV 대역 탐지기)를 가정하고, 각 기구의 유효 면적, 에너지 해상도, 포인트 소스 감도 등을 파라미터화하였다. GeV‑대역(0.1–5 GeV)에서는 10³–10⁴개의 MSP를 검출할 수 있으며, 이들 중 약 30–40 %는 현재 라디오 PTA에 포함되지 않은 새로운 대상이다. 검출된 MSP들의 TOA(도착 시간) 정밀도는 평균 0.5 µs 수준으로, 이는 현재 라디오 PTA의 최상위 MSP와 동등하거나 약간 우수하다.

GW 감도 평가에서는 Hellings‑Downs 각도 상관 함수를 이용해 전체 배열의 스티프니스(stiffness)를 계산하였다. 결과는 A_gwb ≈ 2 × 10⁻¹⁵(현재 라디오 PTA가 보고한 값) 수준의 stochastic GW 배경을 5년 내에 3σ 이상 검출 가능함을 보여준다. 특히, 버즈 MSP가 추가될 경우, GCE와 연관된 GW 신호(예: 중앙 초대질량 블랙홀의 저주파 파동)와의 상관관계를 구분할 수 있는 통계적 파워가 크게 향상된다.

컴프톤/MeV 대역에서는 MSP 스펙트럼이 급격히 감소하고, 배경(특히 Galactic diffuse emission) 대비 S/N이 0.1 이하로 떨어져 실질적인 PTA 구축이 어려움이 확인되었다. 따라서 차세대 PTA는 GeV 대역에 초점을 맞추는 것이 최적임을 결론짓는다.

마지막으로, γ‑ray PTA와 라디오 PTA의 시너지 효과를 논의한다. γ‑ray PTA는 전파 전파 매질(DM) 변동에 전혀 영향을 받지 않으며, 관측 공백이 거의 없고, 동일한 기구와 데이터 파이프라인으로 전천후 모니터링이 가능하다. 이는 라디오 PTA의 시스템atics(예: 전파 간섭, 관측 스케줄 불균형)를 보완해 전체 배열의 민감도를 20 % 이상 향상시킬 수 있다.