Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://uranus.sai.msu.ru/~seif/apj_706_2_1463.pdf Äàòà èçìåíåíèÿ: Fri Nov 13 16:07:34 2009 Äàòà èíäåêñèðîâàíèÿ: Mon Oct 1 19:54:08 2012 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: ï ï ð ð ð ð ð ð ï ï ï ï ï ï ï ï

The Astrophysical Journal, 706:1463­1483, 2009 December 1
C

doi:10.1088/0004-637X/706/2/1463

2009. The American Astronomical Society. All rights reserved. Printed in the U.S.A.

DISCOVERY OF PHOTON INDEX SATURATION IN THE BLACK HOLE BINARY GRS 1915+105
Lev Titarchuk1
1

,2 ,3 ,4 ,5

and Elena Seifina6

6

` Dipartimento di Fisica, Universita di Ferrara, Via Saragat 1, I-44100 Ferrara, Italy; titarchuk@fe.infn.it 2 ICRANET Piazzale d. Repubblica 10-12 65122 Pescara, Italy 3 George Mason University Fairfax, VA 22030, USA 4 US Naval Research Laboratory, Code 7655, Washington, DC 20375, USA; Lev.Titarchuk@nrl.navy.mil 5 Goddard Space Flight Center, NASA, Code 663, Greenbelt MD 20771, USA Moscow State University/Sternberg Astronomical Institute, Universitetsky Prospect 13, Moscow, 119992, Russia; seif@sai.msu.ru Received 2009 February 8; accepted 2009 October 19; published 2009 November 12

ABSTRACT We present a study of the correlations between spectral, timing properties, and mass accretion rate observed in X-rays from the Galactic black hole (BH) binary GRS 1915+105 during the transition between hard and soft states. We analyze all transition episodes from this source observed with the Rossi X-ray Timing Explorer, coordinated with Ryle Radio Telescope observations. We show that broadband energy spectra of GRS 1915+105 during all these spectral states can be adequately presented by two bulk motion Comptonization (BMC) components: a hard component (BMC1, photon index 1 = 1.7­3.0) with turnover at high energies and soft thermal component (BMC2, 2 = 2.7­4.2) with characteristic color temperature 1 keV, and the redskewed iron-line (LAOR) component. We also present observable correlations between the index and the normalization of the disk "seed" component. The use of "seed" disk normalization, which is presumably proportional to mass accretion rate in the disk, is crucial to establish the index saturation effect during the transition to the soft state. We discovered the photon index saturation of the soft and hard spectral components at values of 4.2 and 3, respectively. We present a physical model which explains the index­seed photon normalization correlations. We argue that the index saturation effect of the hard component (BMC1) is due to the soft photon Comptonization in the converging inflow close to the BH and that of soft component is due to matter accumulation in the transition layer when mass accretion rate increases. Furthermore, we demonstrate a strong correlation between equivalent width of the iron line and radio flux in GRS 1915+105. In addition to our spectral model components we also find a strong feature of "blackbody (BB)-like" bump whose color temperature is about 4.5 keV in eight observations of the intermediate and soft states. We discuss a possible origin of this "BB-like" emission. Key words: accretion, accretion disks ­ black hole physics ­ radiation mechanisms: non-thermal ­ stars: individual (GRS 1915+105) Online-only material: color figures

1. INTRODUCTION The study of the characteristic changes in spectral and variability properties of X-ray binaries is proved to be a valuable source of information on the physics governing the accretion processes and on the fundamental parameters of black holes (BHs). The simultaneous study of the spectral and timing evolution of a BH source during a state transition has been a subject of many investigations (see references in a review by Remillard & McClintock 2006). Fender & Belloni (2004), hereafter FB04, introduced a classification of the spectral states in GRS 1915 +105 and studied the spectral state evolution. Using X-ray colors (hardness ratio) they introduced three spectral states. State A: in which the strong blackbody (BB)-like component of color temperature 1 keV dominates in the overall spectrum and little time variability is detected. State B: similar to state A but substantial red-noise variability on scales >1 s occurs in this state. State C: the spectra are harder than those in states A and B. Photon indices of the power-law components vary from 1.8 to 2.5. White­red-noise (WRN) variability on scales >1stakes place in this state. Furthermore, FB04 discussed the connection between states A, B, ­C observed in GRS 1915 with the three "canonical" states in black hole candidates (BHCs) also identified by their timing and spectral properties. 1463

At a low-luminosity state, the energy spectrum is dominated by a hard Comptonization component combined (convolved) with a weak thermal component. The spectrum of this low (luminosity) hard state (LHS) is presumably a result of Comptonization (upscattering) of soft photons, originated in a relatively weak accretion disk, off electrons of the hot ambient plasma (see, e.g., Sunyaev & Titarchuk 1980). Variability in LHS is high (fractional rms variability is up to 40%) and presented by a flat-top broken power-law (WRN) shape, accompanied by quasi-periodic oscillations (QPOs) in the range of 0.01­30 Hz, observed as narrow peaks in the power density spectrum (PDS). In high soft state (HSS), a photon spectrum is characterized by a prominent thermal component which is probably a signature of a strong emission coming from a geometrically thin accretion disk. A weak power-law component is also present at the level of not more than 20% of the total source flux. In the HSS, the flat-top variability ceases, QPOs disappear, and the PDS acquires a pure power-law shape. The total variability in HSS is usually about 5% fractional rms. The intermediate state (IS) is a transitional state between LHS and HSS. Note that in addition to LHS, IS, and HSS sometimes very soft state (VSS) is observed in which the BB component is dominant and the power-law component is either very weak or absent at all. The bolometric luminosity in VSS is a factor of 2­3 lower than that in HSS. FB04 concluded that probably all three states A, B, C of GRS 1915+105 are instances of something similar to the HSS/


1464

TITARCHUK & SEIFINA

Vol. 706

IS observed in other BHC systems, associated to the high accretion rate value for this source, although during the hardest intervals LHS might be sometimes reached. We come to the similar conclusions analyzing spectral and timing data from GRS 1915+105 obtained by the Rossi X-ray Timing Explorer (RXTE; see below). Close correlations of nearly periodic variability (QPOs) observed during low-hard and ISs with the photon index of the Comptonization spectral component have been reported in multiple state transitions observed from accreting BHs (see Vignarca et al. 2003; Shaposhnikov & Titarchuk 2006, 2007, 2009, hereafter V03, ST06, ST07, and ST09, respectively). The ubiquitous nature of these correlations suggests that the underlying physical processes which lead to the observed variability properties are closely tied to the corona; furthermore, they vary in a well-defined manner as the source makes a transition between spectral states. The fact that the same correlations are seen in many sources, which vary widely in both luminosity (presumably with mass accretion rate) and state, suggests that the physical conditions controlling the index and the low-frequency QPOs (LFQPOs) are characteristics of these sources. Moreover, they may be an universal property of all accreting compact systems, including neutron sources too (see Titarchuk, & Fiorito 2004, hereafter TF04 and Titarchuk & Shaposhnikov 2005). When a BH is in LHS, radio emission is also detected and a jet is either seen or inferred (Fender 2001). Several models are successful in reproducing the energy spectrum from the radio domain to the hard X-rays (see, e.g., Markoff et al. 2001; Vadawale et al. 2001; Corbel & Fender 2002; Markoff et al. 2003; and Giannios 2005). The multiplicity of models that can fit well the time average spectrum of Galactic BHs indicates that this alone is not enough to distinguish the most realistic one among them. X-ray timing features can be the key features to finally finding the common physical connection between the corona, the accretion disk, and the jet radio emission in BHs. There is a big debate in the literature on the origin of QPO frequencies (see, e.g., Remillard & McClintock 2006) and its connection with the radio emission. Migliari et al. (2005) reported on correlation between radio luminosity and X-ray timing features in X-ray binaries containing a number of low-magnetic-field neutron stars and one BH GX 339-4. They showed that in the LHS, radio luminosity is correlated with the LFQPO. Note that ST09 demonstrated that in the LHS of Galactic BHs, LFQPO changes by order of magnitude, from 0.2 to 2 Hz whereas the photon index has almost the same value of 1.5. Below we show that in GRS 1915+105 the photon index monotonically increases with LFQPO and disk mass accretion rate, although the radio luminosity does not correlate with LFQPO and X-ray luminosity in the whole range of spectral states, from low hard to high soft through ISs. Recently, Kylafis et al. (2008) suggested a model which explains how the QPO phenomenon is related to an appearance of radio flares (jets). Below (see Section 3) we present details of our observational study of the QPO connection with the X-ray and radio flaring activity in GRS 1915+105. In LHS and IS which we consider in our study, only a small part of the disk emission component is seen directly. The energy spectrum is dominated by a Comptonization component presented by a power law. To calculate the total normalization of the "seed" disk BB component we model the spectrum with a generic Comptonization model (BMC XSPEC model, see details in Titarchuk et al. 1997) which consistently convolves a disk BB with a Green's function of the Compton corona to

produce the Comptonization component. We argue that the disk emission normalization calculated using this approach produces a more accurate correlation with respect to the correlation with the direct disk component which was obtained using the additive model, multicolor disk plus power law (see, e.g., McClintock & Remillard 2006). This paper is a continuation of the study of index­QPO and index­seed photon normalization correlations in BH sources started in ST07 and ST09. Particularly, here we present a study of the index­seed photon normalization (disk flux) correlation observed from GRS 1915+105 when it evolves from LHS to HSS. The description of RXTE data set used is given in Section 2. We have analyzed a broader sample of state transitions from GRS 1915+105 and we found a diverse phenomenology for index evolution through a transition. In Section 3, we provide a detailed description of state transitions analyzed in this study. In Section 4, we discuss and interpret the results of our observational study. Specifically in Section 4, we consider the effect of the bulk motion Comptonization (BMC) in the inner part of the accretion flow on the index evolution during a state transition. Also we show that the index saturation effect is a direct consequence of the existence of this inner BM region and, therefore, can be considered as an observational signature of the converging flow (BH). Furthermore, in Section 4 we discuss the TF04 model and the Monte Carlo simulations by P. Laurent & L. Titarchuk (2009, in preparation) in which the observable index evolution with m has been already predicted. Conclusions follow in Section 5. 2. OBSERVATIONS AND DATA REDUCTION In the present Paper, we have used publicly available data of the RXTE Observatory obtained from 1997 January to 2006 April. In total, our study includes 107 observations made at different BH spectral states (LHS, IS, HSS) of the system. Data sets were selected to represent a complete rise-middledecay track of bright X-ray activity episodes behavior along bright radio flaring events (S15 GHz 250 mJy). Therefore, we have chosen powerful ( 250 all-sky monitor counts s-1 ) flaring episodes of GRS 1915+105 with a good coverage of simultaneous radio/X-ray observation. In the past, some of these data of spectral transitions in GRS 1915+105 were analyzed by Trudolyubov et al. (1999), Trudolyubov (2001), Muno et al. (1999), Reig et al. (2000), ST07, and Rodriguez et al. (2008) for the 1997­1998 and 2005­2006 transitions, respectively. Standard tasks of the LHEASOFT/FTOOLS 5.3 software package were utilized for data processing using methods recommended by RXTE Guest Observer Facility according to The RXTE Cook Book (http://heasarc.gsfc.nasa.gov/docs/ xte/recipes/cook_book.html). For spectral analysis, we used Proportional Counter Array (PCA) Standard 2 mode data, collected in the 3­20 keV energy range. The standard dead-time correction procedure has been applied to the data. To construct broadband spectra, data from High-Energy X-Ray Timing Experiment (HEXTE) detectors have also been used. We subtracted background corrected in off-source observations. Only HEXTE data in the 20­150 keV energy range were used for the spectral analysis in order to exclude the channels with largest uncertainties. The HEXTE data have been re-normalized based on the PCA. The data are available through the Goddard Space Flight Center (GSFC) public archive (http://heasarc.gsfc.nasa.gov). In Tables 1­6, we list groups of observations covering the complete dynamical range LHS­(IS)­HSS­(IS)­LHS of the source evolution during flaring events. We present here period ranges


No. 2, 2009

DISCOVERY OF PHOTON INDEX SATURATION IN THE BB BINARY GRS 1915+105

1465

Table 1 Best-fit Parameters of Spectral Analysis of PCA and HEXTE Observation of GRS 1915+105 in 3­150 keV Energy Range During Rise 1997 Transitiona Observational ID 20402-01-11-00 20402-01-12-00 20402-01-13-00 20402-01-16-00 20402-01-19-00 20402-01-20-00 20402-01-21-00 20402-01-21-01 20402-01-23-00 20402-01-24-00 20402-01-28-00 20402-01-29-00 20402-01-30-01 20402-01-31-00 20402-01-31-02 20402-01-33-00 20402-01-34-01 20402-01-35-00 20402-01-36-00 20402-01-37-01 20402-01-38-01 20402-01-39-02 20402-01-41-01 20402-01-41-02 20402-01-43-00 20402-01-45-00 MJD (day) 50462.06 50471.06 50477.87 50501.88 50517.04 50524.92 50533.83 50534.90 50548.45 50561.13 50586.68 50589.49 50596.20 50602.53 50604.61 50617.54 50621.80 50636.62 50639.62 50641.48 50649.42 50658.51 50679.30 50679.37 50688.24 50698.65 1 = 1 - 1 0.867(9) 0.844(7) 0.92(2) 0.83(1) 0.62(1) 0.85(1) 0.94(1) 0.95(1) 0.91(3) 0.72(3) 1.59(2) 0.95(2) 1.59(1) 1.79(2) 1.80(5) 1.99(6) 1.92(4) 1.85(6) 1.99(9) 1.90(9) 1.9(1) 2.0(1) 2.0(1) 2.00(9) 2.00(5) 1.90(1) log (A1 )
b

Nbmc1 , 2 L39 /d10 0.115(3) 0.1105(2) 0.116(1) 0.1014(5) 0.101(1) 0.092(2) 0.100(1) 0.101(1) 0.112(2) 0.116(5) 0.2001(2) 0.159(1) 0.200(8) 0.223 (5) 0.218(7) 0.22(4) 0.22(1) 0.24(4) 0.20(1) 0.245(9) 0.370(9) 0.26(2) 0.417(6) 0.459(2) 0.342(5) 0.231(2)

2 = 2 - 1 1.84(7) 1.9(2) 1.9(3) 2.0(4) 1.8(2) 1.9(3) 2.0(4) 2.0(2) 2.0(1) 2.0(3) 2.34(3) 2.0(2) 2.3(1) 2.6(1) 2.6(1) 3.2(2) 3.1(1) 3.10(9) 3.2(1) 3.09(9) 3.20(5) 3.2(3) 3.15(5) 3.20(1) 3.20(7) 3.20(2)

log (A2 )

b

Nbmc2 , 2 L39 /d10 0.009(3) 0.009(2) 0.002(5) 0.008(2) 0.002(1) 0.01(2) 0.009(5) 0.009(2) 0.01(3) 0.01(4) 0.109(1) 0.01(3) 0.110(2) 0.118(5) 0.127(6) 0.24(1) 0.17(1) 0.21(5) 0.256(5) 0.24(1) 0.30(1) 0.303(2) 0.422(8) 0.429(1) 0.236(3) 0.236(2)

EW (eV) ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... 63 (8) 61(10) 78(10) 82(10) 124(14) 151(13)

Ela or (keV) ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... 6.1 ± 6.1 ± 6.1 ± 6.09 ± 6.12 ± 5.78 ±

Fluxc 13.59 13.28 13.12 12.01 12.16 11.51 11.34 11.30 11.55 14.27 18.59 18.00 20.45 23.56 23.67 30.21 27.18 31.54 32.24 34.38 38.59 38.15 46.47 49.90 41.87 31.22



2 red

(dof) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73)

0.106(4) 0.102(3) 0.030(6) 0.100(6) 0.137(2) 0.077(8) 0.050(3) 0.025(3) -0.082(8) 0.095(4) -0.178(8) 0.033(3) 0.32(2) 0.07(4) -0.11(8) 0.13(6) 0.49(9) 0.030(3) -0.40(1) -0.26(9) -1.66(8) -0.58(7) -1.82(9) -1.40(8) -0.19(6) -0.22(9)

2.0 2. 0 2. 0 2. 0 2. 0 2. 0 2. 0 2. 0 2. 0 2. 0 2.0 2. 0 2. 0 0.09(6) 0.8(1) 0.18(1) -0.18(2) 2.0 2.0 2.0 2.0 2. 0 0.7(1) 2.0 2.0 0.47(5)

0. 7 1. 5 1. 0 0. 8 0. 6 1. 4

1.09 1.06 1.01 1.07 1.02 1.05 1.09 1.06 1.09 1.08 1.07 1.07 1.03 1.01 1.08 1.01 1.07 0.97 0.97 1.03 1.10 1.12 1.10 1.08 1.01 0.96

Notes. Parameter errors (put in parentheses) correspond to 1 confidence level. a The spectral model of the continuum is wabs (bmc + bmc highecut + laor ). b When the parameter log (A) > 1.0, this parameter is fixed to a value 2.0 (see comments in the text). c Spectral flux in the 3­150 energy range in units of â10-9 erg s-1 cm-2 . Table 2 Best-fit Parameters of Spectral Analysis of PCA and HEXTE Observation of GRS 1915+105 in 3­150 keV Energy Range During Middle 1997 Transition Observational ID 20402-01-45-02 20402-01-45-03 20402-01-46-00 20402-01-48-00 20402-01-50-00 20402-01-50-01 20402-01-51-00 20402-01-52-01 20402-01-52-02 20402-01-53-00 20402-01-53-02 20402-01-54-00 20402-01-55-00 MJD (day) 50696.21 50700.25 50703.41 50720.59 50735.54 50737.40 50743.29 50751.68 50751.75 50752.01 50756.41 50763.20 50769.22 1 = 1 - 1 1.84(1) 1.96(4) 2.00(7) 2.00(7) 1.74(7) 1.720(2) 1.704(2) 2.00(4) 2.00(8) 2.00(2) 2.00(2) 1.99(8) 0.8(3) log (A1 )
a

Nbmc1 , 2 L39 /d10 0.33(1) 0.23(4) 0.25(2) 0.240(6) 0.222(8) 0.236(1) 0.226(1) 0.240(6) 0.320(2) 0.32(1) 0.39(1) 0.412(6) 0.01(6)

2 = 2 - 1 2.98(2) 3.19(6) 3.1(2) 3.09(3) 3.2(2) 1.9(4) 2.0(3) 3.2(1) 3.2(1) 3.00(7) 3.20(1) 3.20(3) 3.15(3)

log (A2 )

a

Nbmc2 , 2 L39 /d10 0.379(2) 0.245(8) 0.13(2) 0.24(1) 0.323(9) 0.02(2) 0.04(5) 0.271(6) 0.482(5) 0.43(1) 0.32(2) 0.369(7) 0.426(4)

EW (eV) 133(5) 192(8) 87(10) 164(13) 125(19) 102(16) 102(13) 348(19) 342(20) 315(10) 53(12) 58(16) 101(13)

Ela or (keV) 6.31 6.10 6.25 6.46 6.48 6.47 6.47 6.31 6.44 6.10 6.18 6.19 6.25 ± ± ± ± ± ± ± ± ± ± ± ± ± 0.08 0.08 0.08 0.05 0.01 0.04 0.03 0.01 0.01 0.01 0.07 0.08 0.09

Fluxb 39.45 33.49 26.96 39.44 20.69 19.19 20.56 34.60 49.68 48.71 43.76 45.09 35.88



2 red

(dof) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73)

0.35(9) -0.1(2) 0.02(1) 0.77(4) 0.16(5) 2.0 1.9(6) 0.15(9) 0.06(1) -0.46(1) -0.45(3) 2.0 0.08(6)

0.56(1) 0.49(2) 0.8(4) 0.15(2) 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 0.322(1)

1.05 1.01 1.09 1.08 1.02 1.02 1.09 1.02 1.04 0.99 1.01 1.02 1.02

Notes. Parameter errors (put in parentheses) correspond to 1 confidence level. a When the parameter log (A) > 1.0, this parameter is fixed to a value 2.0 (see comments in the text). b Spectral flux in the 3­150 energy range in units of â10-9 erg s-1 cm-2 .

MJD = 50462­51081 and MJD = 53382­53852, as different types (samples) of bright X-ray activity, with transitions between hard and soft states. Two selected data sets have different patterns of radio/X-ray behavior and of light-curve shapes. We also use public GRS 1915+105 data from the ASM on-board RXTE (Swank 1999). The ASM light curves (2­ 12 keV energy range) were retrieved from the public RXTE/ ASM archive at HEASARC (http://xte.mit.edu/ASM_lc.html).

The monitoring Ryle Radio Telescope (RT; 15 GHz) data during the 1997­2006 period were kindly provided by Guy Pooley. The technical details of the RT are described by Pooley & Fender (1997). 2.1. Spectral Analysis
2.1.1. BMC and Iron Line Components of the Model Spectrum

The broadband source spectra were modeled in XSPEC with an additive model consisting of two BMC: a BMC with


1466

TITARCHUK & SEIFINA

Vol. 706

Table 3 Best-fit Parameters of Spectral Analysis of PCA and HEXTE Observation of GRS 1915+105 in 3­150 keV Energy Range During Decay 1997 Transition Observational ID 30184-01-01-00 30703-01-14-00 30402-01-09-00c 30402-01-10-00 30402-01-11-00 30703-01-15-00 30703-01-16-00 30402-01-12-01 30182-01-02-00 30182-01-04-01 30402-01-17-00 30703-01-33-00 30703-01-35-00 MJD (day) 50908.00 50909.87 50912.88 50914.39 50923.26 50925.88 50931.67 50945.01 51003.21 51006.21 51067.62 51071.90 51081.81 1 = 1 - 1 1.81(2) 1.74(1) 1.50(2) 1.50(1) 2.00(3) 1.79(5) 1.79(4) 1.69(2) 1.801(2) 1.0(4) 1.520(3) 1.141(7) 0.88(3) log (A1 )
a

Nbmc1 , 2 L39 /d10 0.260(1) 0.236(1) 0.20(3) 0.1949(8) 0.24(1) 0.22(3) 0.21(3) 0.206(4) 0.223(2) 0.169(8) 0.200(1) 0.172(1) 0.152(2)

2 = 2 - 1 2.0(3) 2.0(5) 2.09(5) 2.09(2) 3.2(1) 2.1(1) 2.10(9) 2.10(8) 2.0(3) 2.49(6) 2.1(3) 2.0(4) 2.1(2)

log (A2 )

a

Nbmc2 , 2 L39 /d10 0.009(1) 0.007(2) 0.11(1) 0.111(9) 0.23(3) 0.10(1) 0.09(2) 0.094(2) 0.009(1) 0.125(2) 0.009(2) 0.002(3) 0.009(1)

EW (eV) 284(18) 317(20) 290(18) 411(30) 308(20) 316(10) 326(8) 290(10) 434(11) ... ... ... ...

Ela or (keV) 5.75 6.39 5.7 5.7 5.6 6.4 6.4 5.7 6.41 ± ± ± ± ± ± ± ± ± ... ... ... ... 1. 2 0.05 0. 3 0. 4 0. 2 0. 2 0. 3 0. 5 0.01

Fluxb 24.90 22.85 23.89 22.92 31.28 23.79 23.52 22.17 26.11 24.23 21.15 19.52 18.18



2 red

(dof) (73) (73) (71) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73)

0.32(2) 0.55(2) -0.32(9) -0.32(8) 1.0(2) 0.37(6) 0.37(6) 0.42(3) 2.0 -1.953(8) 0.6(1) 0.19(1) 0.18(4)

2. 0 2. 0 2.0 2.0 2. 0 2. 0 2.0 -0.5(1) 2. 0 2.0 2. 0 2. 0 2. 0

1.09 1.30 1.04 1.02 1.23 1.01 1.07 1.06 0.99 1.03 0.98 1.07 1.02

Notes. Parameter errors (put in parentheses) correspond to 1 confidence a When the parameter log (A) > 1.0, this parameter is fixed to a value 2.0 b Spectral flux in the 3­150 energy range in units of â10-9 erg s-1 cm-2 c the data are fitted with the wabs*(bmc + bmc*highecut + laor + bbody)

level. (see comments in the text). . model; see values of the best-fit BB color temperature and EW in Table 7.

Table 4 Best-fit Parameters of Spectral Analysis of PCA and HEXTE Observation of GRS 1915+105 in 3­150 keV Energy Range During Rise 2005 Transition Observational ID 80701-01-48-00 90701-01-46-00 80701-01-37-00 90701-01-49-00 90105-05-03-00 90105-05-03-04 90105-05-03-05c 91701-01-04-00c 90105-07-01-00 90105-07-02-00 90105-07-03-00 91701-01-08-00 91102-01-01-00 91412-01-01-00 91701-01-10-00 91701-01-10-01 90105-08-02-00 90105-08-03-00 91701-01-11-00 91701-01-12-00 91701-01-13-00 91701-01-17-00 91701-01-19-00 91701-01-20-00 91701-01-24-00 91701-01-25-00 90105-04-01-00 90105-04-03-00 90105-04-03-01 MJD (day) 53382.39 53400.35 53416.34 53422.30 53442.98 53444.01 53444.08 53456.28 53472.92 53473.05 53473.97 53486.12 53488.72 53500.32 53501.06 53501.12 53503.87 53504.71 53508.06 53515.07 53520.05 53547.93 53562.92 53570.78 53600.89 53606.84 53640.38 53641.43 53641.51 1 = 1 - 1 1.78(5) 2.00(1) 2.00(4) 2.00(1) 1.9(1) 1.981(6) 1.982(4) 1.981(5) 1.980(6) 1.981(5) 1.98(1) 1.96(6) 1.60(5) 1.99(2) 2.0(1) 1.9(1) 1.90(3) 2.00(1) 1.70(3) 0.8(1) 0.8(1) 0.7(1) 0.7(1) 0.7(2) 0.7(1) 0.7(2) 0.7(1) 1.8(6) 1.8(7) log (A1 )
a

Nbmc1 , 2 L39 /d10 0.125(1) 0.292(4) 0.164(1) 0.152(4) 0.204 (3) 0.191(3) 0.187(1) 0.185(2) 0.144(1) 0.147(5) 0.147(2) 0.16 (3) 0.123(3) 0.261(1) 0.338(2) 0.210(4) 0.358(1) 0.242(2) 0.343(2) 0.001 (5) 0.05 (6) 0.06 (2) 0.09 (5) 0.05 (3) 0.123 (6) 0.09 (5) 0.04 (5) 0.215(4) 0.215(4)

2 = 2 - 1 2.1(4) 2.90(1) 2.92(3) 3.00(2) 2.92(4) 2.7(0.1) 2.79(6) 2.79(9) 2.85(6) 2.80(6) 2.80(8) 3.1(2) 3.0(1) 3.10(1) 3.20(2) 3.20(2) 3.20(1) 3.20(1) 3.10(1) 3.197(7) 3.200(3) 2.75(3) 3.20(1) 2.57(1) 2.75(9) 2.75(7) 2.16(2) 3.00(1) 3.00(1)

log (A2 )

a

Nbmc2 , 2 L39 /d10 0.009(2) 0.335(1) 0.157(4) 0.158(4) 0.1574(6) 0.111(2) 0.112(1) 0.110(1) 0.112(1) 0.133(3) 0.135(8) 0.38(3) 0.4705(5) 0.3400(9) 0.456(1) 0.433(4) 0.4901(6) 0.278(1) 0.451(1) 0.3480(1) 0.430(1) 0.193(1) 0.4633(7) 0.151(1) 0.124(1) 0.132(1) 0.1335(4) 0.282(3) 0.284(3)

EW (eV) ... 91(10) 594(10) 626(15) 435(15) 402(12) 378(10) 256(15) 55.6(10) 75.3(6) 210(8) 132(11) 74(12) 50(16) ... 68(9) 104(10) 180(12) 69(10) 90(10) 120(10) 95(10) 100(10) 125(10) 100(10) 100(10) 127(5) 114(12) 108(11)

Ela or (keV) 6.35 5.5 6.56 6.45 6.45 6.46 6.33 6.57 6.55 6.39 6.13 6.32 6.19 6.10 6.3 6.11 6.13 6.4 6.4 7.1 7.0 7.0 7.0 7.1 6.40 6.11 6.32 ... ± ± ± ± ± ± ± ± ± ± ± ± ± ... ± ± ± ± ± ± ± ± ± ± ± ± ± ± 0.01 0. 2 0.02 0.01 0.01 0.01 0.01 0.01 0.07 0.01 0.04 0.08 0.09 0.04 0. 1 0.01 0.04 0. 1 0. 1 0. 2 0. 1 0. 3 0. 3 0. 2 0.09 0.07 0.01

Fluxb 6.25 29.94 24.03 22.30 24.37 23.08 22.67 22.69 22.07 22.36 22.61 40.68 38.69 39.35 51.99 45.30 57.53 37.13 56.21 24.18 37.02 11.24 23.69 69.91 5.90 6.22 6.64 33.96 32.26



2 red

(dof) (73) (73) (73) (73) (73) (73) (71) (71) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73)

-0.96(1) -0.21(2) 2. 0 0.26(9) 2.0 2.0 2.0 2.0 2.0 1.0(2) 2. 0 2. 0 -0.47(7) -1.10(8) -1.11(9) -1.06 (8) -0.86 (6) -0.86(7) -1.7(9) 0.3(1) 0.9(4) 2.0 -0.9(5) -0.9(3) -0.95(5) -0.7(3) 0.4(1) 0.70(4) 0.48(4)

2.0 1.1(4) 2. 0 0.75(5) 2. 0 0.880(2) 0.88(5) 0.884(1) 0.251(1) 0.272(3) 2. 0 -0.02(3) -0.0030(5) 0.235(9) 0.34(1) 2. 0 0.802(8) 0.8(2) 2. 0 2.0 2.0 -1.106(1) 0.256(7) -0.537(8) -0.627(1) -0.626(1) -0.530(4) 1.1(2) 1.1(1)

1.05 0.96 1.03 1.01 1.06 1.09 0.96 0.97 0.99 0.99 0.96 1.01 1.08 1.02 1.11 1.39 1.09 1.02 1.04 1.09 1.10 0.99 1.02 1.18 1.24 1.25 1.35 1.08 1.01

Notes. Parameter errors (put in parentheses) correspond to 1 confidence level. a When the parameter log (A) > 1.0, this parameter is fixed to a value 2.0 (see comments in the text). b Spectral flux in the 3­150 energy range in units of â10-9 erg s-1 cm-2 . c These data are fitted with the wabs*(bmc*highecut + bmc + laor + bbody) model; see values of the best-fit parameters in Table 7.

high-energy cutoff (BMC1 component) and BMC2 component: wabs (bmc + bmc highecut). We also use a multiplicative wabs model taking into account an absorption by neutral material. The wabs model parameter is an equivalent hydrogen column NH .A

systematic error of 1% has been applied to the analyzed X-ray spectra. The BMC model describes the outgoing spectrum as a convolution of the input "seed" BB-like spectrum, whose


No. 2, 2009

DISCOVERY OF PHOTON INDEX SATURATION IN THE BB BINARY GRS 1915+105

1467

Table 5 Best-fit Parameters of Spectral Analysis of PCA and HEXTE Observation of GRS 1915+105 in 3­150 keV Energy Range During Middle 2005 Transition Observational ID 91701-01-31-00 91701-01-33-01 91701-01-33-00 91701-01-34-00 91701-01-34-01 91701-01-35-00 90105-06-03-01 90105-06-03-00 90105-06-03-02 91701-01-38-00 91701-01-38-01 91701-01-39-00 91701-01-39-01 92092-01-01-01 92092-02-01-01 92092-03-01-01 91701-01-39-01 91701-01-41-01 91701-01-42-00 91701-01-42-01c 91701-01-43-00c MJD (day) 53646.76 53659.98 53661.70 53669.69 53669.75 53674.59 53694.90 53695.03 53695.30 53696.56 53696.63 53703.57 53703.63 53704.94 53706.00 53707.68 53711.55 53718.56 53723.47 53723.54 53730.41 1 = 1 - 1 2.00(2) 1.60(1) 1.99(1) 1.42(1) 1.7(2) 1.50(5) 1.70(3) 1.89(2) 1.7(1) 1.80(8) 2.07(3) 1.70(2) 1.734(7) 1.74(3) 1.91(2) 1.7(3) 1.8(1) 1.77(2) 1.701(5) 1.91(1) 1.91(4) log (A1 )
a

Nbmc1 , 2 L39 /d10 0.17(1) 0.087(3) 0.112(2) 0.055(1) 0.084(2) 0.116(2) 0.345(1) 0.261(4) 0.282(6) 0.35(1) 0.685(4) 0.159(3) 0.170(1) 0.195(3) 0.116(3) 0.36(2) 0.23(2) 0.184(1) 0.111(2) 0.129(3) 0.15(2)

2 = 2 - 1 3.0(1) 2.90(6) 2.78(1) 2.51(1) 1.98(4) 3.00(2) 3.200(1) 3.20(1) 3.18(3) 3.12(2) 2.7(2) 2.90(1) 2.98(3) 2.9(1) 2.9(1) 3.00(6) 2.9(1) 2.1(5) 2.90(6) 2.70(6) 2.7(2)

log (A2 )

a

Nbmc2 , 2 L39 /d10 0.176(8) 0.143(2) 0.403(1) 0.189(1) 0.185(1) 0.3127(6) 0.552(1) 0.417(5) 0.389(5) 0.479(1) 0.290(3) 0.273(4) 0.163(1) 0.256(2) 0.229(1) 0.49(1) 0.26(2) 0.009(1) 0.110(1) 0.1644(5) 0.14(9)

EW (eV) ... 110(13) 92(8) 135(10) ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...

Ela or (keV) ... 6.35 ± 0.07 6.30 ± 0.01 6.32 ± 0.01 ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...

Fluxb 22.77 25.55 36.88 21.21 21.14 23.31 62.96 49.38 48.40 51.76 57.06 31.30 27.10 28.45 22.75 56.72 34.48 19.95 20.57 20.57 20.29



2 red

(dof) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (71) (71)

2.0 2.0 2.0 2.0 -0.57(7) -0.55(2) -1.4(2) -0.56(7) -0.70(8) -1.1(1) 0.68(8) -0.22(9) 0.35(2) 0.05(1) 2.0 -1.6(1) -0.33(5) 0.85(3) 0.8(1) 2.0 2.0

2. 0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 0.235(1) 2. 0 2.0 0.471(1) 2. 0 -0.363(1) 0.470(1) 2. 0 2. 0 2.0 2.0 2. 0

1.03 1.03 0.99 1.08 1.04 1.01 1.02 1.01 1.07 0.99 0.98 1.11 1.05 1.05 0.95 1.08 0.96 1.06 1.13 1.05 0.91

Notes. Parameter errors (put in parentheses) correspond to 1 confidence level. a When the parameter log (A) > 1.0, this parameter is fixed to a value 2.0 (see comments in the text), b Spectral flux in the 3­150 energy range in units of â10-9 erg s-1 cm-2 . c These data are fitted with the wabs*(bmc*highecut + bmc + laor + bbody) model; see values of the best-fit parameters in Table 7. Table 6 Best-fit Parameters of Spectral Analysis of PCA and HEXTE Observation of GRS 1915+105 in 3­150 keV Energy Range During Decay 2005­2006 Transition Observational ID 91701-01-46-00c 91701-01-49-00c 91701-01-50-00c 91701-01-22-00 92702-01-01-00 92702-01-01-01 92702-01-02-01 92702-01-03-00 92702-01-05-00 92082-01-05-00 92702-01-07-00 92702-01-08-00 92702-01-08-01 92702-01-08-02 90105-02-04-00 MJD (day) 53753.43 53771.51 53778.30 53794.31 53803.28 53803.35 53809.30 53815.19 53829.20 53834.84 53844.16 53851.10 53851.17 53851.23 53852.15 1 = 1 - 1 1.70(1) 1.69(1) 1.65(1) 1.68(2) 1.40(2) 1.50(3) 1.353(9) 1.49(2) 0.75(3) 1.03(1) 1.10(2) 0.94(3) 1.13(1) 1.04(4) 1.09(2) log (A1 )
a

Nbmc1 2 L39 /d10 0.284(7) 0.221(8) 0.143(2) 0.143(2) 0.12(1) 0.11(1) 0.085(2) 0.094(4) 0.125(1) 0.140(2) 0.163(1) 0.127(2) 0.132(1) 0.128(1) 0.125(1)

2 = 2 - 1 3.00(2) 2.41(3) 2.54(7) 2.00(7) 1.8(2) 1.8(1) 1.74(5) 1.21(6) 1.9(5) 2.0(6) 1.9(4) 2.0(3) 2.0(4) 1.9(3) 1.8(6)

log (A2 )

a

Nbmc2 , 2 L39 /d10 0.3807(8) 0.0978(2) 0.1055(5) 0.081(2) 0.083(4) 0.081(7) 0.0798(4) 0.063(2) 0.005(3) 0.009(2) 0.005(3) 0.01(2) 0.008(3) 0.01(3) 0.009(2)

EW (eV) 108(10) ... ... ... ... ... ... ... ... ... ... ... ... ... ...

Ela or (keV) 5.15 ± 0.06 ... ... ... ... ... ... ... ... ... ... ... ... ... ...

Fluxb 38.60 16.00 17.48 18.60 19.35 19.11 16.64 17.46 16.98 16.05 14.92 14.50 14.74 14.03 14.28



2 red

(dof) (71) (71) (71) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73) (73)

0.6(1) 2.0 2.0 2.0 2.0 2.0 2.0 2.0 0.14(1) 0.13(1) -0.15(2) 0.06(1) 0.09(1) 0.05(2) 0.15(1)

2. 0 2. 0 0.38(1) -0.79(1) 0.29 (1) 2.0 2. 0 -0.20(1) 2.0 2.0 2.0 2.0 2.0 2.0 2.0

0.93 1.06 0.95 0.96 1.02 1.04 1.22 1.02 1.02 0.98 1.06 1.07 0.99 0.95 0.98

Notes. Parameter errors (put in parentheses) correspond to 1 confidence level. a When the parameter log (A) > 1.0, this parameter is fixed to a value 2.0 (see comments in the text). b Spectral flux in the 3­150 energy range in units of â10-9 erg s-1 cm-2 . c These data are fitted with the wabs*(bmc*highecut + bmc + laor + bbody) model; see values of the best-fit parameters in Table 7.

normalization is Nbmc and color Comptonization Green's function. XSPEC model, the normalization (disk) luminosity to the square of Nbmc = L 1039 erg s

temperature is kT , with the Similar to the ordinary bbody Nbmc is a ratio of the source the distance: 10 kpc d
2

-1

.

(1)

The resulting spectrum is characterized by the parameter log (A) related to the Comptonized fraction f as f = A/(1 + A) and

a spectral index = - 1. There are several advantages of using the BMC model with respect to other common approaches used in studies of X-ray spectra of accreting compact objects, i.e., sum of BB/multi-color-disk and power-law/ thermal Comptonization. First, the BMC, by the nature of the model, is applicable to the general case where there is an energy gain through not only thermal Comptonization but also via dynamic (bulk) motion Comptonization (see Shaposhnikov & Titarchuk 2006, for details). Second, with respect to the phenomenological powerlaw model, the BMC spectral shape has an


1468

TITARCHUK & SEIFINA

Vol. 706

Figure 1. Schematic view of the proposed geometry for thermal and bulk Comptonization regions in the source hosting a BH with power-law-like emission at high energies. The thermal plus bulk Comptonization spectrum (thermal plus bulk BMC1) arises in the innermost part of the TL, where the disk BB-like seed photons are (thermally and dynamically) Comptonized by the infalling material. Whereas the thermal Comptonization spectrum (thermal BMC2) originates in the outer part of the TL region. (A color version of this figure is available in the online journal.)

appropriate low-energy curvature, which is essential for a correct representation of the lower energy part of the spectrum. Our experience with powerlaw components shows that the model fit with this component is often inconsistent with the NH column values and produces an unphysical component "conspiracy" with the highecut part. Specifically, when a multiplicative component highecut is combined with BMC, the cutoff energies Ecut are in the expected range of 20­30 keV, while in a combination with powerlaw, Ecut often goes below 10 keV, resulting in unreasonably low values for the photon index. As a result, the implementation of the phenomenological powerlaw model makes much harder, or even impossible to correctly identify the spectral state of the source, which is an imminent task for our study. Third, and even a more important property of the BMC model, it calculates consistently the normalization of the original "seed" component, which is expected to be a correct mass accretion rate indicator. Note that the Comptonized fraction is properly evaluated by the BMC model. We consider a scenario related to our spectral model (see Figure 1) where the Compton cloud (CC) along with converging flow (CF) is located in the innermost part of the source and the Keplerian disk is extended from the CC to the optical companion (see, e.g., TF04). An iron K -line (la or ) component (Laor 1991) was included in our model spectrum. To summarize the spectral model parameters are the equivalent hydrogen absorption column density NH ; spectral indices 1 , 2 (photon index = + 1); color temperatures of the BB-like photon spectra kT1 , kT2 ; log (A1 ), log (A2 ) related to the Comptonized fractions f1 , f2 [f = A/(1 + A)]; normalizations of the BB-like components Nbmc1 , Nbmc2 for the BMC1 and BMC2 components of the resulting spectrum, respectively. We find that color temperatures kT1 and kT2 are about 1 keV for all available data and thus we fix values of kT1 and kT2 at 1 keV. An equivalent hydrogen absorption column density was fixed at the level of NH = 5 â 1022 cm-2 (Trudolyubov 2001). When the parameter log (A) 1, we fix log (A) = 2 (see Tables 1­6), because the Comptonized fraction f = A/(1 + A) 1 and variation of A do not improve the fit quality any more. During LHS, the BMC2 component is often very low or barely detectable (see Table 1 (MJD = 50462-50561), Table 3 (MJD = 51067-51081), Table 6 (MJD = 53829-53852), and

Figure 2 (panel "1S")). This observational fact is in agreement with scenario of BH spectral state transition (TF04). During LHS, the spectrum is characterized by a strong hard power-law component. In other words, the energy spectrum is dominated by a Comptonized component seen as a power-law hard emission in the energy range from 10 to 70 keV while the disk emission remains weak (LHS, IS), because only a small fraction of the disk emission component (1 - f ) is directly seen (Figure 2, panel "1S"). Although during IS (Figure 2, panel "2S"), the contributions of two BMC components to the overall spectrum are of the same order that sometimes we can barely identify one of these components (see these cases in Table 2 (MJD = 50737, 50743), Table 3 (MJD = 50908, 50909, 51003), Table 5 (MJD = 53718) when Nbmc2 Nbmc1 ). On the other hand, the model with two BMC components is really needed in the most cases of the IS and HSS. In Figure 3, we demonstrate that the fit qualities are unacceptable when the only one BMC component is included in the spectrum. Specifically, for IS­HSS observation 2 91701-01-11-00 on 2005 May 18 the values of red = 12.3 for 75 degrees of freedom (dof; bottom left panel) for the wa bs bmc hig hecut model. However, 2 is significantly improved when the second BMC component is included in the 2 model. For wa bs (bmc + bmc hig hecut ) red = 3.28 for 76 dof (see the related count spectrum along with the model and the related residual in the central bottom panel of Figure 3). Ultimately, we achieve a remarkable agreement with the data using the wa bs (bmc + bmc hig hecut + la or ) model for 2 which red = 1.04 for 73 dof. In Figure 3 (top and right bottom panels), we show the data along with the best-fit spectra and their components for our two-component BMC model (see Table 4 for the best-fit parameter values). In HSS and VSS (MJD = 53515­53640, Table 4), the soft luminosity is high and the spectrum is dominated by a thermal component (2 > 3.7). Note that the hard power-law component is barely seen in VSS (see Figure 2, panels "4S"). We find a broad emission line between 6 and 7 keV in the IS and HSS spectra. Up to the date, the iron K emission line in IS and HSS of GRS 1915+105 was detected with ASCA, Chandra, XMMNewton, and BeppoSAX by Kotani et al. (2000), Martocchia et al. (2002), and Miller & Homan (2005). However, the determination of the iron line shape with the RXTE is a problem because of the low-energy resolution. As a first trial, we added a Gaussian component to fit the spectrum varying the width and normalization of the line. Fits using the Gaussian always produce residuals around 7 keV while fits with the XSPEC Laor model do not have such a problem. Thus through the paper we incorporate the Laor-line profile to fit the line component. The line feature has a statistical significance of (3­10) depending on the spectral states. This line feature is variable with a time average intensity of 2â10-11 and exhibits equivalent width (EW) in the range of 50­600 eV across the data. We found that adding the Laor-line component significantly improves the fit quality of IS and HSS spectra. Fitting an IS spectrum (e.g., for 90105-05-03-05 observation) 2 without the line components leads to red = 1.24. When the line component (Laor) is included, the fit quality becomes much 2 better, red = 1.01. The fit of the HSS spectrum (90105008-02-00) without the iron-line component is unacceptable: 2 2 red = 2.01 and red = 1.24 when the line component is included. The best-fit parameters of the source spectrum and 2 values of red including degree of freedom are presented in Tables 1­6.


No. 2, 2009

DISCOVERY OF PHOTON INDEX SATURATION IN THE BB BINARY GRS 1915+105

1469

Figure 2. Evolution of spectum shape of GRS 1915+105 during LHS, IS, HSS, and VSS spectral states. Data are taken from RXTE observations 20402-01-11-00 (S1, 1 = 1.8, LHS), 91701-01-33-00 (S2, 1 = 2.9, 2 = 3.7, IS), 91701-01-11-00 (S3, 1 = 2.7, 2 = 4.1, HSS), and 91701-01-19-00 (S4, 2 = 4.2, VSS). Here data are denoted by red points; the spectral model presented with components is shown by blue, black, and dashed purple lines for BMC1, BMC2, and laor components, respectively. (A color version of this figure is available in the online journal.)

2.1.2. Observational Evidence of "BB-like" Component Peaked at 20 keV in Eight IS Spectra

The adopted spectral model shows a very good performance for 99 cases among 107 spectra used in our analysis. Namely, 2 the value of reduced 2 statistic red = 2 /Ndof , where Ndof is a number of degree of freedom for a fit, is less or around 1.0 for most observations. However, for eight observations of IS the fit of the data with the model wabs*(bmc*highecut+bmc+ laor) is 2 not so good; red reaches 1.5 and even higher. We found that in the residual of data versus model there is a characteristic bump around 20 keV (see the left bottom panel of Figure 4) which can be fitted by BB-like shape of color temperature about 4.5 keV (see the right bottom and top panels of Figure 4 and Table 7 for values of the best-fit parameters). This "high-temperature BB" component is strong in each of the eight observations, and its EW varies from 400 to 700 eV. We discuss a possible origin of this "BB" component in Section 4. 2.2. Timing Analysis The RXTE light curves were analyzed using the powspec task. For the timing analysis in the 2­30 keV energy range, we use

the RXTE/PCA data in the binned and event modes containing X-ray events below 13/15 keV and above 13/15 keV for 1997/ 2005 data sets, respectively. Specifically, depending on RXTE epoch, the channel ranges (0-35) for binned and (36­255) for event modes relate to energy bands 1.94­12.99 keV (binned) and 13.36­103.19 keV (event) for epoch 3 (1997­1998 data set), and relate to 2.06­14.76 keV (binned) and 15.18­117.86 keV (event) for epoch 5 (2005­2006 data set). The time resolutions for event and binned modes are 1.52 â 10-5 s and 8 â 10-3 s, respectively. The observational exposition time periods vary from 1.5 to 10 ks. Thus for all of these observations we can obtain power spectra in the wide frequency range (roughly from 0.001 Hz up to 100/10000 Hz for binned and event modes, respectively). These frequency ranges allow us to produce power spectra for all studied cases in the 0.1­ 100 Hz frequency range. We subtracted the contribution due to Poissonian statistics and Very Large Event Window for all PDSs. The data analysis of the PDSs was performed using a simplified version of the diffusion model (see Titarchuk et al. 2007, hereafter TSA07, and Titarchuk & Shaposhnikov 2008) in which the PDS continuum shape at frequencies below


1470

TITARCHUK & SEIFINA

Vol. 706

Figure 3. Best-fit spectrum of GRS 1915+105 during HSS events of 2005 rise transition (observation 91701-01-11-00) in E F (E ) units (top) along with the count spectrum and (bottom panels). Bottom panels (left): the best-fit spectrum and for the model fit, wa bs bmchighecut, which include a single BMC component 2 without an iron-line component (red = 12.3 for 75 dof), center: the model fit, wa bs (bmc + bmchighecut), which includes two BMC components without an 2 iron-line component (red = 3.28 for 76 dof), and right: the same as the latter one but adding of the Laor line component, wa bs (bmc + bmc hi g hecut + laor ), 2 (red = 1.04 for 73 dof, see Table 4). Blue, black, purple lines stand for BMC1, BMC2, and Laor components, respectively. (A color version of this figure is available in the online journal.)

the driving frequency can be approximated by an empirical model PX (1.0 + (x/x )2 )-in (KING model in QPD/PLT). Following TSA07, the break frequency found in the PDS is related to a diffusion time of perturbation propagation while the QPO low frequency is an eigenfrequency of the volume (magnetoacoustic) oscillation of the medium (in our case it is a CC). Note, TSA07 demonstrated that the diffusion model as a linear superposition of Lorentzians related to the eigenvalues of the diffusion problem can also be presented by the continuous shape which is flat below the break frequency and the power law at frequencies above the break. Given these asymptotic forms of the PDS at low- and high-frequency limits they named their diffusion model as a WRN model. For the quasi-uniform perturbation source distribution the slope of the PDS power-law part depends on the law of viscosity in the corona or in the disk

(see details in TSA07). The parameters of this WRN diffusion model are the break (diffusion) frequency and the index of the power-law distribution of the viscosity over the radius. To fit the QPO features, we use the Lorentzian shape. We quote the Lorentzian centroid as a QPO frequency. 3. OBSERVATIONAL RESULTS 3.1. Evolution of Spectral Properties During State Transitions Observations of Galactic BH X-ray binaries reveal diverse spectral and dynamic phenomenologies. The evolution of a BH binary is usually described in terms of spectral states. There are several flavors of BH state definitions in the literature, which slightly differ in BH state definitions and terminology (see, for example, Remillard & McClintock 2006; Belloni et al.


No. 2, 2009

DISCOVERY OF PHOTON INDEX SATURATION IN THE BB BINARY GRS 1915+105

1471

Figure 4. Best-fit spectrum during IS-LHS of the 2005­2006 decay transition in EF (E ) units (top) and in normalized counts units (bottom panel) with for the 2 91701-49-00 observation. Left: fitting without modeling the high-temperature "bbody" component (red = 1.43 for 73 dof) and right: the best-fit spectrum and , 2 when the bump in residuals at 20 keV is modeled by a "high-temperature bbody-like" component with red = 1.06 for 70 dof. On the top panel, the data are denoted by red points; the spectral model presented with components i shown by blue, black, purple, and green lines for BMC1, BMC2, laor, and "high-temperature" bbody components, respectively. (A color version of this figure is available in the online journal.)

2000, 2005; Klein-Wolt & van der Klis 2008). To distinguish different states, the properties observed in the energy spectrum and Fourier PDS are utilized. As we have already emphasized in the introduction section, we use, in our study, the general state classification for four major BH states: LHS, IS, HSS, and VSS. The general picture of LHS­IS­HSS transition is illustrated in Figure 5 where we bring together spectra of LHS, IS, HSS, and VSS to demonstrate the source spectral evolution from lowhard to soft states. We should emphasize different shapes of the spectra for the different spectral states. In the LHS spectrum, the Comptonization component is dominant and the BB component is barely seen in the 3­150 keV energy range. The IS and HSS spectra are characterized by a strong soft BB component and a power law extended up to 150 keV. In VSS, the soft

BB component is dominant and the power-law component is relatively weak with respect of this in IS and HSS. In the RXTE data of GRS 1915+105 observations, there are long periods when the photon index 1 and normalization Nbmc1 of the hard BMC monotonically increase (or decrease) with time. We call these periods as long transition periods. The days, when the source X-ray flux starts to increase while it is still in the LHS, can be considered as a beginning of the rise transition. In these times, the energy spectrum is characterized by low index values 1 1.7 and the thermal component is at a low level or undetectable at all. In Figure 6, from top to bottom we show an evolution of flux density S15 GHz at 15 GHz (Ryle Telescope), RXTE/ASM count rate, BMC normalization, and photon index during the 1997 rise transition of GRS 1915+105 (MJD 50500­50700). Red/black points (for two low


1472

TITARCHUK & SEIFINA
Table 7 Parameters of the Model: wabs*(bmc*highecut+bmc+laor+bbody) for Selected IS Data

Vol. 706

Model bmc1

Parameter 30402-01-09-00 90105-05-03-05 91701-01-04-00 91701-01-42-01 91701-01-43-00 91701-01-46-00 91701-01-49-00 1 = 1 - 1 kT1 (keV) Nbmc1 1.50 ± 0.02 1.982 ± 0.004 1.981 ± 0.005 0. 9 ± 0 . 1 1. 1 ± 0 . 1 1.09 ± 0.1 0.20 ± 0.03 0.187 ± 0.001 0.185 ± 0.002 1.91 ± 0.01 0.99 ± 0.05 0.129 ± 0.003 1.91 ± 0.04 1.00 ± 0.01 0.154 ± 0.015 2. 7 ± 0 . 2 0.99 ± 0.09 0.14 ± 0.09 5.39 ± 0.01 0.06 ± 0.01 180 ± 10 4.68 ± 0.09 0.009 ± 0.001 390 ± 30 1.70 ± 0.01 1.69 ± 0.01 1.02 ± 0.02 0.95 ± 0.06 0.284 ± 0.007 0.221 ± 0.008

91701-01-50-00 1.65 ± 0.01 0.89 ± 0.09 0.143 ± 0.002

bmc2 2 = 2 - 1 kT2 (keV) Nbmc2 laor (keV) 5. 7 ± 0 . 3 6.46 ± 0.01 6.33 ± 0.02 Nla or 0.073 ± 0.02 0.062 ± 0.001 0.028 ± 0.002 EWla or (eV) 290 ± 18 378 ± 10 256 ± 15 (eV) E
la or

2.09 ± 0.05 2.79 ± 0.06 2.79 ± 0.09 2.70 ± 0.06 1. 0 ± 0 . 1 1. 0 ± 0 . 1 1. 1 ± 0. 1 1. 0 ± 0. 2 0.11 ± 0.01 0.112 ± 0.001 0.110 ± 0.001 0.1644 ± 0.0005 5.39 ± 0.01 0.06 ± 0.01 150 ± 15 4. 4 ± 0. 2 0.041 ± 0.005 710 ± 50

3.00 ± 0.02 2.41 ± 0.03 2.54 ± 0.07 1. 0 ± 0 . 2 0.99 ± 0.08 1. 1 ± 0 . 3 0.3807 ± 0.008 0.0978 ± 0.0002 0.1055 ± 0.0005 5.15 ± 0.06 6.40 ± 0.09 0.022 ± 0.002 0.016 ± 0.004 108 ± 10 100 ± 10 4.77 ± 0.04 4.47 ± 0.09 6.40 ± 0.02 0.016 ± 0.009 100 ± 10 4. 5 ± 0 . 1 0.048 ± 0.004 630 ± 40

"bbody" T"bbody " (keV) N"bbody " EW"bbody (eV) Fluxa 2 (dof) S15
GHz b

4. 1 ± 0. 2

4. 5 ± 0 . 3

4. 3 ± 0. 3

0.04 ± 0.01 0.011 ± 0.006 0.036 ± 0.006 550 ± 30 430 ± 20 490 ± 40

0.024 ± 0.006 0.038 ± 0.003 470 ± 30 540 ± 30

23.89 1.04 (71) 88 ± 2

22.67 0.96 (71) 156 ± 2

22.69 0.97 (71) 140 ± 2

20.57 1.01 (71) 41 ± 2

20.29 0.92 (71) 51 ± 2

38.60 0.93 (71) 65 ± 2

16.00 1.06 (71) 40 ± 2

17.48 0.95 (71) 20 ± 2

Notes. See the full list of observations in Tables 1­6. Parameter errors (put in parentheses) correspond to 1 confidence level. a Spectral flux in the 3­150 energy range in units of â10-9 erg s-1 cm-2 . b Spectral flux density in radio band centered at 15 GHz in units of mJy.

panels) correspond to hard/soft components with 1 and 2 , respectively. In the bottom, we plot the photon index versus the BMC normalization (left) and Comptonized fraction (right) for this transition. Here red triangles/black circles correspond to hard/soft components with 1 and 2 , correspondingly. One can see that in the beginning of this transition the resulting spectrum consists of one Comptonization component whose photon index 1 steadily increases from 1.7 toward the softer states and finally saturates at the value of 3. The Comptonization fraction f = A/(1 + A) of the hard component, related to the index 1 , shows a sign of decreasing toward the softer state. When the RXTE/ASM count rate exceeds 50 counts s-1 the soft Comptonized component appears with a weight which is comparable to that of the hard component. The photon index of the soft component 2 saturates to the level of 4.2 when the BMC normalization (disk flux) increases. The Comptonization fraction f of the soft component is about 0.5 and higher. As seen from Figure 6, the start of this rise transition coincides with the active phase of X-rays, and of radio emissions which exceed 10 ASM counts s-1 and 50 mJy levels, respectively, in the 1997 rise transition. Around MJD 50580 day, the source reaches the HSS (when 1 3). Then a long HSS period from MJD 50580 to 50700, when 1 stays almost the same, is followed by the state transition to IS during which 1 decreases to 2.5 (see Figures 7 and 8). We see a similar behavior of X-ray, radio fluxes, and photon indices during the 2005 bright X-ray episode. The only difference of that with the 1997 episode was that the 2005 rise started at the IS and went very quickly to the HSS (see Figures 9 and 10). After MJD 53800 the source came back to IS­LHS when 1 2 (see Figure 10 and Table 6). Note typical X-ray and radio fluxes during IS are about 40­60 ASM counts s-1 and 50 mJy, respectively.

Figure 5. Four representative EFE spectral diagrams during LHS, IS, HSS, and VSS spectral states of GRS 1915+105. Data are taken from RXTE observations 20402-01-11-00 (blue, LHS), 91701-01-33-00 (red, IS), 91701-01-11-00 (green, HSS), and 91701-01-19-00 (purple, VSS). (A color version of this figure is available in the online journal.)

3.2. Observational (Correlated and Non-correlated) Characteristics of X-ray and Radio Emissions In fact, we do not find any correlation of X-ray and radio fluxes when the source changes its spectral states. Also we do not find a correlation of radio activity with the X-ray photon index (see Figure 11). However, we find a strong correlation


No. 2, 2009

DISCOVERY OF PHOTON INDEX SATURATION IN THE BB BINARY GRS 1915+105

1473

Figure 6. From top to bottom: evolutions of the flux density S15 GHz at 15 GHz (Ryle Telescope), RXTE/ASM count rate, BMC normalization, and photon index in the beginning of the 1997 rise transition of GRS 1915+105 (MJD 50460-50700). Red/black points (for two last panels) correspond to hard/soft components with 1 and 2 , respectively. Bottom: spectral index plotted vs. BMC normalization (left) and Comptonized fraction (right) for this transition. Here the red triangles/black circles correspond to hard/soft components with 1 and 2 , correspondingly. Note that in most cases the normalization of the hard BMC component (BMC1) is higher than that of the soft component (BMC2; see red points vs. black points in the BMC normalization time panel). (A color version of this figure is available in the online journal.)

of the iron-line EW with radio flux density S15 GHz at 15 GHz (see Figure 12). In Figure 12, we also include points which have recently been reported by Neilsen & Lee (2009) who have analyzed archival High Energy Transmission Grating Spectrometer (HETGS) observations of GRS 1915+105 from the Chandra X-ray Observatory. The prominent HSS events were observed during 2005­ 2006 observations around MJD 53490 and MJD 53690 (see Figures 9 and 10). The 2005­2006 observations confirm the index evolution versus BMC normalization (disk flux) found in 1997­1998. The indices of the hard and soft components indeed increase and then saturate at values of 3 and 4.2, respectively (see Figure 10). The index 1 started a saturation at lower values of BMC1 normalization (presumably proportional to disk mass accretion rates) than that were in 1997. It is also worth noting a so-called "pivoting" effect, i.e., when inequality Nbmc1 > Nbmc2 switches to Nbmc2 > Nbmc1 or vice versa, which is seen in the 1997­1998 and 2005­2006 observations. One can see this composite pivoting picture combining

Figures 6­10. In fact, these pivoting points correspond to the spectral transitions between adjacent states LHS­IS to HSS and vice versa. In Figure 13, we collect all data points for the indexnormalization correlation for rise and decay stages. We do not find much differences in the correlation patterns related to the rise and decay transitions (compare the left and right panels) in contrast to that ST09 found in other BHs. We also find that the photon index of the X-ray spectrum is tightly correlated with the QPOs frequency (see Figure 14) which can be considered as a strong argument that QPOs and X-ray Comptonization spectrum emerge from the same geometrical configuration (CC). However, the flux density S15 GHz and QPO frequency are not correlated with each other when the source changes its spectral states. In Figure 15, we show an evolution of the flux density S15 GHz at 15 GHz (Ryle Telescope), RXTE/ASM count rate and QPO during 1997 (left column) and 2005 (right column) rise transitions. The left column panel demonstrates the presence of QPO during a low


1474

TITARCHUK & SEIFINA

Vol. 706

Figure 7. From top to bottom: the same types of evolutions which are presented in Figure 6 but now that for the middle of the 1997 rise transition of GRS 1915+105 (MJD 50650-50760). In most cases, the normalization of the soft BMC component (BM C 2) is higher than that of the hard component (BMC1; see black points vs. red points in the BMC normalization­time panel and compare with that in Figure 6). The pivoting effect occurs (i.e., Nbmc1 > Nbmc2 switches to Nbmc2 > Nbmc1 around MJD 50690. (A color version of this figure is available in the online journal.)

radio activity (<30 mJy). The right column panel shows an example of the presence of QPO when the radio flux is high (100­200 mJy). Given that the QPOs of X-ray flux are present independently of the radio flux level, we can conclude that the radio appearances and QPO phenomenon are not closely related and probably the radio and X-ray (oscillating) emission areas have different origins in the source. 3.3. Evolutions of Energy and Power Spectra During a Minor X-ray/Radio Flares In Figure 16, we show the details of a typical evolution of X-ray timing and spectral characteristics for minor X-ray/radio flares. In the top panels of Figure 16, we show the flux density S15 GHz at 15 GHz (Ryle Telescope) and the RXTE/ASM count rate during the 2005 rise transition stage (see also Figure 9). Red points A, B, and C on the panel of the RXTE/ASM count rate versus time correspond to moments at MJD = 53416, 53422, and 53442 (before, during, and after the minor X-ray/radio flare), respectively. Points A and C were chosen as the nearest possible points to point B (taking into account the time table of

archive data). Point B corresponds to the maximum of a radio flux of 300 mJy and a EW of 600 eV. Note that QPO centroid frequency before the flare (at point A) is at 1.8 Hz and shifts to 0.9 Hz (point C) after the flare. PDSs (left bottom column) are plotted versus the energy spectrum (right bottom column) for three points A (top), B (middle), and C (bottom) of the X-ray light curve. There are QPOs at A and C points (A1, C1 panels) but there is none at B point (B1 panel), at the X-ray flare peak. For the photon spectra (right bottom column), red points stand for observational data, while the model is shown by components with a blue line for BMC1, a black line for BMC2, and a dashed purple line for the laor line component, respectively. Note that the spectral characteristics undergo noticeable changes during X-ray/radio flare. Specifically, at the flare peak (point B) the total flux increases at least by a factor of 1.5 with respect to that before the flare, although the photon index of the BMC1 component 1 changes from 2.9 (A and C points) to 3.0 (B point), respectively. We also studied the energy dependence of the PDS shape and integrated power variability as a function of the photon energy.


No. 2, 2009

DISCOVERY OF PHOTON INDEX SATURATION IN THE BB BINARY GRS 1915+105

1475

Figure 8. From top to bottom: the same types of evolutions which are presented in Figures 6 and 7 but now that for the 1997­1998 decay transition of GRS 1915+105. In most cases, the normalization of the hard BMC component (BM C 1) is higher than that of the soft component (BM C 2; see red points vs. black points in the BMC normalization­time panel). The pivoting effect occurs around MJD 50910. (A color version of this figure is available in the online journal.)

In Figure 16 (left bottom column pannel), we show two power spectra for two energy bands 2­15 keV (red) and 15­30 keV (blue). One can see that PDSs weakly depend on the energy band. In particular, a value of the LFQPO QPO is the same for the low-energy and high-energy PDSs. 4. INTERPRETATION AND DISCUSSION OF OBSERVATIONAL RESULTS Before to proceed with the interpretation of the observations, let us briefly summarize them as follows. (1) The spectral data of GRS 1915+105 are well fit by two (soft and hard) BMC components for most of analyzed IS and HSS spectra (see Figure 3) while LHS spectra essentially require only one BMC component: the soft BMC component is very weak (see Tables 1, 3, 4, 6, and panel S1 in Figure 2). (2) In addition to two BMC components, eight IS­HSS spectra require an extra component which can be fitted by the "high-temperature BB-like" profile (see Figure 4 and Table 7). (3) The Green's function indices of each of these components rise and saturate with an increase of the BMC normalization (disk flux). The photon index saturation levels of the soft and hard components are about 4.2 and 3,

respectively (see Figure 13). (4) We also find a tight positive correlation of QPO frequencies with the index (see Figure 14) and consequently that with the disk flux. (5) The iron-line EW correlates with the radio flux (see Figure 12). (6) QPO appearances and their frequency values are not correlated with radio flux when the source undergoes the spectral changes from IS to HSS (see Figure 15). (7) We also do not find any correlation between X-ray and radio fluxes, and X-ray power-law index (see Figure 11). (8) Although we find changes of power and energy spectra during a minor X-ray/radio flares when QPO features disappear in the PDS and the energy spectrum becomes softer than that was before and after the flare (see Figure 16). 4.1. Index­QPO and Index­m Correlations: Index Saturation We confirm the index­QPO correlation in GRS 1915+105 previously found by V03 and ST07. This correlation was indeed predicted by Titarchuk et al. (1998), hereafter TLM98, who argued that the transition layer (CC) formed between the Keplerian disk and the central object (NS or BH) contracts and becomes cooler when the disk mass accretion rate m increases. The observational effect of the CC contraction was


1476

TITARCHUK & SEIFINA

Vol. 706

Figure 9. From top to bottom: the same types of evolutions which are presented in Figures 6­8 but now that for the 2005 rise transition of GRS 1915+105. The pivoting effect occurs around MJD 53475. (A color version of this figure is available in the online journal.)

later demonstrated by ST06, TSA07, TS08, and Montanari et al. (2009) in Cyg X-1 and XTE J1650-500, respectively. As a result of the transition layer (TL) contraction, the QPO low-frequency L , which is presumably the TL's normal mode oscillation frequency, increases with m given that L is inversely proportional to the TL size. On the other hand, the index monotonically increases when the TL (CC) cools down. TF04 provided the details of the index­QPO correlation model where they pointed out that this correlation is a natural consequence of the spectral state transition. In this paper, we have firmly established the index correlation with L along with the index saturation versus the BMC normalization Nbmc (Equation (1)) for the soft and hard Comptonized components of the X-ray spectra of GRS 1915+105 (see Figure 13). Below we show that Nbmc is actually proportional to mass accretion rate in the disk. Namely, the disk flux L (as a source of soft photons for Comptonization, see, e.g., Figure 1 for the geometry of soft photon illumination of the Comptonized region) can be represented as GMbh M = (r )md LEd , R

where R = r RS is an effective radius where the main energy release takes place in the disk, RS = 2GM /c2 is the Schwarzschild radius, = 1/(2r ), md = Md /Mcrit is the dimensionless mass accretion rate in units of the critical mass accretion rate Mcrit = LEd /c2 , and LEd is the Eddington luminosity. On the other hand, L i.e., L
Ed Ed

=

4GM mp c , T

(3)

M and thus using Equations (2)­(3) we obtain that L (r )md m. (4)

L=

(2)

In HSS when the inner disk radius R reaches its lowest value R 3RS , the efficiency of the gravitational energy release (r ) reaches its highest value, and thus the disk flux increases only when the disk mass accretion rate increases (see Equation (4)). Given that BMC normalization Nbmc is proportional to md in HSS, the observational effect of the index saturation with Nbmc is translated to the index saturation with md .


No. 2, 2009

DISCOVERY OF PHOTON INDEX SATURATION IN THE BB BINARY GRS 1915+105

1477

Figure 10. From top to bottom: the same types of evolutions which are presented in Figures 6­9 but now that for the 2005­2006 middle and decay transitions. The pivoting effect occurs between MJD 53750 and MJD 53775. (A color version of this figure is available in the online journal.)

First we interpret the index saturation related to the hard Comptonization component (BMC1). We suggest that this BMC1 component of the emergent spectrum is presumably originated in the CF onto a compact object, in our case to the BH (see Figure 1). In fact, in HSS the plasma temperature of the accretion flow is comparable with the color temperature of the disk photons (see TF04). Thus, in order to explain the high-energy tail observed in HSS of BH sources, one should assume either an unknown source of high-energy nonthermal electrons (see, e.g., Coppi 1999) or consider effects of energy transfer from the CF electrons to the photons emitted from the innermost part of the accretion flow. Optical depth of the CF is proportional to md if one assumes that disk accretion flow continuously goes to the CF and there are no other components in the accretion flow (see, e.g., a model of two-component accretion flow by Chakrabarti & Titarchuk 1995, hereafter CT95). This effect of the index saturation versus optical depth of the bulk flow (BM) was first predicted by Titarchuk & Zannias (1998) and then it was subsequently reproduced in Monte Carlo simulations by Laurent & Titarchuk (1999), hereafter LT99.

It is worth noting that the index saturation effect is an intrinsic property of the BM onto a BH given that the spectral index = - 1 is a reciprocal of the Comptonization parameter Y (see this proof in ST09 and Bradshaw et al. 2007) which saturates when the BM optical depth, or M , increases. In fact, the Y-parameter is a product of the average photon energy exchange per scattering and the mean number of photon scattering Nsc , i.e., Y = Nsc . For the thermal Comptonization case, Y (4kT / me c2 ) 2 given that in this case = 4kT / me c2 and Nsc 2 for 1 (see, e.g., Rybicki & Lightman 1979) and, thus, the thermal Comptonization spectral index is [(4kT / me c2 ) 2 ]-1 . (5)

In the case of CF, the preferable direction for upscattered photons is the direction of BM onto the BH, i.e., along the radius. Note that the fractional photon energy change is E/E = (1 - 1 VR /c)/(1 - 2 VR /c), where 1 and 2 are the cosines of the angles between the direction of the electron velocity n = VR /VR and the


1478

TITARCHUK & SEIFINA

Vol. 706

Figure 12. Equivalent width of the iron line in eV as a function of flux density S15 GHz at 15 GHz (Ryle Telescope) in mJy for 1997/2005 data sets (blue/red points). Here we also include black points which have been recently found by Neilsen & Lee (2009) analyzing archival HETGS observations of GRS 1915+105 from the Chandra X-ray Observatory. (A color version of this figure is available in the online journal.)

Figure 11. BMC photon indices (top), normalization (middle), and X-ray flux (bottom) vs. flux density S15 GHz . Top and middle: red and black points stand for hard and soft spectral components, respectively, for the 1997­1998 active episode. Crimson and blue points correspond to hard and soft spectral components respectively for the 2005­2006 episode. Bottom: blue circles and red triangles stand for 1997­1998 and 2005­2006, respectively. (A color version of this figure is available in the online journal.)

direction of incoming and outcoming (scattered) photons, respectively. The number of scatterings of the up-Comptonized photons Nsc can be estimated as a ratio of the radial characteristic size of the CF L and the free path l in the direction of motion, namely Nsc L/ l = given that E/E has a maximum at 2 = 1for given 1 and VR . On the other hand, the efficiency per scattering for BM flow 1/ when 1 (Laurent & Titarchuk 2007, hereafter LT07). Hence for BMC, the Y-parameter does not depend on for high values of or dimensionless mass accretion rate m. Thus one can conclude that the Comptonization parameter Y = Nsc and hence the energy index = Y -1 saturate to a constant value when optical depth (or mass accretion rate) of the BM flow increases. However, the index saturation value is determined by the plasma temperature during a transition (see LT99). The plasma temperature strongly depends on the mass accretion rate in the BM region Mbm and its illumination by the disk photons L

Figure 13. Photon index plotted vs. BMC normalization for rise (left) and decay (right) transitions. Red/black points stand for hard/soft spectral components for the 1997 activity episode. Crimson and blue points correspond to hard and soft spectral components respectively for the 2005 episode. Points marked with rose oreol correspond to IS­HSS spectra fitted by the model which includes "high-temperature bbody" component (for details see Figure 4 and Table 7). (A color version of this figure is available in the online journal.)


No. 2, 2009

DISCOVERY OF PHOTON INDEX SATURATION IN THE BB BINARY GRS 1915+105

1479

In fact, the index should saturate with mass accretion rate in the CF mbm which is a sum of the disk mass accretion rate md and that in sub-Keplerian flow (CT95). Hence one can argue that this lower value of md at which the index saturates can be a sign of the presence of an extra (sub-Keplerian) component in the accretion flow onto a BH in GRS 1915+105. P. Laurent, & L. Titarchuk (2009, in preparation), hereafter LT09, studied the index­m correlation and also a modification of the disk BB spectrum due to Comptonization in the optically thick CC, which is formed due to the accumulation of accretion matter in the TL. They indeed explained the saturations of the indices of the soft and hard components of the resulting spectrum. Specifically, LT09 showed that gravitational energy of the accretion flow is released in the optically thick and relatively cold TL when the mass accretion rate md is higher than 1. The level of the index saturation depends on the radial velocity in the TL. LT09 also showed that the observable saturation index of the soft BMC component 2 4.2 can be reproduced in their Monte Carlo simulations for values of the TL radial velocity 0.05 c.
Figure 14. Photon index 1 (hard component) plotted vs. QPO centroid for the 1997 rise transition from GRS 1915+105. (A color version of this figure is available in the online journal.)

4.2. Physical Origin of "High-temperature BB-like" Component"? In eight IS-HSS spectra, we find observational evidence of the bump around 20 keV which can be fitted by a "4.5 keV BB-like" profile (see Figure 4 and Table 7). One can argue that this observable bump at 20 keV is a signature of the Compton reflection bump (see, e.g., Basko et al. 1974, ST80, CT95, and Magdziarz & Zdziarski 1995). But this interpretation encounters difficulties given that the hard power-law tails of these spectra are too steep to form the Compton bump. Indeed, ST80 and later LT07 demonstrated that the Compton bump as a result of a photon accumulation due to downscattering of hard photons

(see TLM98 and TF04). For higher Mbm and L the plasma temperature is lower. The level of index saturation decreases when the plasma temperature in the BM increases (TF04). Thus the index saturation levels can be different from source to source depending on the strength of the disk. Looking at Figure 13 one can also note that the index 1 starts its saturation at different values of BMC normalization ( md ) for different types of active episodes.

Figure 15. Evolution of flux density S15 GHz at 15 GHz (Ryle Telescope), RXTE/ASM count rate, QPO during 1997 (left column), and 2005 (right column) rise transitions of GRS 1915+105. Here QPO stands for the centroid frequency of the fundamental QPO. The left column panel demonstrates the presence of QPO even when for radio flux is low (<30 mJy). The right column panel shows an example of the presence of QPOs during prominent radio flux events (100­200 mJy). These panels demonstrate that the QPO appearances are independent of the radio flux during LHS­HSS transitions. (A color version of this figure is available in the online journal.)


1480

TITARCHUK & SEIFINA

Vol. 706

Figure 16. Top: evolution of the flux density S15 GHz at 15 GHz (Ryle Telescope), RXTE/ASM count rate during the 2005 rise transition. Top lower panel: red points A, B, and C mark moments at MJD = 53416, 53422, and 53442 (before, during, and after radio flares), respectively. Bottom: PDSs (left column) are plotted along with energy spectral diagram EF (E ) (right column) for three points A (top), B (middle), and C (bottom) of the X-ray light curve. Power density spectra for the soft (red, 3­15 keV) and hard (blue, 15­30 keV) energy bands are presented. There are QPOs at A and C points (A1, C1 panels) but there is none at B point (B1 panel), at the X-ray flare peak. Here data are denoted by red points; the spectral model presented with components is shown by blue, black, and dashed purple lines for BMC1, BMC2,and laor components, respectively. (A color version of this figure is available in the online journal.)

in the cold medium (for example, disk) cannot be produced if the photon index of the incident hard photon spectrum > 2. In fact, as one can see from Table 7, that in all spectra where we detect this 20 keV feature the index of the hard BMC component 1 > 2 (the indices vary between 2.5 and 3). In principle, this bump may also be a result of photoelectric absorption of the photons below 10 keV in the cold medium (disk) even if the incident spectrum is very steep. The photoelectric absorption cross section ph (7.8 keV/E )3 T , where E is the photon energy and T is the Thomson cross section (e.g., ` CT95). However, Laming & Titarchuk (2004) and Rozanska & Madej (2008) showed that the ionization of such a disk by the intensive X-ray radiation during IS­HSS invalidates the basic assumptions of the presence of the cold material in the innermost part of the source. Note that the hard tail of the X-ray spectrum in IS­HSS is formed in the CF, i.e., in the innermost part of the accretion flow, because we see the CF signature (the index saturation) when the source goes to IS­HSS.

On the other hand, Titarchuk (2002) argued that the specific spectral and timing features of X-ray radiation could be seen in BH sources only. Particularly, he stated that the photon­photon interaction of the effectively upscattered photons results in the powerful pair production near a BH horizon. Indeed, a large fraction of the upscatterred photons going inward are deflected by the relativistic free-fall electron in the outward direction (the aberration effect of light; see, e.g., Rybicki & Lightman 1979 and Appendix A). These diverted upscattered photons of energy Eup interact with incoming photons of energy Ein flux and ultimately this interaction leads to the pair creation if the condition Eup Ein (me c2 )2 is satisfied. Note that free-fall BM with the Lorentz factor 1 should be very close to the horizon, i.e., R 1 2, RS (6)


No. 2, 2009

DISCOVERY OF PHOTON INDEX SATURATION IN THE BB BINARY GRS 1915+105

1481

where R = R - RS is the radial distance to the horizon. Thus the created positrons extensively interact with accreting electrons there and therefore the annihilation line photons are created and distributed over the relatively narrow shell near the BH horizon. Specifically, R 3 â 104 10
2

m cm. 10

(7)

The proper energy (in the comoving flow frame) of annihilation line photons E511 should be seen by the Earth observer (in the zero frame) at the redshifted energy E0 = (1 - RS /R )
1/2

E511

E511 ,

(8)

where (1 - RS /R )1/2 (R/RS )1/2 = 1/ (see accordingly Equation (6)). In other words, the line energy displacement due to gravity as viewed by a far-away observer in free space is z + 1 = 1/(1 - RS /R )1/2 = . A significant fraction of these annihilation line photons strongly gravitationally redshifted can be directly seen by the Earth observer as a bump located at 20 keV and presumably related to the representative value of 20. P. Laurent & L. Titarchuk (2009, in preparation) made extensive Monte Carlo simulations of the X-ray spectral formation in the CF taking into account photon­electron, photon­photon, and pair­electron interactions. These simulations confirm our expectations that in some cases the emergent spectra of IS and HSS consist of the redshifted annihilation line located at 20 keV, which can be fitted by the "high-temperature BB-like" profile, and also the simulated spectra are extended to energies of the order of a few MeV (see Figure 4 and Grove et al. 1998, for details of IS­HSS spectral components). 4.3. Radio­X-ray Connection Migliari et al. (2005), hereafter MFK05, reported on correlations between radio luminosity and X-ray timing (QPO) features in X-ray binary systems containing low-magnetic-field neutron stars and BHs. The MFK05 conclusions on the radio­QPO correlation based on observations of seven neutron stars and one BH GX 339-4. For GX 339-4 they used data only in LHS before and after outburst. Titarchuk & Shaposhnikov (2005) based on the analysis of RXTE data from NS 4U1728-34 confirmed a correlation of X-ray and radio emissions with LFQPOs through all spectral states for this particular NS source. However, we do not find a real LHS in GRS 1915+105 for which the photon index 1 should be about 1.5, as that in GX 339-4 (see ST09), and therefore we cannot confirm or refute the radio­QPO correlation in GRS 1915+105 LHS similar to that found by MFK05 in GX 339-4. On the other hand, we find that LFQPOs do not correlate with the radio flux while they correlate with the hard-component photon index 1 through IS and HSS (see Figures 15 and 14, respectively). The absence of correlation between radio luminosity and QPO can be explained by the different origins of these quantities. While the QPO phenomenon is probably related to the TL oscillations (see Section 4.1), it is confirmed by the index­QPO correlation; the radio emission is presumably originated in the wind or wide open jet. Furthermore, because the radio flux and iron-line EW are strongly correlated (see Figure 12) one can conclude that the iron line is also formed

in the wind (see more on the line formation origin in LT07; Shaposhnikov et al. 2009; and Titarchuk et al. 2009). It is also worth noting that the X-ray flux and photon index do not correlate with the radio flux (see Figure 11). It can be explained by the different mechanisms of the energy releases in X-ray and radio. X-ray radiation is presumably formed in the innermost part of the disk and in the TL (or CC) while the radio emission is presumably formed in the jet or winds which are probably launched at the outskirts of the disk (see, e.g., Meyer-Homeister & Meyer 1999, 2003; Meyer et al. 2000; and TSA07). Using the aforementioned correlations of EW with radio flux we can suggest that some fraction of accretion flow may go to the outflow. Probably, the powerful outflow is launched at outer parts of the accretion disk and it is not by chance that we can see this strong correlation of the iron-line EW with the radio flux. Thus there are two ways for the matter to proceed: (1) in the outflow if the local mass accretion rate in the disk exceeds the critical value (which is proportional to radius; see TSA07), (2) in the disk where the matter proceeds and it ultimately converges onto a BH. This final stage of the accretion we observe as a saturation of index with mass accretion flow (CF signature). Whereas we do not find any correlation between radio and X-ray fluxes during the spectral transition, although, probably, we see an indication of X-ray­radio connection during a minor flare event. In Figure 16, we show the spectral and timing properties of X-ray emission of a typical minor X-ray/radio flare. As one can see from this figure (see the low panels there), the PDS and energy spectra are different at the peak of the flare from those before and after the flare. Specifically, the soft component is more pronounced but QPO features are not seen at the peak of the flare while they are present before and after the flare. Moreover, the flat part (white noise) of the peak PDS is extended to higher frequencies (break frequency b 10 Hz at the peak vs. b 2 Hz before and after the flare). In Section 3.3, we reported the results of the study of the energy dependence of the PDS shape and integrated power variability as a function of the photon energy. In Figure 16, we show two power spectra for two energy bands 2­15 keV (red) and 15­30 keV (blue). One can see that PDSs weakly depend on the energy band. We can suggest that the sizes of the photon emission areas LCC related to these two energy bands are also the same. In fact, Titarchuk et al. (2001) argued that QPO is proportional to the ratio of magnetoacoustic (plasma) velocity VMA and CC size LCC and hence one can conclude that the emission areas are the same because QPO are the same for these two energy bands. Lightman & Eardley (1974), hereafter LE74, suggested that the thin disk is always unstable in the inner region when radiation pressure dominates gas pressure. Using numerical simulations Lightman (1974) found that the innermost region of a disk around a BH is secular unstable against clumping of the gas into rings which observational appearances can be seen by the Earth observer as X-ray­radio flares. We can speculate that an increase of the soft and hard components in the observable X-ray spectrum at the flare peak (compare panels B2 and A2, and C2 of Figure 16) can be a sign that the radiation pressure gets to dominate in the inner disk region. On the other hand, the sign of the destruction of some part of the innermost part of the disk, as an effect of the high pressure instability, should be seen in the PDS. TSA07 argued that the break frequency b in the PDS is proportional to the diffusion frequency d = 1/tvisc /R 2 ^ where is the viscosity, tvisc is the viscous timescale, and R is ^ the radial size of the innermost part of the disk. Given that b


1482

TITARCHUK & SEIFINA

Vol. 706

increases at the peak with respect to that before and after the flare (compare panels B1 and A1, and C1 of Figure 16) it can imply that the size R decreases when b increases, i.e., some part of the innermost part of the disk is probably destroyed. In terms of the diffusion theory, this disk instability arises in the inner region where the viscous stress W is a decreasing function of the surface density and thus an effective diffusion coefficient of the nonlinear equation for becomes negative there (see Equations (4)­(5) in LE74). On the other hand, A. Makeev & L. Titarchuk (2009, in preparation), hereafter MT09, obtain this disk instability as a solution of the linear diffusion equation for , and they do not specify any (ad hoc) assumption about the nature of the disk viscosity (cf. LE74). They just assume the power-law viscosity distribution over the disk and they use the TLM98 angular velocity distribution in the TL. MT09 studied the TSA07 model of the PDS formation and they found that this model predicts the existence of the two distinct zones within the TL of the BH or neutron star, with oppositely different types of diffusion of perturbation taking place in each zone. A simple fact of the change of sign of the angular velocity derivative in the linear diffusion equation for at the critical radius Rmax results in a negative diffusion coefficient of the equation in the interval Rin < r < Rmax . Moreover, the change of sign of the diffusion coefficient at Rmax leads to a turnover of the angular momentum transfer, changing it toward the central body of the accreting system, instead of being pushed outward as in the Keplerian disk. One of the considered scenarios implies an unstable diffusion of the perturbations in the TL inner zone which might be an indication of the development of a X-ray flare followed by a radio flare. The absence of any QPOs in the TL and raising b are possible indications of this instability. 5. CONCLUSIONS We concentrate our efforts on the study of the correlation between the spectral index and the accretion disk luminosity. We argue that the shape of the correlation pattern can contain the direct BH signature. Namely, we show both observationally and theoretically that the index saturates with mass accretion rate which is a signature of a CF. This index saturation effect can exist only in BH sources. Also this correlation pattern carries the most direct information on the BH mass and the distance to the source (see ST09). We compiled the state transition data from GRS 1915+105 collected with the RXTE mission. We examined the correlation between the photon index of the Comptonized spectral component, its normalization, and the QPO frequency (see Figures 6­ 10 and 13­14). The spectral data of GRS 1915+105 are well fitted by two (soft and hard) BMC components for most of analyzed IS and HSS spectra (see Figure 3) while LHS spectra essentially require only one BMC component. In addition to two BMC components, eight IS-HSS spectra require an extra component which can be fitted by a "high-temperature BB-like" profile. We suggest this "BB" component is probably a signature of the redshifted annihilation line formed in the very narrow shell near a BH horizon due to high photon compactness taking place during intermediate and HSSs (see Figure 4 and Table 7). A remarkable result of our study is that the index­ normalization (mass accretion rate) correlation seen in GRS 1915+105 is predicted by the theory of the CF. We demonstrate that a strong index saturation versus disk flux seen in the index­disk flux correlation (see Figure 13) is an observational

signature of the presence of the CF, which should only exist in the BH sources. In other words, this index saturation effect provides robust observational evidence for the presence of BH in GRS 1915+105. We also find a tight positive correlation of QPO frequencies with the index (see Figure 14) and consequently that with the disk flux. Our comprehensive analysis of X-ray and radio emissions in GRS 1915+105 shows that QPOs are seen independently of radio activity of the source during the spectral transition from low-hard to HSS. Specifically, these QPO features have been detected at any level of the radio flux and even when the radio emission is at the noise level (see the left and right panels in Figure 15, correspondingly). We also do not find any correlation between X-ray and radio fluxes and X-ray power-law index (see Figure 11). However, we establish a strong correlation between EW of the iron line and radio flux in binary GRS 1915+105 (see Figure 12). We are grateful to the referee whose constructive suggestions help us to improve the paper quality. We acknowledge productive discussion with Nikolai Shaposhnikov and Ada Paizis, and we also thank Guy Pooley who kindly provided us Ryle Radio Telescope data. APPENDIX THE LIGHT ABERRATION EFFECT For a given electron velocity v moving inward along radius v = -v er the directions of the ray in the observed laboratory and comoving electron frames n0 and n are related to the aberration formula (see, e.g., Rybicki & Lightman 1979): cos = v/c - cos 0 , (v/c) cos 0 - 1 (A1)

where = [1 - (v/c)2 ]-1/2 is the Lorentz factor, cos 0 = -(n0 ·er ) and cos = -(n · er ) are for the laboratory (zero) and electron rest frames, respectively. For 1, we can write that 1 - v/c = 1/(2 2 ) and thus we can rewrite Equation (A1) as follows: cos -2sin2 (0 /2) + 1/(2 2 ) . 2sin2 (0 /2) + cos 0 /(2 2 ) (A2)

If 0 b/ then for b > 1 cos 1 - b2 < 0, 1+ b2 (A3)

i.e., the ray is directed against the electron direction. In fact, 1/ . On the other hand, cos -1 for all0 cos = -(n·er ) > 0 if b < 1. REFERENCES
Basko, M. M., Sunyaev,, & Titarchuk, L. G. 1974, A&A, 31, 249 ´ Belloni, T., Klein-Wolt, M., Mendez, M., van der Klis, M., & Paradijs, J. 2000, A&A, 355, 271 Belloni, T., et al. 2005, A&A, 440, 207 Bradshaw, C. F., Titarchuk, L., & Kuznetsov, S. I. 2007, ApJ, 663, 1225 Chakrabarti, S. K., & Titarchuk, L. 1995, ApJ, 455, 623 (CT95)

(A4)


No. 2, 2009

DISCOVERY OF PHOTON INDEX SATURATION IN THE BB BINARY GRS 1915+105

1483

Coppi, P. K. 1999, in ASP Conf. Ser. 161, High Energy Processes in Accreting Black Holes, ed. J. Poutanen & R. Svensson (San Francisco, CA: ASP), 375 Corbel, S., & Fender, R. P. 2002, ApJ, 573, L35 Fender, R. P. 2001, MNRAS, 322, 31 Fender, R. P., & Belloni, T. 2004, ARA&A, 42, 317 (FB04) Giannios, D. 2005, A&A, 437, 1007 (Paper III) Grove, J. E., et al. 1998, ApJ, 500, 899 Klein-Wolt, M., & van der Klis, M. 2008, ApJ, 675, 1407 Kotani, T., et al. 2000, ApJ, 539, 413 Kylafis, N. D., Papadakis, I. E., Reig, P., Giannios, D., & Pooley, G. G. 2008, A&A, 489, 481 Laming, J. M., & Titarchuk, L. 2004, ApJ, 615, L121 Laor, A. 1991, ApJ, 376, 90 Laurent, P., & Titarchuk, L. 1999, ApJ, 511, 289 (LT99) Laurent, P., & Titarchuk, L. 2007, ApJ, 656, 1056 (LT07) Lightman, A. P. 1974, ApJ, 194, 429 Lightman, A. P., & Eardley, D. M. 1974, ApJ, 187, L1 Magdziarz, P., & Zdziarski, A. A. 1995, MNRAS, 273, 837 (MZ95) Markoff, S., Falcke, H., & Fender, R. 2001, A&A, 372, L25 Markoff, S., Nowak, M., Corbel, S., Fender, R., & Falcke, H. 2003, A&A, 397, 645 Martocchia, A., Matt, G., Karas, V., Belloni, T., & Feroci, M. 2002, A&A, 387, 215 McClintock, J., & Remillard, R. 2006, in Compact Stellar X-ray Sources, ed. W. H. G. Lewin & M. van der Klis (Cambridge: Cambridge Univ. Press), 157 Meyer, F., Liu, B., & Meyer-Homeister, E. 2000, A&A, 361, 175 Meyer-Homeister, E., & Meyer, F. 1999, A&A, 348, 154 Meyer-Homeister, E., & Meyer, F. 2003, A&A, 402, 1013 Migliari, S., Fender, R. P., & van der Klis, M. 2005, MNRAS, 363, 112 Miller, J. M., & Homan, J. 2005, ApJ, 618, 107 Montanari, E., Titarchuk, L., & Frontera, F. 2009, ApJ, 692, 1597

Muno, M. P., Morgan, E. H., & Remillard, R. A. 1999, ApJ, 527, 321 Neilsen, J., & Lee, J. C. 2009, Nature, 458, 481 Pooley, G., & Fender, R. 1997, MNRAS, 292, 925 Reig, P., Belloni, T., van der Klis, M., & Mendez, M. 2000, ApJ, 541, 883 Remillard, R. A., & McClintock, J. E. 2006, ARA&A, 44, 49 Rodriguez, J., et al. 2008, ApJ, 676, 1436 ` Rozanska, A., & Madej, J. 2008, MNRAS, 386, 1872 Rybicki, G. B., & Lightman, A. P. 1979, Radiative Processes in Astrophysics (New York: Wiley Interscience) (RL79) Shaposhnikov, N., & Titarchuk, L. 2006, ApJ, 643, 1098 (ST06) Shaposhnikov, N., & Titarchuk, L. 2007, ApJ, 663, 445 (ST07) Shaposhnikov, N., & Titarchuk, L. 2009, ApJ, 699, 453 (ST09) Shaposhnikov, N., Titarchuk, L., & Laurent, P. 2009, ApJ, 699, 1223 Sunyaev, R. A., & Titarchuk, L. G. 1980, A&A, 86, 121 Swank, J. H. 1999, Nucl. Phys. B Proc. Suppl., 69, 12 Titarchuk, L. 2002, in Meeting Abstract 4th COSPAR Scientific Assembly, The Second World Space Congress Titarchuk, L., Bradshaw, C. F., & Wood, K. S. 2001, ApJ, 560, L55 Titarchuk, L. G., & Fiorito, R. 2004, ApJ, 612, 988 (TF04) Titarchuk, L., Lapidus, I. I., & Muslimov, A. 1998, ApJ, 499, 315 (TLM98) Titarchuk, L., Mastichiadis, A., & Kylafis, N. D. 1997, ApJ, 487, 834 Titarchuk, L., & Shaposhnikov, N. 2005, ApJ, 626, 298 Titarchuk, L., & Shaposhnikov, N. 2008, ApJ, 678, 1230 Titarchuk, L., Shaposhnikov, N., & Arefiev, V. 2007, ApJ, 660, 556 (TSA07) Titarchuk, L., Shaposhnikov, N., & Laurent, P. 2009, ApJ, 700, 1831 Titarchuk, L., & Zannias, T. 1998, ApJ, 493, 863 (TZ08) Trudolyubov, S. 2001, ApJ, 558, 276 Trudolyubov, S., Churazov, E., & Gilfanov, M. 1999, Astron. Lett., 25, 718 Vadawale, S. V., Rao, A. R., & Chakrabarti, S. K. 2001, A&A, 372, 793 Vignarca, F., Migliari, S., Belloni, T., Psaltis, D., & van der Klis, M. 2003, A&A, 397, 729 (V03)