🔭 astrophysics

Radiative cooling effects on black hole hot accretion flows around the sub-Bondi radius

This paper presents two-dimensional magnetohydrodynamic simulations demonstrating that increasing radiative cooling in hot accretion flows around the sub-Bondi radius thins the accretion disk, marginally stabilizes it against convection, and significantly suppresses wind-driven mass loss by shifting the primary mechanism for inward mass inflow reduction from winds to turbulence driven by magnetorotational instability and convection.

Original authors: Mu-Qing Liu, Xiao-Hong Yang, De-Fu Bu

Published 2026-04-16

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Mu-Qing Liu, Xiao-Hong Yang, De-Fu Bu

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a supermassive black hole as a giant, cosmic vacuum cleaner sitting at the center of a galaxy. Usually, we think of it just sucking everything in, but in reality, the gas swirling around it is a chaotic, swirling mess of hot plasma, magnetic fields, and winds.

This paper is like a high-speed, 3D weather simulation, but instead of clouds and rain, it's modeling the "atmosphere" around a black hole. The scientists wanted to answer a specific question: What happens when this hot gas gets hot enough to glow and lose energy (cool down) as it falls in?

Here is the story of their findings, broken down with some everyday analogies.

The Setup: The Cosmic Bathtub

Think of the black hole as the drain in a giant bathtub. The water (gas) is swirling around the drain.

The Bondi Radius: This is the edge of the bathtub where the water starts to feel the pull of the drain.
The Problem: It's hard to simulate the whole bathtub and the drain at the same time because the scale is so huge. Most previous studies only looked at the water right near the drain. This team decided to zoom out and look at the whole "sub-Bondi" area (the region just outside the immediate drain).
The Variable: They ran simulations with different amounts of gas. More gas means more friction and heat, which leads to more "radiative cooling" (the gas glowing and losing energy, like a hot pan cooling down on a stove).

The Big Discovery: The "Shrink Wrap" Effect

When the gas is thin and hot (low accretion rate), it behaves like a puffy, fluffy cloud. It's thick and puffed up.

But as they added more gas (increasing the accretion rate), the gas started to glow brighter and cool down faster.

The Analogy: Imagine a puffy marshmallow. If you put it in a freezer (radiative cooling), it loses its air and shrinks.
The Result: The "hot accretion flow" didn't stay puffy. It got squeezed flat. The disk of gas became much thinner and denser, like a flat pancake instead of a fluffy cloud.

The Wind vs. The Turbulence

One of the biggest mysteries in black hole physics is: Why doesn't all the gas just fall in?
Usually, the answer is "winds." The gas gets so hot and energetic that it blows away, like steam escaping a kettle. This prevents the black hole from eating everything.

The "Weak Cooling" Scenario: When the gas is puffy and hot, strong winds blow out. These winds are the main reason gas doesn't fall in. It's like a strong gale pushing the water away from the drain.
The "Strong Cooling" Scenario: When the gas cools down (due to high density), something surprising happens. The winds die down. The gas loses the energy needed to blow away.
The New Culprit: If the winds stop blowing, why doesn't the gas fall in? The scientists found that turbulence takes over. Think of it like a whirlpool in a sink. Even if there's no wind blowing the water away, the swirling, chaotic motion (turbulence) keeps the water from sliding smoothly down the drain. The gas gets stuck in a chaotic dance of magnetic fields and swirling eddies.

The Magnetic "Rubber Bands"

The black hole is surrounded by invisible magnetic fields, which act like tangled rubber bands.

When the gas cools and shrinks (the "Shrink Wrap" effect), it squeezes these magnetic rubber bands tighter.
This makes the magnetic pressure inside the disk much stronger.
However, this strong magnetic pressure actually makes it harder for the gas to move inward efficiently. It's like trying to push a crowd through a doorway that is being squeezed shut by elastic bands. The gas gets stuck, and the "accretion" (the eating process) slows down significantly.

The "Höiland" Stability Check

The scientists also checked if the gas was stable or if it wanted to churn up like boiling water.

They used a rule called the Höiland criterion (a fancy way of checking if a fluid is stable).
They found that the gas is "marginally stable." It's not perfectly calm, but it's not boiling over either. It's in a delicate balance, about 55–60% stable, meaning it's just barely holding it together.

The Bottom Line: Cooling Slows the Feast

The most counter-intuitive finding is this: Cooling the gas actually stops the black hole from eating.

In previous studies using simpler models, scientists thought cooling would help gas fall in faster. But this study, which included complex magnetic fields, showed the opposite:

Cooling makes the gas shrink.
Shrinking makes the magnetic fields stronger.
Strong magnetic fields create turbulence that blocks the flow.
Result: The black hole gets less food, and the powerful winds that usually blow the gas away disappear.

In a nutshell: When the cosmic gas around a black hole gets too hot to handle and starts to cool down, it doesn't fall in faster. Instead, it gets squeezed flat, tangled up in magnetic knots, and the whole system slows down, leaving the black hole with a much smaller meal than expected.

1. Problem Statement

The accretion processes around supermassive black holes (SMBHs) in low-luminosity active galactic nuclei (LLAGNs) are generally modeled as hot, radiatively inefficient accretion flows (ADAFs). While previous numerical simulations have extensively studied the region close to the black hole ( $r < 10^3 r_s$ ), they have largely neglected the dynamics at larger scales extending to the Bondi radius ( $r_B$ ).

The Gap: It is computationally challenging to simulate the vast dynamical range from the Bondi radius to the event horizon simultaneously. Consequently, the behavior of hot accretion flows at the "sub-Bondi" scale ( $r < r_B$ ) and the specific impact of radiative cooling on wind production and accretion efficiency in this region remain unclear.
Key Questions: How does radiative cooling affect the structure, stability, and wind production of hot accretion flows at the sub-Bondi radius? Does strong cooling promote or suppress accretion and winds in the presence of magnetic fields?

2. Methodology

The authors performed two-dimensional (2D) axisymmetric Magnetohydrodynamic (MHD) simulations using the PLUTO code.

Physical Model:
- Equations: Solved the full set of MHD equations including mass, momentum, energy, and induction equations.
- Gravity: Used the pseudo-Newtonian potential (Paczyński–Wiita) to approximate the Schwarzschild black hole gravity.
- Cooling: Included bremsstrahlung cooling ( $Q^-_{brem}$ ) as a source term in the energy equation. The cooling rate scales with density squared ( $n^2$ ) and temperature.
- Heating: Heating was calculated based on the work done by Maxwell and Reynolds stresses (turbulent dissipation), consistent with conservative numerical schemes.
- Equation of State: Ideal gas with adiabatic index $\gamma = 5/3$ .
Setup:
- Domain: Radial range from $2 \times 10^{-2} r_B$ to $5 r_B$ ; Angular range $0$ to $\pi$ .
- Initial Conditions: A rotating torus with constant specific angular momentum embedded in a low-density ambient medium.
- Magnetic Field: A weak, purely poloidal seed field was introduced. The topology was designed to allow the Magnetorotational Instability (MRI) to grow. The initial plasma $\beta$ (ratio of gas to magnetic pressure) was set to 1000.
- Resolution: High-resolution grids (256 radial zones, ~20,000 angular zones) ensuring the MRI quality factors ( $Q_r, Q_\theta$ ) exceeded 10, guaranteeing adequate resolution of turbulent modes.
Simulation Runs:
- Four main runs were conducted by varying the initial density of the torus ( $\rho_0$ ) to simulate different accretion rates.
- R22, R23, R24: Runs with radiative cooling ON (increasing density from $10^{-24}$ to $10^{-22}$ g cm $^{-3}$ ).
- N24: A control run with the same density as R24 but radiative cooling OFF.

3. Key Contributions

Scale Bridging: Successfully simulated the dynamics of hot accretion flows specifically in the sub-Bondi radius region, a scale previously under-explored in MHD studies.
MHD vs. HD Comparison: Unlike previous Hydrodynamic (HD) studies that used $\alpha$ -viscosity prescriptions, this study uses self-consistent MHD turbulence (MRI) to drive angular momentum transport, providing a more realistic physical picture.
Cooling Mechanism Analysis: Detailed the transition from advection-dominated flows to cooling-dominated flows and quantified the non-linear feedback of cooling on magnetic field amplification and wind suppression.

4. Key Results

A. Disk Structure and Contraction

Contraction: As the accretion rate increases (from R24 to R22), radiative cooling strengthens, causing the accretion flow to contract vertically.
Scale Height: The ratio of disk scale height to radius ( $H/r$ ) decreases significantly in the strong cooling run (R22), dropping below 0.2 at $r > 0.1 r_B$ , compared to 0.3–0.6 in weak cooling runs.
Density Profile: Strong cooling leads to a flatter density profile in the inner region ( $\rho \propto r^{-p}$ with $p \approx 0.05$ for R22 vs. $p \approx 0.37$ for R24).

B. Magnetic Field Amplification

Mechanism: In the strong cooling run (R22), the magnetic field strength near the equatorial plane increases significantly.
Cause: This is primarily driven by vertical contraction of the gas due to cooling (compressing field lines) rather than magnetic flux advection by converging flows.
Dominance: The toroidal magnetic field dominates the poloidal field in most regions, accounting for 60–90% of the magnetic energy.

C. Convective Stability

Høiland Criterion: Analysis using the Høiland criterion reveals that the accretion flows are marginally convectively stable.
Stability Fraction: Approximately 55–62% of the area within the density scale height is convectively stable. This contrasts with some BH-proximal flows which are often found to be fully stable, suggesting the sub-Bondi region promotes convective instability.

D. Wind Production and Mass Outflow

Wind Definition: Winds are defined as outflows with a positive Bernoulli parameter ($Be > 0$).
Cooling Suppression:
- Weak Cooling (R23, R24, N24): Winds constitute 50–70% of the total mass outflow rate and are the primary driver for the inward decrease of the mass inflow rate.
- Strong Cooling (R22): Radiative cooling significantly suppresses wind production. The wind contribution drops below 50% (and <15% outside the disk at large radii).
Driving Mechanism: The winds are driven by a "micro" Blandford–Payne mechanism, where centrifugal force, combined with gas and magnetic pressure gradients, accelerates matter. In strong cooling runs, the net outward force is significantly weaker.

E. Accretion Efficiency

Paradoxical Suppression: Contrary to some HD simulations (which suggested cooling promotes accretion by reducing outflows), the MHD results show that strong radiative cooling suppresses the net accretion rate.
Reason: Strong cooling reduces the gas temperature, which in turn reduces the specific angular momentum flux driven by Maxwell stresses. This leads to a significant drop in the mass inflow rate (by nearly an order of magnitude in R22 compared to the non-cooling case).
Viscosity: The effective viscosity parameter ( $\alpha$ ) drops significantly in the strong cooling run at $r > 0.08 r_B$ , inhibiting angular momentum transport.

5. Significance and Implications

Revising Accretion Models: The study challenges the notion that radiative cooling always enhances accretion in hot flows. In MHD contexts, strong cooling can inhibit accretion by weakening the turbulent stresses required for angular momentum transport.
Wind Feedback: The results suggest that in systems with high accretion rates (where cooling is strong), the feedback from winds may be significantly weaker than predicted by non-radiative or HD models. This has implications for AGN feedback mechanisms in galaxy evolution.
Observational Context: The findings help interpret the properties of LLAGNs (like Sgr A*), suggesting that the interplay between cooling, magnetic fields, and turbulence at the Bondi scale is critical for determining the actual mass supply to the black hole and the strength of outflows.

In summary, this paper demonstrates that radiative cooling acts as a double-edged sword: it compresses the disk and strengthens magnetic fields locally, but ultimately suppresses the global accretion rate and wind production by reducing the efficiency of MRI-driven angular momentum transport.