2.1 The wind

In the wind, the magnetic field is dominated by the azimuthal component, whose magnitude varies as

(2.1) $$\begin{eqnarray}B_{\unicode[STIX]{x1D719}}(r,\unicode[STIX]{x1D703})=B_{\text{lc}}\left(\frac{r_{\text{lc}}}{r}\right)\sin ^{n}\unicode[STIX]{x1D703},\end{eqnarray}$$

where $r_{\text{lc}}=c/\unicode[STIX]{x1D6FA}$ is the radius of light cylinder, $B_{\text{lc}}$ is the typical value of $B$ at the light cylinder and $\unicode[STIX]{x1D6FA}$ is the angular velocity of the pulsar. In the monopole model of the pulsar magnetic field $n=1$ but there is no easy way of finding its value in the more realistic dipolar model. Recent numerical simulations suggest that it could be $n=2$ (Tchekhovskoy, Philippov & Spitkovsky Reference Tchekhovskoy, Philippov and Spitkovsky2016). This magnetic field is unidirectional in the polar region $\unicode[STIX]{x1D703}<\unicode[STIX]{x1D703}_{p}$ but in the equatorial region it comes in the form of stripes with opposite magnetic directions. In the monopole model, $\unicode[STIX]{x1D703}_{p}=\unicode[STIX]{x03C0}/2-\unicode[STIX]{x1D6FC}$ , where $\unicode[STIX]{x1D6FC}$ is the angle between the magnetic and rotational axes of the pulsar but in the dipolar model the poloidal field lines are not straight and this relation may no longer hold.

For $r\gg r_{\text{lc}}$ the wind becomes relativistic and radial and its electric field is dominated by the polar component

(2.2) $$\begin{eqnarray}E_{\unicode[STIX]{x1D703}}=(v/c)B_{\unicode[STIX]{x1D719}}=B_{\unicode[STIX]{x1D719}}.\end{eqnarray}$$

The corresponding potential drop across the flow in the polar zone is

(2.3) $$\begin{eqnarray}\unicode[STIX]{x1D6F7}(\unicode[STIX]{x1D703})=\int _{0}^{\unicode[STIX]{x1D703}}E_{\unicode[STIX]{x1D703}}r\,\text{d}\unicode[STIX]{x1D703}=\unicode[STIX]{x1D6F7}_{0}f(\unicode[STIX]{x1D703}),\end{eqnarray}$$

where

(2.4) $$\begin{eqnarray}\unicode[STIX]{x1D6F7}_{0}=B_{\text{lc}}r_{\text{lc}}\end{eqnarray}$$

and

(2.5) $$\begin{eqnarray}f(\unicode[STIX]{x1D703})=\left\{\begin{array}{@{}ll@{}}1-\cos (\unicode[STIX]{x1D703}) & \text{for }n=1\\ 0.5(\unicode[STIX]{x1D703}-0.5\sin 2\unicode[STIX]{x1D703}) & \text{for }n=2.\end{array}\right.\end{eqnarray}$$

In the striped wind zone one has to take into account the alternation of magnetic field direction and use the stripe-averaged magnetic field, which is smaller than that in (2.1) (see Komissarov Reference Komissarov2013). As a result, inside the striped wind zone the potential grows much slower (see figure 1).

The corresponding Poynting vector in the wind is

(2.6) $$\begin{eqnarray}\boldsymbol{S}_{P}=\frac{1}{4\unicode[STIX]{x03C0}}E_{\unicode[STIX]{x1D703}}B_{\unicode[STIX]{x1D719}}c=\frac{c}{4\unicode[STIX]{x03C0}}\left(\frac{\unicode[STIX]{x1D6F7}_{0}}{r}\right)^{2}\sin ^{2n}\unicode[STIX]{x1D703},\end{eqnarray}$$

which leads to the total wind power

(2.7) $$\begin{eqnarray}L_{w}=a_{n}c\,\unicode[STIX]{x1D6F7}_{0}^{2},\end{eqnarray}$$

where $a_{1}=2/3$ and $a_{2}=8/15$ . Thus one can find $\unicode[STIX]{x1D6F7}_{0}$ given the measurement of $L_{w}$ as the spindown power of the pulsar. For the Crab pulsar $L_{w}\approx 4\times 10^{38}~\text{erg}~\text{s}^{-1}$ which leads to $\unicode[STIX]{x1D6F7}_{0}\approx 42~\text{PeV}$ , the corresponding electron Lorentz factor

(2.8) $$\begin{eqnarray}\unicode[STIX]{x1D6FE}_{\text{max}}\approx 8.2\times 10^{10}\end{eqnarray}$$

and the wind magnetic field

(2.9) $$\begin{eqnarray}B_{\unicode[STIX]{x1D719}}\approx 5\times 10^{-2}r_{,\text{ld}}^{-1}\sin ^{n}\unicode[STIX]{x1D703}\,\text{G},\end{eqnarray}$$

where $r_{,\text{ld}}$ is the distance measured in light days.

At the light cylinder, the Goldreich–Julian number density of charged particles

(2.10) $$\begin{eqnarray}n_{GJ}=\frac{1}{4\unicode[STIX]{x03C0}e}\,\text{div}\,E\approx \frac{1}{4\unicode[STIX]{x03C0}e}\frac{B_{\text{lc}}}{r_{\text{lc}}}=\frac{1}{4\unicode[STIX]{x03C0}e}\frac{\unicode[STIX]{x1D6F7}_{0}}{r_{\text{lc}}^{2}},\end{eqnarray}$$

whereas the actual density is governed by the pair-production rate $n=\unicode[STIX]{x1D706}n_{GJ}$ , where $\unicode[STIX]{x1D706}$ is called the multiplicity factor. The corresponding total $e_{\pm }$ -particle flux of the wind

(2.11) $$\begin{eqnarray}{\dot{N}}\approx 4\unicode[STIX]{x03C0}r_{\text{lc}}^{2}nc=\frac{\unicode[STIX]{x1D6F7}_{0}c}{e}\approx 8.7\times 10^{37}\unicode[STIX]{x1D706}_{4}~\text{s}^{-1},\end{eqnarray}$$

where $\unicode[STIX]{x1D706}_{4}=10^{-4}\unicode[STIX]{x1D706}$ . This corresponds to the particle density in the wind

(2.12) $$\begin{eqnarray}n\approx 3.4\times 10^{-5}\unicode[STIX]{x1D706}_{4}r_{,\text{ld}}^{-2}~\text{cm}^{-3}.\end{eqnarray}$$

Along the particle flux tube of a steady-state wind, the particle flux $S_{\pm }$ and the total energy flux $S_{\text{tot}}=S_{P}+m_{e}c^{2}\unicode[STIX]{x1D6FE}_{w}hS_{\pm }$ , where $\unicode[STIX]{x1D6FE}_{w}$ is the wind Lorentz factor and $h=w/{\tilde{n}}m_{e}c^{2}$ is the relativistic enthalpy per unit particle rest mass, are conserved. This leads to the integration constant along streamlines

(2.13) $$\begin{eqnarray}\unicode[STIX]{x1D70E}_{M}=\frac{S_{\text{tot}}}{S_{\pm }m_{e}c^{2}}=\unicode[STIX]{x1D6FE}_{w}h(1+\unicode[STIX]{x1D70E}),\end{eqnarray}$$

where $\unicode[STIX]{x1D70E}=S_{P}/S_{p}$ is the ratio of the Poynting flux and the particle energy flux $S_{p}=S_{\pm }m_{e}c^{2}h\unicode[STIX]{x1D6FE}_{w}$ . In the theory of pulsar winds, $\unicode[STIX]{x1D70E}_{M}$ is known as Michel’s $\unicode[STIX]{x1D70E}$ (Michel Reference Michel1969) but in the theory of magnetically dominated jets it is simply denoted as $\unicode[STIX]{x1D707}$ (e.g. Komissarov et al. Reference Komissarov, Barkov, Vlahakis and Königl2007). Replacing $S_{P}$ and $S_{\pm }$ with their densities, $B_{\unicode[STIX]{x1D719}}^{2}c/4\unicode[STIX]{x03C0}$ and $nc$ respectively, we obtain

(2.14) $$\begin{eqnarray}\unicode[STIX]{x1D70E}=\frac{B_{\unicode[STIX]{x1D719}}^{2}}{4\unicode[STIX]{x03C0}m_{e}c^{2}hn\unicode[STIX]{x1D6FE}_{w}}=\frac{\tilde{B}_{\unicode[STIX]{x1D719}}^{2}}{4\unicode[STIX]{x03C0}m_{e}c^{2}h{\tilde{n}}},\end{eqnarray}$$

where $\tilde{B}$ and ${\tilde{n}}$ are measured in the wind frame. Using the parameters at the light cylinder to estimate $\unicode[STIX]{x1D70E}_{M}$ , we find

(2.15) $$\begin{eqnarray}\unicode[STIX]{x1D70E}_{M}=\frac{\unicode[STIX]{x1D6F7}_{0}e}{m_{e}c^{2}\unicode[STIX]{x1D706}}\approx 8.2\times 10^{6}\unicode[STIX]{x1D706}_{4}^{-1}.\end{eqnarray}$$

For a cold wind $h=1$ and hence

(2.16) $$\begin{eqnarray}\unicode[STIX]{x1D6FE}_{w}=\frac{\unicode[STIX]{x1D70E}_{M}}{1+\unicode[STIX]{x1D70E}}.\end{eqnarray}$$

Thus the total conversion of the Poynting flux into the kinetic energy of the bulk motion yields the wind speed

(2.17) $$\begin{eqnarray}\unicode[STIX]{x1D6FE}_{w}^{\text{max}}=\unicode[STIX]{x1D70E}_{M}.\end{eqnarray}$$

However in the limit of ideal MHD, reaching this asymptotic value is quite problematic. Indeed, the wind must first become super-fast magnetosonic. At the transonic point

(2.18) $$\begin{eqnarray}\unicode[STIX]{x1D6FE}_{w}=\sqrt{\unicode[STIX]{x1D70E}}\end{eqnarray}$$

(e.g. Komissarov Reference Komissarov2012). Combining this result with (2.16) we find that at the fast magnetosonic point

(2.19a,b ) $$\begin{eqnarray}\unicode[STIX]{x1D6FE}_{w,f}=\unicode[STIX]{x1D70E}_{M}^{1/3}=200\,\unicode[STIX]{x1D706}_{4}^{-1/3}\quad \text{and}\quad \unicode[STIX]{x1D70E}_{f}=4\times 10^{4}\unicode[STIX]{x1D706}_{4}^{-2/3}.\end{eqnarray}$$

Beyond this point the wind becomes causally disconnected, which strongly reduces the efficiency of the ideal collimation–acceleration mechanism (e.g. Beskin, Kuznetsova & Rafikov Reference Beskin, Kuznetsova and Rafikov1998; Komissarov et al. Reference Komissarov, Vlahakis, Königl and Barkov2009). Thus we do not expect the wind’s Lorentz factor to exceed the value given in (2.20a,b ) by more than a factor of ten. Using the fast magnetosonic Mach number $M\approx \unicode[STIX]{x1D6E4}_{w}/\sqrt{\unicode[STIX]{x1D70E}}$ , the wind parameters at the termination shock can be estimated as

(2.20a,b ) $$\begin{eqnarray}\displaystyle \unicode[STIX]{x1D6FE}_{w,u}=M^{2/3}\unicode[STIX]{x1D70E}_{M}^{1/3}=200\,M^{2/3}\unicode[STIX]{x1D706}_{4}^{-1/3}\quad \text{and}\quad \unicode[STIX]{x1D70E}_{u}=(\unicode[STIX]{x1D70E}_{M}/M)^{2/3}=4\times 10^{4}M^{-2/3}\unicode[STIX]{x1D706}_{4}^{-2/3}. & & \displaystyle \nonumber\\ \displaystyle & & \displaystyle\end{eqnarray}$$

Strictly speaking, these results do not apply to the striped wind zone where the magnetic dissipation of the stripes may come into play and allow some additional acceleration of the flow. Even in the polar zone things can be more complicated if a substantial fraction of the energy is carried away by fast magnetosonic waves emitted into the zone because of the misalignment of the magnetic and rotational axes of the Crab pulsar. Lyubarsky (Reference Lyubarsky2003a ) argues that for $\unicode[STIX]{x1D706}\gtrsim 10^{7}$ these waves steepen into shocks and dissipate before reaching the termination shock and the generated heat is quickly converted into the bulk kinetic energy of the wind particles. This process would set an upper limit on the wind magnetisation $\unicode[STIX]{x1D70E}$ at the termination shock equal to the inverse fraction of the energy flux carried by the waves and the Poynting flux of the wind. For example, if the waves contribute only one per cent to the wind energy budget then $\unicode[STIX]{x1D70E}_{u}<100$ . For smaller multiplicity, the waves would dissipate only after crossing the termination shock and the dissipation rate increases sharply in the sub-fast magnetosonic regime (Lyubarsky Reference Lyubarsky2003a ).

As to the magnetic dissipation in the striped wind zone, its rate can be estimated as follows (Arons Reference Arons2012). In the frame of Earth, the typical thickness of the stripes is $l_{s}=cP/2$ , where $P$ is the pulsar period. In the wind frame, it is $l_{s}^{\prime }=\unicode[STIX]{x1D6FE}_{w}cP/2$ and hence the shortest time of stripe dissipation is $t_{s}^{\prime }=\unicode[STIX]{x1D6FE}_{w}P/2$ . In the Earth frame, the time dilation yields $t_{s}=\unicode[STIX]{x1D6FE}_{w}^{2}P/2$ . Once the dissipation is completed $\unicode[STIX]{x1D6FE}_{w}=\unicode[STIX]{x1D6FE}_{w}^{\text{max}}$ , which yields the distance to this point as

(2.21) $$\begin{eqnarray}r_{d}={\textstyle \frac{1}{2}}\unicode[STIX]{x1D70E}_{M}^{2}Pc\approx 1.3\times 10^{7}\unicode[STIX]{x1D706}_{4}^{-2}~\text{light days}.\end{eqnarray}$$

Unless $\unicode[STIX]{x1D706}_{4}>300$ this is above the radius of the wind termination shock $r_{\text{ts}}\approx 170$ light days (if we identify it with the inner X-ray ring of the Crab Nebula) and the stripes survive all the way up to the termination shock, where $\unicode[STIX]{x1D6FE}_{w}\ll \unicode[STIX]{x1D70E}_{M}$ . However, Lyubarsky (Reference Lyubarsky2003b ) has shown that in this case the stripes will dissipate at the termination shock and the post-shock state will be the same as in the case where the stripe dissipation is completed in the wind (see also Sironi & Spitkovsky Reference Sironi and Spitkovsky2012). In both these cases, the magnetisation of the post-shock state is determined by the stripe-averaged magnetic field, which vanishes only in the equatorial plane but is the same as the wind field at the boundary with the polar zone. In the transition between the polar and striped zones the magnetisation parameter drops rapidly from the high value of the polar zone (see (2.20)) to $\unicode[STIX]{x1D70E}\approx 1$ and then keeps decreasing with the polar angle, but at a lower rate (Komissarov Reference Komissarov2013).

2.2 The nebula

The nebula can be described as a prolate spheroid with the minor axis $a\approx 9.5\,\text{ly}$ and the major axis $b\approx 14\,\text{ly}$ (Hester Reference Hester2008). The most spectacular structure of the nebula is its network of thermal filaments. These are associated with the stellar material ejected during the supernova explosion which left behind the Crab pulsar. The non-thermal emission is more diffusive but it also exhibits many features on different length scales. In particular, the radio maps of the nebula show a filamentary structure which is almost identical to network of optical thermal filaments.

The X-ray images of the inner nebula reveal an equatorial torus-like structure of the radius of approximately 500 light days and a polar jet of a similar size. Inside the torus, there is an oval structure which looks like a beaded neckless around the pulsar, which is called the ‘X-ray inner ring’. Its radius is $r_{ir}\approx 170$ light days. The ring is often interpreted as the termination shock of the pulsar wind. Since the pulsar wind is anisotropic, its termination shock is not spherical and it is most extended in the equatorial plane, where it become transverse and forms the so-called Mach belt (Komissarov & Lyubarsky Reference Komissarov and Lyubarsky2004). On purely geometrical grounds, one is tempted to identify the ring with the Mach belt, but the nature of its emission is not clear. Hester et al. (Reference Hester, Mori, Burrows, Gallagher, Graham, Halverson, Kader, Michel and Scowen2002) state that the beads (knots) of the ring ‘form, brighten, fade, dissipate, move about, and occasionally undergo outbursts, giving rise to expanding clouds of nebulosity’, though no quantitative parameters of this variability are given. Just the naked eye inspection of the Chandra images indicates that their size is approximately one arcsec, which is about the angular resolution of Chandra, corresponding to the linear scale ${\approx}10~\text{ld}$ . One cannot exclude that some knots of the ring are even smaller. The ring has no clear identification with features observed in the optical band with HST. However, the HST images reveal a bright compact knot at approximately 0.65 arcsec from the pulsar. This knot is resolved and has the size similar to that of the X-ray knots. The knot is not seen in the radio and X-ray bands (Rudy et al. Reference Rudy, Horns, DeLuca, Kolodziejczak, Tennant, Yuan, Buehler, Arons, Blandford and Caraveo2015). In addition to these features, there are also fine bright arc-like filaments (‘wisps’) originating in the vicinity of the inner ring and moving away from it with speeds up to $0.7\,c.$ These are seen in the radio, optical and X-ray bands, though the wisps seen in different bands are sometimes displaced relative to each other or do not even have a counterpart in the other band (Bietenholz et al. Reference Bietenholz, Hester, Frail and Bartel2004; Hester Reference Hester2008).

The observations of both the synchrotron and inverse-Compton emission of the same population of relativistic leptons allows a more-or-less reliable measurement of the nebula energetics. Using a one-zone model of the nebula, Meyer, Horns & Zechlin (Reference Meyer, Horns and Zechlin2010) arrive at a mean magnetic field strength of $B\approx 0.12~\text{mG}$ , with the corresponding magnetic energy density $e_{m}\approx 6\times 10^{-10}~\text{erg}~\text{cm}^{-3}$ . The lepton population is split into the low-energy sub-population responsible for the radio synchrotron emission and the high-energy population responsible for the optical through the X-ray to the gamma-ray synchrotron emission. The radio leptons have a hard spectrum,

(2.22) $$\begin{eqnarray}\text{d}n/\text{d}\unicode[STIX]{x1D6FE}\propto \unicode[STIX]{x1D6FE}^{-p_{r}},\quad p_{r}=1.6,~20<\unicode[STIX]{x1D6FE}<2\times 10^{5}\end{eqnarray}$$

and the optical-X-ray leptons have a soft spectrum,

(2.23) $$\begin{eqnarray}\text{d}n/\text{d}\unicode[STIX]{x1D6FE}\propto \unicode[STIX]{x1D6FE}^{-p_{o}},\quad p_{o}=3.2,~4.4\times 10^{5}<\unicode[STIX]{x1D6FE}<3\times 10^{8}.\end{eqnarray}$$

The total energies of the low- and high-energy populations are comparable, $E_{l}\approx 3.1\times 10^{48}~\text{erg}$ and $E_{h}\approx 2.3\times 10^{48}~\text{erg}$ respectively. Given the total volume of the nebula, $V\approx 5.7\times 10^{56}~\text{cm}^{3}$ , the corresponding total energy density of the leptons is $e_{p}\approx 9.4\times 10^{-9}~\text{erg}~\text{cm}^{-3}$ and the total pressure (leptons $+$ magnetic field) is $p_{\text{tot}}\approx 3.6\times 10^{-9}$ . This total pressure can be balanced by the pressure of magnetic field $B=0.3~\text{mG}$ , which is close to the standard equipartition value obtained for the Crab Nebula (Hillas et al. Reference Hillas, Akerlof, Biller, Buckley, Carter-Lewis, Catanese, Cawley, Fegan, Finley and Gaidos1998).

The synchrotron cooling of relativistic electrons leads to a break in the particle spectrum at

(2.24) $$\begin{eqnarray}\unicode[STIX]{x1D6FE}_{b}=(c_{2}m_{e}c^{2}B_{\bot }^{2}t_{n})^{-1}=1.6\times 10^{6}B_{-4}^{-2}t_{n,3}^{-1},\end{eqnarray}$$

where $t_{n,3}$ is the nebula age in units of $10^{3}\,$ years and $c_{2}=2e^{4}/3m_{e}^{4}c^{7}=0.00237$ in cgs units (Pacholczyk Reference Pacholczyk1970). The corresponding energy of synchrotron photons is

(2.25) $$\begin{eqnarray}{\mathcal{E}}_{\unicode[STIX]{x1D708}}=c_{1}B_{\bot }{\mathcal{E}}^{2}=4.4B_{-4}^{-3}t_{n,3}^{-2}~\text{eV},\end{eqnarray}$$

where $c_{1}=3eh/4\unicode[STIX]{x03C0}m_{e}^{3}c^{5}=4.14\times 10^{-8}$ cgs units (Pacholczyk Reference Pacholczyk1970). Thus the spectrum of the optical-X-ray leptons is subject to steepening with the spectral index increasing by unity compared to that of the injection spectrum (Kardashev Reference Kardashev1962), which has to be $a_{i}=2.2$ .

Since the nebula has a rich structure on different scales, it is natural to expect strong variations of these parameters from feature to feature. Such variations are observed in the MHD simulations of pulsar wind nebulae (PWN) (Porth et al. Reference Porth, Komissarov and Keppens2013; Olmi et al. Reference Olmi, Del Zanna, Amato, Bucciantini and Mignone2016), with the total pressure exceeding the mean value by up to 5 times at the base of its polar jet and the magnetic field near the termination shock exceeding its mean value by an order of magnitude. Based on these results one may expect the magnetic field in the very inner parts of the Crab Nebula to be as high as one mG.

While it is widely accepted that the high-energy electrons are constantly supplied into the nebula by the pulsar wind, the nature of its radio electrons is still debated. The observed continuity of the synchrotron spectrum at the transition between the populations has been used to argue that the radio electrons are also supplied by the wind (Bucciantini, Arons & Amato Reference Bucciantini, Arons and Amato2011; Arons Reference Arons2012). However, this requires a much more prolific particle-production process in pulsar magnetospheres than the one currently accepted. Given the results by Meyer et al. (Reference Meyer, Horns and Zechlin2010), the high- and low-energy populations contain $N_{h}\approx 3.8\times 10^{48}$ and $N_{l}\approx 3.4\times 10^{51}$ leptons respectively. The corresponding number densities are $n_{h}\approx 6.7\times 10^{-9}~\text{cm}^{-3}$ and $n_{l}\approx 6\times 10^{-6}~\text{cm}^{-3}$ . Given the known age of the nebula $t_{n}=963~\text{yr}$ , the corresponding mean injection rates are ${\dot{N}}_{h}\approx 1.2\times 10^{38}~\text{s}^{-1}$ and ${\dot{N}}_{l}\approx 1.1\times 10^{41}~\text{s}^{-1}$ respectively. Comparing these data with (2.11) one finds that if only the high-energy population leptons are injected by the wind then the required multiplicity parameter $\unicode[STIX]{x1D706}\approx 10^{4}$ which is consistent with the theory of pair production in pulsar magnetospheres, whereas if the radio leptons are injected as well then $\unicode[STIX]{x1D706}\approx 10^{7}$ (Shklovsky Reference Shklovsky1970), which requires a grand revision of this theory (Arons Reference Arons2012). Using (2.12) for the mean number density of particles in the wind one can check the consistency of these estimates of $\unicode[STIX]{x1D706}$ with the size of the termination shock. Indeed, at the termination shock the number density of wind particles as measured in the laboratory (Earth’s) frame does not change significantly and should match the observed densities in the nebula. Matching it with $n_{l}$ gives the distance

(2.26) $$\begin{eqnarray}r\approx 70\,\unicode[STIX]{x1D706}_{7}^{1/2}~\text{ld},\end{eqnarray}$$

whereas matching with $n_{h}$ gives

(2.27) $$\begin{eqnarray}r\approx 70\,\unicode[STIX]{x1D706}_{4}^{1/2}~\text{ld}.\end{eqnarray}$$

These results are consistent with the observations, which give the equatorial radius of the termination shock $r_{\text{ts}}\approx 170~\text{ld}$ .

Using $\unicode[STIX]{x1D706}=10^{7}$ , we find that upstream of the polar zone of the termination shock $\unicode[STIX]{x1D6FE}_{u}\approx 20$ and $\unicode[STIX]{x1D70E}_{u}\approx 400$ . This value of $\unicode[STIX]{x1D6FE}_{u}$ is approximately the same as the minimum value $\unicode[STIX]{x1D6FE}_{r,\text{min}}$ of Crab’s radio leptons (Meyer et al. Reference Meyer, Horns and Zechlin2010), which provides some supports the hypothesis that the radio electrons are supplied by the wind. In the absence of stripes one can rule out the shock acceleration mechanism and hence some other acceleration process has to operate downstream of the shock. According to the results derived in Komissarov & Lyutikov (Reference Komissarov and Lyutikov2011), downstream of a highly magnetised shock

(2.28a,b ) $$\begin{eqnarray}\unicode[STIX]{x1D6FE}_{d}\approx \sqrt{\unicode[STIX]{x1D70E}_{u}}/\sin \unicode[STIX]{x1D6FF}\quad \text{and}\quad \unicode[STIX]{x1D70E}_{d}\approx \sqrt{\unicode[STIX]{x1D70E}_{u}}\unicode[STIX]{x1D6FE}_{u}\sin \unicode[STIX]{x1D6FF},\end{eqnarray}$$

where $\unicode[STIX]{x1D6FF}$ is the angle between the upstream velocity and the shock plane. If the effect of magnetosonic waves is ignored then the upstream parameters can be estimated using (2.20), which yields $\unicode[STIX]{x1D70E}_{u}\approx 400M^{-2/3}\text{ and }\unicode[STIX]{x1D6FE}_{u}\approx 20M^{2/3}$ and hence

(2.29a,b ) $$\begin{eqnarray}\unicode[STIX]{x1D6FE}_{d}=20M^{-1/3}/\sin \,\unicode[STIX]{x1D6FF}\quad \text{and}\quad \unicode[STIX]{x1D70E}_{d}=400M^{1/3}\sin \unicode[STIX]{x1D6FF}.\end{eqnarray}$$

We note that $\unicode[STIX]{x1D6FE}_{d}$ and $\unicode[STIX]{x1D70E}_{d}$ depend rather weakly on the wind Mach number, a parameter which we do not know well but do not expect to exceed $M\approx 10$ by much. The low wind Lorentz factor obtained in this model is consistent with the radio electrons, which also have low Lorentz factors, being supplied by the wind. The high upstream magnetisation rules out particle acceleration at the termination shock and shifts the focus to the processes occurring downstream of the shock instead. The hard spectrum of radio electrons hints at the magnetic reconnection mechanism instead but for $\unicode[STIX]{x1D70E}_{d}>100$ this mechanism yields spectra (Sironi & Spitkovsky Reference Sironi and Spitkovsky2014; Guo et al. Reference Guo, Liu, Daughton and Li2015) which are even harder, with $p\approx 1$ , than that of the radio electronsFootnote ¹ . As we have discussed earlier, for this value of the multiplicity parameter the fast magnetosonic waves produced by the pulsar rotation are expected to dissipate upstream of the termination shock (Lyubarsky Reference Lyubarsky2003a ) and hence to reduce $\unicode[STIX]{x1D70E}_{u}$ and hence $\unicode[STIX]{x1D70E}_{d}$ , provided they carry a fraction of the total energy flux that exceeds $1/\unicode[STIX]{x1D70E}_{u}$ .

For $\unicode[STIX]{x1D706}=10^{7}$ the stripes of the pulsar wind get erased before the termination shock, see (2.21), and hence in the equatorial zone the wind entering the shock should have a unidirectional magnetic field, with the magnetisation $\unicode[STIX]{x1D70E}<1$ (Komissarov Reference Komissarov2013) and the Lorentz factor $\unicode[STIX]{x1D6FE}_{w}=\unicode[STIX]{x1D70E}_{M}\approx 8\times 10^{3}$ . Because these parameters are very different from those of the polar zone one would expect a different particle acceleration regime (or mechanism) as well, leading to different spectral properties of the accelerated particles. This hints a possibility of explaining the existence of two different non-thermal populations of the Crab Nebula – the radio electrons are supplied via the polar zone whereas the high-energy electrons via the equatorial striped zone (cf. Porth et al. Reference Porth, Komissarov and Keppens2013; Olmi et al. Reference Olmi, Del Zanna, Amato and Bucciantini2015). However the details are not clear. For example, according to the PIC simulations of relativistic shocks in a pair plasma, the magnetisation has to be as low as $\unicode[STIX]{x1D70E}=10^{-3}$ for the traditional Fermi mechanism to operate. Such a low magnetisation is expected only in a very narrow equatorial part of the equatorial zone, with the half-opening angle approximately one degree for the magnetic inclination angle $\unicode[STIX]{x1D6FC}=45^{\circ }$ (Komissarov Reference Komissarov2013). For the rest of it, the shock is expected only to thermalise the wind particles leading to a Maxwellian-like spectrum with $\unicode[STIX]{x1D6FE}_{t}\approx \unicode[STIX]{x1D6FE}_{w}\approx 8\times 10^{3}$ . This would imply strong emission at wavelengths of 1–10 mm associated with the termination shock, which is not seen (Gomez et al. Reference Gomez, Krause, Barlow, Swinyard, Owen, Clark, Matsuura, Gomez, Rho and Besel2012). Finally, the very idea that the two non-thermal populations could be supplied via different sections of the pulsar wind with very different physical parameters, and presumably shaped via different acceleration mechanisms, reintroduces the difficulty of explaining the continuity of the spectrum of the nebula emission. Moreover, if the radio electrons are supplied via the polar zone then the polar zone should carry about a thousand times more particles than the striped wind zone, which is very hard to explain.

For the more traditional multiplicity $\unicode[STIX]{x1D706}=10^{4}$ the radio electrons cannot be supplied by the wind and must have a different origin. For example, the observed filamentary structure of the synchrotron radio emission suggests that it may be closely related to the supernova ejecta shocked by the pulsar wind (Komissarov Reference Komissarov2013). However the hard spectral index of the radio emission disfavours the diffusive shock acceleration mechanism. As to the termination shock, (2.20) shows that in the polar zone the wind parameters upstream of the shock are $\unicode[STIX]{x1D6FE}_{u}\approx 10^{3}M^{2/3}$ and $\unicode[STIX]{x1D70E}_{u}\approx 10^{4}M^{-2/3}$ whereas (2.28) shows that downstream of it

(2.30a,b ) $$\begin{eqnarray}\unicode[STIX]{x1D6FE}_{d}\approx 10^{3}M^{-1/3}/\sin \unicode[STIX]{x1D6FF}\quad \text{and}\quad \unicode[STIX]{x1D70E}_{d}\approx 10^{4}M^{1/3}\sin \unicode[STIX]{x1D6FF}.\end{eqnarray}$$

These parameters are even more extreme that for $\unicode[STIX]{x1D706}=10^{7}$ . Moreover, in this case the fast magnetosonic waves cannot dissipate upstream of the termination shock and hence cannot significantly reduce the wind’s $\unicode[STIX]{x1D70E}$ .

As to the striped wind zone, (2.21) tells us that in this case the stripes survive all the way to the termination shock. As they dissipate at the shock, the spectrum of energised particles depends on the value of the parameter

(2.31) $$\begin{eqnarray}\unicode[STIX]{x1D701}=4\unicode[STIX]{x03C0}\unicode[STIX]{x1D706}\frac{r_{\text{lc}}}{r_{\text{ts}}}\end{eqnarray}$$

(Sironi & Spitkovsky Reference Sironi and Spitkovsky2011). For $\unicode[STIX]{x1D701}<10$ it is a relativistic Maxwellian-like spectrum whereas a broad power-law spectrum is produced only for $\unicode[STIX]{x1D701}>100$ . For the Crab pulsar $r_{\text{lc}}\approx 6\times 10^{-8}\,$ ld and $r_{\text{ts}}\approx 170\,$ ld so

(2.32) $$\begin{eqnarray}\unicode[STIX]{x1D701}\approx 4\times 10^{-4}\unicode[STIX]{x1D706}_{4}.\end{eqnarray}$$

Thus we expect a Maxwellian-like spectrum peaking at $\unicode[STIX]{x1D6FE}_{p}\approx \unicode[STIX]{x1D70E}_{M}\approx 8\times 10^{6}\unicode[STIX]{x1D706}_{4}^{-1}$ . The corresponding synchrotron frequency

(2.33) $$\begin{eqnarray}\unicode[STIX]{x1D708}=5\times 10^{16}B_{-3}\unicode[STIX]{x1D706}_{4}^{-2}~\text{Hz}.\end{eqnarray}$$

The inner knot of the Crab Nebula is currently associated with the emission of shocked striped wind (Lyutikov et al. Reference Lyutikov, Komissarov and Porth2016). The knot is seen only in the optical and near IR bands. Given that the theory predicts a Doppler shift of the knot emission by a factor of a few towards higher energies, this requires $\unicode[STIX]{x1D706}\approx 10^{5}$ , which seems quite realistic. More observational data on the knot emission are required to clarify the nature of its emission.

The above analysis shows a glaring conflict between the observations and the current understanding of relativistic magnetised shocks as particle accelerators. Either the current shock theory based on PIC simulations is completely wrong or the observed spectrum is shaped by other processes. These could be a second-order Fermi mechanism of particle acceleration by plasma turbulence or magnetic reconnection, both occurring inside the nebula.

4.1 Peak energy

The fact that the peak energy of the Crab flares is located about the radiation-reaction limit for the synchrotron mechanism suggests that the observed radiation is produced in this regime, whereby the energy gains due to the acceleration mechanism are balanced by the energy losses via the synchrotron radiation. Accelerated particles are confined to the current sheet where the magnetic field is weaker compared to the reconnecting magnetic field $B$ outside. Moreover only the component of magnetic field perpendicular to the particle momentum contributes to the synchrotron energy loss. Denoting this component as $B_{a,\bot }=\unicode[STIX]{x1D705}B$ , where $0<\unicode[STIX]{x1D705}<1$ , we can write the energy loss rate by a particle of energy ${\mathcal{E}}$ as

(4.1) $$\begin{eqnarray}\frac{\text{d}{\mathcal{E}}}{\text{d}t}=-c_{2}B_{a,\bot }^{2}{\mathcal{E}}^{2},\end{eqnarray}$$

where $c_{2}=2e^{4}/3m_{e}^{4}c^{7}=0.00237$ in cgs units. The energy gain is via the electrostatic acceleration by the reconnection electric field

(4.2) $$\begin{eqnarray}\frac{\text{d}{\mathcal{E}}}{\text{d}t}=eE_{a}c,\end{eqnarray}$$

where $E_{a}=\unicode[STIX]{x1D702}B$ is the electric field and $0<\unicode[STIX]{x1D702}<1$ is the reconnection rate. The balance between the two yields the electron energy

(4.3) $$\begin{eqnarray}{\mathcal{E}}_{\text{max}}=\frac{1}{\unicode[STIX]{x1D705}}\left(\frac{\unicode[STIX]{x1D702}ec}{c_{2}B}\right)=1.53\frac{\unicode[STIX]{x1D702}^{1/2}}{\unicode[STIX]{x1D705}}B_{-3}^{-1/2}~\text{PeV},\end{eqnarray}$$

where $B_{-3}$ is the magnetic field in mG. The corresponding Lorentz factor of flare-producing leptons is

(4.4) $$\begin{eqnarray}\unicode[STIX]{x1D6FE}=3\times 10^{9}\frac{\unicode[STIX]{x1D702}^{1/2}}{\unicode[STIX]{x1D705}}B_{-3}^{-1/2}.\end{eqnarray}$$

This is approximately 30 times smaller than the maximal possible Lorentz factor based on the electric potential drop across the pulsar wind given by (2.8), thus supporting the possibility of electrostatic acceleration. Using (2.3)–(2.5) one can estimate the minimum size of the polar wind zone delivering the required potential drop. For $\unicode[STIX]{x1D6FE}=3\times 10^{9}$ this is

(4.5) $$\begin{eqnarray}\unicode[STIX]{x1D703}_{p}^{\text{min}}=\left\{\begin{array}{@{}ll@{}}15^{\circ } & \text{for }n=1;\\ 27^{\circ } & \text{for }n=2.\end{array}\right.\end{eqnarray}$$

These are reasonably small compared to $\unicode[STIX]{x1D703}_{p}\approx 45^{\circ }$ estimated for the Crab pulsar (Harding et al. Reference Harding, Stern, Dyks and Frackowiak2008; Lyutikov et al. Reference Lyutikov, Komissarov and Porth2016).

The energy of synchrotron photons emitted by the particle of energy ${\mathcal{E}}$ in the magnetic field $B$ is

(4.6) $$\begin{eqnarray}{\mathcal{E}}_{\unicode[STIX]{x1D708}}=c_{1}B_{\bot }{\mathcal{E}}^{2},\end{eqnarray}$$

where $c_{1}=3eh/4\unicode[STIX]{x03C0}m_{e}^{3}c^{5}=4.14\times 10^{-8}$ cgs units (Pacholczyk Reference Pacholczyk1970). It is observed in the simulations (Cerutti et al. Reference Cerutti, Werner, Uzdensky and Begelman2013) that the highest-energy photons are emitted when the accelerated particle is ejected from the current layer into the reconnecting magnetic field. With this understanding of $B$ in (4.6) we find the maximum energy of photons emitted from the reconnection region as

(4.7) $$\begin{eqnarray}{\mathcal{E}}_{\unicode[STIX]{x1D708},\text{max}}=\frac{c_{1}}{c_{2}}\frac{\unicode[STIX]{x1D702}}{\unicode[STIX]{x1D705}^{2}}ec=158\left(\frac{\unicode[STIX]{x1D702}}{\unicode[STIX]{x1D705}^{2}}\right)~\text{MeV}.\end{eqnarray}$$

In papers I and II we show that for the reconnection driven by macroscopic stresses in a highly magnetised plasma, the reconnection speed approaches the speed of light and hence $\unicode[STIX]{x1D702}\rightarrow 1$ (papers I and II). Thus $\unicode[STIX]{x1D705}\simeq 0.5$ is sufficient to explain the peak photon energy of Crab flares.

The acceleration time required to reach ${\mathcal{E}}_{\text{max}}$ is

(4.8) $$\begin{eqnarray}t_{\text{acc}}=\frac{{\mathcal{E}}_{\text{max}}}{\unicode[STIX]{x1D702}eBc}=2\frac{1}{\unicode[STIX]{x1D705}\unicode[STIX]{x1D702}^{1/2}}B_{-3}^{-3/2}~\text{days},\end{eqnarray}$$

whereas the gyration period of an electron with this energy in the reconnecting magnetic field is

(4.9) $$\begin{eqnarray}t_{g}=2\unicode[STIX]{x03C0}\unicode[STIX]{x1D702}t_{\text{acc}}.\end{eqnarray}$$

The corresponding gyration radius is

(4.10) $$\begin{eqnarray}r_{g}=\unicode[STIX]{x1D702}ct_{\text{acc}}.\end{eqnarray}$$

4.2 Size and energetics of the flaring site

The acceleration time sets the lower limit on the size of the ‘accelerator’ (the current layer)

(4.11) $$\begin{eqnarray}l_{\text{min}}=ct_{\text{acc}}=2\frac{1}{\unicode[STIX]{x1D705}\unicode[STIX]{x1D702}^{1/2}}B_{-3}^{-3/2}~\text{light days}.\end{eqnarray}$$

Interestingly, this differs from the gyration radius of the electrons with energy ${\mathcal{E}}_{\text{max}}$ only by the factor $1/\unicode[STIX]{x1D702}$ . Assuming that the reconnection region is not moving with relativistic speed, the actual size of the active region should fall within the bounds

(4.12) $$\begin{eqnarray}l_{\text{min}}<l<ct_{f},\end{eqnarray}$$

where $t_{f}$ is the flare duration. Because $t_{f}$ is of the order of a few days, the last two equations imply that the reconnecting magnetic field cannot be significantly below $1\,$ mG. As we discussed in § 2.2, $B\approx 1~\text{mG}$ is too high for most of the nebula but it is not unreasonable for the very inner part of the Crab Nebula. In fact, direct application of (2.9) shows that the wind magnetic field has such a strength at the distance $r=50\sin ^{n}\unicode[STIX]{x1D703}$ light days, which is significantly smaller compared to the equatorial radius of the termination shock. For $B=1~\text{mG}$ the gyration of electrons with energy ${\mathcal{E}}_{\text{max}}$ is of the order of the size of the active region. This implies that they cannot be trapped by the magnetic islands (ropes) as their size is much smaller compared to the length of the current sheet.

The total magnetic energy of the flare region is

(4.13) $$\begin{eqnarray}E_{m}=7\times 10^{38}B_{-3}^{2}l_{\text{ld}}^{3}~\text{erg},\end{eqnarray}$$

where $l_{\text{ld}}$ is the size of the region in light days. If a fraction $\unicode[STIX]{x1D716}_{m}$ of this energy is converted into emission with ${\mathcal{E}}_{\unicode[STIX]{x1D708}}>100~\text{MeV}$ then the corresponding mean luminosity of the flare is

(4.14) $$\begin{eqnarray}L=8\times 10^{33}\unicode[STIX]{x1D716}_{m}B_{-3}^{2}\,l_{\text{ld}}^{3}\,t_{f,d}^{-1}~\text{erg}~\text{s}^{-1},\end{eqnarray}$$

where $t_{f,d}$ is the flare duration in days. Since $l$ cannot exceed $ct_{f}$ , the highest possible luminosity is

(4.15) $$\begin{eqnarray}L_{\text{max}}=8\times 10^{33}\unicode[STIX]{x1D716}_{m}B_{-3}^{2}\,t_{f,d}^{2}~\text{erg}~\text{s}^{-1}.\end{eqnarray}$$

For the April 2011 flare $t_{f,d}=4$ and hence $L_{\text{max}}=1.5\times 10^{35}\unicode[STIX]{x1D716}_{m}B_{-3}^{2}~\text{erg}~\text{s}^{-1}$ whereas the observed isotropic luminosity above 100 MeV is $L_{\text{ob}}\approx 2\times 10^{36}~\text{erg}~\text{s}^{-1}$ . This suggests that either the magnetic field is a few times over 1 mG or the flare emission is beamed with the beaming factor $f_{b}=\unicode[STIX]{x0394}\unicode[STIX]{x1D6FA}/4\unicode[STIX]{x03C0}\approx L_{\text{max}}/L_{\text{ob}}\lesssim 0.08\unicode[STIX]{x1D716}_{m}B_{-3}^{2}$ , where $\unicode[STIX]{x0394}\unicode[STIX]{x1D6FA}$ is the beam solid angle. The angular distribution of particles accelerated in the current layers is not isotropic and the degree of anisotropy increases with the particle energy (Cerutti et al. Reference Cerutti, Werner, Uzdensky and Begelman2012b , Paper I). This kinetic beaming of accelerated particles is reflected in the anisotropy of their synchrotron emission with the beaming factor near the peak of the emission spectrum $f_{b}\approx 0.1$ , which is consistent with our estimate for the April 2011 flare.

If the reconnecting region is moving relative to the Earth with a relativistic speed then its observed emission can also be beamed due to the Doppler effect. The half-opening angle of the Doppler beam corresponding to the solid angle $\unicode[STIX]{x0394}\unicode[STIX]{x1D6FA}/4\unicode[STIX]{x03C0}\approx 0.08$ is $\unicode[STIX]{x1D703}_{d}\approx 35^{\circ }$ , which corresponds to the bulk speed $v\approx 0.5c$ . Such speeds are typical for the inner Crab Nebula. Even higher speeds are expected in the polar region at the base of the Crab’s jet on theoretical grounds (see § 2.2).

4.3 Spectrum of emitting particles

One can estimate the number of gamma-ray-emitting leptons in the flaring region as $N_{\unicode[STIX]{x1D6FE}}\approx Lt_{\text{cool}}/{\mathcal{E}}_{\text{max}}$ . Assuming that the cooling time $t_{\text{cool}}=t_{\text{acc}}$ and using (4.8) and (4.15), we find

(4.16) $$\begin{eqnarray}N_{\unicode[STIX]{x1D6FE}}\approx 5.6\times 10^{35}\frac{\unicode[STIX]{x1D716}_{m}}{\unicode[STIX]{x1D702}}B_{-3}t_{f,d}^{2}.\end{eqnarray}$$

(Comparing this with the pair-production rate of the pulsar, equation (2.11), one can see that it takes a fraction of a second to generate the particles producing the flares.) The corresponding volume density

(4.17) $$\begin{eqnarray}n_{\unicode[STIX]{x1D6FE}}\approx 3\times 10^{-11}B_{-3}t_{f,d}^{-1}~\text{cm}^{-3}\end{eqnarray}$$

is significantly smaller than the number densities of both the high-energy and low-energy lepton populations in the nebula. This shows that only a tiny fraction of particles in the flaring region are accelerated to the energies required to produce its gamma-ray emission.

Most studies of particle acceleration in magnetic reconnection come up with a single power law for the non-thermal particles. Our study also indicates a clear possibility of two populations of accelerated particles: a less energetic population with a softer spectrum which is made out of particles trapped inside the magnetic islands (ropes) and a more energetic population with a very hard spectrum made out of particles which are free inside the current sheet (paper I). Only the latter one would produce the flare emission. When the synchrotron radiation losses are taken into account one can also expect a pile up of particles near the radiation-reaction limit (Clausen-Brown & Lyutikov Reference Clausen-Brown and Lyutikov2012). Keeping these caveats in mind we will use here a simple model which involves a single power-law spectrum truncated at both ends of the energy domain

(4.18) $$\begin{eqnarray}\frac{\text{d}n}{\text{d}\unicode[STIX]{x1D6FE}}=A\unicode[STIX]{x1D6FE}^{-p},\quad \unicode[STIX]{x1D6FE}_{\text{min}}\leqslant \unicode[STIX]{x1D6FE}\leqslant \unicode[STIX]{x1D6FE}_{\text{max}}.\end{eqnarray}$$

The corresponding spectrum of the synchrotron photons is also a power law $\text{d}n_{\unicode[STIX]{x1D708}}/\text{d}{\mathcal{E}}\propto {\mathcal{E}}^{-a}$ , where

(4.19) $$\begin{eqnarray}p=1+2(a-1).\end{eqnarray}$$

For the April 2011 flare this gives $p=1.54\pm 0.24$ . Since $1<p<2$ one can estimate the total number density as

(4.20) $$\begin{eqnarray}n_{\text{tot}}=A\unicode[STIX]{x1D6FE}_{\text{min}}^{1-p}/(p-1)\end{eqnarray}$$

and the density of gamma-ray-emitting leptons as

(4.21) $$\begin{eqnarray}n_{\unicode[STIX]{x1D6FE}}=A\unicode[STIX]{x1D6FE}_{\text{max}}^{1-p}.\end{eqnarray}$$

Combining the two, we find

(4.22) $$\begin{eqnarray}\frac{n_{\text{tot}}}{n_{\unicode[STIX]{x1D6FE}}}=\frac{1}{p-1}\left(\frac{\unicode[STIX]{x1D6FE}_{\text{max}}}{\unicode[STIX]{x1D6FE}_{\text{min}}}\right)^{p-1}.\end{eqnarray}$$

The corresponding mean Lorentz factor

(4.23) $$\begin{eqnarray}\langle \unicode[STIX]{x1D6FE}\rangle =\unicode[STIX]{x1D6FE}_{\text{max}}\frac{p-1}{2-p}\left(\frac{\unicode[STIX]{x1D6FE}_{\text{max}}}{\unicode[STIX]{x1D6FE}_{\text{min}}}\right)^{1-p},\end{eqnarray}$$

has to be equal to $\unicode[STIX]{x1D70E}_{M}$ . From the last two equations we find the simple equation

(4.24) $$\begin{eqnarray}\unicode[STIX]{x1D6FE}_{\text{max}}=(2-p)\unicode[STIX]{x1D70E}_{M}\left(\frac{n_{\text{tot}}}{n_{\unicode[STIX]{x1D6FE}}}\right),\end{eqnarray}$$

which allows us to check the self-consistency of the model. Substituting the expressions for $\unicode[STIX]{x1D70E}_{M}$ from (2.15), $n_{\unicode[STIX]{x1D6FE}}$ from (4.17), using $t_{d}=4$ and assuming that the wind supplies only the high-energy leptons ( $n_{\text{tot}}\approx n_{h}$ ), we obtain

(4.25) $$\begin{eqnarray}\unicode[STIX]{x1D6FE}_{\text{max}}\approx 3.4\times 10^{9}\unicode[STIX]{x1D706}_{4}^{-1}B_{-3}^{-1},\end{eqnarray}$$

which is very close to the estimate (4.4). If the wind supplies the radio electrons as well then $n_{\text{tot}}\approx n_{l}$ and

(4.26) $$\begin{eqnarray}\unicode[STIX]{x1D6FE}_{\text{max}}\approx 3\times 10^{9}\unicode[STIX]{x1D706}_{7}^{-1}B_{-3}^{-1},\end{eqnarray}$$

which is an equally good result. Next we can use (4.22) to find

(4.27) $$\begin{eqnarray}\unicode[STIX]{x1D6FE}_{\text{min}}=\unicode[STIX]{x1D6FE}_{\text{max}}\left(\frac{n_{\unicode[STIX]{x1D6FE}}}{n_{\text{tot}}}\frac{1}{p-1}\right)^{1/(p-1)}.\end{eqnarray}$$

In this equation $\unicode[STIX]{x1D6FE}_{\text{min}}$ is very sensitive to the power index. In the model where the wind supplies only the high-energy population, we find that for $p=1.54\pm 0.24$ $\unicode[STIX]{x1D6FE}_{\text{min}}$ varies in the range $\unicode[STIX]{x1D6FE}_{\text{min}}\in (3\times 10^{3},5\times 10^{6})$ with the median value of $5\times 10^{5}$ , corresponding to the optical electrons. If the wind supplies the radio electrons as well then $\unicode[STIX]{x1D6FE}_{\text{min}}\in (4\times 10^{-7},750)$ with a median value of $1.6$ , showing that only $p\leqslant 1.54$ are acceptable. Interestingly, the observed spectral index of the radio electrons $p_{r}=1.54\pm 0.08$ (Kovalenko, Pynzar’ & Udal’Tsov Reference Kovalenko, Pynzar’ and Udal’Tsov1994; Bietenholz et al. Reference Bietenholz, Kassim, Frail, Perley, Erickson and Hajian1997) is very close to that of the gamma-ray flare electrons.

4.4 Plasma parameters

The observed spectra of flares provide some information on the plasma magnetisation of the flaring region. According to the PIC simulations (Sironi & Spitkovsky Reference Sironi and Spitkovsky2014; Guo et al. Reference Guo, Liu, Daughton and Li2015, Papers I and II,), $p\approx 1.5$ requires the plasma magnetisation $\unicode[STIX]{x1D70E}$ to be of the order of a few tens, with $\unicode[STIX]{x1D70E}>100$ leading to very hard spectra with $p\rightarrow 1$ and $\unicode[STIX]{x1D70E}<10$ to soft spectra $p>2$ .

Another important parameter is the electron skin depth $\unicode[STIX]{x1D6FF}_{s}=c/\unicode[STIX]{x1D714}_{p}$ , where $\unicode[STIX]{x1D714}_{p}^{2}=4\unicode[STIX]{x03C0}e^{2}n/(m_{e}\unicode[STIX]{x1D6FE}_{t})$ is the plasma frequency and $\unicode[STIX]{x1D6FE}_{t}$ is the mean Lorentz factor of particles in the reconnection inflow. Assuming that the flare occurs in a plasma whose bulk motion Lorentz factor is of the order of unity but $\unicode[STIX]{x1D70E}\gg 1$ , one can estimate $\unicode[STIX]{x1D6FE}_{t}\approx \unicode[STIX]{x1D70E}_{M}/\unicode[STIX]{x1D70E}$ and hence

(4.28) $$\begin{eqnarray}\unicode[STIX]{x1D714}_{p}^{2}\approx \frac{4\unicode[STIX]{x03C0}e^{2}}{m_{e}}\frac{n}{\unicode[STIX]{x1D70E}_{M}}\unicode[STIX]{x1D70E},\end{eqnarray}$$

which is sensitive to what we assume about the origin of radio electrons. If they are not supplied by the wind then $n\approx n_{h}$ and

(4.29) $$\begin{eqnarray}\unicode[STIX]{x1D6FF}_{s}\approx \frac{7\times 10^{-3}}{\sqrt{\unicode[STIX]{x1D706}_{4}\unicode[STIX]{x1D70E}}}~\text{ld}\end{eqnarray}$$

and if they are then $n\approx n_{l}$ and

(4.30) $$\begin{eqnarray}\unicode[STIX]{x1D6FF}_{s}\approx \frac{8\times 10^{-6}}{\sqrt{\unicode[STIX]{x1D706}_{7}\unicode[STIX]{x1D70E}}}~\text{ld}.\end{eqnarray}$$

Thus, the size of the flaring region $ct_{f}\approx \text{few}\times 10^{3}\unicode[STIX]{x1D6FF}_{s}$ and $\text{few}\times 10^{6}\unicode[STIX]{x1D6FF}_{s}$ respectively. The corresponding plasma gyration radius

(4.31) $$\begin{eqnarray}r_{g}=\frac{m_{e}c^{2}\unicode[STIX]{x1D6FE}_{t}}{eB}=\unicode[STIX]{x1D6FF}_{s}/\unicode[STIX]{x1D70E}\end{eqnarray}$$

is even smaller than $\unicode[STIX]{x1D6FF}_{s}$ . Thus, the Crab flares undoubtedly involve macroscopic scales.

As the wind flow is brought to a virtual halt in the flaring site, $\unicode[STIX]{x1D6FE}_{w}\approx 1$ and

(4.32) $$\begin{eqnarray}\unicode[STIX]{x1D70E}_{M}=h\unicode[STIX]{x1D70E}=\frac{B^{2}}{4\unicode[STIX]{x03C0}m_{e}c^{2}n}.\end{eqnarray}$$

This yields an independent estimate of the reconnecting magnetic field strength, namely

(4.33) $$\begin{eqnarray}B=(4\unicode[STIX]{x03C0}m_{e}c^{2}n\unicode[STIX]{x1D70E}_{M})^{1/2}\approx 0.7~\text{mG}.\end{eqnarray}$$

The result is not sensitive to the assumption on the origin of the radio electrons (If the radio electrons are supplied by the wind then we should use $n=n_{l}$ and $\unicode[STIX]{x1D706}=10^{7}$ and if they are not then $n=n_{h}$ and $\unicode[STIX]{x1D706}=10^{7}$ , which leaves $n\unicode[STIX]{x1D70E}_{M}$ unchanged.). This is quite close to the lower limit on the magnetic field strength obtained in § 4.2.

4.5 Potential–luminosity relationship and the Alfvén current

There is an important, yet somewhat subtle, relationship between the parameters of the flares that serve as a confirmation of the basic principles of our approach, as we discuss next. In addition to particle acceleration, the reconnecting electric field also drives (a charge-separated) current. This induced current creates its own magnetic field; particle gyration in the induced magnetic field can inhibit the beam propagation. These considerations place an upper limit on the total charger-separated current, known as the Alfvén current. The relativistic Alfvén current (Alfvén Reference Alfvén1939; Lawson Reference Lawson1973) is

(4.34) $$\begin{eqnarray}I_{A}=\unicode[STIX]{x1D6FE}\frac{m_{e}c^{3}}{e}.\end{eqnarray}$$

In the present model of Crab flares, the reconnection events result in a population of gamma-ray-emitting particles. Those particles create electric current, which in turn creates its own magnetic field. For the consistency of the model it is required that the Larmor radius of the particles in their own magnetic field is sufficiently large. Next, using the total flare luminosity and the estimates of the energy of emitting particles we can verify that the current carried by the flare-emitting particles does not violate the Alfvén limit.

Qualitatively, the electromagnetic power of relativistic outflows can be estimated as (Blandford & Znajek Reference Blandford and Znajek1977; Blandford Reference Blandford, Gilfanov, Sunyeav and Churazov2002)

(4.35) $$\begin{eqnarray}L_{EM}\approx V^{2}/{\mathcal{R}},\end{eqnarray}$$

where $V$ is a typical values of the electric potential produced by the central engine and ${\mathcal{R}}\approx 1/c$ is the impedance of free space. In case of Crab flares, the available potential,

(4.36) $$\begin{eqnarray}V\sim \sqrt{L/c}\approx 1.7\times 10^{15}~\text{V},\end{eqnarray}$$

corresponds to $\unicode[STIX]{x1D6FE}\approx 3.4\times 10^{9}$ – very close to the estimate (4.4).

Similarly, the electric current associated with flares

(4.37) $$\begin{eqnarray}I\approx \sqrt{Lc}\approx 5\times 10^{22}~\text{cgs units}=1.7\times 10^{13}~\text{A}\end{eqnarray}$$

is similar to the Alfvén current (4.34) for the required Lorentz factor (4.4):

(4.38) $$\begin{eqnarray}I_{A}=1.7\times 10^{23},~\text{cgs units}=6\times 10^{13}~\text{A}.\end{eqnarray}$$

This is important: the Lorentz factor that we derived from radiation modelling (from the peak energy and duration of flares) nearly coincides with the expected Lorentz factor of particles derived from the assumption that the radiated power is similar to the intrinsic electromagnetic power of the flares. In a related ‘coincidence’, the current required to produce a given luminosity nearly coincides with the corresponding Alfvén current, calculated using the Lorentz factor of the emitting particles.

These coincidences, in our view, are not accidental: they imply that flares are produced by charge-separated flow. This is indeed expected in the acceleration models based on DC-type acceleration. This is also consistent with our PIC simulation that shows the effects of kinetic beaming – a highly anisotropic distribution of the highest-energy particles.

5.1 Magnetisation in different parts of the nebula

After reviewing the theory and observations of the Crab Nebula and the analysis of its gamma-ray flares as magnetic reconnection events, we are now ready to identify the locations in the nebula where this flaring emission may be coming from. In fact, it is quite clear that for most of the Crab Nebula the physical conditions are not favourable for the flares. Firstly, the mean magnetic field of the nebula is approximately an order of magnitude below the minimum strength of one mG required for the accelerated particles to reach the radiation-reaction limit on the typical time scale of the flares (see § 4.2). Secondly, the mean magnetic energy density is approximately an order of magnitude below that of the synchrotron electrons, which implies that the mean magnetisation parameter can be at most $\unicode[STIX]{x1D70E}\approx 0.1$ . For such a low magnetisation, the reconnection rate can be at most $\unicode[STIX]{x1D702}\approx 0.1$ . This not only increases the required acceleration time even further, see (4.11), but also reduces the photon energy in radiation-reaction limit by an order of magnitude (see (4.7)). Finally, for $\unicode[STIX]{x1D70E}<1$ the spectrum of accelerated particles becomes too soft, with the power index $p>4$ (Sironi & Spitkovsky Reference Sironi and Spitkovsky2014; Guo et al. Reference Guo, Liu, Daughton and Li2015), whereas the observations suggest a hard power-law spectrum with $p\approx 1.5$ (Buehler et al. Reference Buehler, Scargle, Blandford, Baldini, Baring, Belfiore, Charles, Chiang, D’Ammando and Dermer2012) or even a monoenergetic spectrum (Striani et al. Reference Striani, Tavani, Vittorini, Donnarumma, Giuliani, Pucella, Argan, Bulgarelli, Colafrancesco and Cardillo2013). If the flares do occur inside the nebula then these must be rather special cites with much stronger magnetic field and much higher magnetisation. In principle, strong shocks may amplify the magnetic field. The RMHD computer simulations of PWN show that secondary shocks are indeed emitted by the highly variable termination shock of the pulsar wind but they are rather weak (Camus et al. Reference Camus, Komissarov, Bucciantini and Hughes2009; Porth et al. Reference Porth, Komissarov and Keppens2013). Moreover, RMHD shocks cannot change the plasma magnetisation from $\unicode[STIX]{x1D70E}\ll 1$ to $\unicode[STIX]{x1D70E}\gg 1$ . The simulations also indicate that near the termination shock the physical conditions can be more extreme so we now turn our attention to the very vicinity of the shock.

In the analysis of the termination shock and physical parameters of the plasma injected by the pulsar wind into the nebula one should differentiate between its polar section, terminating the unstriped polar part of the wind, and its equatorial section, terminating the striped equatorial part of the wind. As Lyutikov et al. (Reference Lyutikov, Komissarov and Porth2016) discuss, the modelling of the inner knot of the Crab Nebula requires that the equatorial flow is weakly magnetised. Thus, the physical conditions downstream of the striped section are still not quite suitable for a flaring site. Indeed, whether the radio electrons are supplied by the wind or not, the downstream magnetisation is very low, with $\unicode[STIX]{x1D70E}\ll 1$ . If the radio electrons are supplied by the wind, then its stripes dissipate and its magnetisation drops to very low levels before the termination shock. If the radio electrons are not supplied by the wind then the stripes dissipate at the termination shock. As we have discussed, for the parameters of the Crab wind this dissipation leads only to plasma heating (§ 2.2) and the downstream magnetisation is too low to make gamma-ray flares via reconnection of the residual magnetic field. The observations of the inner knot of the Crab Nebula, which is best explained as a Doppler-boosted emission of the shocked striped wind (Lyutikov et al. Reference Lyutikov, Komissarov and Porth2016), agree with this verdict – its optical and IR emission does not correlate with the gamma-ray flares (Rudy et al. Reference Rudy, Horns, DeLuca, Kolodziejczak, Tennant, Yuan, Buehler, Arons, Blandford and Caraveo2015).

What remains is the poleward zone of the nebula which is filled with the plasma supplied by the unstriped polar section of the pulsar wind. As we discussed in § 2.2, here the magnetisation can be very high $\unicode[STIX]{x1D70E}\gg 1$ leading to the reconnection rate $\unicode[STIX]{x1D702}\approx 1$ and very hard spectra of accelerated particles. Moreover, the polar region is over-pressured compared to the rest of the nebula due the hoop stress of the azimuthal magnetic field freshly supplied by the pulsar wind and hence its total pressure and the magnetic field strength can be much higher. Early axisymmetric numerical solutions, based on models of weakly magnetised pulsar winds, exhibited such magnetic pinch all the way to the outer edge of the nebula (e.g. Camus et al. Reference Camus, Komissarov, Bucciantini and Hughes2009). However in the recent 3-D RMHD simulations of PWN (Porth et al. Reference Porth, Komissarov and Keppens2013; Olmi et al. Reference Olmi, Del Zanna, Amato, Bucciantini and Mignone2016) such a strong compression is found only at the base of the polar jet, on a scale not exceeding the equatorial size of the termination shock. Further out the z-pinch appears to be destroyed by instabilities.

The total potential drop across the polar zone of the Crab’s pulsar wind is also sufficiently high to accommodate particle acceleration up to $\unicode[STIX]{x1D6FE}\approx 3\times 10^{9}$ , required by the observations, unless the opening angle of the polar zone is below $15^{\circ }$ for $n=1$ and $27^{\circ }$ for $n=2$ .

As the polar angle decreases below the value of $\unicode[STIX]{x1D703}_{p}$ separating the striped and polar wind zones, the termination shock shrinks dramatically in reaction to the much higher magnetisation of the polar zone (Lyubarsky Reference Lyubarsky2012; Lyutikov et al. Reference Lyutikov, Komissarov and Porth2016). Using the asymptotic solution given in Lyutikov et al. (Reference Lyutikov, Komissarov and Porth2016), we find that in the polar zone the shock’s spherical radius

(5.1) $$\begin{eqnarray}R(\unicode[STIX]{x1D703})\simeq \frac{R_{o}}{\sqrt{2\unicode[STIX]{x1D70E}_{u}}}\left\{\begin{array}{@{}ll@{}}\unicode[STIX]{x1D703}^{2}/2 & \text{for }n=1\\ \unicode[STIX]{x1D703}^{3}/3 & \text{for }n=2,\end{array}\right.\end{eqnarray}$$

where $R_{o}=(S_{p}/4\unicode[STIX]{x03C0}cp_{p})^{1/2}$ and $p_{p}$ is the total pressure of the polar region. The equatorial radius of the termination shock can be estimated as $R_{\text{ts}}=(S_{p}/4\unicode[STIX]{x03C0}cp_{n})^{1/2}$ , where $p_{n}$ is the mean total pressure of the nebula. Given $p_{p}\approx 10p_{n}$ and $\unicode[STIX]{x1D70E}_{u}\approx 4\times 10^{4}M^{-2/3}\unicode[STIX]{x1D706}_{4}^{-2/3}$ predicted by the ideal MHD model of the pulsar wing, the maximum distance of the termination shock from the pulsar can be well below one light day. Even if fast magnetosonic waves are responsible for a non-negligible fraction of the initial energy transport in the polar zone and dissipate before reaching the termination shock, the distance is still rather short. A two per cent fraction leads to $\unicode[STIX]{x1D70E}_{u}=50$ and $R<3~\text{ld}$ . So we may be dealing with a rather compact region here.

For $\unicode[STIX]{x1D70E}_{u}\gg 1$ the shock is weak, the speed of the downstream flow is very close to the speed of light and it is highly magnetised. It has to somehow decelerate further in order to match the slow expansion rate of the nebula. Lyubarsky (Reference Lyubarsky2012) analysed a steady-state axisymmetric model of this polar flow assuming that it is confined by the mean pressure of the nebula, using the approximation of $\unicode[STIX]{x1D703}_{p}\ll 1$ and assuming $n=1$ . He concluded that the flow converges on the scale

(5.2) $$\begin{eqnarray}R_{c}=\frac{\unicode[STIX]{x03C0}}{\sqrt{6}}\unicode[STIX]{x1D703}_{p}^{2}R_{\text{ts}}\end{eqnarray}$$

and argued that at this point it will be destroyed by instabilities. For $\unicode[STIX]{x1D703}_{p}=30^{\circ }$ , equation (5.2) yields $R_{c}\approx 0.35R_{\text{ts}}$ . Given that the external pressure in the polar region can exceed the mean pressure of the nebula by an order of magnitude, $R_{c}$ can be reduced by a factor of three, leading to $R_{c}\approx 0.1R_{\text{ts}}\approx 17~\text{ld}$ in the case of the Crab.

5.2 Current filamentation

Recent 3-D RMHD simulations of PWN produced by winds with distinctly different polar and equatorial zones (Porth, Komissarov & Keppens Reference Porth, Komissarov and Keppens2014) offer further insight into the dynamics of the polar region. As the mildly relativistic post-shock flow decelerates, it comes into a global causal contact. Once causally connected, the almost purely toroidal magnetic configuration of the nebula flow is liable to current-driven instabilities Begelman (Reference Begelman1998), Mizuno et al. (Reference Mizuno, Lyubarsky, Nishikawa and Hardee2011). As pointed out by Lyubarsky (Reference Lyubarsky2012), in a highly magnetised expanding post-shock flow, the flow lines are in fact focussed to the axis via the magnetic hoop stress which can act as an additional trigger for instability. Thereby the innermost flow reaches the recollimation point first and is expected to experience violent instability and magnetic dissipation. Lyubarsky (Reference Lyubarsky2012) hence argues that the gamma-ray flares originate in the polar region at the base of the observed jet/plume. On the one hand, the simulations of Porth et al. (Reference Porth, Komissarov and Keppens2014) are in agreement with this dynamical picture and show development of modes with azimuthal number $m\sim \text{few}$ , that lead to the formation of current filaments at the jet base. On the other hand, for numerical stability, Porth et al. (Reference Porth, Komissarov and Keppens2014) could only study moderate polar magnetisations of $\unicode[STIX]{x1D70E}\leqslant 3$ which disfavour rapid particle acceleration. Even without this technical limitation, for non-vanishing wind power, the axis itself cannot carry a DC Poynting flux. Because of this, we are interested in flow lines from intermediate latitudes that acquire high $\unicode[STIX]{x1D70E}$ due to flow expansion.

To visualise the current in the simulations, we now cut the domain by spherical surfaces and investigate current magnitude and magnetisation in the $\unicode[STIX]{x1D719}\unicode[STIX]{x1D703}$ -plane. In general, the current exhibits a quadrupolar morphology, it is distributed smoothly over the polar regions and returns near the equator. When the test surface crosses over the termination shock, we observe an increased return current due to the shock jump of $B_{\unicode[STIX]{x1D719}}$ . Although the field is predominantly toroidal in the high $\unicode[STIX]{x1D70E}$ downstream region, perturbations from the unstable jet occasionally lead to current filaments and even reversals just downstream of the termination shock.

With increasing distance from the termination shock, the flow becomes more and more turbulent which leads to increased filamentation and copious sign changes of the radial current. In the underlying simulations, the magnetisation in the turbulent regions has dropped already below unity which disfavours rapid particle acceleration at distances larger than a few termination shock radii. It is interesting to point out that although the magnetic field injected in the model decreases gradually over a large striped region of $45^{\circ }$ , an equatorial current sheet develops nonetheless in the nebula. This can be seen in figure 2(c,d) which corresponds to a surface outside of the termination shock at $r=3\times 10^{17}~\text{cm}$ .

Figure 2.

Distribution of the radial current, $j_{r}\sin \unicode[STIX]{x1D703}$ (a,c) and magnetisation $\unicode[STIX]{x1D70E}$ (b,d). Dashed contours indicate the sign changes of the current and solid contours indicate levels of $\unicode[STIX]{x1D70E}=4$ . (a,b) Sphere with $r=10^{17}~\text{cm}$ , one can see formation of current filaments (e.g., regions $\unicode[STIX]{x1D703}\approx 2.5$ and $\unicode[STIX]{x1D719}\approx -1.5,-0.5$ ). The current direction can even reverse (the region $\unicode[STIX]{x1D703}\approx 2.7$ , $\unicode[STIX]{x1D719}\approx -1$ ). (c,d) Sphere with $r=3\times 10^{17}~\text{cm}$ . Here the current is highly filamentary and an equatorial current sheet has developed. The average magnetisation in the current sheet and in the turbulent flow is ${<}1$ but magnetically dominated regions can still be obtained ( $\unicode[STIX]{x1D70E}\approx 6$ at $\unicode[STIX]{x1D703}\approx 1.75$ , $\unicode[STIX]{x1D719}\approx -1.5$ ).

Figure 3 shows the results for a wind with $\unicode[STIX]{x1D703}_{p}=45^{\circ }$ and $\unicode[STIX]{x1D70E}_{u}\approx 1$ in the polar zone. As one can see, the polar flow does indeed converge towards the polar axis with $R_{c}\approx 0.5R_{\text{ts}}$ . Moreover, a significant fraction of the shocked equatorial part of the wind also converges, in agreement with the earlier 2-D simulations of weakly magnetised pulsar winds. Furthermore, the polar flow does not look at all like a quasi-stationary axisymmetric jet. It is highly disturbed already on the scales $R\ll R_{c}$ , with even a negative radial velocity in one region.

Figure 3.

Magnetisation $\unicode[STIX]{x1D70E}$ and instantaneous streamlines near the termination shock in RMHD simulations of the Crab Nebula (Porth et al. Reference Porth, Komissarov and Keppens2014). The dot-dashed straight line shows the separation of the polar and striped zones of the pulsar wind. The dashed straight line is the line of sight. The solid red line shows the termination shock and the solid blue loop between the dashed and dot-dashed lines shows the region of Doppler-beamed emission associated with the inner knot of the nebula. The polar beam corresponds to the streamlines originating from the inner part of the termination shock located to the left of the intersection with the dot-dashed line.

To better understand the geometry of the current and magnetic field, we display a representative volume rendering of the polar region in figure 4, which shows the complicated structure of field lines in the polar zone on the scales of $r\in (0.1,0.3)R_{\text{ts}}$ . In the rendering in figure 4 one can see the violently unstable polar beam embedded in the more regular high- $\unicode[STIX]{x1D70E}$ region comprised of toroidal field lines. The plume forms downstream of this structure and is also strongly perturbed. A part of the disrupted plume approaches the termination shock as a flux tube. The presence of such a configuration where two flux tubes can come together very close to the high $\unicode[STIX]{x1D70E}$ region lets us speculate that Crab flares might originate when the right geometry (e.g. parallel flux tubes) coincides with high magnetisation as present in the nebula even for moderate wind magnetisation. For higher magnetisations of the polar beam, the mechanism described by Lyubarsky (Reference Lyubarsky2012) could directly act also without first having to rely on enhancement of $\unicode[STIX]{x1D70E}$ via flow expansion. Extrapolating from the moderate $\unicode[STIX]{x1D70E}$ simulations where the polar beam is highly unstable and forms a filamented current, this seems feasible at the very least.

The overall impression is that the converging part of equatorial flow creates a highly dynamic ‘wall’ that almost stops and crushes the polar flow, which loses its regular structure and develops numerous current sheets well before the convergence point. This must be followed by a rapid forced magnetic reconnection akin the one observed in the simulations of shocks in plasma with magnetic stripes (Sironi & Spitkovsky Reference Sironi and Spitkovsky2011) and in our simulations of magnetic reconnection driven by macroscopic stresses (papers I and II). The conversion of magnetic energy into the energy of particles during the magnetic reconnection increases the plasma compressibility and this may be the mechanism ultimately responsible for the observed slowing down of the polar beam.

Figure 4.

Three-dimensional volume rendering showing current filamentation of the polar beam just downstream of the termination shock. The shock surface is indicated as the orange plane and we draw field lines shaded from white $(\unicode[STIX]{x1D70E}=0)$ to black ( $\unicode[STIX]{x1D70E}=5$ ). One clearly sees two current filaments producing structures similar to magnetic flux tubes. As discussed in Porth et al. (Reference Porth, Komissarov and Keppens2014), streamlines from intermediate latitudes reach the axis behind this inner violently unstable region and form a plume-like outflow of moderate velocity $v\approx 0.7c$ .

Summarising, the poleward regions downstream of the shock terminating the unstriped section of the pulsar wind seems to be the only promising location for the gamma-ray flares of the Crab Nebula. This is a rather compact central region and in the maps of the Crab Nebula it occupies the area inside the X-ray ring of the Crab Nebula. The flare region should be located at (or extend up to) intermediate polar angles in the range of $10^{\circ }{-}30^{\circ }$ . The required acceleration potential within the region is also a substantial fraction (a few per cent) of the total potential of the wind. The flare region has a macroscopic size, orders of magnitude larger than the corresponding skin depth.