별질량 블랙홀의 집단역학과 중력파 검출 전망

초록

이 논문은 별집단(플럼머 구상) 내에서 별질량 블랙홀(BH)이 급속히 핵으로 침강해 밀집한 BH‑코어를 형성하고, 3체 상호작용으로 BH‑BH 이중성을 만들며, 이들 이중성이 초탄성 충돌을 통해 경화돼 결국 중력파 방출로 합병하거나 클러스터 밖으로 탈출한다는 과정을 N‑body 시뮬레이션(NBODY6)으로 조사한다. 5 × 10⁴ 이상의 초기 별수를 가진 중간 연령의 대형 클러스터가 현재 관측 가능한 BH‑BH 합병을 가장 많이 생산하며, 예상 검출률은 Advanced LIGO 기준 연간 약 30건이다. 또한 금속량이 낮을수록 BH 질량이 커져 코어 가열이 강해져 BH‑BH 합병률이 상승한다는 금속성 의존성도 제시한다.

상세 분석

본 연구는 별질량 블랙홀(BH)의 동역학적 진화를 직접 N‑body 시뮬레이션(NBODY6)으로 정밀하게 탐구한다. 초기 조건은 플럼머 구형 프로파일을 갖는 반지름 r_h ≤ 1 pc, 별 수 N₀ ≤ 10⁵, 질량 0.5–1 M☉의 저질량 별들 그리고 동일 질량(10 M☉)의 BH를 포함한다. BH 보유 비율은 전부 보존하거나 절반만 보존하는 두 경우를 고려하였다. 시뮬레이션 결과, BH는 질량 구분 불안정(Spitzer instability)으로 약 50 Myr 이내에 핵으로 급속 침강해 밀도 10⁴ M☉ pc⁻³ 수준의 BH‑코어를 형성한다. 이 코어 내에서 3체 상호작용에 의해 BH‑BH 이중성이 지속적으로 생성되며, Heggie 법칙에 따라 초탄성 충돌을 반복하면서 점차 경화된다. 경화된 이중성은 두 가지 운명을 갖는다. 첫째, 경화가 충분히 진행되어 중력파 방출에 의한 합병 시간 T_merge이 수십 Myr 이하가 되면, 클러스터 내부에서 직접 합병한다. 둘째, 충돌에 의해 얻은 반동 속도가 탈출 속도(v_esc)보다 커지면 이중성이 클러스터를 탈출하고, 외부에서 교란이 거의 없으므로 Hubble 시간 이내에 합병할 가능성이 높다.

시뮬레이션에서 N₀ ≥ 5 × 10⁴인 경우, 평균적으로 클러스터당 0.4건(내부 합병)과 0.9건(탈출 후 합병)의 BH‑BH 합병이 발생한다. 이는 클러스터당 연간 1.3건의 합병률로 환산된다. 은하 내 중간 연령 대형 클러스터(초기 질량 M_cl ≈ 3 × 10⁴ M☉ 이상, 연령 ≈ 수 Gyr)의 공간 밀도 ρ_cl ≈ 3.5 Mpc⁻³을 적용하면, Advanced LIGO가 감지할 수 있는 합병 사건은 연간 약 30건에 달한다. 이는 기존의 원시 이중성(프리미털 바이너리) 모델이 예측하는 검출률보다 10배 이상 높은 수치이다.

또한 금속성(Z)의 영향을 추가 실험하였다. Z = 0.02와 Z = 0.001(=0.01 Z☉) 두 경우를 동일 초기 조건으로 진행했을 때, 저금속 클러스터는 BH 평균 질량이 약 14 M☉(고금속은 ≈ 7 M☉)로 두 배 이상 커졌다. 무거운 BH는 코어 가열 효과를 증폭시켜 코어 반경 r_c가 더 크게 팽창하고, BH‑코어의 밀도가 상승해 3체 이중성 형성 및 경화 속도가 가속된다. 결과적으로 저금속 클러스터는 탈출 BH‑BH 이중성 중 Hubble 시간 이내에 합병 가능한 경우가 3건, 고금속은 1건에 불과했다. 이는 금속성이 BH‑코어 동역학과 궁극적인 GW 발생률에 중요한 파라미터임을 시사한다.

본 연구는 (1) BH‑코어 형성 및 3체 이중성 생성 메커니즘, (2) 경화와 탈출 과정, (3) 내부·외부 합병 비율, (4) 검출률 추정, (5) 금속성 의존성 등 다각적인 측면에서 별집단 내 BH 동역학을 체계적으로 정량화하였다. 향후 연구에서는 초기 질량 함수, 원시 이중성 포함, 은하 조석장 효과 등을 추가해 보다 현실적인 모델링이 필요하다.