A Joint Optimal Search for Gravitational Waves from Resolved and Unresolved Supermassive Binary Black Holes with Pulsar Timing Arrays

기술 요약: 분해된 및 분해되지 않은 초거대 질량 블랙홀 쌍성계로부터의 중력파에 대한 공동 최적 탐색

문제 정의
현재의 펄서 타이밍 어레이(PTA) 분석은 초거대 질량 블랙홀 쌍성계(SMBHB)로부터 발생하는 중력파 배경(GWB)과 개별적으로 분해 가능한 연속파(CW) 신호를 서로 분리되고 단절된 문제로 취급한다. GWB는 일반적으로 확률론적이고 등방적인 가우시안 과정으로 모델링되는 반면, 개별 소스 탐색에는 변형 진폭(strain amplitude)에 대해 무정보 사전 분포(uninformative priors)를 사용한다. 이러한 파편화는 GWB와 가장 밝은 개별 소스 모두를 생성하는 천체물리학적 집단 특성(수밀도 및 질량 규모 등)에 대한 통합된 추론을 방해한다. 또한, 표준 모델들은 종종 가우시안 대수 통계에서 포아송 통계로의 전이를 고려하지 못하는데, 이는 SMBHB 소스의 이산적인 본질에 있어 매우 중요하다.

방법론
저자들은 분해되지 않은 GWB와 가장 밝은 분해 가능한 SMBHB 소스를 동시에 설명하는, 제1원리로부터 유도된 합동 계층 모델(joint hierarchical model)을 도입한다.

• 이론적 프레임워크: Sato-Polito와 Zaldarriaga [15]의 연구를 바탕으로, 이 모델은 특성 소스 수 $N_c$와 특성 진폭 $h_c$를 천체물리학적 하이퍼파라미터로 활용한다. 모델은 전체 특성 진폭($h_t$)의 확률 밀도 함수(PDF)를 가장 밝은 소스($h_{cw}$)와 나머지 집단($h_{t-1}$)의 중첩으로 정의한다.
• 통계적 접근법: 가장 밝은 소스는 결정론적인 방향성 신호(CW)로 취급되는 반면, 나머지 소스들은 등방적인 확률론적 배경을 형성한다. 모델은 중력파 간섭 및 등방성 편차를 마진화(marginalize)하기 위해 가장 밝은 소스의 PDF를 포함한다.
• 데이터 분석: 이 방법은 데이터로부터 추론된 GWB 파워 스펙트럼을 다시 피팅하는 기존 방식과 달리, 모든 신호 및 노이즈 파라미터에 대해 PTA 가능도(likelihood)의 전역 피팅(global fit)을 수행한다. 저자들은 단절된 사전 분포 대신 계층적 모델 $\pi(h_t, h_{cw} | N_c, h_c)$를 사용하여 정보가 배경 추론과 개별 연속파 추론 사이에서 흐를 수 있도록 한다.
• 검증: 방법론은 알려진 노이즈, 관측치 및 추론된 GWB 파와 스펙트럼을 재현한 NANOGrav 15년 시뮬레이션 데이터를 통해 테스트된다. 시뮬레이션에는 15.06년 동안 100 ns의 타이밍 정밀도를 가진 20개의 펄서가 포함된다.

주요 기여

1. 통합 모델: 본 논문은 단일 가능도 프레임워크 내에서 완전히 분해된 소스와 완전히 분해되지 않은 SMBHB 소스를 공동으로 모델링하는 최초의 확정적이고 검증 가능한 모델을 제시한다.
2. 검출 통계량: 저자들은 GWB의 천체물리학적 기원을 위한 주요 검출 통계량으로 특성 소스 수 $N_c$를 제안한다. $N_c \ll 1000$인 경우 포아송/비가우시안 영역을 나타내며, $N_c \geq 1000$인 경우 표준 가우시안 근사에 해당한다.
3. 개선된 제약 조건: 무정보 사전 분포 대신 천체물리학적으로 동기화된 계층적 사전 분포를 사용함으로써, 개별적으로 분해 가능한 SMBHB의 진폭에 대해 더 엄격한 제약을 제공한다.
4. 집단 추론: 이 프레임워크를 통해 GWB 파라미터($N_c, h_c$)를 물리적 집단 특성인 SMBHB 질량 규모($M_{peak}$) 및 풍부도($\rho_{BH}$)로 직접 재구성할 수 있다.

결과
시뮬레이션된 NANOGrav 15년 데이터에 적용한 결과, 다음과 같은 결과가 도출되었다:

• GWB의 본질: 시뮬레이션된 데이터에서 가우시안 GWB($N_c \approx 1000$)라는 귀무가설을 기각할 만한 충분한 증거가 없다. 데이터는 SMBHB 특유의 포아송 영역에 대해 민감하지 않다.
• 개별 소스: 개별적으로 분해 가능한 CW 신호는 검출되지 않았다. 가장 밝은 소스의 주파수($f_{cw}$)에 대한 사후 분포는 전체 사전 범위에 걸쳐 있으며, 사후 예측 샘플은 특별히 밝은 이상치를 요구하지 않고도 데이터와 일치한다.
• AGN 후보와의 긴장 관계: 새로운, 더 엄격한 CW 진폭 상한선은 활동 은하 핵(AGN) 관측으로부터 얻은 114개의 SMBHB 후보 중 21개를 NANOGrav 데이터와 모순(tension)되게 만든다. 반면, 기존 분석(무정보 사전 분포 사용)에서는 단 하나의 후보만이 모순되는 것으로 나타났다.
• 검출 확률:
  • 15년 데이터에서 SNR $\geq 5$인 CW를 검출할 확률은 **2%**로 추정된다.
  • 예상되는 20년 데이터의 경우, 이 확률은 **5%**로 증가한다.
  • 그러나 20년 데이터에서 낮은 유의성을 가진 이상치(SNR $\geq 2$)를 발견할 확률은 **40%**로 추정된다.

의의 및 주장
본 논문은 이 공동 모델이 GWB의 천체물리학적 성격을 식별하기 위한 "결정적 증거(smoking gun)" 접근법이라고 주장한다. 순수하게 가우시안 배경이라는 가정을 넘어, 이 모델은 신호의 이산적 SMBHB 기원을 확인할 수 있는 소스 수 변동(비가우시안성)의 검출을 가능하게 한다.

저자들은 이들의 방법이 파라미터 추정의 퇴화(degeneracy)를 피하기 위해 전역 피팅을 수행하므로, 기존의 인수분해된 가능도(factorized likelihood) 접근 방식보다 더 완전하고 정확하다고 강조한다. 현재의 시뮬레이션 데이터는 분해 가능한 CW를 보여주지 않지만, 이 프레임워크는 사후 분포와 사후 예측 샘플 간의 일관성, 가장 밝은 소스의 주파수에 대한 넓은 지지, 검출 임계값을 초과할 낮은 확률 등을 포함하는 미래의 검출을 위한 엄격한 체크리스트를 구축한다.

궁극적으로, 본 논문은 미래의 PTA 데이터가 단순히 GWB 진폭을 측정하는 것을 넘어, $(\rho_{BH}, M_{peak})$ 평면에서 천체물리학적으로 기대되는 영역과 더 작은 $N_c$ 값 쪽으로 사후 분포를 좁힘으로써, 유한한 수의 초거대 질량 블랙홀 집단으로부터 기대되는 가우시안에서 포아송 통계로의 전이를 확인하고 SMBHB 해석을 검증할 것이라고 주장한다.