Free Access
Issue
A&A
Volume 496, Number 3, March IV 2009
Page(s) 751 - 758
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/200810871
Published online 09 February 2009

Multi-waveband photons and high energy neutrinos from the Crab nebula

X. C. Yang - L. Zhang

Department of Physics, Yunnan University, 650091 Kunming, PR China

Received 28 August 2008 / Accepted 12 January 2009

Abstract
We study multi-waveband photon emission properties and possible high-energy neutrino emission from the Crab nebula produced by pulsar's accelerated heavy nuclei. Heavy nuclei are accelerated in the outer magnetosphere of the Crab pulsar, and then suffer partial photo-disintegration in the non-thermal radiation fields of the pulsar's outer magnetosphere. The products (protons and neutrons) of the photo-disintegration and surviving heavier nuclei are injected into the expanding supernova remnant in which the surviving heavier nuclei are further accelerated in the pulsar wind zone. In such a case, relativistic hadrons and leptons are produced; the hadrons consist of direct protons and the protons produced through neutron decay in nuclei disintegration process, the leptons consist of positrons accelerated in a realistic magneto-sonic lepton-heavy iron plasma shock, electrons from the neutron decay, direct electrons, and secondary electron/positron pairs in a proton-proton interaction. We calculate the photo-disintegration rate for a nucleus of atomic number A releasing i nucleons and the hadronic and leptonic spectra in the Crab nebula. We then calculate the multi-waveband non-thermal photon spectrum and possible high-energy neutrino spectrum of the Crab nebula. Our results indicate that the photon spectrum up to $\sim$10 TeV has a leptonic origin, but it is difficult to distinguish leptonic and hadronic origins of the spectrum with energy greater than 10 TeV.

Key words: radiation mechanisms: non-thermal - ISM: individual objects: Crab nebula - gamma rays: theory

1 Introduction

Pulsar wind nebulae (PWNe) are generally believed to be a kind of Galactic cosmic ray source. The Crab nebula, being a well established gamma-ray source with a complex spectrum extending up to at least $\sim$80 TeV energies, is a possible object in which the particles can be accelerated up to $\sim$ 1015-16 eV (e.g. de Jager & Harding 1992; Aharonian & Atoyan 1995). For the Crab nebula, however, there are two possible kinds of models to explain the observed TeV emission properties: leptonic origin models (e.g. Jager & Harding 1992; Atoyan & Aharonian 1996; Zhang et al. 2008) and hadronic origin models (e.g. Cheng et al. 1990; Atoyan & Aharonian 1996; Bednarek & Prothero 1997; Bednarek & Bartosik 2003; a review see Bednarek 2007).

Photon emission with energy up to several tens of TeV from the Crab nebula has been usually interpreted with the framework of synchrotron self-Compton (SSC) model. In such a model, relativistic leptons inside the nebula produce soft synchrotron photons, and these soft photons are next up-scattered by the same leptons to gamma-ray energies through an inverse Compton (IC) process, possible soft photons also include microwave background radiation (MBR), infrared or the optical background inside the nebula. This model has been successfully applied to the best studied TeV gamma-ray source-the Crab nebula and other TeV PWNe (de Jager & Harding 1992; Atoyan & Aharonian 1996; Zhang et al. 2008). On the other hand, the self-consistent, time dependent hadronic-leptonic model has been proposed to account for the high energy processes inside PWNe (e.g. Bednarek & Prothero 1997; Bednarek & Bartosik 2003, 2005). In such a model, heavy nuclei are accelerated in the pulsar magnetosphere, and then suffer partial photo-disintegration in the non-thermal radiation fields of pulsar's outer magnetosphere, leading to the production of relativistic protons and neutrons. These protons and neutrons, as well as surviving heavier nuclei, are injected into the expanding supernova remnant. Inside the remnant, relativistic hadrons consist of direct protons and the protons produced through neutron decay in the nuclei disintegration process, the leptons are excited by the relativistic nuclei and consist of electrons from the neutron decay, positrons accelerated in the realistic magneto-sonic lepton-heavy iron plasma shock, and secondary electron/positron pairs in the proton-proton interaction.

In this paper we study the production of non-thermal photons and high-energy neutrinos in the Crab nebula based on the model of Bednarek & Protheroe (1997) and Bednarek & Bartosik (2003). In our treatment, we consider the relativistic proton spectrum from both direct proton production and indirect proton production from neutron decay through the disintegration process. We also consider four kinds of leptons: the positrons accelerated in a realistic magneto-sonic lepton-heavy iron plasma shock (Hoshino & Arons 1991; Hoshino et al. 1992; Gallant & Arons 1994), electrons accompanying neutron decay, direct electrons, and secondary e$^{\pm}$ pairs in proton-proton (pp) collisions (We call them nuclei-excited leptons except for the direct electrons, which accompany the nuclei). Different from the work in Bednarek & Bartosik (2003), we include the direct proton contribution which dominates the high-energy end of the proton spectrum, and use the complete form of the positron spectrum (i.e. a Maxwellian distribution plus a power-law tail) given by Hishino et al. (1992). Using spectra of both hadrons and leptons, we calculate the non-thermal spectrum of the Crab nebula, where the inclusive cross sections of stable secondary particles produced by proton-proton interaction, in parameterized formulae given by Kelner et al. (2006), are used. The organization of this paper is as follows. We describe the model in Sect. 2, and apply the model to the Crab nebula in Sect. 3. We give our conclusions and discussion in Sect. 4.

2 Leptons and hadrons produced in the nebula

We describe the basic physics processes which produce non-thermal emission from the Crab nebula in the framework of the hadronic-leptonic model.

2.1 Nucleon injection in the pulsar magnetosphere

Rapidly rotating neutron stars are known to convert a significant part of their rotation energy into radiation and the emission is observable as a periodic signal. It is possible that rotating magnetospheres of a neuron star can accelerate not only leptons but also nuclei extracted from the positively charged polar cap region (Bednarek & Bartosik 2004). The presence of heavy nuclei can explain the high energy phenomena observed around a pulsar (Gil et al. 2003). The amount of irons accelerated in the inner magnetosphere and pulsar wind zone, and injected in the pulsar wind nebula, can be written as (Bednarek & Protheroe 1997)

\begin{displaymath}\dot {N}_{\rm Fe}=\xi L/Z\Phi(B,P)
\end{displaymath} (1)

where $\xi$ is the fraction of total dissipated power taken by relativistic nuclei accelerated in the outer gap, Z=26 is the charge number of iron, B is the surface magnetic field, P is the pulsar period, L is the rate of kinetic rotational energy dissipated

\begin{displaymath}L=4\pi^2 I\frac {\dot p}{p^3} ,
\end{displaymath}

and

\begin{displaymath}\Phi(B,P)\approx5\times10^{16}\Big(\frac {B}{4\times10^{12}~\rm
G}\Big)\Big(\frac {P}{\rm s}\Big)^{4/3} V
\end{displaymath} (2)

is the potential difference across the outer gap (Cheng et al. 1990), where I is the neutron star moment of inertia which is assumed to be $1\times 10^{45}$ g cm2, B is the surface magnetic field strength which is inferred from the pulsar period P and period derivative $\dot P$. The spin period evolves from P0 as

\begin{displaymath}P=P_{0}\Big(1+\frac {t}{\tau_{0}}\Big)^{\frac{1}{n-1}} ,
\end{displaymath} (3)

where n is the braking index (e.g. Gaensler & Slane 2006).

2.2 Photodisintegration of nuclei in the pulsar soft photon field

It has been proposed that there is the possibility of photodisintegration of Fe nuclei during their propagation through the soft photon field from the Crab pulsar (Bednarek & Protheroe 1997). The Fe nuclei escape from the polar cap surface and move along magnetic field lines to enter the outer gap where they can be accelerated in the outer gap potential, as in the model of Cheng et al. (1986a). Since nuclei have much larger Larmor radii and much lower synchrotron and curvature radiations, they easily propagate out across the field lines along the outer gap. The nuclei suffer disintegration in collisions with the nonthermal radiation field in the outer gap produced by accelerated electron and positron cascades; the dimension of the outer gap is of the order of the radius of the light cylinder, namely, $l_{\rm
gap}\approx1.5\times10^8$ cm (Cheng et al. 1986b). A single proton or neutron emits in each nucleus disintegration, with the same Lorentz factor of the nucleus. The nucleons injected into the nebula, both the directly emitted protons and protons from neutron decay, are trapped and accumulate in the nebula.

 \begin{figure}
\par\includegraphics[width=6cm,clip]{0871fig1.eps}
\end{figure} Figure 1:

Reciprocal mean free path for 2656Fe in the Crab magnetosphere.

Open with DEXTER

The photo-disintegration rate for a nucleus of atomic number Areleasing i nucleons can be written as

\begin{displaymath}{R}_{{\rm A},i}=\frac{1}{2\gamma^2}\int_0^\infty\epsilon^{-2}...
...on^\prime
\epsilon^\prime\sigma_{{\rm A},i}(\epsilon^\prime) ,
\end{displaymath} (4)

where \(\gamma\) is the Lorentz factor of the nucleus, and \(\epsilon\), \(\epsilon^\prime\) the photon energies in the observer and nucleus rest frames respectively. $n(\epsilon)$ is the soft photon field in the pulsar magnetosphere and for the Crab pulsar is given by (Cheng et al. 1986)

\begin{displaymath}n(\epsilon)=n_0\epsilon^{-\nu} ,
\end{displaymath} (5)

where n0 is estimated using $n(\epsilon=m_{\rm e}c^2)\approx
3\times 10^{15}$ cm-3, $\nu=2.2$ for $h\nu\ge 1$, $\nu=1.8$ for $3\times 10^{-3}\le h\nu\le 1$, $\nu=1.25$ for $3\times 10^{-6}\le
h\nu\le 3\times 10^{-3}$, and $\nu=0.66$ for $h\nu\le 3\times
10^{-6}$. In Eq. (4), $\sigma_{{\rm A},i}(\epsilon^\prime)$ is the cross section of the photodisintegration; here we use the Gaussian parameterization for the cross section given by Puget et al. (1976) (also see the revision of Stecker & Salamon 1999). In such a soft photon field, we calculate a reciprocal mean free path $\lambda=c/R_{{\rm A},i}$ and show the reciprocal mean free path of 2656Fe in Fig. 1. The result is similar to that of the Lorentian curve approximation (Karakula & Tkaczyk 1993; Bednarek & Protheroe 1997).

2.3 Proton spectrum

When the photodisintegration rate for a nucleus of atomic number Areleasing i nucleons is given, the number of the i nucleons per unit energy can be estimated. For the treatment of entangled acceleration and disintegration process, we assume the nuclear Lorentz factor at position s is $\gamma_{\rm A}(s)$, where s is distance relative to the inner point of outer gap. For simplicity, the parallel component of the electrical field is assumed constant, as $E_{\parallel}\approx \Phi(B,P)/l_{\rm gap}$, so we have $\gamma_{\rm A}(s)\approx (ZE_{\parallel} s)/(A m_{\rm p}c^{2})$ and ${\rm d}s/{\rm d}\gamma_{\rm A}\approx (A m_{\rm p}c^{2} l_{\rm
gap})/(Z\Phi(B,P))$. Therefore, the number of i nucleons at energy Ei per unit energy per nucleus is approximated as

\begin{displaymath}N_i(\gamma_{i}(t), t)=\frac{{\rm d}N}{{\rm d}\gamma_i}=\frac{...
...mbda}=\frac{Al_{\rm gap}m_{\rm p}cR_{{\rm
A},i}}{Z\Phi(B,P)} ,
\end{displaymath} (6)

where i=n, p, and c is the light speed.

We now can construct and compute the total proton spectrum present in the Crab nebula including both the contribution of neutrons and protons. The present spectrum of protons from accelerated nuclei disintegrating in the nebula is the result of trapping and accumulation from the birth of pulsar to the present time and depends on the pulsar evolution since its birth, collisions and absorption by matter; this is estimated by optical depth. Using initial period P0=19 ms and braking index n=2.5 for the Crab pulsar (Lyne et al. 1993), the pulsar period evolves as Eq. (3). During the evolution of the Crab pulsar, the energy losses are mainly due to pp collision and adiabatic expansion. For the present Lorentz factor $\gamma_{\rm p}$ at time  $t_{\rm CN}$, its Lorentz factor at time t can be estimated as (Bednarek & Protheroe 1997)

\begin{displaymath}\gamma_{\rm p}(t)\approx\gamma_{\rm p}\frac{t+t_{\rm
CN}}{2tK^{\tau_{\rm pp}(t)}} ,
\end{displaymath} (7)

where the $\tau_{\rm pp}(t)$ is the optical depth for the collision of protons with matter between t and $t_{\rm CN}$ which can be approximated by $\tau_{\rm
pp}(t)\approx1.3\times10^{17}M_{1}v_{8}^{-3}(t^{-2}-t_{\rm
CN}^{-2})$, and K is the inelasticity coefficient in proton-proton collisions. In such a case, the energy spectra of protons extracted from Fe nuclei at an early time t are given by ${\rm d}N_{\rm
p}(\gamma_{\rm p}, t_{\rm CN})/{\rm d} E=[{\rm d}E(t)/{\rm d}E][{\rm d}N_{\rm p}/{\rm d}E(t)]$, where ${\rm d}E(t)/{\rm d}E=\gamma_{\rm p}(t)/\gamma_{\rm p}$ is estimated using Eq. (7). The optical depth due to absorption by matter can be expressed as $\tau_{\rm nH}\approx\sigma_{\rm pp}n_{\rm
H}r\approx8.6\times10^{14}M_{1}v_{8}^{-2}t^{-2}$, where $M=M_{1}~M_{\odot}$ is the mass ejected during the supernova explosion in units of solar mass, $n_{\rm H}$ is the number density of matter, and v=108v8 cm s-1 is the expansion velocity of the nebula; v=2000 km s-1 is used here. Following Bednareck & Protheroe (1997), the proton spectrum from neutron decay is
    $\displaystyle \Bigg[\frac{{\rm d}N_{\rm p}(\gamma_{\rm p}, t_{\rm CN})}{{\rm d}...
...^2}\int_{0}^{t_{\rm CN}}{\rm d}t\dot
{N}_{\rm Fe}(t){\rm e}^{-\tau_{\rm nH}(t)}$  
    $\displaystyle \quad\times\Bigg[N_{\rm n}(\gamma_{\rm p}(t),t)\gamma_{\rm p}(t)[1-\exp(vt/c\gamma_{\rm p}(t)\tau_{\rm n})]$  
    $\displaystyle \quad+\int_{t}^{t_{\rm CN}}{\rm d}t^{\prime}N_{\rm n}(\gamma_{\rm...
...(t^{\prime})\tau_{\rm n})}{c\gamma_{\rm p}(t^{\prime})\tau_{\rm n}}\Bigg] \cdot$ (8)

The first term is the contribution from neutrons initially decaying inside the nebula, and the second term gives the contribution of neutron initially decaying outside the nebula.

The direct contribution of proton emission is different from the disintegrated neutron component; all successive protons emitted from nuclei disintegrating and accelerating in the outer gap could be accelerated continuously and achieve almost the same energy determined by the outer gap potential $E_{\rm
p}(t)=Z/A\Phi(B,P(t))$, therefore ultimately the direct proton number dN with energy $E_{\rm p}(t)$ at pulsar period P is the sum of the differential energy spectra of disintegrated protons from E=0 to $E=E_{\rm p}(t)$, ${\rm d}N=\dot N_{\rm Fe}{\rm d}t\int_{0}^{E_{\rm
p}(t)}{\rm d}E{'}N_{\rm p}(E{'})$, i.e.

\begin{displaymath}\frac{{\rm d}N}{{\rm d}E_{\rm p}(t)}=\dot N_{\rm Fe}\frac{{\r...
...{\rm
p}(t)}\int_{0}^{E_{\rm p}(t)}{\rm d}E{'}N_{\rm p}(E{'}) .
\end{displaymath} (9)

Taking account of the energy loss, we reach the direct proton contribution of the energy spectrum as
$\displaystyle \Bigg[\frac{{\rm d}N_{\rm p}(\gamma_{\rm p}, t_{\rm CN})}{{\rm d}...
...{\rm nH}(t)}
\int_{0}^{\gamma_{\rm p}(t)}{\rm d}\gamma{'}N_{\rm p}(\gamma{'}) ,$     (10)

where $N_{\rm p}(\gamma{'})$ is given by Eq. (6), ${\rm d}P/{\rm d}E_{\rm p}(t)$ and ${\rm d}t/{\rm d}P$ can be derived using Eqs. (2) and (3), which are ${\rm d}P/{\rm d}E_{\rm
p}(t)=(3/20)\times10^{-16}(A/Z)(4\times10^{12}{\rm G}/B)(P/\rm
s)^{-1/3}$ and ${\rm d}t/{\rm d}P=(\tau_{0}(n-1)/P_{0})(1+t/\tau_{0})^{1-1/(n-1)}$.

The present proton energy spectrum has contributions both of neutron and proton components, the former is result of accumulation since the pulsar birth and the latter is the contribution at given time (or P). Therefore, the total proton spectrum can be written as

\begin{displaymath}\frac{{\rm d}N_{\rm p}}{{\rm d}E}=\Bigg[\frac{{\rm d}N_{\rm p...
...N_{\rm p}(\gamma_{\rm p}, t_{\rm
CN})}{{\rm d} E}\Bigg]_2\cdot
\end{displaymath} (11)

In Fig. 2, the total proton energy spectrum (solid line) is shown, which is the sum of the proton spectrum (dashed line) from neutron decay and the direct proton spectrum (dotted line). Although the cross section of neutrons and protons is similar, the present proton spectrum at low energy is mainly the contribution of low energy neutron decay accumulated in the pulsar evolution period, but at higher energy the direct proton emission contribution is greater than that of neutron decay. The results show that neglecting direct proton emission is not accurate. Following Bednarek & Protheroe (1997), we can check whether the protons remain in the nebula or diffuse into the galactic environment; the typical diffusion distance to the present time is $x_{\rm dif}\approx(cr_{\rm L}t_{\rm
CN}/3)^{1/2}$, where $r_{\rm L}$ is the Larmor radius of the protons, which is $r_{\rm L}\approx10^{-3}(m_{\rm
p}c^{2}/940~\mbox{MeV}) (\gamma_{\rm
p}/10^{5})(B/10^{-5}~\mbox{G})^{-1}~\mbox{pc}$. For the protons with $\gamma_{\rm p}=10^5$ in a magnetic field of $B\sim 10^{-5}$ G, $x_{\rm dif}\approx 0.36$ pc, comparable to the radius of the nebula. Indeed the magnetic fields of the inner part of the Crab nebula near 0.2 pc is even greater than $2\times10^{-3}$ G measured in the synchrotron wisps, so protons with Lorentz factors $\gamma <5\times 10^{5}$ could be trapped and accumulated in the Crab nebula.

 \begin{figure}
\par\includegraphics[width=6.5cm,clip]{0871fig2.eps}
\par\end{figure} Figure 2:

Present proton spectrum in the Crab nebula. The solid line represents the total proton spectrum, dashed and dotted lines represent the spectra of the protons from neutron decay (Eq. (8)) and the direct protons (Eq. (10)).

Open with DEXTER

2.4 Leptonic spectra excited by nuclei

In a PWN, leptons are excited by relativistic nuclei and consist of four components: positrons accelerated in a realistic magneto-sonic lepton-heavy iron plasma shock, electrons from neutron decay, direct electrons, and secondary electron/positron pairs in a proton-proton interaction (Bednarek & Bartosik 2003). We now give the spectra of the four kinds of leptons and show that the last three components are negligible compared to the first component.

For the first component, according to Hoshino et al. (1992), preferential non-thermal heating of positrons occurs in the realistic, magnetosonic, electron-positron-heavy ion shock. The gyrating reflected heavy ions dissipate their energy in the form of left-handed magnetosonic waves which are resonantly absorbed by the positron behind the iron reflection region. The positron energy spectrum can be described by a Maxwellian distribution with a supra-thermal power-law tail, i.e. (Hoshino et al. 1992)

\begin{displaymath}\left[\frac{{\rm d}N}{{\rm d}E_0{\rm d}t}\right]_1=\frac{N^*}...
...ac{m_{\rm e}c^2}{T_{\rm
e}}(\gamma-1)\right]+F(\gamma)\right\}
\end{displaymath} (12)

with

\begin{displaymath}F(\gamma)=n_{\rm
tail}\big(\frac{\gamma}{\gamma_1}\big)^{-\al...
...\;\;\mbox{for}~\gamma_{\rm
min}\le \gamma\le\gamma_{\rm max} ,
\end{displaymath} (13)

and $G(\gamma)=1-\exp\left[-(\frac{\gamma-1}{\gamma_{\rm min}})^r\right],
\mbox{for}~r\gg 1$. The cutoff function becomes a step function at $\gamma=\gamma_{\rm min}$, where $\gamma_{\rm min}=\gamma_{\rm Fe}$and $\gamma_{\rm max}\approx\gamma_{\rm Fe}(A/Z)(m_{\rm i}/m_{\rm
e})$, $m_{\rm e}$ and $m_{\rm i}$ are the electron and nucleon mass (Bednarek 2003). $n_{\rm tail}=0.076$ for $\alpha=2.4$ (Hoshino et al. 1992). To fit the observed radio emission we adopt $\gamma_1\approx 0.14\gamma_{\rm Fe}$, and $ T_{\rm
e}/\gamma_1m_{\rm e}c^2\approx0.14$. N* is the normalization coefficient and can be estimated by using

\begin{displaymath}\int [{\rm d}N/{\rm d}E_0{\rm d}t]_1E_0{\rm d}E_0=\chi\dot N_{\rm Fe}E_{\rm Fe}
\end{displaymath} (14)

where $\chi$ is the fractional energy of Fe transferred to positrons.

For the second component, electrons are produced through the neutron decay disintegrated from nuclei ( $n\rightarrow
p+{\rm e}^{-}+\overline{\nu}_{\rm e}$), where their Lorentz factors are assumed to be equal to those of their parent nuclei. In this case, the electron spectrum can be approximated as

\begin{displaymath}\left[\frac{{\rm d}N}{{\rm d}E_0{\rm d}t}\right]_2\approx\fra...
...\rm n}(E_{\rm n})[1-\exp(vt/c\gamma_{\rm
n}(t)\tau_{\rm n})] ,
\end{displaymath} (15)

where $E_0=\gamma_{\rm n} m_{\rm e}c^2$ is the energy of injected electrons at time t, $N_{\rm n}(E_{\rm n})$ is the emission neutron energy spectrum with Lorentz factor $\gamma_{\rm n}$, and the coefficient is because of the electron energy spectrum transformed from the decayed neutron spectrum (i.e. the same number electrons with the same Lorentz factor of a decayed neutron $\frac{{\rm d}E_{\rm
n}}{{\rm d}E_{\rm e}}=\frac{m_{\rm n}}{m_{\rm e}}$).

According to Bednarek (2003), we consider the direct electron spectrum. We assume that the spin down energy transfer $1-\xi$ to monoenergy electrons with primary energy $E_0=\gamma_{\rm Fe}m_{\rm
e}c^2$, therefore the direct electron spectrum is given by

\begin{displaymath}\left[\frac{{\rm d}N}{{\rm d}E_0{\rm d}t}\right]_3\approx\frac{L(1-\xi)}{Z_{\rm
e}\Phi(B,P)}\delta(E-E_0),
\end{displaymath} (16)

where $Z_{\rm e}$ is the charge number of the electron. We also consider the spectrum of the secondary electron/positron pairs produced in proton-proton interaction ( $pp\rightarrow\pi^{\pm}\rightarrow \mu^{\pm} \rightarrow {\rm e}^{\pm}$), but its contribution to the leptonic spectrum is negligible (Fig. 3).

 \begin{figure}
\par\includegraphics[width=6.3cm,clip]{0871fig3.eps}
\end{figure} Figure 3:

Present spectra of leptons (solid line) in the Crab nebula. For comparison, we plot the spectra of the Maxwellian distribution (dash-dash dotted line) and the power-law tail (dotted line); dashed, dot-dot dashed, and dot dashed lines represent the second to fourth component spectra.

Open with DEXTER
Therefore, when we consider the leptons which consists of the first, second and third components, the leptonic spectrum is

\begin{displaymath}\frac{{\rm d}N}{{\rm d}E_0{\rm d}t}=\left[\frac{{\rm d}N}{{\r...
...ht]_2+\left[\frac{{\rm d}N}{{\rm d}E_0{\rm d}t}\right]_3 \cdot
\end{displaymath} (17)

The leptons with an energy E0 at an earlier time t are injected into the nebula, and the energy E of the leptons inside the nebula at the present time  $t_{\rm obs}$ is determined by the energy loss rate and expansion of the nebula. Leptons injected into the nebula suffer mainly the following energy losses: bremsstrahlung, synchrotron, and inverse Compton losses as well as adiabatic loss due to the expansion of the nebula. The rate of their energy losses reads as

\begin{displaymath}-\frac{{\rm d}E}{{\rm d}t}=\dot{E}_{\rm brem}+\dot{E}_{\rm ad}+\dot{E}_{\rm
syn}+\dot{E}_{\rm comp},
\end{displaymath} (18)

where $\dot{E}_{\rm brem}\approx m_{\rm p}N/(m_{\rm e}X)E$ describes the bremsstrahlung loss, where X is the radiation length and Nis the density of the medium (Bednarek & Bartosik 2003); $\dot{E}_{\rm ad}=E\dot{R}_{\rm neb}/R_{\rm neb}$ is the adiabatic loss; $\dot{E}_{\rm syn}=(4/3)c\sigma_{\rm T}\gamma_{\rm
e}^2\beta_{\rm e}^2B_{\rm neb}^2/8\pi$ represents the synchrotron energy loss, where $\sigma_{\rm T}$ is the Thomson cross section and $B_{\rm neb}$ is the magnetic field strength. For the inverse Compton scattering (ICS), leptons will lose energy by Compton upscattering ambient soft photon fields. Soft photons include the synchrotron radiation created by the leptons in the magnetic field of the nebula, the cosmic microwave background field, and the infrared photons produced by the the dust inside the nebula. The energy losses of leptons during the ICS in the ``Klein-Nishina regime'' can be safely neglected relative to the synchrotron energy losses at the energy range in our paper. For a realistic electron or positron traversing a blackbody photon gas with temperature T, the energy loss rate is described by Blumenthal & Gould (1970) in the Thomson limit with a lowest-order correction with realistic quantum electrodynamics (Klein-Nishina formula); it is applicable in the Thomson and near Klein-Nishina regimes. We adopted the formula given by Gould (1975) $\dot{E}_{\rm comp}=4/3\sigma_{\rm
T}acT^2\gamma_{\rm e}^2(1-24.15\gamma_{\rm e}kT/(m_{\rm e}c^2)$, where $a=7.5678\times10^{-15}$ erg cm-3 K-4; the last factor in the equation accounts for deviation from the Thomson limit when the deviation is small.

The leptonic spectrum at the present time is obtained by summing over the spectra from leptons injected at an early time and over all the time up to the time of the present observation, which is given by Bednarek & Bartisik (2003) (also see Bednarek 2007):

\begin{displaymath}\frac{{\rm d}N(t_{\rm obs})}{{\rm d}E}=\int_{0}^{t_{\rm
obs}}J(t')\frac{{\rm d}N}{{\rm d}E_{0}{\rm d}t}{\rm d}t
\end{displaymath} (19)

where ${\rm d}N/{\rm d}E_{0}{\rm d}t$ is the injection spectrum (Eq. (17)) of the leptons at time t, $t'=t_{\rm obs}-t$ and J(t')=E0/E. In Fig. 3, we show the current spectra of the leptons inside the nebula after the birth of the Crab pulsar, where the solid line represents the total spectrum. For comparison, we plot the spectra of the Maxwellian distribution (dash-dash dotted line) and the power-law tail (dotted line). In Fig. 3, dashed, dot-dot dashed, and dot dashed lines represent the second to fourth component spectra. It can be seen that the Maxwell component dominates in the lower energy range and the power law component dominates in higher energy range.

3 Photon emission inside the nebula

Having calculated the spectra of both leptons and hadrons, which depend on the time of evolution after the supernova explosion, we can compute the photon spectra produced by the leptons in different radiation processes and the hadrons in pp interactions inside the nebula. The main leptonic processes which produce photons inside the nebula include synchrotron radiation, bremsstrahlung, and inverse Compton (IC) scattering. It is possible that photons produced by inverse Compton scattering of the high-energy leptons on synchrotron photons, microwave background photons, and infrared photons reach to very high energies. Despite different accuracies of IC cross section determination, as well as different distributions of magnetic fields and photon fields in nebula used in the past by different authors, the conclusion is unambiguous that the lepton component has a significant contribution to TeV $\gamma $ rays (e.g. Gould 1965; Jones 1968; Rieke & Weekes 1969; Kennel et al. 1984; de Jager & Harding 1992; Atoyan & Aharonian 1996; Zhang et al. 2008). On the other hand, hadrons produce high-energy photons and neutrinos in p-p interactions. Below, we briefly describe the relevant processes and give our results for the Crab nebula.

3.1 Photons produced by leptonic processes

For the synchrotron radiation, the emissivity is given by

$\displaystyle Q_{\rm syn}(E_{\gamma},t)=\frac{2\sqrt{3}{\rm e}^{3}B}{\hbar
E_{\...
...}}^{E_{\rm e,max}}{\rm d}E_{\rm e}J_{\rm e}(E_{\rm e},t)F(E_{\gamma}/E_{\rm c})$     (20)

where $J_{\rm e}(E_{\rm e})=(c\beta/4\pi){\rm d}n(t_{\rm obs})/{\rm d}E$, where ${\rm d}n(t_{\rm obs})/{\rm d}E$ is the differential particle density over the nebula volume, $\theta$ is the electron pitch angle, $E_{\rm
c}=(3/2)\gamma^2(heB\sin\theta/m_{\rm e}c$, B is the local magnetic field strength, and $F(x)=x\int_{x}^{\infty}K_{5/3}(t){\rm d}t$with $x=E_{\gamma}/E_{\rm c}$, where K5/3(t) is the modified Bessel function of order 5/3 (Blumenthal 1970).

The electron-nucleon bremsstrahlung emissivity can be written as

\begin{displaymath}Q_{\rm brem}(E_{\gamma},t)=4\pi n_{\rm H}\int_{E_{\rm
e,min}}...
...)\Big(\frac{{\rm d}\sigma}{{\rm d}E_{\gamma}}\Big)_{\rm e,p} ,
\end{displaymath} (21)

where $({\rm d}\sigma/{\rm d}E_{\gamma})_{\rm e,p}$ is the differential electron-proton bremsstrahlung cross section, which is (Koch & Motz 1959)
                    $\displaystyle \Bigg(\frac{{\rm d}\sigma}{{\rm d}E_{\gamma}}\Bigg)_{\rm e,p}(E_{\rm e},E_{\gamma})$ = $\displaystyle \frac{r_{0}^{2}\gamma_{f}\beta_{f}}{137E_{\gamma}\gamma_{i}\beta_...
...2}\beta_{f}^{2}}{\gamma_{i}^{2}\beta_{i}^{2}\gamma_{f}^{2}\beta_{f}^{2}}
\Bigg)$  
    $\displaystyle +\frac{\alpha_{i}\gamma_{f}}{\gamma_{i}^{3}\beta_{i}^{3}}+\frac{\...
...a_{f}^{3}}
-\frac{\alpha_{i}\alpha_{f}}{\gamma_{i}\beta_{i}\gamma_{f}\beta_{f}}$  
    $\displaystyle +L\Bigg\{\frac{8}{3\beta_{i}\beta_{f}}+\frac{\epsilon_{\gamma}^{2}(1+\beta_{i}^{2}\beta_{f}^{2})}{\gamma_{i}\beta_{i}^{3}\gamma_{f}\beta_{f}^{3}}$  
    $\displaystyle +\frac{\epsilon_{\gamma}}{2\gamma_{i}\beta_{i}\gamma_{f}\beta_{f}...
...ac{\alpha_{f}(\gamma_{i}+\gamma_{f}\beta_{f}^{2})}{\gamma_{f}^{2}\beta_{f}^{3}}$  
    $\displaystyle +\frac{2\epsilon_{\gamma}}{\gamma_{i}\beta_{i}^{2}\gamma_{f}\beta_{f}^{2}}\Bigg]\Bigg\}\Bigg)$ (22)

with

\begin{displaymath}\alpha_{i(f)}=\ln\left[\gamma_{i(f)}^{2}(1+\beta_{i(f)})^{2}\right]
\end{displaymath} (23)

and

\begin{displaymath}L=2\ln\left(\frac{\gamma_{i}\gamma_{f}+\gamma_{i}\beta_{i}\gamma_{f}\beta_{f}-1}{\epsilon_{\gamma}}\right) ,
\end{displaymath} (24)

where $\gamma_{i(f)}$ is the initial(final) electron Lorentz factor, the dimensionless photon energy $\epsilon_{\gamma}=E_{\gamma}/m_{\rm e}c^{2}$, and $r_{0}={\rm e}^{2}/m_{\rm e}c^{2}$ is the classical electron radius.

For the inverse Compton process, we use the full Klein-Nishina cross section for relativistic electrons to calculate the emissivity. The emissivity can be written as

$\displaystyle Q_{\rm {comp},~j}(E_\gamma, t)=4\pi\int^{\infty}_{0}{\rm d}\epsil...
...sh}}
{\rm d}E_{\rm e}J_{\rm e}(E_{\rm e},t)F(\epsilon, E_{\gamma}, E_{\rm e}) ,$     (25)

where $E_{\rm e, thresh}=[E_{\gamma}+(E^2_{\gamma}+E_{\gamma}(m_{\rm
e}c^2)^2/\epsilon)^{1/2}]/2$ is the lowest energy that electrons can scatter a target photon with energy $\epsilon$ to energy $E_{\gamma}$, $\epsilon$ is the target photon energy, and $n_{\rm {j}}$ is the number density of the soft photon component j. The energy distributions of the components of cosmic microwave background (CMB), IR and star light can be given by

\begin{displaymath}n_{\rm {j}}(\epsilon)=\frac {15U_{\rm {j}}}{(\pi k
T_{\rm {j}})^4}\frac{\epsilon^2}{\exp(\epsilon/kT_{\rm {j}})-1}\cdot
\end{displaymath} (26)

The energy density of the synchrotron photons is (Atoyan & Aharonian 1996)

\begin{displaymath}n_{\rm {syn}}(r, \epsilon)=\frac{Q_{\rm {syn}}(\epsilon)}{4\pi
r_{0}^2c}U(x) ,
\end{displaymath} (27)

with

\begin{displaymath}U(x)=\frac{3}{2}\int_0^1\frac{y}{x}\ln\frac{x+y}{\vert x-y\vert}{\rm d}y,
\end{displaymath} (28)

where $x\equiv r/r_0$, $Q_{\rm {syn}}(\epsilon)$ is the synchrotron emissivity. In Eq. (27), it is assumed that the synchrotron emission is homogenous in the nebula with an average emission volume $V_{0}=\frac{4}{3}\pi r_{0}^3$. Function $F(\epsilon, E_{\gamma},
E_{\rm e})$ is given by
$\displaystyle F(\epsilon, E_{\gamma}, E_{\rm e})={3\sigma_{\rm T}\over 4 (E_{\r...
...psilon}\Big[2q\ln q +(1+2q)(1-q)
+{(\Gamma q)^2(1-q)\over 2(1+\Gamma q)}\Big] ,$     (29)

with $\Gamma=4\epsilon(E_{\rm e}/m_{\rm e}c^2)/m_{\rm e}c^2$, $q=E_1/\Gamma (1-E_1)$ with $E_1=E_{\gamma}/E_{\rm e}$ and $1/4(E_{\rm e}/m_{\rm e}c^2) < q <1$.

3.2 Photons and neutrinos produced by hadronic processes

Table 1:   Main parameters for calculating photon and neutrino production.

The protons accumulated in a PWN can generate high-energy photons and neutrinos in p-p interactions ( $pp\rightarrow pp
+n_{\pi}\pi^{\pm,0}$, where $n_{\pi}$ is the pion multiplicity). In this interaction process, high-energy photons are produced through the decay of neutral pions. Kelner et al. (2006) presented a simple analytical approximation with very good accuracy for the energy distribution of secondary products in the energy rangy above 0.1 TeV and proposed an approximate procedure for the calculations toward lower energies. Here we use energy-dependent parameterizations of the inelastic total cross section and differential $\pi ^0$production for proton-proton interactions (Kelner et al. 2006) to calculate the differential photon flux. For the higher energy region of $E_{\rm p} > 0.1$ TeV and $E_{\gamma}/E_{\rm p}\ge 10^{-3}$, the photon emissivity can be expressed as
$\displaystyle Q^{\rm pp}_{\pi}(E_{\gamma})=4\pi cn_{\rm H}\int_{E_{\gamma}}^{\i...
...rac{E_{\gamma}}{E_{\rm p}},E_{\rm p}\right)\frac{{\rm d}E_{\rm p}}{E_{\rm p}} ,$     (30)

where $J_{\rm p}(E_{\rm p}, t)=(c\beta/4\pi){\rm d}N_{\rm p}(E_{\rm
p})/{\rm d}E_{\rm p}$, $n_{\rm H}$ is the density of the ambient hydrogen, $\sigma_{\rm inel}(E_{\rm p})$ is the inelastic part of the total cross section of p-p interactions, which is described as

\begin{displaymath}\sigma_{\rm inel}(E_{\rm
p})=34.3+1.88L+0.25L^{2}\times\left[...
...(\frac{E_{\rm th}}{E_{\rm
p}}\right)^{4}\right]^2\;\mbox{mb} ,
\end{displaymath} (31)

where $E_{\rm th}$ is the threshold energy of production of neutral pions, where $L=\ln(E_{\rm p}/1~{\rm TeV})$. $F_{\gamma}(x,E_{\rm p})$ is the number of photons in the interval $(x,x+{\rm d}x)$ per collision, which is given in Eq. (58) of Kelner et al. (2006), where $x=E_{\gamma}/E_{\rm p}$. For the lower energy region of $1~\mbox{GeV}\ll E_{\rm p} \le 100$ GeV, Aharonian & Atoyan (2000) proposed that the photon spectrum can be treated using a $\delta$-functional approximation and the photon emissivity in this energy region can be written as

\begin{displaymath}Q^{\rm pp}_{\pi}(E_{\gamma})=2\int^{\infty}_{E_{\pi, \rm
min}...
..._{\pi}\frac{q_{\pi}(E_{\pi},t)}{\sqrt{E^2_{\pi}-m_{\pi}c^2}} ,
\end{displaymath} (32)

where $E_{\pi, \rm min}=E_{\gamma}+m^2_{\pi}c^4/4E_{\gamma}$, and
$\displaystyle q_{\pi}(E_{\pi},t)=\tilde{n}\frac{4\pi n_{\rm H}}{K_{\pi}}
\sigma...
...2+E_{\pi}/K_{\pi}\right)
J_{\rm p}\left(m_{\rm p}c^2+E_{\pi}/K_{\pi},t\right) ,$     (33)

where $\tilde{n}$ and $K_{\pi}$ are parameters that depend weakly on the proton energy. Kelner et al. (2006) pointed out that the photon emissivity in the energy range $E_{\rm p}\ge 0.1$ TeV is calculated using Eq. (30), and at lower energy using Eq. (33) with $K_{\pi}=0.17$. The parameter $\tilde{n}$ is treated as a free parameter determined from the conditions that match the emissivity given by Eq. (30) at $E_{\rm
p}=0.1$ TeV.

In the p-p collision, the pion emissivities satisfy the following relationship: $Q_{\pi^{0}}\approx Q_{\pi^{+}}\approx Q_{\pi^{-}}$. The muonic neutrinos are produced at the decay of charged pions $\pi\rightarrow\mu\nu_{\mu}$, and the corresponding third particle decay yields muonic and electronic neutrinos as the mode $\mu\rightarrow e\nu_{\rm e}\nu_{u}$. Taking neutrino oscillations into account, all muonic neutrinos and electronic neutrinos are incorporated and the emissivity of the muonic neutrinos is

\begin{displaymath}Q_{\nu_{\mu}}=P_{\mu\mu}Q_{\nu_{\mu}}+P_{{\rm e}\mu}Q_{\nu_{\rm e}},
\end{displaymath} (34)

where $P_{\mu\mu}=0.39$ and $P_{\mu\mu}=0.22$ are the oscillation probability and the emissivity of muonic neutrinos and electronic neutrinos are calculated by using the formulae given by Kelner et al. (2006). In our calculations, we did not consider the neutrino from neutron decay $p\rightarrow p+{\rm e}^-+\overline{\nu}_{\rm e}$ since the laboratory energy is low(estimated by $E_{\nu}=\gamma\Delta$, $\Delta=m_{\rm n}-m_{\rm p}-m_{\rm e}=0.71$ MeV), it mainly contributes to the energy below several ten GeV.

3.3 Calculation results

We now calculate the photon spectrum and muonic neutrino flux for the Crab nebula. For the photons produced in the nebula, the spectrum observed on the Earth can be expressed as

$\displaystyle \frac{{\rm d}N}{{\rm d}E_{\gamma}}=\frac{V_{\rm neb}}{4\pi d^2}[Q...
...\rm brem}(E_\gamma)
+Q_{\rm {comp},~j}(E_\gamma)+Q^{\rm pp}_{\pi}(E_{\gamma}] ,$     (35)

where $V_{\rm neb}$ is the volume of the nebula and d is the distance from the earth to the nebula. The spectrum of synchrotron and Compton radiation integrated over the nebula volume is independent of $R_{\rm neb}$, but the bremsstrahlung and pp production is effected by $n_{\rm H}=M/V_{\rm neb}m_{\rm p}$; radio and optical observations indicate that the nebula has an angular size of $5\times7$ arcmin2, which corresponds to a spatial size of $3\times4$ pc2 for d=2 kpc. Protons are efficiently trapped by the dense filaments as suggested by Atoyan & Aharonian (1996) and filaments with a density $\sim$500 cm-3 exist in the Crab nebula (Davidson & Fesen 1985). Therefore we use $n_{\rm H}=38$ cm-3 which is slightly greater than the averaged density in our calculation. For the muonic neutrinos produced in the nebula, the the neutrino flux observed on the Earth is

\begin{displaymath}F_{\nu_{\mu}}(E)=\frac{V_{\rm neb}}{4\pi d^2}Q_{\nu_{\mu}} .
\end{displaymath} (36)

We calculated the multiwaveband photon and high energy neutrino spectra of the Crab nebula and the main parameters used here are listed in Table 1. In Fig. 4, we show the comparison of the predicted photon spectrum with the observed data for the Crab nebula. The radio data are taken from Baars & Hartsuijker (1972), millimeter and IR data from Strom & Greidanus (1992), Bandiera et al. (2002), and Green et al. (2004), soft X-ray to $\gamma $-ray data from Kuiper et al. (2001), and VHE $\gamma $-ray data from Aharonian et al. (2004, 2006) and Albert et al. (2008). In Fig. 5, we address the comparison of the predicted result with observed high-energy data. It can be seen that our results can reproduce the observed multi-waveband spectrum well. Our results indicate that the photon spectrum up to $\sim$10 TeV has a leptonic origin, but the spectrum with energy greater than 10 TeV consists of both hadronic and leptonic components.

 \begin{figure}
\par\includegraphics[width=6.7cm,clip]{0871fig4.eps}
\end{figure} Figure 4:

Comparison of the calculated spectrum with observed data for the Crab nebula. The radio data are taken from Baars & Hartsuijker (1972), millimeter and IR data from Strom & Greidanus (1992), Bandiera et al. (2002), and Green et al. (2004), soft X-ray to $\gamma $-ray data from Kuiper et al. (2001), and VHE $\gamma $-ray data from Aharonian et al. (2004, 2006) and Albert et al. (2008). Dash-dotted, dot-dot dashed, and dash-dash dotted curves represent the synchrotron, bremsstrahlung, and self synchrotron Compton emission spectra of the leptons, curves 1, 2, and 3 are the spectra of the inverse Compton scattering of the leptons with the microwave background, infrared photons, and optical photons, the dashed curve represents the $\pi ^0$-decay spectrum in the pp interaction, and the solid line is total emission.

Open with DEXTER

 \begin{figure}
\par\includegraphics[width=6.7cm,clip]{0871fig5.eps}
\end{figure} Figure 5:

Comparison of the calculated spectrum with observed high-energy data for the Crab nebula. Symbols are the same as in Fig. 4.

Open with DEXTER

In Fig. 6, we show the predicted muonic neutrino flux observed on the Earth and the atmospheric neutrino background (Lipari 1993) within $1^{\circ }$ as indicated by the shaded band. Compared to the result given by Bednarek (2003), our result is slightly lower. Bednarek (2003) estimated that the number of $\mu$neutrinos observed by the 1 km2 detector during 1 year is about 1.3 for the Crab nebula.

 \begin{figure}
\par\includegraphics[width=6.5cm,clip]{0871fig6.eps}
\end{figure} Figure 6:

Predicted spectrum of muonic neutrinos observed on the Earth. The atmospheric neutrino background within $1^{\circ }$ of the source is indicated by the shaded band.

Open with DEXTER

4 Discussion and conclusion

We have studied multi-waveband photon emission and possible high-energy neutrino emission from the Crab nebula in the framework of the hadronic-leptonic model. In this model, we considered nucleon injection in the pulsar magnetosphere, proton production by the photo-disintegration process, and leptonic production excited by the nuclei in detail. The proton spectrum consists of two components: the proton component from neutron decay and the direct emitting proton component. The first component dominates over the second component in the lower energy region, but the second component has its main contribution in the higher energy region (see Fig. 2). It should be noted that previous works (e.g. Bednarek & Prothero 1997; Bednarek & Bartosik 2003) did not consider the contribution of the second component. For the leptons excited by the nuclei, the spectrum consists of four components: positrons accelerated in realistic magneto-sonic lepton-heavy iron plasma shocks, electrons from neutron decay, direct electrons, and secondary electron/positron pairs in the pp interaction (see Fig. 3), it can be seen from this figure that the first component dominates over other components, which is different from that given by Bednarek & Bartosik (2003) since the complete form of the positron spectrum given by Hoshino et al. (1992) is included in our calculation.

After determining the spectra of both hadrons and leptons, we calculated the multi-waveband photon spectrum of the Crab nebula (see Figs. 4 and 5). The calculated spectrum up to TeV energy could be explained by upscattering soft photon fields by realistic leptons accompanying nuclei, where the soft photon fields mainly include synchrotron radiation, microwave background, infrared and optical photons. However our predicted spectrum reveals that the hadronic component has an important contribution to gamma rays in the energy range above several TeV. Bednarek & Bartosik (2003) found that the emission from the Crab nebula can be well fitted by the composition of the $\gamma $-ray emission produced by leptons (<10 TeV) and nuclei (>10 TeV). Compared to the results of Bednarek & Bartosik (2003), our results show that it is difficult to clearly distinguish leptonic and hadronic origins in the $\sim$10 TeV energy range (see Fig. 5). For the $\mu$ neutrino emission, our predicted flux (see Fig. 6) is lower than that given by Bednarek (2003).

We now discuss the possible influence of our treatments on the predicted photon and neutrino fluxes. In this paper, we have given the treatment of disintegration as accelerating along the magnetic line with a non-perpendicular electric field in the outer gap. The nucleon energy spectrum from nucleus disintegration is calculated simply by Eq. (6) assuming a constant electronic field. Because the real $\frac{{\rm d}s}{{\rm d}\gamma_{A}}$ near the outer range of the gap should be larger than our average value, the calculation will mainly underestimate the the proton energy spectrum at high energy. If detailed information about $E_{\parallel}(s)=f(s)$ is known, the nucleon energy spectrum could be calculated precisely as could its more distinct contribution to the very high photon spectrum. It should be noted that the direct proton contribution at high energy is significant irrespective of the calculation technique. Another factor which could effect the nucleon energy spectrum and thus the maximum positron energy is that the kinetic energy of nuclei entering the outer gap is assumed to be zero; this could lead to a decrease in the disintegration rate and thus nuclei maximum energy (or maximum positron energy), resulting in the decrease of fluxes at high energy of pp photon emission and synchrotron radiation.

Acknowledgements

This work is partially supported by the National Natural Science Foundation of China under grant 10778702 and a 973 Program (2009CB824800).

References

All Tables

Table 1:   Main parameters for calculating photon and neutrino production.

All Figures

  \begin{figure}
\par\includegraphics[width=6cm,clip]{0871fig1.eps}
\end{figure} Figure 1:

Reciprocal mean free path for 2656Fe in the Crab magnetosphere.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=6.5cm,clip]{0871fig2.eps}
\par\end{figure} Figure 2:

Present proton spectrum in the Crab nebula. The solid line represents the total proton spectrum, dashed and dotted lines represent the spectra of the protons from neutron decay (Eq. (8)) and the direct protons (Eq. (10)).

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=6.3cm,clip]{0871fig3.eps}
\end{figure} Figure 3:

Present spectra of leptons (solid line) in the Crab nebula. For comparison, we plot the spectra of the Maxwellian distribution (dash-dash dotted line) and the power-law tail (dotted line); dashed, dot-dot dashed, and dot dashed lines represent the second to fourth component spectra.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=6.7cm,clip]{0871fig4.eps}
\end{figure} Figure 4:

Comparison of the calculated spectrum with observed data for the Crab nebula. The radio data are taken from Baars & Hartsuijker (1972), millimeter and IR data from Strom & Greidanus (1992), Bandiera et al. (2002), and Green et al. (2004), soft X-ray to $\gamma $-ray data from Kuiper et al. (2001), and VHE $\gamma $-ray data from Aharonian et al. (2004, 2006) and Albert et al. (2008). Dash-dotted, dot-dot dashed, and dash-dash dotted curves represent the synchrotron, bremsstrahlung, and self synchrotron Compton emission spectra of the leptons, curves 1, 2, and 3 are the spectra of the inverse Compton scattering of the leptons with the microwave background, infrared photons, and optical photons, the dashed curve represents the $\pi ^0$-decay spectrum in the pp interaction, and the solid line is total emission.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=6.7cm,clip]{0871fig5.eps}
\end{figure} Figure 5:

Comparison of the calculated spectrum with observed high-energy data for the Crab nebula. Symbols are the same as in Fig. 4.

Open with DEXTER
In the text

  \begin{figure}
\par\includegraphics[width=6.5cm,clip]{0871fig6.eps}
\end{figure} Figure 6:

Predicted spectrum of muonic neutrinos observed on the Earth. The atmospheric neutrino background within $1^{\circ }$ of the source is indicated by the shaded band.

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.