[논문 리뷰] Mass-redshift dependency of Supermassive Black Hole Binaries for the Gravitational Wave Background

Musfar Muhamed Kozhikkal, Siyuan Chen외 3명|arXiv (Cornell University)|2023. 05. 29.

Galaxies: Formation, Evolution, Phenomena인용 수 4

한 줄 요약

이 논문은 적색편이에 따라 변하는 블랙홀-버진 질량 관계를 초거대 블랙홀 이원성(SMBHB)에서 기인하는 중력파 배경(GWB)과 연결하는 분석적 천체물리 모델을 개발한다. 은하의 항성 질량 함수, 쌍의 비율, 합쳐짐 시간 상수, 그리고 천체역학 시뮬레이션(Illustris, SIMBA 등)을 조합하여 주파수의 함수로 중력파 배경 강도를 예측하고, 베이지안 추론을 통해 BH-버진 매개변수를 제약한다. 그 결과, 주요 SMBHB 기여는 z ~ 1.5–2.5에서 발생하며, PTA 관측은 BH-버진 질량 관계의 적색편이에 따른 진화를 강력하게 제약함을 밝혀낸다.

ABSTRACT

Studying how the black hole (BH) - (galaxy) bulge mass relation evolves with redshift provides valuable insights into the co-evolution of supermassive black holes and their host galaxies. However, obtaining accurate measurement of BH masses is challenging due to the bias towards the most massive and luminous galaxies. Instead we focus on the BH and bulge masses as they vary with redshift using the EAGLE, Illustris, TNG100, TNG300, Horizon-AGN and SIMBA large-scale cosmological simulations. We use an analytical astrophysical model with galaxy stellar mass function, pair fraction, merger timescale and BH-bulge mass relation extended to include redshift evolution. The model can predict the intensity of the gravitational wave background (GWB) produced by a population of supermassive black hole binary (SMBHB) as a function of the frequency. This allows us to compare the predictions of this model with the constraints of Pulsar Timing Array observations. Here, we employ Bayesian analysis for the parameter inference. We find that all six simulations are consistent $\leq 3.5σ$ with a range of simulated GWB spectra. By fixing the BH-bulge mass parameters to the simulations we analyze the changes in the constraints on the other astrophysical parameters. Furthermore, we also examine the variation in SMBHB merger rate with mass and redshift between these large-scale simulations.

연구 동기 및 목표

• 초거대 블랙홀 이원성(SMBHB)에서 기인하는 확률적 중력파 배경(GWB)에 영향을 미치는 블랙홀-버진 질량 관계의 적색편이 진화를 모델링하는 것.
• 대규모 천체역학 시뮬레이션(Illustris, SIMBA 등)에서 유도된 적색편이에 따라 변화하는 BH-버진 질량 관계를 포함하여 SMBHB 합쳐짐 비율의 천체물리적 모델을 향상시키는 것.
• 최근에 관측된 저주파수 공통 모드 신호를 포함하여, 이론적 GWB 예측을 펄서 타이밍 어레이(PTA) 관측과 비교하는 것.
• 베이지안 추론을 통해 시뮬레이션 데이터에 기반하여 BH-버진 질량 관계의 정규화 및 적색편이 진화를 포함한 핵심 천체물리적 매개변수를 제약하는 것.
• 다양한 시뮬레이션 출력에서의 비선형성, 체르프 질량, 적색편이 분포에 따른 GWB 진폭 예측의 민감도를 평가하는 것.

제안 방법

• 모델은 은하 항성 질량 함수(GSMF), 이원성에서의 차별적 쌍 비율, 적색편이에 따라 변화하는 합쳐짐 시간 상수를 조합하여 SMBHB 합쳐짐 비율을 계산한다.
• 천체역학 시뮬레이션(EAGLE, Illustris, TNG100, TNG300, Horizon-AGN, SIMBA)에서 유도된 γ*와 β*로 매개변수화된 적색편이 진화를 거친 BH-버진 질량 관계를 통합한다.
• SMBHB 합쳐짐 비율과 중력파 방출의 에너지 스펙트럼을 사용하여 주파수의 함수로 중력파 배경 강도를 계산한다.
• 베이지안 추론을 적용하여 SMBHB 인구의 시뮬레이션 데이터에 기반해 BH-버진 질량 관계 매개변수(γ*, β*)를 피팅하고, GWB 진폭 기대치로부터 매개변수 제약을 가능하게 한다.
• 특성 스트레인 스펙트럼(곡선형 및 평탄한 형태)을 예측하고, PTA 관측 제약과 비교하며, 특히 h_c = 2×10⁻¹⁵ 및 4×10⁻¹⁵의 경우를 중심으로 분석한다.
• 체르프 질량과 적색편이에 따른 합쳐짐 비율을 분석하며, 원형 및 비원형 이원성에 대해 별도로 다루어 동역학적 영향을 평가한다.

실험 결과

연구 질문

• RQ1 블랙홀-버진 질량 관계의 적색편이 진화는 초거대 블랙홀 이원성(SMBHB)에서 기인하는 확률적 중력파 배경(GWB)의 예측 진폭에 어떤 영향을 미치는가?
• RQ2 천체역학 시뮬레이션에서 유도된 GWB 예측과 함께, 베이지안 추론을 통해 BH-버진 질량 관계 매개변수(γ*, β*)에 대한 제약는 무엇인가?
• RQ3 SMBHB가 GWB에 가장 크게 기여하는 적색편이는 어느 정도이며, 이는 이원성의 비선형성과 체르프 질량에 따라 어떻게 달라지는가?
• RQ4 다른 천체역학 시뮬레이션(Illustris 대비 SIMBA) 간의 GWB 예측은 합쳐짐 비율과 진폭 스펙트럼 측면에서 어떻게 비교되는가?
• RQ5 이 모델은 관측된 PTA 공통 모드 신호를 대규모 SMBHB 인구와 조화시킬 수 있는가? 이는 z < 5에서 BH-버진 질량 관계에 대해 어떤 함의를 갖는가?

주요 결과

• 중력파 배경(GWB)의 주요 기여는 적색편이 z ~ 1.5–2.5에서 발생하며, 더 높은 적색편이에서는 작은 비율이지만 무시할 수 없는 기여를 한다.
• h_c = 4×10⁻¹⁵인 경우 원형 이원성은 z < 1.0에서 거의 달성되지 않으며, 시뮬레이션 샘플에서 이 정도 수준에 도달하는 비원형 이원성 역시 극히 소수에 불과하다.
• h_c = 2×10⁻¹⁵일 경우, Illustris 시뮬레이션에서 적색편이에 관계없이 합쳐짐 비율이 거의 일정해지며, SMBHB 인구의 포화 또는 안정된 기여를 나타낸다.
• 모델은 β* (BH-버진 질량 관계의 정규화)에 대해 γ* (적색편이 진화 지수)보다 더 강한 제약를 제공함을 발견하여, PTA 관측이 질량 관계의 정규화에 가장 민감함을 시사한다.
• 제안된 적색편이에 따라 변화하는 BH-버진 질량 관계는 이전의 관계(Venemans et al. 2016 등)보다 시뮬레이션(Illustris, SIMBA 등) 간 산란이 훨씬 작음을 보여주며, 시뮬레이션 데이터와의 일관성을 향상시킨다.
• 특히 2-시그마 하한부에서 원형 및 비원형 이원성 간 체르프 질량 분포에 미세하지만 측정 가능한 차이가 있음을 보여주며, 이는 비선형성이 GWB 스펙트럼 형상에 영향을 미친다는 것을 시사한다.