[논문 리뷰] Beyond Gaps and Bumps: Spectral Siren Cosmology with Non-Parametric Population Models

I. Magaña Hernandez, Anarya Ray|arXiv (Cornell University)|2024. 04. 03.

Cosmology and Gravitation Theories인용 수 6

한 줄 요약

본 논문은 비모수적, 구간화된 가우시안 프로세스(population model)를 사용하여 GWTC-3 BBH 탐지를 바탕으로 스펙트럴 사이렌 우주론을 수행하고, GW 데이터만으로 BBH 질량 스펙트럼과 허블 상수(Hubble constant)를 공동으로 제약합니다. GW170817과 결합 시 GW170817만으로 얻은 제약보다 더 촘촘한 H0 제약을 얻습니다.

ABSTRACT

Gravitational wave standard sirens typically require electromagnetic (EM) data to obtain redshift information to constrain cosmology. Difficult to find EM counterparts for bright sirens and galaxy survey systematics for dark sirens make cosmological constraints with spectral sirens, a gravitational wave data-only approach, extremely appealing. In this work, we use the GWTC-3 BBH detections as spectral sirens to constrain the BBH population and the underlying cosmological expansion with a flexible model for the black hole mass spectrum. We use a binned Gaussian process to model the BBH mass distribution in the source frame without any astrophysical assumptions on the shape and or inclusion (or lack of) features that drive the cosmological constraints as the redshifted detector frame masses become consistent with the underlying astrophysical mass distribution features. For GWTC-3 we find a measurement on the Hubble constant of $H_0=73.0^{+13.3}_{-7.7} \ { m{km \ s^{-1} \ Mpc^{-1}}}$ at $68\%$ C.L. when combined with that obtained from the bright standard siren analysis with GW170817 and its associated host galaxy NGC 4993. We find an improved estimate for the Hubble constant of around a factor of 1.4 times better than the GW170817 measurement alone. We validate our nonparametric spectral siren approach with simulations and benchmark its scalability and constraining performance when compared with parametric methods.

연구 동기 및 목표

• BBH 질량 특징에 대한 천체물리학 가정을 최소화하기 위해 비모수적이고 데이터 기반의 인구 모델을 동력으로 삼아 스펙트럼 사이렌 우주론을 모티브로 삼는다.
• BBH 질량 스펙트럼, 합병률 진화, 그리고 우주상수를 공동으로 추정하는 계층적 베이지안 프레임워크를 개발한다.
• 정해진 형태 가정 없이 소스 프레임에서 BBH 질량 분포를 모델링하기 위해 구간화된 가우시안 프로세스 사전(prior)을 구현한다.
• 감지 프레임에서의 관측치와 선택 효과를 포함하여 GW 데이터만으로 우주론을 추론하는 방법(스펙트럴 사이렌)을 도입한다.
• 시뮬레이션으로 방법을 검증하고 파라메트릭 인구 모델과의 성능을 비교한다.

제안 방법

• CBC 인구를 탐지 프레임 질량과 적색편차에서의 구간화된 비상관 함수로 모델링한다(식(1)).
• 탐지 프레임 관측치를 우주론 의존적 적색편차를 사용하여 소스 프레임 물리量로 변환하고, 우주론과 인구를 공동 추론하는 것을 가능하게 한다(식(2)).
• 탐지를 균일하지 않은 포아송 과정으로 취급하고 탐지 가능 구간 부피와 재가중된 사건 사후확률을 포함하는 우도(likelihood)를 쓴다(식(3)).
• 권선향 편향(Malmquist bias) 보정을 포함하여 detector 잡음에 대한 주입(injections)을 통해 예상 탐지 수를 추정한다(식(4)).
• GP 평균, 길이 스케일, 진폭에 대한 하이퍼파라미터를 가진 가우시안 프로세스 사전(prior)을 설정하고 해밀토니언 몬테카를로(HMC)로 샘플링한다(식(5-6)).
• GWTC-3 BBH 이벤트(69건)와 GW170817를 사용하여 공동 H0 제약을 얻고, 파라메트릭 접근법과의 벤치마크를 수행한다.

Figure 1: Inferred constraints for both the primary and secondary mass distributions using the 69 BBH mergers from GWTC-3 and our BGP spectral siren model. We show our results as the binned piecewise function plotted in purple and for comparison, we also show the corresponding results obtained by th

실험 결과

연구 질문

• RQ1 비-parametric하고 구간화된 가우시안 프로세스 인구 모델이 특정 스펙트럴 특징을 가정하지 않고도 BBH 질량 분포를 복원할 수 있는가?

주요 결과

• GWTC-3 BBH 탐지에서 GW 데이터만으로 BGP 모델을 사용할 때 H0 = 82.2^{+51.7}_{-33.7} km s^{-1} Mpc^{-1} (68% C.L.)로 측정된다.
• GW170817과의 결합 분석은 H0 = 73.0^{+13.3}_{-7.7} km s^{-1} Mpc^{-1} (68% C.L.)로 제약이 약 1.5배 향상된다( GW170817만 사용했을 때에 비해).
• 추정된 Omega_m은 0.30^{+0.10}_{-0.10}이며, GP 길이 스케일은 l_m = 2.1^{+1.0}_{-1.1}, l_z = 2.1^{+1.0}_{-1.1} (68% C.L.)이다.
• 시뮬레이션 결과 BGP 스펙트럼 사이렌은 fiducial 질량 분포와 우주론을 재현할 수 있으며, 고질량 영역에서 희소한 이벤트로 인한 약간의 불일치가 나타난다.
• 이 방법은 fiducial Powerlaw+Peak 질량 모델, 질량 비율, 적색편차 진화하는 합병률을 특징으로 하는 시뮬레이션에서 편향 없는 추론을 보여주며, 불확실성 내에서의 추론이 가능하다는 것을 시사한다.
• 본 연구는 비파라메트릭 접근법이 더 큰 카탈로그로 확장될 수 있으며, 미래 연구에서 질량-적색편차의 비상관 가정 해제가 이점이 될 수 있음을 시사한다.

Figure 2: Corner plot showing the joint marginalized posterior distribution on the Hubble constant $H_{0}$ , the matter density parameter $\Omega_{m}$ , and the inferred length scales $l_{m},l_{z}$ from our BGP spectral siren model and the 69 GWTC-3 BBH events.