NICER X‑레이 펄서 타이밍으로 보는 나노헤르츠 중력파 배경 제한

초록

NICER를 이용해 6개의 밀리초 펄서를 6년간 관측한 결과, 나노헤르츠 대역의 중력파 배경(GWB)에 대한 95% 상한을 log₁₀A₍gwb₎ < ‑13.4 로 설정하였다. Hellings‑Downs 각도 상관관계는 S = 2.5 수준으로 약한 신호가 보였지만 통계적으로 유의미하지 않다. 현재 라디오·γ‑레이 PTA보다 약간 약하지만, X‑레이 타이밍이 GWB 탐색에 실현 가능함을 입증한다.

상세 분석

본 연구는 NICER(Neutron Star Interior Composition Explorer)의 X‑레이 펄서 타이밍 데이터를 활용해 나노헤르츠(GW) 대역의 stochastic GWB를 탐색한 최초 사례 중 하나이다. 분석 대상은 PSR J1939+2134, J1824‑2452A, J0437‑4715, J0030+0451, J0218+4232, J2124‑3358 등 6개의 밀리초 펄서이며, 2017년 6월부터 2023년 9월까지 약 6년간 누적 관측 시간을 확보하였다.

데이터 전처리는 기존 NICER 파이프라인에 따라 광자 이벤트를 필터링하고, 국제 펄서 타이밍 배열(IPTA) DR2 에피머리스를 이용해 태양계 중심(Solar System Barycenter) 기준으로 TOA(Time of Arrival)를 보정하였다. 이후 30일 구간별로 5×10⁴–2×10⁵개의 광자를 모아 평균 펄스 프로파일을 생성하고, 표준 프로파일과의 교차상관을 통해 TOA를 추출하였다. TOA 오차는 프로파일의 Gaussian 샘플링으로 추정하였다.

노이즈 모델링은 두 가지 구성요소로 나뉜다. (1) White noise: 측정 오차 σ_TOA에 대한 EFAC·EQUAD 파라미터를 도입해 통계적 및 시스템적 백그라운드 노이즈를 보정하였다. (2) Red noise: 각 펄서마다 파워 스펙트럼 P_I(f)=A_I²/(12π²)(f/f_c)^(-γ_I)·yr⁻³ 형태의 적색 잡음을 가정하고, γ_I를 2–7 범위의 사전분포로 제한하였다. 전체 공분산 행렬 C_n은 모든 펄서의 white·red noise 행렬을 블록 대각선 형태로 결합하였다.

GWB 모델은 SMBH 바이너리들의 집합적 기여를 가정하고, 스펙트럼 지수 γ_gwb=13/3(=4.33…)을 고정하였다. 공통 적색 잡음(CRN) 형태의 GWB는 Hellings‑Downs(HD) 각도 상관관계 Γ_IJ(θ_IJ)와 결합해 C_gwb=Γ·ṼC_gwb 로 구성하였다. Bayesian 프레임워크는 Enterprise와 Eryn을 이용해 MCMC 샘플링을 수행했으며, 파라미터는 (i) 각 펄서의 red noise A_I,γ_I, (ii) 공통 GWB 진폭 A_gwb, (iii) 고정된 white noise·타이밍 모델 파라미터(사전값 고정) 로 설정하였다. 사전분포는 기존 라디오·γ‑레이 PTA 결과를 반영한 정규·균등분포를 사용하였다.

결과적으로 로그 진폭 로그₁₀A_gwb의 사후분포는 -14 부근에 약간의 피크를 보였지만, Bayes factor ΔlogZ=1.47으로 GWB 존재를 뒷받침하기엔 부족했다. 95% 신뢰구간 상한은 log₁₀A_gwb < ‑13.0(논문에서는 -13.4) 로 제시되었다. HD 상관관계 검증을 위해 펄서 쌍 간 교차상관 계수를 구하고, Jenet et al. (2004) 방식의 통계량 S=2.5를 얻었다. 이는 약 98.8% 신뢰수준이지만 3σ(≈3) 미만이므로 검출이라기보다 우연 가능성이 큰 수준이다. Bayes factor(HD 포함 vs 미포함) ΔlogZ=0.369 역시 무시할 정도다.

비교 분석에서는 현재 라디오 PTA(NANOGrav, PPTA, CPTA 등)와 γ‑레이 PTA(Fermi‑LAT)보다 약 1~3배 느슨한 제한을 보였으며, 이는 관측 시간(T)과 펄서 수(N)의 제한, 그리고 X‑레이 TOA 정밀도가 아직 라디오에 비해 낮은 점이 원인이다. 그러나 X‑레이는 전파 전파 매질(DM) 변동에 민감하지 않아 장기적인 시스템atics를 크게 줄일 수 있다는 장점이 있다. 향후 NICER와 같은 고정밀 X‑레이 타이밍 미션이 늘어나고, 더 많은 MSP와 긴 관측 기간을 확보한다면, 현재 제한을 1–2 옥타브 정도 개선할 수 있을 것으로 전망된다. 또한, 차세대 X‑레이 관측소(Lynx, STROBE‑X 등)와 결합하면, GWB 탐지에 필요한 민감도에 도달할 가능성이 있다.

결론적으로, 본 연구는 NICER X‑레이 타이밍이 GWB 탐색에 실용적임을 최초로 입증했으며, 현재는 제한적인 감도지만 향후 미션 설계와 데이터 축적 전략에 중요한 기준점을 제공한다.