[论文解读] Jet Formation in Black Hole Accretion Systems II: Numerical Models

Jonathan C. McKinney|arXiv (Cornell University)|Jun 16, 2005

Gamma-ray bursts and supernovae被引用 3

一句话总结

本文通过引入正负电子对产生过程，对黑洞吸积中的广义相对论磁流体动力学（GRMHD）模拟，旨在模拟伽马射线暴（GRBs）、活动星系核（AGN）和X射线双星系统中的喷流形成。结果表明，在塌缩星模型中，由于磁场能量的提取，Poynting主导喷流可达到Γ ∼ 100–1000的终端洛伦兹因子；而含有重子污染的Poynting-重子喷流则达到Γ ∼ 1.5–3，喷流加速发生在大半径区域，其时间变异性由环向磁场不稳定性驱动。

ABSTRACT

In a companion theory paper, we presented a unified model of jet formation. We suggested that primarily two types of relativistic jets form near accreting black holes: a potentially ultrarelativistic Poynting-dominated jet and a Poynting-baryon jet. We showed that, for the collapsar model, the neutrino-driven enthalpy flux (classic fireball model) is probably dominated by the Blandford-Znajek energy flux, which predicts a jet Lorentz factor of $Γ\sim 100-1000$. We showed that radiatively inefficient AGN, such as M87, are synchrotron-cooling limited to $Γ\sim 2-10$. Radiatively efficient x-ray binaries, such as GRS1915+105, are Compton-drag limited to $Γ\lesssim 2$, but the jet may be destroyed by Compton drag. However, the Poynting-baryon jet is a collimated outflow with $Γ\sim 1-3$. Here we present general relativistic hydromagnetic simulations of black hole accretion with pair creation used to simulate jet formation in GRBs, AGN, and x-ray binaries. Our collapsar model shows the development of a patchy ``magnetic fireball'' with typically $Γ\sim 100-1000$ and a Gaussian structure. Temporal variability of the jet is dominated by toroidal field instabilities for $\gtrsim 10^2$ gravitational radii. A broader Poynting-baryon jet with $Γ\sim 1.5$ could contribute to a supernova.

研究动机与目标

• 统一解释伽马射线暴（GRBs）、活动星系核（AGN）和X射线双星系统等多样化黑洞吸积系统中相对论性喷流形成的物理机制。
• 基于最少的物理原理，确定喷流能量来源、组分、准直化以及终端洛伦兹因子的起源。
• 通过包含正负电子对产生和中子扩散的数值模拟，检验McKinney（2005b）的理论框架。
• 评估康普顿拖曳和同步辐射冷却在辐射效率高系统中对喷流洛伦兹因子的限制作用。
• 识别喷流结构中的普遍特征，如跨音速流和自相似行为，以供未来建模使用。

提出的方法

• 数值模拟求解包含正负电子对产生有效模型及中子Fickian扩散的广义相对论磁流体动力学（GRMHD）方程，以模拟重子污染。
• 模拟采用Kerr度规时空，并求解理想MHD方程，电磁场的应力-能量张量表示为$ T^{ m EM}_{ ueta} = \frac{b^2}{2}(u_\nu u_\beta + P_{\nu\beta}) - b_\nu b_\beta $，其中$ b^\mu $为磁场四维矢量。
• 模型引入GRBs中的“磁火球”情景，即在喷流区发生正负电子对产生，其过程受限于中子扩散，从而确定终端洛伦兹因子。
• 通过阿尔文面、快磁面和驻点面等特征面分析喷流动力学，仅当阿尔文点位于能层内时，能量提取才可能发生。
• 模拟追踪喷流从内吸积盘到大半径区域的演化过程，识别外快磁面之后的加速区域。
• 对喷流结构应用曲线拟合，提取分段自相似特征，以供未来喷流动力学建模使用。

实验结果

研究问题

• RQ1 什么物理机制决定了黑洞吸积系统中相对论性喷流的终端洛伦兹因子？
• RQ2 正负电子对产生和中子扩散如何影响GRBs中Poynting主导喷流的组分及最大洛伦兹因子？
• RQ3 环向磁场不稳定性在半径大于100倍引力半径区域的相对论性喷流时间变异性中起什么作用？
• RQ4 康普顿拖曳和同步辐射冷却如何限制X射线双星和AGN等辐射效率高系统中的洛伦兹因子？
• RQ5 GRBs、AGN和X射线双星中的喷流结构在多大程度上表现出自相似或普遍特征？

主要发现

• 在塌缩星模型中，Poynting主导喷流达到终端洛伦兹因子Γ ∼ 100–1000，与短时硬伽马射线暴的观测一致。
• 喷流加速发生在广阔的径向范围，而非靠近黑洞，外快磁面标志着主要加速区域。
• 在半径大于∼100倍引力半径处，喷流的时间变异性主要由环向磁场不稳定性主导。
• Poynting-重子喷流形成准直喷流，终端洛伦兹因子为Γ ∼ 1.5，与超新星驱动喷流的观测一致。
• 喷流结构表现出分段自相似性，已提供曲线拟合结果，可用于未来喷流动力学建模。
• 该模型证实，康普顿拖曳可在辐射效率高的系统中破坏Poynting主导喷流，使Γ ≲ 2，而同步辐射冷却将AGN喷流的洛伦兹因子限制在Γ ∼ 2–10之间。