Free Access
Issue
A&A
Volume 507, Number 3, December I 2009
Page(s) 1217 - 1224
Section Astrophysical processes
DOI https://doi.org/10.1051/0004-6361/200912494
Published online 01 October 2009

A&A 507, 1217-1224 (2009)

One-dimensional pair cascade emission in gamma-ray binaries

An upper-limit to cascade emission at superior conjunction in LS 5039

B. Cerutti - G. Dubus - G. Henri

Laboratoire d'Astrophysique de Grenoble, UMR 5571 CNRS, Université Joseph Fourier, BP 53, 38041 Grenoble, France

Received 14 May 2009 / Accepted 17 September 2009

Abstract
Context. In gamma-ray binaries such as LS 5039, a large number of electron-positron pairs are created by the annihilation of primary very high-energy (VHE) gamma rays with photons from the massive star. The radiation from these particles contributes to the total high-energy gamma-ray flux and can initiate a cascade, decreasing the effective gamma-ray opacity in the system.
Aims. The aim of this paper is to model the cascade emission and investigate whether it can account for the VHE gamma-ray flux detected by HESS from LS 5039 at superior conjunction, where the primary gamma rays are expected to be fully absorbed.
Methods. A one-dimensional cascade develops along the line-of-sight if the deflections of pairs induced by the surrounding magnetic field can be neglected. A semi-analytical approach can then be adopted, including the effects of the anisotropic seed radiation field from the companion star.
Results. Cascade equations are numerically solved, yielding the density of pairs and photons. In LS 5039, the cascade contribution to the total flux is large and anti-correlated with the orbital modulation of the primary VHE gamma rays. The cascade emission dominates close to superior conjunction but is too strong to be compatible with HESS measurements. Positron annihilation does not produce detectable 511 keV emission.
Conclusions. This study provides an upper limit to cascade emission in gamma-ray binaries at orbital phases where absorption is strong. The pairs are likely to be deflected or isotropized by the ambient magnetic field, which will reduce the resulting emission seen by the observer. Cascade emission remains a viable explanation for the detected gamma rays at superior conjunction in LS 5039.

Key words: radiation mechanisms: non-thermal - stars: individual: LS 5039 - gamma rays: theory - X-rays: binaries

1 Introduction

The massive star in gamma-ray binaries plays a key role in the formation of very high-energy (VHE, >100 GeV) radiation. The large seed-photon density provided by the O or Be companion star, contributes to the production of gamma rays via inverse Compton scattering on ultra-relativistic electrons accelerated in the system (e.g. in a pulsar wind or a jet). The same photons annihilate with gamma rays, leading to electron-positron pairs production $\gamma+\gamma\rightarrow {\rm e}^{+}+{\rm e}^{-}$. In some tight binaries such as LS 5039, this gamma-ray absorption mechanism is strong if the VHE emission occurs close to the compact object. Gamma-ray absorption can account for an orbital modulation in the VHE gamma-ray flux from LS 5039, as observed by HESS (Dubus 2006; Bednarek 2006; Böttcher & Dermer 2005).

A copious number of pairs may be produced in the surrounding medium as a by-product of the VHE gamma-ray absorption. If the number of pairs created is large enough and if they have enough time to radiate VHE photons before escaping, a sizeable electromagnetic cascade can be initiated. New generations of pairs and gamma rays are produced as long as the secondary particles have enough energy to boost stellar photons beyond the pair production threshold energy. Because of the anisotropic stellar photon field in the system, the inverse Compton radiation produced in the cascade has a strong angular dependence. The cascade contribution depends on the position of the primary gamma-ray source with respect to the massive star and a distant observer.

The VHE modulation in LS 5039 was explained in Dubus et al. (2008) using phase-dependent absorption and inverse Compton emission, ignoring the effect of pair cascading. This model did not predict any flux close to superior conjunction, i.e. where the massive star lies between the compact object and the observer. This is contradicted by HESS observations (Aharonian et al. 2006a). Interestingly, this mismatch intervenes at phases where $\gamma\gamma$-opacity is known to be high $\tau_{\gamma\gamma}\gg 1$. The development of a cascade could contribute to the residual flux observed in the system, with secondary gamma-ray emission filling in for the highly absorbed primary gamma rays. This possibility has been proposed to explain this discrepancy (Aharonian et al. 2006a) and is quantitatively investigated in this article.

The ambient magnetic field strength has a critical impact on the development of pair cascading. If the magnetic field strength is low enough to neglect the induced deflections on pair trajectories then the cascade develops along the line of sight joining the primary source of gamma rays and a distant observer. The particles do not radiate synchrotron radiation. Cascade calculations are then reduced to a one-dimension problem. Such a situation would apply in an unshocked pulsar wind where the pairs are cold relative to the magnetic field carried in the wind. This paper explores the development of an one-dimensional pair cascade in a binary and its implications.

Previous computations of cascade emission in binary environment were carried out by Bednarek (2007); Zdziarski et al. (2009); Bednarek (2006,1997); Sierpowska-Bartosik & Torres (2008); Aharonian et al. (2006b); Orellana et al. (2007); Sierpowska & Bednarek (2005); Khangulyan et al. (2008). Except for Aharonian et al. (2006b), all these works are based on Monte Carlo methods. One peculiarity of the gamma-ray binary environment is that the source of seed-photons for pair production and inverse Compton emission is the high luminosity companion star. This study proposes a semi-analytical model for one-dimensional cascades calculations, taking into account the anisotropy in the seed-photon field. The aim of the paper is to investigate and compute the total contribution from pair cascading in the system LS 5039, and see if it can account for the measured flux close to superior conjunction. The next section presents the main assumptions and equations for cascade computations. The development and the anisotropic effects of pair cascading in compact binaries are investigated. The density of escaping pairs and their rate of annihilation are also calculated in this part. The cascade contribution along the orbit in LS 5039 is computed and compared with the available observations in Sect. 3. The last section concludes on the implications of one-dimensional cascades in gamma-ray binaries. More details about pair production are available in the appendices.

2 Anisotropic pair cascading in compact binaries

2.1 Assumptions

This part examines one-dimensional cascading in the context of binary systems. The massive star sets the seed-photon radiation field for the cascade. For simplicity, the massive star is assumed point-like and mono-energetic. This is a reasonable approximation as previous studies on absorption (Dubus 2006) and emission (Dubus et al. 2008) have shown. The effects of the magnetic field and pair annihilation are neglected (see Sect. 2.5). Triplet pair production (TPP) due to the high-energy electrons or positrons propagating in a soft photon field ( $\gamma+{\rm e}^{+-} \rightarrow {\rm e}^{+-} + {\rm e}^+ +{\rm e}^-$, Mastichiadis 1991) is not taken into account here. The cross section for this process becomes comparable to inverse Compton scattering when $E_{\rm e} \epsilon_0\ga 250 (m_{\rm e} c^2)^2$ that is for electron energies $E_{\rm e}\ga 6$ TeV interacting with $\epsilon_0\approx 10$ eV stellar photons. With a scattering rate of about $\sim$10-2  $\rm {s^{-1}}$, only a few pairs can be created via TPP by each VHE electron, before it escapes or loses its energy in a Compton scattering. The created pairs have much lower energy than the primary electrons. TPP cooling remains inefficient compared to inverse Compton for VHE electrons with energy $\la$PeV. HESS observations of LS 5039 show a break in the spectrum at a few TeV so few electrons are expected to interact by TPP in the cascade. Observations of other gamma-ray binaries also show steep spectra but this assumption will have to be revised if there is significant primary emission beyond $\approx$10 TeV. Pair production due to high-energy gamma rays interacting with the surrounding material is also neglected. This occurs for $\gamma$-rays >1 MeV and the cross-section is of order $0.04 \sigma_{\rm T} Z^2 ~\rm {cm}^{2}$ (see e.g. Longair 1992), with $\sigma_{\rm T}$ the Thomson cross-section. Since the measured $N_{\rm H}$ is at most 1022 cm-2 in gamma-ray binaries, pair production on matter will not affect the propagation of gamma rays towards the observer.

Due to the high velocity of the center-of-mass (CM) frame in the observer frame, the direction of propagation of pairs created by $\gamma\gamma$-absorption is boosted in the direction of the initial gamma ray. For a gamma ray of energy $\epsilon_1=1~$TeV, the Lorentz factor of the CM to the observer frame transform is $\gamma'\sim \epsilon_1/2 m_{\rm e} c^2=10^6\gg 1$ (see the Appendix, Eq. (A.2)). Pairs produced in the cascade are ultra-relativistic with typical Lorentz factor $\gamma_{\rm e}\sim 10^6\gg 1$. Their emission is forward boosted within a cone of semi-aperture angle $\alpha\sim 1/\gamma_{\rm e}\ll 1$ in the direction of electrons. The deviations on the electron trajectory due to scattering in the Thomson regime are $\sim$ $\epsilon_0/m_{\rm e} c^2 \ll 1/\gamma_{\rm e}$. In the Klein-Nishina regime most of the electron energy is given to the photon. It is assumed here that electrons and photons produced in the cascade remain on the same line, a good approximation since $\gamma'$ and $\gamma_{\rm e}\gg 1$. This line joins the primary gamma-ray source to a distant observer (Fig. 1).

Pair cascading is one-dimensional as long as magnetic deviations of pairs trajectories along the Compton interaction length  $\lambda_{\rm ic}$ remain within the cone of emission of the electrons. This condition holds if $\lambda_{\rm ic}/(2R_{\rm L})<1/\gamma_{\rm e}$, with $R_{\rm L}$ the Larmor radius. For a typical interaction length $\lambda_{\rm ic}\sim 1/(n_{\star}\sigma_{\rm ic})\sim 10^{11}~$cm for TeV pairs in LS 5039, the ambient magnetic field must be lower than $B\la 10^{-8}~$G. If the magnetic field strength is much greater, pairs locally isotropize and radiate in all directions. In between, pairs follow the magnetic field lines and the dynamics of each pairs must be followed as treated in Sierpowska & Bednarek (2005). The above limit may appear unrealistically stringent. However, since deviations and isotropization will dilute the cascade flux, the one-dimensional approach can be seen as maximizing the cascade emission. More exactly, this redistribution induced by magnetic deflections would decrease the cascade flux at orbital phases where many pairs are produced to the benefit of phases where only a few are created. Hence, the one-dimensional approach gives an upper limit to the cascade contribution at phases where absorption is strong. If the flux calculated here using this assumption is lower than required by observations then cascading will be unlikely to play a role. Finally, one-dimensional cascading should hold in the free pulsar wind as long as the pairs move strictly along the magnetic field. In Sierpowska & Bednarek (2005) and Sierpowska-Bartosik & Torres (2008), the cascade radiation is computed up to the termination shock using a Monte Carlo approach. Sierpowska & Bednarek (2005) also include a contribution from the region beyond the shock. The cascade electrons in this region are assumed to follow the magnetic field lines (in contrast with the pulsar wind zone where the propagation is radial). There is no reacceleration at the shock and synchrotron losses are neglected. In the method expounded here, the cascade radiation is calculated semi-analytically from a point-like gamma-ray source at the compact object location up to infinity, providing the maximum possible contribution of the one-dimensional cascade in gamma-ray binaries.

\begin{figure}
\par\includegraphics[width=8.5cm,clip]{12494fg1.eps}
\end{figure} Figure 1:

This diagram describes the system geometry. A gamma-ray photon of energy $\epsilon _1$ from the primary source (compact object) interacts with a soft photon of energy $\epsilon _0$ at a distance r from the source and R from the massive star (assumed point-like and mono-energetic), producing a pair $\rm e^+/e^-$ boosted toward a distant observer. The system is seen at an angle $\psi $.

Open with DEXTER

2.2 Cascade equations

In order to compute the contribution from the cascade, the radiative transfer equation and the kinetic equation of the pairs have to be solved simultaneously.

The radiative transfer equation for the gamma-ray density $n_{\gamma}\equiv {\rm d}N_{\gamma}/{\rm d}t {\rm d}\epsilon_1 {\rm d}\Omega$ at a distance r from the source is

\begin{displaymath}\frac{{\rm d}n_{\gamma}}{{\rm d}r}=-n_{\gamma}\left(\frac{{\r...
... d}N}{{\rm d}t{\rm d}\epsilon_1}\ n_{\rm e}\ {\rm d}E_{\rm e},
\end{displaymath} (1)

where $n_{\rm e}\equiv {\rm d}N_{\rm e}/{\rm d}r {\rm d}E_{\rm e} {\rm d}\Omega_{\rm e}$ is the electrons distribution, $n_{\star}$ the seed-photon density from the massive star and ${\rm d}N/{\rm d}t{\rm d}\epsilon_1$ the Compton kernel. The kernel is normalized to the soft photon density and depends on the energy $E_{\rm e}$ of the electron and the angle between the photon and the direction of motion of the electron (Dubus et al. 2008). In the mono-energetic and point-like star approximation the stellar photon density can be estimated as $L_{\star}/4\pi c R^2 \bar{\epsilon_0}$, where $L_{\star}$ is the stellar luminosity, $\bar{\epsilon_0}\approx 2.7 k T_{\star}$ the mean thermal photon energy and R the distance to the massive star (see Fig. 1). The absorption rate ${\rm d}\tau_{\gamma\gamma}/{\rm d}r$ is given by Eq. (B.8), convoluted to the soft photon density.

The kinetic equation for the pairs is given by the following integro-differential equation for $\gamma_{\rm e}\gg 1$ (D'Avezac et al. 2007; Blumenthal & Gould 1970; Zdziarski 1988)

                                      $\displaystyle \frac{{\rm d}n_{\rm e}}{{\rm d}t}$ = $\displaystyle -n_{\rm e}(E_{\rm e})\int_{m_{\rm e} c^2}^{E_{\rm e}}\mathcal{P}\left(E_{\rm e},E'_{\rm e}\right){\rm d}E'_{\rm e}$  
    $\displaystyle + \int_{E_{\rm e}}^{+\infty}n_{\rm e}(E'_{\rm e})\ \mathcal{P}\le...
...}E'_{\rm e} + 2\int n_{\star}\ g_{\gamma\gamma}\ n_{\gamma}\ {\rm d}\epsilon_1,$ (2)

where $\mathcal{P}(E_{\rm e},E'_{\rm e})$ is the transition rate for an electron of energy $E_{\rm e}$ down-scattered at an energy $E'_{\rm e}\le E_{\rm e}$ at r. The first two terms on the right side of the equation describe the inverse Compton cooling of pairs, taking into account catastrophic losses in the deep Klein-Nishina regime. In this case, most of the electron energy is lost in the interaction and the scattered photon carries away most of its energy since $\epsilon_1=E_{\rm e}-E'_{\rm e}\approx E_{\rm e}$. A continuous losses equation inadequately describes sizeable stochastic losses in a single interaction (Zdziarski 1989; Blumenthal & Gould 1970).

Since the inverse Compton kernel gives the probability per electron of energy $E_{\rm e}$ to produce a gamma ray of energy $\epsilon _1$, the scattering rate can be rewritten as

\begin{displaymath}\mathcal{P}(E_{\rm e},E_{\rm e}')=n_{\star}(r)\frac{{\rm d}N}{{\rm d}t{\rm d}E'_{\rm e}}\cdot
\end{displaymath} (3)

The expression of ${\rm d}N/{\rm d}t{\rm d}E'_{\rm e}$ is the same as the Compton kernel as described before but gives the spectrum of scattered electrons instead of the outcoming photon. The first integral in Eq. (2) is the inverse Compton scattering rate and can be analytically expressed as

\begin{displaymath}\int_{m_{\rm e} c^2}^{E_{\rm e}}\mathcal{P}\left(E_{\rm e},E'...
... ic}\ c\ n_{\star}(r)\left(1-\beta_{\rm e}\cos\theta_0\right),
\end{displaymath} (4)

where $\beta_{\rm e}$ is the electron velocity in the observer frame and $\sigma_{\rm ic}$ is the total inverse Compton cross-section (for the full expression see e.g. Rybicki & Lightman 1979, Eq. (7.5)). The last term in the kinetic equation is a source of pairs from $\gamma\gamma$-absorption coupled with the photon density (see the appendices). The pair production kernel $g_{\gamma \gamma }$ is normalized to the soft photon density.

The anisotropic cascade can be computed by inserting the anisotropic kernels for inverse Compton scattering (see Eq. (A.6) in Dubus et al. 2008) and for pair production obtained in Eq. (B.5) in Eqs. (1), (2). The following sections present cascade calculations applied to compact binaries, using a simple Runge-Kutta 4 integration method. It is more convenient to perform integrations over an angular variable rather than r. Here, calculations are carried out using $\psi_r$, the angle between the line joining the massive star to the observation point and the line of sight (see Fig. 1).

2.3 Cascade growth along the line of sight

Figure 2 presents cascade calculations for different distances r from the primary gamma-ray source. For illustrative purpose, the source is assumed isotropic and point-like, injecting a power-law distribution of photons with an index -2 at r=0 but no electrons. The calculations were carried out for a system like LS 5039 and for a viewing angle $\psi =30\hbox {$^\circ $ }$. In this geometric configuration, absorption is known to be strong ( $\tau_{\gamma\gamma}\approx 40$ for 200 GeV photons) and a significant fraction of the total absorbed energy is expected to be reprocessed in the cascade, inverse Compton scattering being also very efficient in this configuration.

Close to the source ($r\la d$ with d the orbital separation), absorption produces a sharp and deep dip in the spectrum (light dashed line) but the cascade starts to fill the gap (black solid line). The angle $\psi_r$ increases with the distance r to the primary source. Hence, the threshold energy for pair production increases as well. Cascading adds more flux to higher energy gamma rays where absorption is maximum. The cascade produces an excess of low energy gamma rays below the minimum threshold energy $\epsilon_1\approx 30$ GeV. Because these new photons do not suffer from absorption, they accumulate at lower energies. This is a well-known feature of cascading.

2.4 Anisotropic effects

\begin{figure}
\par\includegraphics[width=8cm,clip]{12494fg2.eps}
\end{figure} Figure 2:

Cascade development along the path to the observer. The primary source of photons, situated at the location of the compact object, has a power law spectral distribution with photon index -2 (dotted line). Spectra are computed using the parameters appropriate for LS 5039 at superior conjunction ( $d\approx 2 R_{\star }$, $R_{\star}=9.3~R_{\odot}$, $T_{\star}=39~000$ K) for $\psi =30\hbox {$^\circ $ }$. The transmitted spectrum, including cascade emission, is shown at various distances from the primary source: $r=R_{\star }/4$ (black dashed line), $R_{\star }/2$, $R_{\star }$, $2 R_{\star }$ (solid lines) and $r=+\infty $ (dotted-dashed line). Pure absorbed spectra are shown for comparison (light dashed line).

Open with DEXTER
\begin{figure}
\par\includegraphics[width=15.7cm,clip]{12494fg3.eps}
\end{figure} Figure 3:

Spectra as seen by an observer at infinity, taking into account the effect of cascading. Calculations are applied to LS 5039 at periastron for different viewing angle $\psi =30\hbox {$^\circ $ }$, $60\hbox {$^\circ $ }$, $90\hbox {$^\circ $ }$, $120\hbox {$^\circ $ }$ and $150\hbox {$^\circ $ }$. Left panel: complete spectra (solid line) are compared to the pure absorbed (light dashed line) and injected (dotted line) spectra. The contribution from the cascade is presented in the right panel.

Open with DEXTER
This section investigates anisotropic effects in the development of the cascade as seen by a distant observer. Cascades are computed for different viewing angle $\psi $ at infinity, assuming an isotropic power-law spectrum for the primary gamma rays.

The left panel in Fig. 3 shows the complete spectrum taking into account cascading (solid line) compared to the pure absorbed power-law (dashed line). Due to the angular dependence in the pair production process, higher viewing angles shift the cascade contribution to higher energies and decrease its amplitude (Fig. 3, right panel). The cascade flux is low enough to be ignored for $\psi\ga 150\hbox{$^\circ$ }$.

Three different zones can be distinguished in the cascade spectra. First, below the pair production threshold energy, photons accumulate in a low energy tail (photon index $\approx$-1.5) produced by inverse Compton cooling of pairs. For $\psi\la 90\hbox{$^\circ$ }$, a low energy cut-off is observed due to the pairs escaping the system (Cerutti et al. 2008; Ball & Kirk 2000). This low energy cut-off is at about 0.1 GeV for $\psi =30\hbox {$^\circ $ }$. The cutoff occurs when the cascade reaches a distance from the primary source corresponding to $\psi_r\approx 90\hbox{$^\circ$ }$. Then, the electrons cannot cool effectively because the inverse Compton interaction angle diminishes and the stellar photon density decreases as they propagate. For $\psi\ga 90\hbox{$^\circ$ }$, particles escape right away from the vicinity of the companion star and no tail is produced. Second, above the threshold energy, there is a competition between absorption and gamma-ray production by reprocessed pairs, particularly for low angles where both effects are strong. Even if cascading increases the transparency for gamma rays, absorption still creates a dip in the spectrum. Third, well beyond the threshold energy, absorption becomes inefficient. Fewer pairs are created, producing a high-energy cut-off ($\approx$10 TeV, for $\psi =30\hbox {$^\circ $ }$). Klein-Nishina effects also contribute to the decrease of the high-energy gamma-rays production.

2.5 Escaping pairs

The spectrum of pairs produced in the cascade as seen at infinity is shown in Fig. 4. The density depends strongly on the viewing angle as expected, but the mean energy of pairs lies at very high energies ( $\left\langle E_{\rm e}\right\rangle\ga 100$ GeV, see Table 1). The accumulation of very high-energy particles can be explained by two concurrent effects. Far from the massive star ($r\gg d$), most of the pairs are created at very high energy due to the high threshold energy (almost rear-end collision). The second effect is that inverse Compton losses are in deep Klein-Nishina regime for high-energy electrons. The cooling timescale increases and becomes longer than the propagation timescale of electrons close to the companion star, producing an accumulation of pairs at very high energies.

The distribution of pairs allows to assess the fraction of the total absorbed energy escaping the system in the form of kinetic energy in the pairs. This non-radiated power $P_{\rm e}$ can be compared to the radiated power released in the cascade $P_{\rm r}$. Energy conservation yields the total absorbed power $P_{\rm a}=P_{\rm e}+P_{\rm r}$.

\begin{figure}
\par\includegraphics[width=7.7cm,clip]{12494fg4.eps}
\end{figure} Figure 4:

Distribution of escaping pairs seen by a distant observer, depending on the viewing angle $\psi =30\hbox {$^\circ $ }$ (dashed line), $60\hbox {$^\circ $ }$, $90\hbox {$^\circ $ }$, $120\hbox {$^\circ $ }$ and $150\hbox {$^\circ $ }$ (dotted line). The binary parameters are the same as in Fig. 3.

Open with DEXTER

Table 1:   Mean energy of escaping pairs and radiated power efficiency of the cascade.

The asymptotic radiated power reached by the cascade is compared to the total absorbed power integrated over energy in Table 1. The fraction of lost energy increases with the viewing angle. In fact, for $\psi>90\hbox{$^\circ$ }$ most of the power remains in kinetic energy. Once the electrons are created, only a few have time to radiate through inverse Compton interaction. Below ( $\psi<90\hbox{$^\circ$ }$), the radiative power dominates and the cascade is very efficient (recycling efficiency up to 80% for $\psi =30\hbox {$^\circ $ }$). The cascade is fully linear, since the power re-radiated remains much lower than the star luminosity $P_r\ll L_{\star}$ (Svensson 1987). Self-interactions in the cascade are then negligible. This is also a consequence of Klein-Nishina cascading (Zdziarski 1988). In addition, interactions between particles in the cascade would be forcedly rear-end, hence highly inefficient.

The created positrons will annihilate and form a 511 keV line. However, the expected signal is very weak. The annihilation cross-section is $\sigma\sim \sigma_{\rm T} \log \gamma / \gamma$ (see e.g. Longair 1992). The escaping positrons have a very high average Lorentz factor $\gamma \ga 10^5$ (Table 1) so they are unlikely to annihilate within the system. They will thermalize and annihilate in the interstellar medium. Escaping positrons from gamma-ray binaries are unlikely to contribute much to the diffuse 511 keV emission. The average number of pairs created along the orbit in LS 5039 (based on the results to be discussed in the following section) is $\mathcal{N}_{\rm e}\sim 5\times 10^{35}~\rm {s}^{-1}$. This estimate does not take into account contributions from triplet pair production or from the pulsar wind (for a pulsar injecting pairs with $\left\langle\gamma_{\rm e}\right\rangle\sim 10^5$ and a luminosity of $10^{36}~\rm {erg/s}$, about $10^{36}~\rm {s^{-1}}$ pairs are produced). Gamma-ray binaries have short lifetimes and it is unlikely there is more than a few hundred currently active in the Galaxy. Hence, the expected contribution is orders-of-magnitude below the positron flux required to explain the diffuse 511 keV emission ($\sim$ $10^{43}~\rm {s}^{-1}$, Knödlseder et al. 2005). Even if the positrons thermalize close to or within the system (because magnetic fields contain them, see Sect. 5) then, following Guessoum et al. (2006), the expected contribution from a single source at 2 kpc would be at most $\sim$10-9 ph cm-2 s-1, which is currently well below detectability.

3 Cascading in LS 5039

\begin{figure}
\par\includegraphics[width=7.7cm,clip]{12494fg5.eps}
\end{figure} Figure 5:

Orbit-averaged spectra in LS 5039 at INFC ( $0.45<\phi <0.9$, grey lines) and SUPC ($\phi <0.45$ or $\phi >0.9$, black lines) and comparisons with EGRET (dark) and HESS (light) bowties (Aharonian et al. 2006a; Hartman et al. 1999). Dotted-dashed lines represent the primary source of gamma rays with pure absorption, injected at $r\equiv 0$, computed with the model described in Dubus et al. (2008) for a mono-energetic and point-like star. Dashed lines show the contribution from the cascade and thick solid lines the sum of the primary absorbed source and the cascade contributions.

Open with DEXTER

\begin{figure}
\par\includegraphics[width=7.7cm,clip]{12494fg6.eps}
\end{figure} Figure 6:

Computed light-curves along the orbit in LS 5039, in the HESS energy band (flux $\ge $100 GeV). The cascade contribution (dashed line) is compared to the primary pure absorbed source (dotted-dashed line) and HESS observations. The thick solid line shows the sum of both components.

Open with DEXTER

LS 5039 was detected by HESS (Aharonian et al. 2005) and the orbital modulation of the TeV gamma-ray flux was later on reported in Aharonian et al. (2006a). Most of the temporal and spectral features can be understood as a result of anisotropic gamma-ray absorption and emission from relativistic electrons accelerated in the immediate vicinity of the compact object, e.g. in the pulsar wind termination shock (Dubus et al. 2008). However, this description fails to explain the residual flux observed close to superior conjunction where a significant excess has been detected (6.1$\sigma$ at phase $0.0\pm0.05$). The primary gamma rays should be completely attenuated. The aim of this part is to find if cascading can account for this observed flux. The cascade is assumed to develop freely from the primary gamma-ray source up to the observer. The contribution of the cascade as a function of the orbital phase is also investigated.

The primary source of gamma rays now considered is the spectrum calculated in Dubus et al. (2008). Figure 5 shows phase-averaged spectra along the orbit at INFC (orbital phase $0.45<\phi <0.9$) and SUPC ($\phi <0.45$ or $\phi >0.9$) for the primary source, the cascade and the sum of both components. The orbital parameters and the distance (2.5 kpc) are taken from Casares et al. (2005) for an inclination $i=60\hbox{$^\circ$ }$ so $\psi $ varies between $30\hbox{$^\circ$ }{-}150\hbox{$^\circ$ }$. The cascade contribution is highly variable along the orbit and dominates at SUPC for $\epsilon_1\ga 30$ GeV, where a high pair-production rate is expected. At INFC, cascading is negligible compared with the primary flux. With pair cascading the spectral differences between INFC and SUPC are very small at VHE, contrary to what is observed by HESS. In the GeV band, cascades contribute to a spectral hardening at SUPC close to 10-30 GeV.

Orbital light-curves in the HESS energy band give a better appreciation of the contribution from both components (Fig. 6). The contribution from cascading is anti-correlated with the primary absorbed flux. The cascade light-curve is minimum at inferior conjunction ( $\phi\approx 0.72$). The non trivial double peaked structure of the lightcurve at phases 0.85-0.35 is due to competition in the cascade between absorption and inverse Compton emission. Absorption has a slight edge at superior conjunction ( $\phi\approx 0.06$), producing a dip at this phase. Elsewhere, the primary contribution dominates over the cascade emission. At lower energies ( $\epsilon_1<10~$GeV), the cascade contribution is undistinguishable from the primary source.

In this configuration, the cascade does add VHE gamma-ray emission close to superior conjunction but the expected contribution overestimates HESS observations. Decreasing the inclination of the system does not help: the cascade flux in the TeV energy band increases, since the primary source is on average more absorbed along the orbit (see Sect. 3 in Dubus 2006). For $i\la 30\hbox{$^\circ$ }$, the cascade contribution dominates the primary flux at every orbital phases in the VHE band. One-dimension cascades can be ruled out by the current HESS observations of LS 5039.

4 Conclusion

This paper explored the impact of one-dimensional pair cascading on the formation of the very high-energy radiation from gamma-ray binaries in general, LS 5039 specifically. A significant fraction of the total absorbed energy can be reprocessed at lower energy by the cascade, decreasing the global opacity of the primary source. Anisotropic effects also play a major role on the cascade radiation spectrum seen by a distant observer.

A large contribution from cascading is expected in LS 5039, large enough that it significantly overestimates the flux observed by HESS. One-dimensional cascading is too efficient in redistributing the absorbed primary flux and can be ruled out. However, the fact that it overestimates the observed flux means a more general cascade cannot be ruled out (it would have been if the HESS flux had been underestimated). If the ambient magnetic field is high enough ( $B\gg 10^{-8}~$G) the pairs will be deflected from the line-of-sight. For $B\ga 10^{-3}$G the Larmor radius of a TeV electron becomes smaller than the LS 5039 orbital separation and the pairs will be more and more isotropized locally. All of this will tend to dilute cascade emission compared to the one-dimensional case, which should therefore be seen as an upper limit to the cascade contribution at orbital phases where absorption is strong, particularly at superior conjunction. The initiated cascade will be three-dimensional as pointed out by Bednarek (1997). Each point in the binary system becomes a potential secondary source able to contribute to the total gamma-ray flux at every orbital phases. Cascade emission can still be sizeable all along the orbit in LS 5039, yet form a more weakly modulated background in the light-curve on account of the cascade radiation redistribution at other phases. The strength and structure of the surrounding magnetic field (from both stars) has a strong influence on the cascade (Bosch-Ramon et al. 2008b; Sierpowska & Bednarek 2005; Bosch-Ramon et al. 2008a). More realistic pair cascading calculations cannot be treated with the semi-analytical approach exposed here. Complementary investigations using a Monte Carlo approach are needed to better appreciate the cascade contribution in gamma-ray binaries.

Finally, the cascade will be quenched if the created pairs lose energy to synchrotron rather than inverse Compton scattering. This requires ambient magnetic fields $B \ga 5$ G, as found by equating the radiative timescales for a 1 TeV electron at periastron in LS 5039. Such ambient magnetic field strengths could be reached close to the companion star. In this case an alternative explanation is needed to account for the flux at superior conjunction. A natural one to consider is that the primary gamma-ray source is farther from the massive star. The VHE source would not be coincident with the compact object location anymore and would suffer less from absorption. In the microquasar scenario, Bednarek (2007) can account for consistent flux with HESS observations at superior conjunction if some electrons are injected well above the orbital plane (jet altitude $z>10~R_{\star}$). In addition to LS 5039, this possibility was also considered for the system Cyg X-1 by Bosch-Ramon et al. (2008b) and Zdziarski et al. (2009).

In practice, reality may consist of a complex three-dimensional cascade partly diluted and partly quenched depending upon position, angle and magnetic field configuration.

Acknowledgements
G.D. thanks A. Mastichiadis for discussions of triplet pair production. This work was supported by the European Community via contract ERC-StG-200911.

Appendix A: Pair production

The main equations for the pair production process are briefly presented here. Detailed calculations can be found in Gould & Schréder (1967), Bonometto & Rees (1971) and Böttcher & Schlickeiser (1997).

A.1 Kinematics and cross-sections

The interaction of a gamma-ray photon of energy $\epsilon _1$ and a soft photon of energy $\epsilon _0$in the observer frame leads to the production of an electron-positron pair if the total available energy in the center-of-mass (CM) frame is greater than the rest mass energy of the pair

\begin{displaymath}2\epsilon_1\epsilon_0\left(1-\cos\theta_0\right)\ge 4m_{\rm e}^2 c^4,
\end{displaymath} (A.1)

where $m_{\rm e}$ is the electron mass and $\theta_0$ the angle between the two incoming photons in the observer frame. It is useful to define the Lorentz invariant $s=\epsilon_1\epsilon_0\left(1-\cos\theta_0\right)/2$. Pairs are produced if $s\ge m_{\rm e}^2 c^4$ and the velocity $\beta$ of the electron-positron pair in the CM frame is $\beta=(1-m_{\rm e}^2 c^4/s)^{1/2}$.

The differential cross-section ${\rm d}\sigma_{\gamma\gamma}/{\rm d}(\beta\cos\theta'_1)$ in the CM frame depends on $\beta$ and the angle $\theta'_1$ between the outcoming electron-positron pair and the incoming photons. The full expression can be found in e.g. Bonometto & Rees (1971), Eq. (2.7). The differential cross-section presents a symmetric structure, peaked at $\cos\theta'_1=\pm 1$ and minimum for $\cos\theta'_1=0$. Electrons are mostly created in the same and opposite direction with respect to the incoming hard photon direction in the CM frame. The double peaked structure is enhanced with increasing energy ( $s\gg m_{\rm e}^2 c^4$) and becomes less pronounced close to the threshold ( $s\sim m_{\rm e}^2 c^4$). The integration over the angles gives the total pair production cross-section $\sigma_{\gamma\gamma}$, maximum close to the threshold (see Eq. (1) in Gould & Schréder 1967).

The construction of the CM frame with respect to the observer frame can be simplified if one of the incoming photons carries most of the energy. This case is appropriate in the present context. For $\epsilon_1\gg\epsilon_0$, the CM frame can be considered as propagating along the same direction as the high-energy photon. The velocity of the CM frame in the observer frame can be expressed as

\begin{displaymath}\beta'=\left(1-\frac{4s}{\epsilon^2_1}\right)^{1/2}\cdot
\end{displaymath} (A.2)

The total energy of say the electron $E_{\rm e}$ in the observer frame can then be formulated using the Lorentz transform from the CM to the observer frames

\begin{displaymath}E_{\rm e}=\gamma'\left[s^{1/2}+\beta'\left(s-m_{\rm e}^2 c^4\right)^{1/2}\cos\theta'_1\right],
\end{displaymath} (A.3)

providing a relation between $E_{\rm e}$ and $\cos\theta'_1$.

A.2 Rate of absorption and pair spectrum kernels

A gamma-ray photon going through a soft photon gas of density ${\rm d}n/{\rm d}\epsilon {\rm d}\Omega$ is absorbed at a rate per unit of path length l

\begin{displaymath}\frac{{\rm d}\tau_{\gamma\gamma}}{{\rm d}l}=\iint\frac{{\rm d...
...heta\right)\sigma_{\gamma\gamma}{\rm d}\epsilon {\rm d}\Omega.
\end{displaymath} (A.4)

The absorption rate gives the probability for a gamma ray of energy $\epsilon _1$ to be absorbed but does not give the energy of the pair created in the interaction.

Following Bonometto & Rees (1971), the probability for a gamma ray of energy $\epsilon _1$ to be absorbed between l and $l+{\rm d}l$ yielding an electron of energy between $E_{\rm e}$ and $E_{\rm e}+{\rm d}E_{\rm e}$ (with a positron of energy $E_{\rm e^+}\approx\epsilon_1-E_{\rm e}$ for $\epsilon_1\gg\epsilon$) is

\begin{displaymath}g_{\gamma\gamma}=\iint\frac{{\rm d}n}{{\rm d}\epsilon {\rm d}...
...\gamma\gamma}}{{\rm d}E_{\rm e}}{\rm d}\epsilon {\rm d}\Omega.
\end{displaymath} (A.5)

As with anisotropic inverse Compton scattering (Dubus et al. 2008), it is useful to consider the case of a monoenergetic beam of soft photons. The normalized soft photon density in the observer frame is

\begin{displaymath}\frac{{\rm d}n}{{\rm d}\epsilon {\rm d}\Omega}=\delta\left(\e...
...(\cos\theta-\cos\theta_0\right)\delta\left(\phi-\phi_0\right),
\end{displaymath} (A.6)

where $\delta$ is the Dirac distribution. Injecting Eqs. (A.6) into (A.5) gives the anisotropic pair production kernel, a convenient tool for spectral computations. The detailed calculation is presented in Appendix B and the complete expression given in Eq. (B.5). The pair production kernel has a strong angular dependence and a symmetric structure, centered at $E_{\rm e}=\epsilon_1/2$ and peaked at $E_{\rm e}=E_{\pm}$ (see Appendix B, Fig. B.1). The effect of the angle $\theta_0$ is reduced close to the threshold where the particles share equally the energy of the primary gamma-ray photon $E_{\rm e}\approx E_{\rm e^+}\approx\epsilon_1/2$. Far from the threshold, one particle carries away almost all the available energy $E_{\rm e}\approx\epsilon_1$.

The anisotropic kernel integrated over all the pitch angles, in the case of an isotropic gas of photons, is consistent with the kernel found by Aharonian et al. (1983). Note that a general expression for the anisotropic kernel valid beyond the approximation $\epsilon_1\gg\epsilon_0$ is presented in Böttcher & Schlickeiser (1997).

A.3 Pair density

The number of pair created per unit of length path and electron energy depends on the probability to create a pair and on the probability for the incoming gamma ray to remain unabsorbed up to the point of observation so that

\begin{displaymath}\frac{{\rm d}N_{\rm e}}{{\rm d}l {\rm d}E_{\rm e}}=\left\{g_{...
...n_1-E_{\rm e}\right)\right\}{\rm e}^{-\tau_{\gamma\gamma}(l)}.
\end{displaymath} (A.7)

Because of the symmetry in $g_{\gamma \gamma }$ and since electrons and positrons cannot be distinguished here, $g_{\gamma\gamma}\left(\epsilon_1-E_{\rm e}\right)=g_{\gamma\gamma}\left(E_{\rm e}\right)$. The integration over electron energy yields

\begin{displaymath}\frac{{\rm d}N_{\rm e}}{{\rm d}l}=2\left(\int g_{\gamma\gamma...
...au_{\gamma\gamma}}{{\rm d}l}{\rm e}^{-\tau_{\gamma\gamma}(l)}.
\end{displaymath} (A.8)

The total number of pairs produced by a single gamma ray bathed in a soft radiation along the path l up to the distance r is then

\begin{displaymath}N_{\rm e}(r)=2\left(1-{\rm e}^{-\tau_{\gamma\gamma}(r)}\right).
\end{displaymath} (A.9)

For low opacity $\tau_{\gamma\gamma}\ll 1$, pair production is inefficient and the number of particles produced tends to $\approx$ $2\tau_{\gamma\gamma}$. For high opacity $\tau_{\gamma\gamma}\gg 1$, a pair is always created.

Appendix B: Anisotropic pair production kernel

\begin{figure}
\par\includegraphics[width=7.5cm,clip]{12494fg7.eps}
\end{figure} Figure B.1:

Anisotropic pair production kernel $g_{\gamma \gamma }$ with $\epsilon _0$ set at 1 eV for a head-on collision ( $\theta _0=\pi $). The kernel is computed for $\epsilon _1=265$ GeV (dotted line), 300 GeV, 500 GeV, 1 TeV and 10 TeV (dashed line). The yielding of pairs occurs for $\epsilon _1\ge 260$ GeV.

Open with DEXTER
This section is dedicated to the calculation of the pair energy spectrum produced in the interaction between a single gamma-ray photon of energy $\epsilon _1$ and a mono-energetic beam of soft photons. The general expression in Eq. (A.5) can be reformulated using the relativistic invariant s

\begin{displaymath}g_{\gamma\gamma}=\frac{4}{\epsilon^2_1}\iiint \frac{s}{\epsil...
...amma}}{{\rm d}E_{\rm e}} {\rm d}\epsilon {\rm d}s {\rm d}\phi.
\end{displaymath} (B.1)

Combining the expression of $\beta$ with the equations Eqs. (A.2)-(A.3) and defining $x\equiv\gamma'^2$, the differential cross-section variables can be written as

\begin{displaymath}\beta(x)=\left(1-\frac{4 m_{\rm e}^2 c^4 x}{\epsilon^2_1}\rig...
...}-\epsilon_1}{\epsilon_1\left(1-\frac{1}{x}\right)^{1/2}}\cdot
\end{displaymath} (B.2)

The differential cross-section can then be expressed as
                   $\displaystyle \frac{{\rm d}\sigma_{\gamma\gamma}}{{\rm d}E_{\rm e}}$ = $\displaystyle \frac{{\rm d}\sigma_{\gamma\gamma}}{{\rm d}\left(\beta\cos\theta'_1\right)}\frac{{\rm d}\left(\beta\cos\theta'_1\right)}{{\rm d}E_{\rm e}}$  
  = $\displaystyle \frac{2}{\epsilon_1\left(1-\frac{1}{x}\right)^{1/2}}\frac{{\rm d}\sigma_{\gamma\gamma}}{{\rm d}\left(\beta\cos\theta'_1\right)}\cdot$ (B.3)

The complete general formula to compute the spectrum of the pair for a non-specified soft radiation field is

\begin{displaymath}g_{\gamma\gamma}=\frac{\epsilon_1}{4}\iiint\frac{1}{\epsilon^...
...m d}(\beta\cos\theta'_1)}{\rm d}\epsilon {\rm d}x {\rm d}\phi,
\end{displaymath} (B.4)

corresponding to Eq. (2.14) in Bonometto & Rees (1971). The injection of a mono-energetic and unidirectional soft photon density (Eq. (A.6)) in this last equation yields

\begin{displaymath}g_{\gamma\gamma}=\frac{2\left(1-\mu_0\right)}{\epsilon_1 \lef...
...ta\left(x_0\right),\beta\cos\theta'_1\left(x_0\right)\right\},
\end{displaymath} (B.5)

where $\mu_0\equiv\cos\theta_0$ and

\begin{displaymath}x_0=\frac{\epsilon_1}{2\epsilon_0\left(1-\mu_0\right)}\cdot
\end{displaymath} (B.6)

This expression is valid for $\epsilon_1\gg\epsilon_0$ and $s\ge m_{\rm e}^2 c^4$. The minimum E- and maximum E+ energy reached by the particles is set by the kinematics of the reaction and given by

\begin{displaymath}E_{\pm}=\frac{\epsilon_1}{2}\left[1\pm\left(1-\frac{1}{x_0}\r...
...{4 m_{\rm e}^2 c^4 x_0}{\epsilon^2_1}\right)^{1/2}\right]\cdot
\end{displaymath} (B.7)

Figure B.1 presents the pair production kernel for different incoming gamma-ray energy $\epsilon _1$.

Note that a kernel can be calculated as well for the absorption rate. Injecting Eqs. (A.6) into (A.4) is straightforward and gives

\begin{displaymath}\frac{{\rm d}\tau_{\gamma\gamma}}{{\rm d}l}=\left(1-\cos\theta_0\right)\sigma_{\gamma\gamma}\left(\beta\right).
\end{displaymath} (B.8)

References

All Tables

Table 1:   Mean energy of escaping pairs and radiated power efficiency of the cascade.

All Figures

  \begin{figure}
\par\includegraphics[width=8.5cm,clip]{12494fg1.eps}
\end{figure} Figure 1:

This diagram describes the system geometry. A gamma-ray photon of energy $\epsilon _1$ from the primary source (compact object) interacts with a soft photon of energy $\epsilon _0$ at a distance r from the source and R from the massive star (assumed point-like and mono-energetic), producing a pair $\rm e^+/e^-$ boosted toward a distant observer. The system is seen at an angle $\psi $.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=8cm,clip]{12494fg2.eps}
\end{figure} Figure 2:

Cascade development along the path to the observer. The primary source of photons, situated at the location of the compact object, has a power law spectral distribution with photon index -2 (dotted line). Spectra are computed using the parameters appropriate for LS 5039 at superior conjunction ( $d\approx 2 R_{\star }$, $R_{\star}=9.3~R_{\odot}$, $T_{\star}=39~000$ K) for $\psi =30\hbox {$^\circ $ }$. The transmitted spectrum, including cascade emission, is shown at various distances from the primary source: $r=R_{\star }/4$ (black dashed line), $R_{\star }/2$, $R_{\star }$, $2 R_{\star }$ (solid lines) and $r=+\infty $ (dotted-dashed line). Pure absorbed spectra are shown for comparison (light dashed line).

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=15.7cm,clip]{12494fg3.eps}
\end{figure} Figure 3:

Spectra as seen by an observer at infinity, taking into account the effect of cascading. Calculations are applied to LS 5039 at periastron for different viewing angle $\psi =30\hbox {$^\circ $ }$, $60\hbox {$^\circ $ }$, $90\hbox {$^\circ $ }$, $120\hbox {$^\circ $ }$ and $150\hbox {$^\circ $ }$. Left panel: complete spectra (solid line) are compared to the pure absorbed (light dashed line) and injected (dotted line) spectra. The contribution from the cascade is presented in the right panel.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=7.7cm,clip]{12494fg4.eps}
\end{figure} Figure 4:

Distribution of escaping pairs seen by a distant observer, depending on the viewing angle $\psi =30\hbox {$^\circ $ }$ (dashed line), $60\hbox {$^\circ $ }$, $90\hbox {$^\circ $ }$, $120\hbox {$^\circ $ }$ and $150\hbox {$^\circ $ }$ (dotted line). The binary parameters are the same as in Fig. 3.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=7.7cm,clip]{12494fg5.eps}
\end{figure} Figure 5:

Orbit-averaged spectra in LS 5039 at INFC ( $0.45<\phi <0.9$, grey lines) and SUPC ($\phi <0.45$ or $\phi >0.9$, black lines) and comparisons with EGRET (dark) and HESS (light) bowties (Aharonian et al. 2006a; Hartman et al. 1999). Dotted-dashed lines represent the primary source of gamma rays with pure absorption, injected at $r\equiv 0$, computed with the model described in Dubus et al. (2008) for a mono-energetic and point-like star. Dashed lines show the contribution from the cascade and thick solid lines the sum of the primary absorbed source and the cascade contributions.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=7.7cm,clip]{12494fg6.eps}
\end{figure} Figure 6:

Computed light-curves along the orbit in LS 5039, in the HESS energy band (flux $\ge $100 GeV). The cascade contribution (dashed line) is compared to the primary pure absorbed source (dotted-dashed line) and HESS observations. The thick solid line shows the sum of both components.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=7.5cm,clip]{12494fg7.eps}
\end{figure} Figure B.1:

Anisotropic pair production kernel $g_{\gamma \gamma }$ with $\epsilon _0$ set at 1 eV for a head-on collision ( $\theta _0=\pi $). The kernel is computed for $\epsilon _1=265$ GeV (dotted line), 300 GeV, 500 GeV, 1 TeV and 10 TeV (dashed line). The yielding of pairs occurs for $\epsilon _1\ge 260$ GeV.

Open with DEXTER
In the text


Copyright ESO 2009

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.