Particle transport in magnetized media around black holes and associated radiation
Key Words.:X-rays: binaries – radiation mechanisms: non-thermal – gamma-rays: general – neutrinos
Context:Galactic black hole coronae are composed of a hot, magnetized plasma. The spectral energy distribution produced in this component of X-ray binaries can be strongly affected by different interactions between locally injected relativistic particles and the matter, radiation and magnetic fields in the source.
Aims:We study the non-thermal processes driven by the injection of relativistic particles into a strongly magnetized corona around an accreting black hole.
Methods:We compute in a self-consistent way the effects of relativistic bremsstrahlung, inverse Compton scattering, synchrotron radiation, and the pair-production/annihilation of leptons, as well as hadronic interactions. Our goal is to determine the non-thermal broadband radiative output of the corona. The set of coupled kinetic equations for electrons, positrons, protons, and photons are solved and the resulting particle distributions are computed self-consistently. The spectral energy distributions of transient events in X-ray binaries are calculated, as well as the neutrino production.
Results:We show that the application to Cygnus X-1 of our model of non-thermal emission from a magnetized corona yields a good fit to the observational data. Finally, we show that the accumulated signal produced by neutrino bursts in black hole coronae might be detectable for sources within a few kpc on timescales of years.
Conclusions:Our work leads to predictions for non-thermal photon and neutrino emission generated around accreting black holes, that can be tested by the new generation of very high energy gamma-ray and neutrino instruments.
Some physical processes that take place near a black hole can be inferred from X-ray observations of accreting binaries. There is strong evidence that the X-ray emission in some of these systems is powered by accretion onto a black hole. Accretion processes around compact objects are modeled by hydrodynamic equations of viscous differentially-rotating fluids. The standard disk model is the most famous model based on this set of equations, and was developed by Shakura & Sunyaev (1973) and Novikov & Thorne (1973). The model predicts that the accreting gas forms a geometrically thin and optically thick disk, producing a quasi-blackbody spectrum due to thermal emission. The effective temperature of the accreting gas is in the range K, depending on the mass of the compact object and the accretion rate.
The observed X-ray spectrum, however, is too hard in many cases to have been produced by a standard disk alone, since the hard X-ray component corresponds to a temperature of K.
To explain the complete X-ray spectrum of Cygnus X-1 –which is arguably the most well-studied black hole candidate in the Galaxy– an extra component is usually added to the disk, the so-called corona (e.g., Dove et al. 1997; Gierlinski et al. 1997; Poutanen 1998). In this context, the corona is coupled to the disk by the magnetic field, and the plasma is heated by reconnecting magnetic loops emerging from the disk (Galeev et al., 1979).
As a result, the soft photons emitted by the disk gain energy by successive Compton upscatterings in the corona. This process is known as Comptonization, and it is the most accepted mechanism to explain the broad-band hard X-ray spectra of Galactic black holes. Another feature that can be explained by the presence of the corona, is a hardening of the spectra at keV, which is attributed to the Compton reflection of hard radiation (emitted by the corona) from a cold material, i.e. the disk (White et al., 1988; George & Fabian, 1991).
Since their early applications to Cygnus X-1 (Shapiro et al., 1976; Ichimaru, 1977), different models have attempted to explain the complete X-ray spectrum of black holes in binary systems. The observational data of the spectra are used to set some constraints on the geometry of the source. For example, the covering fraction of the hot cloud as viewed from the soft photon source can be estimated from the observed spectral slopes, and the reflection bump limits the solid angle subtended by the disk around the corona (Poutanen et al., 1997). There are two specific geometries that seem to reproduce well most of the observed spectra: the disk + corona ‘sombrero’ model and the advection dominated accretion flow model (ADAF).
In the ‘sombrero’ model, the corona is added to the disk. It is usually assumed that the corona can be represented by a homogeneous spherical cloud of radius around the compact object; the cold disk (with an inner radius ) is truncated at a certain distance from the compact object in the hard state of Galactic black holes (Dove et al. 1997; Poutanen 1998). The ratio depends on the amount of Compton reflection observed, and can be obtained from the energy balance and the electron/positron pair balance (Dove et al. 1997; Poutanen et al. 1997). In the case of Cygnus X-1, this ratio is estimated to be , which means that only a short fraction of the disk is within the hot cloud (Poutanen et al., 1997).
On the other hand, the ADAF model is a self-consistent solution of the hydrodynamic equations of viscous rotating flows (see, e.g., Ichimaru 1977; Narayan & Yi 1994, 1995a, 1995b; Abramowicz et al. 1995). The main characteristic of this solution is that most of the energy is accreted onto the compact object and the gas is unable to cool efficiently, mainly because of the low density. Since the plasma is a poor radiator, the viscous energy is stored in the gas as thermal energy and the gas temperature becomes high. This causes the accreting gas to swell, hence the geometry of ADAFs is quasi-spherical. This geometry is similar to the geometry of the corona+disk ‘sombrero’ model already mentioned above.
Besides the X-ray spectra observed in X-ray binaries (XRBs), sources such as Cygnus X-1 produce steady emission up to a few MeV (McConnell et al., 2000; Cadolle Bel et al., 2006), that is indicative of a non-thermal contribution to the spectral energy distribution (SED). Li et al. (1996) treated simultaneously the transport and acceleration mechanisms for leptons obtaining good fits to the spectra of Cygnus X-1 observed by COMPTEL. In addition, Li & Miller (1997) demonstrated that under the physical conditions expected in the surroundings of accreting black holes, it is possible to accelerate electrons out of a Maxwellian distribution, resulting in a non-thermal tail.
More recently, several works have been devoted to studying the effects of non-thermal populations of electrons in Galactic black hole coronae (e.g., Belmont et al. 2008; Malzac & Belmont 2009; Vurm & Poutanen 2009), as well as in accreting supermassive black holes (e.g., Belmont et al. 2008; Veledina et al. 2011). These works, however, have not studied the effects of a hadronic component in the source.
In Romero et al. (2010c), it was shown that the presence of hadrons can also explain the non-thermal emission detected in Galactic black holes. In the present work, we show the results of a theoretical study of the effects of the injection of non-thermal particles, both electrons and protons, in a magnetized corona around a black hole, that has many refinements compared to previous works. In particular, a self-consistent treatment of photon and particle transport is now presented. We solve the set of coupled kinetic equations for all types of particles, including photons, hence the treatment of absorption is straightforward. We attempt to estimate the SED produced by these relativistic particles, and explain the origin of the non-thermal tail observed in some XRBs (McConnell et al., 2000; Cadolle Bel et al., 2006; Jourdain et al., 2012).
The completion of IceCube opens new possibilities for neutrino detection, so we also study the neutrino production in the source.
The structure of this paper is as follows. In the next section, we present the basic model used to obtain the values of the relevant parameters, and fix the characteristics of the medium where the non-thermal particles are injected. In Sec. 3, we discuss the treatment of the physical processes, along with the details of the numerical method used to solve the corresponding equations. We first study the system in a steady state, to test our code and reproduce the results obtained in Romero et al. (2010c). In Sec. 4, we then consider the case of a transient event, including the time-dependence of the transport equations. Finally, we estimate the electromagnetic emission as well as the neutrino flux.
2 Basic model
In accordance with estimates for Cygnus X-1, we assume a black hole of mass (Orosz et al., 2011). We study the system in the low-hard state, which is the state where the X-ray spectrum is dominated by the coronal emission.
The size of the region where the hard radiation is produced, is limited by the variability observed in the spectra of Galactic black holes. The minimum variability timescale is on the order of milliseconds, hence in the hard state the corona lies within (Poutanen, 1998), where is the gravitational radius (). We considered a spherical corona with a size of (see Fig. 1), and assumed that the luminosity of the corona is 1 % of the Eddington luminosity (Esin et al., 1997), which results in erg s.
In ADAF models, ions and electrons interact only through Coulomb collisions and there is no coupling between the two species. In this case, the plasma has two temperatures, with the ion temperature ( K) being much higher than the electron temperature ( K) (e.g., Narayan & Yi 1995a, b).
We initially consider that the system is in steady state, thus the equipartition of energy between the different components of the system is a reasonable assumption (Esin et al., 1997). The power of the jet observed in the low-hard state of X-ray binaries is related to the magnetic field, which in turn is coupled to the corona. We then assume equipartition between the magnetic energy density and the bolometric photon density of the corona (as done by, e.g., Bednarek & Giovannelli 2007). At the same time, since the magnetic field launches the plasma into a jet, we also assume equipartition between the magnetic and kinetic energy densities (see e.g., Zdziarski 1998). These conditions allow us to obtain the values of the main parameters. Specifically, the conditions are
where is the value of the mean magnetic field in the corona and , are the electron and ion plasma densities, respectively.
The X-ray emission of the corona is characterized by a power law in photon energy with an exponential cut-off at high energies (e.g. Romero et al. 2002)
We adopt and keV, which are typical of Cygnus X-1 (Poutanen et al. 1997). The photon field of the accretion disk is modeled as a blackbody of temperature keV. Both the X-ray emission of the corona and the radiation field of the disk are considered as seed photon sources for Compton scattering and photomeson production in relativistic particle interactions.
In ADAF models, particles are advected onto the compact object at a mean velocity of . The accretion timescale is about a second, which is considerably shorter than the cooling rate for hadrons (Romero et al., 2010c). As a result, only high energy protons are able to lose a significant amount of energy, and most of the power injected in terms of relativistic hadrons falls into the black hole. Since most of the electron/positron pairs are produced by hadronic interactions, in advective-dominated coronae the high energy emission is likely produced in the jet (Markoff et al., 2001; Vila & Romero, 2010; Vila et al., 2012).
In Romero et al. (2010c), a static corona model was shown to be capable of producing the non-thermal emission of Cygnus X-1. In this work, we consider only this possibility. In a static corona, which is supported by magnetic pressure, the relativistic particles can escape mainly by diffusion (Beloborodov, 1999; Nayakshin, 1999). In the Bohm regime, the diffusion coefficient is , where is the gyroradius of the particle. The diffusion rate is
We consider the injection of the non-thermal particle distributions of both electrons and protons in this corona. The mechanism of particle acceleration in black hole coronae is likely related to magnetic reconnection, which is essentially a topological reconfiguration of the magnetic field caused by a change in the connectivity of the field lines. Several works have been published on particle acceleration through magnetic reconnection (e.g. Schopper et al. 1998; Zenitani & Hoshino 2001; de Gouveia Dal Pino et al. 2010) The basic idea is that a first-order Fermi mechanism takes place within the reconnection zone caused by two converging magnetic fluxes of opposite polarity that move toward each other with a velocity (de Gouveia dal Pino & Lazarian, 2005). The resulting injection function of relativistic particles is a power-law with an index .
The detailed analysis of Drury (2012) supports the idea that magnetic reconnection can lead to type I Fermi acceleration in a similar way as diffusive acceleration is driven by shocks. According to Drury, the spectral index does not have a universal value of ; it instead lies somewhere in the range . Despite this caveat, and given the uncertainties in the acceleration models, we adopt the value of , which is consistent with all simulations implemented so far.
In our model, the acceleration mechanism is not included as a term of the energy gain in the transport equation, but used to fix the injection function of primary electrons and protons and determine the maximum energies that relativistic particles can achieve.
As for standard first order Fermi acceleration, the acceleration rate for a particle of energy in a magnetic field is given by
In the corona, the Alfvén speed is , yielding an acceleration efficiency of .
The power available in the magnetized system can be estimated by (del Valle et al., 2011)
where . This yields a power available for non-thermal processes of % . The total power injected into relativistic protons and electrons, , is assumed to be a fraction of the luminosity of the corona, , with . The way in which energy is divided between hadrons and leptons is unknown. To deal with this uncertainty, it is useful to define the parameter as the ratio of the power injected in protons to the one injected in electrons, . Following Romero et al. (2010c), we consider models with (proton-dominated case, as for Galactic cosmic rays) and , that is, models with the same power injected in both hadrons and leptons. In our model, the injection function is assumed to be both homogeneous and isotropic.
We model the escape of radiation from the region using the treatment described in Coppi (1992)
, is the Klein-Nishina cross-section, is the number density of pairs defined by , and represents an average over the particle distribution such that photons with energy escape out of the source in a time . Table 1 summarizes the values of the relevant parameters in our model.
|: black hole mass ||Value for Cygnus X-1 (Orosz et al., 2011).|
|: corona radius |
|: inner disk/corona ratio||This is the typical value in models where the main source of seed photons for Comptonization are provided by the disk (Poutanen, 1998; Haardt & Maraschi, 1993). In our model, there are other relevant sources of photons that can relax the condition ; however, we only use this parameter to constrain the size of the corona, which for Cygnus X-1 cannot differ significantly from the one adopted in our work (Poutanen, 1998)|
|: electron temperature [K]|
|: ion temperature [K]|
|: X-ray spectrum cut-off [keV]|
|: X-ray spectrum power-law index|
|: acceleration efficiency|
|: magnetic field [G]|
|: plasma density [cm]|
|: disk characteristic temperature [keV]|
3 Treatment of radiative processes
Computing the SEDs of black hole coronae, as well as other magnetized plasmas, is a complex task, since it must include, among other issues, a detailed knowledge of the plasma characteristics and the microphysical processes. The first method used to treat this problem was a Monte Carlo simulation (e.g., Aharonian et al. 1985; Stern et al. 1995; Pilla & Shaham 1997). The main problem for this approach is the small number of high-energy photons, which leads to low quality photon statistics. On the other hand, it is usually easy to model the radiative transfer processes (see Pellizza et al. 2010 for a three-dimensional code).
A second method for estimating the spectra of compact sources involves solving the kinetic equations (e.g., Lightman & Zdziarski 1987; Coppi & Blandford 1990; Coppi 1992). The different interactions of particles with the fields of the source ensure that it is unavoidable to study a wide range of particle energies, hence the computation of radiative processes is quite complicated using this approach. The main advantage of this method is that the transport of photons is solved self-consistently. Since, in general, these equations are solved numerically, the availability of computational resources over the past decade has allowed the significant improvement of this approach (see, e.g., Aharonian & Plyasheshnikov 2003; Malzac & Belmont 2009; Poutanen & Vurm 2009; Vurm & Poutanen 2009).
To complete the treatment initiated in Romero et al. (2010c), we take the second approach. As mentioned in the previous section, we are interested in the study of the injection of non-thermal particle distributions of electrons and protons in the system. Once protons are injected into the corona, they interact with both the photon and matter fields, producing pions. In addition, the charged pions decay producing muons, so we also take into account the presence of these transient particles.
An accurate description of a hot, magnetized plasma such as the corona should also treat the processes of pair production and annihilation, hence we include electron/positron pairs.
The main channel for secondary pair production in our model is photon-photon annihilation. The most important background photon field for pair creation is the thermal X-ray radiation of the corona. The energy spectrum of pairs has been studied, for example, by Aharonian et al. (1983) and Böttcher & Schlickeiser (1997). Under the conditions , the pair emissivity (in units of erg s cm) can be approximated by the expression
Here is the Lorentz factor of the electron, , and are the dimensionless photon energies.
Another important source of electron/positron pairs when protons are present is the Bethe-Heitler process. To estimate this contribution, we use the treatment described by Chodorowski et al. (1992).
3.1 Radiative losses
Loss terms in our equations include synchrotron radiation, inverse Compton (IC) scattering, and relativistic bremsstrahlung for electrons and muons. We also consider photon production by pair annihilation. For protons, the relevant mechanisms are synchrotron radiation, photomeson production, and hadronic inelastic collisions.
A complete discussion of cooling times due to synchrotron radiation, IC scattering, and hadronic interactions can be found, for example, in Vila & Aharonian (2009) and Romero et al. (2010c). The injection of secondary particles, such as pions and muons, is also discussed in these works.
The number of pairs is modified by both creation and annihilation. A useful approximation to the annihilation rate is given by Coppi & Blandford (1990),
where represents the electron/positron distribution (in units of erg cm).
In Fig. 2, we show the cooling rates together with the diffusion, decay, and acceleration rates of all particle species. The value adopted for is twice the rest-frame energy of each type of particle. The maximum energy for electrons and protons can be inferred using a balance between the acceleration rate, given by Eq. 5, and the cooling rate. This yields GeV for electrons, and TeV for protons. Particles of such energies satisfy the Hillas criterion, and can be contained within the source.
The cooling rates for the IC scattering and photomeson production shown in the figures are the result of the interactions of particles with two target photon fields: the X-ray emission of the corona and the emission from the disk.
Cygnus X-1 is a binary system with a massive star that produces an intense radiation field. This photon field could be considered as an additional target for IC scattering since it dominates the bolometric luminosity of the source. However, its effect on particle energy losses is negligible. This can be shown by a simple analysis: in the Thompson regime for IC scattering, the cooling rate () is proportional to the energy density of the target photon field, defined by
This magnitude for both the accretion disk and the corona X-ray photon fields is on the order of magnitude erg cm, whereas for the stellar field at the location of the corona it is erg cm. The photon field provided by the star is then an ineffective seed target for IC processes compared with the local X-ray field.
As can be noted from the figures, no unique mechanism clearly dominates the energy losses for a given particle species. For electrons and muons, there are two relevant radiative processes, IC scattering and synchrotron radiation, whereas for protons and charged pions hadronic interactions are also important. Since IC scattering and photomeson production are processes that depend on the radiation field and, at the same time, the photon field is affected by all the interactions of particles with the different fields, the particle cooling times and distributions are strongly coupled with the transport of photons. This point implies that it is essential to improve the model presented in Romero et al. (2010c), solving the system of coupled kinetic equations for all type of particles.
3.2 Set of coupled equations
We determine the relativistic particle and photon distributions solving the set of coupled transport equations in the steady state and assuming spatial homogeneity and isotropy. The set of kinetic equations are
Transport of electron/positron pairs and protons:
where , , .
Transport of charged pions and muons:
where , , , .
Transport of photons:
Here, represents the steady state of each particle distribution (in units of erg cm), includes all radiative losses for a given type of particle, is the timescale over which relativistic particles escape from the system, is the mean decay time for transient particles (pions and muons), and is the injection function.
In Eq. (17), the term represents photon injection due to several radiative processes
where , , and give the contribution from synchrotron radiation, IC scattering, and neutral pion decay to photon injection, respectively.
for or , or in any other case:
and and . We refer to Romero & Paredes (2011) and references therein for formulae on radiative processes.
3.3 Numerical method
We use an Adams-Moulton method (see, e.g., Press et al. 1992) to solve the differential equations in Eqs. (15)-(16). This is an implicit multi-step integration method that can reach higher orders than other numerical algorithms; we use in particular a second order method.
Following the scheme described in Vurm & Poutanen (2009), we define an equally spaced grid on a logarithmic scale for the energy of particles
We then obtain a system of linear algebraic equations of the form
with the boundary condition that , which represents . The matrix contains the particle losses, whereas particle injection is included in the vector
where is the energy step.
We first solve the transport equations and obtain the particle distributions. These are used to estimate to first order the non-thermal luminosity. Once we know the non-thermal photon injection, we solve Eq. 17. An important property of the photon transport equation (17) is its nonlinearity. This is because the cross-section of photopair production () depends explicitly on the photon distribution. We use the approximation discussed in Poutanen & Vurm (2009), which consists in taking the photon distribution from a previous step, , to obtain the current injection (step ) of electron/positron pairs. Since we firstly consider a steady state, the way to solve the photon transport equation is then reduced to a simple iterative scheme given by
The updated photon distribution is then added to the background photon fields (corona power-law plus emission of the disk)222As we have mentioned, at the location of the corona the stellar photon field can be considered as negligible. to compute the IC scattering and hadronic interactions. We calculate the radiative losses and injection of particles where appropriate. The transport equations of massive particles are then solved, and the new distributions are used to compute the luminosity to second order. The process is repeated until all particle distributions converge to a stationary value.
3.4 Photon absorption in the stellar radiation field
The binary system Cygnus X-1 is composed of a massive star and a compact object. The massive star is an O9.7 Iab star of (Orosz et al., 2011). The orbit of the system is circular, with a period of 5.6 days and an inclination of between and (Orosz et al., 2011).
The star produces an intense radiation field that can absorb gamma rays by pair creation within the binary system. The photon field of the star is anisotropic, because its intensity depends on the position of the black hole in its orbit (but outside the corona). The gamma ray absorption in X-ray binaries with a massive companion star has been studied, for example, by Herterich (1974), Carraminana (1992), Bednarek (1993, 2000), Böttcher & Dermer (2005), Dubus (2006), Zdziarski et al. (2009), and Romero et al. (2010a). To estimate the observable spectrum and compare with the available data, it is necessary to include an appropriate treatment of the absorption in the stellar field. We then use the approach described in Romero et al. (2010a), where the case of Cygnus X-1 is considered.
The star has a radius cm, and we assume a blackbody radiation density of temperature K. The orbital radius is cm. Table 2 lists the values of the companion star and orbit parameters (Orosz et al., 2011).
|: Star mass |
|: Star radius |
|: Star temperature [K]|
|: Orbital period [days]|
|: Semi-major axis [cm]|
|: inclination angle |
In Fig. 3, we show the opacity map produced by the photon absorption in the stellar field, in the relevant energy range and along the complete orbit ( is the orbital phase in units of ; we note that corresponds to the compact object at opposition, i.e. superior conjunction). As can be seen in the figure, except for close to the inferior conjunction (compact object in front of the star, ), the emission is completely suppressed by the stellar photon field at energies in the range GeV GeV. To better illustrate this effect, Fig. 4 shows the modulation with the orbital phase of gamma ray emission of the corona in a steady state at a given energy. The stellar photon field is almost transparent when the compact object passes in front of the star (), but almost opaque when the compact object passes behind (). These results agree with those obtained by Romero et al. (2010a).
3.5 Spectral energy distribution
Figure 5 shows the main non-thermal contributions to the total luminosity. The internal absorption is not included in this plot. The synchrotron radiation of electron/positron pairs dominates the spectrum for MeV. We note that given the small size of the corona, the synchrotron radiation below eV is self-absorbed. All radio emission of the source comes from the jet (Stirling et al., 2001).
The relativistic protons injected in the corona collide with thermal protons through different channels ( and are the multiplicities)
and with the photon field
The main electromagnetic result of these interactions is the neutral pion decay
This is the most relevant source of photons in the high-energy gamma-ray band (see Fig. 5).
One assumption of our model is the equipartition of energy between the magnetic energy density and the photon energy density of the corona. It is then expected that the contribution of synchrotron radiation and IC scattering to the total luminosity be comparable. In the analysis of Fig. 2, we showed that IC scattering and synchrotron radiation are the radiative processes dominating the energy losses for electrons at low energies. At higher energies, the Klein-Nishina effect becomes important, the IC cross-section decreases, and this leaves synchrotron as the main mechanism causing electron energy loss. The Klein-Nishina effect is also responsible for the diminution of the IC radiation with respect to the Thompson regime; the gamma-ray flux is proportional to the number of interactions, hence when the IC cross-section decreases, so does the gamma ray emission. This explains why in Fig. 5 synchrotron radiation dominates the luminosity of the source at low energies. Nevertheless, the IC emission of electron/positron pairs is comparable to the synchrotron radiation at eV.
In Fig. 6, we show the total photon flux produced in the corona in two different orbital phases, superior and inferior conjunction. These positions correspond to the maximum and minimum absorption in the stellar photon field, respectively. We also show the SEDs obtained for two different values of the parameter . We compare our results with observations of Cygnus X-1 made by COMPTEL (McConnell et al., 2000), obtaining good agreement. The best-fit model is obtained with for a corona dominated by protons () and for ; both values are significantly lower that the maximum energy available for particle acceleration.
Jourdain et al. (2012) reported observations of Cygnus X-1 made with the SPI instrument onboard the INTEGRAL satellite. Despite no significant emission above 1 MeV being detected by INTEGRAL, the upper limits are consistent with the non-thermal tail observed at several MeV by COMPTEL, and in close agreement with our model.
The gap observed in the energy range keV is produced by the internal absorption in the corona and accretion disk fields. To quantify this effect, we show in Fig. 7 the opacity as a function of the photon energy at different depths inside the corona. Given the high values of the opacity, the emission is completely suppressed. This result is in accordance with the non-detection of Cygnus X-1 in a steady state by Fermi. As pointed out in Romero et al. (2010c), all emission detected in this energy range should be produced in the jet (e.g., Bosch-Ramon et al. 2008).
The absorption in the stellar field partially suppresses the high-energy bump at eV, which makes it difficult to detect this source using either the MAGIC or VERITAS Cherenkov telescopes. The high-energy emission may be detectable by future instruments with higher sensitivity and wider energy ranges, such as the Cherenkov Telescope Array (CTA).
4 Transient episode or flare
It is well-known that X-ray binaries undergo transient radiative episodes. In 1999, the BATSE instrument detected an increase in the luminosity above keV in Cygnus X-1, of an order of magnitude (Stern et al., 2001). Moreover, between 1995 and 2003, seven outbursts have been reported at the location of this source with a significance of 3 or more (Mazets et al., 1996; Romero et al., 2002; Golenetskii et al., 2003). The luminosities above 15 keV of the outbursts were in the range erg s, which are much higher than the typical thermal luminosity in the hard state.
More recently, the claimed detection of Cygnus X-1 during a flaring state by Albert et al. (2007) and the suggested detection by the AGILE satellite (Tavani et al., 2009) constitute the first presumed evidence of very high-energy gamma ray emission produced around a Galactic black hole.
Another example is the X-ray binary Cygnus X-3, from which four gamma-ray flares were detected by AGILE satellite (Tavani et al., 2009). Variable emission was indeed detected by Fermi (Abdo et al., 2009). This emission, however, is likely produced in the jet (Bednarek, 2010; Araudo et al., 2010; Cerutti et al., 2011).
With the aim of studying transient events in the corona, we apply our model to the case of a non-thermal flare, i.e. a flare produced by changes in the total power injected into relativistic particles (). The dynamical element during these type of outbursts is the magnetic field; a sudden injection might be the result of both fast and large-scale reconnection events. This assumption is supported by observations of solar flares that suggest that magnetic reconnection can trigger diffusive acceleration without the requirement of strong shock formation (Tsuneta & Naito, 1998; Lin, 2008; Kowal et al., 2011). The overall thermal luminosity can remain constant during these episodes.
4.1 Particle injection
Although the light curves of Galactic black holes during outbursts or transient episodes can vary from source to source, there are some common features. The rise time tends to be much shorter than the decay time (Grove et al., 1998), thus the light curves are usually called FRED (Fast Rise and Exponential Decay). A simple analytic expression that can represent this behavior is given by (Romero et al., 2010b)
where , , and are the rise, decay, and plateau duration, respectively. We adopt min, h, and h for a rapid flare. The power-law has the standard index of . The normalization constant can be obtained from the total power injected into relativistic protons and electrons, . This power is assumed to be a fraction of the luminosity of the corona . As mentioned in the previous section, in the steady state the best fit to the observations is obtained with for . During the flare, the number of relativistic particles increases. In our model, the power injected into the flare doubles that injected in the steady state, but larger flares are quite possible, as observed in the Sun (Lin, 2008). It is assumed that the thermal corona remains unaffected during the event.
4.2 Spectral energy distribution
In Galactic black hole coronae, cooling timescales are significantly shorter than the flare timescales, which are typically of hours or even days (Malzac & Jourdain, 2000). Consequently, the transport equation could be equally solved assuming a steady state, and considering changes in the luminosity as the flare evolves. Since one of the aims of this work had been to develop a code that can be applied to different environments, we include the time-dependent term in the kinetic equations. The transport equations then have the following form (Ginzburg & Syrovatskii, 1964)
where as before
We solve the set of coupled equations using the treatment described in Sec. 3.3, but now including the time dependence so the computing time significantly increases.
Figure 8 shows the evolution of the electromagnetic emission during a day. In this figure, we do not include the absorption in the photon field of the star, since it depends on the orbital phase where the flare is produced. For this purpose, in Fig. 9 we show the absorption coefficient related to the stellar field. We estimate the opacity for flares occurring at different orbital phases, and we conclude that for some values of the absorption of the star is almost negligible. For example, flares at energies above GeV can be detected at phases , with instruments such as CTA.
4.3 Neutrino emission
The electromagnetic flare may be absorbed above MeV. This is not true for a neutrino burst. If a sudden injection of relativistic protons occured, the neutrino flux produced could be detectable by instruments such as IceCube.
We consider neutrino production by two main channels: charged pion decay
and muon decay
with , and
In this latter expression, , , and . Our calculations take into account pion and muon losses, as in Reynoso & Romero (2009).
The differential flux of neutrinos arriving at the Earth can be obtained as
This quantity, weighted by the squared energy, is shown in Fig. 10, which also shows the IceCube sensitivity for one year of operation. Assuming that the duty cycle of flares in Galactic black holes is around %, the IceCube detector will be able to detect neutrinos from a source at kpc (Reid et al., 2011) after ten years of observations.
Nonetheless, a variability search with carefully binned time spans might yield positive results long before (Vieyro et al. 2012, in preparation), avoiding the disadvantage of the smoothing cause by the averaging of the observations.
We have developed a model to deal self-consistently with the non-thermal emission from a magnetized corona. The assumption of equipartition among the magnetic energy density, the bolometric photon energy density of the corona, and the kinetic energy density of the plasma is used to estimate the value of the most relevant parameters. This is supported by the X-ray binaries spending a significant amount of time in the low-hard state, permitting in turn a significant amount of time for field-particle interactions. Under this assumption, the model presented in this work is self-similar, because if a value of one parameter is changed, the other parameters can be re-scaled, producing no significant differences to the results and predictions.
On the other hand, we found that if we consider a sub-equipartition magnetic field, the radiative output can be very different. The role of the magnetic field in the cooling of the different types of relativistic particles in the corona is very important. If lower values of the magnetic field are considered, the main radiative losses will be produced by the interaction with matter (e.g., relativistic bremsstrahlung), hence the spectra of relativistic particles will be modified.
The hadronic content of the plasma, given by the parameter , is unknown. We have adopted two values: , which represents a proton-dominated corona, and , i.e. a corona with equal contributions from both protons and electrons. Both models are capable of reproducing the non-thermal spectrum observed by COMPTEL and INTEGRAL in Cygnus X-1. The main difference is caused by the energy injected into relativistic particles: the lower the hadronic content, the higher the power injected in relativistic particles. Neutrino production increases with .
A significant amount of flux in the range 10 MeV - 1 TeV is not expected because of absorption in the thermal photon field. We instead predict the existence of a bump at very high energies ( TeV). This high-energy emission may be detectable in the future from different sources. Since it is of hadronic origin, detections or upper limits can be used to place constraints on the number of relativistic protons in the corona. Orbital modulation can be important owing to the variable absorption (Romero et al., 2010a).
For sources where it is difficult to detect an electromagnetic flare, the neutrino production can nevertheless yield detectable events. Our results show that there may be instances in which a neutrino flare could be detected, but the gamma-ray counterpart is not.
Petropoulou & Mastichiadis (2011, 2012) studied the temporal signatures of electron and proton injection in a compact, magnetized source. Both photon quenching and synchrotron emission are responsible for the dynamical non-linear behavior of the system. Although these results cannot be directly extrapolated to our corona model, since interactions and background thermal photon fields cannot be ignored, future research may identify evidence of some kind of dynamical cycle in sources such as Cygnus X-1.
We have illustrated that a consistent treatment of the non-thermal emission from a magnetized corona can be implemented by solving the set of coupled differential equations for all particle species.
Our application to Cygnus X-1 provided both a good fit to the observational data and interesting predictions for very high energy and neutrino instruments.
In the future, we will explore flare episodes taking place in low-mass X-ray binaries. These are attractive objects for the application of our model, because in these systems the companion star has a weak radiation field. Since the absorption will be negligible, both the electromagnetic and the neutrino flares may be detectable.
We thank María Victoria del Valle for help in several aspects of this work. F.L.V. also thanks Abel Valente and Gabriela Vila for their help in performing the numerical simulations and optimization of the computational implementations. This work was partially supported by the Argentine Agencies CONICET (PIP 0078) and ANPCyT (PICT 2007-00848), as well as by Spanish Ministerio de Ciencia e Innovación (MICINN) under grant AYA2010-21782-C03-01.
- Abdo et al. (2009) Abdo, A. A., Ackermann, M., Ajello, M., et al. 2009, Science, 326, 1512
- Abramowicz et al. (1995) Abramowicz, M. A., Chen, X., Kato, S., Lasota, J.-P., & Regev, O. 1995, ApJ, 438, L37
- Aharonian et al. (1983) Aharonian, F. A., Atoyan, A. M., & Nagapetian, A. M. 1983, Astrofizika, 19, 323
- Aharonian et al. (1985) Aharonian, F. A., Kririllov-Ugrimov, V. G., & Vardanian, V. V. 1985, Ap&SS, 115, 187
- Aharonian & Plyasheshnikov (2003) Aharonian, F. A. & Plyasheshnikov, A. V. 2003, Astroparticle Physics, 19, 525
- Albert et al. (2007) Albert, J., Aliu, E., Anderhub, H., et al. 2007, ApJ, 665, L51
- Araudo et al. (2010) Araudo, A., Bosch-Ramon, V., & Romero, G. E. 2010, PoS tsra2010, 184, 1
- Bednarek (1993) Bednarek, W. 1993, A&A, 278, 307
- Bednarek (2000) Bednarek, W. 2000, A&A, 362, 646
- Bednarek (2010) Bednarek, W. 2010, MNRAS, 406, 689
- Bednarek & Giovannelli (2007) Bednarek, W. & Giovannelli, F. 2007, A&A, 464, 437
- Belmont et al. (2008) Belmont, R., Malzac, J., & Marcowith, A. 2008, A&A, 491, 617
- Beloborodov (1999) Beloborodov, A. M. 1999, in Astronomical Society of the Pacific Conference Series, Vol. 161, High Energy Processes in Accreting Black Holes, ed. J. Poutanen & R. Svensson, 295
- Bosch-Ramon et al. (2008) Bosch-Ramon, V., Khangulyan, D., & Aharonian, F. A. 2008, A&A, 489, L21
- Böttcher & Dermer (2005) Böttcher, M. & Dermer, C. D. 2005, ApJ, 634, L81
- Böttcher & Schlickeiser (1996) Böttcher, M. & Schlickeiser, R. 1996, A&AS, 120, C575
- Böttcher & Schlickeiser (1997) Böttcher, M. & Schlickeiser, R. 1997, A&A, 325, 866
- Cadolle Bel et al. (2006) Cadolle Bel, M., Sizun, P., Goldwurm, A., et al. 2006, A&A, 446, 591
- Carraminana (1992) Carraminana, A. 1992, A&A, 264, 127
- Cerutti et al. (2011) Cerutti, B., Dubus, G., Malzac, J., et al. 2011, A&A, 529, A120
- Chodorowski et al. (1992) Chodorowski, M. J., Zdziarski, A. A., & Sikora, M. 1992, ApJ, 400, 181
- Coppi (1992) Coppi, P. S. 1992, MNRAS, 258, 657
- Coppi & Blandford (1990) Coppi, P. S. & Blandford, R. D. 1990, MNRAS, 245, 453
- de Gouveia dal Pino & Lazarian (2005) de Gouveia dal Pino, E. M. & Lazarian, A. 2005, A&A, 441, 845
- de Gouveia Dal Pino et al. (2010) de Gouveia Dal Pino, E. M., Piovezan, P. P., & Kadowaki, L. H. S. 2010, A&A, 518, A5
- del Valle et al. (2011) del Valle, M. V., Romero, G. E., Luque-Escamilla, P. L., Martí, J., & Ramón Sánchez-Sutil, J. 2011, ApJ, 738, 115
- Dove et al. (1997) Dove, J. B., Wilms, J., Maisack, M., & Begelman, M. C. 1997, ApJ, 487, 759
- Drury (1983) Drury, L. O. 1983, Reports on Progress in Physics, 46, 973
- Drury (2012) Drury, L. O. 2012, [arXiv1201.6612O]
- Dubus (2006) Dubus, G. 2006, A&A, 451, 9
- Esin et al. (1997) Esin, A. A., McClintock, J. E., & Narayan, R. 1997, ApJ, 489, 865
- Galeev et al. (1979) Galeev, A. A., Rosner, R., & Vaiana, G. S. 1979, ApJ, 229, 318
- George & Fabian (1991) George, I. M. & Fabian, A. C. 1991, MNRAS, 249, 352
- Gierlinski et al. (1997) Gierlinski, M., Zdziarski, A. A., Done, C., et al. 1997, MNRAS, 288, 958
- Ginzburg & Syrovatskii (1964) Ginzburg, V. L. & Syrovatskii, S. I. 1964, The Origin of Cosmic Rays, ed. Macmillan, New York
- Golenetskii et al. (2003) Golenetskii, S., Aptekar, R., Frederiks, D., et al. 2003, ApJ, 596, 1113
- Grove et al. (1998) Grove, J. E., Johnson, W. N., Kroeger, R. A., et al. 1998, ApJ, 500, 899
- Haardt & Maraschi (1993) Haardt, F. & Maraschi, L. 1993, ApJ, 413, 507
- Herterich (1974) Herterich, K. 1974, Nature, 250, 311
- Ichimaru (1977) Ichimaru, S. 1977, ApJ, 214, 840
- Jourdain et al. (2012) Jourdain, E., Roques, J. P., & Malzac, J. 2012, ApJ, 744, 64
- Kowal et al. (2011) Kowal, G., de Gouveia Dal Pino, E. M., & Lazarian, A. 2011, ApJ, 735, 102
- Lazarian & Vishniac (1999) Lazarian, A. & Vishniac, E. T. 1999, ApJ, 517, 700
- Li et al. (1996) Li, H., Kusunose, M., & Liang, E. P. 1996, ApJ, 460, L29
- Li & Miller (1997) Li, H. & Miller, J. A. 1997, ApJ, 478, L67
- Lightman & Zdziarski (1987) Lightman, A. P. & Zdziarski, A. A. 1987, ApJ, 319, 643
- Lin (2008) Lin, R. P. 2008, in American Institute of Physics Conference Series, Vol. 1039, American Institute of Physics Conference Series, ed. G. Li, Q. Hu, O. Verkhoglyadova, G. P. Zank, R. P. Lin, & J. Luhmann , 52–62
- Lipari et al. (2007) Lipari, P., Lusignoli, M., & Meloni, D. 2007, Phys. Rev. D, 75, 123005
- Malzac & Belmont (2009) Malzac, J. & Belmont, R. 2009, MNRAS, 392, 570
- Malzac & Jourdain (2000) Malzac, J. & Jourdain, E. 2000, A&A, 359, 843
- Markoff et al. (2001) Markoff, S., Falcke, H., & Fender, R. 2001, A&A, 372, L25
- Mazets et al. (1996) Mazets, E. P., Aptekar, R. L., Frederiks, D. D., et al. 1996, in American Institute of Physics Conference Series, Vol. 384, American Institute of Physics Conference Series, ed. C. Kouveliotou, M. F. Briggs, & G. J. Fishman, 492–496
- McConnell et al. (2000) McConnell, M. L., Ryan, J. M., Collmar, W., et al. 2000, ApJ, 543, 928
- Narayan & Yi (1994) Narayan, R. & Yi, I. 1994, ApJ, 428, L13
- Narayan & Yi (1995a) Narayan, R. & Yi, I. 1995a, ApJ, 444, 231
- Narayan & Yi (1995b) Narayan, R. & Yi, I. 1995b, ApJ, 452, 710
- Nayakshin (1999) Nayakshin, S. 1999, in Astronomical Society of the Pacific Conference Series, Vol. 162, Quasars and Cosmology, ed. G. Ferland & J. Baldwin, 87
- Novikov & Thorne (1973) Novikov, I. D. & Thorne, K. S. 1973, in Black Holes (Les Astres Occlus), ed. C. Dewitt & B. S. Dewitt, 343–450
- Orosz et al. (2011) Orosz, J. A., McClintock, J. E., Aufdenberg, J. P., et al. 2011, ApJ, 742, 84
- Pellizza et al. (2010) Pellizza, L. J., Orellana, M., & Romero, G. E. 2010, Int. J. Mod. Phys. D, 19, 671
- Petropoulou & Mastichiadis (2011) Petropoulou, M. & Mastichiadis, A. 2011, A&A, 532, A11
- Petropoulou & Mastichiadis (2012) Petropoulou, M. & Mastichiadis, A. 2012, MNRAS, [arXiv1201.2091]
- Pilla & Shaham (1997) Pilla, R. P. & Shaham, J. 1997, ApJ, 486, 903
- Poutanen (1998) Poutanen, J. 1998, in Theory of Black Hole Accretion Disks, ed. M. A. Abramowicz, G. Bjornsson, & J. E. Pringle, Cambridge Contemporary Astrophysics (Cambridge University Press, Cambridge, UK), 100
- Poutanen et al. (1997) Poutanen, J., Krolik, J. H., & Ryde, F. 1997, 292, L21
- Poutanen & Vurm (2009) Poutanen, J. & Vurm, I. 2009, ApJ, 690, L97
- Press et al. (1992) Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical recipes in FORTRAN. The art of scientific computing, ed. Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P.
- Reid et al. (2011) Reid, M. J., McClintock, J. E., Narayan, R., et al. 2011, ApJ, 742, 83
- Reynoso & Romero (2009) Reynoso, M. M. & Romero, G. E. 2009, 493, 1
- Romero et al. (2010a) Romero, G. E., Del Valle, M. V., & Orellana, M. 2010a, 518, A12
- Romero et al. (2002) Romero, G. E., Kaufman Bernadó, M. M., & Mirabel, I. F. 2002, 393, L61
- Romero & Paredes (2011) Romero, G. E. & Paredes, J. M. 2011, Introducción a la Astrofísica Relativista, 1st edn. (Publicacions i Edicions de la Universitat de Barcelona, Barcelona)
- Romero et al. (2010b) Romero, G. E., Reynoso, M. M., & Christiansen, H. R. 2010b, A&A, 524, A4
- Romero et al. (2010c) Romero, G. E., Vieyro, F. L., & Vila, G. S. 2010c, 519, A109
- Schopper et al. (1998) Schopper, R., Lesch, H., & Birk, G. T. 1998, A&A, 335, 26
- Shakura & Sunyaev (1973) Shakura, N. I. & Sunyaev, R. A. 1973, 24, 337
- Shapiro et al. (1976) Shapiro, S. L., Lightman, A. P., & Eardley, D. M. 1976, 204, 187
- Stern et al. (2001) Stern, B. E., Beloborodov, A. M., & Poutanen, J. 2001, ApJ, 555, 829
- Stern et al. (1995) Stern, B. E., Poutanen, J., Svensson, R., Sikora, M., & Begelman, M. C. 1995, ApJ, 449, L13
- Stirling et al. (2001) Stirling, A. M., Spencer, R. E., de la Force, C. J., et al. 2001, MNRAS, 327, 1273
- Svensson (1982) Svensson, R. 1982, ApJ, 258, 335
- Tavani et al. (2009) Tavani, M., Bulgarelli, A., Piano, G., et al. 2009, Nature, 462, 620
- Tsuneta & Naito (1998) Tsuneta, S. & Naito, T. 1998, ApJ, 495, L67
- Veledina et al. (2011) Veledina, A., Vurm, I., & Poutanen, J. 2011, MNRAS, 414, 3330
- Vieyro et al. (2012) Vieyro, F. L., Sestayo, Y., Romero, G. E., & Paredes, J. M. 2012, in preparation
- Vila & Aharonian (2009) Vila, G. S. & Aharonian, F. A. 2009, in Compact Objects and their Emission, ed. G. E. Romero & P. Benaglia, Asociación Argentina de Astronomía, Book Series (Editorial Paideia, La Plata, Argentina), 1–38
- Vila & Romero (2010) Vila, G. S. & Romero, G. E. 2010, MNRAS, 403, 1457
- Vila et al. (2012) Vila, G. S., Romero, G. E., & Casco, N. A. 2012, A&A, 538, A97
- Vurm & Poutanen (2009) Vurm, I. & Poutanen, J. 2009, ApJ, 698, 293
- White et al. (1988) White, T. R., Lightman, A. P., & Zdziarski, A. A. 1988, ApJ, 331, 939
- Zdziarski (1998) Zdziarski, A. A. 1998, MNRAS, 296, L51
- Zdziarski et al. (2009) Zdziarski, A. A., Malzac, J., & Bednarek, W. 2009, MNRAS, 394, L41
- Zenitani & Hoshino (2001) Zenitani, S. & Hoshino, M. 2001, ApJ, 562, L63