Universal behaviour of $\alpha$-viscosity in black hole accretion discs

This paper proposes a new formula for the $\alpha$-viscosity coefficient in black hole accretion discs, derived from GRMHD simulations and expressed in terms of the "gyration radius," to accurately capture its universal behavior and improve analytic disc models.

The Gist

Here is an explanation of the paper, translated into everyday language with some creative analogies.

The Big Picture: Fixing the "Traffic Rules" of Black Holes

Imagine a black hole as a giant, cosmic whirlpool. Matter (gas and dust) swirls around it like water going down a drain, forming a spinning disk called an accretion disk. As this matter spirals inward, it gets super hot and glows brightly, which is how we see black holes.

For decades, scientists used a simple rule to predict how this disk behaves. They assumed that the "friction" inside the disk (which causes the matter to lose energy and fall inward) was constant everywhere. Think of it like driving a car and assuming the friction of the tires on the road is exactly the same whether you are on a highway or in a tight parking lot.

This paper says: "That assumption is wrong."

Using powerful supercomputer simulations, the authors discovered that the friction inside a black hole's disk changes dramatically depending on how close you are to the black hole. They have proposed a new, more accurate "rulebook" for how this friction works.

---

The Problem: The Old Rule Was Too Simple

In the 1970s, two scientists named Shakura and Sunyaev created a famous model. They introduced a number called $\alpha$ (alpha) to represent the "viscosity" or internal friction of the disk.

• The Old Idea: They assumed $\alpha$ was a constant number, like a fixed setting on a thermostat.
• The Reality: Modern supercomputer simulations (called GRMHD) show that $\alpha$ is not constant. It's more like a dimmer switch that changes brightness based on your location.
  • Far away: The friction is low.
  • Very close to the black hole: The friction spikes up.
  • Right at the edge (the Horizon): The friction drops to zero.

The Discovery: A Universal Pattern

The authors looked at data from several different teams of scientists who ran these complex simulations. Despite using different computers and methods, they all saw the same pattern. It's as if they were looking at different cars, but they all had the same engine behavior.

They noticed three key features:

1. The Zero Point: At the very edge of the black hole (the Event Horizon), the stress/friction is zero. Nothing can push back against the black hole once it crosses that line.
2. The Peak: Just outside the horizon, near the path where light itself would orbit the black hole (the "photon orbit"), the friction hits its maximum. It's like a traffic jam forming right before the exit ramp.
3. The Tail: Far away from the black hole, the friction settles down to a low, steady value.

The Solution: A New Formula

The authors propose a new mathematical formula to describe this behavior. Instead of a single number, they use a formula that changes based on your distance from the black hole.

To explain this, they use a concept called the "Gyration Radius."

• The Analogy: Imagine a figure skater spinning. If they pull their arms in, they spin faster. The "gyration radius" is a measure of how spread out their mass is.
• In the Paper: The authors found that the friction ($\alpha$) is directly linked to this "gyration radius." As you get closer to the black hole, the geometry of space-time changes this radius, which in turn changes the friction.

Their new formula is like a smart thermostat that automatically adjusts the friction based on how close you are to the black hole, rather than keeping it fixed.

Why Does This Matter?

You might ask, "Why do we need a new formula? The old one worked fine."

1. It's More Realistic: The old model was a great simplification, but it breaks down when the black hole is eating a lot of food (high accretion rates). The new formula accounts for the weird physics of Einstein's General Relativity.
2. Better Predictions: By using this new rule, scientists can build better models of what these disks actually look like. This helps astronomers interpret the light we see from black holes in the real universe.
3. Understanding the "Plunge": The area right before the black hole (the "plunging region") is where the physics gets weird. The old models couldn't explain why the friction behaves the way it does there. This new formula captures that behavior perfectly.

The "Why" Behind the Magic

The paper also hints at why this happens, though they promise a deeper math explanation later:

• Why Zero at the Horizon? Once matter crosses the event horizon, it can't go back. There is no "upstream" flow to create friction against. It's like a waterfall; once you go over the edge, you can't push back against the water falling down.
• Why the Peak? The peak happens near the "photon orbit." In this zone, the rules of rotation flip. Usually, things spin faster as they get closer to the center. But near a black hole, the geometry is so twisted that stability requires things to behave differently, causing a massive buildup of stress (friction) right there.

Summary

Think of this paper as updating the GPS navigation system for black holes.

• Old GPS: "Drive at a steady speed." (Too simple, leads to wrong turns).
• New GPS: "Slow down here, speed up there, and stop completely at the edge." (Accurate, based on real-time data).

The authors have found a universal "traffic law" for black hole disks that works for everyone, from the smallest stellar black holes to the giants at the center of galaxies. This will help us understand the universe a little bit better.

Technical Summary

Here is a detailed technical summary of the paper "Universal behaviour of $\alpha$-viscosity in black hole accretion discs" by Abramowicz et al.

1. Problem Statement

The classic Shakura-Sunyaev (1973) model of accretion discs relies on the ansatz that the viscous stress $T$ is proportional to the total pressure $p$ via a constant coefficient, $\alpha = T/p = \text{const}$. While this assumption allowed for the development of foundational analytic models, modern General Relativistic Magnetohydrodynamic (GRMHD) simulations have revealed that $\alpha$ is not constant.

• The Discrepancy: Simulations show that $\alpha$ varies by approximately an order of magnitude within the disc. It is significantly smaller at large radii and increases dramatically in the "plunging region" (between the sonic radius and the event horizon).
• The Gap: Current analytic models (including "slim disc" models) still rely on the constant $\alpha$ assumption or lack a realistic prescription for the stress in the plunging region. This limits the accuracy of calculating radiation flux and the observational appearance of black hole accretion, particularly for high-luminosity discs where radial heat advection is significant.
• The Goal: To derive a universal, analytic formula for $\alpha(r)$ that accurately reflects the behavior observed in GRMHD simulations, specifically capturing the stress's behavior near the event horizon and the photon orbit.

2. Methodology

The authors employed a comparative analysis of existing GRMHD simulation data to identify universal patterns and then constructed a semi-analytic model based on those patterns.

• Data Sources: The study analyzed results from three independent GRMHD simulation sets:
  • L19 (Lančová et al. 2019): Radiative GRMHD simulations of a "puffy" disc around a non-rotating (Schwarzschild) black hole. This is the most physically realistic dataset, including radiation pressure and energy exchange.
  • P13 (Penna et al. 2013): Non-radiative GRMHD simulations.
  • R25 (Rule et al. 2025): Non-radiative simulations focusing on Maxwell stress in the local fluid frame.
• Stress Calculation: The authors examined the total stress $T$ (sum of Reynolds and Maxwell components) calculated in the comoving tetrad frame. They noted that while simulation methods differ (e.g., inclusion of radiation pressure), the radial profile of $\alpha$ exhibits a consistent universal shape.
• Theoretical Framework: The authors utilized the concept of the gyration radius ($\tilde{r}$), a geometric quantity defined in the Kerr metric by the ratio of the Killing vectors for axial symmetry ($\xi^\mu$) and time symmetry ($\eta^\mu$):
  $$\tilde{r} = \sqrt{\frac{-(\xi\xi)}{(\eta\eta)}}$$
  In the Schwarzschild limit, this relates the specific angular momentum $j$ and angular velocity $\Omega$ via $j = \Omega \tilde{r}^2$.

3. Key Contributions

The primary contribution is the proposal of a new, universal $\alpha$-prescription formula that replaces the constant $\alpha$ assumption.

The Proposed Formula:
The authors propose a coordinate-independent formula for $\alpha$ that depends on the gyration radius and the black hole horizon radius ($r_H$):
$$\alpha = \left[ \alpha_P \left( \frac{r_H}{\tilde{r}} \right)^2 + \alpha_\infty \right] (1 - \frac{r_H}{r})$$
Or, expressed using Killing vectors for the Schwarzschild case:
$$\alpha = \alpha_P \left( \frac{r_H}{\tilde{r}} \right)^2 + \alpha_\infty (\eta\eta)$$

Parameters:

• $\alpha_P \approx 4.71$: A phenomenological constant dominating near the black hole.
• $\alpha_\infty \approx 0.01$: A constant representing the asymptotic value at large radii.
• The term $(1 - r_H/r)$ ensures the stress vanishes at the event horizon.

Key Physical Insights:

1. Zero Stress at Horizon: The formula enforces $T=0$ at $r=r_H$, consistent with the physical argument that no upstream mass flux exists across the horizon to generate turbulent stress.
2. Peak at Photon Orbit: The formula naturally produces a maximum in $\alpha$ near the circular photon orbit ($r_{ph} = 1.5 r_H$). The authors hypothesize this is linked to the reversal of the radial direction's dynamical meaning at the photon sphere (where "outward" points toward the black hole).
3. Asymptotic Behavior: At large radii ($r \gg r_H$), $\alpha$ approaches a constant value related to the ratio of shear ($\sigma$) to vorticity ($\omega$), consistent with Magnetorotational Instability (MRI) turbulence theory.

4. Results

• Universal Fit: The proposed formula provides an excellent fit to the $\alpha$ profiles derived from all three simulation sets (L19, P13, R25), despite their differences in physics (radiative vs. non-radiative) and numerical methods.
• Stress Profile Confirmation: The simulations confirm three distinct features:
  1. $\alpha = 0$ at the horizon.
  2. A sharp peak in $\alpha$ near the photon orbit ($r_{ph}$).
  3. A significant drop in $\alpha$ in the outer disc regions compared to the plunging region.
• Validation of "Puffy" Discs: The L19 radiative simulations show that even in geometrically thick, radiation-pressure-dominated discs, this universal stress behavior holds. The "puffy" structure (dense core + puffy region + funnel) does not invalidate the underlying stress mechanism.
• Analytic Utility: The formula is short, explicit, and expressed in terms of coordinate-independent geometric quantities, making it immediately applicable to analytic models.

5. Significance

• Improving Analytic Models: This work bridges the gap between complex numerical simulations and tractable analytic models. By replacing the constant $\alpha$ with the proposed function, researchers can construct improved "thin" and "slim" disc models that accurately account for stress in the plunging region.
• Observational Accuracy: Accurate modeling of the stress in the innermost regions is crucial for calculating the radiation flux $F(r)$ and the spectral energy distribution of black hole accretion discs. This is particularly important for high-accretion-rate systems (near Eddington limits) where the plunging region contributes significantly to the emission.
• Fundamental Physics: The paper suggests that the stress behavior is not merely a result of local turbulence but is fundamentally constrained by the global geometry of the black hole spacetime (specifically the horizon and photon orbit). This points toward a deeper connection between General Relativity and magnetohydrodynamic turbulence.
• Future Directions: The authors plan to use this prescription to develop stationary and non-stationary slim disc models, potentially explaining the long-term evolution and observational appearance of black hole systems more accurately than current models allow.