[논문 리뷰] Numerical simulations of black hole accretion flows

Agnieszka Janiuk, Kostas Sapountzis외 2명|arXiv (Cornell University)|2018. 05. 29.

Gamma-ray bursts and supernovae참고 문헌 25인용 수 3

한 줄 요약

이 논문은 감마선 폭발(GRBs)과 블라자르와 같은 극한 천체물리 환경에서 블랙홀 순환 원반과 제트 형성 모델링을 위한 일반 상대성 이론 기반 자기유체역학(GRMHD) 수치 시뮬레이션 프레임워크를 제시한다. 이는 초순환 원반에서 유도된 자기장 기반 바람, 중성자율 과정, 그리고 Blandford–Znajek 메커니즘이 초상대론적 제트를 공급하고, 킬로노바 방출을 설명할 수 있음을 시사하며, 러스트 인자(Lorentz factor)는 약 100에 이르며, $^{44}$Ti와 $^{57}$Co와 같은 방사성 동위원소로부터 기인한 X선 방출이 관측됨을 보여준다.

ABSTRACT

We model the structure and evolution of black hole accretion disks, and their neighboring regions, using numerical simulations. The numerics is governed by the equations of general relativistic magneto-hydrodynamics (GRMHD). In particular, such disks and outflows can be found at the base of very energetic ultra-relativistic jets produced by cosmic explosions, so called gamma-ray bursts (GRBs). Another, more persistent type of the jet phenomena, are blazars, emitted from the centers of galaxies. Long-lasting, detailed computations are essential to properly determine the physics of these explosions, and confront the theoretical models with any potential observables. From the point of view of numerical methods and computational techniques, three ingredients need to be considered. First, the numerical scheme must work in a conservative manner, which is achieved by solving a set of non-linear equations at each time-step, to advance the conserved quantities from one time step to the next. Second, the efficiency of computations intrinsically depends on the code parallelization methods, which may use various techniques. Third, the analysis of results is possible via the post-processing of the computed time-dependent physical quantities, and visualization of the flow properties. This is done via implementing various packages and libraries that are standardized in the field of computational astrophysics and supported by community developers. In the present paper, we discuss the physical picture of the cosmic sources which are modeled using numerical framework. We also describe several technical issues, in the particular context of our own experience with the performance of the GRMHD code which we develop. We also present a suite of performance tests, done on the High-Performance Computer cluster (HPC) in the Center for Mathematical Modeling of the Warsaw University.

연구 동기 및 목표

• 감마선 폭발(GRBs)과 블라자르와 같은 극한 환경에서 블랙홀 주위의 자화된 순환 원반의 구조와 진화를 모델링하기 위해.
• 일반 상대성 이론 기반 자기유체역학(GRMHD)이 블랙홀 엔진으로부터 초상대론적 제트를 유도하는 데 기여하는 방식을 조사하기 위해.
• 미세물리학, 특히 핵반응, 상태방정식(EOS), 중성자율 수송이 원반 구조와 분출 물질 조성에 미치는 영향을 검토하기 위해.
• 모의 원반 성질과 관측 가능한 서명, 즉 방사성 붕괴로부터 기인한 X선 방출과 킬로노바 광도 곡선을 연결하기 위해.
• GW170817과 GRB 130603B를 포함한 다중 메시저 관측 결과와의 일치를 통해 디스크-제트 시스템 성질을 제약하기 위해.

제안 방법

• 곡선 시공간에서 전체 상대론적 MHD 방정식을 해결하는 보존형 유한체적 GRMHD 코드를 사용하여 수치 시뮬레이션을 수행한다.
• GRB 및 킬로노바 환경에 적합한 고온, 고밀도, 디제너레이트 플라즈마에 대한 현실적인 상태방정식(EOS)을 프레임워크에 통합한다.
• 비평형 조건 하에서 r-과정 동위원소, 즉 $^{44}$Ti, $^{57}$Co, $^{65}$Zn의 합성 계산을 위해 핵반응 네트워크를 구현한다.
• 후처리 도구를 사용하여 방사성 붕괴로부터 유도된 합성 X선 및 가시광선 방출을 계산하여 관측 결과와의 비교를 가능하게 한다.
• MPI를 사용하여 코드를 병렬화하고 고성능 계산(HPC) 클러스터에 최적화하였으며, 성능 벤치마크 결과를 보고한다.
• 장기적인 난류성, 자기장이 가득한 순환 원반의 정확한 진화를 확보하기 위해 경계 조건과 수치적 안정성에 주의 깊게 대응한다.

실험 결과

연구 질문

• RQ1 초순환 원반의 GRMHD 시뮬레이션은 짧은 감마선 폭발과 킬로노바의 관측된 성질을 어떻게 재현하는가?
• RQ2 GRB 엔진에서 초상대론적 제트를 공급하는 데 Blandford–Znajek 메커니즘과 중성자율 기반 바람 중 어느 것이 더 기여하는가?
• RQ3 핵반응 네트워크와 상태방정식(EOS)을 포함한 미세물리학이 원반 분출 물질의 조성과 에너지 예산에 얼마나 큰 영향을 미치는가?
• RQ4 방사성 동위원소로부터 기인한 모의 X선 방출이 GRB 후광에서 관측된 광도와 일치하는가, 예를 들어 GRB 130603B의 경우?
• RQ5 원반 난류의 시간스케일과 제트의 변동성은 블라자르 및 GRB 방출에서 관측된 짧은 시간스케일 변동성과 어떻게 관련되는가?

주요 결과

• 시뮬레이션 결과, 원반에서 유도된 자기장 기반 바람, 특히 낮은 전자율($Y_{\rm e}$) 바람이 킬로노바 방출을 공급하며, GW170817의 관측 결과와 일치한다.
• 제트의 최종 러스트 인자는 관측 결과와 일치하는 약 100 수준에 이른다.
• 12–80 keV 에너지 대역에서 $^{44}$Ti, $^{57}$Co, $^{65}$Ga와 같은 방사성 동위원소로부터 기인한 X선 방출이 예측되며, 이는 GRB 130603B에서 관측된 초과 NIR 방출과 일치한다.
• Blandford–Znajek 메커니즘과 중성자율 과정이 GRB 엔진에서 제트 에너지 공급에 비슷한 기여를 한다고 밝혀졌다.
• 원반에서 기인한 역동적 분출 물질($M_{\rm ej} \sim 0.01 M_{\odot}$)은 약 1주일 동안 옵티컬 및 가까운 적외선 대역에서 $10^{40}-10^{41}$ erg/s의 방출을 가능하게 한다.
• GW170817의 전자기 대응체 광도 곡선은 모델에 의해 잘 재현되며, r-과정 핵합성과 일치하는 빠른 퇄소와 빨간빛을 띤 방출 특성을 보인다.