The Stochastic Siren: Astrophysical Gravitational-Wave Background Measurements of the Hubble Constant

기술 요약: 스토캐스틱 사이렌(The Stochastic Siren): 허블 상수의 천체물리학적 중력파 배경 측정

문제 제기
본 논문은 초기 우주 측정치($H_0 \approx 67\text{--}68 \text{ km s}^{-1} \text{ Mpc}^{-1}$)와 후기 우주 측정치($H_0 \approx 72\text{--}74 \text{ km s}^{-1} \text{ Mpc}^{-1}$) 사이의 유의미한 불일치($>5\sigma$)인 "허블 긴장(Hubble tension)"을 다룬다. 중력파(GW) "표준 사이렌"(개별적으로 분해된 병합 사건)은 $H_0$를 측정하는 독립적인 방법을 제공하지만, 스펙트럼 사이렌 방법(질량 스펙트럼 특징 활용)을 사용하는 현재의 LIGO-Virgo-KAGRA(LVK) 협업체의 제약은 여전히 광범위하다($H_0 = 46^{+49}_{-26} \text{ km s}^{-1} \text{ Mpc}^{-1}$). 나아가, 미분해된 항성 질량 블랙홀 쌍성(BBH) 병합으로부터 발생하는 스토캐스틱 중력파 배경(GWB)은 아직 검출되지 않았다. 저자들은 이 배경의 '비검출(non-detection)'이 여전히 우주의 팽창 역사에 대한 의존성을 활용하여 $H_0$에 대한 우주론적 제약을 제공할 수 있는지 조사한다.

방법론
저자들은 분해된 전경(foreground, FG) BBH 병합의 관측 가능성(likelihood)과 미분해 배경(background, BG) GWB 신호의 가능성을 결합한 "스토캐스틱 사이렌"이라 명명된 새로운 방법을 제안한다.

1. GWB의 우주론적 의존성: 무차원 중력파 에너지 밀도 $\Omega_{\text{gw}}$는 $H_0$에 반비례하여 유도된다. 구체적으로, (고정된 병합율 밀도를 가정할 때) $\Omega_{\text{gw}} \propto H_0^{-3}$이다. 낮은 $H_0$는 더 큰 공변 부피(comoving volumes)를 의미하며, 이는 더 높은 통합 병합 횟수와 그에 따른 더 강한 $\Omega_{\text{gw}}$로 이어진다. 반대로, 높은 $H_0$는 더 약한 신호를 초래한다. 따라서, GWB의 비검출(즉, $\Omega_{\text{gw}}$에 대한 상한선)은 검출 가능한 배경을 예측하는 낮은 $H_0$ 값을 배제할 수 있다.
2. 가능도(Likelihood) 구성:
   • 전경($L_{\text{FG}}$): GWTC-3에서 관측된 42개의 분해된 BBH 카탈로그의 가능도는 인구 통계적 하이퍼파라미터($\vec{\lambda}$) 및 $H_0$에 의존하는 불균질 포아송 분포(inhomogeneous Poisson distribution)를 사용하여 모델링된다.
   • 배경($L_{\text{BG}}$): 교차 상관 통계량 $\hat{C}$(GWB 탐색에 사용됨)의 가능도는 가우시안 분포로 모델링된다. 이 가능도는 예측된 $\Omega_{\text{gw}}(\vec{\lambda}, H_0)$에 의존한다. 저자들은 $\Omega_{\text{gw}}$에 대한 균등 사전 확률(uniform prior)을 가정하면서, 기존의 LVK 탐색으로부터 얻은 $\Omega_{\text{gw}}$ 상한선을 사용하여 이 가능도를 구축한다.
   • 결합 가능도(Joint Likelihood): 전경 신호와 배경 신호 사이의 독립성을 가정할 때(분해된 병합이 전체 $\Omega_{\text{gw}}$ 상한선에 기여하는 바가 매우 미미하므로($\sim0.3\%$), 이는 정당화된다), 결합 가능도는 $L_{\text{joint}} = L_{\text{FG}} L_{\text{BG}}$가 된다.
3. 인구 모델: 분석에는 질량 분포를 모델링하기 위한 8개의 파라미터를 가진 PowerLaw+Peak 모델과 적색편이 진화를 위한 Madau-Dickinson 함수 형태를 가진 4개의 파라미터를 포함하여 총 12개의 하이퍼파라미터를 사용하는 BBH 인구 모델이 채택된다.
4. 추론: 저자들은 15개의 파라미터(12개 인구 파라미터 + 3개 우주론적 파라미터: $H_0, \Omega_{m,0}, w_0$)에 대해 중첩 샘플링(nested sampling, bilatby 내의 dynesty 샘플러 사용)을 통해 베이지안 추론을 수행한다. 전경 가능도를 효율적으로 평가하기 위해 이전 GWTC-3 인구 분석에서 도출된 커널 밀도 추정기(KDE)를 활용한다.

주요 결과

• 정밀도 향상: 결합된 "스토캐스틱 사이렌" 측정값은 스펙트럼 사이렌 방법 단독 사용 시보다 더 타이트한 $H_0$ 제약을 제공한다. 결합 측정의 최대 사후 확률(MAP)은 $H_0 = 72^{+44}_{-37} \text{ km s}^{-1} \text{ Mpc}^{-1}$로, 전경 전용 분석의 $57^{+43}_{-35} \text{ km s}^{-1} \text{ Mpc}^{-1}$ 및 스펙트럼 사이렌의 $46^{+49}_{-26} \text{ km s}^{-1} \text{ Mpc}^{-1}$와 비교된다.
• 사후 분포의 이동: GWB 비검출을 포함하면 사후 분포가 더 높은 $H_0$ 값 쪽으로 이동하며, 이는 후기 우주의 전자기학적 측정치(SH0ES)와 더 잘 일치하게 되면서도, 큰 불확실도 내에서 초기 우주 측정치(Planck)와는 여전히 일관성을 유지한다.
• 하한선 메커니즘: 저자들은 GWB의 비검출이 효과적으로 $H_0$에 대한 하한선을 설정함을 입증한다. 관측 시간이 증가하고 GWB가 계속 검출되지 않으면, 이 하한선은 점진적으로 상승하여 측정의 정밀도를 높일 것이다.
• 파라미터 제약: GWB 비검출을 포함하는 것이 $H_0$ 제약을 개선하기는 하지만, 현재 데이터의 민감도를 고려할 때 다른 인구 하이퍼파라미터(질량 분포)나 우주론적 파라미터($\Omega_{m,0}, w_0$)에 대한 제약을 유의미하게 강화하지는 못한다.

의의 및 주장
본 논문은 스토캐스틱 중력파 배경을 이용한 최초의 $H_0$ 측정을 제시한다고 주장한다. 본 연구의 주요 의의는 검출된 배경의 '부재'가 강력한 우주론적 탐사 도구가 될 수 있음을 보여준 데 있다. "스토캐스틱 사이렌"은 다음과 같은 이유로 허블 긴장을 해결하는 데 독특한 역할을 한다:

1. 전자기학적 관측 및 우주 거리 사다리와 완전히 독립적이다.
2. 배경이 직접 검출되기 전에도, 비검출이 지속됨에 따라 $H_0$를 아래에서부터 점진적으로 제약할 수 있는 메커니즘을 제공한다.
3. 기존의 "표준 사이렌" 및 "스펙트럼 사이렌" 접근법을 보완한다.

저자들은 실제 배경이 다른 소스(예: 중성자별 쌍성)를 포함하여 $\Omega_{\text{gw}}$를 증가시키고 $H_0$ 제약을 더욱 강화할 수 있음에도 불구하고, 오직 BBH 기여만을 고려했기에 본 결과가 보수적이라고 언급한다. 향-미래의 전망에 따르면, 검출기 업그레이드(예: A#)와 함께 GWB가 1년 미만의 시간 내에 신호 대 잡음비(SNR) 8로 검출될 수 있으며, 이 시점에서 스토캐스틱 사이렌 방법은 훨씬 더 강력한 우주론적 제약을 제공할 것이다.