[논문 리뷰] Dynamics and spin alignment in massive, gravito-turbulent circumbinary discs around supermassive black hole binaries

Martin A. Bourne, Davide Fiacconi외 3명|arXiv (Cornell University)|2023. 11. 28.

Astrophysical Phenomena and Observations인용 수 4

한 줄 요약

이 연구는 초거대 브라운홀 이중성계 주위의 거대하고 중력 난류가 있는 원형이중행성원반(CBDs)에 대한 3차원 유체역학 시뮬레이션을 수행하며, 슈퍼라그랑주안의 정밀화와 미세구역 셰카라-수냐에프 원반 모델을 사용하여 브라운홀 질량과 스핀의 자발적 진화를 추적한다. 이전 연구들과 달리, 보조 브라운홀로의 선호적 물질 붕괴가 없음을 발견하여 질량 비율이 일정하고 중력파 배경이 더 약해지며, 스핀 정렬은 낮은 질량 비율 시스템에서도 빠르게 발생하지만, 혼란스러운 붕괴 시나리오에서는 기울어진 미니원반으로 인해 정렬이 방해받을 수 있음을 시사한다.

ABSTRACT

Parsec-scale separation supermassive black hole binaries in the centre of gas-rich galaxy merger remnants could be surrounded by massive circumbinary discs (CBDs). Black hole mass and spin evolution during the gas-rich binary inspiral are crucial in determining the direction and power of relativistic jets that radio observations with LOFAR and SKAO will probe, and for predicting gravitational wave (GW) emission that IPTA and LISA will measure. We present 3D hydrodynamic simulations capturing gas-rich, self-gravitating CBDs around a $2 imes 10^6$ M$_{\odot}$ supermassive black hole binary, that probe different mass ratios, eccentricities and inclinations. We employ a sub-grid Shakura-Sunyaev accretion disc to self-consistently model black hole mass and spin evolution together with super-Lagrangian refinement techniques to resolve gas flows, streams and mini-discs within the cavity, which play a fundamental role in torquing and feeding the binary. We find that higher mass ratio and eccentric binaries result in larger cavities, while retrograde binaries result in smaller cavities. All of the simulated binaries are expected to shrink with net gravitational torques being negative. Unlike previous simulations, we do not find preferential accretion onto the secondary black hole. This implies smaller chirp masses at coalescence and hence a weaker GW background. Critically this means that spin-alignment is faster than the binary inspiral timescale even for low mass ratios. When considering initially misaligned systems, the orientation of the mini-discs around each black hole can vary significantly. We discuss the implications of this behaviour for black hole spin alignment and highlight the need for broader parameter space studies of misaligned systems to understand the impact on black hole recoil velocities.

연구 동기 및 목표

• 은하단 융합의 가스 풍부 단계 동안 거대하고 자가중력 작용이 있는 원형이중행성원반(CBDs)에서 초거대 브라운홀 이중성계의 역학과 스핀 진화를 이해하기 위해.
• 이중성계 질량 비율, 이심률, 기울기가 원형이중행성원반의 공동 구조, 기체 유입, 붕괴 패tern에 어떤 영향을 미치는지 조사하기 위해.
• 이전 연구에서 제안된 바와 같이 붕괴가 보조 브라운홀로 선호적으로 유입되는지 확인하고, 이를 초음파질량과 중력파 방출에 미치는 영향를 평가하기 위해.
• 특히 기울어진 시스템에서 브라운홀 스핀 정렬의 시간스케일이 이중성계의 수렴 시간스케일과 비교하여 어떻게 되는지 평가하고, 이를 통해 재진입 속도와 은하계 이동에 미치는 영향을 분석하기 위해.
• 미니원반과 유입 흐름의 역학이 브라운홀 스핀을 어떻게 토크하여 장기적인 이중성계 진화에 영향을 미치는지 평가하기 위해.

제안 방법

• 유동, 유입 흐름, 미니원반을 모두 해상도를 높여 해석할 수 있도록 초래그랑주안 정밀화를 사용한 3차원 유체역학 시뮬레이션을 수행하며, 메쉬 정밀화 코드를 사용한다.
• 점성 토크를 통해 브라운홀 질량과 스핀을 자발적으로 진화시키기 위해, 미세구역 셰카라-수냐에프 붕괴 원반 모델을 통합한다.
• 기체 열역학과 원반 안정성을 유지하기 위해, 느린 β 냉각을 포함한 등온 상태 방정식을 사용한다.
• 이중성계 매개변수를 체계적으로 변화시켜 질량 비율(q = 1, 1/3, 1/10), 이심률(e = 0, 0.5), 기울기(i = 0°, 45°, 180°)를 포함한 총 12개의 시뮬레이션 설정을 수행한다.
• 이중성계에 작용하는 중력 토크를 다양한 반경(R < a 및 R > a)에서 기체로부터 추적하여 순수한 수렴 속도를 결정한다.
• 이중성계 수렴 시간스케일과 비교하여 미니원반의 형성, 크기, 정렬 시간스케일을 분석하여 스핀 진화를 평가한다.

실험 결과

연구 질문

• RQ1 거대하고 중력 난류가 있는 원형이중행성원반에서 붕괴가 일부 이론 모델이 예측한 것처럼 보조 브라운홀로 선호적으로 공급되는가?
• RQ2 이중성계 질량 비율, 이심률, 기울기가 원형이중행성원반의 공동 크기와 구조, 기체 유입에 어떤 영향을 미치는가?
• RQ3 이중성계에 작용하는 순수 중력 토크는 무엇이며, 이는 수렴 시간스케일에 어떤 영향을 미치는가?
• RQ4 브라운홀 스핀은 시스템의 전체 운동량 방향과 얼마나 빨리 정렬되는가? 특히 초기에 기울어진 구성에서 이 정렬이 보장되는가?
• RQ5 기울어진 미니원반은 브라운홀 스핀 진화와 합성 후 잠재적인 재진입 속도에 어떤 영향을 미치는가?

주요 결과

• 높은 질량 비율과 이심률이 큰 이중성계는 더 큰 원형이중행성원반 공동을 형성하지만, 반대 방향 운동을 하는 이중성계는 더 강한 기체 토크로 인해 더 작은 공동을 형성한다.
• 모든 구성에서 이중성계에 작용하는 순수 중력 토크는 음수이며, 이는 수렴을 유도한다. 낮은 질량 비율과 반대 방향 이중성계는 더 빠르게 수축한다.
• 붕괴가 보조 브라운홀로 선호적으로 공급되지 않으며, 질량 비율은 약간의 변동만을 보이며, 이는 합성 시 더 낮은 초음파질량과 더 약한 중력파 배경을 의미한다.
• 미니원반의 효율적인 토크 작용으로 인해, 낮은 질량 비율 시스템에서도 이중성계 수렴 시간스케일과 유사하거나 이를 빠르게 정렬하는 스핀 정렬이 발생한다.
• 초기 기울어진 구성에서, 미니원반은 빠르게 정렬하지 않으며, 결과적으로 브라운홀 스핀은 합성 시 기울어진 채로 남을 수 있어, 상당한 재진입 속도와 은하계 이동 위험을 증가시킨다.
• 보조 브라운홀로의 선호적 붕괴가 없는 것은 등온 상태 방정식, 느린 β 냉각, 낮은 점성도, 기존의 싱크 파article 붕괴 모델과의 차이에서 기인한다.