EDP Sciences
Free Access
Issue
A&A
Volume 591, July 2016
Article Number A24
Number of page(s) 6
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/201628294
Published online 06 June 2016

© ESO, 2016

1. Introduction

In black-hole (BHXB) and neutron-star (NSXB) X-ray binaries, the transfer of mass from the optical companion to the compact star is mediated by an optically thick, geometrically thin, accretion disk. Because the mass-to-radius ratio in neutron stars and black holes is similar1, the characteristic velocities and dynamical time scales near the compact object are also similar (v ~ 0.5c, t ≲ 1 ms). Thus, it is not surprising that BHXB and NSXB show very similar X-ray spectral and temporal characteristics (van der Klis 1994, 2005; Maccarone & Coppi 2003; Barret 2004; Di Salvo et al. 2006; Muñoz-Darias et al. 2014).

It is generally accepted that the X-ray spectral continuum of BHXB and NSXB results from the sum of two main components: one thermal component which follows a blackbody distribution and another non-thermal component that is best described as a power law. The thermal component dominates at lower energies, hence it is referred to as the soft component, while the power-law component extends up to a few hundred keV and it is often known as the “hard tail” (Barret 2004; McClintock & Remillard 2006; Lin et al. 2007; Done et al. 2007; Church et al. 2012).

In BHXB, the soft thermal component is attributed to the emission from the innermost, hotter regions of the accretion disk, while the hard power law results from Comptonization of low-energy photons by energetic electrons. The physical origin of the Comptonizing medium could be an advection-dominated accretion flow (Narayan & Yi 1994; Esin et al. 1997), a low angular momentum accretion flow (Ghosh et al. 2011; Garain et al. 2012), or a radio jet (Band & Grindlay 1986; Georganopoulos et al. 2002; Reig et al. 2003; Giannios 2005; Markoff et al. 2005).

In NSXB, the emission spectrum is made more complex by the presence of a hard surface. Thus below ~10 keV, NSXB may display two soft components: a single-temperature blackbody from the neutron star surface or from the boundary layer between the disk and the neutron star and a multi-temperature blackbody from the accretion disk (Mitsuda et al. 1984; Gilfanov et al. 2003). The temperature of the neutron star photons is expected to be higher (kT ~ 2−3 keV) than that of the disk soft photons (kT ~ 1 keV), because the effective area of the neutron star surface is more compact than that of the accretion disk (Farinelli et al. 2007). Above ~20 keV, the most prominent spectral feature in the X-ray spectrum of NSXB is a power-law distribution (the hard tail), that usually extends up to 200300 keV without evidence of a cutoff (Di Salvo et al. 2000, 2001; Iaria et al. 2001; D’Amico et al. 2001; Farinelli et al. 2005; Migliari et al. 2007; Ding et al. 2011).

Both kinds of binary systems exhibit distinct spectral states in response to changes in the mass accretion rate, for example in the course of an outburst. The main states are generally referred to as soft and hard, depending on whether the bulk of the luminosity is radiated below or above a few keV, respectively.

Another property common to both BHXB and NSXB is the presence of radio jets (Fender et al. 2004; Migliari & Fender 2006). Optically thick, compact, steady radio emission is detected during the low (in the 120 keV band) X-ray lumonisity hard state, while optically thin radio flares occur during transitions from the hard to the soft sate. The jet is quenched in the soft state, when the X-ray luminosity approaches the Eddington limit (Paizis et al. 2006; Migliari et al. 2007; Miller-Jones et al. 2010).

A useful way to describe the rich phenomenology exhibited by these systems in the X-ray band is the hardness-luminosity diagram (HLD). In the HLD, BHXB and NSXB display distinct curves, along which the sources move smoothly as the mass accretion rate varies. In Kylafis & Belloni (2015), a physical inerpretation was offered for the phenomenology of BHXB along the q-shaped curve in the HLD. In NSXB, the shape of the curve in the HLD devides them into Z sources and atoll sources (Hasinger & van der Klis 1989; van der Klis 2006). The three branches that form the Z-shaped HLD are called horizontal (HB), normal (NB), and flaring (FB) branches. The HB corresponds to the hardest X-ray state, while the FB is the soft state. In the case of atoll NSXB, the branches are called island and banana branches, with subcategories such as extreme island and lower and upper banana. The extreme island branch corresponds to the hardest spectral sate. The hard tail has been detected in almost all the currently known Z sources (Di Salvo et al. 2002) and many of the atoll sources (Paizis et al. 2006).

The association of the compact radio jet and the X-ray hard tail with a specific region of the source in the HLD is well documented in BHXB (Fender et al. 2009) and NSXB (Paizis et al. 2006; Migliari & Fender 2006). Both the radio emission and the strength of the hard tail become weaker at higher accretion rates. Radio and hard X-rays show the strongest intensity in the hard states of BHXB, the horizontal branch of Z-NSXB (Hjellming et al. 1990; Migliari et al. 2007; D’Aí et al. 2007), and the island or lower banana states in atoll-NSXB (Migliari et al. 2003). In BHXB, radio emission is absent in the soft sate. In some NSXB, radio emission may be still detected in the softer states, albeit highly reduced (Homan et al. 2004).

While there is a general consensus that the radio emission is produced by a compact jet, the physical origin of the hard tail is not well understood. It naturally arises from Comptonization, however, the details of how the electrons acquire their energy (thermal, non-thermal, or bulk-motion) and the source of seed photons (neutron-star surface, inner accretion disk, synchrotron photons) remain unclear (Barret et al. 2000; Di Salvo et al. 2006; Paizis et al. 2006; Lin et al. 2007; Markoff et al. 2005; Farinelli et al. 2009; Revnivtsev et al. 2014). A weak hard tail contributing a small percentage of the total flux (≲ 3%) is also detected in the soft state, when the jet is quenched or absent (see, e.g., McClintock & Remillard 2006, and references therein). In this state, the power-law tail is thought to come from Comptonization in non-thermal flares above and below the thermal disk (Poutanen et al. 1997; Gierliński et al. 1999).

In a series of papers (Reig et al. 2003; Giannios et al. 2004; Giannios 2005; Kylafis et al. 2008; Reig & Kylafis 2015), we showed that Compton upscattering of soft photons from the accretion disk in the jet can explain a number of observational relations between the spectral and timing parameters in the hard state of BHXB. Our results clearly demonstrate that jets play a central role in all the observed phenomena, not only in the radio emission. Motivated by the similarities of the X-ray spectral continuum in BHXB and NSXB and especially by the clear connection of the presence of a hard tail with a radio-loud state of the source, we investigate whether our model can also account for the spectral properties of NSXB.

Our objective in this paper is to demonstrate that emission from the neutron-star surface plus Comptonization in a jet can reproduce the observed spectra of bright NSXB, namely a soft thermal component that is well described by a blackbody distribution with kT = 2−3 keV and a hard tail that follows a power-law with photon index in the range 1.8–3 with no evidence for a cutoff up to 200300 keV.

2. The model

The model that we have used in this work is identical to that used in Reig & Kylafis (2015). For our Monte Carlo code we follow Pozdnyakov et al. (1983). Photons from a blackbody distribution of temperature kTbb are injected at the base of the jet with an upward isotropic distribution. As the photons travel through the medium, they experience Compton scatterings with the spiraling electrons. In each scattering, the photons on average gain energy from the bulk motion (i.e., v and v) of the electrons. Comptonization can occur everywhere in the jet. The optical depth to electron scattering, the energy change and the new direction of the photons after scattering are computed using the corresponding relativistic expressions. Each model is run for 107 photons. More details of how the code works can be found in Kylafis et al. (2008).

Since we are not interested in reproducing the radio spectrum (see Giannios 2005), we have assumed mono-energetic electrons in the jet with Lorentz factor equal to the smallest in the distribution, namely (1)where v = v0 = const. is the terminal velocity of the jet and v is the smallest peprendicular velocity of the electrons in the lab frame.

The flow velocity in the jet is given by (2)where z1 and p are parameters. In other words, the jet has an acceleration region of thickness z1, beyond which the flow has constant velocity v0.

For a parabolic jet, i.e., one whose radius at height z is R(z) = R0(z/z0)1/2, the electron density is obtained from the continuity equation and for z>z1 it is inversely proportional to z.

The fixed parameters of our models and their reference values are: the radius R0 = 50 RNS of the base of the jet, where RNS = 1.25 × 106 cm is the radius of the neutron star, the distance z0 = 1 RNS of the bottom of the jet from the neutron star center, the height H = 105RNS of the jet, the terminal velocity v0 = 0.8c of the jet, the thickness z1 = 5 RNS of the acceleration zone, the exponent p = 1/2, and the temperature kTbb = 2.5 keV of the soft-photon input.

In our model, the parameters that we have varied are the Thomson optical depth τ along the axis of the jet and the minimum Lorentz factor γmin. Since v0 is a constant in our models, the variation of γmin is equivalent to a variation in v.

Because the jet is relativistic, the results also depend on the angle θ of observation with respect to the jet axis. Here we consider an intermediate range of observing angles 0.2 < cosθ< 0.6. Practically, for the Monte Carlo simulation this means that we count only photons that leave the jet in this range of angles.

thumbnail Fig. 1

Emergent spectra from the Monte Carlo simulations. The models plotted correspond to a fixed γmin = 3.34 (v = 0.52c) and optical depth τ = 0.1,0.25,0.5,0.75,1.0,1.25,1.5. The input source of photons follows a blackbody distribution with kT = 2.5 keV. The normalization of the input blackbody was chosen so that the hard component contributes about 5% of the 0.1300 keV luminosity. In the energy range covered by current missions (< 300 keV) the hard tail shows no break.

Open with DEXTER

thumbnail Fig. 2

Photon index and cutoff energy as functions of τ and γmin. In the left-hand panels, each point represents a calculation with the same γmin = 3.34 (v = 0.52c) and different τ. In the right-hand panels, each point represents a calculation with the same τ = 1 and different γmin. The scale of the Y-axis was left the same to faciliate the comparison.

Open with DEXTER

Table 1

Monte Carlo simulation results.

3. Results

The spectral parameters relevant to the present study are the photon index of the hard-tail component and the cutoff energy. We have investigated the dependence of these two quantities on the choice of two parameters of our model, namely the optical depth and the electron velocity (see Table 1).

Figure 1 shows the emerging spectra that result from our simulations using the reference values of the parameters reported in Sect. 2 with γmin = 3.34 (v = 0.52c) and various optical depths τ = 0.1, 0.25, 0.5, 0.75, 1.0, 1.25, and 1.5. As seen in the figure, although the spectra roll over above a certain energy, this energy is well above 200300 keV. Thus, up to 200300 keV, the spectra are essentially pure power-laws with photon index between 1.83, in agreement with the observations. Figure 2 shows the photon index Γ and the cutoff energy Ec as functions of optical depth τ and γmin (or equivalently v). The scale of the vertical axis was kept the same in the two plots to better assess the differences. The cutoff energy has a weak dependence on τ for the range of optical depths considered, but it is strongly dependent of γmin. This is expected because the cutoff is mainly determined by the energetics of the electrons (Giannios 2005). In contrast, the photon index strongly depends on τ, but only weakly on γmin. The reason is that the slope of the spectrum is very sensitive to the number of scatterings, which increases with optical depth.

Because of the significant flow velocity in the jet, the up-scattered photons escape preferentially in the forward direction. Thus we expect the results to also depend on the angle θ of observation with respect to the jet axis. The hardest spectra are found at small observation angles, as these are the photons that have suffered the most energetic scatterings. The cutoff energy increases as the observation angle decreases, but stabilizes at intermediate angles. However, once the angle range is fixed, the parameters follow a similar trend with optical depth and Lorentz gamma factor, irrespective of the range considered.

We also computed the time lag of hard photons in the energy range 7–40 keV with respect to softer photons in the range 27 keV. The time-lag spectrum, that is the time lag as a function of Fourier frequency, is shown in Fig. 3. The spectra follow a power-law dependence on Fourier frequency tlag ~ νβ, with β = 1.0 ± 0.1. This slope is somewhat steeper than the average slope measured in black-hole binaries, tlag ~ να with α = 0.7 ± 0.1 (Nowak et al. 1999; Pottschmidt et al. 2003; Cassatella et al. 2012), but agrees with the results reported by Olive & Barret (2001). We have not found a significant dependence of the shape of the lag spectrum with either the optical depth or γmin.

thumbnail Fig. 3

Time-lag spectrum that results from the model with τ = 1 and γmin = 3.34. All the remaining parameters of the model were fixed at the reference values (see Sect. 2).

Open with DEXTER

4. Discussion

We ran Monte Carlo simulations of the Compton upscattering of soft photons in an outflow moving at relativistic velocity (representing a radio jet) and computed the emerging X-ray spectrum. The input parameters of the model have been selected so that the results match what is seen in bright NSXB. Assuming a blackbody input source of photons with kT = 2.5 keV, our model reproduces very well the observed 1300 keV X-ray spectrum in the harder states, a power law continuum with a photon index Γ in the range 1.8–3 without evidence for a cutoff up to 200300 keV.

The main difficulty of the models that seek to reproduce the high-energy X-ray spectral continuum in NSXB lies in explaining the absence of a cutoff, because the spectrum resulting from Comptonization is expected to roll over at a certain energy (an effect known as Compton recoil), when the energy of the photons approaches that of the electrons. Various models have been put forward to explain the absence of a cutoff in the hard tails: bulk-motion Comptonization (Bradshaw et al. 2003; Paizis et al. 2006; Farinelli et al. 2007, 2008), Comptonization by a hybrid–thermal-non thermal corona (Gierliński et al. 1999; Özel et al. 2000; Zdziarski et al. 2001; Farinelli et al. 2005; Revnivtsev et al. 2014), or synchrotron emission from the electrons of a jet (Markoff et al. 2005). Because of the limited sensitivity of current detectors at energies above ~100 keV, statistically significant detections are given up to ~200 keV, typically. The spectral continuum could in principle still roll over at a certain energy provided that this cutoff energy is larger than the limit energy of the detectors2. Recently, Revnivtsev et al. (2014) reported the highest lower limit on the cutoff energy that has been detected so far: Ec> 330 keV, at a 2σ significance level.

This relatively high lower limit of the cutoff energy was used by Revnivtsev et al. (2014) to rule out the bulk-motion Comptonization model. The reason for this is that the cutoff energy in thermal and bulk-motion Comptonization depends on the temperature and the velocity of the Comptonizing electrons, respectively. For realistic values of kTe< 50 keV and assuming free-fall velocity vff ≈ 0.5c onto a neutron star of canonical mass (1.4 M) and radius (106 cm), the cutoff energy is expected to be below 200 keV. Moreover, the radiation pressure expected to result from the emission of the neutron-star surface may reduce the infall velocity to less than 0.2c, decreasing the bulk-motion Comptonization effect (Farinelli et al. 2008).

Observations collected over the past decade show that there exists a clear connection between the X-ray spectral states and the radio emission in NSXB. Similarly to BHXB, the NSXB are radio-loud in the hard state, usually displaying a flat spectrum indicative of a relativistic jet. The radio emission is quenched when the source enters the soft state in most (but not all) cases (Penninx et al. 1988; Hasinger et al. 1990; Hjellming et al. 1990; Oosterbroek et al. 1994; Migliari & Fender 2006; Migliari et al. 2007; Miller-Jones et al. 2010).

Motivated by the disputed origin of the hard X-ray tail in NSXB and by the clear relationship between the presence of this hard component and radio emission from the source in the form of a jet, we have investigated whether the observed X-ray spectrum of NSXB can be reproduced by our jet model. Our goal has been to reproduce the observed spectrum above ~20 keV, namely, a power-law distribution with a photon index in the range 1.83 and the absence of a high-energy cutoff up to an observed energy of 300 keV. In Sect. 3, we demonstrated that such a range of photon indices can be obtained with an optical depth varying in the range 0.1 ≤ τ ≤ 1.5 (Table 1). The lack of a cutoff in the observed spectra translates into cutoff energies above 200 keV in our models. These energies are obtained when the electrons move with γmin ≳ 3.2 (v ≳ 0.52c, for v = 0.8c). We note that, in principle, and contrary to the bulk-motion Comptonization model, there is no restriction on the velocity components of the electrons in the jet, provided that the parallel component is significantly larger than the perpendicular component. Therefore, although we selected models that give cutoff energies above but not too far from 300 keV, we could increase our limit on Ec simply by increasing γmin (see Fig 2). We also note that such high γmin is possible, since Fomalont et al. (2001), who studied the evolution of the radio emission in Sco X–1, found that v> 0.95c.

4.1. Comparison with black-hole binaries

There are two major differences between the jets in NSXB and BHXB: i) the radio luminosities of neutron-star jets are typically ~30 times lower than those of black-hole jets at comparable X-ray luminosities (Fender & Kuulkers 2001; Migliari et al. 2003); and ii) radio emission is not completely quenched in the soft state of some NSXB (Migliari et al. 2004), as opposed to their non-detection in BHXB, when they are in the same state.

Both of these differences have been addressed successfully in Kylafis et al. (2012), where the formation and the destruction of jets is explained using the Cosmic Battery (Contopoulos & Kazanas 1998; Contopoulos et al. 2006; Christodoulou et al. 2008).

That the radio luminosity is significantly lower in neutron-star jets than in black-hole ones can be understood as follows: while the radio luminosity scales as M for both types of compact object (here, is the mass accretion rate in units of the Eddington value and M is the mass of the compact object), the X-ray luminosity due to the ADAF-like inner accretion flow scales as M for neutron stars and as 2M for black holes (Narayan et al. 1997; Migliari & Fender 2006; Abramowicz & Fragile 2013). Thus, for comparable X-ray luminosities, the radio luminosity of neutron-star jets is ~ times smaller than that of black-hole ones.

The radio emission is quenched in the soft state of BHXB, because the Cosmic Battery works very inefficiently when the accretion flow is in the form of a Shakura-Sunyayev type (Shakura & Sunyaev 1973). On the other hand, in low-mass NSXB with weak magnetic field (B ~ 108−109 G), the Cosmic Battery can work moderately efficiently due to the “spreading layer” (Inogamov & Sunyaev 2010) that forms on the surface of the neutron star, provided that there is a significant difference between the spin frequency of the neutron star and the Keplerian frequency of the disk.

The findings from our calculations are consistent with the characteristics of neutron-star and black-hole jets discussed above. We have found that in NSXB the required optical depths are in the range 0.1 ≤ τ ≤ 1.5, while in BHXB they are in the range 1 ≤ τ ≤ 10 (Reig et al. 2003; Kylafis et al. 2008; Reig & Kylafis 2015). Because the optical depth is proportional to density and so is the radio emission, the lower optical depth in NSXB is consistent with the observational fact that the radio emission of the jet in NSXB is significantly smaller than that in BHXB.

To explain the fact that BHXB exhibit a high-energy cutoff in the X-ray spectrum while NSXB do not, we were forced to assume that γmin in the jets of NSXB is larger than in BHXB. In particular, we found from our calculations that for BHXB γmin ≲ 2.4 (Reig & Kylafis 2015), while for NSXB γmin ≳ 3. It will be interesting to see if this is confirmed by observations.

In a recent paper (Koljonen et al. 2015), an important anti-correlation was found between the break frequency νb in the jet spectrum (spectral energy distribution) and the power-law index Γ of the hard X-ray spectrum for BHXB and AGN. If the X-ray spectrum is produced in the jet by inverse Compton scattering (Reig et al. 2003; Giannios et al. 2004; Giannios 2005; Kylafis et al. 2008; Reig & Kylafis 2015), then this anticorrelation can be understood at least qualitatively. As the jet weakens, τ decreases, and as a result Γ increases (Reig & Kylafis 2015) and νb decreases. If our model is correct, that is to say if the hard X-ray spectrum is produced in the jet by inverse Comptonization, then we predict that NSXB as a class will exhibit lower values of νb than BHXB.

In our model, time or phase lags result from the random walk of the photons in the scattering medium. Hard photons scatter on average more times than softer ones before escaping the Comptonizing medium (the jet), hence they are delayed with respect to the softer ones. While the amplitude of the lags depends on the optical depth and dimensions of the jet (Giannios et al. 2004; Kylafis et al. 2008), the slope of the time-lag spectrum appears to be rather insensitive to those parameters or the Lorentz factor. In practice, the slope is determined for given energy ranges of the soft and hard photons and calculated for a given frequency range. Although most of our models produce slopes of ~1, smaller indices, similar to those of BHXB, could be obtained by different combinations of frequency and energy ranges. Given the lack of publications reporting Fourier time-lag spectra of NSXB, we cannot conclusively establish whether the different slopes that our calculations suggest corresponds to a distinguishing feature between NSXB and BHXB.

5. Conclusion

We ran Monte Carlo simulations of the Compton upscattering of soft photons in an outflow moving at relativistic velocity (representing a radio jet) and computed the emerging X-ray spectrum of bright NSXB. Our Monte Carlo simulations reproduce the observed spectral features: a power-law energy distribution at energies above 20 keV with photon index in the range 1.8–3 and the absence of a high energy cutoff up to ~300 keV. We demonstrate that the hard tail detected in these systems can be explained by Comptonization in the jet. The observed connection between the hard X-ray tail and the radio emission in low-mass X-ray binaries with neutron star companions supports this scenario. We explain the differences in the spectral and timing parameters between BHXB and NSXB by assuming that the electron population in the jets of NSXB have higher velocities and lower densities than in BHXB.


1

A 5M black hole has a Schwarzchild radius of ~15 km.

2

From the point of view of the energetics, if the slope of the power law is Γ ≤ 2, then the spectrum must show a rollover at a certain energy to prevent a diverging luminosity.

Acknowledgments

This research has been supported in part by the “RoboPol” project, which is implemented under the “ARISTEIA” action of the “Operational Program Education and Lifelong Learning” and is co-funded by the European Social Fund (ESF) and National Resources.

References

All Tables

Table 1

Monte Carlo simulation results.

All Figures

thumbnail Fig. 1

Emergent spectra from the Monte Carlo simulations. The models plotted correspond to a fixed γmin = 3.34 (v = 0.52c) and optical depth τ = 0.1,0.25,0.5,0.75,1.0,1.25,1.5. The input source of photons follows a blackbody distribution with kT = 2.5 keV. The normalization of the input blackbody was chosen so that the hard component contributes about 5% of the 0.1300 keV luminosity. In the energy range covered by current missions (< 300 keV) the hard tail shows no break.

Open with DEXTER
In the text
thumbnail Fig. 2

Photon index and cutoff energy as functions of τ and γmin. In the left-hand panels, each point represents a calculation with the same γmin = 3.34 (v = 0.52c) and different τ. In the right-hand panels, each point represents a calculation with the same τ = 1 and different γmin. The scale of the Y-axis was left the same to faciliate the comparison.

Open with DEXTER
In the text
thumbnail Fig. 3

Time-lag spectrum that results from the model with τ = 1 and γmin = 3.34. All the remaining parameters of the model were fixed at the reference values (see Sect. 2).

Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.