[논문 리뷰] Accretion onto a Supermassive Black Hole Binary Before Merger

Mark Avara, Julian H. Krolik외 4명|arXiv (Cornell University)|2023. 05. 29.

Black Holes and Theoretical Physics인용 수 8

한 줄 요약

이 논문은 equal-mass, non-spinning supermassive binary black holes의 inspiraling 과정에서 20M에서 ~9M까지 진행되는 최초의 relativistic 3D GRMHD 시뮬레이션을 제시하며, 지속적인 accretion이 약간의 감소를 보이고 풍부한 전자기 서명을 나타냅니다.

ABSTRACT

While supermassive binary black holes inspiral toward merger they may also experience significant accretion of matter from a surrounding disk. We study the dynamics of this system, simultaneously describing the evolving spacetime and magnetized plasma, and present the first relativistic calculation simulating two equal-mass, non-spinning black holes as they inspiral from a $20M$ ($G=c=1$) initial separation almost to merger, $\simeq 9M$ ($M$=binary mass). Our dynamical results imply important observational consequences: for instance, the accretion rate $\dot M$ onto the black holes first decreases and then reaches a plateau, dropping by only a factor of $\sim 3$ despite the rapid inspiral. An estimated bolometric light curve thus suggests some merging SMBBHs may be quite luminous past the predicted decoupling from the circumbinary disk. The minidisks through which the accretion reaches the black holes are very non-standard: Reynolds, not Maxwell, stresses dominate, and they oscillate between two states. In one part of the cycle, ``sloshing" streams transfer mass from one minidisk to the other through the L1 point at a rate $\sim 0.1 imes$ the accretion rate, carrying kinetic energy at a rate that can be as large as the peak minidisk bolometric luminosity. We also discover that episodic accretion drives minidisks with time-varying tilts. The unsigned poloidal magnetic flux on the black hole event horizon is roughly constant at a dimensionless level $ϕ\sim 2-3$, but doubles just before merger; if the black holes had significant spin, this flux indicates the potential for powerful jets with variability driven by binary dynamics, another prediction of potentially unique EM signatures. This simulation is the first to employ our multipatch infrastructure \pwmhd, decreasing computational expense to $\sim 3\%$ of conventional single-grid methods' cost.

연구 동기 및 목표

• 합병을 향해 inspiral하는 이진 블랙홀 시스템에 대한 가스 역학 및 물질 흡수를 조사한다.
• 써클럼 바이너리 디스크가 미니디스크 및 이진체에 어떻게 공급되는지, 시간 의존적 물질 흡수율을 포함하여 결정한다.
• 후기 단계의 inspiral 및 approaching merger와 연관된 관찰 가능한 전자기 서명을 식별한다.
• 흡수 및 방출에 대한 자기장, 토크, 그리고 디스크 불안정성의 역할을 탐구한다.

제안 방법

• 일반 상대론적 MHD 방정식을 이진체에 대한 2.5PN/3.5PN 시공에 맞춰 Schwarzschild 패치를 매칭하는 시간 의존적 시공에서 해결한다.
• 좌표 특이점을 피하고 계산 비용을 줄이기 위해 두 패치(외부 구면 CBD 영역과 내부 직교 좌표 미니디스크 영역)를 갖는 새로운 다중 패치 인프라스트럭처(PatchworkMHD)를 사용한다.
• 이상 MHD와 목표 엔트로피를 기반으로 하는 냉각 규칙으로 방출되는 열을 방사하는 magnetized 가스를 진화시킨다.
• Constrained transport(FluxCT) 및 패치 경계에서의 발산 제거 루틴으로 자기장 발산 오차를 최소화한다.
• 사전 평형 상태의 RunSE 상태에서 circumbinary disk를 초기화하고, 두 패치 구성을 맞추기 위해 내부 영역을 수정하며 ~36 개의 이진 궤도를 실행한다(20M에서 ~9M).
• 수렴성을 테스트하고 단일 m mesh 실행과의 비교를 위해 여러 해상도(PM.IN20s 및 PM.IN20sHR)를 제공한다.

실험 결과

연구 질문

• RQ1 20M에서 ~9M으로 축약되면서 이진 억지에서 circumbinary disk가 두 미니디스크로, 가스가 어떻게 흐르는가를 측정한다?
• RQ2 각 블랙홀 및 이진 전체에 대한 시간 의존적 흡수율은 어떻게 되고, 분리 감소와 함께 어떻게 변하는가?
• RQ3 후기 단계의 inspiral 동안 미니디스크, 흘러가는 스트림, 디스크 경사에 의해 어떤 전자기 서명이 발생하는가?
• RQ4 블랙홀 지평선의 자기 플럭스가 어떻게 진화하며, 합병 중 제트 생성에 어떤 영향을 주는가?

주요 결과

• 이진체에 대한 총 흡수율은 분리 거리가 줄어들면서 약 3–4배 감소하지만 최종 단계까지는 대략 일정하게 유지된다.
• 볼로메트릭 광도 곡선은 흡수율과 동일한 패턴을 따르는 경향이 있어 이중체의 진화에도 불구하고 합병 근처에서 상당한 EM 방출이 가능하다는 시사점을 제공한다.
• 미니디스크는 Maxwell 응력보다 Reynolds 응력에 의해 지배되며 순환 상태 변화와 주기적 상태 변화를 보이고, 약 10%의 흡수율로 질량을 갭 너머로 운반하는 ‘sloshing’ 스트림이 발생하여 상당한 운동에너지를 전달한다.
• 슬로싱에 의한 질량 전달은 충격에 에너지를 침착시키는 역할을 할 수 있으며, 미니디스크 광도에 필적하는 소멸률을 만들어낸다.
• 주기적 흡수가 미니디스크를 궤도면에 대해 시간에 따라 기울기를 갖도록 만들어 주기적 광도 곡선 특징에 기여한다.
• 지평선의 unsigned poloidal 자기 플럭스는 φ≈2–3을 유지하다가 합병 직전에 두 배로 증가하며, 회전이 크다면 강력한 제트를 만들 가능성을 시사한다.