CDAF: Convective-Dominated Accretion Flow

• CDAF is a radiatively inefficient accretion flow dominated by large-scale convection that traps gas in quasi-stationary eddies.
• The model exhibits a flat density profile (ρ ∝ r⁻¹/2) with mass inflow linearly decreasing with radius, leading to a suppressed net accretion rate.
• While pure CDAF dynamics occur in non-magnetized flows, realistic MHD simulations suggest ADIOS-type outflows often prevail.

A convection-dominated accretion flow (CDAF) is a class of radiatively inefficient, hot accretion solution in which large-scale convective motions govern both energy and angular-momentum transport. In this regime, convection dominates over both advective inflow and outflow, trapping most of the accreting gas in quasi-stationary eddies and substantially suppressing the net inward accretion rate at small radii. This theoretical framework was originally developed in the context of black hole accretion at low accretion rates, particularly for systems with low Bondi mass supply and minimal radiative cooling. CDAFs produce accretion structures that are geometrically thick, optically thin, and inefficiently radiating, with inflow rates falling sharply toward the black hole, explaining the observed low luminosities of supermassive black holes such as Sgr A*, M31*, and M87* (Inayoshi et al., 2017, Yuan et al., 2012, Yuan et al., 2012).

1. Theoretical Foundations and Key Physical Principles

The CDAF solution arises in axisymmetric, radiatively inefficient accretion flows with low viscosity, where hydrodynamic or MHD turbulence fails to transport all of the liberated angular momentum outward, and the entropy gradient (generated by viscous heating exceeding cooling) produces an inward-increasing entropy profile. This triggers convective instability, as indicated by the Høiland or Solberg–Høiland criteria, resulting in radial and vertical convective eddies that redistribute thermal energy and angular momentum. In the CDAF regime, net inward angular-momentum transport by convection nearly cancels the outward viscous (Shakura–Sunyaev α prescription) transport, yielding a low net accretion rate. The mass inflow circulates in large eddies; only a small fraction “leaks” to the hole via the turbulent cascade (Inayoshi et al., 2017, Yuan et al., 2012, Ghasemnezhad, 2017).

The CDAF model contrasts with the classical advection-dominated accretion flow (ADAF), which assumes a constant mass accretion rate, and the adiabatic inflow-outflow solution (ADIOS), which attributes declining $\dot{M}(r)$ to mass loss in outflows. In CDAF, the suppression is due to gas being “locked” in convective cells, not expelled (Yuan et al., 2012, Yuan et al., 2012, Ghasemnezhad, 2017).

2. Scaling Laws and Flow Structure

The radial profiles in a CDAF are characterized by distinct self-similar scalings. The density falls as $\rho(r) \propto r^{-1/2}$, much shallower than the ADAF scaling $\rho \propto r^{-3/2}$. The mass inflow (angle-integrated, not net) decreases linearly with radius, $\dot{M}_{\mathrm{in}}(r)\propto r$, as a result of the reduced net radial velocity $v_r \sim \nu/r \propto r^{-1/2}$. The net accretion rate (mass actually reaching the innermost radii) is nearly independent of radius and is suppressed by several orders relative to the Bondi rate: $$\dot{M}/\dot{M}_\mathrm{B}\propto (\alpha/0.01)^{0.6}$$, typically yielding $\dot{M}/\dot{M}_\mathrm{B}\sim 10^{-3}$ to $10^{-2}$ (for $\alpha\sim10^{-2}$). This suppression operates in the low-luminosity, radiatively inefficient regime $\dot{M}_\mathrm{B}/\dot{M}_\mathrm{Edd}\ll10^{-3}$ (Inayoshi et al., 2017).

The steady-state radial structure in global CDAF solutions divides into three zones:

• At large radii ($r\gtrsim R_B$), a rotational equilibrium zone with $\rho(r)\sim (1+R_B/r)^{3/2}$, negligible radial motion.
• In the intermediate zone ($R_C\lesssim r\lesssim 2R_C$), a geometrically thick torus with subsonic circulations and large-scale convection, with outflows near the poles.
• Inside the centrifugal radius ($r\lesssim R_C$), CDAF scalings hold: $\rho(r)\propto r^{-1/2}$ and $\dot{M}_{\mathrm{in}}(r)\propto r$ (Inayoshi et al., 2017).

3. Angular Momentum, Energy Transport, and Dynamics

Angular momentum in CDAFs is transported both outward by viscosity and inward by convection. The classic “convective envelope” is established when the convective angular-momentum flux (parameterized via a convective α$_c$) nearly balances the outward viscous flux. In self-similar models, this can be expressed as $J_{\mathrm{con}} = r\nu_c\Sigma\,d(r^2\Omega)/dr$ (for inward transport). The balance leads to a suppression of the radial inflow and yields the characteristic flat density profile and sluggish rotation ($\Omega<\Omega_K$), with disk thickness $H/r\sim0.1-1$, increasing for stronger convection and toroidal magnetic field (Ghasemnezhad, 2017, Ghasemnezhad et al., 2017).

Thermal energy generated by viscous dissipation ($Q^+_{\mathrm{vis}}$) is largely transported outward by convection, as quantified by the convective luminosity $$L_\mathrm{conv}\simeq\eta_\mathrm{conv}\dot{M}_\mathrm{B}c_\infty^2$$, with $\eta_\mathrm{conv}\simeq0.2$ (Inayoshi et al., 2017). At small radii ($r\lesssim 10^{-3}R_B$), thermal conduction becomes dominant, flattening the temperature profile and terminating the CDAF regime (Inayoshi et al., 2017). At these radii, the flow transitions into a thin, optically thin disk structure.

4. MHD Effects, Outflows, and the CDAF–ADIOS Distinction

Analysis of global, two-dimensional hydrodynamical (HD) and magnetohydrodynamical (MHD) simulations demonstrates that the classic CDAF solution is realized only in non-magnetized, convectively unstable flows (Yuan et al., 2012, Yuan et al., 2012). In MHD flows, the Høiland criteria generally indicate convective stability throughout most of the domain, so the CDAF mechanism is not operative. Instead, systematic inflow and outflow are observed: in HD, buoyancy drives hot, entropy-rich parcels outward; in MHD, angular momentum is efficiently transported outward by tangled magnetic fields, producing centrifugal “micro–Blandford–Payne” outflows (Yuan et al., 2012).

Simulations reveal systematic thermodynamic and kinematic differences between inflow and outflow, inconsistent with purely convective (CDAF) dynamics. Temperature of HD outflows exceeds that of inflow at a given radius, while in MHD, outflows carry significantly higher specific angular momentum. This supports an ADIOS framework (real outflows removing mass) rather than a CDAF one (mass “locked” in eddies). CDAFs may coexist with ADIOS-type outflows in some regimes, but do not globally dominate the dynamics in MHD (Yuan et al., 2012, Yuan et al., 2012).

5. Self-Similar and Resistive Solutions with Magnetic Fields and Outflows

Extensions to the CDAF framework incorporate large-scale toroidal and vertical magnetic fields and explicit magnetic resistivity. In height-integrated (1.5D) models, the presence of ordered fields modifies the disk geometry, with toroidal fields thickening and vertical fields thinning the disk via pressure and tension gradients, respectively (Ghasemnezhad, 2017, Ghasemnezhad et al., 2017). The generalized mass accretion rate is modeled as $$\dot{M}(r)=\dot{M}_\mathrm{out}(r/r_\mathrm{out})^s$$, $s>0$, to capture the influence of outflows.

For strong outflows (large $s$), the disk becomes hotter, thicker, rotates more slowly, and accretes more rapidly at outer radii, but the net accretion to small radii is diminished (Ghasemnezhad et al., 2017). Increasing the magnetic field strength further thickens, heats, and accelerates the flow, while resistivity effects are subdominant, mainly modifying infall speed and surface density by modest factors.

Thermal bremsstrahlung is the principal cooling mechanism in these regimes. The broadband emergent spectrum is controlled by the surface density, temperature, and aspect ratio of the flow, with the spectra peaking at $h\nu \sim 40\,\mathrm{keV}$ for typical parameters (Ghasemnezhad et al., 2017).

6. Astrophysical Implications and Observational Diagnostics

The CDAF solution yields very low black hole feeding rates and luminosities, with dimensionless accretion rates $\dot{M}/\dot{M}_\mathrm{Edd}\sim10^{-7}$–$10^{-6}$ and Eddington ratios $L_\mathrm{bol}/L_\mathrm{Edd}\lesssim10^{-7}$, matching the quiescent emission of galactic nuclei such as Sgr A*, M31, and M87. This provides a self-consistent explanation for observed sub-Eddington states without necessitating feedback regulation (Inayoshi et al., 2017).

Key observational signatures of CDAFs (especially magnetized CDAFs) include:

• Very flat radial density profiles
• Large geometric thickness ($H/r\sim0.5$)
• Low radiative efficiency and hard X-ray spectra (high Compton $y$-parameters)
• Weaker variability and suppressed short-timescale modulations due to slow accretion
• Potential for reduced reflection features and spectral flattening due to energy transport by convection (Ghasemnezhad, 2017).

In practice, the observed properties of accreting black holes are often best matched when ADIOS-like outflows are combined with convective suppression of inflow, with pure CDAF dynamics being subdominant—especially in the presence of dynamically significant magnetic fields (Yuan et al., 2012, Yuan et al., 2012).

7. Limitations, Open Problems, and Future Directions

The viability of CDAFs as the global structure of hot accretion flows depends critically on the competition between convection, outflows, and magnetic stresses. Simulations show CDAFs may describe the innermost regions of purely hydrodynamic, low-$\alpha$, radiatively inefficient flows; however, in realistic astrophysical contexts where MHD turbulence and outflows are ubiquitous, the ADIOS scenario is typically favored.

Ongoing challenges include:

• Establishing the conditions under which CDAFs can form and persist in the presence of tangled or large-scale magnetic fields.
• Quantifying the transition between CDAF- and ADIOS-dominated regions, especially in 3D MHD simulations with realistic microphysics.
• Distinguishing CDAF and ADIOS observational diagnostics (spectral, timing, and outflow properties).
• Further modeling the impact of thermal conduction and resistivity at the smallest radii.

Future work will require extended simulations with larger dynamic range, improved turbulence models, and direct comparison with multiwavelength observational data to clarify the role of CDAFs within the broader landscape of radiatively inefficient accretion physics (Inayoshi et al., 2017, Yuan et al., 2012, Ghasemnezhad, 2017, Ghasemnezhad et al., 2017, Yuan et al., 2012).