A&A 420, 1-16 (2004)
DOI: 10.1051/0004-6361:20035741
B. Czerny1 - A. Rózanska 1 - M. Dovciak2,3 - V. Karas2,3 - A.-M. Dumont4
1 - Copernicus Astronomical Center, Bartycka 18, 00-716 Warsaw, Poland
2-
Astronomical Institute, Academy of Sciences of the Czech Republic,
Prague, Czech Republic
3 -
Faculty of Mathematics and Physics, Charles University,
Prague, Czech Republic
4 -
Observatoire de Paris-Meudon, DAEC, Meudon, France
Received 25 November 2003 / Accepted 10 February 2004
Abstract
We discuss a model of the X-ray variability of active galactic
nuclei (AGN).
We consider multiple spots that originate on the surface of an
accretion disk following
intense irradiation by coronal flares.
The spots move with the disk around the central black
hole and eventually decay while new spots continuously emerge.
We construct time sequences
of the spectra of the spotted disk and compute the corresponding
energy-dependent
fractional variability amplitude.
We explore the dependence on the disk inclination and other model parameters.
AGN seen at higher inclination with respect to the observer, such as
Seyfert 2 galaxies, are expected to have a
fractional variability amplitude of the direct emission
that is by a factor of a few higher
than objects seen face on, such as Seyfert 1s.
Key words: radiative transfer - accretion, accretion disks - atomic processes - galaxies: active - X-rays: galaxies - radiation mechanisms: thermal
Indeed, AGN are variable in the X-ray band (e.g. Lawrence et al. 1987; Taylor et al. 2003) but the origin of the variability cannot be examined directly because the relevant central regions still remain unresolved by direct imaging. However, studies of the X-ray spectral variability offer a direct insight into the structure of accretion flows onto a central black hole, which is assumed to power AGN and determine their spectra.
In the flare model, sudden dissipation occurs in very localized regions above the disk surface (coronal loops, in analogy with the solar corona). In this case the local irradiation flux can be orders of magnitude higher then the steady level of the disk flux itself, but the irradiation does not last very long and only a small fraction of the disk surface is irradiated at any moment.
Such an arrangement of X-ray emitting region was suggested for the first time by Galeev et al. (1979) and was developed in many papers (e.g. Abramowicz et al. 1991; Haardt et al. 1994; van Oss et al. 1993; Poutanen & Fabian 1999; Zycki 2002). The development of a flare leads to a burst of "primary emission'' as well as to the formation of a hot spot underlying the flare where roughly half of the X-ray flux is reprocessed by the disk. Irradiation of the disk surface in hydrostatic equilibrium leads to the formation of a strongly stratified medium - a hot fully ionized skin covering a cooler, more neutral zone - and so the X-ray spectrum resulting from reprocessing should contain signatures that are characteristic of multi-temperature gas (Nayakshin et al. 2000; Ballantyne et al. 2001; Rózanska et al. 2002). The irradiating flux is locally very large, considerably exceeding the stationary energy flux that is dissipated inside the disk (e.g. Nayakshin 2000; Ballantyne et al. 2001). The resulting ionized skin is relatively optically thick (Collin et al. 2003).
The flare/spot model is a possible (although not unique) explanation of
narrow emission features, which have been reported in 5-6 keV
X-ray spectra of several AGNs and interpreted in terms of
localized iron-line emission (Turner et al. 2002, 2004; Guainazzi 2003;
Yaqoob et al. 2003; Dovciak et al. 2004).
In the present paper we test the flare/spot model by analyzing its predictions for the fractional variability amplitude in the X-ray band. We model the local spot/flare spectrum as a sum of flare (primary) emission of a power law shape and a spot (reflected) emission. The spot emission is determined by using the coupled TITAN/NOAR codes (Dumont et al. 2000) to solve the radiative transfer and the code of Rózanska et al. (1999) to calculate the hydrostatic equilibrium. We assume a random distribution of spots and flares across the disk and we account for their motion during the time-integrated observation. We consider both non-rotating (Schwarzschild) and rapidly rotating (Kerr) black holes, and we apply general relativity corrections using the KY code of Dovciak et al. (2003).
In our approach, a sequence of solutions is specified by fixing
statistical properties that the flare/spot distribution is required
to obey. Given a particular solution, the corresponding energy-dependent
fractional variability amplitude
is calculated.
The model is described in Sect. 2. Results are
given in Sect. 3. The relevance of the
model to explain the spectrum and variability of a typical
Seyfert 1 galaxy is discussed in Sect. 4.
Reprocessing of the incident flux is computed with the code TITAN (Dumont et al. 2000). The radiative transfer is solved both for the continuum and for lines, because the escape probability approach does not give reliable results for the lines (Dumont et al. 2003). The Comptonization is taken into account by coupling with the Monte Carlo code NOAR (also Dumont et al. 2000).
The vertical structure of the underlying disk is calculated with the
code of Rózanska et al. (1999). As our principal model, we consider
an irradiated disk in hydrostatic equilibrium, which requires iterations
between the disk structure and the radiation transfer computations.
These iterations are carried out as described by Rózanska
et al. (2002). The computations of equilibrium are extremely time-consuming
because of coupling between the equations of radiation
transfer for the lines and for continuum. Therefore, in the present paper
we performed computations at a certain representative radius in the disk
(
in terms of gravitational radii,
)
with the aim of obtaining a typical reflection spectrum,
.
Local reflection spectra at other radii were obtained with the help of
scaling.
The fundamental parameters of the model are the mass of the black hole, M,
and the mean X-ray luminosity,
.
This X-ray luminosity
is assumed to originate from
flares per second,
on average.
The actual number of flares fluctuates at any given moment
around this mean value according to a Poisson's distribution.
Flares are assumed to co-rotate with the underlying disk.
A fraction of the radiation from each flare reaches the observer as the direct component. A part of flare photons is reprocessed by the disk surface, producing the reflection component.
In the model, a flare occurs at height h above the disc, it
irradiates the disk surface and creates a hot spot.
Irradiation is greatest and is strictly perpendicular just
below the flare (in the centre of the spot where it first appears).
Away from the centre, the irradiation slowly decreases over
distance up to
from the central point
(with the inclination angle of the incident X-rays increasing)
and it drops fast further out. We start by simplifying
our analysis and assuming that the irradiation does not show
any gradient across the spot radius up to
.
Then it
drops sharply to zero. Linked with the spatial dependence of
illumination across the spot surface there should be also
a corresponding change in the incident angle. However, this incident
angle is about 60
at most of the spot surface so
we assume illumination to be isotropic.
In this way, we consider the reprocessed radiation as coming from
uniformly radiating circular spots.
Given the incident radiation flux of a single flare,
,
and the spot spectrum
,
the reflection is calculated through simple scaling,
![]() |
(1) |
The nine input parameters
of the model, M,
,
,
,
,
F0,
,
,
and
fully define the flare
distribution, including the total mean luminosity of the source,
.
In practice it may be more convenient to use
as a free parameter instead of, for example,
.
Prediction of the observed spectrum involves two additional
parameters which are inherent to any model involving an accretion disk
around a rotating black hole: inclination angle of the observer, i,
and dimension-less angular momentum of the black hole, a. These
parameters span the range
(
the disk axis, while
is the disk plane), and
(a=0 refers to a non-rotating black hole, while a=1 corresponds
to a Kerr black hole with maximum rotation).
We first derive the mean luminosity,
.
Uniform distribution of
flares over the disk surface means that probability of a flare
occuring at a given radial and azimuthal position, (
),
is equal to
![]() |
(6) |
![]() |
(8) |
![]() |
(10) |
![]() |
(12) |
The mean number of flares during an observation of duration
is given by
Since some of the flares existing during the observation were born before
the beginning of the observation or they will disappear after the end of
observation, the mean observing time of any single flare is
clearly shorter than the flare duration time. We therefore introduce
the effective
life-time,
,
which is a combination of the two
timescales:
![]() |
(14) |
Therefore, each simulated observation will show clear
dispersion of the total luminosity around the assumed mean value,
.
We can calculate this dispersion analytically,
deriving the mean number of flares during observation
from Eq. (13). Luminosity of a single flare
is specified by Eq. (2) and combined with
Eqs. (4) and (11), assuming
random fluctuations around the mean value of flare number,
as well as random variations in the distribution of
flares and spots over the disk surface.
The time-integrated luminosity is given by
![]() |
(15) |
The mean (expected) value of this quantity is defined in the same way as
in Eq. (8). This expression simplifies to
![]() |
(16) |
![]() |
(17) |
![]() |
= | ![]() |
|
= | ![]() |
||
![]() |
|||
![]() |
(18) |
![]() |
(19) |
![]() |
(21) |
![]() |
(22) |
We stress that Eq. (23) for the normalized variance does not include relativistic corrections. It can be applied to the case of a non-rotating Schwarzschild black hole with acceptable accuracy, but it turns out to be quite inaccurate for a maximally rotating Kerr solution.
Naturally, spots should only be generated between the adopted
inner radius of the disk,
,
and the outer radius,
.
To model a source with a prescribed
luminosity, non-zero spot size has to be taken into account. A correction
factor is obtained numerically through iterations.
The spot radius
is constrained to be much smaller
than
,
and also it must be
smaller than
.
![]() |
(25) |
![]() ![]() |
(26) |
An example of the generated spot distribution is shown in Fig. 1.
Each spot is localized at the disk surface at the moment of its
generation. Spots occurring closer
to the black hole are brighter than those which happen to be born farther out,
because we assumed the power law scaling with
in Eq. (4).
![]() |
Figure 1:
Illustration of the spot distribution created by flares that
illuminate the underlying disk. The plane of the disk is shown, as viewed
from the top, along the symmetry axis; (x,y) coordinates are scaled with the black
hole gravitational radius,
![]() |
Open with DEXTER |
The belts are relevant for observational properties inferred from the model,
because actual observations have, indeed, some finite (and non-negligible)
duration,
.
If the integration time is long enough, some
spots may complete one whole revolution, or even more
full orbits, and so the length of the belt can be much longer
than the circumference of a circle at the corresponding radius. This is
particularly relevant for spots in the innermost part of the disk where
motion is fast. The right panel of Fig. 1 shows an example
of the resulting coverage of the disk surface by radiation belts.
Computation of the time integrated spectrum is therefore performed
as the integration of the local intrinsic emissivity of the belts.
The locally emitted flux must be renormalized to take into
account that the emission now originates from the belt instead of a
spot. Also, belt size is a function of radius
.
We find
![]() |
(28) |
The routine employs pre-stored data tables to look up and interpolate four
relevant quantities. These are (i) the energy shift of photons,
(ii) relative delay of their arrival time, (iii) magnification of
radiation (lensing effect) in the Kerr metric, and (iv) the local emission
angle (to account for directional anisotropy of the
emission). Given the
local (intrinsic) emissivity in the disk plane, the predicted (observed) spectrum can
be integrated. To achieve higher accuracy than it was possible
with previous routines, the resolution in energy and time can be controlled
and the grid covering the disk plane can be adjusted. Light rays are
integrated in Kerr ingoing coordinates, so that the frame-dragging effect
is tracked even very close to a fast rotating black hole (Misner et al. 1973). The equation of geodesic deviation is integrated
separately to ensure precise determination of the flux magnification due
to lensing near caustics.
In general, the disk emissivity is characterized by radiation energy
and the related photon flux produced
in the equatorial plane as function of emission energy
,
As mentioned above, the influence of relativistic corrections on the local spectrum
can be parameterized by the angular momentum of the black hole, a,
and the inclination angle of an observer, i. On the other hand, by scaling
lengths with
and time intervals with
,
one
conveniently ensures that graphs of predicted spectra do not explicitly show
dependence on the black hole mass.
Relativistic effects, particularly in the case approaching a maximally rotating Kerr solution, significantly affect the shape of the radiation spectrum (especially at large inclination angle). Also the magnitude of the observed luminosity as well as the observed variance turn out to be more sensitive to inclination if relativistic effects are taken into account.
![]() |
Figure 2:
Local reflection spectrum of an individual hot spot,
which arises due to irradiation of the
disk surface by a flare. Left: broad energy range is shown
(log-log scale) with the intrinsic spectral resolution
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
We computed several radiation spectra representing different realizations of the same mean spot distribution. In this way we simulated independent time-integrated observations of a given source. We generate each distribution independently which means that we assume the lapse of time among simulated observations to be larger than the life-time of flares.
Fractional variability amplitude was then calculated in the standard way
(see e.g. Edelson et al. 2002), as the energy-dependent variance of
radiation flux normalized by the mean flux at a given energy.
The variance was determined numerically, including relativistic
corrections, from a finite number of realizations, K. Then the
relative error of the normalized dispersion
was estimated as
![]() |
(32) |
Radiative transfer computations were performed for flare parameters
similar to those used by Ballantyne et al. (2001). We adopt the black
hole mass of
.
The irradiating X-ray flux,
erg cm-2 s-1, extends from 1 eV to 100 keV
in the form of a power law (photon index
). Iterations between the
vertical disk structure (in hydrostatic equilibrium) and radiative transfer
are performed untill convergence is achieved (Rózanska et al. 2002).
We choose the dimensionless accretion rate of
,
(
is scaled to the Eddington rate,
g s-1), and
radius
as typical values where local computations are
carried out.
Hence, the ratio of the X-ray incident flux to the
disk flux due to the internal dissipation for an assumed incident flux
and accretion rate is equal to 144. The disk thickness is
determined self-consistently by proceeding from the disk
surface down to the equatorial plane and adopting the diffusion approximation
deeply below the disk atmosphere (for Thomson optical depth 4). The intrinsic spectrum of a spot (i.e. the reflection component) is
shown in Fig. 2. Notice that the disk surface is
highly ionized, and so the reflection is efficient also at low energies.
Sulphur and silicate lines are the strongest. Several other lines from highly ionized species, such as oxygen and carbon, are also clearly visible.
Table 1:
Equivalent widths of the most intensive
emission lines from the local spot model (this paper),
measured (i) with respect to the reflected continuum,
,
and (ii) to the sum of the
incident plus reflected continuum,
.
Only lines
with EWs higher than 1 eV are listed.
An expanded region of the iron line (in linear scale) is shown in
the right panel of Fig. 2.
The iron line is strong and multi-component.
This emission is mostly due to helium-like and hydrogen-like ions.
The two-component structure
is clearly visible, although the 6.7 keV component slightly dominates.
Narrow-line components are accompanied by a broad shoulder due to
Comptonization in the disk surface layers. The iron edge is quite deep,
as is usual in the case of a highly ionized medium, and complex, as we
actually see two separate edges due to Fe XXV and Fe XXVI,
at 8.83 and 9.28 keV
correspondingly.
These features are also considerably
smeared by Compton scattering.
![]() |
Figure 3:
The reflection spectrum from the spotted disk.
The local intrinsic emissivity decreases (
![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
There is a clear distinction between the relativistic smearing of a monochromatic line and the shape of the relativistically smeared realistic reflection component. From Fig. 3 we see that the sharp blue edge, which is characteristic for a monochromatic line, is replaced by a much smoother feature. This is due to the presence of a double-peaked and Compton-broadened iron line in the realistic spectrum, as well as due to the contribution of the complex and Compton-smeared iron edge.
Weak emission lines of other elements are also noticeable in smeared spectra. Therefore it is not surprising that traces of these relativistically broadened emission lines are noticeable in the observational data (e.g. Rózanska et al. 2004 for Ton S180; Kaastra et al. 2002 for NGC 5548). Soft X-ray lines as strong as those claimed to be detected in the data by Branduardi-Raymond (2001) and Mason et al. (2003) are not expected within the framework of our model.
This result clearly supports the argument (Zycki et al. 1997; see also Bao et al. 1998; Young et al. 1998; Martocchia et al. 2000; Gondoin et al. 2002; Ballantyne et al. 2003) that proper interpretation of the observational data requires consideration of the relativistic broadening of the entire reflected component, i.e. the reflected continuum together with the iron line, instead of treating the two components separately.
![]() |
Figure 4: The dependence of the luminosity of the uniformly radiating disk belts on the inclination angle for the maximally rotating Kerr black hole a=0.998. Power law spectrum and no radial dependence of the intensity were assumed. |
Open with DEXTER |
The radiation originating in the innermost part of the disk shows
significantly different dependence on the inclination angle than the emission
of the outer parts.
We illustrate this by assuming, as a special example, that the disk surface
is radiating uniformly in the local frame co-rotating with the disk.
No separate spots and no radial dependence of
luminosity are present. We consider two examples of such spectra:
(A) the local emission is a pure power law with an energy index
;
(B) the local emission is a pure reflection component as described
in Sect. 3.1. In both cases, we consider the emission coming from
three detached regions: (1) an inner ring,
,
(2) an intermediate ring,
,
and (3) an outer ring,
.
We assume no limb darkening (locally isotropic emission).
In Fig. 4 we consider a strictly power-law emissivity
profile and we show the corresponding behaviour of radiation flux at
6 keV on the inclination angle. We notice that the emission from the outer
ring roughly follows a purely geometrical effect (trivial in the flat space),
i.e. the decrease proportional to ,
where i is the inclination angle
of observation. However, the emission from the innermost part is more complicated:
moving the observer from the disk axis towards its plane (i.e. i increasing),
the observed flux first brightens, peaks at some intermediate inclination, and
finally it decreases when the disk is inclined very strongly. This behaviour
is due to competition between the relativistic boosting and reduction
of the projected radiating area. The plot does not depend on the energy
at which the monochromatic luminosity is measured, because the shape of the
observed spectrum does not depend on inclination in the absence of
spectral features.
In the case of the reflection component the results are practically identical outside the region of the intense iron line. On the other hand, there is a more complex behaviour within the line range, as can be seen from Fig. 3.
The assumption of radius-independent emissivity (between
and
)
and zero emission beyond the outer edge is rather
artificial. Therefore we also show the results for the disk whose radial
emissivity normalization decays with the radius as
(
). In Fig. 5 observed luminosity
is plotted as a function of inclination. We see that the
dependence is now even stronger than for the emission from
the inner ring in Fig. 4, because with the adopted
steep emissivity law most of the radiation is generated below
.
A decline in luminosity is observed only at very large inclination, i.e.
close to the disk equatorial plane (
,
probably unrealistic
for real objects). Additionally, at such inclinations very close to 90
multiple images appear (Zakharov & Repin 2003).
![]() |
Figure 5:
The dependence of the luminosity of uniformly radiating disk
on the inclination angle with isotropic local emission, limb brightening
(
![]() ![]() ![]() |
Open with DEXTER |
Table 2:
Exemplary flare distribution properties for the black hole/disk system
seen at a moderate inclination,
.
Parameters above the double line are assumed while the parameters below the
double line are calculated as described in Sect. 2.2.
This dependency on inclination angle is only seemingly in contradiction with several other examples of such a trend published in the literature (see e.g. Gierlinski et al. 2001; Ebisawa et al. 2001; Bhattacharyya et al. 2001). There are three reasons why our results are different.
First, the emissivity law is usually assumed to be less concentrated towards the center. When the emission from the multi-black-body Keplerian disk around a low mass object is calculated (as in all papers mentioned above) the emissivity goes down to zero at the marginally stable orbit due to the assumed boundary conditions. This is in accordance with the standard thin disk model. However, dissipation in the hot corona does not need to follow this law. Also the presence of the large-scale magnetic field may eventually lead to more dissipation close to the marginally stable orbit (or even below it) than predicted in models based on zero-torque condition (see Afshordi & Paczynski 2003 for the discussion). Models of the broad iron line actually indicate a high concentration of the generated radiation flux towards the center (e.g. Wilms et al. 2001; Martocchia et al. 2002).
Second, many computations have assumed some kind of limb darkening.
For illustration we include in Fig. 5 the same
limb-darkening profile as adopted by Laor (1991),
.
We see that in this case, under the assumed radial
emissivity law, the flux still rises
with the inclination but more slowly.
The angle-dependent computations of the Compton reflection component show
this effect clearly. The increase of emissivity with inclination strongly
depends on the ionization state of the reflecting material. For the incident
angle of
of the incoming beam the ratio of local intensity at
to that at
at
6 keV is equal to 9
(for a weakly ionized
medium) and 5 (for a strongly ionized medium; cf. Zycki & Czerny 1994;
George & Fabian 1991; Martocchia et al. 2000). Therefore,
in Fig. 5 we show
an example of the flux dependence with limb brightening
roughly representing
realistic reflection by the strongly ionized medium. In this case
the flux rises almost by 50% with the change of inclination angle
from
to
.
However, in the total spectrum the
primary emission dominates over the reflected component. The angular
distribution of this emission is not known so in further considerations
we assume isotropic emission.
The dependence of the observed radiation on the inclination angle has two important consequences. First, we recognize the tendency of the spectral shape on the inclination angle, as already illustrated in Fig. 3. Second, at low inclination angles the emission from the innermost part of the disk is under-represented in comparison with the emission from the outer parts, while at higher inclinations the opposite tendency appears. Since in the model outlined in our paper most of the variability is due to rare but luminous flares from the inner region, we find strong dependence of the variability properties on the inclination angle. We will discuss this in more detail in Sect. 3.3.
The model depends on a number of parameters but most of them are well constrained by observations, and hence the actual freedom is not so large.
The mean number of flares at any moment,
,
cannot be less than
10 because observed lightcurves do not show
strong dips even on short timescales. Neither it can be very large, because
in that case the variability amplitude would be far too small.
To quantify variability, we use a standard rms as its measure.
The typical normalized rms of AGN is known to be of
order of 0.3 (see e.g. Abrassart & Czerny
2000; Uttley et al. 2002; Markowitz et al. 2003).
The parameter
is expected to be
3
on the basis of theoretical arguments as well as successful fits of some
iron
profiles in such sources, including MCG-6-15-30
(e.g. Iwasawa et al. 1996; Sulentic et al. 1998; Nandra et al. 1997;
Guainazzi et al. 1999), although more complex distributions were also advocated
(see Merloni & Fabian 2003 and references therein).
In Table 2 we give several combinations of parameters
that we considered together with the properties of the chosen flare distribution.
In most cases, we fix the mass of a black hole at a representative value
,
for which the local computations of the disc structure and radiation
reprocessing were performed. The inclination angle of
was
taken as a representative value (the view of the nucleus is frequently obscured
at high inclinations by the presence of a molecular or dusty torus, so
the average inclination in type 1 objects must be significantly lower than
). The opening angle of the dusty torus is statistically constrained
by the ratio of type 2 to type 1 objects. It was estimated to
be most probably between 2 and 4 (see Krolik 1999; Veron-Cetty & Veron
2000). Therefore, the torus opening angle is of the order of
,
a typical type 1 object is seen at
,
and a typical type 2 object at
.
![]() |
Figure 6: Two random realizations of the same spot distribution (thin lines) and the original local emission (thick line histogram) for model A; the ratio of the two random realizations are shown below. The local input spectrum is the sum of the primary (flare) emission and the reflected (spot) component, without any limb-darkening. |
Open with DEXTER |
An example of the X-ray spectra from two random realizations of
model A is shown in Fig. 6. The intrinsic
spectrum of spots/flares was assumed to be a combination of the
intrinsic flare emission and reflected/reprocessed spot emission.
The observed spectral features are now weaker than the example in
Fig. 2
where only the spot emission was taken into account. The two
spectra of statistically identical distributions of spots differ in
normalization, but otherwise their form is very similar. Variations
are slightly enhanced at the energy where the strong iron line occurs.
The tails of the line vary more significantly because this emission comes
from the several most luminous spots at the innermost part of the disk,
and so random variations in the position of those few spots
are relatively frequent and important.
![]() |
Figure 7:
Fractional variability amplitude calculated as a function of energy.
The model A (
![]() ![]() ![]() |
Open with DEXTER |
This effect is best seen in plots where the energy-dependent fractional variability amplitude is calculated from many random realizations. Such a plot (for model A) is shown in Fig. 7.
In the case of a Schwarzschild black hole,
variations are almost independent of energy. There is only
marginal enhancement of variability in the red wing of the iron
line, around 5 keV. There is also a small dip in the plot
at 7 keV. However, the overall energy dependence
is not strong, and so the value of the dispersion at 1 keV
is representative for the entire spectrum (Table 2).
The energy dependence of
is stronger in the case of
emission from the disk surrounding a fast-rotating Kerr black hole.
We show the corresponding plot for model D in the lower panel of
Fig. 7. Here we adopted a higher value of the mean number
of clouds, because otherwise the obtained variance would be
unrealistically high in the context of AGN. We see in this case strong
variations in the region of the red wing of the iron line, around 5 keV,
and again a small dip in the plot at
7 keV. Still, the
overall variability is well characterized by the dispersion at 1 keV, as
given in Table 2.
The fractional variability amplitude in this case must be calculated with
the relativistic corrections.
The effect depends also on the inclination angle. An example for the
model D but computed
this time at a much higher inclination angle of
is also shown
in Fig. 7. At such a high inclination, line wings are still
broader than at
,
influencing the entire energy range
between 4 and 8 keV.
There is a noticeable difference between the variability level
of the source observed at
versus the case of
(the value typically expected for broad-line galaxies on
the basis of an AGN unification scheme). The fractional variability amplitude
is enhanced by 47%. Further increase, by a factor of 3.1-4.6 in the
variability level, is expected for objects observed at
-
.
Such inclinations are typical of Seyfert 2 galaxies.
In these objects, in the soft X-ray band, we do not see the direct emission from
the nucleus. Instead, in some of them the direct emission is visible in hard
X-rays (Done et al. 2003), which allows us to test our prediction
regarding the expected level of variability.
![]() |
Figure 8: The dependence of the normalized dispersion measured at 1 keV on the inclination angle of observation for model A (upper panel, continuous line), model J (upper panel, dashed line) and model L (lower panel, middle curve). Two other curves in the lower panel show the result of increasing the mean number of flares to 3000 and decreasing the mean number of flares to 300. Number of exposures was 100. |
Open with DEXTER |
Several authors have recognized the importance of general relativity effects for the observed dependence of the variability properties on the inclination. Abramowicz & Bao (1994) noticed the flattening of the power spectrum slope with an increase of the inclination angle. Their approximate results were in general confirmed by Xiong et al. (2000) for a Schwarzschild geometry. Karas (1997) discussed Fourier-phase analysis of the spotted disk model.
The quantitative effect of the inclination depends on details of the
assumed variability model. We show this dependence for model J, which is
characterized by steeper emissivity (
)
than
adopted in other cases. We see that the rise of the fractional variability
amplitude with inclination is faster than for
because more radiation comes from the innermost region where the relativistic
effects are important.
The value of the fractional variability amplitude given by
Eq. (23) is
independent of both the energy and inclination angle. In
Table 2 we give these analytical values as
.
We also calculate semi-analytical values
taking into account the finite size of the spot but
neglecting relativistic corrections;
a new radius is
calculated iteratively and inserted into Eq. (23) with
.
We see that the analytical formula for
gives a rather good approximation
of the exact value in the model K, which involves the Schwarzschild solution
and relatively small size of spots (less than
).
For other Schwarzschild models
with larger sizes of the spots, the simplest approximation is no longer
satisfactory, but
still provides a reasonable
approximation. Relativistic effects are
important but the formulae neglecting them can be used for
qualitative estimates.
For Kerr models, the analytical value given by Eq. (23) over-predicts the variability level at low inclination angle by a factor up to 3, but it represents roughly the correct values at large inclination angles. In this case, analytical estimates give large errors, and so accurate numerical energy-dependent and inclination-dependent results with general relativity effects are needed.
Increasing the assumed radius of the disk also increases the variance
because, in our model, larger
translates
into a smaller number of flares generated in the innermost part of the disk
where flares are the most energetic. This is clearly related to
our assumption of the uniform covering of the disk and the power-law radial
decay of flare luminosity, as given by
.
Equivalently, we can modify the
dispersion by changing the slope
:
a lower value (model F)
reduces variations, while a higher value increases it (model J).
The energy-dependent fractional variability amplitude was studied for this source
by Matsumoto et al. (2003) on the basis of ASCA data. They showed that
the overall 2-10 keV amplitude of 18-20% is enhanced in the
region of the red wing of the iron line (modeled as a broad feature around 5.4 keV)
up to 40% (but with a large error), and somewhat less at its blue wing
(around 6.4 keV). A detailed analysis was performed by
Markowitz et al. (2003) on the basis of several RXTE observations.
They found that the overall variability decreases with
energy, going from
20% at 4 keV down to
13-14% above
10 keV. There is also a slight relative enhancement at
5 keV and
at
8 keV, and a dip at
6.5 keV. A similar trend was noticed by
Fabian et al. (2002) in their long XMM observations,
but in this case the variability level was found to be slightly higher.
We model the observed fractional variability assuming parameters appropriate
for this source.
The mass of the black hole is rather small on the basis of X-ray power
spectrum analysis:
follows from the normalization method
(Czerny et al. 2001), while
has been derived from the
high-frequency break model (Uttley et al. 2002). However, in view
of the uncertainties, a mass as large as
is still compatible.
No results have been obtained from reverberation studies for this source.
The mass was estimated to be
(Czerny et al. 2001)
using the disk method. Reynolds (2000) also argues for this higher
value on the basis of the bulge-black hole mass relation (Magorrian et al. 1998).
We thus adopt
in our computations.
X-ray luminosity varies significantly, but taking the typical
2-10 keV flux to be
erg s-1 cm-2
(Weaver et al. 2001) we can estimate the total X-ray flux (1 eV-100 keV)
to be of the order of
erg s-1(assuming the Hubble constant H0=70 km s-1 Mpc-1).
The bolometric luminosity of this source can be estimated either from
the 2-10 keV luminosity, taking into account the bolometric correction of
27 (Padovani & Rafanelli 1988), or directly from the broad-band
data corrected for extinction. The first method gives
erg s-1 while the second one gives
erg s-1 (Niko
ajuk,
in preparation). This means that for the adopted mass of the black hole
the source radiates at
10% of the Eddington luminosity.
We considered three values of the mean number of flares: 300, 1000, and 3000. The corresponding models are denoted L, P and Q in Table 2. We choose the duration of the observations to be 103 s, the typical duration of continuous observation in the data, but the result does not depend strongly on this value. The adopted life-time of flares, 105 s, is chosen to reproduce the knee in the power spectrum of this source, as determined by Uttley et al. (2002). Therefore, we model the day-to-day variability (and the result does not apply to longer timescales).
Comparison of the observed
in short (day-to-day)
timescales (see Fig. 5 of Markowitz et al. 2003) with theoretical
level of variability indicates that the most appropriate mean number
of flares is between 300 and 1000.
There is evident similarity in the energy dependence of
between observed data and the theoretical diagram:
variability in the region of the iron line red wing
(
5 keV) is enhanced, and a dip occurs at about 7 keV.
However, the dependence on energy is generally much weaker
in the model than in the data. Also, the model does not show the
observed strong rise of the
towards low energies.
This presence of strong variation at 3-4 keV in the data and
the lack of such variations in the model may indicate that
this energy
band is still affected by varying conditions in the warm absorber
(Inoue & Matsumoto 2003).
The flare/spot model is an attractive explanation of the X-ray spectra and
variability of AGN. In this scenario X-ray emission is generated both in hot
magnetic loops above an accretion disks and in the bright spots created under
the loops by strong irradiation. In the present paper we tested this model by
analyzing the predicted fractional variability amplitude,
.
We derived simple analytical formulae which allow us to estimate the level of variability from the assumed mean number of flares, flare duration, integration time of a single observation, the ratio of inner to outer disk radius and the radial dependence of the flare luminosity. This energy-independent and inclination-independent formula roughly applies to the case of a non-rotating black hole and small size of the spots. The formula is based on an assumption of a uniformly covered disk surface but it can be readily generalized to another kind of distribution.
Larger spots and/or a fast rotating black hole, with an inner disk radius close to
the marginally stable orbit require a numerical approach, and we show the
dependence of the resulting
both
on energy and inclination angle.
This leads to a firm prediction of our model
which can be used to test the basic scenario. The model shows that
if the disk extends
close enough to the black hole, the general relativity effects lead to
significant
enhancement of the variability at large inclination angles of observation.
Therefore, Seyfert 2 galaxies, if intrinsically identical to Seyfert 1
objects
but viewed at a larger angle, should exhibit statistically higher
when measured at the same energy. A factor of 3 difference
between Seyfert 2 and Seyfert 1 galaxies is expected
for a Kerr black hole and a factor of 1.4 for a Schwarzschild black hole.
The observational evidence for a trend in this direction is scarce
at present, but it
may be used to support or reject our view more firmly in the future.
However, some values can be given already.
The mean variability level of Seyfert 1 objects is 18%
on the short timescales of 1 day in the sample of Markowitz et al. (2003).
The Seyfert 2 galaxy NGC 4945 is seen through the torus and it
displays variability at the level of
40%
in the hard X-ray band. Another example,
NGC 7582 (Mihara et al. 2000), has
revealed the normalized variability amplitude in the hard X-ray band
to be about
30%.
Studies of more Seyfert 2 objects are clearly needed.
The overall intra-day variability level of Seyfert galaxies
is well explained by the flare model if the mean number
of flares is of the order of 30-100 for an assumed non-rotating black hole,
or 300-1000 for a fast rotating black hole. The number of flares
required is much larger than the usual expectation of 10 flares.
Such a small number of flares is predicted when an assumption is made that
all flares have the same luminosity. In our more realistic model we allow
for the flare luminosity to depend on the occurrence radius. When we assume that
the flare luminosity scales with radius as
r-3, most of the
source X-ray luminosity comes from a few flares generated in the innermost
part of the disk which enhances the variability. This trend is clearly seen
from our analytical expression (23) for the normalized variance,
which, for
,
and
,
reduces to
Eq. (24). Therefore, the total number of flares,
can still be quite high for a source with moderate variance.
The energy dependence of
is generally weaker in the model than in the data. Observed
variations show trends with energy in 1-10 keV band. For example, in
Markowitz et al. (2003) Akn 564 varies at the level of 18-22%, IC 4329A
at the level of 11-14% and MCG-6-30-15 of 14-21%, depending on the
considered energy. In our models the trends with energy in
never exceed 2%.
This discrepancy possibly can be solved by relaxing several simplifications:
We have considered a few special cases.
A change of the life-time
of a flare from radius-independent to scaled with a Keplerian timescale,
without a change of other parameters, resulted in a fractional variability
amplitude even less dependent on the energy than previously. This is
caused by the fact that in the case of such a scaling relatively more energy
is dissipated in the outer region, so the most relativistically broadened and
variable inner region contributes less to the total lightcurve. The overall
normalization of the
depends on the proportionality constant
between the life times and the Keplerian timescale (with other parameters
fixed).
If we assume that the overall radial dependence of the dissipation should not be modified, we can consider two representative examples.
In the first case we assume that the life-time of a flare scales with the Keplerian timescale but the probability of a flare appearing at a given radius scales inversely with the local Keplerian timescale. We have therefore fewer long-living flares in the outer region and more short-living flares in the inner region, with the scaling of a single flare luminosity with radius unchanged. In this case again the fractional variability amplitude is less dependent on the energy than in our basic model. Having more flares localized in the inner region leads to reduction in the fluctuations in the relativistically smeared red wing.
In the second case we again assume that the life-time of a flare scales with
the Keplerian timescale, we still adopt the uniform distribution of the flares
across the disk surface but this time we assume that flare luminosity
decreases with radius even more strongly (
instead of
the usually adopted
)
to compensate for
an increase of the
flare life-time. Such a solution leads to slightly enhanced dependence of
the fractional variability amplitude on the energy but the effect is not
strong.
Complex dependence of the spectral shape of the reflected component on the
disk radius may introduce significant modification to the predicted energy
dependence. A clear suggestion of what is needed can be obtained
from the recent
analysis of MCG-6-30-15 by Vaughan & Fabian (2003). A constant
component with complex energy dependence is apparently needed to
formally model the fractional variability amplitude in this source (see their
Figs. 11 and 16, top panel). This component may perhaps be understood as a
reflection component (above 1 keV) and a contribution from
emission/scattering by some extended medium (below
1 keV). However,
it is not clear whether there is any possibility finding a parameter range
that would satisfy the two "opposing'' trends seen in these data: we need strong
reflection from the distant disk region to reproduce the constant
component but we need strong reflection from the innermost region to
explain the strongly relativistically broadened iron line profile also seen
in these data. Attempts by Zycki & Rózanska (2001) and Ballantyne et al.
(2003) were not successful.
Flares are not the only possibility to explain the X-ray emission of AGN. Other scenarios include the lamp-post (standing shock) model (e.g. Henri & Pelletier 1991; Malzac et al. 1998), model of gradual or rapid disk evaporation and its replacement by the hot flow (e.g. Narayan & Yi 1994; Liu et al. 2002), possibly with an outflow (e.g. Blandford & Begelman 1999), and the cloud model (e.g. Collin et al. 1996; Karas et al. 2000). Further work, taking into account variability issues, is needed to determine whether the flare/spot model is the most satisfactory.
Acknowledgements
We thank Piotr Zycki for very helpful discussions and we are grateful to Sergei Nayakshin, the referee, for comments which helped to improve the manuscript. Part of this work was supported by grant 2P03D00322 of the Polish State Committee for Scientific Research and by Jumelage/CNRS No. 16 "Astronomie France/Pologne''. V.K. and M.D. acknowledge support from grants GAUK 188/2001 and GACR 205/03/0902 and M.D. acknowledges support from the grant GACR 202/02/0735 in the Czech Republic.