{'fraction_non_alphanumeric': 0.032105263157894734, 'fraction_numerical': 0.025146198830409357, 'mean_word_length': 4.5325137495956005, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Swift and RXTE observations of the black hole transient GX 339-4 in the hard state between outbursts

We use simultaneous Swift and RXTE observations of the black hole binary GX 339-4 to measure the inner radius of its accretion disk in the hard state down to 0.4% L_{Edd} via modeling of the thermal disk emission and the relativistically broadened iron line. For the luminosity range covered in this work, our results rule out a significantly truncated disk at 100-1000 R_g as predicted by the advection-dominated accretion flow paradigm. The measurements depend strongly on the assumed emission geometry, with most results providing no clear picture of radius evolution. If the inclination is constrained to roughly 20 degrees, however, the measurements based on the thermal disk emission suggest a mildly receding disk at a luminosity of 0.4% L_{Edd}. The iron abundance varies between roughly 1-2 solar abundances, with the i=20 degrees results indicating a negative correlation with luminosity, though this is likely due to a change in disk illumination geometry.


INTRODUCTION
Stellar mass black holes are typically observed in binary orbits with a companion star.
Accretion of material from the companion star onto the black hole forms an accretion disk, which converts gravitational energy into radiation primarily in the X-ray band. Accretion in black hole binaries (BHBs) can be a sporadic process, resulting in X-ray outbursts separated by periods of low luminosity. As a BHB evolves over time, it can be observed to transition between three canonical emission states: thermal, hard, and steep power law (Remillard & McClintock 2006). The thermal state is generally characterized by high luminosity, a thermally dominated spectrum, and low levels of X-ray variability. In the hard state, sources have spectra dominated by a power law component, relatively high X-ray variability, and are frequently accompanied by compact radio jets (Fender 2006).
While the hard state is typically seen at low levels of luminosity at the beginning or end of an outburst, a source can remain in the hard state even at high luminosities (Homan et al. 2001;Fender et al. 2004;Tomsick et al. 2005;Dunn et al. 2010). The steep power law is not well understood, occurring at high luminosities with strong variability in the form of quasi-periodic oscillations (van der Klis 2006) and a steep power law index. Note that BHBs are also frequently classified into states on the basis of model independent parameters such as spectral hardness and intensity (e.g. Belloni 2010).
The common framework for understanding BHB state transitions is described in Esin et al. (1997). Within this model, BHBs are usually found in the hard state at low accretion rates (ṁ 0.01M Edd ). The inner part of the accretion disk is optically thin and is characterized by a low density, possibly caused by radiative evaporation (Meyer et al. 2000).
Most of the viscously released energy is contained in the ions, which are unable to radiate efficiently in an advection-dominated accretion flow (ADAF). This results in energy being advected toward the black hole rather than being radiated away. Additionally, the flow becomes spherical due to the high pressures expected in the ion gas. At these low accretion rates, the thermally-dominated accretion disk is predicted to be truncated at an inner radius on the order of 100-1000 R g (R g = GM/c 2 ). At higher mass accretion rates (thermal state), the accretion disk reaches all the way down to the innermost stable circular orbit (ISCO). The precise luminosity and/or mass accretion rate at which the truncation occurs, if at all, is still a matter of debate (Miller et al. 2006a;Rykoff et al. 2007;Tomsick et al. 2008;Reis et al. 2010).
One way to test for the presence of an ADAF involves measuring the inner radius of the thermally-dominated accretion disk of BHBs at low luminosities. The thermal emission of the disk is dependent on both its inner temperature and inner radius, and modeling the shape of the thermal component in low luminosity BHB spectra can provide measurements of both. This is complicated, however, due to the power law domination of BHB spectra in the hard state; the low normalization and temperature of the thermal component are difficult to constrain. There are also uncertainties in corrections required by spectral hardening via scattering in the corona and general relativistic boundary conditions at the inner radius (Kubota et al. 1998;Shimura & Takahara 1995). Another method of radius measurement uses the shape of the relativistically broadened iron line (Fabian et al. 1989).
The effect of relativistic broadening as a function of disk geometry has been calculated by Laor (1991), but models based on these calculations fail to include thermal broadening inherent to the line. A model has been incorporated into XSPEC that includes thermal broadening (reflionx; Ross & Fabian (2005)) and also calculates the emission lines and continuum due to reflection of emission off of the disk in a physically self-consistent manner.
GX 339-4 is a galactic BHB with a persistent iron line in its various emission states.
The distance has been estimated to be 6 kpc < D < 15 kpc by Hynes et al. (2004) via analysis of the Na and Ca lines in the optical spectrum of the companion star. Hynes et al. (2003) calculated the mass function to be 5.8 M ⊙ using radial velocity data of the companion star. There is debate in the literature regarding the inclination, with some results at ∼ 20 − 30 • (Miller et al. 2006b;Done & Diaz Trigo 2010), and a more recent study favoring ∼ 50 • (Shidatsu et al. 2011). There is also debate on the inner accretion disk radius of GX 339-4 at low luminosities. Tomsick et al. (2009) have measured the radius at the lowest luminosity (0.14% L Edd ) and find a significantly truncated disk at > 65 R g . At slightly higher luminosities (∼ 2 − 4% L Edd ), however, several results indicate a disk radius consistent with the ISCO (Miller et al. 2006b;Reis et al. 2008), while Shidatsu et al. (2011) find a marginally truncated disk at 13.3 +6.0 −6.4 R g . This paper examines the inner disk radius of GX 339-4 using both thermal modeling and relativistic line broadening at a variety of luminosities in the hard state.

OBSERVATIONS AND DATA REDUCTION
We analyzed Swift X-ray Telescope (XRT; Burrows et al. (2005)) and respectively. The data was reduced using HEASOFT v6.11. PCA spectra were extracted using saextract, with marfrmf 3.2.6 used to produce the response matrices and the SkyVLE model to produce the background spectra. XRT spectra were extracted using xselect, with swxwt0to2s6 20010101v013.rmf used as the response and xrtmkarf 0.5.9 to produce the ancillary response file. Due to the proximity of GX 339-4 to the Galactic plane, observations from March 2001, February 2002, and September 2003 to create a quiescent spectrum that was subtracted from all PCA spectra. To improve statistics, data were combined into seven PCA/XRT spectra as indicated by Table 1. The PCA and XRT lightcurves are shown in Fig. 1 with boxes indicating spectral groupings.
Spectral fitting was performed in XSPEC v12.7.1. XRT spectra were restricted to 0.5-10 keV and PCA spectra to 3-25 keV. Systematic errors of 0.5% were added to the PCA spectra.

Spectral Model and Fitting Procedure
To investigate the possible presence of an iron line, we examined the residuals in the 5-8 keV range. Specifically, all seven spectra were fit with an absorbed multicolor blackbody disk model plus a power law component, the standard model for BHBs (Remillard & McClintock 2006). The data/model ratios in Fig. 2 show a clear, broad peak near 6-7 keV in the PCA data. The peak is also evident in the XRT residuals, except in spectra 3 and 4. The PCA residuals above 10 keV show a positive deviation from the power law, a feature typical of disk reflection. Confident in the presence of a broad iron line in all seven spectra, we switched to a model including iron line emission with relativistic blurring and disk reflection.   is evident in all spectra, as well as a reflection continuum above 10 keV.
PCA and XRT spectra were fit simultaneously using an absorbed (phabs) multicolor blackbody disk model (diskbb; Mitsuda et al. (1984)) plus a standard power law component. The Wilms et al. (2000) interstellar medium abundances and Balucinska-Church & McCammon (1992) photoelectric cross sections were used in the phabs component. Spectral features due to reflection of the power law component off of the disk were modeled using reflionx (Ross & Fabian 2005). reflionx assumes a folding energy of 300 keV in the power law, which was included with the multiplicative component highecut. The power law and reflection model components were convolved with kdblur (a relativistic broadening model using the calculations from Laor (1991)) in order to account for relativistic effects from the compact object.
To allow for a difference in PCA and XRT calibrations, we applied a normalization constant to the PCA model and left it free to vary during the fits. All diskbb and powerlaw parameters were left free. The physical reflionx parameters of ionization parameter and and iron abundance were left free, but the photon index was fixed to that of the powerlaw component for consistency. The redshift was set to 0 due to the galactic nature of the source, and the normalization was left variable. The kdblur inner radius parameter was left free, the outer radius was assumed to be 400 R g (R g = GM/c 2 ), and the disk inclination and emissivity index were varied between separate fitting rounds. As discussed in the introduction, there is a discrepancy between prior results on the inclination i of  Table 2. The rows give, in descending order, the inner disk temperature, diskbb normalization, kdblur inner radius, photon index, power law normalization, iron abundance, reflionx ionization parameter, reflionx normalization, best fit χ 2 , and model luminosity.

Inner Radius Determination
Inner radii were determined from the kdblur R in parameter and by using the radial dependence of the diskbb normalization parameter. The normalization parameter is defined as N DBB /cosi = r in (km)/d(10 kpc), where i is the inclination of the disk, r in is the inner radius, and d is the distance to the object. Physical radii obtained from N DBB must be converted into units of gravitational radii. We assume a distance of 8 kpc for convenient comparison to prior studies. The inclination was taken to be either 20 • or 50 • depending on the fitting round. In order to derive a physical inner radius, correction factors due to general relativistic boundary conditions and spectral hardening must be taken into account. Based on Kubota et al. (1998) and Shimura & Takahara (1995), we write the physical radius as R in = ξκ 2 r in , and have adopted a relativistic correction factor ξ = 0.412 and a spectral hardening factor κ = 1.7. Finally, we convert R in to units of gravitational radii R g = GM/c 2 . We take the black hole mass to be 10 M ⊙ , again for the purposes of convenient comparison.
To investigate the evolution of disk radius with luminosity, we calculate luminosities for each observation in units of Eddington luminosity. The Eddington luminosity is estimated using a black hole mass of 10 M ⊙ , and our source luminosity is estimated using the unabsorbed model flux from 1-100 keV and a distance of 8 kpc. Inner radii obtained from kdblur are plotted against source luminosity in Fig. 5 along with the relevant results from Table 6 of Shidatsu et al. (2011), and those obtained from diskbb are shown in Fig. 6.
Error bars are 90% confidence limits. All of these results used the same distance and mass assumptions as above to determine luminosity in units of L Edd .

Correlation Analysis
Two statistical tests were performed on the inner radii obtained from each fitting round in order to quantify the relationship the inner disk radius has with luminosity.
First, a one-tailed test was used to determine whether the radii were consistent with a constant. The data were fit to a constant, and the resulting χ 2 value was compared to the χ 2 distribution with the appropriate degrees of freedom (5 for diskbb, 6 for kdblur).
A χ 2 value in the upper tail of the distribution led to a rejection of the null hypothesis that the radius was constant. Second, the Pearson correlation coefficient distribution was simulated to determine how the radius depended on luminosity. For each observation and fitting round, 10 5 inner radii were generated based on the asymmetric 90% confidence regions and assuming Gaussian statistics. From this Monte-Carlo generated data, 10 5 Pearson correlation coefficients were calculated and binned into a histogram. Because it was noticed that the iron abundance has a large variation in all of the fitting rounds, these tests were also carried out on the iron abundance parameter. The P-value for rejection of the null hypotheses and the average correlation coefficients (ρ) with standard deviations for each fitting round are given in Table 3. To be clear, a low P-value is cause to reject the hypothesis of a constant parameter. Plots of the correlation coefficient distributions are shown in Fig. 7.

DISCUSSION
The inner radii obtained from the kdblur component depend strongly on the assumed disk geometry. It is clear from Fig. 5 that there is a general lack of correlation of the inner radii from our various fitting rounds. We attribute this to insufficient energy resolution and/or count rates to reliably distinguish between Compton/thermal broadening in reflionx and relativistic broadening in kdblur. Because of this, it is difficult to draw conclusions about the inner radius evolution based on the kdblur R in parameter without independent constraints on the inner disk geometry. While the Pearson correlation coefficient distributions based on kdblur R in shown in Fig. 7 for i = 50 • (top left) show strong negative correlation between inner radius and luminosity, suggestive of the beginnings of disk recession, only the i = 50 • , e = 3.0 fitting round has a P-value inconsistent with a constant radius (see Table 3 The stronger reflection at soft energies then suppresses the diskbb component by lowering its normalization, though in observation 5 the inner temperature is driven low instead. The dependence of inner radii obtained from diskbb on assumed disk geometry is much weaker. Our results shown in Fig. 6 show a similar functional form of inner radius between the fitting rounds. Note, however, the higher radius of observation 7 at i = 50 • due to the effect of kdblur described above. Because this observation occurs at the highest luminosity (5.21% L Edd ), it has a high impact on the correlation distributions and causes the i = 50 • , e = 2.1 and i = 50 • , e = 3.0 fitting rounds to be consistent with a constant radius.
The other four fitting rounds, all three i = 20 • rounds in particular, have P-values highly inconsistent with a constant radius and are negatively correlated at 99.9% confidence (none of the 1000 simulated ρ were positive for these three cases). Thus, our results at i = 20 • based on diskbb would constitute strong evidence of a mildly receding inner disk with independent inclination constraints.
The P-values in Table 3 all indicate a varying iron abundance. The strength of an iron emission line depends on the amount of iron in the source, but also on how that iron is illuminated (i.e. strength and incidence angle of the illuminating flux). Relativistic effects at the inner disk could produce a wide range of illumination geometries, including a varying emissivity index as pointed out in §3.1. Indeed, Bhayani & Nandra (2010) have investigated variability of relativistic iron lines in Seyfert galaxies using rms spectra, and frequently find enhanced variability with respect to the continuum in the red wings of iron lines that are thought to originate at the inner disk. The reflionx model does not account for a varying emissivity index, meaning our varying iron abundance could actually be indicating a change in this geometry. The iron abundances obtained in the i = 20 • fits are negatively correlated with luminosity, a phenomenon discovered in Seyferts by Nandra et al. (1997). This monotonic relationship of iron abundance is consistent with the change in disk geometry suggested by our i = 20 • diskbb inner radii. Nandra et al. (1997) suggest a more ionized disk at higher luminosities to explain this relationship, though our reflionx ionization parameter shows no signs of positive correlation with luminosity. At i = 50 • , the relationship of iron abundance with luminosity is weakly and ambiguously correlated.   Table 6 of Shidatsu et al. (2011). Blue points are from Tomsick et al. (2008), yellow points are from Shidatsu et al. (2011), teal points are from Done & Diaz Trigo (2010), the green lower limit is from Tomsick et al. (2009), the orange point is from Miller et al. (2006b), and the purple point is from Reis et al. (2008).