[论文解读] Numerical simulations of black hole accretion flows

Agnieszka Janiuk, Kostas Sapountzis+2 more|arXiv (Cornell University)|May 29, 2018

Gamma-ray bursts and supernovae参考文献 25被引用 3

一句话总结

本文提出了一套广义相对论磁流体动力学（GRMHD）数值模拟框架，用于模拟伽马射线暴（GRBs）和耀变体等极端天体物理环境中黑洞吸积盘及喷流的形成。结果表明，磁驱动的喷流（来自超吸积盘）、中微子过程以及Blandford–Znajek机制共同作用，可驱动相对论喷流并解释千新星辐射，喷流的洛伦兹因子可达约100，X射线辐射来自放射性同位素如$^{44}$Ti和$^{57}$Co。

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）验证模型，以约束吸积盘-喷流系统的物理特性。

提出的方法

• 采用保守型有限体积GRMHD代码，在弯曲时空下求解完整的相对论磁流体动力学方程组进行数值模拟。
• 模拟框架中引入了适用于GRB和千新星环境的热、致密且简并等离子体的真实状态方程（EOS）。
• 实现了核反应网络，以在非平衡条件下计算r过程同位素（包括$^{44}$Ti、$^{57}$Co和$^{65}$Zn）的合成。
• 使用后处理工具计算放射性衰变产生的合成X射线和光学辐射，以便与观测结果比较。
• 采用MPI并行化代码，并针对高性能计算（HPC）集群进行优化，报告了性能基准测试结果。
• 仔细处理边界条件和数值稳定性，以确保湍流磁化吸积盘长期演化模拟的准确性。

实验结果

研究问题

• RQ1 GRMHD模拟在超吸积盘条件下能否再现短伽马射线暴和千新星的观测特征？
• RQ2 在GRB引擎中，Blandford–Znajek机制与中微子驱动喷流对喷流能量的相对贡献如何？
• RQ3 微观物理过程（特别是核反应网络和状态方程）在多大程度上决定了吸积盘喷出物的成分与能量预算？
• RQ4 模拟的放射性同位素X射线辐射能否与GRB余晖（如GRB 130603B）中观测到的通量相匹配？
• RQ5 吸积盘湍流时标与喷流变异性时标如何与耀变体和GRB辐射中观测到的短时标变异性相关联？

主要发现

• 模拟结果表明，来自吸积盘的磁驱动、电子分数较低（$Y_{\rm e}$）的喷流可产生千新星辐射，与GW170817的观测结果一致。
• 喷流的最终洛伦兹因子可达约100量级，与超相对论性GRB喷流的观测结果一致。
• 预测在12–80 keV能量 band内，放射性同位素如$^{44}$Ti、$^{57}$Co和$^{65}$Ga会产生X射线辐射，与GRB 130603B中观测到的红外过量辐射相匹配。
• 发现Blandford–Znajek机制与中微子过程在GRB引擎中对喷流能量的贡献相当。
• 吸积盘的动态喷出物（质量$M_{\rm ej} \sim 0.01 M_{\odot}$）可在约1周的时间尺度内，在光学与近红外波段发射$10^{40}-10^{41}$ erg/s的辐射。
• 该模型能良好再现GW170817电磁对应体的光曲线，其快速衰减、偏红的辐射特征与r过程核合成一致。