[논문 리뷰] Hydrodynamical Evolution of Black-Hole Binaries Embedded in AGN Discs

Rixin Li, Dong Lai|arXiv (Cornell University)|2022. 02. 15.

Astrophysical Phenomena and Observations인용 수 4

한 줄 요약

이 연구는 활성은하핵(AGN) 디스크에 매달린 별질량 블랙홀 이진체(BBHs)의 진화를 모의하기 위해 시어링박스 프레임워크 내에서 고해상도 2차원 유체역학 시뮬레이션을 사용한다. 원형에 가까운 질량비가 유사한 BBHs는 질량 두 배 증가율의 몇 배 수준의 시간스케일로 빠르게 수축하며, 이는 빠른 하드닝을 의미한다. 이 Gauss형 이진체는 이심률 감쇠를 경험하며, 모든 현실적인 매개변수 조건에서 이동보다 이론적으로 더 빠른 궤도 감쇠를 보인다. 또한 질량 공급은 매우 변동성이 크며, 이심률이 증가함에 따라 주기적 변동 주파수는 이진체 궤도 주파수에서 반경 방향 진동 주파수로 이동한다.

ABSTRACT

Stellar-mass binary black holes (BBHs) embedded in active galactic nucleus (AGN) discs are possible progenitors of black-hole mergers detected in gravitational waves by LIGO/VIRGO. To better understand the hydrodynamical evolution of BBHs interacting with the disc gas, we perform a suite of high-resolution 2D simulations of binaries in local disc (shearing-box) models, considering various binary mass ratios, eccentricities and background disc properties. We use the $γ$-law equation of state and adopt a robust post-processing treatment to evaluate the mass accretion rate, torque and energy transfer rate on the binary to determine its long-term orbital evolution. We find that circular comparable-mass binaries contract, with an orbital decay rate of a few times the mass doubling rate. Eccentric binaries always experience eccentricity damping. Prograde binaries with higher eccentricities or smaller mass ratios generally have slower orbital decay rates, with some extreme cases exhibiting orbital expansion. The averaged binary mass accretion rate depends on the physical size of the accretor. The accretion flows are highly variable, and the dominant variability frequency is the apparent binary orbital frequency (in the rotating frame around the central massive BH) for circular binaries but gradually shifts to the radial epicyclic frequency as the binary eccentricity increases. Our findings demonstrate that the dynamics of BBHs embedded in AGN discs is quite different from that of isolated binaries in their own circumbinary discs. Furthermore, our results suggest that the hardening timescales of the binaries are much shorter than their migration timescales in the disc, for all reasonable binary and disc parameters.

연구 동기 및 목표

• LIGO/Virgo 중력파 병합의 잠재적 형성 채널인 기체로 이루어진 AGN 디스크에 매달린 별질량 블랙홀 이진체(BBHs)의 유체역학적 진화를 이해하기 위해.
• 이전 연구에서 BBHs가 AGN 디스크에서 하드닝되는지 확장되는지에 대한 모순된 결과가 있었는데, 이는 수치 해상도와 질량 공급 모델링의 차이 때문이므로 이를 해결하기 위해.
• 이중체 질량비, 이심률, 디스크 특성가 BBHs의 질량 공급률, 토크, 궤도 진화에 미치는 영향을 정량화하기 위해.
• AGN 디스크 내 BBHs의 하드닝 시간스케일이 이동 시간스케일보다 짧은지 여부를 판단하기 위해, 이는 병합 형성 효율성에 있어 핵심적인 질문이다.

제안 방법

• 중심 초거대 블랙홀 주변 AGN 디스크의 국소 영역을 모의하기 위해 시어링박스 근사법을 사용한 고해상도 2차원 유체역학 시뮬레이션을 수행한다.
• γ-법 상태방정식을 사용하고, 이중체에 대한 질량 공급률, 토크, 에너지 이행률을 계산하기 위해 견고한 후처리 방법을 적용한다.
• 작은 중력 연화 길이(~이중체 간격의 0.08배)와 높은 공간 해상도를 사용하여 각 이중체 성분 주변의 원형단일디스크(CSDs)를 정확히 해상화한다.
• 질량비(q = 0.1에서 1.0), 이심률(e = 0.0에서 0.5), 디스크 매개변수(h = 0.01, λ = 2.5)를 변화시킨 이중체를 시뮬레이션하며, 최종 300개 궤도 주기 동안 시간 평균된 양상을 분석한다.
• 이중체 궤도 진화를 반지름, 이심률, 이중체 특별 에너지의 시간 평균 변화율을 통해 평가한다.
• 일관된 후처리 프레임워크를 적용하여 시뮬레이션 출력물에서 질량 공급률과 토크를 추출하고, 장기적 진화에 대한 신뢰성 있는 추정을 보장한다.

실험 결과

연구 질문

• RQ1 AGN 디스크 내 BBHs의 궤도 감쇠 속도가 이동 시간스케일을 초과하는가? 어떤 조건에서 그러한가?
• RQ2 이중체 질량비는 AGN 디스크 내 BBHs의 질량 공급률과 궤도 진화에 어떤 영향을 미치는가?
• RQ3 이중체 이심률이 질량 공급의 주요 변동 주파수와 이심률 감쇠 속도를 결정하는 데 어떤 역할을 하는가?
• RQ4 적절히 해상화된 CSDs를 포함한 고해상도 시뮬레이션은 이전에 보고된 저해상도 연구에서의 BBHs의 확장 현상과 어떻게 다를까?
• RQ5 원형단일디스크의 열역학적 상태(예: 온도)가 순수 토크에 영향을 미치며, 이는 이중체의 하드닝 또는 확장에 영향을 미치는가?

주요 결과

• 원형에 가까운 질량비가 유사한 BBHs는 빠르게 수축하며, 질량 두 배 증가율의 몇 배 수준의 궤도 감쇠 속도를 보여, 빠른 하드닝을 시사한다.
• 이심률이 있는 이중체는 항상 이심률 감쇠를 경험하며, 이심률이 높거나 질량비가 작은 후행 이중체는 궤도 감쇠 속도가 느리며, 일부 경우에서는 궤도 확장까지 관찰된다.
• 시간 평균 질량 공급률은 수급체의 물리적 크기에 의존하며, 더 큰 수거 반경일수록 더 높은 질량 공급률을 보인다.
• 원형 이중체의 경우 질량 공급 변동은 고정된 기준에서의 이중체 궤도 주파수에 의해 지배되지만, 이심률이 증가함에 따라 주파수는 반경 방향 진동 주파수로 이동한다.
• 모든 현실적인 매개변수 조합에서 AGN 디스크 내 BBHs의 하드닝 시간스케일은 이동 시간스케일보다 현저히 짧으며, 이는 빠른 병합 형성에 유리하다.
• 결과는 AGN 디스크 내 BBH 역학이 원형이중체 디스크 내 고립된 이중체와 본질적으로 다르며, 더 강한 기체 토크와 독특한 질량 공급 서명을 가짐을 보여준다.