A&A 397, 711-722 (2003)
DOI: 10.1051/0004-6361:20021079
A. Rau1,2 - J. Greiner1,2
1 - Astrophysical Institute
Potsdam, An der Sternwarte 16, 14482 Potsdam, Germany
2 -
Max-Planck-Institute for extraterrestrial Physics, Giessenbachstrasse, 85748 Garching, Germany
Received 8 March 2002 / Accepted 31 May 2002
Abstract
We present a comprehensive X-ray study of four years of
pointed RXTE observations of GRS 1915+105 in the
-state. We
interpret the behavior of the hard power law tail spectrum as
coming from inverse Compton scattering of soft disk photons on a
thermally dominated hybrid corona above the accretion disk. GRS 1915+105 shows a strong, variable reflection amplitude. As in other BHC and
in Seyfert galaxies, a correlation between the power law slope and
the reflection was found. Also, the radio fluxes at 2.25 GHz and
15 GHz correlate with the power law slope, thus revealing a
connection between the outflowing matter and the comptonizing
region in the
-states.
Key words: X-rays: stars - stars: binaries: close - X-rays: individual: GRS 1915+105
With the discovery of superluminal ejections (Mirabel & Rodriguez
1994) from the galactic transient source GRS 1915+105 and GRO J1655-40, much
observational and theoretical attention has been directed to the field
of microquasars. Microquasars are thought to be downscaled analogs to
quasars exhibiting much smaller time scales and are therefore
potential laboratories for studying accretion and relativistic jets
near black holes (Mirabel et al. 1992). The most prominent object of
this type is the galactic X-ray binary system GRS 1915+105 which shows
dramatic variability (Greiner et al. 1996) in light curve,
quasi-periodic oscillations, phase lags and coherence behavior (Morgan
et al. 1997; Muno et al. 2001). GRS 1915+105 is the most energetic object
known in our galaxy with a luminosity of
1038 erg/s in the
low state.
GRS 1915+105 harbors a black hole of 14
(Greiner et al. 2001b) making
it the most massive stellar black hole known. Being in the galactic
plane and at a distance of
12 kpc the source suffers a very high
extinction in the optical band of the order of 25-30 mag (Greiner
et al. 1994; Chaty et al. 1996). Greiner et al. (2001a) found the donor to be an K-M III late type giant as determined from absorption line
measurements in the near infrared.
Black hole transients generally exhibit different states of intensity and spectrum: the high/soft state with a prominent disk component and a weak or negligible steep power law tail, the low/hard state with negligible accretion disk and flat power law. In the intermediate and very high state both the accretion disk and a power law are seen.
Originally discovered by Granat (Castro-Tirado et al. 1992), GRS 1915+105 was extensively monitored by the RXTE since 1996 and a number of investigations of these data have been published (e.g. Morgan et al. 1997; Belloni et al. 1997, 2000; Muno et al. 1999, 2000).
The source evades simple classification, although it seems to spend
most of the time in the very high state. Several attempts to
categorize the behavior of GRS 1915+105 have been made in the past. Belloni
et al. (2000) defined 12 different states depending on light curve
variation and hardness colors. One of these states, the so-called
-state is characterized by a lack of obvious variations in the
light curve and spectrum and is associated with continuous radio
emission of differing strength.
-states correspond to the
low/hard state of GRS 1915+105, exhibiting relatively low flux from the
accretion disk and showing a hard power law tail. Theses states
resemble long variants of the lulls in
-states, when part of
the inner accretion disk is proposed to be absent. Depending on the
strength of the radio emission, differing phase lag and power spectrum
behavior are seen. During the so-called radio quiet
-states the
phase lag of the hard X-ray photons to the soft X-ray photons is
partly negative, whereas it is always positive when the radio emission
is strong (Muno et al. 2001). Also, the frequency of the 0.5-10 Hz
QPO decreases with increasing radio flux.
In this paper we present a comprehensive study of the X-ray spectral
behavior of GRS 1915+105 in the
-state with over nearly four years of
observation with RXTE. The variation of the spectral properties
is then compared to those of the radio emission. Valuable information
about the disk geometry and properties of the compact object is
derived.
For our investigation of GRS 1915+105 we used public RXTE data from
November 1996 to September 2000 provided by HEASARC. We selected
Proportional Counter Array (PCA) and High Energy X-Ray Timing
Experiment (HEXTE) data of 139
-state observations from 89
different days. The selection of the datasets was based on the
-states defined by Belloni et al. (2000) and on PCA light curves
provided by
E. H. Morgan
. Data
with more than 30 min away from the South Atlantic Anomaly and
without X1908+075 in the HEXTE
background
were selected. The lack of obvious variability in the light curve and
its spectral hardness allows long continuous exposure times and
therefore high signal-to-noise ratios. Typically, the spectra during
this state show blackbody emission arising from an accretion disk and
a dominating power law hard energy tail.
We reduced the RXTE data using the standard reduction script REX
included in the HEAsoft5.04 package. We restricted the analysis to
Standard 2 binned data of all layers of PCU0 of the PCA from
3-25 keV and HEXTE cluster 0 from 20-190 keV only. We used PCARSP
7.10 to produce a particular response matrix for each PCU0 dataset. In
order to account for part of the uncertainties in the PCU0 instrument
we added a systematic error of 1% as recommended
(Remillard
).
The X-ray spectral fitting was done using XSPEC 11.0 (Arnaud 1996). A
consistent model should fit all RXTE
-state
data. Therefore we tested several X-ray radiation models and finally
selected a model consisting of (i) photoelectric absorption (WABS;
Balucinska-Church & McGammon 1992), (ii) a spectrum from an accretion
disk consisting of multiple blackbody components (DISKBB) and (iii) a
power law spectrum reflected from an ionized relativistic accretion
disk (REFSCH; Fabian et al. 1989; Magdziarz & Zdziarski 1995).
Several attempts in the past to fit GRS 1915+105 spectra have shown
complicated residuals in the soft X-ray band suggesting the existence
of emission and/or absorption features near 6.4 keV, the energy of
the Fe K
line (Kotani et al. 2000). Because of the low spectral
resolution of
1 keV of the PCA at this energy an additional
line fit gives no meaningful results. Therefore we ignored all energy
bins from 4.5 to 8.5 keV while fitting our model. We fixed the
hydrogen column density at
cm-2 as
determined by Greiner et al. (1994) with ROSAT. It has been shown
that the disk blackbody + power law assumption strongly overestimates
the flux at lower energies compared to the thermal and non-thermal
comptonization models used by Vilhu et al. (2001) and Zdziarski et al.
(2001), respectively. This explains the smaller
(
cm-2) values they found for
-state
observations. However, because we ignored the energy bins from 4.5 to
8.5 keV the amount of data bins needed to adjust the hydrogen column
density was to small to let
be a free parameter.
We fixed 9 of the combined 16 model parameters
(Table 1), leaving free the accretion disk temperature,
,
and relative normalization,
,
the
power law photon index,
,
with
and
relative normalization,
(photons/keV/cm2/s at 1 keV),
the reflection index, R, the ionization parameter,
,
and a
factor to account for the relative normalization between PCU0 and
HEXTE cluster 0. If not stated otherwise, we plot 1
errors for
each parameter of interest.
| fix. parameter | value |
| N |
|
|---|---|
| cutoff energy | no cutoff |
| redshift | 0 |
| Z>2 element abundances | 1 |
| iron abundance to abundance above | 1 |
| inclination angle | 70 |
| power law index for reflection emissivity | -2 |
| inner disk radius | 6 GM/c2 |
| outer disk radius | 1000 GM/c2 |
The appearance of steady radio emission in
-states suggested us
to search for correlations between RXTE data and radio data at
15 GHz (Ryle Telescope = RT) and 2.25 GHz (Green Bank
Interferometer = GBI). GRS 1915+105 shows variability at all frequencies on
time scales of seconds to hours. For a useful statement in
-states a suitably selection of corresponding datasets is
therefore required.
According to the general interpretation that the radio emission is synchrotron emission from ejected plasma in sporadic or continuous jets (Fender et al. 1995) the radio flux should reach a maximum 15 min after the actual ejection (Mirabel et al. 1997) and therefore after a possible determining X-ray event.
Figure 1 shows the 15 GHz radio flux,
,
from the RT for JD 2450898-2450913 together with two RXTE
observations. The RT observed GRS 1915+105 several times a day with
five minute exposures.
Because the radio exposure is much shorter than the X-ray exposure
(
hour) it is important to select radio data simultaneous with
the RXTE observations. 37 of the 139 analyzed RXTE
observations have simultaneous RT data and 9 have simultaneous
GBI data.
Still, the selection of simultaneous radio observations is
non-trivial. Occasionally the radio emission varies also during a
single
-state X-ray observation, whereas no variability is seen
in X-ray count rate and hardness ratio (5.2-60 keV/2-5.2 keV)
(Fig. 2). But the variation of
during a
-state
observation is negligible compared to the uncertainties of the
individual radio measurements and to the variation between different
RXTE observations. Therefore, the radio fluxes were averaged for
each RXTE observation.
| |
Figure 1:
Variation of the 15 GHz flux
(crosses) between JD 2450898 and JD 2450913. The dotted
lines mark the time of RXTE observations from 29.03.1998 (
|
| Open with DEXTER | |
![]() |
Figure 2:
15 GHz light curve ( RT, upper pannel), 2-60 keV PCA
light curve (middle pannel) and hardness ratio
|
| Open with DEXTER | |
![]() |
Figure 3:
Combined PCU0 and HEXTE cluster 0 spectra of the
|
| Open with DEXTER | |
The spectra of all 139 datasets have been fitted with the
above-described DISKBB+REFSCH model. Two example model spectra are
shown in Fig. 3. These spectra have the extreme values of
and R, revealing the spectral differences between different
-state observations.
The best fit parameters for the
-states of GRS 1915+105 are shown in
Fig. 4 and Table 2, respectively. The upper
two panels in Fig. 4 show the strong variability of the
source in the 1.5-12 keV (RXTE All Sky Monitor (ASM)) and
2.25 GHz (GBI) bands from JD 2450300-2451900. Irregular
outburst and relatively quiet phases alternate in X-ray and radio
without any obvious coupling.
The lower six panels of Fig. 4 show the PCU0 count rate
and the fit parameters of the
-state observations. The red and
black data points mark two groups of observations with different
(
)
behavior (see Fig. 5).
Between the
-states, the PCU0 count rate varies, with some
exceptions, by a factor of 2 at most. Only during the long continuous
-state at
and at JD
2451750
recurrent variability in the X-ray count rate is observed. The model
parameters of DISKBB+REFSCH are variable. The power law slope,
,
varies between 2.4 and 3.5 with a long-term periodicity of
590 days (Rau & Greiner 2002, in preparation). No correlation
of
and PCU0 count rate is seen, except around
and
.
The power law normalization
behaves similarly to the slope.
The reflection component, R, is variable between different
-state observations and shows a long-term variability similar to
the power law slope. It varies between 0 and 10 (upper limit of our model) with rather large
uncertainties when R>4. Except for five observations when R<1,
this rules out an isotropically sandwiching corona above the entire accretion disk.
![]() |
Figure 4:
1.5-12 keV ASM count rate (1st panel from top) and 2.25 GHz GBI
flux (2nd) from November 1996 to September 2000. ( GBI
was off-line some time before
|
| Open with DEXTER | |
| Obs-ID(1) | GD | JD | Exposure(2) | PCU0(3) | HEXTE0(4) |
|
|
R(8) | |
| (-2400000) | [s] | [cts/s] | [cts/s] | [keV] | |||||
| K-02-01 | 14.11.1996 | 50401.12 | 2688 | 995 | 70 | 2.60 +4.81-0.49 | 1.67 +4.90-1.60 | 2.89 +0.05-0.07 | 1.37 +0.55-0.86 |
| K-06-00 | 11.12.1996 | 50428.86 | 8656 | 587 | 60 | 2.50 +0.35-0.30 | 2.52 +1.46-1.45 | 2.69 +0.02-0.02 | 2.40 +0.40-0.44 |
| K-07-00 | 19.12.1996 | 50436.74 | 8992 | 598 | 61 | 2.21 +0.36-0.21 | 3.42 +3.25-2.32 | 2.70 +0.01-0.03 | 2.25 +0.20-0.56 |
| . | . | . | . | . | . | . | . | . | . |
| . | . | . | . | . | . | . | . | . | . |
| P-28-00 | 21.09.2000 | 51808.70 | 2384 | 706 | 69 | 2.70 +1.39-0.73 | 1.14 +8.85-1.01 | 2.84 +0.03-0.06 | 3.11 +1.42-2.64 |
| P-28-01 | 21.09.2000 | 51808.77 | 2400 | 651 | 66 | 3.15 +3.59-1.05 | 0.38 +0.21-0.37 | 2.83 +0.08-0.04 | 2.74 +0.22-1.29 |
| P-28-02 | 21.09.2000 | 51808.83 | 2336 | 645 | 44 | 2.00 +1.72-0.28 | 6.65 +23.43-6.47 | 2.87 +0.05-0.03 | 3.24 +0.24-2.29 |
| Obs-ID |
|
c(11) | ||
| [Pho/keV/cm2/s] | ||||
| K-02-01 | 43.5 +2.0-2.6 | 4999 +1-3521 | 0.82 | 0.69 |
| K-06-00 | 15.6 +1.2-0.3 | 5000 +0-4421 | 0.78 | 1.00 |
| K-07-00 | 16.3 +1.1-0.4 | 996 +4004-338 | 0.79 | 1.01 |
| . | . | . | . | . |
| . | . | . | . | . |
| P-28-00 | 34.1 +0.7-5.0 | 2902 +2098-2902 | 0.82 | 0.94 |
| P-28-01 | 29.8 +0.1-2.1 | 462 +4538-461 | 0.81 | 0.84 |
| P-28-02 | 27.2 +2.8-0.8 | 34 +2900-34 | 0.52 | 0.77 |
The significance of the accretion disk component varies through the
-states. For some observations the disk is more or less absent
(e.g. JD 50400-50600; relatively low disk temperature and large
inner disk radius (small normalization)), for other observations the
disk component provides a non-negligible contribution to the X-ray
flux (e.g. JD 51500-51600; high temperature and small inner
radius). In order to be consistent and to fit all
-state
observations with the same model, we included the DISKBB component in
all of our fits, although it could be excluded for several observations.
![]() |
Figure 5:
Power law slope, |
| Open with DEXTER | |
The accretion disk component shows a variable disk temperature of
1-4 keV. Sometimes, large uncertainties due to the small
contribution of the DISKBB component to the total flux are seen. The
inner disk radius, which can be determined from the disk
normalization, varies between 1 and 20 km. For a non-rotating black
hole of mass 14
(as measured for GRS 1915+105, Greiner et al. 2001b)
the Schwarzschild radius is
40 km. It is known that the DISKBB
model underestimates the inner disk radius by a factor of 1.7-3 due
to Doppler blurring and gravitational redshift (Merloni et al.
2000). Also the neglect of comptonization in the surface layers of the
disk leads to unphysical values when using then DISKBB model
(Zdziarski et al. 2001). But even a maximally rotating black hole
(where the inner disk radius reaches the Schwarzschild radius) cannot
account for the majority of the parameter values. This problem has to
be kept in mind when discussing the absolute values of the parameters.
Another free physical parameter of the REFSCH model is the ionization
parameter,
.
It has huge uncertainties because no Fe K
line could be fitted but has a negligible influence on the hard X-ray
continuum and our model parameters. Therefore we will not plot or
discuss
further. Theoretically, the reflection in the hard
spectrum should depend on the ionization parameter because higher
ionization means lower absorption and therefore higher reflection
probability. But the REFSCH model includes a simple 1-ionization zone
model only, which is very unlikely to be present in the disk and does
not show any dependence of reflection and ionization at all.
The model fits reveal an increasing power law normalization,
,
with a steepening power law component (Fig. 5) suggesting a pivoting behavior. Two branches with different slopes are seen in the correlation. No correlation of the power law slope with the X-ray count rate in ASM and/or PCU is observed.
The strength of the correlation can be tested using a Spearman
rank-order correlation test (Press et al. 1992). Both branches show
strong correlations (steeper:
,
flatter:
). (Note, this statistical method does not take into
account the particular uncertainties of the data points.)
The best descriptions of the correlations are
functions of the type
| (1) |
The pivoting has been tested in more detail by schematically plotting the obtained power law spectra. The two different branches consequently show different pivoting behavior (Fig. 6). Whereas the upper branch from Fig. 5 shows a pivoting energy of 4-8 keV (gray), the lower branch pivots at around 20-30 keV (black).
![]() |
Figure 6: Schematic power law spectra for the upper (top) and lower (bottom) branch of Fig. 5, the red and black points from Fig. 4, respectively, artificially offset by a factor of 100. The pivoting happens at different energies (4-8 keV & 20-30 keV). |
| Open with DEXTER | |
Due to the known uncertainties of the PCA below 6 keV (Jahoda et al. 1996) and the low energy resolution, stringent conclusions about the behavior of an iron line in the RXTE spectra at 6.4-7 keV not can be drawn. For the analysis of the reflection, therefore, only the continuum radiation can be used, manifesting itself in the reflection hump from 10-30 keV.
![]() |
Figure 7:
R( |
| Open with DEXTER | |
The reflection amplitude in GRS 1915+105 as a function of
is shown
in Fig. 7 together with that of three black hole
candidates (Cyg X-1, GX 339-4 and GS 1354-644; Gilfanov et al.
2000). Although R shows large uncertainties for the
-states of
GRS 1915+105, a similar correlation is seen as in the other X-ray
binaries. The steeper the power law component, the higher is the
reflected fraction of photons. Figure 7 contains the
confidence intervalls of R for three observations to clarify the
influence of the model on the R(
)-correlation. Although the
contours are elongated similiar to the correlation (higher
has higher R), they are intrinsically steeper compared to the
overall correlation. Thus, the existence of the
R(
)-correlation in the
-states of GRS 1915+105 is no artifact
of the model.
A group of observations at high
(3.3-3.5) and small R(0-2) behaves remarkably differently. These points belong to five
datasets from JD 2451497 and JD 2451513 showing very high PCU0 count
rates (1300-2100 cts/s) and are short duration
-states between
different high variability states.
The correlation of R and
is tested using a Spearman
rank-order correlation test. The correlation is distinct
(
)
but not strong.
For a quantitative description of the correlation, a phenomenological
function was fitted to the data. The best-fitting model is a power law
| (2) |
No obvious correlation between the soft X-ray component and the radio emission is seen. Neither the soft X-ray flux, the temperature of the accretion disk nor the inner disk radius show a correlation with the 15 GHz radio flux from RT and/or GBI .
An unexpected result is found when plotting
vs.
(Fig. 8). The power law slope correlates positively
with the radio flux at 2.25 GHz and 15 GHz. Observations with high
radio emission show a softer X-ray spectrum (steeper power law
component). Note that no bimodality in the radio emission
exists. Instead, a continuous spread is observed. A strict separation
of radio loud and radio quiet
-states, as done before
(e.g. Muno et al. 2001; Trudolyubov 2001) seems therefore
unsubstantiated.
The correlation of
and FR is most obvious for simultaneous
(
hrs) observations. Datasets with radio observations
5 hrs offset also show the correlation (Fig. 8
upper panel) but with significant scattering. Usually, no statement
about the variability state before and after the RXTE exposure
can be made, because the source may have had several state
alterations. Thus, non-simultaneity is the likely reason for the
strong scattering of the correlation for
hrs in
comparison to
hrs.
![]() |
Figure 8:
|
| Open with DEXTER | |
To determine the strength of the correlation a Spearman rank-order
correlation test was made for the RT data with
h. The correlation is strong (
). No test was made
for GBI due to the small number of simultaneous GBI/RXTE datasets.
The next step was to fit a phenomenological function (dotted line in
Fig. 8) to the correlation. The best model functions
are of the form
| (3) |
Besides the
(
)-correlation we also investigated the
relation between
and the X-ray flux of the power law
component. This introduces some freedom as to at which lower energy
the power law is "chosen''. Figure 9 shows the X-ray flux
from 20-200 keV and 1-200 keV, respectively, over the radio flux.
![]() |
Figure 9: 20-200 keV (left) and 1-200 keV (right) X-ray flux shown over the 15 GHz radio flux ( RT). The 20-200 keV flux has a negative correlation with the radio flux, whereas the 1-200 keV flux does not show such a correlation. |
| Open with DEXTER | |
With increasing radio emission, the X-ray flux in the power law
component between 20-200 keV decreases (neglecting reflection). The
total 1-200 keV X-ray flux in the power law component is more or less
constant. No correlation with the radio flux is seen. This is because
the flux in the 1-20 keV range by far dominates, thus washing out
the correlation. Muno et al. (1999) did not find a correlation of
(50-100 keV) with the 15 GHz radio flux. But they used a model
consisting only of DISKBB and a broken power law (BKNPO) for the X-ray
spectra and RT data with
12 hrs offset.
A Spearman rank-order correlation test for
(20-200) gives
,
thus showing a strong anti-correlation. For
(1-200 keV) the test results in
and ratifies the
lack of a correlation between the 1-200 keV non-thermal X-ray flux
and
.
| correlation between: | type(1) | ref. | |
|
|
+ | ||
| + | |||
| R | + | ||
| - | |||
|
|
- | Muno et al. (2001) | |
|
|
|
+ | Muno et al. (1999) |
|
(1):
direction of correlation ("+" = positive, "-" = negative),
(2): frequency of the 0.5-10 Hz QPO. |
Table 3 summarizes the correlations found here plus 2 relevant correlations found by Muno et al. (1999, 2000). In the following we will discuss implications of these correlations.
The origin of the hard power law component in AGN and X-ray binaries is still under discussion. Three main models exist, where the hard photons originate by inverse Compton scattering of soft disk photons on hot electrons. The main difference between these models is the distribution of the electrons. They can be thermal (Maxwellian), non-thermal (power law like) or free falling from the last stable orbit onto the event horizon of the black hole. The high energy spectrum presents an important test for the distinction between these models. Therefore, the behavior of the power law component in GRS 1915+105 is crucial for the understanding of the electron distribution and for the geometry of the system.
It is still unclear whether or not the hard spectrum of GRS 1915+105 in the
-state extends with or without cutoff up to MeV energies. Due to
the rapid variability of the source and the long required exposure
times it is nearly impossible to get sufficiently accurate
-spectra with the present high energy satellites. For example,
the OSSE spectrum from 14th-20th of May 1997 reveals a power law tail
up to 1 MeV without a cutoff (Iyudin 2000; Zdziarski et al. 2001). But
as several contemporaneous RXTE observations show (which cover
3% of the OSSE exposure), during the observation GRS 1915+105 went
through several
/
state transitions and the PCA photon
count rate from 2-40 keV varied between 5 and 25 kcts/s. Thus, the
OSSE spectrum is the sum of the spectra of the
and
(and possible other) states, with an unkown fraction coming from the
-state. It is therefore premature to unequivocally relate the
hard MeV tail to the
-state.
| |
Figure 10: Bulk motion geometry. Electrons become spherical free-falling between the last stable orbit and the event horizon due to shocking. The photons are comptonized on these accelerated electrons. |
| Open with DEXTER | |
One model for the comptonization is the bulk motion comptonization
(Blandford & Payne 1981; Chakrabarti & Titarchuk 1995)
(Fig. 10). The accretion stream passes the transition
radius,
,
at the last stable orbit around the black hole and
gets shocked. The matter falls spherical onto the event horizon. In
front of the shock, the stream is similar to an optically thick,
geometrically thin accretion disk (as required by the strong observed
blackbody component).
The free falling electrons are accelerated up to the speed of light
and inverse Compton scattering of soft photons on the electrons
provides the observed power law in the hard spectrum. The model
predicts a cutoff depending on the mass accretion rate (Ebisawa et al.
1996) due to inverse Compton scattering and Compton recoil. As stated
above, the existence/absence of a cutoff during the
-states is
still not secured, wherefore not statement about the model can be made
from that issue.
Two major problems exist for the description of
-state data by
the BMC model. The RXTE spectra of GRS 1915+105 partly show a strong
reflection component. This is not compatible with the geometry of a
thin disk outside a central, spherical accretion stream. The accretion
stream is accelerated away from the inner disk edge onto the event
horizon of the black hole. Therefore only a small part of the soft
photons are comptonized back into the disk plane at
.
In
conflict to the observations, no strong reflection component can
emerge. A possible explanation for the observed reflection can then be
a partly covering of the disk by an absorbing medium with
(Zdziarski 2000). Recent Chandra X-ray spectra indeed
show a high X-ray column density and abundance excesses for Si and Fe,
which may be related to material that is associated with the immediate
environment of the source (Lee et al. 2002).
The other problem for the BMC model is the time lag of hard and soft
X-ray photons. Muno et al. (2001) found a phase lag in the 5 Hz QPO
(which are connected to the hard X-ray photons) of 0.5. This
corresponds to a delay of 0.1 s of the hard photons to the soft
photons. The expected delay for scattering on a convergent electron
stream inside
is
0.2 ms, much smaller than
observed. The phase lag is not constant in time and frequency and
seems to vary with the radio flux.
It is obvious that the BMC model not can explain the hard power law
component in the
-states of GRS 1915+105. There are indications that the
BMC model fits other variability states. Shrader & Titarchuk (1998)
found good spectral agreement of the BMC model and RXTE
observations in the
-state. Note, however, that they completely
neglected the time delay of hard and soft photons and the strong
variability of the
-state for their X-ray spectral modelling.
Another possible model for the origin of the hard X-ray component is the disk corona geometry (Fig. 11) (e.g. Haardt & Maraschi 1993; Svensson & Zdziarski 1994). Soft disk photons are inverse Compton scattered on hot (>20 keV) electrons located above the disk.
| |
Figure 11: Disk corona geometry. The disk photons are comptonized in an electron distribution located above the accretion disk. (Note, the shape and size of the electron distribution in this plot has no real physical or observational meaning.) Part of the photons are scattered back into the disk plane and are reflected or reprocessed. The photons reaching the detector are compound of the soft disk photons, the comptonized and the reflected photons. |
| Open with DEXTER | |
Simultaneous CGRO/OSSE and X-ray observations of X-ray binaries with black holes showed that the high energy continuum in the low/hard state originates most likely due to thermal comptonization. The predicted cutoff in the hard X-ray spectrum is seen in Cyg X-1 (Gierlinski et al. 1997), in the LMXB GX 339-4 (Zdziarski et al. 1998) and some other X-ray transients (Grove et al. 1998).
Also the high/soft state in X-ray binaries can be explained by thermal
comptonization but requires high electron temperatures in some systems
to explain the observed unbroken power law up to MeV. In addition
is required in order to keep the spectrum soft (Zdziarski
2000).
The power law slope is a measure of the Compton amplification,
.
is the dissipated energy
in the corona and
the energy of the incoming
photons. Table 4 shows how A depends on
with
![]() |
(4) |
|
|
y | A | |
| 2.4 | 84% | 1.33 | 1.84 |
| 3.0 | 22% | 3.65 | 1.22 |
| 3.5 | 9% | 7.30 | 1.09 |
The minimum energy which a photon receives when passing a thermal
electron distribution depends on the electron temperature
![]() |
(5) |
The radiation processes and geometry can also be described by comptonization of soft disk photons on a hybrid (thermal and non-thermal) electron distribution above an optical thick accretion disk (as in Fig. 11) (Coppi 1992). Part of the electrons in the thermal corona are accelerated, possibly in reconnection events. The disk photons are comptonized on these electrons forming the observed power law component in the hard X-ray spectrum. The energy of the non-thermal electrons is partially transfered to the thermal electrons due to Coulomb scattering. This leads to heating of the thermal electrons above the Compton temperature. Then the thermal electrons play an important part in the up-scattering of the soft disk photons, too.
The compactness (ratio of luminosity and size) is an important parameter of the coronal plasma (Coppi 1999). A large compactness leads to electron positron pair creation due to photon photon collisions resulting in a pair annihilation line at 511 keV in the spectra. When the compactness is low, the loss of energy of the electrons due to Coulomb scattering dominates the loss due to inverse Compton scattering. The ratio of thermal to non-thermal electrons increases and the plasma becomes thermally dominated. Together with the resulting break in the distribution of the non-thermal electrons this results in a cutoff of the photon spectrum.
The rate at which non-thermal electrons appear in the hybrid electron
distribution can be written as a power law (Gierlinski et al. 1999)
![]() |
(6) |
| (7) |
In GRS 1915+105, the observed continuous distribution of
between 2.4
and 3.5 is evidence for a variability in the power law component and
thus in the composition and locus of the electron distribution. From
the observation of the reflection and the 0.5-10 Hz QPO it is
conspicuous that the electron distribution must have a small (relative
to the system) spatial size. As a small compactness, e.g. large
plasma volume, results in a spectral cutoff, for the
-states of
GRS 1915+105 a hybrid electron distribution not can be ruled out. It is
difficult to distinguish spectroscopically between a hybrid and a
thermal plasma if the hybrid plasma is thermally dominated. The only
possibility is the search for the predicted photon excess at high
energies and for the annihilation line (Coppi 1999). For a better
understanding, enhanced high energy detectors, such as on INTEGRAL,
are required. The analysis of such spectra also needs further
developed physical models, instead of a simple power law.
The investigation of the RXTE data of GRS 1915+105 with the
DISKBB+REFSCH model does not allow a definite conclusion about the
origin of the hard spectral component. Only the BMC model is ruled out
for the
-states. Though a clear distinction between a thermal
and a hybrid not can be proposed, it seems clear that the thermal
electrons dominate.
It has been shown above (Fig. 5) that two different states of
comptonization with two different pivoting energies are observed. Very
recently (Zdziarski et al. 2002) pivoting was found also in Cyg X-1. A
varying amount of soft seed photons undergoing the comptonization
process probably accounts for the pivoting of the spectrum. One
possible explanation is a constant disk blackbody and an
expanding/contracting hot coronal plasma. A change in the size of the
corona leads to a change in the fraction of intercepted soft photons.
Assuming constant optical plasma depth,
,
if the soft
luminosity is high, the electrons in the corona are cooled efficiently
and the spectrum of the comptonized photons is soft. On the other
hand, if the soft luminosity decreases, the temperature of the corona
increases and the spectrum hardens). The ocurrance of two differing
pivoting branches suggests the existence of two states with different
and/or different electron composition (thermal/non-thermal).
The reflection of X-rays on an accretion disk manifests itself in two
specific spectral features: (i) the characteristic emission line
spectrum (mainly the K
lines of the most common metals) with
the 6.4 keV iron emission lines as the strongest and (ii) a
characteristic hump arises above 10 keV because of the energy
dependence of the cross section of absorption and Compton scattering.
Observations with ASCA in different variability states indicate
a variable iron absorption or emission in GRS 1915+105 (Ebisawa et al.
1997). The iron emission at 6.4 keV dominates the spectra from
25.10.1996 and 25.04.1997 whereas a distinct absorption at 7 keV and
no emission was seen on 27.09.1994 and 26.04.1995. During the
25.04.1997 observation GRS 1915+105 switched between
- and
-states (based on RXTE data). Therefore, it is not clear
whether the emission line originates during the
- or during the
-interval. No statement can be made about the states during
the other ASCA observations due to lack of corresponding RXTE observations.
Recent Chandra data of GRS 1915+105 obtained in the low hard state
revealed neutral K absorption edges, ionized resonance absorption from
Fe (XXV, XXVI) and possible emission from neutral Fe K
and
ionized Fe XXV (Lee et al. 2002), suggesting conditions favorable for
reflection.
Investigations of Seyfert galaxies and X-ray binaries in the low/hard state show a correlation between the slope of the hard X-ray component and the reflection. Incoming radiation with a steeper slope is more reflected than radiation with a flatter slope.
The reflection is stronger in black hole candidates than in Seyfert
galaxies for a given
(Zdziarski 1999). This can be explained
using an optical thick accretion disk. The different masses of the
black holes imply different maximum energies of the blackbody photons
(Svensson 1996). The hard X-ray photons from the corona can ionize the
upper layers of the accretion disk and an ionized layer above neutral
matter evolves. If
is large the heating of the disk is small
and the cold layer of the disk lies near the disk surface. With
decreasing power law slope the temperature of the upper layer and
therefore the ionization increases, the absorption is smaller, more
photons are scattered and the characteristic reflection hump is
suppressed, sufficient to explain the observed
R(
)-correlation.
The correlation in Seyfert galaxies is generally steeper than that
observed in X-ray binaries (Zdziarski 1999). In Seyferts the
correlation has been explained with a model consisting of a static,
thermal corona above a neutral reflector (Svensson 1996). The
existence of the R(
)-correlation implies a feedback where the
existence of the reflecting matter influences the hardness of the
X-ray spectrum (Böttcher et al. 1998; Zdziarski et al. 1999). Assuming
that the cold matter (the accretion disk) emits soft photons which
become seed photons of the comptonization, than with increasing solid
angle of the reflector (here the accretion disk) the flux of soft
photons increases and the cooling rate of the hot corona
increases. For a thermal plasma the resulting power law component
steepens with increasing cooling rate.
On the other hand, models which are based on non-thermal electrons should not show a dependence of spectral hardness and cooling rate (Lightman & Zdziarski 1987; Zdziarski et al. 1999). Models in which the seed photons are intrinsically produced in the hot plasma (e.g. synchrotron radiation) not can reproduce the observed correlation. The power law slope is then independent of the reflection.
The
-states in GRS 1915+105 show a similiar correlation as seen in
Seyfert galaxies and X-ray binaries (Zdziarski et al. 1999). The
reflection amplitude is strongly effected by the disk model. The
higher the disk temperature, the larger the reflection amplitude. The
disk temperatures in X-ray binaries are higher than in Seyfert
galaxies. In GRS 1915+105 the temperature of the disk is still higher,
explaining why a majority of the GRS 1915+105 data lies above the correlation
function (showing slightly higher R for a given
)
in
Fig. 7.
The reflection amplitude of R> 1 strongly suggests an anisotrophic
inverse Compton process. The dominant fraction of the up-scattered
soft disk photons is directed back into the plane of the accretion
disk. This produces the observed large reflection amplitude in
contrast to an isotrophic scattering where
.
Observations of GRS 1915+105 show a strongly varying
,
therefore the
corona can hardly be static above the reflector. In GRS 1915+105,
increases with increasing radio emission. Higher radio emission may
imply more outflow away from the disk with
.
Thus, the
higher the mass outflow, the lower the electron temperature which is
required to produce the observed amount of comptonization. On the
other hand, with a dynamical, thermal corona the spectrum hardens with
increasing
due to relativistic aberration, giving a flatter
power law with higher outflow velocity (Malzac et al. 2001).
The observation in GRS 1915+105 predicts an at least partly thermal source of the hard X-ray photons and an important contribution of the cooling due to the soft photons. This allows both a thermal and a hybrid electron distribution.
It is more difficult to interpret the behavior of observations with
high
and small R in Fig. 7. Unlike in a
thermal plasma, the Compton scattering in a non-thermal plasma does
not depend on the cooling rate of the soft photons, instead it depends
on the slope of the electron distribution (Poutanen & Coppi 1998). A
change in the properties of the plasma can provide the differing
behavior.
For a better determination of the correlation and for conclusions about the distinct system components, more complex ionization models are needed. They partly exist (Done & Nayakshin 2000) but require too much computing power to be useful for the analysis of real spectra at this stage.
The lack of a correlation of the radio flux with any of the disk
parameters is somewhat unexpected since the (short-duration) low/hard
states have earlier been related to jet ejections with radio emission
by synchrotron radiation of the ejected plasma (Pooley & Fender
1997). Since IR (Eikenberry et al. 2000) and radio observations (Pooley
& Fender 1997; Fender & Pooley 2000) have led to the conclusion that
GRS 1915+105 shows jets on various scales, or even has a continuous
distribution of jet strength, one would have expected that also in
-states the radio flux correlates with the disk
temperature/emissivity.
The matter in the sporadic, relativistic jets may originate from
disruption of the inner part of the accretion disk during
-states (Mirabel et al. 1998). In
-states no correlation
of the radio emission and the inner disk radius is observed. Instead,
the hard spectral component is correlated with the radio flux.
Muno et al. (1999) found a positive correlation of the accretion disk
temperature and the frequency of the 0.5-10 Hz quasi periodic
oscillations, and of the radio flux with the QPO frequency (Muno et al.
2001). With increasing
and decreasing radio emission the QPO
frequency increases. This predicts a negative correlation of
and radio flux. The fact that nothing like this is found in the
present analysis can have several reasons. As mentioned before, Muno
et al. (1999) used the standard model DISKBB+BKNPO for the analysis of
the RXTE spectra. This seems not suitable for the
-states
of GRS 1915+105. Accordingly, the disk temperature and therefore the
correlation of the QPO frequency and the disk temperature have to be
interpreted carefully.
It is generally assumed that both the radio jets (Fendt & Greiner
2001; Fender 2001) and the hard spectral component originate near the
black hole. With increasing outflow of matter (therefore increasing
radio emission) the interaction with the electron distribution becomes
significant. Similar to the interpretation of the pivoting behavior
of the X-ray spectra, an increasing outflow of matter implies an
increasing size of the scattering medium and therefore an increasing
amount of intercepted soft seed photons. This leads to a lower plasma
temperature due to cooling and to a softer X-ray spectrum, which is
seen in the
(
)-correlation. Simultaneously, the
outflowing matter intermingles with the coronal matter and pushes it
away from the accretion disk. The formerly thermally dominated
electron distribution may become non-thermal dominated. This should
result in a shift of the cutoff of the power law component in the
-states with higher radio emission above the HEXTE range to
higher energies.
When the radio emission is low, the observed power law is supposed to
originate in the corona above the accretion disk. Assuming that the
radio quiet state lacks outflow of matter or the effect of the outflow
on the corona is negligible, the hard spectral component possibly
originates due to comptonization of disk photons in the corona. A
thermal electron distribution can then explain the X-ray spectra of
the radio quiet
-states.
The hard X-ray component in
-states therefore comes most likely
from soft disk photons which are inversely scattered on a thermal
dominated electron distribution (when radio flux is low) or on the
base if a continuous jet (when radio flux increases). The suggestion
of the base of the jet as source of the hard X-ray photons was already
made by Fender (2001). Only the increasing reflection with increasing
radio flux is still unresolved.
In conflict with the anti-correlation of
(20-200) and
in
GRS 1915+105 a linear correlation is observed on the LMXB GX 339-4 (Corbel et al.
2000). With increasing radio flux the X-ray flux increases. This is
true for the hard (20-100 keV) and for the soft X-ray flux
(2-12 keV) where GRS 1915+105 shows no correlation. The soft disk component
in GX 339-4 is negligible above 2 keV and the power law component alone
quantifies the flux. Although GRS 1915+105 shows an equal
R(
)-correlation as GX 339-4 the
(
)- and
(
)-correlation is opposite. That is remarkable because the
R(
)-correlation suggests a similar structure near the black
hole (especially the corona) in GRS 1915+105 and GX 339-4. No such correlation of
and
is observed in Cyg X-1 (Brocksopp et al. 1999).
In principle the hard X-ray photons can originate in the continuous jet itself. But following the equations of Marscher (1983) the hard X-ray photons from self comptonization of a compact syncrotron jet are negligible with regard to the coronal contribution. An analogous result was found for Cyg X-1 and GX 339-4. The contribution of thermal X-ray bremsstrahlung in the jet is also too small to explain the observed hard X-ray luminosities (Memola et al. 2002).
The analysis of 139 RXTE observations of GRS 1915+105 in the
-state revealed variable components and parameters of the X-ray
spectra. The structure of the hard X-ray power law depends e.g. on
the radio flux. Further, a two-branch correlation of the power law
slope and the power law normalization was found. Theses branches show
different pivoting behavior.
The most probable geometry is that of a hot corona above the accretion
disk. A bulk motion comptonization seems to be ruled out for the
-states because of the time lags of hard and soft X-ray photons.
The continuous outflow of matter is not correlated with the accretion disk parameters as measured at X-rays. Neither the disk temperature nor the inner disk radius are connected with the radio flux. Therefore, the outflowing matter should not be provided by a disruption of the inner part of the accretion disk, as thought to be the case for the sporadic jets in GRS 1915+105, and a positive correlation between radio emission and power law slope is found. Because both the base of the jet and the corona are probably located near the black hole the correlation indicates an interaction of both structures. Whereas during low radio emission the hard X-ray component originates due to comptonization on a thermally-dominated corona, in radio loud states the comptonization should appear in the outflowing matter.
GRS 1915+105 shows a positive R(
)-correlation as seen in other
X-ray binaries and AGN. The large reflection amplitude suggests a
highly anisotrophic inverse Compton scattering with the dominant part
of the soft photons being scattered back into the disk plane.
Acknowledgements
The authors thank E. H. Morgan (MIT) and G. G. Pooley (MRAO) for providing the ASM and RT data and W. A. Heindl (RXTE) and K. A. Arnaud (XSPEC) for their comments. We thank the referee, A. Zdziarski, for the thorough reading and several suggestions which improved the presentation and T. Belloni for the constructive criticism. The Ryle Telescope is supported by PPARC. The Green Bank Interferometer is a facility of the National Science Foundation operated by the NRAO in support of NASA High Energy Astrophysics programs.