The effect of an offset-dipole magnetic field on the Vela pulsar’s -ray light curves
Over the past six years, the Fermi Large Area Telescope has detected more than 150 -ray pulsars, discovering a variety of light curve trends and classes. Such diversity hints at distinct underlying magnetospheric and/or emission geometries. We implemented an offset-dipole magnetic field, with an offset characterised by parameters epsilon and magnetic azimuthal angle, in an existing geometric pulsar modelling code which already includes static and retarded vacuum dipole fields. We use these different -field solutions in conjunction with standard emission geometries, namely the two-pole caustic and outer gap models (the latter only for non-offset dipoles), and construct intensity maps and light curves for several pulsar parameters. We compare our model light curves to the Vela data from the second pulsar catalogue of Fermi. We use a refined chi-square grid search method for finding best-fit light curves for each of the different models. Our best fit is for the retarded vacuum dipole field and the outer gap model.
The discovery of the first pulsar in 1967 by Bell and Hewish  gave birth to pulsar astronomy. Pulsars are believed to be rapidly-rotating, compact neutron stars that possess strong magnetic, electric, and gravitational fields . They emit radiation across the entire electromagnetic spectrum, including radio, optical, X-ray, and -rays . Since the launch of the Fermi Gamma-ray Space Telescope in June 2008, over 150 -ray pulsars have been detected, of which the Crab and Vela pulsars are the brightest sources. Fermi consists of two parts including the Large Area Telescope (LAT) and the Gamma-ray Burst Monitor. The LAT measures -rays in the energy range between 20 MeV and 300 GeV . Over the past six years, Fermi has released two pulsar catalogues, both describing the light curve profiles and spectral characteristics of -ray pulsars [1, 2]. The light curves show great variety in profile shape, and may be divided into three general classes based on the relative phase differences between the radio and -ray pulses [11, 24]. The light curves also show energy-dependent behaviour. Most of the young and millisecond pulsars exhibit two -ray peaks whereas some pulsars including Vela display three peaks .
Models of pulsar geometry are characterised by the inclination angle () between the rotation () and magnetic () axes, the observer’s viewing angle () between the rotation axis and the observer’s line of sight, and the impact angle (). Geometric models assume the presence of several ‘gap regions’ in the pulsar magnetosphere. These are defined as regions where particle acceleration and emission take place. These geometric models include the two-pole caustic (TPC)  (the slot gap (SG)  model may be its physical representation), outer gap (OG) [6, 22], and pair-starved polar cap (PSPC) model . All of these models are represented in figure 1. The emissivity of high-energy photons within this gap region is assumed to be uniform in the corotating frame for geometric models (for physical models the changes with radial distance when an electric field is assumed). Since the -rays are expected to be emitted tangentially to the local magnetic field in the corotating frame , the assumed magnetic field geometry is very important with respect to the predicted light curves. Several magnetospheric structures have been studied, including the static dipole field , the retarded vacuum dipole field (RVD)  and the offset-dipole -field. The latter is motivated by the fact that retardation of the -field at the light cylinder causes offset of the polar caps (PCs), always toward the trailing edge of the PC () . The offset is characterised by parameters epsilon () and magnetic azimuthal angle () which represents a shift of the PC away from the magnetic axis, with corresponding to the static-dipole case.
In this paper we studied the effect of using different combinations of magnetospheric structures, geometric models, and model parameters on -ray light curves. In Section 2 we discuss the implementation of an offset-dipole solution and associated electric field. Section 3 describes how we matched limiting cases of the low-altitude and high-altitude -fields using a matching parameter . In Section 4 we describe the chi-squared () method we used to search for best-fit light curves. Our results are given in Section 5 and the conclusions follow in Section 6.
2 Implementation of an offset-dipole -field
We implemented  an offset-dipole magnetic field  into an existing geometric pulsar modelling code  which already includes static and RVD fields. The implementation involves a transformation of the -field from spherical to Cartesian coordinates, rotating both the -field components and its Cartesian frame through an angle , thereby transforming the -field from the magnetic frame () to the rotational frame (). We extended the range of for which we could solve the PC rim (for details, see ) by enlarging the colatitude range thought to contain the last open field line (tangent to , see figure 1).
3 Matching parameter
It is important to take the accelerating -field into account (in a physical model) when such expressions are available, since this will modulate the in the gap (as opposed to geometric models where we just assume constant at all altitudes). We use analytic expressions  for a low-altitude and high-altitude SG -field in an offset-dipole magnetic geometry. These are matched at a critical scaled radius to obtain a general -field valid for all altitudes . In previous work we chose . In this paper we solve on each field line. Using the general -field we could solve the particle transport equation [23, 7] (taking only curvature losses into account) to obtain the particle Lorentz factor that is necessary for calculating the curvature emissivity.
4 Finding best-fit light curves
We tested several fitting methods in order to select the most suitable one. We decided to use the standard method. For each combination of -field and geometric model, we calculate for each set of free model parameters , , normalisation () and phase shift (). We assume that for the bright Vela pulsar, with a large amount of counts in each bin, the will be Gaussian distributed with the number of degrees of freedom. For the standard Gaussian distribution, assuming very small Gaussian errors, the reduced values are very large and therefore we needed to scale the values with the optimal and multiply by . The scaled is denoted by , and the reduced scaled will have a value of . We next determined confidence intervals (1, 2, and 3) in (,) space for these model parameters, using . We only fitted the -ray light curve, because we do not want to bias our results with a simplistic radio model.
5.1 Phaseplots and light curves
In figure 3, 3 and 4, we show the intensity maps or phaseplots (emission per solid angle versus and pulse phase ) and their corresponding light curves (i.e., cuts at constant ) for the offset-dipole -field and the TPC model. The dark circle in figure 3 is the non-emitting PC, and the sharp, bright regions are the emission caustics, where radiation is bunched in phase due to relativistic effects. Figure 3 and 3 represent phaseplots for a fixed and , with ranging from 0 to 0.18 with increments of 0.03. We contrasted the cases of constant and variable . We observed a qualitative difference in caustic structure. The caustics seem larger and more pronounced in the constant case. In figure 4 we chose a fixed value of for variable , with ranging from 0 to 90 and from 15 to 90, both with a resolution of 15. This shows examples of various light curves that may be obtained in this model.
5.2 Contours and best-fit light curves
In figure 5 we represent our best fits we obtained using the method. Our overall best fit is for the RVD field and the OG model, with , , , and . For the offset-dipole solution and the TPC model, we have a best fit (assuming constant ) for parameters , , , , and . The best fit for the offset-dipole solution and the TPC model, assuming variable , is for parameters , , , , and . The best-fit parameters for each -field and geometric model combination are summarised in table 1. The table includes the different model combinations, the optimal value (before scaling), the free parameters with errors (found using 3 connected (,) contours), a reference fit found using radio polarisation data, and the comparison between models using the difference between the respective optimal values of , represented by .
|\br Model||Our model parameters||Ref. fit ||Radio pol. |
|\mr Static dipole:|
|\mr Offset dipole - constant :|
|\mr Offset dipole - variable :|
We have studied the effect of implementing the offset-dipole -field on -ray light curves for the TPC geometry. We find an optimal best-fit for Vela for the RVD -field and the OG model. The OG model displays reduced off-peak emission. We note that the best fits for the offset-dipole -field for constant favour smaller values of and for variable larger values. When including an -field the resulting phaseplots becomes qualitatively different compared to constant . In future, we want to continue to produce light curves using improved geometric models and -fields, and also using more data, in order to search for best-fit profiles, thereby constraining the low-altitude magnetic structure and system geometry of several bright pulsars.
This work is supported by the South African National Research Foundation. AKH acknowledges the support from the NASA Astrophysics Theory Program. CV, TJJ, and AKH acknowledge support from the Fermi Guest Investigator Program.
-  Abdo A A et al. 2010 ApJS 187 460–94
-  Abdo A A et al. 2013 ApJS 208 1–59
-  Atwood W B et al. 2009 ApJ 697 1071–102
-  Becker W, Gil J A, and Rudak B 2007 Highlights of Astronomy 14 109–38
-  Breed M, Venter C, Harding A K and Johnson T J 2013 58th Ann. Conf. of the SA Inst. of Phys. in press
-  Cheng K S, Ho C and Ruderman M 1986 ApJ 300 500–39
-  Daugherty J K and Harding A K 1996 ApJ 458 278–92
-  Deutsch A J 1955 Annales d’Astrophysique 18 1–10
-  Dyks J and Rudak B 2003 ApJ 598 1201–6
-  Dyks J, Harding A K and Rudak B 2004 ApJ 606 1125–42
-  Espinoza C M et al. 2013 MNRAS 430 571–87
-  Goldreich P and Julian W H 1969 ApJ 157 869–80
-  Griffiths D J 1995 Introduction to Electrodynamics (San Francisco: Pearson Benjamin Cummings) p 246
-  Harding, A K 2004 22nd Texas Symp. on Relativistic Astrophysics ed P Chen, E Bloom et al. p 40
-  Harding A K, Usov V V and Muslimov A G 2005 ApJ 622 531–43
-  Harding A K and Muslimov A G 2011 ApJ 743 181–96
-  Hewish A, Bell S J, Pilkington J D H, Scott P F and Collins R A 1968 Nature 217 709–13
-  Johnston S 2005 MNRAS 364 1397–412
-  Muslimov A G and Harding A K 2003 ApJ 588 430–40
-  Muslimov A G and Harding A K 2004 ApJ 606 1143–53
-  Pierbattista M et al. 2014 Astronomy & Astrophysics in press
-  Romani R W 1996 ApJ 470 469–78
-  Sturner S J 1995 ApJ 446 292–99
-  Venter C, Harding A K and Guillemot L 2009 ApJ 707 800–22
-  Watters K P, Romani R W, Weltevrede P and Johnston S 2009 ApJ 695 1289–301