Modeling gamma-ray light curves with more realistic pulsar magnetospheres

TL;DR

This paper models gamma-ray light curves of pulsars using realistic magnetosphere simulations, comparing geometric and particle trajectory methods, and finds the particle trajectory approach aligns better with physical conditions and observations.

Contribution

It introduces a detailed particle trajectory method considering conductivity variations, providing more physically consistent gamma-ray light curve modeling for pulsars.

Findings

Low conductivity magnetospheres favor geometric modeling.

High conductivity outside the light cylinder improves particle trajectory modeling.

Particle trajectory method aligns better with observed gamma-ray light curves.

Abstract

We study the gamma-ray emission patterns and light curves in dissipative pulsar magnetospheres. We produce the gamma-ray light curves by using the geometric method and the particle trajectory method. For the geometric method, assuming the gamma-ray emission originates in a finite-width layer along the last closed lines, we generate the gamma-ray light curves based on the uniform emissivity along the magnetic field lines in the comoving frame (CF). For the particle trajectory method, we consider the spatial distribution of conductivity $\sigma$ by assuming a very high conductivity within the light cylinder and a finite conductivity outside the light cylinder (LC). Assuming that all the $\gamma$-ray emission originates in the outer magnetosphere outside the LC, we generate the gamma-ray light curves by computing realistic particle trajectories and Lorentz factors, taking into account both the accelerating electric field and curvature radiation loss. Further, we compare the modeling light curves to the observed light curves at $>0.1\, \rm GeV$ energies for Vela pulsar. Our results show that the magnetosphere with the low $\sigma$ value is more preferred for the geometric method. However, the magnetosphere with a near force-free regime within the LC and a high $\sigma$ value outside the LC is more favored for the particle trajectory method. It is noted that the particle trajectory method uses the parallel electric fields that are self-consistent with the magnetic fields of the magnetosphere. We suggest that the results from the particle trajectory method are more supported on the physical ground.

Equations

$$\frac{1}{c}\frac{\partial \mathbf{B}}{\partial t}$$

$$\frac{1}{c}\frac{\partial \mathbf{E}}{\partial t}$$

$$\nabla \cdot \mathbf{B}$$

$$\nabla \cdot \mathbf{E}$$

$$\mathbf{J}=c{\rho }_{\mathrm{e}}\frac{\mathbf{E}\times \mathbf{B}}{B^2}+\frac{c}{4\pi }\frac{(\mathbf{B}\cdot \nabla \times \mathbf{B}-\mathbf{E}\cdot \nabla \times \mathbf{E})\mathbf{B}}{B^2}.$$

$$\mathbf{J}=c{\rho }_{\mathrm{e}}\frac{\mathbf{E}\times \mathbf{B}}{B^2+E_0^2}+\sigma {\mathbf{E}}_{\Vert },$$

$$B_0^2-E_0^2$$

$$E_0{B}_0$$

$$\mathbf{v}=\left(\frac{\mathbf{E}\times \mathbf{B}}{B^2+E_0^2}+f\frac{\mathbf{B}}{B}\right)c,$$

$${\mu }_{\mathrm{em}}={\beta }_z,{\varphi }_{\mathrm{em}}=\mathrm{atan}\left(\frac{{\beta }_{\mathrm{y}}}{{\beta }_{\mathrm{x}}}\right),$$

$$\Phi ={\varphi }_{\mathrm{rot}}-{\varphi }_{\mathrm{em}}-{\mathbf{r}}_{\mathrm{em}}\cdot {\mathbf{\eta }}_{\mathrm{em}}/r_{\mathrm{L}},$$

$$\frac{d\gamma }{dt}=f\frac{q_{\mathrm{e}}c{E}_{\Vert }}{m_{\mathrm{e}}{c}^2}-\frac{2q_{\mathrm{e}}^2{\gamma }^4}{3R_{\mathrm{CR}}^2{m}_{\mathrm{e}}c},$$

Full text

1. Introduction

The launch of Fermi Gamma-Ray Space Telescope has opened a new era in the study of pulsar physics. In the first year of operation of the Large Area Telescope (LAT), more than 40 new gamma-ray pulsars were discovered (Abdo et al., 2010). To date, more than 150 gamma-ray pulsars have been detected by the LAT (Ackermann et al., 2015), 117 of which are included in the Second Fermi Pulsar Catalog (Abdo et al., 2013). They are divided into three groups: millisecond pulsars, young radio-loud pulsars and young radio-quiet pulsars. High-quality light curves, phase-averaged spectra and phase-resolved spectra obtained by the LAT observations provide an opportunity of understanding the emission sites and magnetic field geometries of these gamma-ray pulsars.

The realistic structures of the pulsar magnetosphere still remain uncertain. Knowledge about the pulsar magnetosphere structures can be used to identify the potential sites of particle acceleration and gamma-ray emission. A vacuum dipole field is generally adopted in the early study of pulsar emission, because it has an exact analytical solution given by Deutsch (1955, hereafter VRD). Based on this field structure, different theoretical models have been developed to explain the observed pulsar emission. In these models, it is widely believed that particles are accelerated in the gap region where an accelerating electric field is created because of the deficit of charges. Gamma-rays emission is produced by curvature or inverse-Compton radiation from high-energy particles accelerated in these gaps. Due to different emission zone locations, standard pulsar radiation models include the polar cap (e.g., Ruderman & Sutherland, 1975; Daugherty & Harding, 1982), the slot-gap (SG) (e.g., Dyks & Rudak, 2003; Dyks et al., 2004; Muslimov & Harding, 2004), and the outer-gap (OG) (e.g., Cheng et al., 1986; Zhang & Cheng, 1997; Cheng et al., 2000; Zhang & Cheng, 2001; Zhang et al., 2004) models. These gap models have achieved great successes in explaining pulsar high-energy emissions and light curves (e.g., Watters el al., 2009; Romani & Watters, 2010). However, the vacuum solution has no plasma, it is not able to reproduce any pulsar phenomenons. It is well known that the pulsar magnetosphere should be filled with plasma (Goldreich & Julian, 1969). In the presence of abundant plasma, all accelerating electric fields can be efficiently screened to form a force-free (FF) magnetosphere. The FF solution for an aligned rotator was first obtained by Contopoulos et al. (1999, hereafter CKF). The CKF solution consists of a closed field line region extending to the LC, an open field line region, and an equatorial current sheet beyond the LC. Moreover, the time-dependent simulations for the FF axisymmetric rotator also confirmed the closed-open CKF solution (e.g., Komissarov, 2006; McKinney, 2006; Timokhin, 2006; Yu, 2012; Parfrey et al., 2012; Cao et al., 2016a; Etienne et al., 2016).

Three-dimensional (3D) structures of the force-free pulsar magnetosphere have been available in recent years (e.g., Spitkovsky, 2006; Kalapotharakos & Contopoulos, 2009; Pétri, 2012). Moreover, these studies have been also extended to the general-relativistic force-free regime that takes time-space curvature and frame-dragging effects into accounts (Pétri, 2016; Carrasco et al., 2018) and the full magnetohydrodynamic (MHD) regime that takes plasma pressures and inertial into account(Tchekhovskoy et al., 2013). The 3D force-free solutions are similar to the CKF solution with an equatorial current sheet outside the LC. These simulations gave an impetus to the study of pulsar high-energy radiation with a more realistic field geometry instead of the vacuum dipole field. Gamma-ray light curves have been modeled using the force-free magnetic field (Bai & Spitkovsky, 2010b; Bogovalov et al., 2018). The force-free solutions provide different pulse profiles due to the increased polar cap size compared to that of the vacuum dipole. It should be noted that the force-free solutions cannot allow any accelerating electric fields. Therefore, they cannot provide any information about the sites of particle acceleration and radiation.

It is well known that the vacuum solution has no any particle distributions, while the force-free solution does not allow particle acceleration along magnetic field lines. Therefore, more realistic pulsar magnetosphere should lie between the vacuum and force-free limits. In fact, the plasma resistivity can produce a none-zero accelerating electric field along the magnetic field, and a sets of resistive solutions that smoothly bridges the gap between the vacuum and force-free solutions have been constructed based on the finite-difference time-domain approach (Li et al., 2012; Kalapotharakos et al., 2012a). Moreover, the resistive pulsar magnetospheres have been used to model the pulsar $\gamma$-ray spectra and light curves (Kalapotharakos et al., 2014; Brambilla et al., 2015). Also, particle-in-cell (PIC) methods with a self-consistent treatment between particles and fields are developed to simulate the structures of the pulsar magnetosphere (Philippov & Spitkovsky, 2014; Chen & Beloborodov, 2014; Belyaev, 2015; Cerutti et al., 2015; Philippov et al., 2015; Kalapotharakos et al., 2018; Brambilla et al., 2018). Very recently, full PIC simulations also start to predict the pulsar light curves by including the radiation reaction (Cerutti et al., 2016; Philippov & Spitkovsky, 2018; Kalapotharakos et al., 2018).

In the previous paper, we presented the structures of oblique pulsar magnetosphere with the uniform conductivity by a pseudo-spectral method (Cao et al., 2016b). We constructed a set of resistive solutions that smoothly bridges the gap between the vacuum and force-free limits. In this paper, we focus on the the influence of the conductivity on gamma-ray light curves in dissipation pulsar magnetospheres. The light curves are produced by the geometric method and the particle trajectory method. As an application, we compare the predicted light curves with the gamma-ray light curves of Vela pulsar observed by Fermi-LAT. The paper is organized as follows: In section 2, we present the basic equations describing the structure of pulsar magnetosphere. In section 3 and 4, we present the method that we have used to produce the gamma-ray light curves. In section 5, we apply the model to Vela pulsar. Finally, a brief discussion and conclusions are given in section 6.

2. Basic Equations

The time-dependent Maxwell equations are given by

[TABLE]

with two initial conditions

[TABLE]

where $\rho_{\rm e}$ is the charge density and ${\bf J}$ is the current density. The structure of pulsar magnetosphere can be determined by the prescription for the current density ${\bf{J}}$. For the force-free magnetosphere, the current density can be derived by force-free condition ($\bf E\cdot\bf B=0$) and the Maxwell equations as (Gruzinov, 1999; Blandford, 2002)

[TABLE]

which implies that the electric field parallel to the magnetic field component ${\bf{E}}_{\|}$ will be zero.

The force-free approximations do not allow for any particle acceleration along the magnetic field, ${\bf{E}}_{\|}=0$. Therefore, they do not account for the particle acceleration and pulsed emission in the magnetosphere. More realistic pulsar magnetosphere should allow for ${\bf{E}}_{\|}\neq 0$ in some regions of the magnetosphere. The resistive magnetosphere can allow for a non-zero parallel electric field component by involving a conductivity parameter $\sigma$. In fact, the resistive pusar plasma is effectively collisionless. Therefore, the resistivity coming from coulomb collisions can be practically neglected. We can define the current density for the resistive magnetosphere by a form of Ohm’s law given by (Kalapotharakos et al., 2014)

[TABLE]

where

[TABLE]

The first term in equation (6) is the drift current component perpendicular to ${\bf{B}}$, while the second term controls ${\bf{E}}_{\|}$ by the conductivity $\sigma$. The $E_{0}$ term in equation (6) ensures the drift current to be subluminal. When the conductivity $\sigma$ increases from $\sigma=0$ to $\sigma\rightarrow\infty$, we expect to obtain a set of solutions that smoothly transition from the vacuum field to the force-free field. It is noted that there is no unique prescription for the current density in the resistive electrodynamics. We checked the field structures by using the prescription for the current density given by Li et al. (2012). Our results are also qualitatively very similar to those of Li et al. (2012).

3. Geometric Method

3.1. Magnetic Field Structures

The time-dependent Maxwell equations are solved by a spectral method in spherical coordinates. The radial coordinate $r$ is expanded into the Chebyshev function. However, the angle coordinates $\theta$ and $\phi$ are expanded into the vector spherical harmonic functions. A spectral filter is used in all directions in order to ensure the stability of the algorithm and increase the convergent rate of the solution. For a detailed discussion about the pseudo-spectral algorithm, see Pétri (2012) and Cao et al. (2016a, b). The computational domain extends from the inner boundary $r_{\rm min}=0.2\,r_{\rm L}$ to the outer boundary $r_{\rm max}=3\,r_{\rm L}$. A resolution of $N_{r}\times N_{\theta}\times N_{\phi}=128\times 32\times 64$ is necessary to get a good accuracy. The magnetic field is initialized to be an rotating vacuum dipolar with magnetic inclination angle $\alpha=\{0^{\circ},15^{\circ},30^{\circ},45^{\circ},60^{\circ},75^{\circ},90^{\circ}\}$. We impose inner boundary condition at the stellar surface with a rotating electric field ${\bf{E}}=-({\bf\Omega}\times{\bf r})\times{\bf B}/c$. A characteristic compatibility method is used in order to prevent the inward reflection from the outer boundary. We let the system evolve for several rotational periods to reach a final steady solution. We perform a series of simulations with uniform conductivities $\sigma=\{0.1,\,0.3,\,1,\,3,\,5,\,10,\,30,\,60\}\,\,\Omega$. Magnetic field structures in the $\Omega-\mu$ plane for magnetic inclination $\alpha=60^{\circ}$ with increasing conductivity are show in figure 1. The normalized Poynting flux $L/L_{\rm aligned}$ as a function of radius $r$ for magnetic inclination $\alpha=60^{\circ}$ with different $\sigma$ values at time $t=3\,P$ are shown in figure 2. We can see that the Poynting flux increases with increasing conductivity. To demonstrate the convergence of the solution, the normalized Poynting flux for $\sigma=30\,\Omega$ at time $t=4\,P$ is also shown as the red dashed line. The Poynting flux is nearly constant when the time $t>3\,P$, the system will relax to the steady-state solution. In figure 3, we show the normalized Poynting flux $L/L_{\rm aligned}$ as a function of radius $r$ for magnetic inclination $\alpha=0^{\circ}$ with different $\sigma$ values at time $t=3\,P$. It is found that our results still have some numerical dissipations due to the low resolution used in the presented 3D simulations. A higher resolution is necessary to resolve more fine plasmoids in the current sheet. Our simulations only present the lower limit on resolution to get convergence. We can see a smooth transition for magnetic field structures and the Poyting flux from the vacuum limit to the force-free limit with increasing conductivity. For a detailed description about the structures of resistive magnetosphere, see Cao et al. (2016b).

The polar cap shape is defined by the footprints of the last closed field lines on the stellar surface. The last closed field lines can be found by checking whether the field lines close inside or outside the LC. It has been known that the polar cap shape in vacuum case is different from that in the force-free case (e.g., Bai & Spitkovsky, 2010b; Kalapotharakos et al., 2014). Here, a third-order Runge-Kutta integration is used to find the last closed field lines in the lab frame, and the bisection method is used to find the magnetic colatitude $\theta^{\rm rim}_{m}$ of the rim of the polar cap for the fixed magnetic azimuth $\phi^{\rm rim}_{m}$. Open volume coordinates are be defined by $(r_{\rm ovc},\phi_{\rm m})$ on the polar cap. $r_{\rm ovc}\geq 1$ represent the closed field lines and $r_{\rm ovc}<1$ represent the open field lines. The polar cap shapes on a sphere of radius 0.2 $r_{\rm L}$ for magnetic inclination $\alpha=60^{\circ}$ with different values of $\sigma$ are shown in figure 4. The polar cap shapes for the low $\sigma$ value are similar to that for the vacuum field, and a notch appears on the polar cap rim. In the resistive solution, on the one hand, the polar shapes are larger than that of vacuum solution because of the larger open magnetic flux. On the other hand, the polar caps are shifted towards the trailing side. The notches on the polar caps disappear at high $\sigma$ value and the polar cap (PC) shapes are more circular.

3.2. Sky Maps and Light Curves

The light curves can be computed by the geometry of pulsar emission models for a given magnetosphere structure. In the SG model, all accelerating electric fields are screened by the pair cascades above the pair formation fronts. There is a narrow gap with $E_{\|}\neq 0$ between the last closed lines and the pair-formation fronts. The particles continue to accelerate and radiate from the neutron star surface to high altitude. The OG model is a vacuum gap that also forms between the null charge surfaces and the last closed field lines. The outer gap extends from the null-charge surfaces to the light cylinder, and the gap width is limited by the screening of $E_{\|}\neq 0$ by pair cascades. After the current through the gap is taken into account, it has been shown that the inner boundary of the outer gap is shifted inward from the null charge surface (Takata et al., 2004; Hirotani, 2006). Therefore, the gaps would exist inside the light cylinder of the pulsar magnetosphere, the details of $E_{\|}$ depend on the pair cascade process which is not taken into account in the simulation of pulsar magnetosphere so far.

To explore how the magnetic field structures and offset PCs influence $\gamma$-ray pulsar light curves, we use a geometric method to produce the $\gamma$-ray light curves based on the uniform $\sigma$ model. We assume that the $\gamma$-ray emission originates in a layer of width $\omega=0.1$ along the last closed field lines between $r_{\rm ovc}=0.9$ and $r_{\rm ovc}=1$. The emission region extends from the neutron star surface to the light cylinder, and the $\gamma$-ray emission is restricted to lie within a cylindrical radius $\rho_{\rm max}\lesssim r_{\rm L}$. We assume the uniform emissivity along the field lines in the CF. The photon direction is assumed to be tangent to the magnetic field in the CF, obtained through a Lorentz transformation from the inertial observer’s frame (Bai & Spitkovsky, 2010a). The emitted photon is collected in the sky map in viewing angle $\zeta$ and the observed phase $\Phi$, taking into account the aberration and time-delay effects. The light curves are then obtained by a cut through this sky map at constant $\zeta$.

In figure 5, we show the sky maps and the corresponding light curves from the geometric method in different inclination angles and view angles with different magnetic field configurations. We see narrow double-peak profiles at relative large $\alpha$ and $\zeta$, which are very similar to the observed $\gamma$-ray light curves. We see that the light curves for the low $\sigma$ value are very close to that of the vacuum dipole, but the peaks are shifted in phase. As the conductivity $\sigma$ increases, the peak phase is shifted to the lager phase relative to the magnetic pole. The similar results are also found by Harding (2011) and Kalapotharakos et al. (2012b). Overall, there is a significant progression in the LC shapes as $\sigma$ increases. This is because that the magnetospheres with high conductivity have more sweepback field lines, which produces the larger shift of PC and thus causes a lager phase lag of the light curve peaks.

4. Particle Trajectory Method

4.1. Real Particle Trajectories

We study the contribution of the curvature radiation to the $\gamma$-ray light curves by including the accelerating electric fields provided by the solutions themselves. We define the trajectory of particles in inertial observer’s frame ( IOF ) by (Kalapotharakos et al., 2014)

[TABLE]

where the first term in the equation (8) is a drift velocity component, while the second term is a component parallel to the magnetic field. By requiring that the condition $v\simeq c$ and that the motion of the particle is outward, the spatial distribution of $f$ can be uniquely determined by equation (9). We assume that the direction of photon emission, ${\bm{\eta}}_{\rm em}$, is along the direction of particle motion ${\bm{\beta}}={\bf v}/c$. Therefore, the direction of photon emission in the IOF is determined by

[TABLE]

where the view angle $\zeta=\rm acos(\mu_{\rm em})$.

The observe phase is determined by include the rotation and time-delay correction

[TABLE]

where $\phi_{\rm rot}=\Omega\,dt$ is the rotation phase, $\phi_{\rm em}$ is the azimuthal angle of the emitted photon, and ${{\bf r}_{\rm em}}$ is the location of the emitted photon.

Under above assumptions, we can calculate the trajectories of radiating particles passing through each magnetospheric point in the computational domain, which also allows the determination of curvature radius $R_{\rm CR}$ along each particle trajectory. We assume that the charge particles $(e^{-}\rm,e^{+})$ are uniformly distributed on the polar cap. The $e^{\pm}$ pairs are then injected from the stellar surface with small Lorentz factor ($\gamma<100$). The Lorentz factor $\gamma$ of radiating particle along each trajectory is calculated by including the influence of the accelerating electric field and curvature radiation loss. Specifically, the $\gamma$ value along each trajectory is integrated by the expression

[TABLE]

where $q_{\rm e}$ and $m_{\rm e}$ are the electron charge and rest mass, respectively. The first term in equation (11) is the energy gain rates of the particles due to the accelerating electric field and the second term is the energy loss rates due to curvature radiation. Then, we calculate the characteristic energy of curvature radiation $E_{\rm c}=\frac{3}{2}c\hbar\frac{\gamma^{3}}{R_{\rm CR}}$ and its bolometric luminosity $L_{\rm bol}=\frac{2}{3}q^{2}_{\rm e}c\frac{\gamma^{4}}{R^{2}_{\rm CR}}$ along each trajectory. We integrated each particle trajectory from the neutron star surface up to $r=2.5\,r_{\rm L}$. We construct the $\gamma$-ray light curve by collecting the bolometric emission with $E_{\rm c}>0.1\,\rm GeV$ from all the emitting particles in sky maps. For a detail description about the trajectory and radiation of particles, see Kalapotharakos et al. (2014).

4.2. Sky Maps and Light Curves

Kalapotharakos et al. (2014) studied the $\gamma$-ray light curves using dissipative pulsar magnetospheres. It is found that the uniform $\sigma$ model cannot well match the distribution of $\gamma$-ray peak separation and the radio lag observed by Fermi-LAT. Further, they found that a significant improve can be achieved by applying a force-free regime inside the LC and a dissipative regime (FIDO) outside the LC. In fact, we also find that the uniform $\sigma$ model cannot explain the pulsar light curves observed by Fermi LAT. Therefore, we adopt a model that is similar to the FIDO model introduced by Kalapotharakos et al. (2014). Kalapotharakos et al. (2014) used a $\sigma$ linear approximation for $E_{\|}$ based on the force-free solution at high $\sigma$ value. In this work, we adopt a more accurate treatment for $E_{\|}$. We employ a very high conductivity ($\sigma=60\,\Omega$) within the LC and a finite conductivity outside the LC, we then remove any $E_{\|}$ solutions within the LC. We assume that all the $\gamma$-ray emission comes from the outer magnetosphere outside the LC. In figure 6, we shown the $E_{\|}$ distribution in the $\Omega-\mu$ plane for magnetic inclination $\alpha=60^{\circ}$ with $\sigma=60\,\Omega$ within the LC and with $\sigma=30\,\Omega$ outside the LC. The field line structure is very similar to that of the force-free solution. We see the strong parallel electric field along the equatorial current sheet outside the LC. We also note that the parallel electric distribution is very similar to that of the global PIC simulation (Kalapotharakos et al., 2018; Philippov & Spitkovsky, 2018).

In the following study, we assume the standard pulsar parameters for pulsar period $P=0.1\,\rm s$ and surface magnetic field $B_{\star}=10^{12}\,\rm G$. We calculate the direction of photon emission and the corresponding bolometric emission along each particle trajectory. This information allow us to produce $\gamma$-ray light curves by collecting all the photons with $E_{\rm c}>0.1\,\rm GeV$ in sky maps. In figure 7$-$8, we show the sky maps and the corresponding light curves in different inclination angles and observer viewing angles for the magnetosphere with a near force-free regime within the LC and with $\sigma=0.3\,\Omega$ and $30\,\Omega$ outside the LC. For $\sigma=0.3\,\Omega$ outside the LC, the light curves generally display only one broad peak. We see the double-peak light curves as $\alpha$ and $\zeta$ increase, but the peaks seem to be broad. For very high $\sigma$ value outside the LC, we see the significant double-peak profiles except for some additional secondary peaks. There are some cases where the light curves are narrow, which are very similar to those observed by Fermi LAT. We find that all the emission originates in a thin layer with $r_{\rm ovc}=0.9-1$ near the equatorial current sheet outside the LC.

5. Application to Vela Pulsar

The observed pulsar light curves potentially provide an important diagnostic for pulsar magnetosphere structure. As an application, we compare the light curves predicted from the two methods with those observed at $>0.1\,\rm GeV$ energies for Vela pulsar. The results are shown in figure 9 and 10. For the geometric method, we find that the predicted light curve with $\sigma=0.3\,\Omega$ provides a better match to the Fermi observed data for Vela pulsar. The double-peak profile with phase separation of $\sim$ 0.43 can be better reproduced by our model with the low value of $\sigma$. In fact, the first peak phase lag of the high $\sigma$ model is too large to explain the observed data for Vela pulsar. We suggest that the magnetosphere with the low $\sigma$ value is more favorable for the geometric method. For the particle trajectory method, we use the measured parameters of Vela pulsar: $P=0.089$ s, $B=3.4\times 10^{12}$ G. We find that the magnetosphere with a near force-free regime within the LC and with $\sigma=30\,\Omega$ outside the LC can better produce the observed data of Vela pulsar. Our results indicate that the magnetosphere with a near force-free regime within the LC and a high $\sigma$ value outside the LC is more favorable for the particle trajectory method.

6. Discussion and Conclusions

In this paper, we explore the $\gamma$-ray emission patterns and light curves in dissipation pulsar magnetospheres. The $\gamma$-ray light curves are produced by the geometric method and particle trajectory method. For the geometric method, we assume that the $\gamma$-ray emission comes from a finite-width layer along the last closed lines, we generate the $\gamma$-ray light curves by assuming the uniform emissivity along the field lines in the CF. We find that there is a significant progression in the $\gamma$-ray peak phase with increasing $\sigma$. For the particle trajectory method, we consider the spatial distribution of $\sigma$ by assuming a very high conductivity within the LC and a finite conductivity outside the LC. Then, we use the field structure of these models to define realistic trajectories of radiated particles. Assuming that all the $\gamma$-ray emission comes from the outer magnetosphere outside the LC, we compute the Lorentz factors of radiating particles and the characteristic energy of curvature radiation along each particle trajectory under the influence of both the accelerating electric field and curvature radiation loss. We produce the $\gamma$-ray light curves by collecting the bolometric emission with $E_{\rm c}>0.1\,\rm GeV$ from all the emitting particles. We find that the light curve shapes are very sensitive to the value of $\sigma$, and all the $\gamma$-ray emission is produced in regions near the equatorial current sheet for very high $\sigma$ value.

As an application, we compare the modeling light curves from the two methods with observed light curves at $>0.1\,\rm GeV$ energies for Vela pulsar. For the geometric method, we find that the magnetosphere with the low $\sigma$ value provides a better match to the observed data for Vela pulsar. For the particle trajectory method, we find that the magnetosphere with a near force-free regime within the LC and a high $\sigma$ value outside the LC can better produce the light curves of Vela pulsar. We note that the particle trajectory method uses the parallel electric fields that are self-consistent with the magnetic fields of the magnetosphere. Therefore, the results from the particle trajectory method are more physically motivated. In fact, the microphysical processes of the pair production and particle acceleration are not self-consistently included in the simulation of the resistive magnetosphere. However, the parallel electric field distribution from the high $\sigma$ model outside the LC is only confined near the equatorial current sheet, which is very similar to that of the global PIC simulation (Kalapotharakos et al., 2018; Philippov & Spitkovsky, 2018). This gives a microphysical explanation for the origin of the parallel electric field from the high $\sigma$ model outside the LC. Our results from the particle trajectory are also similar to those from the PIC simulation. The particle trajectory method also allows the derivation of energy spectrum, which can be compared directly with the Fermi observations. Therefore, we will explore the $\gamma$-ray energetic and spectral properties in future work. Further, we will improve our code by identifying different field lines at each time-step of simulation. This will allow us apply $\sigma$ along specific magnetic field lines. This technique will be helpful for us to understand the underlying pair production process in the gaps or current sheets.

We thank the anonymous referee for valuable comments and suggestions. We would like to thank Jér$\hat{\rm o}$me Pétri, Gabriele Brambilla, Xuening Bai and Li Zhang for some useful discussions. We acknowledge the financial support from the National Natural Science Foundation of China 11573060 and 11661161010, the National Science Foundation of China 11673060, the National Science Foundation of China 11871418 and Fund for Science and Technology Innovation team in universities at Yunnan. The simulation presented in this work is performed at HPC Center, Yunnan Observatories, CAS, China.