[論文レビュー] Universal behaviour of $α$-viscosity in black hole accretion discs

M. A. Abramowicz, Axel Brandenburg+5 more|arXiv (Cornell University)|Mar 11, 2026

Astrophysical Phenomena and Observations被引用数 0

ひとこと要約

論文はブラックホール付着円盤のGRMHDに基づく普遍的な α(r) 粘性プロファイルを導出し、地平線でゼロ、光子軌道付近で最大となる性質を利用して薄い円盤およびスリム円盤モデルを改善可能にする。

ABSTRACT

The Shakura-Sunyaev $α$-viscosity coefficient, defined as the ratio of total stress to total pressure, $α= \mathbb{T}/p$, played an important role in the development of the accretion disc theory in the early 1970s. The origin of turbulence that causes the stress $\mathbb{T}$ was unknown at that time. Shakura and Sunyaev assumed $α=$ const. Today we know that this was not quite realistic - the modern general relativistic magneto-hydrodynamic simulations (GRMHD) of black hole accretion discs revealed that $α$ changes by about an order of magnitude within the disc, being smaller far away from the black hole and larger in the plunging region close in. It was found that the behaviour of $α$ reflects some underlying, fundamental properties of the stress $\mathbb{T}$ itself. In particular, as argued by several authors, the stress must be zero at the black hole horizon. We notice that the stress calculated in GRMHD simulations by different authors, including us, has a maximum rather close to the location of the circular photon orbit. We propose a formula that accurately describes this universal behaviour of $α$ in terms of the "gyration radius'', a physical characteristic of rotation well known in Newtonian dynamics and in the black hole case uniquely defined by the Kerr space-time geometry. Analytic and semi-analytic models of black hole accretion discs provide an invaluable insight into fundamental physics, and the GRMHD simulations do not aspire to replace them. Rather, simulations could help to improve analytic models by making them more realistic. For example, our $α$-formula, deduced from the GRMHD simulations, may be handy in the construction of improved versions of thin and slim disc models.

研究の動機と目的

• Shakura–Sunyaev 円盤理論における一定の α を超え、GRMHD シミュレーションで観測される現実的なストレス挙動を取り入れる必要性を動機づける。
• Kerr/シュワルツシルト時空における主要なストレス特性を捉える解析的 α(r) 処方を提案する。
• α-プロファイルを回転対称時空の測地ベクトル幾何学とジラレーション半径（gyration radius）に結びつけ、座標不変の定式化を行う。
• GRMHD 結果（薄・スリム円盤モデルを改善するための半解析モデルへ適用可能なフィット）を提供する。

提案手法

• 提案された α(r) 処方を導入: α = αP (rH/ r̃)^2 + α∞ (ηη), ここで r̃ は gyration radius、αP、α∞ は現象的定数。
• r̃ を Kerr 極座標の Killing ベクトルから表現: r̃^2 = - (ξξ)/(ηη) および j = Ω r̃^2 との関係。
• シュワルツシルト/ Kerr 時空における Killing ベクトル η および ξ を用いた座標不変形を導出。
• GRMHD シミュレーション（L19, P13, R25）からのストレス Tμν を T = <Tμν> e^(μ)_(r) e^(ν)_(φ) = <T_(r)(φ)> として計算し、α = T/p を定義する（p には利用可能な気体・磁気・放射圧を含む）。
• ストレスが地平線でゼロとなり、円形光子軌道付近でピークし、外径で減衰することを示し、MRI ドリブン乱流と大域的な時空制約と整合する。
• Puffy disk（放射GRMHD）および非放射シミュレーションを用いて処方を検証する。

Figure 1: The $\alpha$ -viscosity coefficient calculated in GRMHD simulations by L19 , Penna et al. ( 2013 ) , denoted by P13 , and Rule et al. ( 2025 ) , denoted by R25 , all for non-rotating black hole accretion discs, fitted to our $\alpha$ -prescription formula ( 13 ). Three characteristic featu

実験結果

リサーチクエスチョン

• RQ1 GRMHD シミュレーションにおけるブラックホ hole 付着円盤の α 粘性係数の半径依存性は何か？
• RQ2 なぜイベントホライズンでストレスがゼロになり、円形光子軌道付近でピークとなるのか、そしてこれを単純な解析形にどう取り込むのか？
• RQ3 提案された α(r) 処方はサブエディントン領域からほぼエディントン領域に至る範囲で半解析的薄円盤・スリム円盤モデルを改善できるのか？
• RQ4 gyration radius と時空幾何学がブラックホール周りの円盤のストレス分布をどのように支配するのか？

主な発見

• GRMHD の独立系グループによるシミュレーションは、提案された形式を介して普遍的な半径パターンを示す。
• α は地平線でゼロ、円形光子軌道付近で最大となる（ r_ph ≈ (3/2) r_H 付近）。
• 大きな半径で α は plunging 領域に比べて非常に小さく、MRI駆動乱流および重力場全体の制約と整合する。
• αP ≈ 4.71 および α∞ ≈ 0.01 のフィット型は、ケース全体でシミュレーション結果を正確に記述する。
• 座標不変の定式化は α を gyration radius および Killing ベクトルに結びつけ、座標選択に依存しない頑健な記述を提供する。
• この枠組みは α-式を用いて、現実的なストレス挙動を組み込んだ薄円盤およびスリム円盤モデルを構築するのに適用可能である。

Figure 2: The structure of a puffy disc consists of three regions: C: dense core, P: puffy region, and F: funnel. Left: The colour map of gas density with the magnetic field lines. Right : Viscous $\alpha$ colour map, and the magnetosonic surface (yellow contour). In both panels, the white dashed li