[논문 리뷰] Parameter estimation with the current generation of phenomenological waveform models applied to the black hole mergers of GWTC-1

M. Mateu-Lucena, S. Husa외 6명|arXiv (Cornell University)|2021. 05. 12.

Pulsars and Gravitational Waves Research인용 수 5

한 줄 요약

이 연구는 최신 세대의 현상학적 웨이브폼 모델—주파수 도메인의 IMRPhenomX와 시간 도메인의 IMRPhenomT—를 사용하여 GWTC-1 카탈로그에 포함된 10개의 이중 블랙홀 병합을 재분석하여, 보조 모드와 프리세션을 포함함으로써 매개변수 추정 정확도를 향상시킨다. 이는 원래의 GWTC-1 결과와 일관됨을 확인한다 (Jensen-Shannon 분산 <0.045 비트) 동시에, 특히 프리세션과 고차 모드에서의 웨이브폼 체계적 오차가, GW170729와 GW190521와 같은 고질량 시스템에서 사후 분포에 영향을 줄 수 있음을 보여준다.

ABSTRACT

We consider the ten confidently detected gravitational-wave signals in the GWTC-1 catalog which are consistent with mergers of binary black hole systems, and perform a thorough parameter estimation re-analysis. This is made possible by using computationally efficient waveform models of the current (fourth) generation of the IMRPhenom family of phenomenological waveform models, which consists of the IMRPhenomX frequency-domain modelsand the IMRPhenomT time-domain models. The analysis is performed with both precessing and non-precessing waveform models with and without subdominant spherical harmonic modes. Results for all events are validated with convergence tests, discussing in particular the events GW170729 and GW151226. For the latter and the other two lowest-mass events, we also compare results between two independent sampling codes, Bilby and LALInference. We find overall consistent results with the original GWTC-1 results, with all Jensen-Shannon divergences between the previous results using IMRPhenomPv2 and our default IMRPhenomXPHM posteriors below 0.045 bits, but we also discuss cases where including subdominant harmonics and/or precession influences the posteriors.

연구 동기 및 목표

• 최신 현상학적 웨이브폼 모델을 사용하여 GWTC-1에 포함된 10개의 확실하게 감지된 이중 블랙홀 사건을 재분석함으로써 물리적 정확도를 향상시킨다.
• 매개변수 추정에서 보조 구면 조화 모드와 스핀 프리세션을 포함함으로써 사후 분포에 미치는 영향를 평가한다.
• 수렴 테스트, 독립적인 샘플링 코드(Bilby 및 LALInference), 체계적 사전 분포 변형을 통해 결과를 검증한다.
• 다양한 질량 영역에서 주파수 도메인(IMRPhenomX)과 시간 도메인(IMRPhenomT) 모델 간의 계산적 및 물리적 상호 교환 조건을 평가한다.
• 현재 웨이브폼 모델의 스트레스 테스트를 수행하고 향후 중력파 천체물리학의 기준을 설정한다.

제안 방법

• 수치相对론에 캘리브레이션된 완전한 프리세션과 보조 모드를 포함한 네 번째 세대의 IMRPhenomXPHM 및 IMRPhenomTPHM 웨이브폼 모델을 사용한다.
• Bilby 및 LALInference 코드의 내장된 네스티드 샘플링 알고리즘을 사용하여 베이지안 매개변수 추정을 수행함으로써 방법론적 일관성을 확보한다.
• 넓은 사전 분포와 수렴 테스트를 통해 다양한 샘플링 구성과 사전 가정 하에서 사후 분포의 강건성을 검증한다.
• 모델에 의존하는 체계적 오차를 평가하기 위해 IMRPhenomXPHM(주파수 도메인)과 IMRPhenomTPHM(시간 도메인) 모델 간의 결과를 비교한다.
• 필요에 따라 공개된 GWOSC 데이터를 사용하고, 추가로 자체 제작한 PSD 및 캘리브레이션 엔velope를 보완한다.
• Jensen-Shannon 분산을 사용하여 사후 분포 간의 차이를 정량화함으로써 원래의 GWTC-1 결과와의 일관성을 평가한다.

실험 결과

연구 질문

• RQ1 보조 구면 조화 모드와 스핀 프리세션의 포함이 GWTC-1 사례의 이중 블랙홀 매개변수 사후 분포에 어떤 영향을 미치는가?
• RQ2 특히 위상 및 프리세션에서의 웨이브폼 모델 체계적 오차가 고질량 또는 고스핀 시스템의 매개변수 추정에 얼마나 큰 영향을 미치는가?
• RQ3 현재 세대의 현상학적 모델에서 도출된 사후 결과는 원래의 GWTC-1 분석과 일관성이 있는가, 특히 GW170729와 GW151226 같은 사례에서?
• RQ4 Bilby와 LALInference와 같은 다른 샘플링 코드와 사전 선택이 추론된 사후 분포의 강건성에 어떤 영향을 미치는가?
• RQ5 GWTC-1의 다양한 질량 영역에서 IMRPhenomXPHM 또는 IMRPhenomTPHM 중 어느 웨이브폼 모델이 더 신뢰할 수 있는 결과를 제공하는가?

주요 결과

• GWTC-1 이중 블랙홀 사례 10건 모두에서 본 연구 결과와 원래의 GWTC-1 분석 결과 사이에 일관성이 있으며, Jensen-Shannon 분산이 0.045 비트 이하이다.
• 총 질량가 120 M☉ 이상인 GW170729의 경우, IMRPhenomTPHM 모델이 IMRPhenomXPHM 모델보다 약간 더 우수한 지지를 보이며 프리세션 처리의 우수성 덕분이다.
• GW190521의 경우, IMRPhenomTPHM 모델이 IMRPhenomXPHM 모델보다 유의미하게 더 좋은 결과를 도출하여, 고질량·프리세싱 시스템에 대해 시간 도메인 모델링의 중요성을 강조한다.
• 보조 모드와 프리세션의 포함은 고질량 또는 고스핀 영역에서 사후 분포에 체계적으로 영향을 미치며, 이는 무시할 수 없는 웨이브폼 체계적 오차를 시사한다.
• 수렴 테스트 및 Bilby와 LALInference 간의 비교를 통해 결과의 강건성이 확인되었으며, 극단적인 사전 선택 조건에서도 다중 모달 사후 분포의 증거는 없었다.
• 질량 및 거리 사전 분포의 변화는 사후 분포에 거의 영향을 주지 않아, 합리적인 사전 분포 변형 하에서도 결과가 안정적임을 시사한다.