Propagating fluctuations in the mass accretion rate of a precessing flow as a power spectral model for black hole binaries
This model and description are due to Adam Ingram.
The power spectrum of black hole binaries (BHBs) in the rise to outburst generally consists of a quasi-periodic oscillation (QPO) and additional band limited noise. This can be phenomenologically modelled using a number of broad Lorentzians for the band limited noise and a narrow Lorentzian for each harmonic of the QPO. Here, we assume a physical origin for this variability in a truncated disc / hot inner flow geometry. Mass accretion rate fluctuations are generated everywhere in the flow (primarily) at the local viscous frequency and propagate inwards towards the black hole. We also assume that the entire flow is precessing due to frame-dragging (Lense-Thirring precession), which gives rise to a QPO. The key assumption is the surface density profile as this sets the precession (and therefore QPO) frequency but also, by mass conservation, sets the viscous frequency as a function of radius. Further details of the model, propfluc, are included in Ingram & Done (2011), particularly in the appendix.
We assume that the surface density is given by a smoothly broken power law, consistent with the results of general relativistic magneto hydrodynamic (GRMHD) simulations. The break occurs at the bending wave radius, rbw, with the power law dependence on radius parametrised by zeta for r>>rbw and by lambda for r<<rbw. The sharpness of the break is given by another parameter, kappa. The QPO frequency is calculated self-consistently and the power spectrum of the QPO is represented by Lorentzians centred at f_QPO, 2f_QPO, 3f_QPO and 1/2f_QPO to represent the 1st, 2nd, 3rd and sub harmonics respectively.
As this is a model for the power spectrum rather than the spectral energy
distribution, the process for loading the data into XSPEC is slightly
different. First of all, a power spectrum can easily be created from a light
curve using powspec
from the XRONOS package (for example). The power spectrum
will then typically be written in the form
Parameters in propfluc:
Download the tar file propfluc.tar. If you are using this without any other models in xspec12 then untar this in a clean directory, fire up xspec12 and type @load. The model should then be set up. Otherwise, build it as normal with other local models.
Keith Arnaud, Lab. for High Energy Astrophysics, NASA/Goddard Space Flight Center
HEASARC Home | Observatories | Archive | Calibration | Software | Tools | Students/Teachers/Public
Last modified: Thursday, 04-Aug-2011 09:24:55 EDT