Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Constraining the mass of the black hole GX 339-4 using spectro-temporal analysis of multiple outbursts

Constraining the mass of the black hole GX 339-4 using spectro-temporal analysis of multiple... We carried out spectro-temporal analysis of the archived data from multiple outbursts spanning over the last two decades from the black hole X-ray binary GX 339-4. In this paper, the mass of the compact object in the X-ray binary system GX 339-4 is constrained based on three indirect methods. The first method uses broadband spectral modelling with a two component flow structure of the accretion around the black hole. The broadband data are obtained from RXTE (Rossi X-ray Timing Explorer) in the range 3.0 to 150.0 keV and from Swift and NuSTAR (Nuclear Spectroscopic Telescope Array) simultaneously in the range 0.5 to 79.0 keV. In the second method, we model the time evolution of Quasi-periodic Oscillation (QPO) frequencies, considering it to be the result of an oscillating shock that radially propagates towards or away from the compact object. The third method is based on scaling a mass dependent parameter from an empirical model of the photon index (Γ) - QPO (ν) correlation. We compare the results at 90 percent confidence from the three methods and summarize the mass estimate of the central object to be in the range 8.28 − 11.89 M . Keywords: Black hole candidate; accretion discs; X-ray outbursts 1. Introduction et al., 2012, and references therein). These states are char- acterized by the variations in spectral and temporal prop- Black Hole (BH) X-ray binaries (XRBs) are very in- erties as the outburst progresses. The Hardness-Intensity teresting objects to study because of their peculiar be- Diagram (HID) (Maccarone and Coppi, 2003; Homan and haviour like sudden increase in intensity leading to an Belloni, 2005; Remillard and McClintock, 2006; Nandi et al., outburst. Such outbursts typically result in the source 2012; Radhika and Nandi, 2014; Nandi et al., 2018) is the undergoing spectral state transitions occupying the low- log-log plot of hardness (ratio of flux in higher energy band hard state (LHS), hard-intermediate state (HIMS), soft- to the flux in lower energy band) versus intensity in the intermediate state (SIMS) and the high-soft state (HSS) total energy band. The HID for an entire outburst usu- during the rising phase and a reverse trend during the de- ally takes the shape of a ‘q’ and hence is often referred cay phase (Homan et al., 2001; Belloni et al., 2005; Nandi to as the ‘q’- diagram. This diagram helps in identifying spectral state transitions of a BH XRB. These state transi- tions give rise to evolution of spectral parameters obtained Corresponding author from various models, which can be studied to get a better Email addresses: hjsreehari@gmail.com (Sreehari H.), understanding of accretion dynamics around a black hole. nirmal.iyer@gmail.com (Nirmal Iyer), radhikad.isac@gmail.com (Radhika D.), anuj@isac.gov.in (Anuj Nandi), It is generally accepted that there is a disc along the sam.cenb@gmail.com (Samir Mandal) Preprint submitted to Advances in Space Research November 13, 2018 arXiv:1811.04341v1 [astro-ph.HE] 11 Nov 2018 equatorial plane of the black hole (Shakura and Sunyaev, 2016). Most of the outbursts of GX 339-4 lasted for a du- 1973) and around the black hole there is a ‘corona’ or ration of a year or more (Sreehari et al., 2018). The evo- ‘Compton cloud’ (Sunyaev and Titarchuk, 1980) which lution of spectral and temporal properties of this source inverse-Comptonize soft radiation from the disc. Com- during the 2002 outburst have been remarkable and paved bining these, several disc-corona models (Svensson and way for having a clear understanding of the HID (Belloni Zdziarski, 1994; Zycki et al., 1995; Done and Kubota, 2006; et al., 2005). Identification and classification of the dif- Ingram and Done, 2011) have been proposed to explain ferent types of QPOs during the RXTE era have already radiations from black hole binary systems. These mod- been done by Belloni et al. (2005); Motta et al. (2011). els have helped to understand in detail the evolution of Corral-Santana et al. (2016) has provided the photo- energy spectral parameters like disc temperature and pho- metric properties determined for this source till date. Though ton index during the outburst. Model with two component many detailed observations have been done on the system accretion flows (Chakrabarti and Titarchuk, 1995; Smith (Samimi et al. (1979), Grebenev et al. (1991)), a constraint et al., 2002; Wu et al., 2002) is also considered often. In on the mass of the central BH is yet to be made. The the two component accretion flow model (Chakrabarti and fact that GX 339−4 lies very much within the galactic Titarchuk, 1995), the BH is considered to be surrounded disc makes it difficult to dynamically estimate its mass, by a post-shock corona (PSC) and apart from the Keple- as even in quiescent state the signature of companion is rian disc there is a sub-Keplerian halo as well. This model hardly detectable. Nevertheless attempts to estimate the is used for broadband spectral modelling to determine ac- mass of the source by Hynes et al. (2003) based on Na- cretion rates and black hole mass (Iyer et al., 2015; Nandi I and Ca-II line emissions resulted in a lower mass limit et al., 2018; Radhika et al., 2018). of 5.8 ± 0.3 M . Further Munoz-Darias et al. (2008) has The temporal properties of BH sources have given clear come up with a tighter bound of minimum 6 M based evidence for the occurrence of low frequency QPOs in the on a ‘stripped giant’ model where in the companion mass range 0.1 - 20 Hz. These are further divided into 3 types is assumed to be ≥ 0.17 M . Later, Shaposhnikov and i.e. A, B and C based on the value of QPO frequency, Titarchuk (2009) (hereafter ST09) has estimated the mass width and amplitude (Casella et al., 2004; Belloni et al., of the black hole in GX 339-4 to be 12.3±1.4 M based on 2005). It has been found that for most of the outbursting a photon index-QPO scaling method. Also, Parker et al. sources the value of QPO frequency increases with time (2016) has estimated the mass of the black hole to be +1.6 during the rising phase of the outburst, gets saturated 9.0 M and a distance of 8.4 ± 0.9 kpc by spectral −1.2 during the SIMS and decreases as the decay phase occurs modelling of XRT and NuSTAR data. But Parker et al. (Belloni et al., 2005; Nandi et al., 2012). C-type QPOs in (2016) have used only the data from the very high state particular go through a monotonic increase in frequency of GX 339-4. Recently, Heida et al. (2017) determined during the rising phase of the outburst. They show the the mass function to be f(M) = 1.91 ± 0.08 M and a reverse trend during the decay phase of the outburst. The mass range of 2.3 − 9.5 M at a binary inclination of 37 C-type QPO frequency is related to the size of the oscil- to 78 . From the above mentioned studies, the mass of lating region (Chakrabarti et al., 2008). From this it is the compact object of GX 339-4 is gauged to be in a very inferred that, during rising phase as the QPO frequency wide range from about 2.3 M to 13.7 M . In this paper, ⊙ ⊙ increases the size of the oscillating region decreases and we attempt to give a better constraint on the mass of the vice versa during the decay phase. The temporal evolu- black hole GX 339-4. tion of C-type QPOs have been modelled by Chakrabarti Using spectral and temporal modelling, mass of the et al. (2008) with black hole mass as a parameter. In- sources IGR J17091−3624 and XTE J1859+226 have been terestingly, a positive correlation has also been observed estimated previously by Iyer et al. (2015) and Nandi et al. for the variation of QPO frequency with the spectral pho- (2018) respectively. In this paper, we explore the spectral ton index (Γ) which keeps increasing in the rising phase and temporal evolution of GX 339−4 during all its out- of the outburst and gets saturated as the source enters its bursts in the RXTE era. We study the variation of QPO SIMS (Titarchuk and Fiorito, 2004). This correlation has frequency with time and the correlation between spectral been empirically modelled by Shaposhnikov and Titarchuk photon index and QPO frequency in order to estimate the (2007) (hereafter ST07) and used to obtain source mass source mass. We also model the X-ray energy spectra from based on a scaling relation. RXTE using the two component flow model. Further as In this paper, we have considered the well known BH RXTE observations are not available after 2012, we anal- source GX 339−4 during its 2002, 2004, 2007 and 2010 yse the broadband spectra from simultaneous Swift-XRT outbursts in the era of RXTE and the 2013 and 2015 and NuSTAR observations for the outbursts in 2013 and 2 3 outburst observed by Swift and NuSTAR . The source 2015. Finally, combining all these results we provide a is also known as 1H J1659−487 (Corral-Santana et al., mass range for the compact object in GX 339-4. The paper is organized as follows: In the next section we discuss the observations and data analysis procedures https://heasarc.gsfc.nasa.gov/docs/xte/XTE.html http://www.swift.ac.uk/ https://www.nustar.caltech.edu/ 2 4 using the software HEAsoft . In § 3, we explain the dif- files. The deadtime correction was done using hxtdead ferent methodologies used to indirectly estimate the mass and the HEXTE response was generated with hxtrsp (see of the compact object in a black hole binary using spectral Radhika et al. (2016) for details). The spectral modelling and temporal data. In § 4, we present the results of ap- was done using XSPEC-12.9.1. For broadband spectral plying these methods on the observational data from the modelling, we used RXTE’s PCA and HEXTE combined source GX 339-4. In § 5 we discuss on results and men- spectrum. We used 1% systematic error for the spectral tion the caveats of the methods used. Finally, we provide analysis with RXTE data. concluding remarks on the results in § 6. For the temporal data, we consider the science event files and made use of the tool seextrct applying a bin time of 8 ms to generate the lightcurves. The Power Density 2. Observation and Analysis Spectrum (PDS) for each data set corresponding to 3.0 - 30.0 keV was obtained from the lightcurves using the We have used High Energy Astrophysics Science Archive powspec tool. These power spectra were modelled with Research Center (HEASARC) archival data to analyse the powerlaws, constants and Lorentzians. We extracted QPO outbursts of GX 339−4 in the period 1997-2015. We anal- parameters for the Lorentzians that had a quality factor ysed the 2002, 2004, 2007 and 2010 outbursts of GX 339-4 ν/FWHM greater than 3, where ν is the Lorentzian cen- using RXTE data and 2013 and 2015 outburst using si- troid and FWHM is the Lorentzian width. The QPO pa- multaneous observations from Swift and NuSTAR. The rameters obtained are consistent with Belloni et al. (2005) 1998-1999 outburst was partially observed by RXTE and and Motta et al. (2011). so the data is insufficient for modelling. The spectral and temporal analysis are performed using HEAsoft v 6.21. 2.2. Swift and NuSTAR 2.1. RXTE As RXTE was decommissioned in 2012, we searched for broadband data from other instruments for the period The RXTE consists of an All Sky Monitor (ASM), Pro- after 2012. So for the 2013 and 2015 outbursts of GX 339- portional Counter Array (PCA) and High-Energy X-ray 4, we used data from Swift-XRT (Evans et al., 2009) and Timing Experiment (HEXTE) instruments. We use data data from NuSTAR observatory to do the simultaneous from PCA and HEXTE for the analysis presented in this broadband spectral analysis. The Swift spectral data were paper. PCA and HEXTE together covers a broad energy obtained using standard procedures in the energy range range extending from 3.0 keV to 150.0 keV. To obtain the 0.5 to 10.0 keV, which is then grouped to a minimum of PCA spectral data, at first we generate the good time 30 counts per bin using the grppha tool. interval (gti) file using maketime, and then apply this gti With the NuSTAR data, we initially ran the nupipeline and filter file information into saextrct. The cmbrightV LE tool to generate the pipeline outputs. The source and background model and the South Atlantic Anomaly (SAA) ′′ background regions each of 30 were chosen from the cleaned history file for PCA are given as inputs to the pcabackest event file using the ds9 tool. Further, we executed the command to estimate background. Further the saextrct nuproducts command to obtain the energy spectra in the command is used to extract the background spectrum for range 3.0 to 79.0 keV, grouped at a minimum of 30 counts PCA. The pcarsp tool is used to generate the response per bin. Combining XRT and NuSTAR for simultaneous matrix. Of the five proportional counter units (PCUs) of observations, we obtain broadband spectra ranging from PCA, we used only PCU2 which was stable and in working 0.5 keV to 79.0 keV. condition over all considered outbursts. The HEXTE instrument has two clusters A and B. We use cluster A data for our analysis. The source spectrum 3. Methodology was generated using saextrct and the background spectrum using hxtback followed by saextrct. For the 2007 outburst, In this paper, we apply three different methods to esti- we used the HEXTE Cluster B data as the cluster A ex- mate the mass of the black holes in X-ray binary systems. perienced a rocking anomaly where the on and off source The first method is based on the two component accretion modulation ceased . For the 2010 outburst as the satel- flow model. The second and third methods involve obser- lite had a rocking problem with cluster B also, the mission vations of QPOs, one based on temporal evolution of QPOs team decided to have cluster B permanently in off-source and the other depends on QPO frequency correlation with position and cluster A in on source position. In this case, the photon index obtained from spectral modelling. the source spectrum was extracted from FS50* files and the corresponding background was generated from FS56* 3.1. Method I: Two Component Accretion Flow files which were further processed with hxtebackest to We consider an accretion disc around a Schwarzschild obtain background spectrum corresponding to the FS50* black hole. The accretion disc consists of two components https://heasarc.nasa.gov/lheasoft/ 6 http://www.swift.ac.uk/user objects/index.php https://heasarc.gsfc.nasa.gov/docs/xte/whatsnew/big.html 3 (Chakrabarti and Titarchuk, 1995), one being the Kep- outbursts after 2012. The results of modelling the energy lerian disc at the equatorial plane and the other a sub- spectra using the two component accretion flow for the dif- Keplerian halo above and below the Keplerian disc. The ferent outbursts of the black hole GX 339-4 are presented in-falling matter necessarily becomes sub-Keplerian (Giri in § 4.1. and Chakrabarti, 2013) near the black hole horizon to sat- isfy the inner boundary conditions. A sub-Keplerian flow 3.2. Method II: QPO evolution with multiple sonic points, becomes supersonic after cross- Low Frequency QPOs (LFQPOs) are often observed ing the outer sonic point and connect the inner sonic point in the power density spectra obtained from fourier trans- through a shock transition (Das et al., 2001) where the forms of X-ray variability. Of the different types of LFQ- centrifugal force becomes comparable to gravity. The Ke- POs, C-type QPOs are usually found to be evolving with plerian disc remains subsonic (Das et al., 2001) through- time. These QPOs are found in the LHS and HIMS of the out and it is truncated at the shock. The region between outburst (Belloni and Motta, 2016). In the rising phase the shock and the black hole horizon is the post shock of the X-ray outburst, the C-type QPO frequencies are corona (hereafter PSC) which is sub-Keplerian in nature. found to increase from around 0.1 Hz to about 20 Hz or At the shock the supersonic flow suddenly becomes sub- so (Casella et al., 2005). During the decaying phase of the sonic and the excess kinetic energy is converted into ther- outburst, QPO frequency evolution exhibits the opposite mal energy. As a result the post-shock region becomes trend starting at a larger frequency and then gradually hotter and denser than the pre-shock region. The soft X- decreasing towards lower frequencies. rays emitted by the X-ray binary system corresponds to The origin of QPOs can be related to oscillating shocks. thermal emission from the Keplerian disc. This can be We consider the Propagating Oscillatory Shock (POS) model modelled with a multi-temperature blackbody like diskbb (Chakrabarti et al., 2008; Nandi et al., 2012; Radhika and model available in XSPEC. Inverse-Comptonization of the Nandi, 2014) that describes the time evolution of C-type soft photons from the optically thick Keplerian disc by the QPO frequency. In POS model, QPO is considered to hot electrons in the post shock corona gives rise to a power- be triggered by an oscillating shock front which propa- law distribution of hard X-rays (Titarchuk and Lyubarskij, gates towards the black hole during the rising phase of an 1995; Chakrabarti and Titarchuk, 1995). In the two com- outburst and moves away from the black hole during the ponent accretion flow model, the contribution from the declining phase of the outburst. In this model the QPO 3/2 soft X-rays are accounted for by the Keplerian-disc ac- frequency is inversely proportional to r (Chakrabarti cretion rate and the contribution from hard X-rays varies and Manickam, 2000), where r is the instantaneous shock with the accretion rate of the sub-Keplerian component. location. In this method, mass of the source is determined by spec- This model is represented by: tral modelling. In calculating the radiation spectrum, the mass of the source decides the soft photon flux from the ν = , (1) qpo Keplerian disc and the number density of electrons in the 2πRr r r − 1 g s s central cloud (Compton corona i.e., PSC) through shock at ut + location. r = r ∓ , (2) s so The two component accretion flow is included in XSPEC r as an additive table model (Iyer et al., 2015). This model where ν is the QPO frequency, R is the shock com- qpo has five parameters namely the black hole mass, halo accre- pression ratio, u is the initial velocity of the shock front tion rate, disc accretion rate, shock location and normal- ρ and a its acceleration. R = is the ratio of density of ization. The accretion rates are in units of Eddington rate, the post-shock region (ρ ) to the density of the pre-shock 2 + shock location is in units of r = 2GM/c (Schwarzschild region (ρ ). The instantaneous shock location, r and the − s radius) where G is the gravitational constant and c is the initial shock location, r are in units of r . The negative so g speed of light and black hole mass is in units of M . The sign in Equation 2 corresponds to the rising phase and the normalization varies proportional to cos(θ)/D where θ is plus sign corresponds to the decaying phase. the angle of inclination of the system in degrees and D is As defined in the previous sub section r depends on the distance to the source in units of 10 kpc. Considering mass of the black hole which is a parameter of this model. a distance of 8.4 kpc and an inclination angle of 37 to So one can obtain mass of the black hole from POS model 78 for GX 339-4, we get a norm value ranging from 0.58 fitting as was done by Iyer et al. (2015); Sreehari et al. to 2.26. For our analysis, we have used the average norm (2018). We apply the POS method on the observed C-type of 1.4. Iyer et al. (2015) and Nandi et al. (2018) has used QPOs during the 2002, 2004, 2007 and 2010 outbursts of this model with mass as a free parameter whereas Debnath the source GX 339-4 to estimate the black hole mass as is et al. (2014) has frozen the mass of the central source. We discussed in § 4.2. follow the procedure adopted by Iyer et al. (2015). We use 3.0 to 150.0 keV RXTE data for the spectral modelling of outbursts before 2012 and 0.5 to 79.0 keV simultane- ous Swift and NuSTAR data for spectral modelling of the 4 3.3. Method III: Photon Index-QPO correlation the photon index of the energy spectra. Using the tempo- ral and spectral parameters obtained from the above men- It has been observed that the QPO frequency and the tioned models fitting, we estimate the mass of the compact photon index of the power-law used to fit the energy spec- object using three indirect methods. trum are well correlated. The correlation is modelled us- All the QPO signatures in power spectra are modelled ing an empirical relation developed by ST07 based on with lorentzians and the lorentzian centroid frequency is Titarchuk and Fiorito (2004). This model is expressed as taken as the QPO frequency. We model the time evolution ν − ν tr of the C-type LFQPO frequencies to obtain the BH source Γ(ν) = A − BD ln exp + 1 , (3) D mass. where A is the index saturation level, B is the initial slope, 4.1. Results from Method I D determines how fast the index goes to saturation and ν tr is the transition frequency. Value of parameter B varies in As mentioned in previous section, we have modelled the proportion to the mass of the compact object. As shown energy spectrum of the source in all the outbursts observed below, for two sources whose correlation plot has similar by RXTE using the two component accretion flow model. saturation the ratio of its mass to B is a constant. Additionally phabs is used to model the interstellar absorp- We model both the target and reference source photon- tion, smedge to model the intrinsic absorption, gaussian to index versus QPO frequency correlation and estimate the model the iron line and a calibration constant is used to parameters A, B, D and ν . Our intention is to estimate tr correct for the offset between PCA and HEXTE spectra. the scaling factor S which is defined as the ratio of mass ν We have used both PCA and HEXTE data in the range of the target black hole to the mass of the reference black 3.0 - 150.0 keV for fitting the model. These broadband hole. In equation 3, if we consider the exponential term to spectra are shown in Figures 1, 2 and 3 for 2002, 2004 and be much larger than unity, we can see that the dimensions 2010 outbursts respectively. The PCA spectra are shown −1 of B corresponds to the unit Hz . Thus, according to in black and the HEXTE spectra are shown in red. We equation 3, B has an inverse relation to QPO frequency. did not use 2007 outburst data for spectral modelling due Also, QPO frequencies are inversely related to mass of the to large offset between PCA and HEXTE spectra. Be- black hole (i.e. M ∝ 1/ν) as is shown in Remillard and sides this, we also fitted the simultaneous Swift-XRT and McClintock (2006) and ST09. This in turn implies that NuSTAR spectra from the 2013 and 2015 outbursts of the mass of the black hole is directly proportional to the value source. This is shown in Figure 4. of parameter B. Thus ν /ν = M /M = S = B /B . We obtain shock location (size of the Compton cloud), r t t r ν t r For instance in ST07, mass of Cyg X-1 is estimated to halo accretion rate, the disc accretion rate and the mass of be 8.7 ± 0.8 from photon index - QPO correlation data the black hole by spectral fitting of each state of the out- of Cyg X-1 spanning over the entire RXTE mission, by bursts using this model. We follow the state classification scaling with respect to the 2005 outburst (decay phase) of done by Motta et al. (2011) except for the 2013 and 2015 source GRO J1655-40 whose dynamically determined mass outbursts. For these cases, we have considered the nature is 6.3± 0.5M (Greene et al., 2001). With GRO J1655-40 ⊙ of the energy spectra based on photon index values and as the reference source and Cyg X-1 as the target source, the contribution of the thermal component to classify the we can write states. The energy spectra corresponding to the four states in M = B . (4) t t 2002 outburst of GX 339-4, modelled with the two com- The results of this scaling method based on photon ponent accretion flow are shown in Figure 1. In the LHS index - QPO frequency correlation for the source GX 339- shown in top-left panel the spectrum extends from 3.0 keV 4 are presented in § 4.3. to 150.0 keV. The shock location corresponding to this spectrum is 274 r which implies a large corona around the 4. Results black hole. The halo accretion rate is 0.83 M showing Edd that the contribution of inverse-Comptonized radiation is The energy spectra from different states of each out- significant. This is because the free-falling electrons in the burst is modelled with the two component accretion flow halo can acquire large kinetic energies as they approach the which has mass as a free parameter. We also modelled black hole. This triggers thermal Comptonization once ki- all outburst energy spectra with pexrav model following netic energy has been converted into thermal energy at the Motta et al. (2011). The reflection parameter of this model shock. The disc accretion rate of 1.70 M indicates the Edd was left free. A multi-colour disc blackbody (diskbb) and presence of a Keplerian disc. As the system transits into a gaussian with centroid around 6.4 keV and width 0.8 the HIMS (top-right panel), the shock location decreases keV was used whenever necessary. The N value of 5 × to 92.4 r , the halo accretion rate changes to 1.12 M g Edd 21 2 10 atoms/cm used for fitting is obtained from NASA’s and the disc accretion rate increases to 4.29 M . This is Edd N value calculator . From this modelling, we obtain evident from the figure where we see comparatively higher photon flux in the 3 to 6 keV range. In the bottom-left https://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 5 panel, we have shown the SIMS spectrum extending from 3 - 100 keV, fitted with the two component flow model. Here MJD 52388.05 0.1 MJD 52410.52 LHS 0.1 HIMS the shock location has reached 28 r and the halo rate has 0.01 0.01 −3 −3 reduced to 0.48 M while the disc rate has become 3.37 10 10 Edd −4 ˙ 10 −4 M . Here also the disc is prominent as is evident from Edd the figure where the photon count rate per unit area is −1 −1 high below 6 keV. And finally in the bottom-right panel, −2 −2 5 10 20 50 100 5 10 20 50 Energy (keV) Energy (keV) we show the fit corresponding to the HSS where the size of the corona has become 25 r , the halo accretion rate de- 0.1 creased to 0.11 M and the disc accretion rate increased Edd MJD 52416.59 0.1 MJD 52687.68 0.01 SIMS HSS to 3.81 M . This HSS spectrum ranges only from 3 - 0.01 Edd −3 −3 −4 60 keV. The mass limits for the BH obtained with this 10 −4 −5 method for the 2002 outburst is 9.04 - 10.61 M . Similarly, we have modelled the other outbursts ob- 0 −1 −2 −2 served with RXTE as well and the results are tabulated in 5 10 20 50 5 10 20 50 Energy (keV) Energy (keV) Table 1. It is observed that, generally the disc accretion rate is low in the hard state and then it increases as the Figure 1: Spectral fitting with two component flow model (Method states evolve reaching a maximum at the soft state while I) for 2002 outburst. PCA spectra is shown in black and HEXTE spectra is shown in red. The hard state spectra extends upto 150.0 the halo accretion rate is maximum at the hard states and keV while the soft state spectra extends only upto 60.0 keV for this reduces as the source enters the soft intermediate and high outburst. soft states. Also the shock location is found to decrease as the states evolve from LHS to HSS suggesting that the corona becomes compact in the soft states. With the Swift-XRT and NuSTAR data of 2013 and 0.1 MJD 53201.32 MJD 53230.95 0.01 LHS HIMS 2015 outburst, we initially applied a phenomenological 0.01 −3 −3 model consisting of diskbb, powerlaw and gaussian. For −4 −4 the 2013 outburst, the powerlaw is dominant with a pho- −5 −5 ton index of 1.43 ± 0.01. The diskbb component accounts 1 1 for only 3 % of the total flux. The source is in LHS dur- −1 −1 −2 −2 5 10 20 50 100 5 10 20 50 100 ing this observation. For the 2015 outburst, the photon Energy (keV) Energy (keV) index is 2.60 ± 0.01 and there is a prominent disc at 0.84 keV implying significant thermal emission. Around 44.8 0.1 MJD 53414.13 MJD 53322.68 0.1 SIMS HSS % of the total flux is contributed by the thermal compo- 0.01 0.01 nent. The source is in its soft intermediate state (SIMS) −3 −3 −4 during this observation. Further, we modelled these two −4 −5 2 2 broadband observations (Figure 4) using two component 1 1 0 0 accretion flow. In the observation corresponding to 2013 −1 −1 −2 −2 5 10 20 50 5 10 outburst the shock location is 137 r , halo rate is 0.54 Energy (keV) Energy (keV) ˙ ˙ M and has a disc rate of 0.13 M . During the obser- Edd Edd Figure 2: Spectral fitting with two component flow model (Method I) vation corresponding to 2015 outburst, the shock location for the 2004 outburst of GX 339-4. The hard state spectrum extends is only 15 r , halo rate is 0.23 M and has a significant g Edd till 150.0 keV while the soft state spectrum extends only up to 15.0 disc rate of 4.92 M . The mass estimates at 90 percent- keV for this outburst. Edd age confidence from all the outbursts analysed are given in Table 1. −2 −1 −1 −2 −1 −1 −2 −1 −1 −2 −1 −1 χ Photons cm s keV χ Photons cm s keV χ Photons cm s keV χ Photons cm s keV −2 −1 −1 −2 −1 −1 −2 −1 −1 −2 −1 −1 χ Photons cm s keV χ Photons cm s keV χ Photons cm s keV χ Photons cm s keV a propagating shock front. So the shock location, shock MJD 55301.79 velocity, shock strength (compression ratio) and mass of 0.1 MJD 55277.47 0.1 HIMS LHS 0.01 the the compact object are the parameters in this model. 0.01 −3 The model is described by the Equations 1 and 2. We −3 −4 10 10 have used the fitted shock locations from two component −5 1 flow model for the first day of QPO of each outburst as −1 −1 the initial guess of shock locations for Method II. −2 −2 5 10 20 50 5 10 20 50 Initially, we have fitted the time evolution of QPO fre- Energy (keV) Energy (keV) quencies with constant velocity, but this gives poor fits. For instance the constant velocity based POS model gave MJD 55315.69 MJD 55354.27 0.1 0.1 SIMS HSS a χ /dof equal to 56/11 for the 2002 outburst QPO evolu- 0.01 0.01 −3 tion. The parameters obtained from this fit are a velocity −3 −4 of 3.7 ± 0.3 m/s and a mass of 8.05 ± 0.45 M . Further, −4 10 ⊙ we introduced an acceleration parameter to the model and 0 0 −1 −2 found that it fits the data properly with a χ close to one. −2 red 5 10 20 50 5 10 Energy (keV) Energy (keV) With the inclusion of the acceleration parameter the fit improved considerably giving a χ /dof equal to 12.44/10. Figure 3: Spectral fitting with two component flow model (Method Here the initial velocity obtained is 8.5 ± 0.65 m/s, the I) for the 2010 outburst of GX 339-4. −06 2 acceleration is −4.46 × 10 m/s and a mass range of 9.97 − 10.81 M . Except for 2004 rising phase data, it 0.1 is observed that the change in velocity is significant and MJD 56528.52 LHS hence the most general model should include shock front 0.01 acceleration as a parameter. The POS model fit results for 2002, 2004, 2007 and 2010 outbursts of GX 339-4 are −3 shown in Figure 5. For the 2002 outburst as QPO frequency increases from −4 0.16 Hz to 5.8 Hz the shock location changes from 278 r to 18.8 r and shock velocity decreases from 8.5 m/s to 6.0 m/s. Similar trend is observed in the rising phases of other outbursts as well. This implies that the shock front slows down as the system transits into SIMS during the −2 rising phase. We also found that during the decay phase 1 10 of 2007 outburst the QPO frequencies evolved from 4.1 Energy (keV) Hz to 0.67 Hz, the corresponding shock location increased 1 MJD 57092.63 from 25.0 r to 82.0 r and the shock velocity increased g g SIMS from 1.52 m/s to 2.00 m/s. So during the decay phase, the 0.1 shock front velocity increases as time progresses. In Table 0.01 2, we present the values of POS model parameters from −3 the QPO frequency evolution fitting. The overall mini- mum to maximum mass range obtained from this method −4 4 is 8.79 − 11.89 M . −2 −4 1 10 Energy (keV) Figure 4: Broadband spectral modelling with Method I for outbursts in 2013 (upper panel) and 2015 (lower panel) showing simultaneous XRT (black) and NuSTAR (red) spectra. The XRT spectral data used, is in the range 0.5 keV to 9.0 keV and the NuSTAR data is from 3.0 keV to 79 keV. 4.2. Results from Method II In the POS model, C-Type QPOs which are generally found in LHS and HIMS states of the rising phase of out- bursts are considered to be formed by the oscillations of −2 −1 −1 −2 −1 −1 χ Photons cm s keV χ Photons cm s keV −2 −1 −1 −2 −1 −1 (data−model)/error Photons cm s keV (data−model)/error Photons cm s keV −2 −1 −1 −2 −1 −1 χ Photons cm s keV χ Photons cm s keV Table 1: Fitted parameters from Two component Flow model. The XRT-NuSTAR spectra are denoted by †. All other observations are from RXTE. ˙ ˙ MJD Outburst State r (r ) m ˙ (M ) m ˙ (M ) Mass (M ) χ /dof s g h Edd d Edd ⊙ 52388.05 2002 LHS 274 ± 15 0.83 ± 0.01 1.70 ± 0.03 10.49 ± 0.12 111.42/70 52410.52 2002 HIMS 92.4 ± 4 1.12 ± 0.02 4.29 ± 0.09 09.59 ± 0.12 88.79/62 52416.59 2002 SIMS 28 ± 1 0.48 ± 0.002 3.37 ± 0.07 09.55 ± 0.07 113.81/62 52687.68 2002 HSS 25 ± 1 0.11 ± 0.001 3.81 ± 0.4 09.49 ± 0.45 63.94/60 53201.32 2004 LHS 136 ± 12 0.38 ± 0.007 0.30 ± 0.02 09.74 ± 0.32 97.30/90 53230.95 2004 HIMS 97.5 ± 1.5 0.27 ± 0.003 2.27 ± 0.04 08.89 ± 0.10 85.89/82 53322.68 2004 SIMS 10.6 ± 0.21 0.12 ± 0.007 6.50 ± 0.23 10.19 ± 0.63 66.27/49 53414.13 2004 HSS 17 ± 0.5 0.09 ± 0.001 6.69 ± 0.25 10.05 ± 0.86 13.95/10 55277.47 2010 LHS 453 ± 34 0.30 ± 0.005 2.9 ± 0.12 10.16 ± 0.16 125.72/117 55301.79 2010 HIMS 43 ± 1.3 0.29 ± 0.003 5.9 ± 0.17 08.97 ± 0.12 167.97/115 55315.69 2010 SIMS 23 ± 0.8 0.18 ± 0.002 9.28 ± 0.38 09.87 ± 0.27 125.15/79 55354.27 2010 HSS 5.7 ± 0.04 0.04 ± 0.001 7.10 ± 0.13 08.60 ± 0.30 14.46/13 56528.52 † 2013 LHS 137 ± 3 0.54 ± 0.003 0.13 ± 0.001 09.89 ± 0.06 1521/1296 57092.63 † 2015 SIMS 15 ± 0.1 0.23 ± 0.001 4.92 ± 0.02 10.17 ± 0.04 2580/1503 Table 2: POS model fitted parameters (Method II) 2 2 Outburst rs (r ) v (m/s) acceleration (m/s ) Mass (M ) χ /dof 0 g 0 ⊙ −06 2002 278 ± 4 8.50 ± 0.65 −4.46 × 10 9.97 − 10.81 12.44/10 −08 2004 143 ± 2 1.21 ± 0.09 +3.55 × 10 9.97 − 10.79 17.28/12 −06 2007 299 ± 7 12.32 ± 0.09 −9.00 × 10 10.04 − 10.49 15.87/20 −07 2007decay 25 ± 0.6 1.52 ± 0.06 +5.00 × 10 9.09 − 10.99 17.92/15 −06 2010 487 ± 22 11.72 ± 2.15 −4.79 × 10 8.79 − 11.89 14.42/13 tra. We fitted this correlation using empirical formulae de- veloped in ST07. We could not apply the Method III for 2002 outburst data of GX 339-4 as a matching reference POS modelling - GX 339-4 outburst is unavailable. So, we used the data from out- bursts in 2004, 2007 and 2010 where QPOs are observed for implementing this method. XTE J1550-564 with a dy- +1.02 namically estimated mass of 10.56 M (Orosz et al. −0.88 (2002)) and GRO J1655-40 with a dynamically estimated mass of 6.3±0.5 M (Greene et al., 2001) are the reference sources used. For the 2004 outburst as the QPO frequency increases from 0.30 Hz to 8.05 Hz the photon index (Γ) value changes from 1.45 to 2.21. The corresponding reference source used is the 2005 outburst of GRO J1655-40. During the 2007 -2 outburst Γ of GX 339-4 evolves from 1.74 to 2.76 as the 0 5 10 15 20 25 30 35 QPO evolves from its lowest to highest value. The ref- Days since QPO started erence source used is the 1998a outburst of XTE J 1550- 564. Finally the range of Γ for the 2010 outburst is 1.50 Figure 5: Evolution of QPO frequencies modelled with Method II to 2.56 which matches the reference outburst 1998b of the (POS model) for the rising phase of the outbursts of GX 339-4. The lower panel shows the residuals obtained after POS model fitting. source XTE J1550-564. ST09 has classified 1998 outburst of XTE J1550-564 into 1998a and 1998b based on photon index evolution during the outburst. 4.3. Results from Method III The fits corresponding to Method III are shown in Fig- As was described earlier the photon index parameter ure 6. For photon index - QPO frequency correlation, we obtained by modelling the energy spectra are correlated have used data from same energy ranges for both reference with the QPO frequencies of the corresponding power spec- χ QPO (Hz) and target. We have used 3 - 30 keV data for the 2004 and 2010 outbursts for both reference and target sources and 3 - 100 keV data for the case of 2007 outburst. This is 2.2 because the photon index values are significantly different over the two energy ranges and a reference outburst with 2.0 matching photon index evolution for the 2007 outburst of GX 339-4 was obtained only from the broadband spectral 1.8 analysis of both the target and reference sources. The re- sults from Method III presented in this paper are based 1.6 on the relations given in ST07 (see Equations 3 and 4). GX 2004 For the reference source GRO J1655-40 during its 2005 1.4 GRO 2005 outburst rising phase the parameter B value obtained from 0.05 1.0 10.0 +0.01 QPO (Hz) the correlation fitting is 0.15 . And for 2004 outburst −0.003 +0.02 of GX 339-4, we have B = 0.22 . Substituting these −0.01 values in Equation 4, we get the mass of GX 339-4 to be 2.8 in the range 8.50−10.85 M . Similarly we have computed the mass of the target source using data from different out- bursts. In Table 3, we have provided the parameter values 2.4 obtained by fitting the Gamma (Γ) - QPO (ν) correlation to different outbursts of GX 339-4 and the reference out- bursts used. Thus by using the constant M/B relation in 2.0 Equation 4, method III gives an over all mass estimate for GX 339-4 in the range 8.28 to 11.81 M . The mass 1.6 GX 2007 estimates for the different outbursts are given in Table 3. XTE 1998a 0.05 1.0 10.0 QPO (Hz) 2.6 2.2 1.8 GX 2010 1.4 XTE 1998b 0.05 1.0 10.0 QPO (Hz) Figure 6: Fitting of the QPO-photon index correlation based on Method III (ST07). The reference source is shown in blue and the target source in green. The reference sources are chosen such that the photon index saturation levels are similar to that of the target source. 4.4. Consolidated Results Most of the outbursting black hole sources show spo- radic jet activities (Fender et al., 2004, 2009; Radhika et al., 2016) and X-ray variabilities (Sreehari et al., 2018; Radhika et al., 2018) in the intermediate states. So we intend to get a tight constraint on the mass of GX 339- 4 using only LHS data. For this we use the results from method I, which is the most robust method among the three that we considered. As broadband LHS data is not available for 2015 outburst, we have used the only available SIMS data for this case. We have converted the outputs of Photon Index Photon Index Photon Index Table 3: Model fitted parameter values (Method III) Source - Outburst A B ν (Hz) D Mass (M ) tr ⊙ +0.34 +0.02 GX 339 - 2004 2.28 0.22 3.0 2.72 ± 0.08 8.50 − 10.85 −0.14 −0.01 +0.08 +0.035 GX 339 - 2007 2.73 0.266 3.8 0.42 ± 0.24 8.58 − 11.81 −0.06 −0.004 +0.09 +0.02 +0.69 GX 339 - 2010 2.51 0.40 2.0 1.81 8.28 − 10.49 −0.06 −0.03 −0.12 Reference - Outburst +0.01 GRO J1655 - 2005 2.29 ± 0.03 0.15 6.0 1.0 6.3 ± 0.5 −0.003 +0.031 +1.02 XTE J1550 - 1998a 2.79 0.285 4.2 1.0 10.56 −0.027 −0.88 +0.03 +0.25 +1.02 XTE J1550 - 1998b 2.35 ± 0.03 0.46 2.0 0.68 10.56 −0.02 −0.09 −0.88 the steppar command of XSPEC into a probability density function (PDF) following the steps in Nandi et al. (2018). The same method has been employed in Iyer et al. (2015) and Radhika et al. (2018) to combine different mass es- timates. Figure 7 shows the PDF corresponding to each outburst and we have also indicated the union and in- tersection of these PDFs. Each PDF is of unit area. The union PDF suggests a mass range of 9.72 - 10.70 M where as the intersection of PDFs gives a tighter constraint of 10.10 - 10.35 M . We consider the results from intersec- tion of PDFs, as per the naive Bayes rule for combining independent probabilities (Iyer et al., 2015). 2002 LHS 2004 LHS 2010 LHS 2013 LHS 4 2015 SIMS Intersection Union Figure 8: Forest plot of the mass estimates from different methods at 90% confidence. The size of the square represents the significance of each result. The maximum and minimum values from all three methods are used to plot the final range. 8.0 8.5 9.0 9.5 10.0 10.5 11.0 11.5 12.0 Mass (M ) As was mentioned earlier, we exclude the 2007 outburst Figure 7: The PDFs corresponding to mass estimates from spectral results from Method I due to large offset between PCA and modelling of each outburst is plotted here. Except for 2015 outburst, HEXTE spectra. For method II, we have used all rising we have used data from LHS state to obtain the PDFs. The yellow phase QPO evolution data available from RXTE. We have shaded portion is the intersection of individual PDFs and the cyan shaded portion is the union of individual PDFs. See text for more considered the decay phase data from 2007 as well for POS details and mass limits. modelling, but did not include the decay phase data from 2010 outburst (Nandi et al., 2012) as the QPO evolution is We also include together the results from all the three not smooth. And finally for method III, we have included methods in the Forest plot in Figure 8. The method and three outbursts for which we could find matching reference year are shown in the first column, the corresponding mass sources. For each individual method, we have considered ranges at 90% confidence are shown in the second column the minimum and maximum mass values as the range of and a plot of the mass range is shown in third column. the mass parameter. Thus the summary measure is given The size of the squares in the plot shows how significant as the minimum to maximum mass values considering all the results are. three methods together. So the mass of the black hole in GX 339-4 at a 90% confidence range is 8.28 − 11.89 M . https://cran.r-project.org/web/packages/forestplot/vignettes/forestplot.html PDF It should be noted that the PDFs corresponding to that the basic POS model with a constant velocity does Methods II and III will be broad (see Iyer et al. (2015)). not fit the temporal evolution of QPO perfectly. Hence This is also evident from the mass estimates in Tables we introduced an acceleration parameter that accounts for 2 and 3 which covers the range obtained from intersec- the change in velocity of the shock front. We found that tion of PDFs from Method I. Hence including PDFs from the shock front decelerates as the system approaches the Methods II and III would not improve the estimate ob- soft intermediate state during the rising phase. We also tained from the joint PDF (intersection) of Method I. The noted that the shock front accelerates as it recedes during mass estimate from the likelihood analysis (both intersec- the decay phase. During the rising phase of the outburst, tion and union of PDFs) lies within the extreme values particularly in LHS and HIMS, the increase of frequency obtained from the Forest plot. of QPO with time demands a forward movement of shock front. The location of shock in the accretion flow is deter- mined by the Rankine-Hugoniot conditions which states a 5. Discussion balance of both energy and pressure of the flow before and Method I is based on broadband spectral modelling us- after the shock. The key physical parameters are mass ing the two component accretion flow. For this we used accretion rate and viscosity of the flow. Accretion rate RXTE data in the range 3.0 - 150.0 keV and combined controls the rate of cooling and hence affect the total pres- XRT - NuSTAR data in the range 0.5 - 79.0 keV. From the sure whereas viscosity decides the transfer of angular mo- data of multiple outbursts, the mass of GX 339-4 estimated mentum. Of course in a magnetized accretion flow, the using this method to be 8.30 - 10.91 M . Though the magnetic pressure also needs to be taken into account. methods we employed are already used by other authors, An increase of accretion rate increases the cooling and this is the first time broadband spectral data has been shock front moves inward. The same trends are expected modelled to constrain the mass of the system. The results with the increases of viscosity as well as strength of mag- of spectral modelling from Method I makes it clear that netic field. In general, the dependency of these effects the error bars on mass parameter from RXTE data analy- (rate of cooling, pressure etc.) with the flow parameters sis (outbursts in 2002 to 2010) is significantly larger than are non-linear and hence an accelerated/decelerated move- the error bars on parameters obtained with XRT - NuS- ment of the shock front may result along with the evolu- TAR data. It can be seen that, the mass range obtained tion of the outburst. A constant speed of shock front is using Method I from 2013 outburst data is 9.82 − 9.95 M , a special case due to internal adjustment of the underly- where as the same method applied to 2015 outburst data ing physical processes. Sukov´a and Janiuk (2015) details gave 10.13 − 10.21 M . This is because these two data are the shock formation and numerical analysis of oscillating obtained when the source was in different states of the out- shock propagation in the scenario of super massive as well burst. The data systematics will be more in the case of the as stellar mass black holes. SIMS as compared to LHS due to presence of other activ- In Method II, there is uncertainty in the initial shock ity like relativistic jet ejections (Fender et al., 2004, 2009). location as one is not sure of the exact time at which the Also, an analytical spectral modelling for optically thin QPO has originated. Based on the observations, we con- (LHS state) or optically thick (HSS state) conditions pro- sider the time at which the lowest frequency is detected duce the required results, but in SIMS the central cloud’s (first instance of detection) as the zero time for modelling (i.e., Compton cloud) Thompson optical depth approaches QPO frequency evolution in the rising phase. However, unity where the analytically calculated spectrum may dif- this is unlikely to be accurate and hence the mass esti- fer from the expected one (say using monte-carlo simu- mate might also have some uncertainty. Moreover, in this lation) and this may also introduce some amount of un- method also, we have considered that the compression ra- certainty in the estimation. Besides this one should also tio is a constant value equal to 3 for all the outbursts. As consider that the RXTE data starts only from 3 keV on- this need not be the actual scenario, some uncertainty in wards, whereas the XRT - NuSTAR data starts from 0.5 the estimated mass is possible. keV onwards and hence we get better constraints from The third method used is the scaling based on pho- 2013 and 2015 outburst analysis. In this method, we have ton index - QPO frequency correlation. The over all mass fixed a constant value of 3 for the compression ratio which range obtained from this method is 8.28 − 11.81 M . The is within the theoretical limits of 1 to 4 for hydrodynamic disadvantage of this method is the requirement of an ap- shocks (Mandal and Chakrabarti, 2005). As compression propriate reference source whose photon index saturation ratio may vary as the shock location advances and the levels match with that of the target source. As even the black hole system transits over different states, it can in- same source does not have similar index saturation levels troduce some uncertainty on the estimated source mass. on different outbursts, it is not possible to find a perfect In the POS model (Method II), the QPO is considered reference source for this method. Moreover, even if the to originate from the oscillations of a propagating shock saturation levels match, the evolution of the photon index front. This model has mass of the black hole as a free pa- is significantly different for target and reference sources. rameter. Fitting the QPO frequency evolution using POS It must be noted that ST09 has obtained a higher mass model gives us a mass range of 8.79−11.89 M . We found estimate with a ΔM/M = 22%. ST09 modified equation 11 3 by introducing a parameter β to get, maximum mass value obtained from this method over all the outbursts considered is 8.30 M and 10.91 M . ⊙ ⊙     Mass determination using the method II, which is based 1 − on the POS model also gave us consistent results. Here we tr     Γ(ν) = A − BD ln exp + 1 . (5)     studied the temporal evolution of the fundamental QPO frequency as a function of the propagating shock location. As per this model, the C-type LFQPO frequency varies −3/2 After getting a good fit using Equation 5 for the refer- approximately as r . The mass range of the black hole ence source correlation, it is multiplied by a scaling factor obtained using this method is 8.79 − 11.89 M . along the frequency axis in order to get a new scaled curve The third method is based on the fact that the photon that matches the photon index-QPO correlation data of index of the energy spectrum has a correlation with the the target source. This scaling factor is further used to es- QPO frequency. As the QPO frequency increases the pho- timate the mass of the target source from the mass of the ton index also increases initially and then saturates. This reference source as mass and QPO frequency are inversely correlation has been empirically modelled by ST07 with an related. But this geometrical scaling of the fitted curve expression having four variables A, B, D and ν . We fit tr of the reference pattern does not give a perfect match on the correlation using Equation 3 for a reference source as the target patterns considered. Moreover, the scaled curve well as for the target source. The reference source should is quite different from the actual fitted curve of the tar- have similar photon index saturation levels to that of the get pattern. As a result, here we presented the results of target source during its outburst. Once the fit parameters method III based on ST07 rather than ST09. for both reference and target source are obtained, we use Using the three methods considered in this paper, we the scaling relation given in Equation 4 to estimate the have estimated the mass of GX 339-4 to be 8.28 − 11.89 M mass of the target source. Using this method, the mass with a ΔM/M of 35%. The ΔM/M for Parker et al. estimate of GX 339-4 at 90 % confidence gives the result (2016) including both negative and positive error is 31%. 8.28 − 11.81 M . But Parker et al. (2016) has used only data from the very The forest plot in Figure 8 shows all the mass estimates high state, while we have considered multiple outbursts we obtained. We quote the mass estimate based on the and all different states of the outbursts. Moreover, based three indirect methods to be in the range 8.28 − 11.89 M on a likelihood analysis we have obtained the most likely by considering the overall maximum and minimum val- mass of 10.10 - 10.35 M and a broader range of 9.72 - ues. Previous estimate by Parker et al. (2016) used only 10.70 M obtained from union of PDFs. Here ΔM/M is the data from very high state of the source and the re- 9.5%. sults partially overlaps with our estimate. Similarly the Corral-Santana et al. (2016); Tetarenko et al. (2016) mass determined by Heida et al. (2017) and Shaposhnikov provides the mass estimates of more than twenty compact and Titarchuk (2009) also partially overlaps with our re- objects in binary systems. Of this only GRS 1915+105, 4U sults. Through a likelihood analysis we also infer that the 1956+350, XTE J1550-564 and 4U 0540-697 have masses most likely mass for the source is 10.10 - 10.35 M . A larger than GX 339-4. So in the distribution of dynamical wider limit of 9.72 - 10.70 M is obtained from the union and candidate BH XRBs, GX 339-4 is a comparatively of PDFs. In method I, we have used data from multi- higher mass black hole. ple outbursts and also analysed different states of the BH X-ray binary system in each outburst. In method II, we have used the time evolution of C-type LFQPO frequen- 6. Conclusion cies in the LHS and HIMS state of the black hole binary. In this paper, we aimed at determining the mass of the In method III, we have used the LFQPOs in each out- Galactic black hole GX 339-4 using spectral and tempo- burst appearing in all the states. Thus using a complete ral techniques. We used three different methods to do the spectro-temporal analysis of the outbursts of GX 339-4 we mass estimation. Together these methods give a reason- have constrained the mass of the black hole in this system. ably narrow range of mass for the source. Among the BHBs whose mass has been estimated so far, The first method we used for mass determination is the GX 339-4 is a comparatively high mass black hole. two component accretion flow model wherein mass of the compact object is a free parameter. The other parame- Acknowledgements ters in this model are disc and halo accretion rates and shock location. We have tabulated the parameter values The authors thank both reviewers for their comments obtained by fitting this model to different outburst data and suggestions that helped to improve the quality of the in Table 1. It is usually observed that, the halo rate is manuscript. This research has made use of data obtained dominant in the hard states and the disc rate is dominant through the HEASARC web service, provided by the NASA/ in the soft states. Also the shock location which indicates Goddard Space Flight Center and also made use of data the size of the PSC is larger in the hard states and it supplied by the UK Swift Science Data Centre at the Uni- decreases as the spectra gets softer. The minimum and versity of Leicester. We thank GD, SAG; DD, PDMSA 12 and Director, URSC for encouragement and support to Giri, K., Chakrabarti, S.K., 2013. Hydrodynamic simulation of two-component advective flows around black holes. Mon. Not. carry out this research. R. Astron. Soc. 430, 2836–2843. doi:10.1093/mnras/stt087, arXiv:1212.6493. Grebenev, S., Sunyaev, R., Pavlinsky, M., Dekhanov, I., 1991. De- References tection of quasi-periodic oscillations of x-radiation from the black hole candidate gx339-4. Pisma v Astronomicheskii Zhurnal 17, Belloni, T., Homan, J., Casella, P., Van Der Klis, M., Nespoli, E., 985–990. Lewin, W., Miller, J., M´endez, M., 2005. The evolution of the Greene, J., Bailyn, C.D., Orosz, J.A., 2001. Optical and infrared timing properties of the black-hole transient gx 339–4 during its photometry of the microquasar gro j1655–40 in quiescence. The 2002/2003 outburst. Astronomy & Astrophysics 440, 207–222. Astrophysical Journal 554, 1290. Belloni, T.M., Motta, S.E., 2016. Transient Black Hole Bi- Heida, M., Jonker, P.G., Torres, M.A.P., Chiavassa, A., 2017. naries, in: Bambi, C. (Ed.), Astrophysics of Black Holes: The Mass Function of GX 339-4 from Spectroscopic Obser- From Fundamental Aspects to Latest Developments, p. 61. vations of Its Donor Star. Astrophysical Journal 846, 132. doi:10.1007/978-3-662-52859-4_2, arXiv:1603.07872. doi:10.3847/1538-4357/aa85df, arXiv:1708.04667. Casella, P., Belloni, T., Homan, J., Stella, L., 2004. A study of the Homan, J., Belloni, T., 2005. The Evolution of Black low-frequency quasi-periodic oscillations in the X-ray light curves Hole States. Astrophysics and Space Science 300, 107–117. of the black hole candidate XTE J1859+226. Astronomy and doi:10.1007/s10509-005-1197-4, arXiv:astro-ph/0412597. Astrophysics 426, 587–600. doi:10.1051/0004-6361:20041231, Homan, J., Wijnands, R., van der Klis, M., Belloni, T., van Paradijs, arXiv:astro-ph/0407262. J., Klein-Wolt, M., Fender, R., M´endez, M., 2001. Correlated Casella, P., Belloni, T., Stella, L., 2005. The ABC of Low-Frequency X-Ray Spectral and Timing Behavior of the Black Hole Candi- Quasi-periodic Oscillations in Black Hole Candidates: Analo- date XTE J1550-564: A New Interpretation of Black Hole States. gies with Z Sources. The Astrophysical Journal 629, 403–407. Astrophys. J. Suppl. Ser. 132, 377–402. doi:10.1086/318954, doi:10.1086/431174, arXiv:astro-ph/0504318. arXiv:astro-ph/0001163. Chakrabarti, S., Titarchuk, L.G., 1995. Spectral Properties Hynes, R.I., Steeghs, D., Casares, J., Charles, P., O’Brien, K., 2003. of Accretion Disks around Galactic and Extragalactic Black Dynamical evidence for a black hole in gx 339–4. The Astrophys- Holes. Astrophysical Journal 455, 623. doi:10.1086/176610, ical Journal Letters 583, L95. arXiv:astro-ph/9510005. Ingram, A., Done, C., 2011. A physical model for the con- Chakrabarti, S.K., Debnath, D., Nandi, A., Pal, P.S., 2008. Evo- tinuum variability and quasi-periodic oscillation in accreting lution of the quasi-periodic oscillation frequency in GRO J1655- black holes. Mon. Not. R. Astron. Soc. 415, 2323–2335. 40 - Implications for accretion disk dynamics. Astronomy and doi:10.1111/j.1365-2966.2011.18860.x, arXiv:1101.2336. Astrophysics 489, L41–L44. doi:10.1051/0004-6361:200810136, Iyer, N., Nandi, A., Mandal, S., 2015. Determination of the Mass of arXiv:0809.0876. IGR J17091-3624 from “Spectro-temporal” Variations during the Chakrabarti, S.K., Manickam, S.G., 2000. Correlation among Onset Phase of the 2011 Outburst. Astrophysical Journal 807, Quasi-Periodic Oscillation Frequencies and Quiescent-State Du- 108. doi:10.1088/0004-637X/807/1/108, arXiv:1505.02529. ration in Black Hole Candidate GRS 1915+105. Astro- Maccarone, T.J., Coppi, P.S., 2003. Hysteresis in the light physical Journal Letters 531, L41–L44. doi:10.1086/312512, curves of soft X-ray transients. Mon. Not. R. Astron. arXiv:astro-ph/9910012. Soc. 338, 189–196. doi:10.1046/j.1365-8711.2003.06040.x, Corral-Santana, J.M., Casares, J., Mun˜oz-Darias, T., Bauer, F.E., arXiv:astro-ph/0209116. Mart´ınez-Pais, I.G., Russell, D.M., 2016. BlackCAT: A cata- Mandal, S., Chakrabarti, S.K., 2005. Accretion shock sig- logue of stellar-mass black holes in X-ray transients. Astronomy natures in the spectrum of two-temperature advective flows and Astrophysics 587, A61. doi:10.1051/0004-6361/201527130, around black holes. Astronomy and Astrophysics 434, 839–848. arXiv:1510.08869. doi:10.1051/0004-6361:20041235. Das, S., Chattopadhyay, I., Chakrabarti, S.K., 2001. Stand- Motta, S., Munoz-Darias, T., Casella, P., Belloni, T., Homan, J., ing Shocks around Black Holes: An Analytical Study. 2011. Low-frequency oscillations in black holes: a spectral-timing Astrophysical Journal 557, 983–989. doi:10.1086/321692, approach to the case of gx 339-4. Mon. Not. R. Astron. Soc. 418, arXiv:astro-ph/0107046. 2292–2307. Debnath, D., Chakrabarti, S.K., Mondal, S., 2014. Implementation Munoz-Darias, T., Casares, J., Mart´ınez-Pais, I., 2008. On the of two-component advective flow solution in XSPEC. Mon. Not. masses and evolutionary status of the black hole binary gx 339-4: R. Astron. Soc. 440, L121–L125. doi:10.1093/mnrasl/slu024, a twin system of xte j1550-564? Monthly Notices of the Royal arXiv:1402.0989. Astronomical Society 385, 2205–2209. Done, C., Kubota, A., 2006. Disc-corona energetics in the very Nandi, A., Debnath, D., Mandal, S., Chakrabarti, S.K., 2012. Ac- high state of Galactic black holes. Mon. Not. R. Astron. cretion flow dynamics during the evolution of timing and spectral Soc. 371, 1216–1230. doi:10.1111/j.1365-2966.2006.10737.x, properties of gx 339-4 during its 2010–11 outburst. Astronomy & arXiv:astro-ph/0511030. Astrophysics 542, A56. Evans, P.A., Beardmore, A.P., Page, K.L., Osborne, J.P., O’Brien, Nandi, A., Mandal, S., Sreehari, H., Radhika, D., Das, S., P.T., Willingale, R., Starling, R.L.C., Burrows, D.N., Godet, Chattopadhyay, I., Iyer, N., Agrawal, V.K., Aktar, R., O., Vetere, L., Racusin, J., Goad, M.R., Wiersema, K., An- 2018. Accretion flow dynamics during 1999 outburst of XTE gelini, L., Capalbi, M., Chincarini, G., Gehrels, N., Kennea, J.A., J1859+226-modeling of broadband spectra and constraining Margutti, R., Morris, D.C., Mountford, C.J., Pagani, C., Perri, the source mass. Astrophysics & Space Science 363, 90. M., Romano, P., Tanvir, N., 2009. Methods and results of an doi:10.1007/s10509-018-3314-1, arXiv:1803.08638. automatic analysis of a complete sample of Swift-XRT obser- Orosz, J.A., Groot, P.J., van der Klis, M., McClintock, J.E., Garcia, vations of GRBs. Mon. Not. R. Astron. Soc. 397, 1177–1201. M.R., Zhao, P., Jain, R.K., Bailyn, C.D., Remillard, R.A., 2002. doi:10.1111/j.1365-2966.2009.14913.x, arXiv:0812.3662. Dynamical evidence for a black hole in the microquasar xte j1550– Fender, R.P., Belloni, T.M., Gallo, E., 2004. Towards a unified 564. The Astrophysical Journal 568, 845. model for black hole X-ray binary jets. Mon. Not. R. Astron. Parker, M.L., Tomsick, J.A., Kennea, J.A., Miller, J.M., Harrison, Soc. 355, 1105–1118. doi:10.1111/j.1365-2966.2004.08384.x, F.A., Barret, D., Boggs, S.E., Christensen, F.E., Craig, W.W., arXiv:astro-ph/0409360. Fabian, A.C., Fu¨rst, F., Grinberg, V., Hailey, C.J., Romano, Fender, R.P., Homan, J., Belloni, T.M., 2009. Jets from black hole X- P., Stern, D., Walton, D.J., Zhang, W.W., 2016. NuSTAR and ray binaries: testing, refining and extending empirical models for Swift Observations of the Very High State in GX 339-4: Weigh- the coupling to X-rays. Mon. Not. R. Astron. Soc. 396, 1370–1382. ing the Black Hole with X-Rays. Astrophys. J. Lett. 821, L6. doi:10.1111/j.1365-2966.2009.14841.x, arXiv:0903.5166. 13 doi:10.3847/2041-8205/821/1/L6, arXiv:1603.03777. servational aspects of outbursting black hole sources: Evolution of Radhika, D., Nandi, A., 2014. spectro-temporalcharacteristics and spectro-temporal features and X-ray variability. Journal of Astro- disk-jet connection of the outbursting black hole source xte j1859+ physics and Astronomy 39, 5. doi:10.1007/s12036-018-9510-0, 226. Advances in Space Research 54, 1678–1697. arXiv:1802.05163. Radhika, D., Nandi, A., Agrawal, V.K., Seetha, S., 2016. ‘Spectro- Sukova´, P., Janiuk, A., 2015. Oscillating shocks in the low temporal’ variabilities and possible physical mechanism for angular momentum flows as a source of variability of accret- jet ejections. Mon. Not. R. Astron. Soc. 460, 4403–4416. ing black holes. Mon. Not. R. Astron. Soc. 447, 1565–1579. doi:10.1093/mnras/stw1239, arXiv:1605.08351. doi:10.1093/mnras/stu2544, arXiv:1411.7836. Radhika, D., Sreehari, H., Nandi, A., Iyer, N., Mandal, S., Sunyaev, R.A., Titarchuk, L.G., 1980. Comptonization of X-rays in 2018. Broad-band spectral evolution and temporal variabil- plasma clouds - Typical radiation spectra. Astronomy and Astro- ity of IGR J17091-3624 during its 2016 outburst: SWIFT physics 86, 121–138. and NuSTAR results. Astrophysics & Space Science 363, 189. Svensson, R., Zdziarski, A.A., 1994. Black hole accretion doi:10.1007/s10509-018-3411-1, arXiv:1808.05556. disks with coronae. The Astrophysical Journal 436, 599–606. Remillard, R.A., McClintock, J.E., 2006. X-Ray Proper- doi:10.1086/174934. ties of Black-Hole Binaries. Annu. Rev. Astron. Astro- Tetarenko, B.E., Sivakoff, G.R., Heinke, C.O., Gladstone, J.C., 2016. phys. 44, 49–92. doi:10.1146/annurev.astro.44.051905.092532, WATCHDOG: A Comprehensive All-sky Database of Galactic arXiv:astro-ph/0606352. Black Hole X-ray Binaries. The Astrophysical Journal Supplement Samimi, J., Share, G.H., Wood, K., Yentis, D., Meekins, J., Evans, 222, 15. doi:10.3847/0067-0049/222/2/15, arXiv:1512.00778. W.D., Shulman, S., Byram, E.T., Chubb, T.A., Friedman, H., Titarchuk, L., Fiorito, R., 2004. Spectral Index and Quasi- 1979. GX339-4 - A new black hole candidate. Nature 278, 434– Periodic Oscillation Frequency Correlation in Black Hole Sources: 436. doi:10.1038/278434a0. Observational Evidence of Two Phases and Phase Transition Shakura, N.I., Sunyaev, R.A., 1973. Black holes in binary systems. in Black Holes. The Astrophysical Journal 612, 988–999. Observational appearance. Astronomy and Astrophysics 24, 337– doi:10.1086/422573, arXiv:astro-ph/0405360. 355. Titarchuk, L., Lyubarskij, Y., 1995. Power-Law Spectra as a Result Shaposhnikov, N., Titarchuk, L., 2007. Determination of black hole of Comptonization of the Soft Radiation in a Plasma Cloud. The mass in cygnus x-1 by scaling of spectral index-qpo frequency Astrophysical Journal 450, 876. doi:10.1086/176191. correlation. The Astrophysical Journal 663, 445, ST07. Wu, K., Soria, R., Campbell-Wilson, D., Hannikainen, D., Harmon, Shaposhnikov, N., Titarchuk, L., 2009. Determination of black hole B.A., Hunstead, R., Johnston, H., McCollough, M., McIntyre, V., masses in galactic black hole binaries using scaling of spectral and 2002. The 1998 Outburst of XTE J1550-564: A Model Based on variability characteristics. The Astrophysical Journal 699, 453, Multiwavelength Observations. The Astrophysical Journal 565, ST09. 1161–1168. doi:10.1086/324328, arXiv:astro-ph/0109222. Smith, D.M., Heindl, W.A., Swank, J.H., 2002. Two Different Zycki, P.T., Collin-Souffrin, S., Czerny, B., 1995. Accretion Discs Long-Term Behaviors in Black Hole Candidates: Evidence for with Accreting Coronae in Active Galactic Nuclei - Part One - Two Accretion Flows? The Astrophysical Journal 569, 362–380. Solutions in Hydrostatic Equilibrium. Mon. Not. R. Astron. Soc. doi:10.1086/339167, arXiv:astro-ph/0103304. 277, 70. doi:10.1093/mnras/277.1.70, arXiv:astro-ph/9505045. Sreehari, H., Nandi, A., Radhika, D., Iyer, N., Mandal, S., 2018. Ob- http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Astrophysics arXiv (Cornell University)

Constraining the mass of the black hole GX 339-4 using spectro-temporal analysis of multiple outbursts

Astrophysics , Volume 2018 (1811) – Nov 11, 2018

Loading next page...
 
/lp/arxiv-cornell-university/constraining-the-mass-of-the-black-hole-gx-339-4-using-spectro-SkSKYfX60Y
ISSN
0273-1177
eISSN
ARCH-3330
DOI
10.1016/j.asr.2018.10.042
Publisher site
See Article on Publisher Site

Abstract

We carried out spectro-temporal analysis of the archived data from multiple outbursts spanning over the last two decades from the black hole X-ray binary GX 339-4. In this paper, the mass of the compact object in the X-ray binary system GX 339-4 is constrained based on three indirect methods. The first method uses broadband spectral modelling with a two component flow structure of the accretion around the black hole. The broadband data are obtained from RXTE (Rossi X-ray Timing Explorer) in the range 3.0 to 150.0 keV and from Swift and NuSTAR (Nuclear Spectroscopic Telescope Array) simultaneously in the range 0.5 to 79.0 keV. In the second method, we model the time evolution of Quasi-periodic Oscillation (QPO) frequencies, considering it to be the result of an oscillating shock that radially propagates towards or away from the compact object. The third method is based on scaling a mass dependent parameter from an empirical model of the photon index (Γ) - QPO (ν) correlation. We compare the results at 90 percent confidence from the three methods and summarize the mass estimate of the central object to be in the range 8.28 − 11.89 M . Keywords: Black hole candidate; accretion discs; X-ray outbursts 1. Introduction et al., 2012, and references therein). These states are char- acterized by the variations in spectral and temporal prop- Black Hole (BH) X-ray binaries (XRBs) are very in- erties as the outburst progresses. The Hardness-Intensity teresting objects to study because of their peculiar be- Diagram (HID) (Maccarone and Coppi, 2003; Homan and haviour like sudden increase in intensity leading to an Belloni, 2005; Remillard and McClintock, 2006; Nandi et al., outburst. Such outbursts typically result in the source 2012; Radhika and Nandi, 2014; Nandi et al., 2018) is the undergoing spectral state transitions occupying the low- log-log plot of hardness (ratio of flux in higher energy band hard state (LHS), hard-intermediate state (HIMS), soft- to the flux in lower energy band) versus intensity in the intermediate state (SIMS) and the high-soft state (HSS) total energy band. The HID for an entire outburst usu- during the rising phase and a reverse trend during the de- ally takes the shape of a ‘q’ and hence is often referred cay phase (Homan et al., 2001; Belloni et al., 2005; Nandi to as the ‘q’- diagram. This diagram helps in identifying spectral state transitions of a BH XRB. These state transi- tions give rise to evolution of spectral parameters obtained Corresponding author from various models, which can be studied to get a better Email addresses: hjsreehari@gmail.com (Sreehari H.), understanding of accretion dynamics around a black hole. nirmal.iyer@gmail.com (Nirmal Iyer), radhikad.isac@gmail.com (Radhika D.), anuj@isac.gov.in (Anuj Nandi), It is generally accepted that there is a disc along the sam.cenb@gmail.com (Samir Mandal) Preprint submitted to Advances in Space Research November 13, 2018 arXiv:1811.04341v1 [astro-ph.HE] 11 Nov 2018 equatorial plane of the black hole (Shakura and Sunyaev, 2016). Most of the outbursts of GX 339-4 lasted for a du- 1973) and around the black hole there is a ‘corona’ or ration of a year or more (Sreehari et al., 2018). The evo- ‘Compton cloud’ (Sunyaev and Titarchuk, 1980) which lution of spectral and temporal properties of this source inverse-Comptonize soft radiation from the disc. Com- during the 2002 outburst have been remarkable and paved bining these, several disc-corona models (Svensson and way for having a clear understanding of the HID (Belloni Zdziarski, 1994; Zycki et al., 1995; Done and Kubota, 2006; et al., 2005). Identification and classification of the dif- Ingram and Done, 2011) have been proposed to explain ferent types of QPOs during the RXTE era have already radiations from black hole binary systems. These mod- been done by Belloni et al. (2005); Motta et al. (2011). els have helped to understand in detail the evolution of Corral-Santana et al. (2016) has provided the photo- energy spectral parameters like disc temperature and pho- metric properties determined for this source till date. Though ton index during the outburst. Model with two component many detailed observations have been done on the system accretion flows (Chakrabarti and Titarchuk, 1995; Smith (Samimi et al. (1979), Grebenev et al. (1991)), a constraint et al., 2002; Wu et al., 2002) is also considered often. In on the mass of the central BH is yet to be made. The the two component accretion flow model (Chakrabarti and fact that GX 339−4 lies very much within the galactic Titarchuk, 1995), the BH is considered to be surrounded disc makes it difficult to dynamically estimate its mass, by a post-shock corona (PSC) and apart from the Keple- as even in quiescent state the signature of companion is rian disc there is a sub-Keplerian halo as well. This model hardly detectable. Nevertheless attempts to estimate the is used for broadband spectral modelling to determine ac- mass of the source by Hynes et al. (2003) based on Na- cretion rates and black hole mass (Iyer et al., 2015; Nandi I and Ca-II line emissions resulted in a lower mass limit et al., 2018; Radhika et al., 2018). of 5.8 ± 0.3 M . Further Munoz-Darias et al. (2008) has The temporal properties of BH sources have given clear come up with a tighter bound of minimum 6 M based evidence for the occurrence of low frequency QPOs in the on a ‘stripped giant’ model where in the companion mass range 0.1 - 20 Hz. These are further divided into 3 types is assumed to be ≥ 0.17 M . Later, Shaposhnikov and i.e. A, B and C based on the value of QPO frequency, Titarchuk (2009) (hereafter ST09) has estimated the mass width and amplitude (Casella et al., 2004; Belloni et al., of the black hole in GX 339-4 to be 12.3±1.4 M based on 2005). It has been found that for most of the outbursting a photon index-QPO scaling method. Also, Parker et al. sources the value of QPO frequency increases with time (2016) has estimated the mass of the black hole to be +1.6 during the rising phase of the outburst, gets saturated 9.0 M and a distance of 8.4 ± 0.9 kpc by spectral −1.2 during the SIMS and decreases as the decay phase occurs modelling of XRT and NuSTAR data. But Parker et al. (Belloni et al., 2005; Nandi et al., 2012). C-type QPOs in (2016) have used only the data from the very high state particular go through a monotonic increase in frequency of GX 339-4. Recently, Heida et al. (2017) determined during the rising phase of the outburst. They show the the mass function to be f(M) = 1.91 ± 0.08 M and a reverse trend during the decay phase of the outburst. The mass range of 2.3 − 9.5 M at a binary inclination of 37 C-type QPO frequency is related to the size of the oscil- to 78 . From the above mentioned studies, the mass of lating region (Chakrabarti et al., 2008). From this it is the compact object of GX 339-4 is gauged to be in a very inferred that, during rising phase as the QPO frequency wide range from about 2.3 M to 13.7 M . In this paper, ⊙ ⊙ increases the size of the oscillating region decreases and we attempt to give a better constraint on the mass of the vice versa during the decay phase. The temporal evolu- black hole GX 339-4. tion of C-type QPOs have been modelled by Chakrabarti Using spectral and temporal modelling, mass of the et al. (2008) with black hole mass as a parameter. In- sources IGR J17091−3624 and XTE J1859+226 have been terestingly, a positive correlation has also been observed estimated previously by Iyer et al. (2015) and Nandi et al. for the variation of QPO frequency with the spectral pho- (2018) respectively. In this paper, we explore the spectral ton index (Γ) which keeps increasing in the rising phase and temporal evolution of GX 339−4 during all its out- of the outburst and gets saturated as the source enters its bursts in the RXTE era. We study the variation of QPO SIMS (Titarchuk and Fiorito, 2004). This correlation has frequency with time and the correlation between spectral been empirically modelled by Shaposhnikov and Titarchuk photon index and QPO frequency in order to estimate the (2007) (hereafter ST07) and used to obtain source mass source mass. We also model the X-ray energy spectra from based on a scaling relation. RXTE using the two component flow model. Further as In this paper, we have considered the well known BH RXTE observations are not available after 2012, we anal- source GX 339−4 during its 2002, 2004, 2007 and 2010 yse the broadband spectra from simultaneous Swift-XRT outbursts in the era of RXTE and the 2013 and 2015 and NuSTAR observations for the outbursts in 2013 and 2 3 outburst observed by Swift and NuSTAR . The source 2015. Finally, combining all these results we provide a is also known as 1H J1659−487 (Corral-Santana et al., mass range for the compact object in GX 339-4. The paper is organized as follows: In the next section we discuss the observations and data analysis procedures https://heasarc.gsfc.nasa.gov/docs/xte/XTE.html http://www.swift.ac.uk/ https://www.nustar.caltech.edu/ 2 4 using the software HEAsoft . In § 3, we explain the dif- files. The deadtime correction was done using hxtdead ferent methodologies used to indirectly estimate the mass and the HEXTE response was generated with hxtrsp (see of the compact object in a black hole binary using spectral Radhika et al. (2016) for details). The spectral modelling and temporal data. In § 4, we present the results of ap- was done using XSPEC-12.9.1. For broadband spectral plying these methods on the observational data from the modelling, we used RXTE’s PCA and HEXTE combined source GX 339-4. In § 5 we discuss on results and men- spectrum. We used 1% systematic error for the spectral tion the caveats of the methods used. Finally, we provide analysis with RXTE data. concluding remarks on the results in § 6. For the temporal data, we consider the science event files and made use of the tool seextrct applying a bin time of 8 ms to generate the lightcurves. The Power Density 2. Observation and Analysis Spectrum (PDS) for each data set corresponding to 3.0 - 30.0 keV was obtained from the lightcurves using the We have used High Energy Astrophysics Science Archive powspec tool. These power spectra were modelled with Research Center (HEASARC) archival data to analyse the powerlaws, constants and Lorentzians. We extracted QPO outbursts of GX 339−4 in the period 1997-2015. We anal- parameters for the Lorentzians that had a quality factor ysed the 2002, 2004, 2007 and 2010 outbursts of GX 339-4 ν/FWHM greater than 3, where ν is the Lorentzian cen- using RXTE data and 2013 and 2015 outburst using si- troid and FWHM is the Lorentzian width. The QPO pa- multaneous observations from Swift and NuSTAR. The rameters obtained are consistent with Belloni et al. (2005) 1998-1999 outburst was partially observed by RXTE and and Motta et al. (2011). so the data is insufficient for modelling. The spectral and temporal analysis are performed using HEAsoft v 6.21. 2.2. Swift and NuSTAR 2.1. RXTE As RXTE was decommissioned in 2012, we searched for broadband data from other instruments for the period The RXTE consists of an All Sky Monitor (ASM), Pro- after 2012. So for the 2013 and 2015 outbursts of GX 339- portional Counter Array (PCA) and High-Energy X-ray 4, we used data from Swift-XRT (Evans et al., 2009) and Timing Experiment (HEXTE) instruments. We use data data from NuSTAR observatory to do the simultaneous from PCA and HEXTE for the analysis presented in this broadband spectral analysis. The Swift spectral data were paper. PCA and HEXTE together covers a broad energy obtained using standard procedures in the energy range range extending from 3.0 keV to 150.0 keV. To obtain the 0.5 to 10.0 keV, which is then grouped to a minimum of PCA spectral data, at first we generate the good time 30 counts per bin using the grppha tool. interval (gti) file using maketime, and then apply this gti With the NuSTAR data, we initially ran the nupipeline and filter file information into saextrct. The cmbrightV LE tool to generate the pipeline outputs. The source and background model and the South Atlantic Anomaly (SAA) ′′ background regions each of 30 were chosen from the cleaned history file for PCA are given as inputs to the pcabackest event file using the ds9 tool. Further, we executed the command to estimate background. Further the saextrct nuproducts command to obtain the energy spectra in the command is used to extract the background spectrum for range 3.0 to 79.0 keV, grouped at a minimum of 30 counts PCA. The pcarsp tool is used to generate the response per bin. Combining XRT and NuSTAR for simultaneous matrix. Of the five proportional counter units (PCUs) of observations, we obtain broadband spectra ranging from PCA, we used only PCU2 which was stable and in working 0.5 keV to 79.0 keV. condition over all considered outbursts. The HEXTE instrument has two clusters A and B. We use cluster A data for our analysis. The source spectrum 3. Methodology was generated using saextrct and the background spectrum using hxtback followed by saextrct. For the 2007 outburst, In this paper, we apply three different methods to esti- we used the HEXTE Cluster B data as the cluster A ex- mate the mass of the black holes in X-ray binary systems. perienced a rocking anomaly where the on and off source The first method is based on the two component accretion modulation ceased . For the 2010 outburst as the satel- flow model. The second and third methods involve obser- lite had a rocking problem with cluster B also, the mission vations of QPOs, one based on temporal evolution of QPOs team decided to have cluster B permanently in off-source and the other depends on QPO frequency correlation with position and cluster A in on source position. In this case, the photon index obtained from spectral modelling. the source spectrum was extracted from FS50* files and the corresponding background was generated from FS56* 3.1. Method I: Two Component Accretion Flow files which were further processed with hxtebackest to We consider an accretion disc around a Schwarzschild obtain background spectrum corresponding to the FS50* black hole. The accretion disc consists of two components https://heasarc.nasa.gov/lheasoft/ 6 http://www.swift.ac.uk/user objects/index.php https://heasarc.gsfc.nasa.gov/docs/xte/whatsnew/big.html 3 (Chakrabarti and Titarchuk, 1995), one being the Kep- outbursts after 2012. The results of modelling the energy lerian disc at the equatorial plane and the other a sub- spectra using the two component accretion flow for the dif- Keplerian halo above and below the Keplerian disc. The ferent outbursts of the black hole GX 339-4 are presented in-falling matter necessarily becomes sub-Keplerian (Giri in § 4.1. and Chakrabarti, 2013) near the black hole horizon to sat- isfy the inner boundary conditions. A sub-Keplerian flow 3.2. Method II: QPO evolution with multiple sonic points, becomes supersonic after cross- Low Frequency QPOs (LFQPOs) are often observed ing the outer sonic point and connect the inner sonic point in the power density spectra obtained from fourier trans- through a shock transition (Das et al., 2001) where the forms of X-ray variability. Of the different types of LFQ- centrifugal force becomes comparable to gravity. The Ke- POs, C-type QPOs are usually found to be evolving with plerian disc remains subsonic (Das et al., 2001) through- time. These QPOs are found in the LHS and HIMS of the out and it is truncated at the shock. The region between outburst (Belloni and Motta, 2016). In the rising phase the shock and the black hole horizon is the post shock of the X-ray outburst, the C-type QPO frequencies are corona (hereafter PSC) which is sub-Keplerian in nature. found to increase from around 0.1 Hz to about 20 Hz or At the shock the supersonic flow suddenly becomes sub- so (Casella et al., 2005). During the decaying phase of the sonic and the excess kinetic energy is converted into ther- outburst, QPO frequency evolution exhibits the opposite mal energy. As a result the post-shock region becomes trend starting at a larger frequency and then gradually hotter and denser than the pre-shock region. The soft X- decreasing towards lower frequencies. rays emitted by the X-ray binary system corresponds to The origin of QPOs can be related to oscillating shocks. thermal emission from the Keplerian disc. This can be We consider the Propagating Oscillatory Shock (POS) model modelled with a multi-temperature blackbody like diskbb (Chakrabarti et al., 2008; Nandi et al., 2012; Radhika and model available in XSPEC. Inverse-Comptonization of the Nandi, 2014) that describes the time evolution of C-type soft photons from the optically thick Keplerian disc by the QPO frequency. In POS model, QPO is considered to hot electrons in the post shock corona gives rise to a power- be triggered by an oscillating shock front which propa- law distribution of hard X-rays (Titarchuk and Lyubarskij, gates towards the black hole during the rising phase of an 1995; Chakrabarti and Titarchuk, 1995). In the two com- outburst and moves away from the black hole during the ponent accretion flow model, the contribution from the declining phase of the outburst. In this model the QPO 3/2 soft X-rays are accounted for by the Keplerian-disc ac- frequency is inversely proportional to r (Chakrabarti cretion rate and the contribution from hard X-rays varies and Manickam, 2000), where r is the instantaneous shock with the accretion rate of the sub-Keplerian component. location. In this method, mass of the source is determined by spec- This model is represented by: tral modelling. In calculating the radiation spectrum, the mass of the source decides the soft photon flux from the ν = , (1) qpo Keplerian disc and the number density of electrons in the 2πRr r r − 1 g s s central cloud (Compton corona i.e., PSC) through shock at ut + location. r = r ∓ , (2) s so The two component accretion flow is included in XSPEC r as an additive table model (Iyer et al., 2015). This model where ν is the QPO frequency, R is the shock com- qpo has five parameters namely the black hole mass, halo accre- pression ratio, u is the initial velocity of the shock front tion rate, disc accretion rate, shock location and normal- ρ and a its acceleration. R = is the ratio of density of ization. The accretion rates are in units of Eddington rate, the post-shock region (ρ ) to the density of the pre-shock 2 + shock location is in units of r = 2GM/c (Schwarzschild region (ρ ). The instantaneous shock location, r and the − s radius) where G is the gravitational constant and c is the initial shock location, r are in units of r . The negative so g speed of light and black hole mass is in units of M . The sign in Equation 2 corresponds to the rising phase and the normalization varies proportional to cos(θ)/D where θ is plus sign corresponds to the decaying phase. the angle of inclination of the system in degrees and D is As defined in the previous sub section r depends on the distance to the source in units of 10 kpc. Considering mass of the black hole which is a parameter of this model. a distance of 8.4 kpc and an inclination angle of 37 to So one can obtain mass of the black hole from POS model 78 for GX 339-4, we get a norm value ranging from 0.58 fitting as was done by Iyer et al. (2015); Sreehari et al. to 2.26. For our analysis, we have used the average norm (2018). We apply the POS method on the observed C-type of 1.4. Iyer et al. (2015) and Nandi et al. (2018) has used QPOs during the 2002, 2004, 2007 and 2010 outbursts of this model with mass as a free parameter whereas Debnath the source GX 339-4 to estimate the black hole mass as is et al. (2014) has frozen the mass of the central source. We discussed in § 4.2. follow the procedure adopted by Iyer et al. (2015). We use 3.0 to 150.0 keV RXTE data for the spectral modelling of outbursts before 2012 and 0.5 to 79.0 keV simultane- ous Swift and NuSTAR data for spectral modelling of the 4 3.3. Method III: Photon Index-QPO correlation the photon index of the energy spectra. Using the tempo- ral and spectral parameters obtained from the above men- It has been observed that the QPO frequency and the tioned models fitting, we estimate the mass of the compact photon index of the power-law used to fit the energy spec- object using three indirect methods. trum are well correlated. The correlation is modelled us- All the QPO signatures in power spectra are modelled ing an empirical relation developed by ST07 based on with lorentzians and the lorentzian centroid frequency is Titarchuk and Fiorito (2004). This model is expressed as taken as the QPO frequency. We model the time evolution ν − ν tr of the C-type LFQPO frequencies to obtain the BH source Γ(ν) = A − BD ln exp + 1 , (3) D mass. where A is the index saturation level, B is the initial slope, 4.1. Results from Method I D determines how fast the index goes to saturation and ν tr is the transition frequency. Value of parameter B varies in As mentioned in previous section, we have modelled the proportion to the mass of the compact object. As shown energy spectrum of the source in all the outbursts observed below, for two sources whose correlation plot has similar by RXTE using the two component accretion flow model. saturation the ratio of its mass to B is a constant. Additionally phabs is used to model the interstellar absorp- We model both the target and reference source photon- tion, smedge to model the intrinsic absorption, gaussian to index versus QPO frequency correlation and estimate the model the iron line and a calibration constant is used to parameters A, B, D and ν . Our intention is to estimate tr correct for the offset between PCA and HEXTE spectra. the scaling factor S which is defined as the ratio of mass ν We have used both PCA and HEXTE data in the range of the target black hole to the mass of the reference black 3.0 - 150.0 keV for fitting the model. These broadband hole. In equation 3, if we consider the exponential term to spectra are shown in Figures 1, 2 and 3 for 2002, 2004 and be much larger than unity, we can see that the dimensions 2010 outbursts respectively. The PCA spectra are shown −1 of B corresponds to the unit Hz . Thus, according to in black and the HEXTE spectra are shown in red. We equation 3, B has an inverse relation to QPO frequency. did not use 2007 outburst data for spectral modelling due Also, QPO frequencies are inversely related to mass of the to large offset between PCA and HEXTE spectra. Be- black hole (i.e. M ∝ 1/ν) as is shown in Remillard and sides this, we also fitted the simultaneous Swift-XRT and McClintock (2006) and ST09. This in turn implies that NuSTAR spectra from the 2013 and 2015 outbursts of the mass of the black hole is directly proportional to the value source. This is shown in Figure 4. of parameter B. Thus ν /ν = M /M = S = B /B . We obtain shock location (size of the Compton cloud), r t t r ν t r For instance in ST07, mass of Cyg X-1 is estimated to halo accretion rate, the disc accretion rate and the mass of be 8.7 ± 0.8 from photon index - QPO correlation data the black hole by spectral fitting of each state of the out- of Cyg X-1 spanning over the entire RXTE mission, by bursts using this model. We follow the state classification scaling with respect to the 2005 outburst (decay phase) of done by Motta et al. (2011) except for the 2013 and 2015 source GRO J1655-40 whose dynamically determined mass outbursts. For these cases, we have considered the nature is 6.3± 0.5M (Greene et al., 2001). With GRO J1655-40 ⊙ of the energy spectra based on photon index values and as the reference source and Cyg X-1 as the target source, the contribution of the thermal component to classify the we can write states. The energy spectra corresponding to the four states in M = B . (4) t t 2002 outburst of GX 339-4, modelled with the two com- The results of this scaling method based on photon ponent accretion flow are shown in Figure 1. In the LHS index - QPO frequency correlation for the source GX 339- shown in top-left panel the spectrum extends from 3.0 keV 4 are presented in § 4.3. to 150.0 keV. The shock location corresponding to this spectrum is 274 r which implies a large corona around the 4. Results black hole. The halo accretion rate is 0.83 M showing Edd that the contribution of inverse-Comptonized radiation is The energy spectra from different states of each out- significant. This is because the free-falling electrons in the burst is modelled with the two component accretion flow halo can acquire large kinetic energies as they approach the which has mass as a free parameter. We also modelled black hole. This triggers thermal Comptonization once ki- all outburst energy spectra with pexrav model following netic energy has been converted into thermal energy at the Motta et al. (2011). The reflection parameter of this model shock. The disc accretion rate of 1.70 M indicates the Edd was left free. A multi-colour disc blackbody (diskbb) and presence of a Keplerian disc. As the system transits into a gaussian with centroid around 6.4 keV and width 0.8 the HIMS (top-right panel), the shock location decreases keV was used whenever necessary. The N value of 5 × to 92.4 r , the halo accretion rate changes to 1.12 M g Edd 21 2 10 atoms/cm used for fitting is obtained from NASA’s and the disc accretion rate increases to 4.29 M . This is Edd N value calculator . From this modelling, we obtain evident from the figure where we see comparatively higher photon flux in the 3 to 6 keV range. In the bottom-left https://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 5 panel, we have shown the SIMS spectrum extending from 3 - 100 keV, fitted with the two component flow model. Here MJD 52388.05 0.1 MJD 52410.52 LHS 0.1 HIMS the shock location has reached 28 r and the halo rate has 0.01 0.01 −3 −3 reduced to 0.48 M while the disc rate has become 3.37 10 10 Edd −4 ˙ 10 −4 M . Here also the disc is prominent as is evident from Edd the figure where the photon count rate per unit area is −1 −1 high below 6 keV. And finally in the bottom-right panel, −2 −2 5 10 20 50 100 5 10 20 50 Energy (keV) Energy (keV) we show the fit corresponding to the HSS where the size of the corona has become 25 r , the halo accretion rate de- 0.1 creased to 0.11 M and the disc accretion rate increased Edd MJD 52416.59 0.1 MJD 52687.68 0.01 SIMS HSS to 3.81 M . This HSS spectrum ranges only from 3 - 0.01 Edd −3 −3 −4 60 keV. The mass limits for the BH obtained with this 10 −4 −5 method for the 2002 outburst is 9.04 - 10.61 M . Similarly, we have modelled the other outbursts ob- 0 −1 −2 −2 served with RXTE as well and the results are tabulated in 5 10 20 50 5 10 20 50 Energy (keV) Energy (keV) Table 1. It is observed that, generally the disc accretion rate is low in the hard state and then it increases as the Figure 1: Spectral fitting with two component flow model (Method states evolve reaching a maximum at the soft state while I) for 2002 outburst. PCA spectra is shown in black and HEXTE spectra is shown in red. The hard state spectra extends upto 150.0 the halo accretion rate is maximum at the hard states and keV while the soft state spectra extends only upto 60.0 keV for this reduces as the source enters the soft intermediate and high outburst. soft states. Also the shock location is found to decrease as the states evolve from LHS to HSS suggesting that the corona becomes compact in the soft states. With the Swift-XRT and NuSTAR data of 2013 and 0.1 MJD 53201.32 MJD 53230.95 0.01 LHS HIMS 2015 outburst, we initially applied a phenomenological 0.01 −3 −3 model consisting of diskbb, powerlaw and gaussian. For −4 −4 the 2013 outburst, the powerlaw is dominant with a pho- −5 −5 ton index of 1.43 ± 0.01. The diskbb component accounts 1 1 for only 3 % of the total flux. The source is in LHS dur- −1 −1 −2 −2 5 10 20 50 100 5 10 20 50 100 ing this observation. For the 2015 outburst, the photon Energy (keV) Energy (keV) index is 2.60 ± 0.01 and there is a prominent disc at 0.84 keV implying significant thermal emission. Around 44.8 0.1 MJD 53414.13 MJD 53322.68 0.1 SIMS HSS % of the total flux is contributed by the thermal compo- 0.01 0.01 nent. The source is in its soft intermediate state (SIMS) −3 −3 −4 during this observation. Further, we modelled these two −4 −5 2 2 broadband observations (Figure 4) using two component 1 1 0 0 accretion flow. In the observation corresponding to 2013 −1 −1 −2 −2 5 10 20 50 5 10 outburst the shock location is 137 r , halo rate is 0.54 Energy (keV) Energy (keV) ˙ ˙ M and has a disc rate of 0.13 M . During the obser- Edd Edd Figure 2: Spectral fitting with two component flow model (Method I) vation corresponding to 2015 outburst, the shock location for the 2004 outburst of GX 339-4. The hard state spectrum extends is only 15 r , halo rate is 0.23 M and has a significant g Edd till 150.0 keV while the soft state spectrum extends only up to 15.0 disc rate of 4.92 M . The mass estimates at 90 percent- keV for this outburst. Edd age confidence from all the outbursts analysed are given in Table 1. −2 −1 −1 −2 −1 −1 −2 −1 −1 −2 −1 −1 χ Photons cm s keV χ Photons cm s keV χ Photons cm s keV χ Photons cm s keV −2 −1 −1 −2 −1 −1 −2 −1 −1 −2 −1 −1 χ Photons cm s keV χ Photons cm s keV χ Photons cm s keV χ Photons cm s keV a propagating shock front. So the shock location, shock MJD 55301.79 velocity, shock strength (compression ratio) and mass of 0.1 MJD 55277.47 0.1 HIMS LHS 0.01 the the compact object are the parameters in this model. 0.01 −3 The model is described by the Equations 1 and 2. We −3 −4 10 10 have used the fitted shock locations from two component −5 1 flow model for the first day of QPO of each outburst as −1 −1 the initial guess of shock locations for Method II. −2 −2 5 10 20 50 5 10 20 50 Initially, we have fitted the time evolution of QPO fre- Energy (keV) Energy (keV) quencies with constant velocity, but this gives poor fits. For instance the constant velocity based POS model gave MJD 55315.69 MJD 55354.27 0.1 0.1 SIMS HSS a χ /dof equal to 56/11 for the 2002 outburst QPO evolu- 0.01 0.01 −3 tion. The parameters obtained from this fit are a velocity −3 −4 of 3.7 ± 0.3 m/s and a mass of 8.05 ± 0.45 M . Further, −4 10 ⊙ we introduced an acceleration parameter to the model and 0 0 −1 −2 found that it fits the data properly with a χ close to one. −2 red 5 10 20 50 5 10 Energy (keV) Energy (keV) With the inclusion of the acceleration parameter the fit improved considerably giving a χ /dof equal to 12.44/10. Figure 3: Spectral fitting with two component flow model (Method Here the initial velocity obtained is 8.5 ± 0.65 m/s, the I) for the 2010 outburst of GX 339-4. −06 2 acceleration is −4.46 × 10 m/s and a mass range of 9.97 − 10.81 M . Except for 2004 rising phase data, it 0.1 is observed that the change in velocity is significant and MJD 56528.52 LHS hence the most general model should include shock front 0.01 acceleration as a parameter. The POS model fit results for 2002, 2004, 2007 and 2010 outbursts of GX 339-4 are −3 shown in Figure 5. For the 2002 outburst as QPO frequency increases from −4 0.16 Hz to 5.8 Hz the shock location changes from 278 r to 18.8 r and shock velocity decreases from 8.5 m/s to 6.0 m/s. Similar trend is observed in the rising phases of other outbursts as well. This implies that the shock front slows down as the system transits into SIMS during the −2 rising phase. We also found that during the decay phase 1 10 of 2007 outburst the QPO frequencies evolved from 4.1 Energy (keV) Hz to 0.67 Hz, the corresponding shock location increased 1 MJD 57092.63 from 25.0 r to 82.0 r and the shock velocity increased g g SIMS from 1.52 m/s to 2.00 m/s. So during the decay phase, the 0.1 shock front velocity increases as time progresses. In Table 0.01 2, we present the values of POS model parameters from −3 the QPO frequency evolution fitting. The overall mini- mum to maximum mass range obtained from this method −4 4 is 8.79 − 11.89 M . −2 −4 1 10 Energy (keV) Figure 4: Broadband spectral modelling with Method I for outbursts in 2013 (upper panel) and 2015 (lower panel) showing simultaneous XRT (black) and NuSTAR (red) spectra. The XRT spectral data used, is in the range 0.5 keV to 9.0 keV and the NuSTAR data is from 3.0 keV to 79 keV. 4.2. Results from Method II In the POS model, C-Type QPOs which are generally found in LHS and HIMS states of the rising phase of out- bursts are considered to be formed by the oscillations of −2 −1 −1 −2 −1 −1 χ Photons cm s keV χ Photons cm s keV −2 −1 −1 −2 −1 −1 (data−model)/error Photons cm s keV (data−model)/error Photons cm s keV −2 −1 −1 −2 −1 −1 χ Photons cm s keV χ Photons cm s keV Table 1: Fitted parameters from Two component Flow model. The XRT-NuSTAR spectra are denoted by †. All other observations are from RXTE. ˙ ˙ MJD Outburst State r (r ) m ˙ (M ) m ˙ (M ) Mass (M ) χ /dof s g h Edd d Edd ⊙ 52388.05 2002 LHS 274 ± 15 0.83 ± 0.01 1.70 ± 0.03 10.49 ± 0.12 111.42/70 52410.52 2002 HIMS 92.4 ± 4 1.12 ± 0.02 4.29 ± 0.09 09.59 ± 0.12 88.79/62 52416.59 2002 SIMS 28 ± 1 0.48 ± 0.002 3.37 ± 0.07 09.55 ± 0.07 113.81/62 52687.68 2002 HSS 25 ± 1 0.11 ± 0.001 3.81 ± 0.4 09.49 ± 0.45 63.94/60 53201.32 2004 LHS 136 ± 12 0.38 ± 0.007 0.30 ± 0.02 09.74 ± 0.32 97.30/90 53230.95 2004 HIMS 97.5 ± 1.5 0.27 ± 0.003 2.27 ± 0.04 08.89 ± 0.10 85.89/82 53322.68 2004 SIMS 10.6 ± 0.21 0.12 ± 0.007 6.50 ± 0.23 10.19 ± 0.63 66.27/49 53414.13 2004 HSS 17 ± 0.5 0.09 ± 0.001 6.69 ± 0.25 10.05 ± 0.86 13.95/10 55277.47 2010 LHS 453 ± 34 0.30 ± 0.005 2.9 ± 0.12 10.16 ± 0.16 125.72/117 55301.79 2010 HIMS 43 ± 1.3 0.29 ± 0.003 5.9 ± 0.17 08.97 ± 0.12 167.97/115 55315.69 2010 SIMS 23 ± 0.8 0.18 ± 0.002 9.28 ± 0.38 09.87 ± 0.27 125.15/79 55354.27 2010 HSS 5.7 ± 0.04 0.04 ± 0.001 7.10 ± 0.13 08.60 ± 0.30 14.46/13 56528.52 † 2013 LHS 137 ± 3 0.54 ± 0.003 0.13 ± 0.001 09.89 ± 0.06 1521/1296 57092.63 † 2015 SIMS 15 ± 0.1 0.23 ± 0.001 4.92 ± 0.02 10.17 ± 0.04 2580/1503 Table 2: POS model fitted parameters (Method II) 2 2 Outburst rs (r ) v (m/s) acceleration (m/s ) Mass (M ) χ /dof 0 g 0 ⊙ −06 2002 278 ± 4 8.50 ± 0.65 −4.46 × 10 9.97 − 10.81 12.44/10 −08 2004 143 ± 2 1.21 ± 0.09 +3.55 × 10 9.97 − 10.79 17.28/12 −06 2007 299 ± 7 12.32 ± 0.09 −9.00 × 10 10.04 − 10.49 15.87/20 −07 2007decay 25 ± 0.6 1.52 ± 0.06 +5.00 × 10 9.09 − 10.99 17.92/15 −06 2010 487 ± 22 11.72 ± 2.15 −4.79 × 10 8.79 − 11.89 14.42/13 tra. We fitted this correlation using empirical formulae de- veloped in ST07. We could not apply the Method III for 2002 outburst data of GX 339-4 as a matching reference POS modelling - GX 339-4 outburst is unavailable. So, we used the data from out- bursts in 2004, 2007 and 2010 where QPOs are observed for implementing this method. XTE J1550-564 with a dy- +1.02 namically estimated mass of 10.56 M (Orosz et al. −0.88 (2002)) and GRO J1655-40 with a dynamically estimated mass of 6.3±0.5 M (Greene et al., 2001) are the reference sources used. For the 2004 outburst as the QPO frequency increases from 0.30 Hz to 8.05 Hz the photon index (Γ) value changes from 1.45 to 2.21. The corresponding reference source used is the 2005 outburst of GRO J1655-40. During the 2007 -2 outburst Γ of GX 339-4 evolves from 1.74 to 2.76 as the 0 5 10 15 20 25 30 35 QPO evolves from its lowest to highest value. The ref- Days since QPO started erence source used is the 1998a outburst of XTE J 1550- 564. Finally the range of Γ for the 2010 outburst is 1.50 Figure 5: Evolution of QPO frequencies modelled with Method II to 2.56 which matches the reference outburst 1998b of the (POS model) for the rising phase of the outbursts of GX 339-4. The lower panel shows the residuals obtained after POS model fitting. source XTE J1550-564. ST09 has classified 1998 outburst of XTE J1550-564 into 1998a and 1998b based on photon index evolution during the outburst. 4.3. Results from Method III The fits corresponding to Method III are shown in Fig- As was described earlier the photon index parameter ure 6. For photon index - QPO frequency correlation, we obtained by modelling the energy spectra are correlated have used data from same energy ranges for both reference with the QPO frequencies of the corresponding power spec- χ QPO (Hz) and target. We have used 3 - 30 keV data for the 2004 and 2010 outbursts for both reference and target sources and 3 - 100 keV data for the case of 2007 outburst. This is 2.2 because the photon index values are significantly different over the two energy ranges and a reference outburst with 2.0 matching photon index evolution for the 2007 outburst of GX 339-4 was obtained only from the broadband spectral 1.8 analysis of both the target and reference sources. The re- sults from Method III presented in this paper are based 1.6 on the relations given in ST07 (see Equations 3 and 4). GX 2004 For the reference source GRO J1655-40 during its 2005 1.4 GRO 2005 outburst rising phase the parameter B value obtained from 0.05 1.0 10.0 +0.01 QPO (Hz) the correlation fitting is 0.15 . And for 2004 outburst −0.003 +0.02 of GX 339-4, we have B = 0.22 . Substituting these −0.01 values in Equation 4, we get the mass of GX 339-4 to be 2.8 in the range 8.50−10.85 M . Similarly we have computed the mass of the target source using data from different out- bursts. In Table 3, we have provided the parameter values 2.4 obtained by fitting the Gamma (Γ) - QPO (ν) correlation to different outbursts of GX 339-4 and the reference out- bursts used. Thus by using the constant M/B relation in 2.0 Equation 4, method III gives an over all mass estimate for GX 339-4 in the range 8.28 to 11.81 M . The mass 1.6 GX 2007 estimates for the different outbursts are given in Table 3. XTE 1998a 0.05 1.0 10.0 QPO (Hz) 2.6 2.2 1.8 GX 2010 1.4 XTE 1998b 0.05 1.0 10.0 QPO (Hz) Figure 6: Fitting of the QPO-photon index correlation based on Method III (ST07). The reference source is shown in blue and the target source in green. The reference sources are chosen such that the photon index saturation levels are similar to that of the target source. 4.4. Consolidated Results Most of the outbursting black hole sources show spo- radic jet activities (Fender et al., 2004, 2009; Radhika et al., 2016) and X-ray variabilities (Sreehari et al., 2018; Radhika et al., 2018) in the intermediate states. So we intend to get a tight constraint on the mass of GX 339- 4 using only LHS data. For this we use the results from method I, which is the most robust method among the three that we considered. As broadband LHS data is not available for 2015 outburst, we have used the only available SIMS data for this case. We have converted the outputs of Photon Index Photon Index Photon Index Table 3: Model fitted parameter values (Method III) Source - Outburst A B ν (Hz) D Mass (M ) tr ⊙ +0.34 +0.02 GX 339 - 2004 2.28 0.22 3.0 2.72 ± 0.08 8.50 − 10.85 −0.14 −0.01 +0.08 +0.035 GX 339 - 2007 2.73 0.266 3.8 0.42 ± 0.24 8.58 − 11.81 −0.06 −0.004 +0.09 +0.02 +0.69 GX 339 - 2010 2.51 0.40 2.0 1.81 8.28 − 10.49 −0.06 −0.03 −0.12 Reference - Outburst +0.01 GRO J1655 - 2005 2.29 ± 0.03 0.15 6.0 1.0 6.3 ± 0.5 −0.003 +0.031 +1.02 XTE J1550 - 1998a 2.79 0.285 4.2 1.0 10.56 −0.027 −0.88 +0.03 +0.25 +1.02 XTE J1550 - 1998b 2.35 ± 0.03 0.46 2.0 0.68 10.56 −0.02 −0.09 −0.88 the steppar command of XSPEC into a probability density function (PDF) following the steps in Nandi et al. (2018). The same method has been employed in Iyer et al. (2015) and Radhika et al. (2018) to combine different mass es- timates. Figure 7 shows the PDF corresponding to each outburst and we have also indicated the union and in- tersection of these PDFs. Each PDF is of unit area. The union PDF suggests a mass range of 9.72 - 10.70 M where as the intersection of PDFs gives a tighter constraint of 10.10 - 10.35 M . We consider the results from intersec- tion of PDFs, as per the naive Bayes rule for combining independent probabilities (Iyer et al., 2015). 2002 LHS 2004 LHS 2010 LHS 2013 LHS 4 2015 SIMS Intersection Union Figure 8: Forest plot of the mass estimates from different methods at 90% confidence. The size of the square represents the significance of each result. The maximum and minimum values from all three methods are used to plot the final range. 8.0 8.5 9.0 9.5 10.0 10.5 11.0 11.5 12.0 Mass (M ) As was mentioned earlier, we exclude the 2007 outburst Figure 7: The PDFs corresponding to mass estimates from spectral results from Method I due to large offset between PCA and modelling of each outburst is plotted here. Except for 2015 outburst, HEXTE spectra. For method II, we have used all rising we have used data from LHS state to obtain the PDFs. The yellow phase QPO evolution data available from RXTE. We have shaded portion is the intersection of individual PDFs and the cyan shaded portion is the union of individual PDFs. See text for more considered the decay phase data from 2007 as well for POS details and mass limits. modelling, but did not include the decay phase data from 2010 outburst (Nandi et al., 2012) as the QPO evolution is We also include together the results from all the three not smooth. And finally for method III, we have included methods in the Forest plot in Figure 8. The method and three outbursts for which we could find matching reference year are shown in the first column, the corresponding mass sources. For each individual method, we have considered ranges at 90% confidence are shown in the second column the minimum and maximum mass values as the range of and a plot of the mass range is shown in third column. the mass parameter. Thus the summary measure is given The size of the squares in the plot shows how significant as the minimum to maximum mass values considering all the results are. three methods together. So the mass of the black hole in GX 339-4 at a 90% confidence range is 8.28 − 11.89 M . https://cran.r-project.org/web/packages/forestplot/vignettes/forestplot.html PDF It should be noted that the PDFs corresponding to that the basic POS model with a constant velocity does Methods II and III will be broad (see Iyer et al. (2015)). not fit the temporal evolution of QPO perfectly. Hence This is also evident from the mass estimates in Tables we introduced an acceleration parameter that accounts for 2 and 3 which covers the range obtained from intersec- the change in velocity of the shock front. We found that tion of PDFs from Method I. Hence including PDFs from the shock front decelerates as the system approaches the Methods II and III would not improve the estimate ob- soft intermediate state during the rising phase. We also tained from the joint PDF (intersection) of Method I. The noted that the shock front accelerates as it recedes during mass estimate from the likelihood analysis (both intersec- the decay phase. During the rising phase of the outburst, tion and union of PDFs) lies within the extreme values particularly in LHS and HIMS, the increase of frequency obtained from the Forest plot. of QPO with time demands a forward movement of shock front. The location of shock in the accretion flow is deter- mined by the Rankine-Hugoniot conditions which states a 5. Discussion balance of both energy and pressure of the flow before and Method I is based on broadband spectral modelling us- after the shock. The key physical parameters are mass ing the two component accretion flow. For this we used accretion rate and viscosity of the flow. Accretion rate RXTE data in the range 3.0 - 150.0 keV and combined controls the rate of cooling and hence affect the total pres- XRT - NuSTAR data in the range 0.5 - 79.0 keV. From the sure whereas viscosity decides the transfer of angular mo- data of multiple outbursts, the mass of GX 339-4 estimated mentum. Of course in a magnetized accretion flow, the using this method to be 8.30 - 10.91 M . Though the magnetic pressure also needs to be taken into account. methods we employed are already used by other authors, An increase of accretion rate increases the cooling and this is the first time broadband spectral data has been shock front moves inward. The same trends are expected modelled to constrain the mass of the system. The results with the increases of viscosity as well as strength of mag- of spectral modelling from Method I makes it clear that netic field. In general, the dependency of these effects the error bars on mass parameter from RXTE data analy- (rate of cooling, pressure etc.) with the flow parameters sis (outbursts in 2002 to 2010) is significantly larger than are non-linear and hence an accelerated/decelerated move- the error bars on parameters obtained with XRT - NuS- ment of the shock front may result along with the evolu- TAR data. It can be seen that, the mass range obtained tion of the outburst. A constant speed of shock front is using Method I from 2013 outburst data is 9.82 − 9.95 M , a special case due to internal adjustment of the underly- where as the same method applied to 2015 outburst data ing physical processes. Sukov´a and Janiuk (2015) details gave 10.13 − 10.21 M . This is because these two data are the shock formation and numerical analysis of oscillating obtained when the source was in different states of the out- shock propagation in the scenario of super massive as well burst. The data systematics will be more in the case of the as stellar mass black holes. SIMS as compared to LHS due to presence of other activ- In Method II, there is uncertainty in the initial shock ity like relativistic jet ejections (Fender et al., 2004, 2009). location as one is not sure of the exact time at which the Also, an analytical spectral modelling for optically thin QPO has originated. Based on the observations, we con- (LHS state) or optically thick (HSS state) conditions pro- sider the time at which the lowest frequency is detected duce the required results, but in SIMS the central cloud’s (first instance of detection) as the zero time for modelling (i.e., Compton cloud) Thompson optical depth approaches QPO frequency evolution in the rising phase. However, unity where the analytically calculated spectrum may dif- this is unlikely to be accurate and hence the mass esti- fer from the expected one (say using monte-carlo simu- mate might also have some uncertainty. Moreover, in this lation) and this may also introduce some amount of un- method also, we have considered that the compression ra- certainty in the estimation. Besides this one should also tio is a constant value equal to 3 for all the outbursts. As consider that the RXTE data starts only from 3 keV on- this need not be the actual scenario, some uncertainty in wards, whereas the XRT - NuSTAR data starts from 0.5 the estimated mass is possible. keV onwards and hence we get better constraints from The third method used is the scaling based on pho- 2013 and 2015 outburst analysis. In this method, we have ton index - QPO frequency correlation. The over all mass fixed a constant value of 3 for the compression ratio which range obtained from this method is 8.28 − 11.81 M . The is within the theoretical limits of 1 to 4 for hydrodynamic disadvantage of this method is the requirement of an ap- shocks (Mandal and Chakrabarti, 2005). As compression propriate reference source whose photon index saturation ratio may vary as the shock location advances and the levels match with that of the target source. As even the black hole system transits over different states, it can in- same source does not have similar index saturation levels troduce some uncertainty on the estimated source mass. on different outbursts, it is not possible to find a perfect In the POS model (Method II), the QPO is considered reference source for this method. Moreover, even if the to originate from the oscillations of a propagating shock saturation levels match, the evolution of the photon index front. This model has mass of the black hole as a free pa- is significantly different for target and reference sources. rameter. Fitting the QPO frequency evolution using POS It must be noted that ST09 has obtained a higher mass model gives us a mass range of 8.79−11.89 M . We found estimate with a ΔM/M = 22%. ST09 modified equation 11 3 by introducing a parameter β to get, maximum mass value obtained from this method over all the outbursts considered is 8.30 M and 10.91 M . ⊙ ⊙     Mass determination using the method II, which is based 1 − on the POS model also gave us consistent results. Here we tr     Γ(ν) = A − BD ln exp + 1 . (5)     studied the temporal evolution of the fundamental QPO frequency as a function of the propagating shock location. As per this model, the C-type LFQPO frequency varies −3/2 After getting a good fit using Equation 5 for the refer- approximately as r . The mass range of the black hole ence source correlation, it is multiplied by a scaling factor obtained using this method is 8.79 − 11.89 M . along the frequency axis in order to get a new scaled curve The third method is based on the fact that the photon that matches the photon index-QPO correlation data of index of the energy spectrum has a correlation with the the target source. This scaling factor is further used to es- QPO frequency. As the QPO frequency increases the pho- timate the mass of the target source from the mass of the ton index also increases initially and then saturates. This reference source as mass and QPO frequency are inversely correlation has been empirically modelled by ST07 with an related. But this geometrical scaling of the fitted curve expression having four variables A, B, D and ν . We fit tr of the reference pattern does not give a perfect match on the correlation using Equation 3 for a reference source as the target patterns considered. Moreover, the scaled curve well as for the target source. The reference source should is quite different from the actual fitted curve of the tar- have similar photon index saturation levels to that of the get pattern. As a result, here we presented the results of target source during its outburst. Once the fit parameters method III based on ST07 rather than ST09. for both reference and target source are obtained, we use Using the three methods considered in this paper, we the scaling relation given in Equation 4 to estimate the have estimated the mass of GX 339-4 to be 8.28 − 11.89 M mass of the target source. Using this method, the mass with a ΔM/M of 35%. The ΔM/M for Parker et al. estimate of GX 339-4 at 90 % confidence gives the result (2016) including both negative and positive error is 31%. 8.28 − 11.81 M . But Parker et al. (2016) has used only data from the very The forest plot in Figure 8 shows all the mass estimates high state, while we have considered multiple outbursts we obtained. We quote the mass estimate based on the and all different states of the outbursts. Moreover, based three indirect methods to be in the range 8.28 − 11.89 M on a likelihood analysis we have obtained the most likely by considering the overall maximum and minimum val- mass of 10.10 - 10.35 M and a broader range of 9.72 - ues. Previous estimate by Parker et al. (2016) used only 10.70 M obtained from union of PDFs. Here ΔM/M is the data from very high state of the source and the re- 9.5%. sults partially overlaps with our estimate. Similarly the Corral-Santana et al. (2016); Tetarenko et al. (2016) mass determined by Heida et al. (2017) and Shaposhnikov provides the mass estimates of more than twenty compact and Titarchuk (2009) also partially overlaps with our re- objects in binary systems. Of this only GRS 1915+105, 4U sults. Through a likelihood analysis we also infer that the 1956+350, XTE J1550-564 and 4U 0540-697 have masses most likely mass for the source is 10.10 - 10.35 M . A larger than GX 339-4. So in the distribution of dynamical wider limit of 9.72 - 10.70 M is obtained from the union and candidate BH XRBs, GX 339-4 is a comparatively of PDFs. In method I, we have used data from multi- higher mass black hole. ple outbursts and also analysed different states of the BH X-ray binary system in each outburst. In method II, we have used the time evolution of C-type LFQPO frequen- 6. Conclusion cies in the LHS and HIMS state of the black hole binary. In this paper, we aimed at determining the mass of the In method III, we have used the LFQPOs in each out- Galactic black hole GX 339-4 using spectral and tempo- burst appearing in all the states. Thus using a complete ral techniques. We used three different methods to do the spectro-temporal analysis of the outbursts of GX 339-4 we mass estimation. Together these methods give a reason- have constrained the mass of the black hole in this system. ably narrow range of mass for the source. Among the BHBs whose mass has been estimated so far, The first method we used for mass determination is the GX 339-4 is a comparatively high mass black hole. two component accretion flow model wherein mass of the compact object is a free parameter. The other parame- Acknowledgements ters in this model are disc and halo accretion rates and shock location. We have tabulated the parameter values The authors thank both reviewers for their comments obtained by fitting this model to different outburst data and suggestions that helped to improve the quality of the in Table 1. It is usually observed that, the halo rate is manuscript. This research has made use of data obtained dominant in the hard states and the disc rate is dominant through the HEASARC web service, provided by the NASA/ in the soft states. Also the shock location which indicates Goddard Space Flight Center and also made use of data the size of the PSC is larger in the hard states and it supplied by the UK Swift Science Data Centre at the Uni- decreases as the spectra gets softer. The minimum and versity of Leicester. We thank GD, SAG; DD, PDMSA 12 and Director, URSC for encouragement and support to Giri, K., Chakrabarti, S.K., 2013. Hydrodynamic simulation of two-component advective flows around black holes. Mon. Not. carry out this research. R. Astron. Soc. 430, 2836–2843. doi:10.1093/mnras/stt087, arXiv:1212.6493. Grebenev, S., Sunyaev, R., Pavlinsky, M., Dekhanov, I., 1991. De- References tection of quasi-periodic oscillations of x-radiation from the black hole candidate gx339-4. Pisma v Astronomicheskii Zhurnal 17, Belloni, T., Homan, J., Casella, P., Van Der Klis, M., Nespoli, E., 985–990. Lewin, W., Miller, J., M´endez, M., 2005. The evolution of the Greene, J., Bailyn, C.D., Orosz, J.A., 2001. Optical and infrared timing properties of the black-hole transient gx 339–4 during its photometry of the microquasar gro j1655–40 in quiescence. The 2002/2003 outburst. Astronomy & Astrophysics 440, 207–222. Astrophysical Journal 554, 1290. Belloni, T.M., Motta, S.E., 2016. Transient Black Hole Bi- Heida, M., Jonker, P.G., Torres, M.A.P., Chiavassa, A., 2017. naries, in: Bambi, C. (Ed.), Astrophysics of Black Holes: The Mass Function of GX 339-4 from Spectroscopic Obser- From Fundamental Aspects to Latest Developments, p. 61. vations of Its Donor Star. Astrophysical Journal 846, 132. doi:10.1007/978-3-662-52859-4_2, arXiv:1603.07872. doi:10.3847/1538-4357/aa85df, arXiv:1708.04667. Casella, P., Belloni, T., Homan, J., Stella, L., 2004. A study of the Homan, J., Belloni, T., 2005. The Evolution of Black low-frequency quasi-periodic oscillations in the X-ray light curves Hole States. Astrophysics and Space Science 300, 107–117. of the black hole candidate XTE J1859+226. Astronomy and doi:10.1007/s10509-005-1197-4, arXiv:astro-ph/0412597. Astrophysics 426, 587–600. doi:10.1051/0004-6361:20041231, Homan, J., Wijnands, R., van der Klis, M., Belloni, T., van Paradijs, arXiv:astro-ph/0407262. J., Klein-Wolt, M., Fender, R., M´endez, M., 2001. Correlated Casella, P., Belloni, T., Stella, L., 2005. The ABC of Low-Frequency X-Ray Spectral and Timing Behavior of the Black Hole Candi- Quasi-periodic Oscillations in Black Hole Candidates: Analo- date XTE J1550-564: A New Interpretation of Black Hole States. gies with Z Sources. The Astrophysical Journal 629, 403–407. Astrophys. J. Suppl. Ser. 132, 377–402. doi:10.1086/318954, doi:10.1086/431174, arXiv:astro-ph/0504318. arXiv:astro-ph/0001163. Chakrabarti, S., Titarchuk, L.G., 1995. Spectral Properties Hynes, R.I., Steeghs, D., Casares, J., Charles, P., O’Brien, K., 2003. of Accretion Disks around Galactic and Extragalactic Black Dynamical evidence for a black hole in gx 339–4. The Astrophys- Holes. Astrophysical Journal 455, 623. doi:10.1086/176610, ical Journal Letters 583, L95. arXiv:astro-ph/9510005. Ingram, A., Done, C., 2011. A physical model for the con- Chakrabarti, S.K., Debnath, D., Nandi, A., Pal, P.S., 2008. Evo- tinuum variability and quasi-periodic oscillation in accreting lution of the quasi-periodic oscillation frequency in GRO J1655- black holes. Mon. Not. R. Astron. Soc. 415, 2323–2335. 40 - Implications for accretion disk dynamics. Astronomy and doi:10.1111/j.1365-2966.2011.18860.x, arXiv:1101.2336. Astrophysics 489, L41–L44. doi:10.1051/0004-6361:200810136, Iyer, N., Nandi, A., Mandal, S., 2015. Determination of the Mass of arXiv:0809.0876. IGR J17091-3624 from “Spectro-temporal” Variations during the Chakrabarti, S.K., Manickam, S.G., 2000. Correlation among Onset Phase of the 2011 Outburst. Astrophysical Journal 807, Quasi-Periodic Oscillation Frequencies and Quiescent-State Du- 108. doi:10.1088/0004-637X/807/1/108, arXiv:1505.02529. ration in Black Hole Candidate GRS 1915+105. Astro- Maccarone, T.J., Coppi, P.S., 2003. Hysteresis in the light physical Journal Letters 531, L41–L44. doi:10.1086/312512, curves of soft X-ray transients. Mon. Not. R. Astron. arXiv:astro-ph/9910012. Soc. 338, 189–196. doi:10.1046/j.1365-8711.2003.06040.x, Corral-Santana, J.M., Casares, J., Mun˜oz-Darias, T., Bauer, F.E., arXiv:astro-ph/0209116. Mart´ınez-Pais, I.G., Russell, D.M., 2016. BlackCAT: A cata- Mandal, S., Chakrabarti, S.K., 2005. Accretion shock sig- logue of stellar-mass black holes in X-ray transients. Astronomy natures in the spectrum of two-temperature advective flows and Astrophysics 587, A61. doi:10.1051/0004-6361/201527130, around black holes. Astronomy and Astrophysics 434, 839–848. arXiv:1510.08869. doi:10.1051/0004-6361:20041235. Das, S., Chattopadhyay, I., Chakrabarti, S.K., 2001. Stand- Motta, S., Munoz-Darias, T., Casella, P., Belloni, T., Homan, J., ing Shocks around Black Holes: An Analytical Study. 2011. Low-frequency oscillations in black holes: a spectral-timing Astrophysical Journal 557, 983–989. doi:10.1086/321692, approach to the case of gx 339-4. Mon. Not. R. Astron. Soc. 418, arXiv:astro-ph/0107046. 2292–2307. Debnath, D., Chakrabarti, S.K., Mondal, S., 2014. Implementation Munoz-Darias, T., Casares, J., Mart´ınez-Pais, I., 2008. On the of two-component advective flow solution in XSPEC. Mon. Not. masses and evolutionary status of the black hole binary gx 339-4: R. Astron. Soc. 440, L121–L125. doi:10.1093/mnrasl/slu024, a twin system of xte j1550-564? Monthly Notices of the Royal arXiv:1402.0989. Astronomical Society 385, 2205–2209. Done, C., Kubota, A., 2006. Disc-corona energetics in the very Nandi, A., Debnath, D., Mandal, S., Chakrabarti, S.K., 2012. Ac- high state of Galactic black holes. Mon. Not. R. Astron. cretion flow dynamics during the evolution of timing and spectral Soc. 371, 1216–1230. doi:10.1111/j.1365-2966.2006.10737.x, properties of gx 339-4 during its 2010–11 outburst. Astronomy & arXiv:astro-ph/0511030. Astrophysics 542, A56. Evans, P.A., Beardmore, A.P., Page, K.L., Osborne, J.P., O’Brien, Nandi, A., Mandal, S., Sreehari, H., Radhika, D., Das, S., P.T., Willingale, R., Starling, R.L.C., Burrows, D.N., Godet, Chattopadhyay, I., Iyer, N., Agrawal, V.K., Aktar, R., O., Vetere, L., Racusin, J., Goad, M.R., Wiersema, K., An- 2018. Accretion flow dynamics during 1999 outburst of XTE gelini, L., Capalbi, M., Chincarini, G., Gehrels, N., Kennea, J.A., J1859+226-modeling of broadband spectra and constraining Margutti, R., Morris, D.C., Mountford, C.J., Pagani, C., Perri, the source mass. Astrophysics & Space Science 363, 90. M., Romano, P., Tanvir, N., 2009. Methods and results of an doi:10.1007/s10509-018-3314-1, arXiv:1803.08638. automatic analysis of a complete sample of Swift-XRT obser- Orosz, J.A., Groot, P.J., van der Klis, M., McClintock, J.E., Garcia, vations of GRBs. Mon. Not. R. Astron. Soc. 397, 1177–1201. M.R., Zhao, P., Jain, R.K., Bailyn, C.D., Remillard, R.A., 2002. doi:10.1111/j.1365-2966.2009.14913.x, arXiv:0812.3662. Dynamical evidence for a black hole in the microquasar xte j1550– Fender, R.P., Belloni, T.M., Gallo, E., 2004. Towards a unified 564. The Astrophysical Journal 568, 845. model for black hole X-ray binary jets. Mon. Not. R. Astron. Parker, M.L., Tomsick, J.A., Kennea, J.A., Miller, J.M., Harrison, Soc. 355, 1105–1118. doi:10.1111/j.1365-2966.2004.08384.x, F.A., Barret, D., Boggs, S.E., Christensen, F.E., Craig, W.W., arXiv:astro-ph/0409360. Fabian, A.C., Fu¨rst, F., Grinberg, V., Hailey, C.J., Romano, Fender, R.P., Homan, J., Belloni, T.M., 2009. Jets from black hole X- P., Stern, D., Walton, D.J., Zhang, W.W., 2016. NuSTAR and ray binaries: testing, refining and extending empirical models for Swift Observations of the Very High State in GX 339-4: Weigh- the coupling to X-rays. Mon. Not. R. Astron. Soc. 396, 1370–1382. ing the Black Hole with X-Rays. Astrophys. J. Lett. 821, L6. doi:10.1111/j.1365-2966.2009.14841.x, arXiv:0903.5166. 13 doi:10.3847/2041-8205/821/1/L6, arXiv:1603.03777. servational aspects of outbursting black hole sources: Evolution of Radhika, D., Nandi, A., 2014. spectro-temporalcharacteristics and spectro-temporal features and X-ray variability. Journal of Astro- disk-jet connection of the outbursting black hole source xte j1859+ physics and Astronomy 39, 5. doi:10.1007/s12036-018-9510-0, 226. Advances in Space Research 54, 1678–1697. arXiv:1802.05163. Radhika, D., Nandi, A., Agrawal, V.K., Seetha, S., 2016. ‘Spectro- Sukova´, P., Janiuk, A., 2015. Oscillating shocks in the low temporal’ variabilities and possible physical mechanism for angular momentum flows as a source of variability of accret- jet ejections. Mon. Not. R. Astron. Soc. 460, 4403–4416. ing black holes. Mon. Not. R. Astron. Soc. 447, 1565–1579. doi:10.1093/mnras/stw1239, arXiv:1605.08351. doi:10.1093/mnras/stu2544, arXiv:1411.7836. Radhika, D., Sreehari, H., Nandi, A., Iyer, N., Mandal, S., Sunyaev, R.A., Titarchuk, L.G., 1980. Comptonization of X-rays in 2018. Broad-band spectral evolution and temporal variabil- plasma clouds - Typical radiation spectra. Astronomy and Astro- ity of IGR J17091-3624 during its 2016 outburst: SWIFT physics 86, 121–138. and NuSTAR results. Astrophysics & Space Science 363, 189. Svensson, R., Zdziarski, A.A., 1994. Black hole accretion doi:10.1007/s10509-018-3411-1, arXiv:1808.05556. disks with coronae. The Astrophysical Journal 436, 599–606. Remillard, R.A., McClintock, J.E., 2006. X-Ray Proper- doi:10.1086/174934. ties of Black-Hole Binaries. Annu. Rev. Astron. Astro- Tetarenko, B.E., Sivakoff, G.R., Heinke, C.O., Gladstone, J.C., 2016. phys. 44, 49–92. doi:10.1146/annurev.astro.44.051905.092532, WATCHDOG: A Comprehensive All-sky Database of Galactic arXiv:astro-ph/0606352. Black Hole X-ray Binaries. The Astrophysical Journal Supplement Samimi, J., Share, G.H., Wood, K., Yentis, D., Meekins, J., Evans, 222, 15. doi:10.3847/0067-0049/222/2/15, arXiv:1512.00778. W.D., Shulman, S., Byram, E.T., Chubb, T.A., Friedman, H., Titarchuk, L., Fiorito, R., 2004. Spectral Index and Quasi- 1979. GX339-4 - A new black hole candidate. Nature 278, 434– Periodic Oscillation Frequency Correlation in Black Hole Sources: 436. doi:10.1038/278434a0. Observational Evidence of Two Phases and Phase Transition Shakura, N.I., Sunyaev, R.A., 1973. Black holes in binary systems. in Black Holes. The Astrophysical Journal 612, 988–999. Observational appearance. Astronomy and Astrophysics 24, 337– doi:10.1086/422573, arXiv:astro-ph/0405360. 355. Titarchuk, L., Lyubarskij, Y., 1995. Power-Law Spectra as a Result Shaposhnikov, N., Titarchuk, L., 2007. Determination of black hole of Comptonization of the Soft Radiation in a Plasma Cloud. The mass in cygnus x-1 by scaling of spectral index-qpo frequency Astrophysical Journal 450, 876. doi:10.1086/176191. correlation. The Astrophysical Journal 663, 445, ST07. Wu, K., Soria, R., Campbell-Wilson, D., Hannikainen, D., Harmon, Shaposhnikov, N., Titarchuk, L., 2009. Determination of black hole B.A., Hunstead, R., Johnston, H., McCollough, M., McIntyre, V., masses in galactic black hole binaries using scaling of spectral and 2002. The 1998 Outburst of XTE J1550-564: A Model Based on variability characteristics. The Astrophysical Journal 699, 453, Multiwavelength Observations. The Astrophysical Journal 565, ST09. 1161–1168. doi:10.1086/324328, arXiv:astro-ph/0109222. Smith, D.M., Heindl, W.A., Swank, J.H., 2002. Two Different Zycki, P.T., Collin-Souffrin, S., Czerny, B., 1995. Accretion Discs Long-Term Behaviors in Black Hole Candidates: Evidence for with Accreting Coronae in Active Galactic Nuclei - Part One - Two Accretion Flows? The Astrophysical Journal 569, 362–380. Solutions in Hydrostatic Equilibrium. Mon. Not. R. Astron. Soc. doi:10.1086/339167, arXiv:astro-ph/0103304. 277, 70. doi:10.1093/mnras/277.1.70, arXiv:astro-ph/9505045. Sreehari, H., Nandi, A., Radhika, D., Iyer, N., Mandal, S., 2018. Ob-

Journal

AstrophysicsarXiv (Cornell University)

Published: Nov 11, 2018

References