펨토헤르츠 중력파 탐색을 위한 펄서 파라미터 드리프트 분석

초록

이 논문은 피코헤르츠(10⁻¹² Hz) 대역의 연속 중력파를 탐지하기 위해 이진 궤도 주기 변화율(˙P_b)과 스핀 주기 2차 미분(¨P)이라는 두 민감한 펄서 타이밍 파라미터의 드리프트를 이용한 최초의 베이지안 검색 방법을 제시한다. 기존 PT A보다 두 배 이상의 펄서를 사용해 감도 10배 향상에 성공했으며, 현재 데이터에서는 신호를 찾지 못했지만 차세대 SKA 관측으로 이 대역의 중력파를 실질적으로 탐지할 가능성을 제시한다.

상세 분석

본 연구는 피코헤르츠(PHz) 대역, 즉 10⁻¹⁸ ~ 10⁻¹⁶ Hz 사이의 초저주파 중력파가 기존 펄서 타이밍 배열(PTA) 분석에 직접적으로 나타나지 못한다는 점에 착안한다. 이 주파수에서는 파동 주기가 수천 년에 이르므로, 수십 년에 불과한 관측 기간 동안은 파동의 위상 변화가 거의 정적이다. 저주파 파동은 펄서와 태양계 중심(Solar System Barycenter) 사이의 유효 속도와 가속도를 미세하게 변형시키며, 이는 관측 가능한 파라미터인 이진 궤도 주기 변화율 ˙P_b와 스핀 주기 2차 미분 ¨P에 누적된 형태로 나타난다. 논문은 이러한 효과를 정량화하기 위해 GW에 의해 유도되는 속도 v_GW, 가속도 a_GW, 그리고 ‘jerk’ j_GW를 도출하고, 각각을 ˙P_b와 ¨P에 연결시킨다(식 3·4).

베이지안 프레임워크는 7개의 신호 파라미터(진폭 h₀, 주파수 f_GW, 그리고 천구좌표·위상·편광을 나타내는 Θ={α,δ,i,ϕ,Ψ₀})를 포함한다. 관측값 y_GW,p={a_GW,p, j_GW,p}와 각 펄서별 불확도 σ_p를 이용해 다변량 정규분포 형태의 우도 함수를 구성하고, 사전분포는 대부분 균등(α,δ 등)이며 h₀에 대해서는 로그 균등과 선형 균등을 상황에 따라 교차 적용한다. MCMC 샘플링을 통해 사후분포를 추정하고, 베이즈 팩터(BF)를 이용해 검출 여부를 판단한다.

펄서 샘플은 기존 14개에서 30개(˙P_b)와 29개(¨P)로 확대했으며, 각 펄서의 은하 가속도와 시프코프스키 효과를 독립적으로 보정한다. 특히, 파동 파장보다 펄서 거리 측정 오차가 작아야 하는 조건 f·Δd_p≤1을 만족하는 절반 이상의 펄서를 선별해 분석에 포함시켰다. 결과적으로, 현재 데이터에서는 BF가 검출 임계값 이하이며, h₀에 대한 상한이 기존 연구보다 약 10배 낮아졌다.

앞으로 SKA가 제공할 수 있는 수십 년에 걸친 장기 관측과 수백 개의 고정밀 MSP 타이밍은 감도를 현재보다 1~2계량 이상 향상시킬 것으로 예측된다. 이는 SMBHB(초대질량 블랙홀 이진) 초기 단계, 즉 1 pc 이하의 분리에서 방출되는 피코헤르츠 연속 파동을 직접 탐지할 수 있는 첫 번째 기회를 제공한다. 또한, 코스믹 스트링 네트워크, 초기에 일어난 1차 상전이, 그리고 인플레이션 전후의 프리히팅 등 새로운 물리학 시그널도 이 대역에서 검출 가능성이 열려 있다.