Idealized Global Models of Accretion Disks with Strong Toroidal Magnetic Fields

技术摘要：具有强环向磁场的理想化全球吸积盘模型

问题陈述
关于吸力盘中角动量输运的本质，特别是对于在没有净垂直磁通量的情况下由强环向磁场主导的吸积盘，目前仍存在争议。虽然零净垂直磁通量（ZNVF）的局部剪切盒模拟已经确定了一个由局部发电机维持的“低 $\beta$ 态”（其中等离子体 $\beta \sim 1$），但全球模拟在历史上产生的磁盘要么是高度磁化的（$\beta \ll 1$），要么因为浮力损失（帕克不稳定性）而无法维持强环向磁场。先前的全球研究（例如 Hopkins 等人 2, 2024a,b；Guo 等人 2024）经常发现 $\beta \ll 1$ 且 $\alpha \gg 1$，而对矮新星和 X 射线瞬变源的观测推断则表明 $\alpha \sim 0.1 - 0.3$，对应于中度磁化的磁盘。局部与全球结果之间定量差异的起源，以及维持 $\beta \sim 1$ 环向主导磁盘所需的条件，尚未得到充分理解。

方法论
作者展示了使用 AthenaK 代码进行的理想化 3D 全球磁流体力学（MHD）模拟。该设置涉及一个固定的质点引力势（牛顿力学，无广义相对论）和一个旋转的磁化气体云，其初始条件为环向磁场。气体在初始状态下并不处于平衡态；它发生流入、圆化并形成磁盘。

• 状态方程： 使用了局部等温状态方程，固定了气体的热标高 $H_{th}$ 相对于半径的关系（$H_{th} = \sqrt{2}hR$）。这使得作者能够控制相对于热标高的分辨率。
• 初始条件： 初始气体云在较大半径处具有均匀的密度和角动量。磁场在圆化半径（$R_{circ} = 16 r_0$）之外被初始化为纯环向磁场，在较大半径处等离子体 $\beta_\phi = 1$。部分运行包括用于对比的净垂直磁场。
• 数值分辨率： 模拟利用了高达 9 层的自适应网格细化（AMR）。一个关键参数是每个热标高的单元数（$h r / \Delta x$）。作者通过改变热标高参数 $h$ 和网格分辨率进行参数研究，以测试收敛性。
• 边界条件： 内真空汇（inner vacuum sink）移除了进入 $r < r_0$ 的气体，而外部边界则固定为初始条件。

主要结果

1. 依赖于分辨率的磁性转变： 最显著的发现是，磁盘的磁性状态取决于气体热标高是否被充分解析。

   • 未解析（$h r / \Delta x \lesssim 3$）： 磁盘保持高度磁化，$\beta \ll 1$ 且高吸积应力 $\alpha \gg 1$。磁盘由磁压在垂直方向上支撑。
   • 已解析（$h r / \Delta x \gtrsim 4$）： 磁盘在 $\sim 50-180$ 个轨道周期内发生转变。它从内向外发生垂直塌缩，形成一个中度磁化的中平面，其 $\beta \sim 1$ 且 $\alpha \sim 0.1$。在这种状态下，中平面由气体压力主导（接近热平衡），而中平面上方和下方的区域仍保持强磁化（$\beta \ll 1$）。

2. 磁盘结构与演化：

   • 系统最初形成一个快速吸积、强磁化的开普勒磁盘。
   • 随着磁盘演化，它过渡到一个准稳态，其中中平面具有跨阿尔芬波（trans-Alfvénic）和跨声速（trans-sonic）的湍流波动（$\delta v^2 \sim \delta v_A^2 \sim v_A^2 \sim c_s^2$）。
   • 角动量输运由平均应力和波动应力共同驱动。在中平面内，波动麦克斯韦应力占主导地位，而在较高高度处，平均应力占主导地位。

3. 维持磁场的机制：

   • 方位角磁场由于浮力作用不断从垂直方向逃逸，但被局部发电机所补充。
   • 在已解析的磁盘中（$R < R_{circ}$），发电机由径向磁场的剪切驱动（该磁场由大规模流入产生），并与湍流扩散达到平衡。径向磁场平流在中盘内部是次要的，但在圆化半径之外变得重要，有助于维持外部区域的 $\beta \ll 1$ 状态。
   • 对于已解析情况的全球结果，在很大程度上与局部 ZNVF 剪切盒模拟（Squire 等人 2024）中发现的“低 $\beta$ 态”一致。

4. 鲁棒性测试：

   • $\beta \sim 1$ 状态对于初始磁场强度的变化、初始湍流的引入以及不同重构方法（PLM 与 PPM）的使用是鲁棒的。
   • 包含净垂直磁场会导致不同的吸积模式，其具有更强磁化的中平面（$\beta \lesssim 1$）以及由平均麦克斯韦应力主导的应力，这与之前的 NVF 研究一致。
   • 网格细化测试确认，向 $\beta \sim 1$ 的转变是一个数值分辨率效应：细化一个未充分解析的磁盘会触发向 $\beta \sim 1$ 态的塌缩，而降低分辨率则会返回 $\beta \ll 1$ 态。

意义与主张
本文声称解决了关于是否存在 $\beta \sim 1$ 环向主导磁盘的局部与全球模拟之间的差异。作者认为，先前的全球模拟未能捕捉到这种状态，可能是因为热标高未得到充分解析。

• 观测相关性： 在已解析的 $\beta \sim 1$ 态中，产生的吸积应力 $\alpha \sim 0.1$ 与从矮新星和 X 射线瞬变爆发中推断出的数值一致，这表明此类磁盘可能自然地解释这些现象，而无需要求净垂直磁通量。
• 物理洞察： 该研究证明，只要热结构被充分解析，强环向场可以通过局部发电机在全域环境下对抗浮力损失得以维持。
• 局限性： 作者指出，结果并非严格收敛；局部剪切盒研究表明，在使用更高分辨率（需要 $\sim 10$ 个单元/标高）时，$\alpha$ 和 $\beta$ 可能仍会继续演化。此外，模拟是理想化的（牛顿力学、等温），忽略了辐射传输、自引力和广义相对论效应，这些因素可能会改变现实天体系统中定量的细节。这些结果与其他全球模拟（例如 Hopkins 等人 2024b）之间的差异尚不明确，这可能归因于分辨率或热力学及流入性质的差异。