The study of the emission patterns of gamma-ray pulsars has progressed ever since their discovery. There have been different models proposed in the past to explain the high energy radiation, i.e. the polar cap model (Ruderman & Sutherland 1975), the outer gap model (Cheng et al. 1986a, 1986b) and the slot gap model (Arons 1983). And, ever since the launch of the
A previous model, the simple two-layer outer gap model (Wang et al. 2010), was used to study the phase-averaged spectrum of pulsars. In this model, the outer gap is considered as a combination of a main acceleration region and a screening region. The charge density distribution of the outer gap was assumed to be a simple step function. The acceleration region has a charge density much lower than the Goldreich-Julian density and is nearly a vacuum, which raises a strong electric field; particles are accelerated and curvature photons with energy in the GeV range are emitted. The particles produced by the pair creation process stop the growth of the acceleration region and form the screening region, which produce photons with energy of around 100 MeV. However, this can only explain the phase-averaged spectra, not the phase-resolved spectra. A previous study has considered the azimuthal structure of the outer gap (Wang et al. 2011) and the phase-resolved spectra and light curves of pulsars. In their work, they tested their model using the Vela pulsar, which has an energy dependent third peak in the light curve. The model successfully explained the phase-averaged spectra, the phase-resolved spectra, and the light curve.
One study has shown that non-thermal X-ray emission from the outer magnetosphere could be generated by Inverse-Compton Scattering between primary particles in the gap and the radio waves coming towards the pulsar (Ng et al. 2014). When the energy of the primary particles is large enough, the Inverse-Compton Scattering is capable of producing GeV photons.
In this paper, we assume, based on the three-dimensional two-layer outer gap model, that the magnetosphere of a pulsar is time-dependent. Instead of solving a time-dependent Poisson equation, we assume that one of the spectral parameters has a power-law distribution. For a certain value of the parameter we assign a possibility. Hence when we consider the distribution of the parameter, we can obtain a weighted spectrum. And, using the Inverse-Compton Scattering component, we can calculate the high energy tail of powerful pulsars.
This paper is organized as follows. In section 2, we present a simple non-stationary two-layer three-dimensional outer gap model. Here, instead of solving a time-dependent Poisson equation, we choose to give a probability distribution to the outer gap structure, showing that the structure of the outer gap changes over time. Then, to explain the slow dropping high energy tail of Geminga pulsar, we discuss the Inverse Compton Scattering component in relation to the high energy part. In section 3, we present the model fitting results of the non-stationary model and the model combining the Inverse-Compton Scattering component. In sections 4 and 5, we present the discussion and conclusions.
2.1 Non-stationary Outer Gap Model
It is known that the outer gap extends from the last open field lines to the region between the null charge surface and the light cylinder. The previous study by Wang et al. (2010, 2011) assumed that the outer gap can be divided into two parts, the acceleration and the screening regions. They assumed that in the trans-field direction the charge density would have a simple step function, as shown below:
In the acceleration region, the charge density is assumed to be around 10% of the Goldreich-Julian density, thus producing a large electric field that accelerates the particles to ultra-relativistic velocity and produces GeV curvature photons. In the screening region, which is around the upper boundary, the charge density is larger than the Goldreich-Julian density due to the pair creation process.
When we assume that the physical quantities change slowly in the azimuthal angle
Here we have five spectral parameters,
However, when we tried to apply this model to the Geminga pulsar, we found that when the energy of the photons exceeded the cutoff energy, the spectrum dropped faster than expected, which means that if this model is used, a flux of above GeV photons will be underestimated . We introduced the non-stationary model to solve this problem. We believe that the structure of the outer gap is not stationary, but changes from time to time. In this case, the Poisson equation needs to add a dimension for time, and by solving the time-dependent Poisson equation, this can be the model we are looking for. However, solving time-dependent Poisson equations is very complicated. In a simple case, we can assume a distribution for one of the spectral parameters proposed above. For one chosen parameter, we can assume a power-law distribution,
2.2 Inverse-Compton Scattering in Outer Gap
Although to some degree the non-stationary model can explain the spectral property of gamma ray pulsars, we noticed that for some very powerful pulsars like the Geminga pulsar there were difficulties. Ng et al. (2014) discovered that besides synchrotron radiation from secondary pairs in the outer gap, there is another way to produce X-ray emission for millisecond pulsars with large magnetic fields: Inverse-Compton Scattering between the radio emission and the primary particles in the outer gap. According to their study, as the radio, X-ray, and gamma ray emission regions are likely to be at the same positions, so the radio waves emitted from the outer gap can, then, scatter with the primary particles with ultra-relativistic energy. To compare this energy with the curvature radiation from the outer gap, we may argue as follows.
We can conclude from this equation that the flux predicted by Inverse-Compton Scattering has about the same magnitude as the curvature radiation, and thus is observable for the pulsars. To calculate the Inverse-Compton Scattering spectrum, Ng et al. (2014) considered the scattering between the outwards primary particles and the outwards radio waves. The radio emission region is assumed to be just above the outer gap; for each particle, the power per unit energy per unit solid angle of the Inverse-Compton scattering is
Here
As we can see from Eq. (8), the spectrum of Inverse-Compton Scattering relies a lot on the radio spectrum. Although this spectrum is currently unknown, it is suggested that at a low frequency of less than 1 GHz, there will be a turnover point in the spectrum at around 100 MHz for millisecond pulsars (Kuzmin & Losovsky 2001); we set this frequency as the turnover value. After comparing the theoretical results with the observations, we here assume a broken power-law radio spectrum:
We also set the spectral indexes
In this section, we apply our model to the Geminga pulsar. We chose to give a power-law distribution to the charge density, in Eq. (5). In the previous study we found that the Geminga pulsar is very powerful, which means that it needs a large outer gap to provide its emission, and thus the fitted value of the fractional gap size (
Parameters used in the non-stationary three-dimensional outer gap model to fit the phase-averaged spectrum, except for the parameter we assigned to the probability distribution, i.e. the charge density
In Fig. 3, the combined spectrum of Inverse-Compton Scattering and the outer gap is presented. When combining the Inverse-Compton scattering spectrum with the curvature radiation, we add only the part above 1 GeV.
In the beginning we believed that the high energy emission of the Geminga pulsar would be explained by the non-stationary model; however, from the results it is obvious that this is far from the reality. Different distribution patterns of selected parameters may lead to even better results; however, if we still want to undertake a detailed study of this problem, we cannot avoid solving a time-dependent Poisson equation.
It is argued that most photons emitted from the gap are from above the null charge surface, as shown in the results above; however, fellow researchers have predicted that the inner boundary of the gap can extend inwards (Takata et al. 2004), which may also be the reason for the lack of high energy (above 3 GeV) photons.
In this paper, we discussed the non-stationary three dimensional two-layer outer gap model and used it to explain the high energy emission (above approximately 3 GeV) of the Geminga pulsar, observed by