Measurement of the Hubble constant using the Dark Energy Survey Year 6 Gold galaxy catalogue and the fourth Gravitational-Wave Transient Catalogue

기술 요약: DES Y6 Gold 및 GWTC-4.0을 이용한 허블 상수 측정

문제 정의
초기 우주(CMB)와 후기 우주(초신성) 측정치 사이의 허블 상수($H_0$) 긴장은 우주론의 중대한 과제로 남아 있다. 중력파(GW) 표준 사이렌은 $H_0$를 측정하기 위한 독립적인 경로를 제공한다. "밝은 사이렌(bright sirens)"(GW170817과 같이 전자기파 대응물이 있는 사건)은 직접적인 적색편이 측정을 제공하지만, 그 빈도가 매우 낮다. 대응물이 없는 "어두운 사이렌(dark sirens)"(event without counterparts)은 은하 카탈로그와 GW 국지화 부피를 통계적으로 상관시켜 적색편이를 추론하는 방식에 의존한다. DES를 사용한 기존 분석은 단일 사건이나 고정된 모집단 모델에 국한되었으며, GLADE+와 같은 기존 카탈로그는 $z > 0.05$에서 불완전성을 겪는다. 본 연구는 우주론적 파라미터와 GW 모집단 파라미터의 공동 추론을 위해 전체 네 번째 중력파 과도 현상 카탈로그(GWTC-4.0)를 활용할 수 있도록 깊고 넓은 시야의 광도 측정 카탈로그를 적응시키는 문제를 다룬다.

방법론
저자들은 gwcosmo 파이프라인을 사용하여 우주론적 파라미터($H_0$)와 컴팩트 이중성 병합(CBC) 모집단 파라미터의 공동 베이지안 추론을 수행한다.

• 데이터 소스:
  • GW 사건: 137개의 이중 블랙홀과 5개의 잠재적 중성자별 성분을 포함한 GWTC-4.0(O4a 런)의 142개 CBC 후보. GW170817은 별도의 밝은 사이렌으로 취급된다.
  • 은하 카탈로그: 4923 deg$^2$를 커버하며 $r$-band 깊이가 23.9 mag인 DES Year 6 (Y6) Gold 광도 측정 카탈로그.
• 카탈로그 준비 및 선택:
  • 객체 선택: 오분류된 밝은 별을 제거하기 위해 품질 컷(FLAGS_FOOTPRINT, FLAGS_FOREGROUND, FLAGS_GOLD)과 별-은하 분리 기법인 XGB_PRED(그래디언트 부스팅 결정 트리)를 적용하였다.
  • 적색편이 범위: $0.05 < z < 0.50$으로 제한하였다. 하한값($z > 0.05$)은 DES의 CCD 포화 효과를 완화하기 위함이며, 상한값($z < 0.5$)은 SDSS 필터에서 유도된 $K$-correction의 유효성을 보장하고 $u$-band 데이터가 누락된 $0.15 < z < 0.35$ 구간에서의 광도 측정 적색편이 편향을 피하기 위함이다.
  • 광도 함수: 카탈로그 완결성을 정량화하기 위해 $1/V_{max}$ 방법을 사용하여 $r$-band Schechter 광도 함수를 측정하였다. 진화 파라미터($M^*(z)$, $\alpha(z)$, $\phi^*(z)$)는 적색편이에 따라 선형적으로 모델링되었다.
  • 가중치 부여: 병합률이 항성 질량에 비례한다고 가정하여, 호스트 은하 확률을 광도($\epsilon=1$)에 따라 가중치를 두었다.
• 시선(LOS) 적색편이 사전 분포(Prior):
  • 카탈로그 내 은하(광도 및 광도 측정 적색편이 불확실성에 의해 가중치가 부여됨)와 미검출 은하를 위한 "빈 카탈로그(empty catalogue)" 항(공변 부피에 대해 균등함)을 결합하여 LOS 사전 분포를 구축하였다.
  • 실제 대규모 구조(예: 은하 벽) 및 광도 측정 적색편이 편향으로 인해 발생하는 최대 ~21%의 편차를 주목하며, 이 사전 분포를 균등 가정에 대해 검증하였다.
• 컴퓨팅 프레임워크:
  • 사후 분포 샘플에 대한 커널 밀도 추정(KDE)을 처리하기 위해 GPU 벡터화(PyTorch 및 Numba/CUDA 활용)를 통한 하드웨어 가속을 구현하였다. 이를 통해 분석 실행 시간을 몇 주에서 몇 시간으로 단축하여 수백 개의 사건을 처리할 수 있게 하였다.

주요 기여

1. DES Y6 Gold의 적응: DES Y6 Gold 카탈로그가 하늘의 약 12%만을 커버함에도 불구하고, 이를 어두운 사이렌 우주론에 성공적으로 적응시켰음을 입증하였으며, 이 카탈로그는 풋프린트 내에서 GLADE+보다 약 100~1000배 더 높은 은하 밀도를 제공한다.
2. 공동 추론: 이전의 DES 분석에서 사용된 고정된 모집단 모델을 넘어, 전체 GWTC-4.0 데이터셋과 깊은 광도 측정 카탈로그를 사용하여 $H_0$와 CBC 모집단 파라미터의 첫 번째 공동 추론을 수행하였다.
3. 편향 완화: 엄격한 적색편이 컷($0.05 < z < 0.5$)을 적용함으로써 DES 적색편이 분포의 특정 편향(저 $z$에서의 포화 및 중간 $z$에서의 광도 측정 적색편이 산란)을 식별하고 완화하였다.
4. 컴퓨팅 효율성: gwcosmo를 통해 GPU 가속으로 1~2 자릿수(order of magnitude)의 속도 향상을 달-성하였으며, 이는 대규모 어두운 사이렌 분석을 가능하게 한다.

결과

• 어두운 사이렌만 사용 시: $0.05 < z < 0.5$ 컷이 적용된 DES Y6 Gold 카탈로그를 사용하여 $H_0 = 71.7^{+22.7}_{-18.3}$ km s$^{-1}$ Mpc$^{-1}$를 측정하였다. 이는 GLADE+를 사용한 이전 O4a 분석에 비해 68.3% 신뢰 구간 폭이 약 11% 감소한 결과이다.
• 결합 (어두운 + 밝은 사이렌): 밝은 사이렌인 GW170817과 어두운 사이렌 결과를 결합하여 $H_0 = 73.1^{+11.9}_{-8.5}$ km s$^{-1}$ Mpc$^{-1}$를 얻었다. 이는 O4a 표준 사이렌 결과보다 정밀도가 약 10% 향상된 것이다.
• 제약 능력: 저자들은 이 데이터셋에서 $H_0$에 대한 제약 능력의 대부분이 은하 카탈로그 자체보다는 "스펙트럼 사이렌(spectral siren)" 측면(질량 분포의 특징)에서 온다는 점에 주목하였다. 은하 카탈로그는 대부분의 GW 사건의 국지화 부피를 덮기에는 너무 좁거나 얕기 때문에, 현재는 "빈 카탈로그" 사례에 비해 미미한 개선만을 제공한다.

의의 및 주장
본 논문은 깊은 은하 탐사가 GW 우주론을 위해 어떻게 적응될 수 있는지를 입증하며, 카탈로그의 불완전성과 적색편이 편향을 처리하는 데 필요한 방법론을 강조한다. 저자들은 기존 카탈로그의 깊이와 하늘 피복 영역이 GW 국지화 부피에 비해 제한적이지만, 입증된 프레임워크가 차세대 탐사(예: LSST, Euclid, DESI)의 잠재력을 극대화하는 데 필수적이라고 주장한다. 이 연구는 정밀 GW 우주론 시대(O4b, O5 등 미래 관측 런 및 더 깊은 카탈로그가 어두운 사이렌을 주요 도구로 만들 시기)를 위한 준비 단계 역할을 한다. 또한, 정밀도 향상의 핵심은 모집단 파라미터의 공동 추론에서 오며, 은하 카탈로그는 점진적이지만 필요한 세부 조정을 제공한다는 점을 강조한다.