Implementation of Two Component Advective Flow Solution in XSPEC

Abstract: Spectral and Temporal properties of black hole candidates can be explained reasonably well using Chakrabarti-Titarchuk solution of two component advective flow (TCAF). This model requires two accretion rates, namely, the Keplerian disk accretion rate and the halo accretion rate, the latter being composed of a sub-Keplerian, low angular momentum flow which may or may not develop a shock. In this solution, the relevant parameter is the relative importance of the halo (which creates the Compton cloud region) rate with respect to the Keplerian disk rate (soft photon source). Though this model has been used earlier to manually fit data of several black hole candidates quite satisfactorily, for the first time, we made it user friendly by implementing it into XSPEC software of GSFC/NASA. This enables any user to extract physical parameters of the accretion flows, such as two accretion rates, the shock location, the shock strength etc. for any black hole candidate. We provide some examples of fitting a few cases using this model. Most importantly, unlike any other model, we show that TCAF is capable of predicting timing properties from the spectral fits, since in TCAF, a shock is responsible for deciding spectral slopes as well as QPO frequencies.

• The paper presents the implementation of the TCAF model in XSPEC to automatically extract physical accretion parameters from X-ray spectra of black hole binaries.
• It constructs a library of roughly 400,000 model spectra by varying key parameters, enabling precise fitting of both thermal and non-thermal emission features.
• The integration bridges spectral and timing analysis, allowing prediction of QPO frequencies and setting the stage for future enhancements including spin and jet effects.

Implementation of the Two-Component Advective Flow (TCAF) Solution in XSPEC

Introduction and Motivation

Understanding the emission mechanisms and accretion dynamics of black hole X-ray binaries (BHXBs) is indispensable for constraining compact object astrophysics. The radiative output from matter infalling onto black holes displays both thermal and non-thermal components, requiring models that can account for soft blackbody emission and the hard, power-law tail often attributed to Comptonization processes. Conventional Shakura-Sunyaev (SS73) or Novikov-Thorne (NT73) models, while effective for thermal disk emission, are inadequate for capturing the full phenomenology of BHXB spectra, especially features dominated by inverse Compton scattering in a hot corona or 'Compton cloud'.

The Chakrabarti-Titarchuk Two-Component Advective Flow (TCAF) model provides a semi-analytic framework that incorporates a Keplerian disk plus a sub-Keplerian, low-angular-momentum halo enveloping the disk. The TCAF scenario naturally produces a Comptonizing shock region (the CENBOL), whose properties set both the geometry and thermodynamics of the hard X-ray source. Previous applications of TCAF required manual fitting and custom code. The key technical development described in this work is the implementation of the TCAF model as an additive local table model in NASA/GSFC's XSPEC spectral analysis package, enabling automated extraction of physical accretion parameters directly from X-ray spectra for any black hole X-ray binary (1402.0989).

Theoretical Framework and Key Model Components

The TCAF model specifies the accretion geometry using five independent parameters: the Keplerian disk accretion rate ($\dot{m}_d$), the sub-Keplerian (halo) accretion rate ($\dot{m}_h$), the central black hole mass ($M_{BH}$), the shock location ($X_s$), and the compression ratio at the shock ($R$). The shock forms in the low-angular-momentum inflow region due to the centrifugal barrier, leading to the formation of the CENBOL. This post-shock, high-temperature region acts as the physical Compton cloud responsible for inverse-Compton upscattering of the soft disk photons.

Density and temperature profiles within the flow are determined self-consistently by solving the two-temperature hydrodynamic equations, incorporating radiative cooling (bremsstrahlung, Comptonization) and heating processes. The model computes the spectral index, electron temperature, and optical depth within the CENBOL, ensuring that the emergent spectrum is linked directly to the physical accretion mechanics.

Adjustments to the original CT95 code are made to accommodate both strong and weak shock cases, with new parameterizations for shock height and temperature. Effects such as spectral hardening and the pseudo-Newtonian potential for the space-time geometry are included, though synchrotron cooling and black hole spin are deferred for future implementation.

XSPEC Integration and Technical Implementation

To embed TCAF in the XSPEC environment, the authors generated a library of model spectra ($\sim$4\times10^5$$entries) by systematically varying the five physical parameters over astrophysically relevant domains:</p> <ul> <li>Keplerian disk rate$$\dot{m}_d$:$0.1 - 12.1$$\dot{M}_{Edd}$</li> <li>Sub-Keplerian halo rate$\dot{m}_h$:$0.01 - 12.01$$\dot{M}_{Edd}$</li> <li>Black hole mass:$5 - 15~M_\odot$</li> <li>Shock location$X_s$:$6 - 456~r_g$</li> <li>Compression ratio$R$:$1 - 4$$</li> </ul> <p>Model spectra for each parameter vector are synthesized via a FORTRAN implementation of the TCAF hydrodynamical + radiative transfer code, producing the additive table model TCAF.fits. XSPEC users can now perform parameter estimation and spectral fitting within a standard analysis workflow, using physically meaningful quantities rather than phenomenological power laws or disk models. The model normalization encodes source distance and disk inclination.</p> <h2 class='paper-heading' id='application-to-observational-data'>Application to Observational Data</h2> <p>Three RXTE/PCA spectra from the rising outburst phases of H~1743-322, GX~339-4, and GRO~J1655-40 are fitted using TCAF within XSPEC. Systematic errors and interstellar absorption are treated consistently, and Fe K$$\alpha$$lines are modeled with Gaussians as needed.</p> <p>The fitted models yield reduced$$\chi^2$values typically$\leq 2$$in hard and intermediate spectral states, indicating statistically acceptable fits with direct mapping to accretion flow parameters. In the soft state, omission of the bulk motion dominated advective flow (BDAF) and pre-jet limits fit quality, but this is a stated limitation to be addressed in future model releases.</p> <p>The TCAF fits not only reproduce the observed spectral shape but also yield unique estimations for accretion rates, shock location, and compression, parameters that are otherwise only indirectly constrained in standard disk or Comptonization models.</p> <h2 class='paper-heading' id='predictive-timing-diagnostics-qpo-frequency-estimation'>Predictive Timing Diagnostics: QPO Frequency Estimation</h2> <p>A critical innovation of TCAF is its <em>predictive</em> link between spectral fit parameters and timing properties, specifically low-frequency <a href="https://www.emergentmind.com/topics/quasi-periodic-oscillations-qpos" title="" rel="nofollow" data-turbo="false" class="assistant-link" x-data x-tooltip.raw="">quasi-periodic oscillations</a> (LFQPOs). In the shock oscillation model (SOM), the LFQPO frequency is inversely proportional to the post-shock infall time:</p> <p>$$\nu_{QPO} \sim \frac{C}{R~X_s(X_s-1)^{1/2}}$</p> <p>where$C$depends only on$M_{BH}$. Thus, knowledge of$R$and$X_s$ from spectral fits suffices to predict the dominant QPO frequency. For the selected spectra, estimated QPO frequencies from TCAF parameters are within the observational uncertainties of the measured values. This cross-modal predictive capability is not present in conventional multi-component phenomenological XSPEC models.

Implications and Future Directions

By introducing TCAF into XSPEC, the authors effectively bridge the gap between semi-analytic accretion theory and observational data analysis for black hole binaries. The most significant practical implication is the capability to extract fundamental physical accretion parameters and QPO diagnostics directly from routinely collected X-ray spectra.

On the theoretical side, TCAF's success in mapping spectral and timing features to properties of advective shocks and disk-halo configurations lends support to models where centrifugal barrier-induced structures (CENBOL, pre-Jet) play a central role in both radiative and timing variability. Future model expansions are anticipated to incorporate the effects of bulk motion, jet emission, and black hole spin, increasing the applicability to soft states and more complex spectral-timing morphologies.

Conclusion

The integration of the Two-Component Advective Flow solution into XSPEC marks a technically robust advance, enabling direct and physically motivated spectral-timing analysis for black hole candidates (1402.0989). The model's capacity to simultaneously fit spectral data and predict QPO frequencies from first principles underscores the physical realism embedded in the two-component advective paradigm. Forthcoming enhancements will further broaden its relevance for advanced black hole astrophysics.