TCAF: Two-Component Accretion Model

• TCAF is a framework that models accretion disks with two interpenetrating flows: a cool, viscous Keplerian disk and a hot, fast-moving sub-Keplerian halo.
• It explains phenomena like shock formation, spectral state transitions, and QPOs by linking accretion rates with Comptonization in the post-shock region (CENBOL).
• Numerical simulations and XSPEC fittings using TCAF yield key physical parameters, offering insights into accretion geometry, viscous timescales, and disk–jet interactions.

The two component advective flow (TCAF) model provides a first-principles, hydrodynamic description of black hole accretion disks, rooted in the existence of two physically and dynamically distinct accretion components—a cool, optically thick, viscous Keplerian disk and a hot, sub-Keplerian, low-viscosity halo. This paradigm successfully unifies spectral and timing phenomena in both stellar-mass black holes and AGN through explicit modeling of accretion rates, shock formation, and Compton cloud properties. Numerical simulations and exhaustive spectral-timing fits have established TCAF as a comprehensive and predictive framework for interpreting state transitions, quasi-periodic oscillations (QPOs), disk–jet coupling, and spectral energy distributions across a diversity of systems (Chatterjee et al., 2016, Giri et al., 2014, Jana et al., 2017, Nandi et al., 2019).

1. Dynamical Foundations and Core Equations

TCAF is based on the decomposition of the inflowing matter into two co-existing, inter-penetrating components:

• Keplerian disk: A high-viscosity, optically thick flow, primarily in the equatorial plane, rotating nearly at the local Keplerian velocity, with small radial velocity. Mass accretion rate $\dot{M}_d$ is regulated by viscous angular momentum transport and radiative cooling—leading to a standard multi-color blackbody emission.
• Sub-Keplerian halo: A low-viscosity, optically thin, hot flow with high radial infall velocity and sub-Keplerian specific angular momentum, denoted $\dot{M}_h$. This component is advective and radiatively inefficient except near the innermost regions.

The set of height-integrated, axisymmetric hydrodynamic equations includes mass continuity, radial momentum (Euler) equation, angular momentum transport (in the disk via Shakura–Sunyaev $\alpha$-viscosity), and energy equations incorporating viscous heating and cooling (Giri et al., 2014, Chatterjee et al., 2016). The equations are generally formulated in the Paczyński–Wiita pseudo-Newtonian potential to capture Schwarzschild geometry effects.

The transition from supersonic to subsonic flows in the sub-Keplerian halo produces a centrifugal barrier, resulting in the formation of a standing or oscillating shock. Rankine–Hugoniot conditions determine (i) the shock location $x_s$ (in units of $r_g=2GM/c^2$), and (ii) the density compression ratio $R=\rho_+/\rho_-$, where "+" and "−" denote post- and pre-shock quantities (Chatterjee et al., 2016).

2. Shock Formation, CENBOL, and Comptonization

When low-angular-momentum matter in the halo approaches the black hole, it encounters a centrifugal barrier and may undergo a shock if the required hydrodynamic conditions are satisfied. The post-shock region, termed the centrifugal pressure-supported boundary layer (CENBOL), is hot (electron temperatures $\sim10^9$ K) and moderately optically thick.

The CENBOL acts as the Compton cloud: soft photons originating from the Keplerian disk are inverse-Comptonized by the hot electrons in CENBOL, producing the power-law hard X-ray tail observed in accreting black holes (Chatterjee et al., 2016, Jana et al., 2017). The analytical expression for QPO frequency in the shock-oscillation model relates directly to the shock’s dynamical properties: $$\nu_{\mathrm{QPO}} = \beta\,\left[x_s(x_s-1)^{1/2}\right]^{-1}$$ where $\beta$ is a function of the compression ratio evolution (Chatterjee et al., 2016).

The optical depth and electron temperature of the CENBOL, as well as its geometry (set by $x_s$ and $R$), govern the magnitude of Comptonization and thus control the spectral slope and hardness.

3. Numerical Simulations and Physical Realization

Viscosity and cooling naturally segregate the accretion flow into two components. Hydrodynamic simulations with spatially dependent $\alpha$-viscosity (maximal on the equatorial plane) and power-law cooling produce a steady configuration in which a cold, dense, nearly Keplerian disk forms in the midplane, surrounded by a hot, rarefied, sub-Keplerian halo (Giri et al., 2014, Giri et al., 2012). Standing shocks arise self-consistently in the inner disk; the shock location and strength depend sensitively on the specifics of the inflow angular momentum and viscosity.

Simulations further show cyclic transitions: strong cooling collapses the CENBOL, reproducing the spectral softening and transitions observed in black hole outbursts. On removal of enhanced viscosity/cooling, the Keplerian disk evaporates, demonstrating accretion hysteresis similar to that seen in X-ray binary outbursts (Giri et al., 2012).

4. Spectral States, State Transitions, and Timing Features

The ratio of the halo to disk accretion rates, termed the accretion rate ratio $\mathrm{ARR} = \dot{M}_h/\dot{M}_d$, is a primary diagnostic of spectral state in the TCAF paradigm (Chatterjee et al., 2016, Mondal et al., 2014). The canonical spectral states are mapped as follows:

• Hard state: $\mathrm{ARR} \gg 1$, strong shock at large $x_s$, extended, hot CENBOL, spectrum dominated by power-law, low/absent QPOs.
• Hard-intermediate state (HIMS): Both accretion rates significant, shock moves inward, CENBOL shrinks, type-C QPOs with rising frequency.
• Soft-intermediate state (SIMS): $\dot{M}_d \gtrsim \dot{M}_h$, weak/close-in shock, thermal disk emission dominant, sporadic type-B/A QPOs.
• Soft state: $\dot{M}_d \gg \dot{M}_h$, shock disappears, only disk present, blackbody-dominated spectrum, QPOs absent.

Transitions occur as the rates evolve, and QPO frequencies trace $x_s$ and $R$ evolution in a quantitative manner (Chatterjee et al., 2016, Mondal et al., 2014, Debnath et al., 2013). The shock oscillation model quantitatively predicts QPO frequencies as the inverse of the post-shock region’s infall time, providing a direct timing–spectral linkage.

5. Observational Implementation: XSPEC Table Model and Parameter Estimation

The TCAF solution, as implemented in a large FITS table model for XSPEC, allows direct extraction of underlying flow parameters from X-ray spectra. The five principal free parameters are black hole mass $M_{BH}$, $\dot{M}_d$, $\dot{M}_h$, shock location $x_s$, and compression ratio $R$ (Chatterjee et al., 2016, Debnath et al., 2014, Nandi et al., 2019). An additional normalization encodes geometric factors (distance, inclination, mass).

The approach is broadly applicable: TCAF-based XSPEC fitting accurately recovers the daily evolution of disk and halo rates, shock properties, and even predicts the timing features (QPO frequencies) directly from spectral fits (Debnath et al., 2013, Mondal et al., 2014, Nandi et al., 2019).

6. Empirical Validation, State Lags, and Broader Implications

Analysis of X-ray variability and timing confirms the physical separation into two viscous components: RXTE/ASM arrival time lag studies reveal delayed soft flux responses in systems with large Keplerian disks, in line with predictions of different viscous timescales for the disk and halo (Ghosh et al., 2018). Typical viscous delays of up to $20$ days are observed in low-mass X-ray binaries (LMXBs), and are essentially zero in wind-fed HMXBs, as anticipated (Ghosh et al., 2018, Jana et al., 2016). Time-resolved studies of accretion rates further yield viscous timescales (e.g., $t_{\mathrm{visc}} \sim 10$ days in MAXI J1836-194) directly from the lag between halo and disk peak rates (Jana et al., 2016).

The TCAF scheme has been applied to a broad range of sources, including AGN such as NGC 4151, where broad-band X-ray fits using TCAF yield both reliable black hole mass estimates and physical characterization of the accretion geometry (Nandi et al., 2019). The normalization parameter in TCAF fits remains stable across spectral states, providing a built-in gauge to detect additional X-ray components (e.g., jet base contributions) (Jana et al., 2017).

7. Limitations, Extensions, and Outlook

Current implementations of TCAF neglect magnetic fields and apply only to Schwarzschild geometry, with outflows and jet contributions either absent or handled phenomenologically (Giri et al., 2014, Nandi et al., 2019). Extensions to include self-consistent mass ejection (JeTCAF), Comptonization in jet bases, black hole spin (Kerr metric), and full radiative transfer are under development (Mondal et al., 2021). Numerical resolution and spectral grid density limit parameter extraction precision, and a full treatment of disc reflection and relativistic lines is necessary for high-fidelity AGN and high-spin X-ray binary fits.

Despite these, the TCAF framework robustly unifies dynamical, spectral, and timing properties in accretion systems, and remains a predictive paradigm for interpreting multi-wavelength observations and probing the physical regimes of relativistic accretion (Chatterjee et al., 2016, Giri et al., 2014, Nandi et al., 2019).