Issue |
A&A
Volume 569, September 2014
|
|
---|---|---|
Article Number | A33 | |
Number of page(s) | 9 | |
Section | Stellar structure and evolution | |
DOI | https://doi.org/10.1051/0004-6361/201323294 | |
Published online | 15 September 2014 |
Research Note
The variability classes in IGR J17091–3624 during the 2011 outburst⋆
1
Key Laboratory for Particle Astrophysics, Institute of High Energy Physics,
Chinese Academy of Sciences (CAS),
19B Yuquan Road, 100049
Beijing
PR China
e-mail:
zhangzhen@ihep.ac.cn
2
University of Chinese Academy of Sciences,
19A Yuquan Road,
100049
Beijing, PR
China
3
Department of Astronomy, Beijing Normal University,
100875
Beijing, PR
China
4
Xinjiang Astronomical Observatory, CAS, 150, Science 1-Street, Urumqi, 830011
Xinjiang, PR
China
Received: 20 December 2013
Accepted: 24 July 2014
Aims. Based on the Rossi X-ray Timing Explorer (RXTE) data, we studied the typical variability classes in IGR J17091−3624 during the 2011 outburst and attempted to understand the physical mechanism that produces the highly structured variability.
Methods. Using the rate-resolved spectroscopy, we produced color-color diagrams (CCDs), hardness-intensity diagrams (HIDs), and rate-resolved energy spectra for each class. The time lags were calculated by means of the cross-power spectrum (CPS) analysis.
Results. Our study provides a comprehensive understanding of the ν, ρ, α, β/λ, μ and γ classes during this outburst. The source shows up with different track directions in the CCD loops and different types of the lag at ~1−6 s which are class dependent. It is notable that each variability class showing hard lags presents an anti-clockwise loop in the CCD, while each variability class showing soft lags presents a clockwise loop in the CCD. We discuss these observational results with the current models proposed for this type of system. The observed hardness ratio shows an increasing trend toward the high disk temperature for different classes, which can be regarded as evidence for thermal-viscous instability.
Key words: stars: individual: IGR J17091 / 3624 / X-rays: binaries / accretion, accretion disks / black hole physics
Tables 3 and 4 are available in electronic form at http://www.aanda.org
© ESO, 2014
1. Introduction
The object IGR J17091−3624 is the second black-hole X-ray binary after GRS 1915+105, which has shown a bewildering variety of X-ray light curve patterns (Altamirano et al. 2011). These patterns are usually called “variability classes” and characterized mostly with the large variabilities in amplitude and structure. Some variability classes have presented a loop-like behavior in their color-color diagrams (CCDs) and hardness-intensity diagrams (HIDs; Belloni et al. 1997a, 2000; Altamirano et al. 2011). The loop-like behaviors need to be studied to resolve the issue of whether there are differences in physical origin between different classes. Using the rate-resolved spectroscopy, Pahari et al. (2013a) extracted energy spectra during the peak and dip intervals of the different classes and found that their spectra are phenomenologically complex. We performed more refined spectral analysis for different classes in IGR J17091−3624. In IGR J17091−3624, the power density spectra (PDS) show a band limited noise component up to a break frequency νb. The frequency νb most plausibly corresponds to a viscous timescale (Ingram et al. 2009), over which flare may occur frequently (Belloni et al. 1997b). Pahari et al. (2013b) found that different types of lags of ~0.1−1.0 s were observed in different classes. The lag behaviors around νb provide a new perspective to our understanding the physical origins of the variability classes.
2. Observations and data analysis
Observational details of different classes used for the analysis.
We examined IGR J17091-3624 in its 2011 outburst using RXTE public archival observations and selected out six typical observations. The observational data was obtained with the Proportional Counter Array (PCA) instruments. These observations represent six variability classes (Altamirano et al. 2011, 2012). They are listed in Table 1. All the opened proportional counter units (PCUs) were used for the data analysis to improve statistics.
For each observation interval, we created a light curve with 1 s bin in three PHA channel intervals A = 0−14 (2.0−6.5 keV), B = 15−35 (6.5−15.0 keV), and C = 36−255 (15.0−60.0 keV). Using the rate resolved spectroscopy, we got the spectral parameters from different rate bands. Following Pahari et al. (2013a), we obtained the spectral parameters from the peak and dip, respectively. For those between the peak and dip, we obtained the spectral parameters from the ascending and descending parts, respectively. Flares were separated into bins of Δ ct/s wide intervals, where Δ is is a given number. From these data, we obtained a total count rate (R = A + B + C) and two hardness values: HR1 = B/A and HR2 = C/A. All the CCDs and HIDs were created by setting Δ = 50.
The PDS and light curves were produced from the PCA using standard techniques (e.g., Belloni et al. 2002; Altamirano et al. 2011). The PDS fitting was carried out in a standard way (Belloni et al. 2002). The errors were within 90% confidence level. The time lags between the high energy intervals (from E1 to E5) and the lowest interval (E0) were quantified by the CPS analysis (Cui et al. 1997). The absolute channels and energy intervals are: E0: 0−8 (2.06−4.09 keV), E1: 9−13 (4.09−6.12 keV), E2: 14−18 (6.12−8.17 keV), E3: 19−23 (8.17−10.22 keV), E4: 24−35 (10.22−15.18 keV) and E5: 36−57 (15.18−24.85 keV).
Rate-resolved energy spectra were extracted for each observation. All the energy spectra were fit in the 3.0−30.0 keV band using XSPEC version 12.7. The energy spectra have two main components: thermal and non-thermal (Belloni et al. 1997b; Pahari et al. 2013a). Different models can be used to well describe the two components (Pahari et al. 2013a). Following Pahari et al. (2013a), we fit these spectra with different disk models, such as diskbb or diskpbb in addition to the powerlaw model. Some spectra required a Gaussian for good fit. The errors were within 90% confidence level. To further improve statistic, we perform spectral analysis by setting Δ = 100. The hydrogen column density NH was fixed to 1.1 × 1022 cm-2 (Rodriguez et al. 2011; Krimm et al. 2011).
3. Results
![]() |
Fig. 1 HR1 vs. Tin for the six typical observations (Tables 3 and 4). Different symbols correspond to different observations. Each color stands for a different disk model. Purple: diskbb. Green: diskpbb. |
In Fig. 1, we plot HR1 versus inner disk temperature Tin. The data points of each class distribute diversely. However, there exists a global relationship between the HR1 and Tin for different classes (Fig. 1). With limited statistics, the HR1 appears to increase with Tin, especially when the Tin is higher than 1.2 keV. However, a similar trend is not observed in a HR2-Tin plot (Tables 3 and 4).
![]() |
Fig. 2 Left panels: CCDs extracted from the six typical observations. Insets show HID1s. Rates are given in per rate band with all the opened PCUs, 2−60 keV, and background subtracted. Middle panels: PDS. Shaded areas mark the frequency bands around νQPO and νmax. Right panels: time lags as a function of energy for various frequency bands vs. the lowest energy interval (2.06−4.09 keV). Different symbols correspond to different frequency bands. |
Summary of PDS and lag properties of the six observations.
Figure 2 shows the CCDs for flares in different classes, where the source does not always follow the same pattern. The source traces out an anti-clockwise loop for the ν and α observations and a clockwise loop for the μ, ρ, and β/λ observations. In Fig. 2, the insets show HIDs of the PCA count rate and the HR1 value (hereafter HID1s). All the HID1 loops are traced in the same direction with the CCD loops. In the γ observation, the source does not trace out a CCD loop. In Fig. 3, we plot PCU2 count rate versus HR2 (hereafter HID2s). A strong anticorrelation between PCU2 count rate and HR2 is observed for different classes and their evolutions follow the remarkable similar pattern. Each class shows a single elongated branch stretching out diagonally.
Figure 2 also shows representative PDS. Each one is dominated by a band-limited noise component, which breaks at a frequency νb, below which the spectrum is so strong that the fractional power density is about 0.1% per hertz (Fig. 2). In the PDS, one observes a ~10−100 mHz QPO (Q = νQPO/full width half maximum >2) from the “flares”. Each quasi-periodic oscillation (QPO) has a centroid frequency νQPO below or approximately equal to νb (Table 2). In the right panels of Fig. 2, we present the energy-lag spectra from the frequency bands around νQPO, νb, and νmax (Table 2). The QPOs with νQPO ≲ νb show obvious time lags. For the γ observation, even if there are no such QPO features visible in the PDS, we still performed the time lag analysis for the strong noise component at frequency below or approximately equal to νb.
In the ν and α observations, the ~10 mHz QPOs show hard lags of ~1−6 s (Fig. 2 and Table 2). In the μ, ρ, and β/λ observations, the ~10−100 mHz QPOs show soft lags of ~0.1−6 s (Fig. 2 and Table 2). In the μ observation, the time-lag spectra show soft lags of ~1−6 s in the frequency band from ~10 mHz to νb (Fig. 2). In the γ observation, the strong noise component does not show any hints of the time lag around 1−6 s (Fig. 2). It is notable that the observations observed with νQPO ≲ νb show obvious lags, while the one with no similar QPO feature behaves with a less prominent time lag (Table 2).
For different classes, the largest absolute value of the observed lags tend to occur in different frequency bands (see νmax values in Table 2). For a fixed frequency band, the time lags appear to vary between different energy intervals. As observed in Fig. 2, for example, the hard lags at νmax tend to increase with energy ranging from ~1 to ~6 s, while the soft lags at νmax tend to decrease with energy ranging from ~−1 to ~−6 s. It is clear that the time lags of ~1−6 s depend strongly on photon energy. The higher the energy, the larger the lag timescale (with the exception of the 45 mHz QPO in the ρ observation, where the lag timescale increases at first and then drops above 12 keV).
We found a strong relation between the CCD and the time lag: the observations showing hard lag features present a CCD loop traced out in an anti-clockwise direction while the observations showing soft lag features present a CCD loop traced out in a clockwise direction. It suggests an intrinsic relation between the track direction of the CCD loop and the sign of the time lag.
![]() |
Fig. 3 HIDs of the PCU2 count rate and the HR2 value. Different symbols correspond to different observations. The anticorrelation is well described by a power law of the form y = ax− b + c, where a = 30.96 ± 16.98, b = 1.05 ± 0.25, c = 13.22 ± 29.37, and |
4. Discussion and conclusions
Belloni et al. (1997b) interpreted the large variability of the X-ray flux of GRS 1915+105 with the rapid removal and replenishment of matter, which forms a hot accretion flow and is driven by the thermal-viscous instability. Since IGR J17091−3624 and GRS 1915+105 were observed with similar behavior, it is reasonable to speculate that both have a similar physics underneath (Capitanio et al. 2012). In this model, the source evolves in form of loops between the lower and upper branches of the S-shaped Ṁ(Σ)-curve (Belloni et al. 1997b). Each loop refers to a single cycle of instability (Janiuk et al. 2005) and, as suggested in Belloni et al. (1997b), should show up in CCD. In IGR J17091−3624, several classes were observed with the loop-like behaviors in the CCD and HID1, which are also due to the cycles of instability. We find for the first time that the track directions of their CCD and HID1 loops are highly class dependent, which is in turn suggestive of having the different intrinsic physical origins.
For different classes, the HR1 shows an increasing trend toward high Tin. This global trend suggests that these classes may share with a common physical mechanism. If the variability classes are caused by thermal-viscous instability, it is usually thought that the disk is truncated at a region where Prad is balanced with Pgas (Belloni et al. 1997b). Here Prad is the radiation pressure and Pgas is the gas pressure. This region is the innermost part of an accretion flow, according to the description by a standard disk model. Thus, the inner disk temperature can be used to describe the region where Prad ~ Pgas. When the luminosity increases, the disk moves inwards and the hot accretion flow extends toward outside to strengthen the corona, resulting in a harder spectrum and hence a larger hardness ratio (Lei et al. 2013). With the Prad increasing, Pgas increases because Pgas ~ Prad. Then Tin increases with the increasing Pgas(=nkTin). This trend can be regarded as an evidence for the thermal-viscous instability.
A similar pattern in the HID2 for the different classes is shown, and those classes showing HID1 loops present no loop in their HID2s. One interpretation is that the thermal-viscous instability become less prominent in the non-thermal dominant inner hot accretion flow with Prad ≫ Pgas. Due to the same reason, the HR2 shows no trend of increasing toward higher Tin for different classes.
In IGR J17091−3624, the observed lags are of ~1−6 s. The timescale of only reflection or the Comptonization process is expected to be about ≲1 s (Lei et al. 2013; Pahari et al. 2013b). Therefore, these lags should be related to other process. As mentioned in Janiuk et al. (2005), the hard lag of ~1−6 s is comparable to the viscous timescale in the hot accretion flow. Such a hard lag is required by the hot accretion flow to adjust to the changing conditions in the disk. However, the model reported in Janiuk et al. (2005) did not explain both energy dependence of the lag and the soft lag. The viscous fluctuation model is also a candidate for an alternative interpretation. In this model, the viscous fluctuations come from outer regions (Pahari et al. 2013b). Since the outer regions typically emits lower energy photons than the inner ones, the expected lag should be at the hard band (Pahari et al. 2013b). When the viscous fluctuations occur at outer regions and propagate inwards, it causes an accretion rate fluctuation in the hot accretion flow, and the expected hard lag from this model should correlate with the energy in a positive manner. However, this viscous propagation model cannot be used to explain the soft lag of ~1 s (Pahari et al. 2013b). Such a soft lag can be explained when strong disturbance (such as hot gas steams) is produced from the hot accretion flow with Prad ~ Pgas (Frank et al. 1987) and propagates to the outer and optically thick disk, hence affecting the soft X-ray emission region. In this scenario, the disturbances come from inner regions, which is different from the viscous fluctuation model. The disturbance from the optically thin hot accretion flow contributes to the decreasing of the scattering depth, so that more photons with lower energies can escape from the outer disk. These soft photons escaping from the outer region have a time lag with respect to those hard photons escaping from the inner region. This disturbance from the hot accretion flow perturbs the outer material and then produces the observed soft lag. From this point of view, the soft lag of 1−6 s can be regarded as the time taken for the disk to adjust to the disturbance in the hot accretion flow. While the disturbance propagation from the hot accretion flow could manifest itself in a lag observationally, any radial movement of the propagation would cause the lag timescale to increase with energy (Cui et al. 1997). However, it remains unclear how to associate the lag signs with the track directions followed along the CCD.
Online material
Spectral parameters for the different rate bands from the six typical RXTE observations, which are fit with diskbb+powerlaw model.
Spectral parameters for the different rate bands from the six typical RXTE observations, which are fit with diskpbb+powerlaw model.
Acknowledgments
We are thankful to the anonymous referee for providing useful comments/suggestions which improve the paper. This work is supported in part by the National Natural Science Foundation of China (NSFC), XTP project XDA04060604, NSFC-11173024, NSFC-11133002, NSFC-11073021, the NSFC for Young Scientists, grant No. 11203064, and the Strategic Priority Research Program “The Emergence of Cosmological Structures” of the Chinese Academy of Sciences, grant No. XDB09000000.
References
- Altamirano, D., & Belloni, T. 2012, ApJ, 747, L4 [NASA ADS] [CrossRef] [Google Scholar]
- Altamirano, D., Belloni, T., Linares, M., et al. 2011, ApJ, 742, L17 [NASA ADS] [CrossRef] [Google Scholar]
- Belloni, T., Méndez, M., King, A. R., van der Klis, M., & van Paradijs, J. 1997a, ApJ, 488, L109 [NASA ADS] [CrossRef] [Google Scholar]
- Belloni, T., Méndez, M., King, A. R., van der Klis, M., & van Paradijs, J. 1997b, ApJ, 479, L145 [NASA ADS] [CrossRef] [Google Scholar]
- Belloni, T., Psaltis, D., & van der Klis, M. 2000, A&A, 355, 271 [NASA ADS] [Google Scholar]
- Belloni, T., Psaltis, D. & van der Kils, M. 2002, ApJ, 572, 392 [NASA ADS] [CrossRef] [Google Scholar]
- Capitanio, F., Del Santo, M., Bozzo, E., et al. 2012, MNRAS, 422, 3130 [NASA ADS] [CrossRef] [Google Scholar]
- Cui, W. 1999, ApJ, 524, L59 [NASA ADS] [CrossRef] [Google Scholar]
- Cui, W., Zhang, S. N., Fock, W., & Swank, J. H. 1997, ApJ, 484, 383 [NASA ADS] [CrossRef] [Google Scholar]
- Frank, J., King, A. R., Lasota, J. P. 1987, A&A, 178, 137 [NASA ADS] [Google Scholar]
- Ingram, A., Done, C., & Fragile, C. P. 2009, MNRAS, 397, L101 [NASA ADS] [CrossRef] [Google Scholar]
- Janiuk, A., & Czerny, B. 2005, MNRAS, 356, 205 [NASA ADS] [CrossRef] [Google Scholar]
- Krimm, H. A., Barthelmy, S. D., Baumgartner, W., et al. 2011, ATel, 3144 [Google Scholar]
- Lei, Y. J., Zhang, H. T., & Zhang, C. M. 2013, ApJ, 146, 60 [NASA ADS] [CrossRef] [Google Scholar]
- Pahari, M., Bhattacharyya, S., Yadav, J. S., & Pandey, S. K. 2012, MNRAS, 422, L87 [NASA ADS] [CrossRef] [Google Scholar]
- Pahari, M., Misra, R., Mukherjee, A., Yadav, J. S., & Pandey, S. K. 2013a, MNRAS, 436, 2334 [NASA ADS] [CrossRef] [Google Scholar]
- Pahari, M., Yadav, J. S., Rodriguez, J., et al. 2013b, ApJ, 778, 46 [NASA ADS] [CrossRef] [Google Scholar]
- Pahari, M., Yadav, R., & Bhattacharyya, S. 2014, ApJ, 783, 141 [NASA ADS] [CrossRef] [Google Scholar]
- Rodriguez, J., Corbel, S., Caballero, I., et al. 2011, A&A, 533, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
Spectral parameters for the different rate bands from the six typical RXTE observations, which are fit with diskbb+powerlaw model.
Spectral parameters for the different rate bands from the six typical RXTE observations, which are fit with diskpbb+powerlaw model.
All Figures
![]() |
Fig. 1 HR1 vs. Tin for the six typical observations (Tables 3 and 4). Different symbols correspond to different observations. Each color stands for a different disk model. Purple: diskbb. Green: diskpbb. |
In the text |
![]() |
Fig. 2 Left panels: CCDs extracted from the six typical observations. Insets show HID1s. Rates are given in per rate band with all the opened PCUs, 2−60 keV, and background subtracted. Middle panels: PDS. Shaded areas mark the frequency bands around νQPO and νmax. Right panels: time lags as a function of energy for various frequency bands vs. the lowest energy interval (2.06−4.09 keV). Different symbols correspond to different frequency bands. |
In the text |
![]() |
Fig. 3 HIDs of the PCU2 count rate and the HR2 value. Different symbols correspond to different observations. The anticorrelation is well described by a power law of the form y = ax− b + c, where a = 30.96 ± 16.98, b = 1.05 ± 0.25, c = 13.22 ± 29.37, and |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.