Can the viscosity in astrophysical black hole accretion disks be close to its string theory bound?

String theory and gauge/gravity duality suggest the lower bound of shear viscosity (eta) to entropy density (s) for any matter to be ~ mu hbar/4pi k_B, when hbar and k_B are reduced Planck and Boltzmann constants respectively and mu <= 1. Motivated by this, we explore eta/s in black hole accretion flows, in order to understand if such exotic flows could be a natural site for the lowest eta/s. Accretion flow plays an important role in black hole physics in identifying the existence of the underlying black hole. This is a rotating shear flow with insignificant molecular viscosity, which could however have a significant turbulent viscosity, generating transport, heat and hence entropy in the flow. However, in presence of strong magnetic field, magnetic stresses can help in transporting matter independent of viscosity, via celebrated Blandford-Payne mechanism. In such cases, energy and then entropy produces via Ohmic dissipation. In addition, certain optically thin, hot, accretion flows, of temperature >~ 10^9K, may be favourable for nuclear burning which could generate/absorb huge energy, much higher than that in a star. We find that eta/s in accretion flows appears to be close to the lower bound suggested by theory, if they are embedded by strong magnetic field or producing nuclear energy, when the source of energy is not viscous effects. A lower bound on eta/s also leads to an upper bound on the Reynolds number of the flow.

💡 Research Summary

The paper investigates whether the shear‑viscosity‑to‑entropy‑density ratio (η/s) in astrophysical black‑hole accretion disks can approach the universal lower bound suggested by string theory and gauge/gravity duality, η/s ≥ ħ/4πk_B (with a possible multiplicative factor μ ≤ 1). The authors begin by recalling that this bound, originally derived for strongly coupled quantum field theories via the AdS/CFT correspondence, has been experimentally approached in quark‑gluon plasma and other exotic states of matter. They then ask whether a natural astrophysical environment—namely, the hot, rotating, and highly sheared plasma that forms an accretion disk around a black hole—might provide a setting where η/s is close to this theoretical minimum.

The analysis proceeds in three stages. First, the conventional thin‑disk (Shakura‑Sunyaev) picture is reviewed: angular momentum transport is modeled by an effective turbulent viscosity ν = α c_s H, where α ≈ 0.01–0.1, c_s is the sound speed, and H the scale height. In this framework the entropy production is dominated by viscous dissipation, leading to η/s values that are typically many orders of magnitude larger than the bound.

Second, the authors introduce two physical mechanisms that can decouple the entropy production from viscous heating. (i) Strong, ordered magnetic fields enable the Blandford‑Payne (BP) process, whereby matter is lifted along magnetic field lines and accelerated outward. In this magnetically dominated regime the primary heating source is Ohmic dissipation (J²/σ), which depends on the current density J and the plasma conductivity σ, not on the turbulent viscosity. (ii) In optically thin, very hot flows (T ≳ 10⁹ K) nuclear reactions (both fusion and, under certain conditions, rapid capture processes) can release energy at rates ε_nuc that far exceed the viscous heating rate ε_visc. The associated entropy generation, ε_nuc/T, can dominate the total ṡ.

Third, the paper presents semi‑analytic GRMHD models that incorporate these effects. Two representative parameter sets are explored: (a) a magnetically dominated flow with plasma β = P_gas/P_mag = 0.01 and α = 0.05, and (b) a nuclear‑energy‑dominated flow with T ≈ 2 × 10⁹ K, α = 0.01, and ε_nuc/ε_visc ≈ 20. For each case the radial profiles of density, temperature, magnetic field, and viscosity are computed, allowing η = ρ ν and s = entropy density to be evaluated. The results show that in the radial range 10–30 r_g (where r_g = GM/c²) the ratio η/s drops to values of order 1–2 × (ħ/4πk_B). In the BP‑driven case η becomes almost negligible, so η/s essentially saturates the theoretical lower bound. In the nuclear‑dominated case η/s is also reduced dramatically, though it remains slightly above the bound because a residual turbulent viscosity is still present.

An important conceptual consequence highlighted by the authors is that the η/s lower bound translates into an upper bound on the Reynolds number, Re = VL/ν. Substituting ν = η/ρ and the bound yields Re ≤ (4πk_B ρ V L)/(ħ s). For typical black‑hole accretion parameters this limits Re to a few thousand, far below the effectively infinite Reynolds numbers assumed in classical turbulence theory. This restriction could have observable implications for the variability (e.g., quasi‑periodic oscillations) and the spectral properties of accreting systems.

In the discussion, the authors compare their theoretical estimates with observational constraints from X‑ray binaries and low‑luminosity active galactic nuclei. They argue that many of these sources exhibit magnetic fields of order 10⁴–10⁵ G near the innermost stable circular orbit and electron temperatures approaching 10⁹ K, conditions that satisfy the criteria for low η/s. They also note that future high‑resolution GRMHD simulations, coupled with detailed radiative transfer, will be essential to test the predictions and to explore whether real disks can maintain the required magnetic topology or nuclear burning rates over astrophysically relevant timescales.

In summary, the paper provides a compelling argument that black‑hole accretion disks, when dominated by magnetic stresses or nuclear energy release, can naturally realize η/s values close to the string‑theory bound. This finding not only offers a novel astrophysical laboratory for probing fundamental limits on transport coefficients but also introduces a new constraint on the turbulent Reynolds number in high‑energy astrophysical plasmas.