핫스팟 중력파 탐지를 위한 교차상관법 연구

초록

본 논문은 교차상관법을 이용해 몇 개의 화소에 걸쳐 있는 작은 천구 영역, 즉 ‘핫스팟’에서 발생하는 확률적 중력파 배경을 탐지하는 방법을 제시한다. 대표적인 목표물로 베르길 군집을 선택하고, LIGO·Virgo·LCGT·ET 등 차세대 탐지기 쌍의 기준선을 이용해 연간 수준의 관측 시간으로 신호‑대‑잡음비(SNR)≈1을 달성할 수 있음을 보인다. 다중 기준선 결합 시 SNR이 더욱 향상된다.

상세 분석

이 연구는 기존의 전천구 스테레오스코픽 맵핑과 회전 중성자별(펄서) 타깃 탐색에 사용된 교차상관 기법을, ‘핫스팟’이라 부르는 소규모 천구 영역에 적용한다는 점에서 새롭다. 화소의 크기는 두 탐지기 사이의 기준선 길이와 탐지기 대역폭에 의해 결정되는 회절 한계(Δθ≈c/(2πfL))에 의해 정의되며, 실제로는 몇 도 정도의 각도 범위가 된다. 베르길 군집은 약 10¹¹개의 중성자별을 포함하고, 그 중 극히 일부가 밀리초 펄서 형태로 고주파(≈100 Hz–1 kHz) 중력파를 지속 방출한다고 가정한다.

논문은 먼저 교차상관 통계량 S를 정의하고, 각 시간 구간 Δt(32–192 s)마다 단일 화소에 대한 신호와 잡음의 평균·분산을 계산한다. 최적 필터 Q는 신호 스펙트럼 H(f)와 두 탐지기의 잡음 PSD P₁, P₂, 그리고 방향 의존 겹침 감소 함수 γ(t,f,Ω̂)를 이용해 Q∝H·γ*/(P₁P₂) 형태로 도출된다. 이때 γ는 탐지기 안테나 패턴과 두 탐지기 사이의 위상 차이를 포함한다.

SNR은 연속적 한계에서 ρ≈2 √(T Δf) · ⟨H²/(P₁P₂)⟩¹ᐟ² · ⟨Γ²⟩¹ᐟ² 로 간단히 표현된다. 여기서 Γ는 하루 평균의 방향 겹침 함수이며, 표 I에 제시된 바와 같이 LIGO-Livingston–LIGO-Hanford, LIGO-Livingston–AIGO 등 특정 기준선이 가장 큰 ⟨Γ²⟩¹ᐟ² 값을 갖는다.

핫스팟 신호 스펙트럼 H(f)는 개별 중성자별의 중력파 진폭 h∝ε I R⁻¹ f²와 별 개수 분포 N(f) 를 곱해 H²∝ε² I² R⁻² f⁴ N(f) 로 모델링한다. 논문은 은하 내 밀리초 펄서 수를 4 × 10⁴개(주파수 50 Hz 이상)로 추정하고, 베르길 군집에 약 10³개의 은하가 존재하므로 N_high≈4 × 10⁷, N_low≈10¹¹을 얻는다. 이 분포를 로그-정규 모델로 근사하고, 주파수 대역 100 Hz–1 kHz에서 평균적으로 화소당 10~100개의 펄서가 겹쳐 있어 확률적 배경으로 간주한다.

시뮬레이션 결과는 고감도 차세대 탐지기(특히 ET)의 PSD와 비교했을 때, 1년 관측 시 ρ≈1–3 수준을 달성한다는 것을 보여준다. 다중 기준선(예: 3~4개의 탐지기 조합) 사용 시 ρ는 √N_baseline 만큼 향상되어, 실질적인 검출 가능성이 크게 높아진다. 또한, 베르길 군집 외에도 은하단, 초대질량 블랙홀 주변 등 다른 고밀도 천체 집합에도 동일한 방법을 적용할 수 있음을 제시한다.

핵심적인 한계점은 (1) 중성자별의 비대칭성 파라미터 ε와 회전 주파수 분포에 대한 불확실성, (2) 실제 탐지기 잡음이 시간에 따라 비정상적일 수 있다는 점, (3) 화소 크기가 기준선·대역폭에 의존하므로, 짧은 기준선(예: LIGO–Virgo)에서는 화소가 넓어 배경 잡음이 증가한다는 점이다. 그러나 이러한 제약에도 불구하고, 교차상관법은 전천구 스테레오스코픽 분석보다 특정 소규모 영역에 대한 감도를 크게 향상시킬 수 있다.