[논문 리뷰] Merging of unequal mass binary black holes in non-axisymmetric galactic nuclei

Peter Berczik, Manuel Arca Sedda외 4명|arXiv (Cornell University)|2020. 08. 10.

Galaxies: Formation, Evolution, Phenomena참고 문헌 55인용 수 8

한 줄 요약

이 연구는 ϕ-GPU 코드를 사용한 고정밀 N체 시뮬레이션을 통해 비축대칭이고 회전하는 은하핵에서 질량이 다른 초거대 블랙홀 이중성(SMBHB)의 역학적 경화 및 융합을 조사한다. 유의미한 정체(즉, '마지막 파리세크 문제' 없음)가 관찰되지 않았으며, 비대칭 이중성은 등질량 이중성보다 최대 10,000배 빠르게 융합되며, 경화는 회전하는 시스템에서 별의 산란에 의해 효율적으로 일어난다.

ABSTRACT

In this work, we study the stellar-dynamical hardening of unequal mass supermassive black hole (SMBH) binaries in the central regions of merging galactic nuclei. We present a comprehensive set of direct $N$-body simulations of the problem, varying both the total mass and the mass ratio of the SMBH binary (SMBHB). Simulations were carried out with the $\varphi-$GPU $N$-body code, which enabled us to fully exploit supercomputers equipped with graphic processing units (GPUs). As a model for the galactic nuclei, we adopted initial axisymmetric, rotating models, aimed at reproducing the properties of a galactic nucleus emerging from a galaxy merger event, containing two SMBHs which were unbound initially. We found no 'final-parsec problem', as our SMBHs tend to pair and shrink without showing significant signs of stalling. This confirms earlier results and extends them to large particle numbers and rotating systems. We find that the SMBHB hardening depends on the binary-reduced mass ratio via a single parameter function. Our results suggest that, at a fixed value for the SMBHB primary mass, the merger time of highly asymmetric binaries is up to four order of magnitudes smaller than the equal-mass binaries. This can significantly affect the population of SMBHs potentially detectable as gravitational wave sources.

연구 동기 및 목표

• 은하 융합의 결과로 발생하는 비축대칭이고 회전하는 은하핵에서 질량이 다른 초거대 블랙홀 이중성(SMBHB)의 별역학적 경화를 조사하기 위해.
• 이전 시뮬레이션을 확장하여 질량 비율(q = 0.01에서 1까지)의 넓은 범위를 탐색하고, 은하핵 모델에 회전을 포함하기 위해.
• 대규모 질량 비대칭이 있는 비축대칭이고 회전하는 시스템에서 '마지막 파리세크 문제'—즉, 경화의 정체 가능성—가 지속되는지 여부를 규명하기 위해.
• 융합 시간스케일이 이중성 질량 비율과 총 질량에 어떻게 의존하는지 정량화하고, 중력파 탐지에 대한 함의를 평가하기 위해.
• 총질량과 감소질량을 포함하는 두 매개변수 스케일링 법칙의 적용 가능성을 검증하고 확장하기 위해.

제안 방법

• ϕ-GPU N체 코드를 사용하여 직접 N체 시뮬레이션을 수행하였으며, 이는 GPU 가속 계산을 활용해 최대 10^6 개의 입자를 포함한 대규모 입자 시스템을 시뮬레이션할 수 있다.
• 은하핵을 초기에는 축대칭이고 회전하는 시스템으로 모델링하였으며, 은하 융합의 결과로 나타나는 두 개의 비결합 초거대 블랙홀을 반영하였다.
• 후뉴턴 보정을 방정정식 운동에 포함시켜, 후기 단계에서 중력파 방출 효과를 자동으로 포함시켰다.
• 다양한 총질량(10^4–10^9 M⊙)과 질량 비율(q = 0.01–1)을 가진 SMBHB를 여러 시뮬레이션에 걸쳐 변화시켜, 경화 시간스케일을 매핑하였다.
• 시간에 따른 이중성 간격, 유도 에너지, 그리고 경화 속도를 추적하여, 시스템 매개변수에 따라 변화하는 융합 시간스케일 T_merge 를 유도하였다.
• 감소질량 µ 와 총질량 M12 를 통해 차원 없는 경화 속도 스케일링을 사용하여, 다양한 질량 비율에서 통합적인 스케일링 관계를 식별하였다.

실험 결과

연구 질문

• RQ1 비축대칭이고 회전하는 은하핵에서 질량이 다른 SMBHB를 지닌 시스템에서 '마지막 파리세크 문제'가 지속되는가?
• RQ2 질량 비율(q = m2/m1)이 밀도가 높고 회전하는 별 시스템에서 SMBHB의 경화 속도와 융합 시간스케일에 어떻게 영향을 미치는가?
• RQ3 은하핵 내의 회전 존재가 SMBHB의 별역학적 경화를 강화하거나 억제하는 정도는 어느 정도인가?
• RQ4 총질량 M12 와 감소질량 µ 를 포함하는 두 매개변수 스케일링 법칙으로 경화 시간스케일을 기술할 수 있는가?
• RQ5 매우 비대칭적인 SMBHB(예: q = 0.01)의 예측 융합 시간스케일은 등질량 이중성(q = 1)과 비교해 어떻게 되며, 이는 중력파 소스로서의 탐지 가능성에 어떤 영향을 미치는가?

주요 결과

• 시뮬레이션에서 유의미한 정체(즉, '마지막 파리세크 문제' 없음)가 관찰되지 않았으며, 비축대칭이고 회전하는 핵에서 SMBHB는 계속해서 효율적으로 경화되고 융합된다.
• 고정된 주요 SMBH 질량에서, 보조 질량 비율 q 가 1에서 0.01으로 감소함에 따라 융합 시간스케일 T_merge 가 최대 4개의 자리수 감소한다.
• 경화 속도는 감소질량 µ 와 총질량 M12 에 따라 스케일링되며, 다양한 질량 비율에서 데이터를 통합하는 데 유용한 두 매개변수 스케일링 관계를 지지한다.
• 10^4–10^9 M⊙ 의 구성 질량과 질량 비율 q = 10^{-4}–1 을 가진 SMBHB 는 효율적인 별역학적 경화로 인해 T_merge ≈ 10^3–10^7 년의 시간스케일로 융합된다.
• 매우 비대칭 이중성(q ≈ 0.01)의 융합 시간스케일은 등질량 이중성(q = 1)보다 최대 10,000배 짧으며, 이는 중력파 소스로서의 탐지 가능성의 급격한 증가를 의미한다.
• 결과는 이전의 발견을 더 큰 입자 수와 회전하는 시스템으로 확장하여, 비대칭이고 회전하는 핵에서도 별경화가 여전히 효과적으로 유지됨을 확인한다.