라이다와 LGWA를 위한 중간·초대질량 블랙홀 합병률 예측

초록

본 연구는 라이다(LISA)와 달 기반 중력파 탐지기(LGWA)의 감도 범위에 맞춰, 가벼운 시드(빛나는 별 붕괴) 모델을 이용해 중간질량(10³–10⁵ M☉) 및 초대질량(>10⁵ M☉) 블랙홀(BH) 이진의 합병률을 반감기와 최악·최선 시나리오를 통해 추정한다. PINOCCHIO 시뮬레이션과 GAEA 기반 반정밀 모델을 결합해 은하·다크물질 병합 트리를 생성하고, Chandrasekhar 동역학 마찰 공식을 이용해 BH 병합까지의 시간을 계산한다. 낙관적 가정에서는 연간 15건 이상의 IMBH 합병이 감지될 수 있으며, 비관적 경우에도 약 5건 이하가 검출될 것으로 예측한다. 고‑z(>8)에서도 일부 합병이 탐지 가능하고, LGWA는 고‑SNR IMBH, LISA는 무거운 SMBH에 더 민감함을 확인한다.

상세 분석

이 논문은 차세대 중력파 관측기인 LISA와 달 기반 LGWA가 탐지할 수 있는 블랙홀 합병 사건의 발생률을, ‘라이트 시드’(light‑seed) 모델을 기반으로 정량화한 최초의 연구 중 하나이다. 저자들은 먼저 PINOCCHIO 코드를 이용해 59.7 Mpc ³ 부피의 다크물질(N‑body) 시뮬레이션을 수행하고, 이를 통해 은하와 다크물질의 병합 트리를 고해상도로 재구성한다. 이어서 GAEA 반정밀 모델(CAM25)을 적용해 각 다크물질 halo에 별 형성, 가스 냉각, 피드백, 그리고 블랙홀 성장 과정을 삽입한다. 핵심은 ‘라이트 시드’ 가정으로, 모든 은하 중심에 10 M☉~10⁵ M☉ 수준의 초기 블랙홀을 배치한다는 점이다. 이는 Pop III 별 붕괴 혹은 조밀성단 내 런어웨이 합병을 반영한 것으로, 초기 질량 분포가 넓어 IMBH와 SMBH가 동시에 형성될 수 있다.

병합 시간 계산은 Chandrasekhar 동역학 마찰 공식에 기반한다. 저자들은 두 가지 경계값을 설정했는데, ‘낙관적’ 경우는 병합 대상 블랙홀 질량을 최대값(즉, 병합 전후 질량 합계)으로 잡아 마찰이 빠르게 작용하도록 하고, ‘비관적’ 경우는 질량을 최소값(주로 작은 위성 블랙홀)으로 가정해 마찰 시간을 길게 만든다. 이를 통해 동역학 마찰 시간의 불확실성을 상한·하한으로 포괄한다.

감도 계산은 LISA와 LGWA 각각의 파워 스펙트럼 밀도와 목표 관측 기간(5년, 10년 등)을 적용해 신호‑대‑잡음비(SNR) 임계값을 8로 설정하였다. 결과적으로, 낙관적 시나리오에서는 연간 15건 이상의 IMBH(10³–10⁵ M☉) 합병이 두 탐지기 모두에서 SNR > 8을 만족한다. 비관적 시나리오에서는 연간 5건 이하가 검출 가능하며, 대부분은 z ≈ 2–4 범위에 집중된다. 고‑z(>8) 합병도 일부 감지될 수 있지만, 사건률은 급격히 감소한다. 질량별로 보면, LISA는 10⁵–10⁷ M☉ SMBH 합병에 최적이며, LGWA는 10³–10⁵ M☉ IMBH 합병에 대해 높은 SNR을 제공한다. 두 탐지기의 파라미터 공간이 겹치는 구간에서는 공동 관측이 가능해, 파라미터 추정 정확도가 크게 향상될 것으로 기대된다.

논문은 또한 기존 문헌과 비교해, 라이트 시드 기반 SAM이 예측하는 합병률이 ‘헤비 시드’ 모델(10⁴–10⁶ M☉ 초기 질량)보다 1–2 dex 높으며, 이는 IMBH 존재 가능성을 크게 확대한다는 점을 강조한다. 그러나 저자들은 모델이 초기 별 형성, 가스 흐름, 그리고 ‘떠도는 블랙홀’(wandering BH) 등을 단순화했으며, 실제 우주에서는 추가적인 물리적 과정(예: 삼중계 상호작용, 가스 디스크 마찰 등)이 합병률을 변동시킬 수 있음을 인정한다.