젊은 대질량 및 개방성 클러스터에서의 블랙홀 병합과 중력파 관측 비교

초록

본 연구는 질량 10⁴–10⁵ M☉, 반지름 1–3 pc, 금속성 0.0002–0.02 범위의 90개 젊은 대질량·개방성 클러스터 모델을 NBODY7로 직접‑포스트뉴턴 연산해, 동역학적으로 형성된 BBH 병합을 조사한다. 모델은 관측된 LVK‑O4a 중력파 사건들의 질량·스핀·비율 분포와 일부 일치하나, 저질량(≈10 M☉) 프라이머리 부족을 보인다. 전체 병합률은 관측치의 25‑33 % 수준이며, 적색이동에 따라 별 형성률보다 빠르게 증가한다. 또한 모델은 Gaia에서 발견된 BH‑MS·NS‑MS 이진과도 매칭되는 파라미터를 만든다.

상세 분석

이 논문은 젊은 대질량 클러스터(YMC)와 그 진화 후 개방성 클러스터(OC)가 중력파(BBH) 발생에 기여할 수 있음을 정량적으로 검증한다. 90개의 모델은 초기 질량을 10⁴, 2×10⁴, 3×10⁴, 5×10⁴, 7.5×10⁴, 10⁵ M☉ 로 설정하고, 반지름 1, 2, 3 pc, 금속성 Z=0.0002, 0.001, 0.005, 0.01, 0.02 로 조합하였다. 각 클러스터는 은하 중심거리 8.5 kpc의 원형 궤도에 놓여 외부 MW 조석장을 받으며, 초기 별 질량은 Kroupa IMF(0.08–150 M☉)와 10 % 원시 이진 비율을 갖는다. 특히 O형 별(≥16 M☉)은 100 % 이진으로 설정해 관측된 고이진율을 반영하였다.

NBODY7은 4차 Hermite 적분과 KS·Chain 정규화를 이용해 별별 중력과 포스트-뉴턴 효과를 직접 계산한다. 별 진화는 BSE와 연계해 지연된 잔여질량 모델을 적용했으며, PPSN·PSN을 포함해 40–120 M☉ 사이의 전통적 질량갭을 재현한다. 초신성(Natal kick)은 Maxwellian σ=265 km s⁻¹에서 시작해 질량 낙하 비율에 따라 감쇠시켜, 8 M☉ 이상 BH는 대다수가 클러스터에 남고, 저질량 BH와 핵융합형 NS는 대부분 탈출한다. 전자포획 초신성(ECS)으로 형성된 NS는 낮은 킥을 받아 보존된다.

시뮬레이션 결과, 대부분의 BBH 병합은 클러스터 내부에서 동역학적으로 조립된 이중 블랙홀이며, 탈출 후 병합은 소수에 불과하고 주로 원시 이진에서 유래한다. 질량비 q는 0.1–0.2까지 내려가지만 평균적으로 q≈0.8에 집중된다. 병합 효율은 클러스터 질량에 비단조적이며, M₀≈7.5×10⁴ M☉에서 최고에 도달하고, M₀≤3×10⁴ M☉에서도 두 번째 피크를 보인다. 이는 중간 질량 클러스터가 충분히 높은 밀도와 충분한 BH 보유량을 동시에 갖추어 동역학적 조립을 촉진하기 때문이다.

모델에서 도출된 병합률은 z≈0.2에서 관측된 전체 BBH 병합률의 25–33 %에 해당한다. 적색이동이 증가함에 따라 병합률은 우주 별 형성률보다 약간 더 가파르게 상승한다는 점에서 LVK가 제시한 “전역적인 성장”과 일치한다. 스핀 분포는 z=0에서 양의 비대칭을 보이며, 고스핀·고질량 BH 쌍에서도 스핀 축이 크게 틀어진(미정렬) 사건을 재현한다. 이는 클러스터 내 다중 상호작용이 스핀 방향을 무작위화하기 때문이다.

또한, 시뮬레이션은 BH‑주계열성(또는 NS‑주계열성) 이진을 생성한다. 이들 이진은 질량, 궤도 주기, 이심률 면에서 Gaia DR3에서 보고된 후보와 일치한다. 따라서 YMC/OC 채널이 중력파 이외에도 관측 가능한 전자기 현상과 연결될 가능성을 제시한다.

한계점으로는 (1) 원시 이진 비율을 전체 10 %로 고정했으나 실제 클러스터에서는 환경에 따라 크게 변동할 수 있다. (2) 금속성 분포를 단순히 5가지 값으로만 샘플링했으며, 우주적 금속성 진화와 연계된 가중치가 부족하다. (3) 시뮬레이션 시간은 대부분 2–4 Gyr에 머물러, 10 Gyr 이상 살아남는 구형 클러스터와의 직접 비교가 제한적이다. 향후 연구에서는 더 넓은 파라미터 공간, 고해상도 Monte‑Carlo 보조, 그리고 관측 기반 금속성·이진 분포를 적용해 모델을 확장할 필요가 있다.