[论文解读] Papaloizou-Pringle instability suppression by the magnetorotational instability in relativistic accretion discs

Matteo Bugli, Jérôme Guilet+4 more|arXiv (Cornell University)|Jul 6, 2017

Astrophysical Phenomena and Observations参考文献 79被引用 42

一句话总结

本研究通过3D广义相对论磁流体动力学（GRMHD）模拟，探究了相对论性、磁化吸积环在广义相对论框架下，Papaloizou-Pringle不稳定性（PPI）与磁旋转不稳定性（MRI）之间的相互作用。结果表明，即使存在微弱的环向磁场，MRI驱动的湍流也会抑制主导的m=1 PPI模态，将能量重新分配至更小尺度，从而抑制大尺度不稳定性。该发现对致密天体并合与潮汐瓦解事件中的角动量输运具有重要意义。

ABSTRACT

Geometrically thick tori with constant specific angular momentum have been widely used in the last decades to construct numerical models of accretion flows onto black holes. Such discs are prone to a global non-axisymmetric hydrodynamic instability, known as Papaloizou-Pringle instability (PPI), which can redistribute angular momentum and also lead to an emission of gravitational waves. It is, however, not clear yet how the development of the PPI is affected by the presence of a magnetic field and by the concurrent development of the magnetorotational instability (MRI). We present a numerical analysis using three-dimensional GRMHD simulations of the interplay between the PPI and the MRI considering, for the first time, an analytical magnetized equilibrium solution as initial condition. In the purely hydrodynamic case, the PPI selects as expected the large-scale $m=1$ azimuthal mode as the fastest growing and non-linearly dominant mode. However, when the torus is threaded by a weak toroidal magnetic field, the development of the MRI leads to the suppression of large-scale modes and redistributes power across smaller scales. If the system starts with a significantly excited $m=1$ mode, the PPI can be dominant in a transient phase, before being ultimately quenched by the MRI. Such dynamics may well be important in compact star mergers and tidal disruption events.

研究动机与目标

• 理解磁旋转不稳定性（MRI）如何影响相对论性、几何厚吸积环中Papaloizou-Pringle不稳定性（PPI）的发展。
• 解决一个开放性问题：磁场是否能够抑制或改变PPI，否则PPI将导致低效的大尺度角动量输运。
• 首次利用解析磁化平衡解作为初始条件，开展数值分析，以实现对MRI-PPI相互作用的受控研究。
• 评估MRI诱导湍流对PPI非线性演化及饱和状态的影响，特别是在黑洞吸积和潮汐瓦解事件等天体物理相关情景中的意义。

提出的方法

• 使用ECHO代码进行3D广义相对论磁流体动力学（GRMHD）模拟，求解Kerr-Schild坐标系下的理想MHD方程。
• 采用解析磁化平衡解，构建具有恒定比角动量、并被微弱环向磁场贯穿的几何厚环结构。
• 采用高分辨率有限体积格式，并结合约束传输法以保持磁通量守恒，确保数值稳定性。
• 实现一种保守算法以恢复原始变量，包括熵守恒，防止反演过程中出现负压。
• 通过方位角模数m的傅里叶分解追踪非轴对称模态的增长与演化，重点关注m=1及更高阶谐波。
• 利用适当的雅可比变换，在Boyer-Lindquist与Kerr-Schild坐标系之间转换流体和电磁量，以保持观测者参考系中测量的一致性。

实验结果

研究问题

• RQ1 在相对论性吸积环中，微弱环向磁场的存在如何影响Papaloizou-Pringle不稳定性（PPI）的增长与饱和？
• RQ2 磁旋转不稳定性（MRI）是否抑制或改变PPI的主导m=1模态？如果是，其物理机制是什么？
• RQ3 MRI驱动的湍流在不同方位模态间如何重新分配能量？是否导致从大尺度动力学向小尺度动力学的转变？
• RQ4 在MRI抑制之前，是否可观测到瞬态的PPI活动？这种现象在何种初始条件下发生？
• RQ5 这些结果如何影响我们对致密天体系统（如黑洞并合与潮汐瓦解事件）中角动量输运与吸积效率的理解？

主要发现

• 在纯流体动力学情况下，PPI选择m=1方位模态作为增长最快且非线性主导的不稳定性，与先前的理论预期一致。
• 当引入微弱环向磁场后，MRI迅速发展并抑制大尺度m=1模态，能量被重新分配至更小的方位模数。
• MRI诱导的湍流有效抑制了PPI，即使m=1模态初始被强烈激发，也会导致一个短暂的过渡阶段后完全被抑制。
• MRI导致能量向更小尺度重新分配，表明系统从相干的大振幅螺旋波向更各向同性、类似湍流的能量分布转变。
• MRI对PPI的抑制意味着磁场可能在接近恒定角动量的环中阻止低效的波驱动吸积，转而促进更高效的MRI驱动吸积。
• 这些结果表明，MRI可能在真实磁化吸积流中主导角动量输运，即使PPI原本可能活跃，对潮汐瓦解事件和致密双星并合中的吸积建模具有重要启示。