Accretion Flow Properties of Swift J1753.5-0127 during its 2005 outburst
Galactic X-ray binary black hole candidate Swift J1753.5-0127 was discovered on June 30 2005 by Swift/BAT instrument. In this paper, we make detailed analysis of spectral and timing properties of its 2005 outburst using RXTE/PCA archival data. We study evolution of spectral properties of the source from spectral analysis with the additive table model fits file of the Chakrabarti-Titarchuk two-components advective flow (TCAF) solution. From spectral fit, we extract physical flow parameters, such as, Keplerian disk accretion rate, sub-Keplerian halo rate, shock location and shock compression ratio, etc. We also study the evolution of temporal properties, such as observation of low frequency quasi-periodic oscillations (QPOs), variation of X-ray intensity throughout the outburst. From the nature of the variation of QPOs, and accretion rate ratios (ARRs=ratio of halo to disk rates), we classify entire 2005 outburst into two harder (hard-intermediate and hard) spectral states. No signature of softer (soft-intermediate and soft) spectral states are seen. This may be because of significant halo rate throughout the outburst. This behavior is similar to a class of other short orbital period sources, such as, MAXI J1836-194, MAXI J1659-152 and XTE J1118+480. Here, we also estimate probable mass range of the source, to be in between to based on our spectral analysis.
Galactic black holes candidates (BHCs) are situated generally in binary systems. Some of these sources are transient in nature. These black hole X-ray binaries (BHXRB) spend most of their lifetime in the quiescent states. Occasionally, they show outbursts in which high energy radiation intensity increases by a factor of hundred or more. A possible explanation is that an outburst is triggered by the sudden rise in viscosity in the accreting matter coming from the companion (Chakrabarti & Titarchuk 1995; Ebisawa et al. 1996). The mass transfer from the companion to the compact object can occur via Roche lobe overflow or by wind accretion depending on the nature of the companion. In general, if the companion is a low mass K or M type star, Roche lobe overflow is expected. If the companion is a mass losing high massive star, such as a type of O, A or B, winds form and accretion of winds are expected. In a BHXRB, the matter does not fall radially on the compact object. Rather, it forms an inward spiraling accretion disk around it. This disk emits electromagnetic radiation in the range from radio to -rays. High energy radiations, such as X-rays, are very important to study the properties of black holes (BHs), as they come from regions close to the BH event horizon and therefore carry accurate information about the compact object.
In general, X-ray spectrum consists of a multi-color black body type component of soft photons and a power-law (PL) type tail of hard photons. The multi-color black body component is believed to come from the standard accretion disk (Shakura & Sunyaev 1973; Novikov & Throne 1973) while the PL component is believed to come from a hot-electron filled Compton cloud region (Sunyaev & Titarchuk 1980, 1985). In the two component advective flow (TCAF) solution Chakrabarti and his collaborators (Chakrabarti & Titarchuk, 1995, hereafter CT95; Chakrabarti, 1997), the so-called Compton cloud is replaced with CENtrifugal pressure supported BOundary Layer (CENBOL). In the TCAF solution, the accretion disk consists of two components of accreting matter - a standard Keplerian disk with temperature modified by the Compton cloud and a sub-Keplerian halo. The Keplerian disk is a high viscous disk with high angular momentum while the sub-Keplerian disk has low viscosity and low angular momentum. The sub-Keplerian halo piles up at the centrifugal barrier and an axisymmetric shock is formed (Chakrabarti 1990). The soft photons from the Keplerian disk are inverse-Comptonized at the post-shock region and become hard photons. CENBOL is also considered to be the base of jet (Chakrabarti 1999). Its oscillation is believed to cause low frequency quasi-Periodic Oscillations (QPOs). QPOs are observed if resonance conditions are satisfied when infall timescale roughly matches with the cooling time scale (Molteni et al. 1996; Chakrabarti et al. 2015) or the Rankine-Hugoniot conditions are not satisfied (Ryu et al. 1997). The observed frequency of the QPOs are inversely proportional to the infall time in the post-shock region (Chakrabarti & Manickam 2000).
Recently, Debnath and his collaborators’ have successfully included this most generalized accretion solution into HEASARC’s spectral analysis software package XSPEC (Arnaud 1996) as an additive table model (Debnath et al. 2014; Mondal et al. 2014, hereafter MDC14; Debnath et al. 2015a,b, hereafter DMC15, DMCM15 respectively; Jana et al. 2016, hereafter JDCMM16; Chatterjee et al. 2016, hereafter CDCMJ16) to find the nature of accretion flow properties of a few transient BH sources during their X-ray outbursts. Apart from the mass of the black hole and constant normalization, one requires four parameters (two types of accretion rates, shock location and compression ratio) to fit a BH spectrum with the TCAF solution in XSPEC. Transient BHXRBs show several spectral states during an entire outburst, which could be well understood by observing variations of these physical flow parameters and nature of QPOs. Four spectral states, namely, hard (HS), hard-intermediate (HIMS), soft intermediate (SIMS) and soft state (SS) are commonly observed in an outburst of a transient BHC. In a typical outburst, BHCs generally stay in the HS at the start of an outburst. After that, it enters into the HIMS, SIMS and SS one by one. The shock moves closer to the black hole as it moves to the softer state and before disappearing altogether in soft/very soft state. The declining phase starts when the matter supply to the inner disk is turned off. During this phase of an outburst, the shock moves away from the black hole and the spectra start to become harder. Generally SIMS, HIMS and HS are observed one after another, to complete the hysteresis loop of the spectral states (see, Debnath et al. 2013 and references therein).
In general, in the harder states (HS and HIMS), the halo rate dominates while in the softer states (SIMS and SS), the disk rate dominates. From the spectral fits of data of each day with TCAF solution, we can also estimate the mass of black hole more accurately using constant normalization method (Molla et al. 2016; hereafter MDCMJ16). According to the TCAF solution, model normalization (N) should not vary on a daily basis, since it only depends on the mass, distance and disk inclination angle of the disk, unless precession in the disk occurs to change the projected surface of the emission area or there are some significant outflow activities which are not included in current TCAF model fits file. This method has been successfully applied to estimate masses of few transient Galactic BHCs, such as, MAXI J1659-152 (MDCMJ16), MAXI J1836-194 (JDCMM16), MAXI J1543-564 (CDCMJ16), H 1743-322 (Bhattacharjee et al. 2017; Molla et al. 2017). This motivated us to estimate the mass of Swift J1753.5-0127 using the same method. Moreover, Mondal et al. (2016) also estimated spin parameter of BHC GX 339-4, from their spectral study of the relativistic broadening of the double Fe-line with combined TCAF plus LAOR model fits using Swift/XRT and NuSTAR/FMA data of the 2013 outburst.
Swift J1753.5-0127 was discovered by Swift/BAT on 2005 June 30 (Palmar et al. 2005) at the sky position of RA, DEC. It was also detected in optical (Halpern et. al 2005) and radio (Fender et al. 2005). BHC Swift J1753.5-0127 has a short orbital period of hrs (Zurita et al. 2007). However Neustroev et al. (2014) suggested the orbital period is hrs. They also reported that the binary has a primary BH mass of less than and companion as and disk inclination angle of . Recently, Shaw et al. (2016) predicted the mass of central BH to be greater than . Cadolle Bel et al. (2007) estimated the distance of the system is likely to be in between kpc. However, from UV spectral study, Froning et al. (2014) calculated the source distance as kpc for inclination of and kpc for inclination of . The secondary companion is also reported to be a K or M type main-sequence star (Bell et. al 2007). Radio jet was also observed during the 2005 outburst of source (Fender et al. 2005; Soleri et al. 2010). The source is still active in HS with observable low luminosity X-ray intensities. In fact, the source has not gone into the quiescent state, since its discovery. The source was active over 11 years before going to quiescence state in November, 2016 (Shaw et al. 2016; Plotkin et al. 2016) for months, before showing new flaring activity in February, 2017 (Kong et al. 2017; Qasim et al. 2017).
Recent studies of different transient BHCs show that outbursts are of two types: type-I or classical type, where all spectral states (HS SS via two intermediate states, HIMS and SIMS) are observed, and type-II or harder type, where SS (sometimes even SIMS) are absent. The later type of outbursts are generally termed as ‘failed’ outbursts. Our recent study of few type-II outbursts of few transient BHCs, such as 2010 outburst of MAXI J1659-152 (Debnath et al. 2015b), 2011 outburst of MAXI J1836-194 (Jana et al. 2016), 2000 outburst of XTE J1118+480 (Chatterjee et al. 2017) confirms that short orbital period BHCs generally show these class of outbursts.
The is organized in following way. In §2, we briefly discuss observation and data analysis procedure. In §3, we present TCAF model fitted spectral analysis result. There we also present timing analysis results. We also estimate mass of the BHC Swift J1753.5-0127 based on TCAF model fitted constant normalization method. Finally in §4, a brief discussions and concluding remarks are presented.
2 Observation and Data Analysis
RXTE/PCA monitored the source two days after its discovery on a daily basis (Morgan et al. 2005) Here we study archival data of observational IDs starting from the first PCA observed day; 2005 July 2 (Modified Julian Day, i.e., MJD=53553.06) to 2005 October 19 (MJD=53662.77) using XSPEC version 12.8. To analyze the data we follow standard data analysis technique of the RXTE/PCA instrument as presented in Debnath et al. (2013, 2015a).
For timing analysis we use the PCA Science Binned mode (FS3f*.gz) data with a maximum timing resolution of to generate light curves for well-calibrated Proportional Counter Unit 2 (PCU2; including all six layers) in keV (0-59 channels) and keV (0-35 channels). To generate the power-density spectra (PDS), “powspec” routine of the XRONOS package is used to compute rms fractional variability on keV light curves of sec time bin. To find centroid frequencies of the observed QPOs () in PDS, we fit each QPO with a Lorentzian profile and use the “fit err” command to get ‘’ error limits.
For spectral analysis we use Standard2 mode Science Data (FS4a*.gz) of the PCA instrument. The keV background subtracted PCU2 spectra are fitted with both the TCAF-based model fits file and combined disk black body (DBB) and PL model components in XSPEC. However, the spectra of the last eight observations, were fitted only with PL, as no significant DBB component was needed. To achieve the best spectral fits, a Gaussian line of peak energy around keV (iron-line emission) is used. Hydrogen column density (N) was kept frozen at 1.010 atoms cm (Kennea et al., 2011) for absorption model wabs. For the entire outburst, we also use a fixed % systematic instrumental error for the spectral analysis. The XSPEC command ‘err’ is used to find 90% confidence ‘’ error values for the model fitted parameters after achieving the best fit based on reduced chi-square value ().
To fit spectra using the TCAF-based model additive table fits file, one needs to supply five model input parameters such as, black hole mass () in solar mass () unit, Keplerian accretion rate ( in Eddington rate , where is the Eddington luminosity and is the velocity of light), sub-Keplerian accretion rate ( in ), location of the shock ( in Schwarzschild radius =) and the ) compression ratio. (R=, where and are post- and pre-shock densities respectively) of the shock. The normalization (N), though simultaneously determined is a constant and once it is determined for one outburst, it remains the same, though it is found to fluctuate within a narrow region due to errors in fits and possible inaccuracies in data. Thus, once the mass () and the normalization (N) are known, only four parameters are needed to fit the data unless there are precession or jet activity. The size of the hot Compton cloud, i.e., CENBOL and its nature (such as, height, temperature, optical depth, etc.) could be calculated from the shock parameters.
Recent studies by our group show that intricacies of accretion dynamics around BHC become transparent when the data is analyzed with the TCAF solution. This motivated us to study the accretion flow properties of the 2005 outburst of the well known BHC Swift J1753.5-0127 under TCAF paradigm. Here, we study observations between 2005 July 2 (MJD=53553.06) and 2005 October 19 (MJD=53662.77) RXTE/PCA data for our analysis. We use RXTE/PCU2 data for our spectral and timing analysis. To extract physical accretion flow parameters during the outburst of the source, keV spectra are fitted with the current version (v0.3) of the TCAF model fits file. We also compare our results with the spectral analysis result using DBB plus PL model. Combined DBB plus PL model provides only gross properties of accretion disk such as disk temperature, photon index and flux from thermal and non-thermal components, but real physical reasons behind the observed spectra come by fitting with TCAF solution. From the spectral fits, we estimate the mass of the BHC, since it is also an input parameter of the present TCAF model fits file.
3.1 Spectral Analysis with the TCAF Solution and with Combined DBB and PL models
We use keV RXTE/PCU2 data for spectral analysis of BHC Swift J1753.5-0127 during its 2005 outburst. Disk black body temperature () and photon index () are obtained from the combined DBB plus PL model fit of the spectra. We also calculated individual model component fluxes (DBB and PL separately) in keV energy band to have a rough estimate about the photon contribution from the thermal (DBB flux) and non-thermal (PL flux) processes around the BH during the 2005 outburst. Spectral fits with the current version of the TCAF solution fits file provide us with two component (Keplerian disk rate and sub-Keplerian halo rate ) accretion rates and shock parameters (location and compression ratio ). These parameters, in turn, provide us with other physical quantities such as the viscosity, temperature, size and the optical depth of the Compton cloud.
In Fig. 1a, we show keV PCU2 count rate variation with day (MJD), which seems to be a sharp rise and exponential decay type. In Figs. 1b & 1d, variation of TCAF model fitted two component accretion rates (Keplerian disk rate , and sub-Keplerian halo rate ) are shown. To compare in Fig. 1c & 1e, variation of individual model components of the combined DBB and PL model fitted spectra are shown.
In Figs. 2a & 2b, we show variations of the combined DBB plus PL model fitted DBB temperature ( in keV) and PL photon index () respectively. TCAF model fitted ratio between the sub-Keplerian halo rate and Keplerian disk rate, i.e., accretion rate ratio (ARR = ) are shown in Fig. 2c. In Fig. 2d, we show the evolution of the observed (primary dominating) QPO frequencies with day. In Fig. 3(a-d), we show variations of the TCAF model fitted shock location (), compression ratio (R), normalization (N), and the derived mass.
Depending on the nature of the variations of ARR, and , we classify the entire outburst into two harder spectral states in both the rising (Ris.) and declining (Dec.) phases, although many previous authors reported the entire outburst phase as the low HS (see, Cadolle Bel et al. 2007; Ramadevi & Seetha 2007). The spectral states are also observed in the sequence: HS (Ris.) HIMS (Ris.) HIMS (Dec.) HS (Dec.). In Fig. 4(a-d), as examples of our fits, TCAF fitted spectra with residuals for 4 observations selected from four different states are shown. The observation Ids are: 91094-01-01-00 (MJD=53553.05), 91094-01-01-03 (MJD=53556.19), 91423-01-02-06 (MJD=53564.91), and 91423-01-09-00 (MJD=53612.75) are from HS (Ris.), HIMS (Ris.), HIMS (Dec.), and HS (Dec.) respectively. In Fig. 5(a-d), we show unabsorbed theoretical TCAF model spectra (online solid black curves) with its two components i.e., blackbody (online dashed blue curves) and Comptonized (online dot-dashed red curves), which are used to fit four different state spectra shown in Fig. 4. Photons in the multi-color blackbody type spectral component are mainly originated from the Keplerian disk, where as photons in the Comptonized part originate from the CENBOL. Since the Keplerian disk in the pre-shock region is also heated up by absorbing hard radiations from the CENBOL, it is not exactly the standard Keplerian disk as explained in CT95. Similarly, the so-called ‘power-law’ component is not really a power-law, since it has cut-off due to recoil effect and a possible component reflected from the CENBOL to the observer (built-it in TCAF procedure). The number of photons under the PL component is those from the pre-shock Keplerian disk intercepted by the CENBOL. Thus, in effect, both the components are inter-linked.
3.2 Evolution of Timings Properties during the Outburst
Evolution of temporal properties of this BHC has been discussed by several authors (Zhang et al. 2007; Ramadevi & Seetha 2007). From the outbursting profile as shown in Fig. 1a, we infer that this outburst is a Fast Rise and Slow Decay (FRSD) type rather than a slow rise slow Decay (SRSD; see Debnath et al. 2010). Low frequency QPOs are observed almost every day of the 2005 outburst of Swift J1753.5-0127. In the rising phase of the outburst, QPOs are observed on all observations, where as in the declining phase, QPOs are observed in 28 observations out of total observations. Most of the observed QPOs are strong type-C QPOs.
A sharp rise in X-ray intensity (see, Fig. 1a) and QPO frequencies (from to Hz) are observed in the initial two days, from 2005 July 2 (MJD=53553.06) to 2005 July 4 (MJD=53555.21), then slow rise until 5th observation i.e., 2005 July 6 (MJD=53557.24) of the outburst. After that, both PCU2 count rates and QPO frequencies are showed a general trend of slow decreasing nature till the end of our observations. 2005 September 17 (MJD=53630.31) is the last day during the outburst, where a prominent QPO of Hz is observed.
3.3 Evolution of Spectral Properties during the Outburst
We use keV RXTE/PCU2 data for the spectral analysis. With TCAF solution, we extract physical parameters such as two types of accretion rates ( and ) in Eddington accretion rate unit, shock location () in () unit and compression ratio () from each spectral fit. From combined DBB and PL model fitted spectra, we obtain the value of PL photon index (), DBB temperatures ( in keV). Based on the variations of the ARRs (=) and nature (shape, frequency, value, rms% etc) of QPOs, we find only two harder (HS and HIMS) spectral states during the 2005 outburst of the BHC Swift J1753.5-0127. Similar to 2011 outburst of MAXI J1836-194 (Jana et al. 2016), these spectral states are observed in the sequence of : HS (Ris.) HIMS (Ris.) HIMS (Dec.) HS (Dec.). The nature of the source during observed spectral states are discussed below.
3.3.1 Hard State - Rising
Initial two observations (from 2005 July 2 to 4; MJD=53553.06 to 53555.21) belong to this spectral state. During this phase of days, a sharp rise in the TCAF model fitted both accretion rates ( and ) as well as individual model flux components of the combined DBB and PL model fits are observed. QPO frequency also rises sharply during these two observations. Similar to other transient BHCs (e.g., H 1743-322, GX 339-4, MAXI J1659-152, MAXI J1836-194), recently studied by our group with the TCAF solution (see, Mondal et al. 2014; Debnath et al. 2015a,b; Jana et al. 2016), a local maximum in ARR values is also observed on the 2nd observation of the current outburst of Swift J1753.5-0127. So, we define this observation as the transition from HS (Ris.) to HIMS (Ris.).
3.3.2 Hard-Intermediate State - Rising
The source has been observed to this spectral state for the next two days (from 2005 July 4 to 6; MJD=53555.21 to 53557.24), three observations. During this phase of the outburst, a slow rise in both disk and halo accretion rates are observed, although rise in the disk rate is slightly higher than halo rate. As a result of the ARR values are observed to decrease slowly, and becomes lowest on the 3rd observation day (MJD=53557.24). We define this day as the transition day from HIMS (Ris.) to HIMS (Dec.). After that, X-ray intensities as well as both the components of accretion rates start to decrease. Also, on this day, maximum values of the PL photon indices (=) and QPO frequencies (= Hz) during entire outburst are observed.
3.3.3 Hard-Intermediate State - Declining
The source was observed in this spectral state for around next days. During this phase of the outburst, X-ray intensity as well as TCAF model fitted two component accretion rates decrease monotonically, although there is a rise in ARR values. This indicates that the rate of fall in Keplerian disk rate is faster than the sub-Keplerian halo rates, which implies that the spectrum becomes harder as day progresses. This hardening of the spectra also could be concluded from the variation of the PL indices. A general decreasing nature of the values are observed during this phase of the outburst.
The shock is found to move slowly away from the BH, while compression ratio (R) stays roughly at a constant value of . QPOs are also observed to decrease with time (day). We define observation on the 2005 August 7 (MJD=53589.49) as the transition day from HIMS (Dec.) to HS (Dec.). On this particular observation, the shock becomes stronger with . After that, ARRs, QPO frequencies, PL indices became roughly constant and the trend continued till the end of our observation. The shock started to recede back with an acceleration with progressively increasing shock strength after the transition day.
3.3.4 Hard State - Declining
The source is found in this spectral state till the end of our observation starting from the HIMS (Dec.) to HS (Dec.) transition day (MJD=53589.49). During this phase of the outburst, both TCAF model fitted accretion rates decrease roughly in the same rate due to lack of supply of matter from the companion. As a result of that, ARRs are observed roughly at a constant value of . QPOs are also observed roughly at a constant frequency (within a short range of Hz), although they have not been detected for all observations during this phase of the outburst. During the initial period of this phase, the shock moves away from the BH with increasing shock strength, before it settles down into a roughly constant value starting from 2005 August 25 (MJD=53607.25). The outburst continued for a long duration even after our observations with low X-ray activity.
3.4 Estimation of the Mass of the Black Hole
Recently, MDCMJ16 used two independent methods, such as, TCAF model fitted constant normalization (N), and propagating oscillatory model (POS) model (see, Chakrabarti et al. 2008; Debnath et al. 2010, 2013; Nandi et al. 2012) fitted QPO frequency evolution, to measure unknown mass of BHC MAXI J1659-152. Furthermore, these methods are successfully applied to measure mass of other BHCs of unknown mass, such as, MAXI J1836-194 (JDCMM16), MAXI J1543-564 (CDCMJ16), H 1743-322 (Bhattacharjee et al. 2017; Molla et al. 2017), which motivated us to measure the mass of the current BHC Swift J1753.5-0127 using constant normalization method.
We believe that TCAF model fitted N values should be constant for a particular source, since it only depends on physical parameters such as the mass, distance, disk inclination angle, etc. Here, during the entire 2005 outburst of Swift J1753.5-0127, TCAF normalization show a constant signature of , except for the high values () for observations in the initial period of HIMS (Dec.). A significant deviation of the N value could be due to high X-ray emissions from the jet (Soleri et al. 2010), whose effects have not been included in the current version of the TCAF model fits file. Emission from the jet will be discussed in another work (Jana et al. 2017). Mass of the BH () is an important parameter of the TCAF model fits file, which was kept free for all observations. From the spectral fit with the TCAF solution, we find that varies in between to . This is our estimated range of the probable mass the source.
4 Discussions and Concluding Remarks
We study spectral and temporal evolutions of BHC Swift J1753.5-0127 during its 2005 outburst to infer about the accretion flow properties of the source. For temporal study, we use keV RXTE/PCU2 data and for spectral analysis keV data of the same instrument for observation Ids are used. Low frequency QPOs and their nature, such as the shape, frequency, value, rms%, etc. are studied during the entire 2005 outburst of the BHC Swift J1753.5-0127. To study spectral properties, spectra are fitted with the current version (v0.3) of the TCAF solution fits file to extract physical accretion flow parameters, such as, Keplerian disk rate (), sub-Keplerian halo rate () in Eddington rate, shock location () in and shock compression ratio (), and mass of the BH. A comparison is made with the spectral analysis result done with the combined DBB and PL models. A combined analysis of DBB and PL model gives us a rough estimation of disk temperature, photon indices and thermal and non-thermal flux contributions, etc. However, a TCAF model fits provide us with the evolution of the physical parameters relating to the accretion flow itself. Variation of TCAF extracted parameters together with ARR, and observed QPOs help us to determine the spectral state of the BHC. We conclude that this source has been in two harder (HS and HIMS) spectral states during the 2005 outburst, although there are many reports from other authors that during the entire outburst phase source was in low HS (see, Cadolle Bel et al. 2007; Ramadevi & Seetha 2007). No signature of softer (SIMS and SS) spectral states are observed by us. The reason may be the short orbital period ( hrs, Zurita et al. 2007) of the binary system which makes the system very compact. There is a general excess of matter surrounding the whole system. So, there may be a constant supply of sub-Keplerian matter from the wind. This is also the reason that the outburst continued till date in HS with low X-ray intensity. Since during the outburst, this source does not follow canonical hysteresis loop of HID or ‘q’-diagram or ARRID (see, JDCMM16 and references therein), as it never visited softer spectral states, the 2005 outburst of Swift J1753.5-0127 may be termed as a ‘failed’ outburst. Several other BHCs, such as, MAXI J1659-152, MAXI J1836-194, XTE J1118+480, Swift J1357.2-0933 also have shown similar behavior during their outbursts.
Similar to MAXI J1836-194 during its 2011 outburst (JDCMM16), the spectral states are observed in the sequence of: HS (Ris.) HIMS (Ris.) HIMS (Dec.) HS (Dec.) to form a hysteresis loop without softer spectral states. X-ray intensities, QPO frequencies as well as two component accretion rates, PL indices rise during initial five days of the outburst. This phase of the outburst is divided into HS (Ris.) and HIMS (Ris.) with a transition on the 2nd observation (2005 July 4; MJD=53557.24), when a sharp rise in above mentioned parameters from the first observation and a local maximum of the ARR are observed. This property is found to be a common feature on HS HIMS transition days, observed for other transient BHCs, such as H 1743-322, GX 339-4, MAXI J1659-152, MAXI J1836-194 during their outbursts (see, MDC14, DMC15, DMCM15, JDCMM16). After that, X-ray intensity, QPO frequency as well as two component accretion rates are observed to decrease with time (day) till the end of our observations. We define this phase of the outburst into two declining harder spectral states, i.e., HIMS (Dec.) and HS (Dec.). The transition between these spectral states is defined on 2005 August 7 (MJD=53589.49), since after that QPO frequencies, ARRs, PL indices become roughly constant. Our definition of 2005 July 6 (MJD=53557.24) observation as the transition between two hard-intermediate spectral states is clearly justified as on this particular observation maximum values of X-ray intensities, QPO frequencies, Keplerian rates are observed.
QPO frequencies monotonically increase during the rising phase of the outburst, whereas during the declining phase of the outburst a general trend of decreasing nature of QPO frequencies is observed though on a daily basis the decrease is not monotonic. There may be some turbulence present in the flow. According to the TCAF solution, origin of these low frequency QPOs could easily be explained by the resonance oscillation (type-C), weak oscillation (type-B) of the Compton cloud boundary (i.e., shock) or even of the shock-free centrifugal barrier (type-A). So, we suspect that during the non-QPO days (in the declining HS), resonance condition between cooling and infall time of the post-shock region was not satisfied by a large margin (see, Chakrabarti et al., 2015 for more details).
Through TCAF fits, we find that an outburst of a transient BHC is manifested due to a large supply of Keplerian and sub-Keplerian matter close to the inner disk. It is reasonable to assume that this is due to rapid rise in viscosity, perhaps of magnetic origin, owing to enhanced magnetic activity in the companion or convective turbulence at the outer disk (Chakrabarti 2013). Similarly, declining phase of an outburst starts when supply from the companion is turned off, i.e., due to low supply of the high viscous matter. Although the 2005 outburst of Galactic BHC Swift J1753.5-0127 started in the canonical way, it did not move to the burst-off phase for over a decade now. This may be due to the constant supply of matter from the companion star, most probably in the form of wind flow of low angular momentum low viscous sub-Keplerian halo matter. For this reason, sometimes even for high accretion rate in the Keplerian component could not able to cross the comparative limiting value, such that source reaches soft spectral state. The binary system, Swift J1753.5-0127 seems to belong to a type-II class of outbursting sources along with MAXI J1836-194, MAXI J1659-152, XTE J1118+480, etc. of shorter orbital periods, where softer spectral states are also missing during their outbursts.
For comparison of parameters obtained from TCAF fits with those we have obtained from DBB & PL model, we compute the , the inner edge of the so-called truncated standard disk. In TCAF fits, is the shock location which is the inner boundary of a modified standard disk and the outer boundary of the CENBOL. (in ) computed from DBB normalization from the relation , where, is in and is inclination angle, is given in the Table in Appendix -I. The normalization is used for the whole spectrum to match photons emitted at the disk frame and those observed by us. , calculated from DBB norm, is not really a physical parameter. For instance, we see from the Appendix Table that on MJD=53583, is less than 30 km. For a BH mass of , goes inside the event horizon which is not possible. Similar value of was found in seven more observations. In contrast, normalization in TCAF is constant across the spectral states and the size of the Compton cloud could be tens to hundreds of Schwarzschild radii away. In other word, TCAF solves for totally different physical parameters of the system than those from the DBB & PL model.
It is important to note that according to TCAF, all the four parameters (other than the mass and the Normalization which do not change in principle on a daily basis, unless there are addition sources of photons from jets), are essential to fit any general data from a stellar mass black hole having no outflows or magnetic field. In certain period of time, one or more parameters may be constant or nearly constant as special case. For instance in soft/very soft states, only accretion rate of the disk would be required, and two components are not necessary. Fitting with DBB plus PL models often requires an additional ‘reflection’ component to take care of a ‘bump’ like features. In TCAF this is automatically produced when the Keplerian disk is heated up by the fraction of incident Comptonized photons from the CENBOL. Thus we do not require this additional component. It is not unlikely that this reprocessing may not fully account for the reflection component, especially when the disk is partly ionized. However, in our case the Keplerian disk is flanked by a hot advective flow, and thus this problem of partial ionization, especially closer to the black hole, does not arise.
All the four parameters play crucial role in fitting a spectrum even it appears to have simple features. The coupled differential equations given in CT95 produce a spectrum in a self-consistent way. Keplerian disk component, depending on its accretion rate , emits seed photons, a fraction of which are intercepted by the CENBOL. The fraction depends on the optical depth, height (related to the CENBOL temperature) of the CENBOL, and hence depends on all four parameters. The degree of Comptonization of the seed photons will decide the PL slope of the spectrum. If the CENBOL size () gets reduced, it it expected that it should be easier to cool it since the optical depth is much higher in a denser region. However, the reduction in size reduces the path length as well as the number of intercepted seed photons from the Keplerian flow. The Compression ratio () increases the density in CENBOL by this factor, causing more Comptonization. The internal density of the CENBOL depends on the sum of the two rates ( and ) of the accretion flow. As the CENBOL cools down with higher , the possibility of thermal pressure driven jets from CENBOL becomes lesser. Thus in TCAF scenario, Jets and outflows are possible in harder states.
In sort, all the four parameters play their roles in a complex and non-linear way. It is true that when the disk rate is very small, as shown in CT95, the spectral slope is fairly insensitive to the rate, till the time the disk rate becomes comparable to the halo rate to cool CENBOL down. So one would imagine that the would be degenerate. This is not true, since the same rate also is responsible for the multi-color blackbody bump. On the other hand, the resulting spectral slope would be different when the mass changes, since changes the seed photon temperature, without changing the advective flow properties. Hence the information of mass enters through the spectral properties of the softer component. Since TCAF obtains the entire spectrum at a time at the rest frame of the emitter, only one normalization is needed to map the whole spectrum with the observed spectrum. There could be fluctuations due to observational errors, but our model normalization is not expected to change violently for the same object unless there are addition source of X-rays from the base of the outflows and jets. Thus the average normalization obtained from one outburst is expected to fit the data in other outbursts as well.
Using TCAF, we also estimated probable range of the mass of the BHC from our TCAF model fitted spectral study. Mass of the BH is an important parameter of the current TCAF model fits file, which was kept free while fitting spectra. We assumed the model normalization to remain constant (as in MDCMJ16) to estimate mass of the source to be about to . TCAF model normalization (N) does not vary significantly for a particular BHC even when it passes through different outbursts and spectral states, unless there is a precession in the disk which may change the projected surface of the emission area or there are some significant outflow activities which we missed to include in current model fits file. During the entire 2005 outburst of Swift J1753.5-0127, model normalization is found to be constant at with a little variation in between . A large deviation (N ) is observed in five observations during the HIMS (Dec.), and this may be due to high Jet emissions (Soleri et al. 2010) during that period of observations, whose effects were not included in the current version of the TCAF model fits file. In future, we will calculate X-ray flux coming from the jet and study their variations during the outburst using the variation of N during an outburst. This will be published elsewhere. A broad band X-ray study of the source starting from July 2005 till date, using archival data of RXTE (PCA and HEXTE instruments) and SWIFT (XRT and BAT instruments) satellites will also be done and published elsewhere.
A.J. and D.D. acknowledge support from ISRO sponsored RESPOND project fund (ISRO/RES/2/388/2014-15). D.C. and D.D. acknowledge support from DST sponsored Fast-track Young Scientist project fund (SR/FTP/PS-188/2012). SM acknowledge supports FONDECYT post-doctoral fellowship grand (# 3160350).
- Arnaud (1996) Arnaud, K.A., 1996, ASP Conf. Ser., Astronomical Data Analysis Software and Systems V, ed. G.H. Jacoby & J. Barnes, 101, 17
- Cadolle Bel et al. (2007) Cadolle Bel, M., Ribó, M., Rodriguez, J., et al., 2007, ApJ, 659, 549
- Bhattacharjee et al. (2017) Bhattacharjee, A., Banerjee, I., Banerjee, A., et al., 2017, MNRAS, 466, 1372
- Chakrabarti (1990) Chakrabarti, S. K., 1990, “Theory of Transonic Astrophysical Flows”, World Scientific (Singapore)
- Chakrabarti & Titarchuk (1995) Chakrabarti, S.K., & Titarchuk, L.G., 1995, ApJ, 455, 623 (CT95)
- Chakrabarti (1997) Chakrabarti, S.K., 1997, ApJ, 484, 313
- Chakrabarti (1999) Chakrabarti, S.K., 1999, A&A, 351, 185
- Chakrabarti et al. (2008) Chakrabarti, S.K., Debnath, D., & Nandi, A., et al., 2008, A&A, 489, L41
- Chakrabarti (2013) Chakrabarti, S.K., 2013, ASI Conf. Series, Eds. S. Das, A. Nandi and I. Chattopadhyay, 8, 1
- Chakrabarti et al. (2015) Chakrabarti, S.K., Mondal, S., & Debnath, D., 2015, MNRAS, 452, 3451
- Chatterjee et al. (2016) Chatterjee, D., Debnath, D., Chakrabarti, S.K., et al. 2016, ApJ, 827, 88 (CDCMJ16)
- Chatterjee et al. (2017) Chatterjee, D., Debnath, D., Chakrabarti, S.K., et al. 2017 (in preparation)
- Corbel et al. (2005) Corbel, S., Kaaret, P., & Fernder, R.P., et al., 2005, ApJ, 632, 504
- Debnath et al. (2010) Debnath, D., Chakrabarti, S.K., & Nandi, A., 2010, A&A, 520, A98
- Debnath et al. (2013) Debnath, D., Chakrabarti, S.K., & Nandi, A., 2013, AdSpR, 52, 2143
- Debnath, Chakrabarti & Mondal (2014) Debnath, D., Chakrabarti, S.K., & Mondal, S., 2014, MNRAS, 440, L121
- Debnath, Mondal & Chakrabarti (2015a) Debnath, D., Mondal, S., & Chakrabarti, S.K., 2015a, MNRAS, 447, 1984 (DMC15)
- Debnath et al. (2015b) Debnath, D., Molla, A.A., Chakrabarti, S.K., & Mondal, S., 2015b, ApJ, 803, 59 (DMCM15)
- Ebisawa,Titarchuk & Chakrabarti (1996) Ebisawa, K., Titarchuk, L.G., & Chakrabarti, S.K., 1996, PASJ, 48, 59
- Fender et al. (2005) Fender, R. P., Garrington, S., Muxlow, T., 2005, ATel, 558, 1
- Froning et al. (2014) Froning, C.S., Maccarone, T.J., France, K., et al., 2014, ApJ, 780, 48
- Halpern (2005) Halpern, J.P., 2005, ATel, 549, 1
- Jana et al. (2016) Jana, A., Debnath, D., Chakrabarti, S.K., et al., 2016, ApJ, 803, 107 (JDCMM16)
- Jana, Chakrabarti & Debnath (2017) Jana, A., Chakrabarti, S.K., Debnath, D., 2017, ApJ (submitted)
- kong et al. (2017) Kong, A. K. H., 2017, ATel, 10081, 1
- Molla et al. (2016) Molla, A. A., Debnath, D., & Chakrabarti, S. K., et al., 2016, MNRAS, 460, 3163 (MDCMJ16)
- Molla et al. (2017) Molla, A. A., Chakrabarti, S. K., Debnath, D., & Mondal, S., 2017, ApJ, 834, 88
- Molteni, Sponholz & Chakrabarti (1996) Molteni, D., Sponholz, H., & Chakrabarti, S.K., 1996, ApJ, 457, 805
- Mondal, Debnath & Chakrabarti (2014) Mondal, S., Debnath, D., & Chakrabarti, S.K., 2014, ApJ, 786, 4 (MDC14)
- Mondal, Chakrabarti & Debnath (2016) Mondal, S., Chakrabarti, S.K., & Debnath, D., 2016, ApSS, 361, 309
- Nandi et al. (2012) Nandi, A., Debnath, D., Mandal, S., & Chakrabarti, S.K., 2012, A&A, 542, A56
- Novikov & Thorne (1973) Novikov, I., & Thorne, K.S., 1973, in Black Holes, Ed. C. DeWitt & B.S. DeWitt (New York: Gordon & Breach), 343
- (Neustroev et al. 2014) Neustroev, V.V., Veledina, A, Poutanen, J., et al. 2014, MNRAS, 445, 2424
- Palmer et al. (2005) Palmer, D.M., Barthelmey, S.D. & Cummings, J.R., et al. 2005, ATel, 546, 1
- Plotkin et al. (2016) Plotkin, R. M., Miller-Jones, J. C. A., Jonker, P. G., et al., 2016, ATel, 9765, 1
- Qasim et al. (2017) Qasim, A. A., AlMannaei, Aisha., Russell, D. M., et al. 2017, ATel, 10075, 1
- Ramadevi & Seetha (2007) Ramadevi, M.C., Seetha, S., 2007, MNRAS, 378, 182
- Ryu, Chakrabarti & Molteni (1997) Ryu, D., Chakrabarti, S.K., & Molteni, D., 1997, ApJ, 474, 378
- Shakura & Sunyaev (1973) Shakura, N.I., & Sunyaev, R.A., 1973, A&A, 24, 337
- Shaw et al. (2016) Shaw, A.W., Charles, P.A., Casares, J., et al., 2016, MNRAS, 463, 1314
- Shaw et al. (2017) Shaw, A. W., Tomsick, J. A., Bahramian, A., et al., 2016, ATel, 9735, 1
- (Soleri et al. 2010) Soleri, P., Fender, R.P., Tudose, V., et al. 2010, MNRAS, 406, 1471
- Sunyaev & Titarchuk (1980) Sunyaev, R.A., & Titarchuk, L.G., 1980, ApJ, 86, 121
- Sunyaev & Titarchuk (1985) Sunyaev, R.A., & Titarchuk, L.G., 1985, A&A, 143, 374
- Zhang et al. (2007) Zhang, G. B., Qu, J. L., Zhang, S., et al. 2007, ApJ, 659, 1115
- Zurita et al. (2007) Zurita, C., Torres, M.A.P., Durant, M., et al., 2007, ATel, 1130, 1
2.5-25 keV Combined DBB plus PL Model and TCAF Model Fitted Spectral Parameters with QPOs
|30||Y-07-01||53597.30||– –||– –||– –||46.49/41|
|31||Y-08-01||53603.38||– –||– –||– –||53.12/43|
|32||Y-08-02||53605.47||– –||– –||– –||33.77/40|
|33||Y-08-03||53607.25||– –||– –||– –||44.67/43|
|34||Y-09-00||53612.75||– –||– –||– –||56.21/41|
|35||Y-10-00||53619.69||– –||– –||– –||34.92/41|
|36||Y-11-00||53623.76||– –||– –||– –||48.11/40|
|37||Y-12-00||53630.31||– –||– –||– –||43.00/44|
|38||Y-12-01||53634.24||– –||– –||– –||64.99/44|
|39||Y-13-00||53638.65||– –||– –||– –||69.31/44|
|40||Y-13-01||53641.60||– –||– –||– –||66.99/43|
|41||Y-14-00||53646.98||– –||– –||– –||63.53/44|
|42||Y-15-00||53652.67||– –||– –||– –||58.29/44|
|43||Y-15-01||53654.63||– –||– –||– –||49.88/44|
|44||Y-16-01||53662.77||– –||– –||– –||63.57/44|
X=91094-01 and Y=91423-01 are the prefixes of observation Ids. Blank lines mark transitions between different spectral states.
, and values indicate combined DBB and PL model fitted DBB temperatures in keV and PL photon indices respectively.
DBBf, PLf represent combined DBB and PL model fitted fluxes in 2.5-25 keV band for DBB and PL model components respectively in units of .
(in Eddington rate unit), (in Eddington rate unit), (in units of Schwarzschild radius), and are TCAF fitted Keplerian disk, sub-Keplerian halo,
shock location and compression ratio values respectively. Norm. and are TCAF fitted normalization and BH mass values respectively.
Frequencies of the dominating QPO in Hz are mentioned. DOF means degrees of freedom of the spectral model fits.
The values of , and DOF of TCAF model fitted spectra are mentioned in Col. 17. The ratio is the .
Note: average values of 90% confidence error values obtained using ‘err’ task in XSPEC, are placed as superscripts of fitted parameters.
Also note that the accretion rate is sometimes higher than the Eddington rate. Even though the luminosity may be sub-Eddington.
This is because TCAF self-consistently computes the efficiency factor for each data and the observed luminosity .