[논문 리뷰] The erratic path to coalescence of LISA massive black hole binaries in sub-pc resolution simulations of smooth circumnuclear gas disks

Rafael Souza Lima, Lucio Mayer외 3명|arXiv (Cornell University)|2020. 03. 30.

Astrophysical Phenomena and Observations참고 문헌 1인용 수 4

한 줄 요약

이 연구는 부드러운 중심핵 가스 디스크(CND) 내에서 거대 블랙홀(MBH) 이중성의 궤도 감쇠를 고해상도 유체역학 시뮬레이션을 통해 조사하며, 디스크 토크에 의해 유도되는 불규칙한 감쇠 과정을 규명한다. 주요 발견은 하위-파세크(sub-pc) 감쇠가 비효율적이고 취약하며, 마지막 느린 단계가 존재함으로써 원형 이중성 디스크(CBD) 형성이 고해상도 시뮬레이션에서만 가능하다는 점으로, 이는 가스가 풍부한 환경에서 융합을 위해 항성 3체 상호작용이 필수적일 수 있음을 시사한다.

ABSTRACT

We report on high-resolution simulations that explore the orbital decay of massive black hole (MBH) pairs with masses between $10^5$ and $10^7 M_{\odot}$ embedded in a circumnuclear gas disk (CND). An adiabatic equation of state is adopted, with a range of adiabatic indices, which maintains a smooth flow. Mergers between MBHs in this mass range would be detectable by the upcoming Laser Inteferometer Space Antenna (LISA). The orbital evolution is followed from the CND scale ($100$~pc) down to separations of $0.1$--$0.01$~pc at which a circumbinary disk (CBD) could form. The decay is erratic and strongly dependent on the gas flow within the disk, that ultimately determines the net torques experienced by the sinking MBH. Overall, we can identify three different evolutionary stages: (i) an initially slow decay that leads to no significant change in the orbital angular momentum, resulting in some circularization; (ii) a fast migration phase in which the orbital angular momentum decreases rapidly; and (iii) a final, very slow decay phase, in which orbital angular momentum can even increase, and a CBD can form. The fast migration phase owes to disk-driven torques originating primarily from the co-orbital region of the secondary MBH, at a distance of 1--3 Hill radii. We find strong analogies with fast Type III migration for massive planets in protoplanetary disks. The CBD forms only when the decay rate becomes small enough to allow it enough time to carve a cavity around the primary MBH, at scales $\lesssim 1$~pc; when this happens, the MBH separation nearly stalls in our higher-resolution run. We suggest an empirically modified gap-opening criterion that takes into account such timescale effects as well as other deviations from standard assumptions made in the literature. [Abriged]

연구 동기 및 목표

• 100 pc에서 0.01 pc까지의 중심핵 가스 디스크(CND) 내 거대 블랙홀(MBH) 이중성의 궤도 감쇠를 조사하기 위해.
• 하위-파세크 분리에서 MBH 이중성 이동의 효율성과 메커니즘에 영향을 미치는 가스 디스크 구조와 해상도의 영향을 규명하기 위해.
• 중력파 방출의 핵심 전구단계인 원형 이중성 디스크(CBD) 형성이 가능한 조건을 평가하기 위해.
• 고해상도, 매끄러운 CND 시뮬레이션에서 표준 갭 열림 기준의 한계를 평가하고, 경험적으로 수정된 기준을 제안하기 위해.
• '최종 파세크 문제'에 대한 함의와 항성 3체 상호작용과 같은 추가적인 경화 메커니즘의 역할을 탐색하기 위해.

제안 방법

• 10⁵에서 10⁷ M☉ 사이의 질량을 가진 매끄럽고 압력 지지 CND 내 MBH 이중성의 고해상도 스무스한 입자 유체역학(SPH) 시뮬레이션을 수행하였다.
• 매끄러운 가스 유동을 유지하고 디스크 온도 및 압력을 제어하기 위해 다양한 다이아바틱 지수를 가진 등온 상태 방정식을 사용하였다.
• 약 100 pc에서 0.01 pc까지 궤도 진화를 추적하며, 토크 기반 이동과 각운동량 진화에 집중하였다.
• 이동 시간스케일과 비이상적 조건을 고려하기 위해 경험적으로 수정된 기준을 위해 현상학적 보정 인자 g를 도입하였다.
• 고해상도 및 저해상도 런 간의 결과를 비교하여, 공동 및 CBD 형성에 영향을 미치는 해상도 의존성 효과를 분리하였다.
• 특히 1–3 힐 반경 범위 내 공궤도 영역에서의 토크 기여를 분석하여 빠른 이동 단계의 기원을 규명하였다.

실험 결과

연구 질문

• RQ1 매끄러운 CND 내 MBH 이중성의 궤도 감쇠는 약 100 pc에서 하위-파세크 스케일로 어떻게 진화하며, 각 단계에서 감쇠를 이끄는 원리는 무엇인가?
• RQ2 어떤 조건에서 원형 이중성 디스크(CBD)가 형성되며, 해상도는 이 과정에 어떻게 영향을 미치는가?
• RQ3 왜 하위-파세크 분리에서 감쇠 속도가 비효율적이고 불규칙해지며, 마지막 느린 단계 또는 정지 상태를 유도하는 물리적 메커니즘은 무엇인가?
• RQ4 특히 제2 MBH의 공궤도 영역에서 기인하는 디스크 기반 토크는 어떻게 빠른 이동 단계를 촉진하는가?
• RQ5 고해상도, 매끄러운 CND 시뮬레이션에서 표준 갭 열림 기준을 신뢰성 있게 적용할 수 있는가, 아니면 경험적으로 조정이 필요한가?

주요 결과

• 궤도 감쇠는 세 단계로 나뉘어 진행된다: 초기 느린 단계(원형화), 1–3 힐 반경 범위 내 공궤도 토크에 의해 유도되는 빠른 이동 단계, 그리고 각운동량가 오히려 증가할 수 있는 마지막 느린 단계.
• 빠른 이동 단계는 주로 제2 MBH의 공궤도 영역 내 가스에서 기인하는 토크에 의해 유도되는 타입 III 행 星 이동과 유사하다.
• 원형 이중성 디스크(CBD)는 감쇠 속도가 충분히 느려져 주요 MBH 주변에 ≤1 pc 스케일에서 공동 형성이 가능해지는 고해상도 시뮬레이션에서만 형성된다.
• 저해상도 런에서는 빠른 감쇠로 인해 갭 및 CBD 형성이 억제되어 더 빠른 궤도 감쇠와 정지 없이 진행된다.
• 표준 갭 열림 기준은 시간스케일 효과로 인해 이러한 시뮬레이션에서 실패한다; 경험적으로 수정된 기준(보정 인자 g를 포함)은 감쇠 속도가 느려질 경우 CBD 형성을 성공적으로 예측한다.
• CBD 및 공동이 존재하는 동안에도 MBH 이중성은 비제로 이심률을 유지하며, 항성 3체 상호작용은 하위-파세크 스케일에서 여전히 타당한 경화 메커니즘이다.