Gamma-ray emission from millisecond pulsars - an Outergap perspective
Gamma-ray emission from millisecond pulsars - an Outergap perspective
Journal of Astronomy and Space Sciences. 2013. Sep, 30(3): 153-158
Copyright ©2013, The Korean Space Science Society
This is an open Access article distributed under the terms of theCreative Commons Attribution Non-Commercial License ( which premits unrestrictednon-commercial use, distribution, and reproduction in any medium,provided the original work is properly cited.
  • Received : November 11, 2012
  • Accepted : December 12, 2012
  • Published : September 15, 2013
Export by style
Cited by
About the Authors
Kwong Sang Cheng
In this review paper we explain the following gamma-ray emission features from the millisecond pulsars. (1) Why is the dipolar field of millisecond pulsars so weak but the magnetic pair creation process may still be able to control the size of the outergap? (2) A sub-GeV pulse component could occur in the vicinity of the radio pulse of millisecond pulsars. (3) Orbital modulated gamma-rays should exist in the black widow systems for large viewing angle.
Pulsars are known to be rapidly spinning and strongly magnetized neutron stars so they can behave like a unipolar inductor. Without the charge screening young pulsars can easily develop huge potential drop along the open field lines (cf. Manchester & Taylor 1977, Lyne & Graham-Smith 1998 or a general review). Such large potential drop can accelerate charged particles to extremely relativistic speed. Since the magnetic field in the open field lines are still sufficiently strong so charged particles are confined to move along the field lines, hence high energy curvature photons are emitted. It has been shown that pulsar magnetosphere is filled with charge separated plasma, whose charge distribution is given by the Goldreich-Julian charge density
Lager Image
, where
Lager Image
is the angular velocity vector and
Lager Image
is the local magnetic field vector (Goldreich & Julian 1969).
Lager Image
is perperdicular to
Lager Image
, the charge density is zero and such region is called null charge surface and charged carriers on two sides of the null charge surface are opposite. With this charge distribution the electric field along the magnetic field is screened out. However Cheng et al. (1986a) argued that when the current pass through the null charge surface, opposite charged carriers are removed from the vicinity of null charge surface and hence vacuum region will be formed. Charged particles in this vacuum region will be accelerated to extremely relativistic speed so this region is also called “Outergap” accelerator. The question is if the current continue to pass through the null charge region will the vacuum region to grow unlimitedly. The answer is “no” because the energy of the curvature photons is proportional to the size of vacuum region. When the energy of curvature photons is energetically enough, curvature photons and soft photons from the stellar surface can be converted into electron/positron pairs, which can limit the size of the Outergap. Zhang & Cheng (1997) argued that the soft photons from the surface should be dominated by thermal photons emitted by the hot polar cap, which is heated by the return current from the Outergap. In this particular model, the fractional size of the Outergap ( fpp ) limited by the photonphoton pair creation process is completely determined by the global parameters of pulsars, i.e. rotation period P and the dipolar magnetic field B , and it is given by
Lager Image
where B 12 is the dipolar magnetic field in units of 10 12 G and the total gamma-ray luminosity emitted by the Outergap is given by
Lager Image
where fgap is the fractional size of the Outergap and Lsd is the spin-down power of the pulsar.
It has been shown that many observed features from gamma-ray pulsars can be explained in terms of the Outergap models (e.g. Cheng et al. 1986b, 2000, Romani 1996, Zhang & Cheng 1997, Takata et al. 2004, Hirotani 2006, 2008, Wang et al. 2010)
Takata et al. (2010) suggest that if the surface multiple field is much stronger than the dipolar field near the stellar surface in particular they can bend those open field lines connecting to the Outergap sideward, then the curvature photons emitted by the return current from the Outergap can become pairs near the surface and their pitch angles can be larger than 90 degrees(cf. Fig. 2 of Takata et al. 2010). In this situation, the magnetic pairs can stream back to the outermagnetosphere and restrict the size of the Outergap. The Outergap size restricted by the magnetic pair creation process is given by
Lager Image
where P -1 is the rotation period in units of 10 -1 s, K ( Bm , s) depends on the local multiple magnetic field ( Bm ) and the local curvature radius ( s ). By fitting the first Fermi pulsar catalog (Abdo et al. 2010), they find that K is roughly a constant but the value is unity for canonical pulsars and is about 10 for millisecond pulsars respectively. These two values correspond to quadruple field strength Bm ~10 13 G for canonical pulsars and Bm ~10 11 G for millisecond pulsars respective. They also show that the magnetic pair creation process is more important for pulsars with spin-down power lower than 10 36 erg/s. Although the existence of multiple field about 10 times of the dipolar field is possible (e.g. Arons 1983), why could old neutron stars like millisecond pulsars, which is likely older than billion years old, can still maintain a surface multiple field hundred to even thousand times stronger than the dipolar field? Ruderman (1991) argue that in the core of the neutron star the interpinning between the quantized magnetic flux tubes and the quantized rotation vortex lines is very strong, the spin-down and spin-up of the neutron star can cause the vortex lines to move outward and inward respectively and hence the flux tubes move together with the vortex lines. During the accretion spinup epoch of the millisecond pulsars the core magnetic field, which cannot decay due to interpinning with the vortex lines, will be squeezed in the rotation axis due to the dragging of vortex lines. Due to the conservation of magnetic flux the local magnetic field strength will increase but this squeezed magnetic field cannot be stronger than the critical field Bc ~10 15 G otherwise the quantized flux tube structure will break down then the magnetic field lines will be reconnected to reduce the local magnetic field strength to below the critical field strength. The structure of this squeezed magnetic field can be illustrated as follow. The vortex lines will drag the north (south) poles of the flux tubes arrive the rotation axis first then the south (north) poles will follow (cf. Fig. 2 of Chen & Ruderman 1993). But these north and south poles cannot be reconnected because underneath these two magnetic poles they are separated by the vortex lines. Therefore the north will cluster near the rotation axis surrounded by the south poles in the perimenter, a quadruple magnetic field structure is formed. We can estimate the strength of quadruple moment as Q ~ Bc a 4 , where Bc is the critical field strength and a is the characteristic size of the quadruple moment, which can be estimated by the conservation of vortex lines
Lager Image
, where Ω is the angular velocity, r is the radius from the rotation axis and κ is the quantized vorticity. We obtain
Lager Image
where R 6 is the core radius in units of 10 6 cm, Pi is the initial period before spin-up and Pf ,-3 is the final period after spinup in units of 10 -3 s. It is very important to note that this quadruple moment is located at the boundary between the core and the crust. The magnetic field produced by this quadruple moment on the surface of the neutron star is given by
Lager Image
where l is the thickness of the stellar crust.
Cheng et al. (2000) have shown that most pairs should be created around the null charge surface. In this Outergap model, the current flow inside the Outergap is dominated by the outflow current from the null charge surface toward the light cylinder whereas the inflow current dominates from the null charge surface toward the star. It has been shown by Hirotani (2005) that the inner boundary of the Outergap is not located at the null charge surface when the outergap current is not zero. Roughly speaking the location of the inner boundary can be estimated as
Lager Image
where rin and rnull are the distances to the inner boundary and null charge surface respectively, and Jgap and JGJ are the gap current and the Goldreich-Julian Current respectively. The electric field between the null charge surface and the inner boundary can be approximated (Hirotani 2005, Tang et al. 2008) by
Lager Image
where E ( rnull ) is the electric field at the null charge surface. Since the electric field from the null charge surface to the light cylinder is roughly constant whereas the electric field from the null surface to the inner boundary decreases quadratically therefore the characteristic curvature photon energy and the radiation power by the inward current flow are weaker than those of the outflow. The average electric field in the inflow region is
Lager Image
, which gives the characteristic curvature photon energy of inflow as
Lager Image
where Eϒ ( in ) and Eϒ ( out ) are the characteristic curvature photon energies emitted by the inflow and outflow charged particles respectively. We can also estimate the inflow radiation power as
Lager Image
where r L is the light cylinder radius and we have used
Lager Image
, where α is the inclination angle. Since the caustic effect of inflow radiation is small, this soft gammaray component should occur very near the radio pulse emitted from the polar cap. Unlike in the case of millisecond pulsars this soft gamma-ray component is very difficult to be observed in canonical pulsars because most of inflow curvature photons will be converted into pairs by the strong magnetic field (cf. Cheng & Zhang 1999, Wang et al. 2013a). Unless the inclination angle and the viewing are both small most inflow curvature photons will be reprocessed into hard X-rays and PSR1509-58 is the representative example (Wang et al. 2013a). The detail fitting of energy dependent light curves of millisecond pulsars can be found in Wang et al. (2013b).
It is believed that most of pulsar spin-down power eventually wi l l be converted from low frequency electromagnetic dipole radiation into the particle kinetic energy of the pulsar wind. However the exact conversion process from EM wave energy into particle energy is still unclear. In the Crab nebula in a distance ~3 × 10 17 from the pulsar almost all EM wave energy becomes particle kinetic energy (Kennel & Coroniti 1984ab). On the other hand, by fitting the pulsed TeV data of the Crab pulsar/nebula detected by MAGIC Aharonian et al. (2012) have argued that only a few times of the light cylinder radii from the star a good fraction of spin-down power has already been in the kinetic energy of particles. It is still controversial if the pulsed TeV gamma-rays from the direction of the Crab pulsar are emitted in the magnetosphere or outside the light cylinder. PSR B1259-63/LS 2883 is one of most studied gamma-ray binaries, which is a binary system in which a 48 ms pulsar orbits around a Be star in a high eccentric orbit with a long orbital period of about 3.4 yr. It is special for having asymmetric two-peak profiles in both the X-ray and TeV light curves (Johnston et al. 1994, 1996, 2005, Aharonian et al. 2005, 2009, Chernyakova et al. 2006, 2009, Uchiyama et al. 2009). Recently, an unexpected GeV flare has also been detected by the Fermi gamma-ray observatory several weeks after the last periastron passage (Abdo et al. 2011, Tam et al. 2011). Although X-rays and TeV gamma-rays are generally expected to be emitted from the shock region (Tavani & Arons 1997, Takata & Taam 2009), its aysmmetric light curves are very difficult to be explained. In order to explain its asymmetric two-peak multi-wavelength light curves, Kong et al. (2011, 2012) argue that one of important factors to cause such asymmetric light curve result from the fact that the particle kinetic energy of the pulsar wind is position dependent. They argue that the particle kinetic energy of the pulsar wind should gradually increase and can be approximately described by a power law as
Lager Image
where σ ( r ) is called the magnetization parameter and is defined as the ratio of the magnetic energy density and the particle kinetic energy density in the pulsar wind at position r , σ L is the magnetization parameter at the light cylinder and β is a fitting parameter of order of unity. They estimate the magnetization parameter at the light cylinder as
Lager Image
where B L is the magnetic field at the light cylinder, r L is the radius of the light cylinder,
Lager Image
is the e ± multiplicity and GJ ~ 5.26 × 10 31 ( B / 3×10 11 G)( P / 47.762 ms) -2 s -1 is the Goldreich-Julian particle flow at the light cylinder. Replacing the characteristic values of PSR B1259-63 by the values of millisecond pulsars σ L is roughly 10 4 . Most black widow systems have an orbital radius of order of 10 11 cm and a light cylinder radius of order of 10 7 cm, which imply that a good fraction of pulsar spin-down energy has already been converted into the particle kinetic energy of the pulsar wind in a distance of 10 11 cm. The TeV gamma-rays of PSR B1259- 63 are produced by inverse Compton scattering. Although the companion stars of black widow systems are white dwarfs, intense optical emission from the companion stars are observed and it is believed to result from the irradiation of pulsars. Therefore it is interesting to ask if high energy photons can be produced via inverse Compton (IC) process and such emission should exihibit orbital modulation. The characteristic energy of IC photons depends on the bulk Lorentz factor of the pulsar wind and the characteristic energy of the soft photons from the companion star.
Cheng et al. (2010) have estimated bulk Lorentz factor of the pulsar wind from millisecond pulsars γω , which is given by
Lager Image
where L 34 is the pulsar spin-down power in units of 10 34 erg/s,
Lager Image
is number ratio between e ± pairs and protons. For the characteristic values of millisecond pulsars, i.e. P ~3 ms and B ~3 × 10 8 G, which gives γω ~ 3×10 4 . Takata et al. (2010, 2012) have suggested that the companion star in the black widow systems will emit optical due to the irradiation of gamma-rays from pulsar. The optical luminosity due to reprocessing the pulsar wind/gammarays is given by
Lager Image
where Rro and Rb are the radius of the Roche lobe and the binary system respectively. The effective temperature of the companion star is given by
Lager Image
The inverse Compton scattering between the optical photons from the companion star and the pulsar wind can produce high energy photons with the characteristic energy given by
Lager Image
The luminosity of the inverse Compton scattering is given by
Lager Image
Lager Image
is the optical photon density at distance r from the companion star. If the line of sight is sufficiently closed to the companion star, i.e. r ~ Rro , the maximum inverse Compton luminosity occurs when the companion star is between the observer and the pulsar and it could be about a few percents of the spin-down power and this luminosity falls off as 1/ r . This component sensitively depends on the viewing angle. Recently the first black widow system, PSR B1957+20 has been shown some evidence of an orbital modulated gamma-ray component at energy >3 GeV (Wu et al. 2012)
We argue that the reason why millisecond pulsars can maintain a surface quadruple field with strength ~10 11 G results from interpinning of quantized flux tubes and quantized vortex lines. The magnetic flux tubes are dragged toward the spin-axis during the accretion spin-up phase. Consequently an extremely strong quadruple can be formed in the boundary between the inner crust and the core. This multiple field is very important to ensure the magnetic pair creation process can occur even in millisecond pulsars with dipolar field ~10 8 G. In the Outergap models we predict that a sub-GeV component should exist in the vicinity of radio pulse, which is emitted by the inflow current. But the characteristic energy and luminosity of this component are expected to be lower than the main outflow components. We also predict that an orbital modulated gamma-ray component could be produced by the inverse Compton scattering between the pulsar wind and the optical photons from the companion star.
We thank D. Hui, A. Kong, M. Ruderman, J. Takata, T. Tamand Y. Wang for useful discussion. This work is supported bya GRF grant of Hong Kong Government under 700911P.
Abdo AA , Ackermann M , Ajello M , Atwood WB , Axelsson M (2010) The First Fermi Large Area Telescope Catalog ofGamma-rays Pulsars ApJS 187 460 - 494
Abdo AA , Ackermann M , Ajello M , Atwood WB , Axelsson M (2011) Discovery of High-energy Gamma-ray Emissionfrom the Binary System PSR B1259-63/LS 2883 aroundPeriastron with Fermi ApJ 736 L11 - L16
Aharonian F , Akhperjanian AG , Anton G , Barres de Almeida U , Bazer-Bachi AR (2009) Very high energyγ-ray observations of the binary PSR B1259-63/SS2883around the 2007 Periastron A&A 507 389 - 396
Aharonian F , Aharonian AG , Aye KM , Bazer-Bachi AR , Beilicke M (2005) Discovery of the binary pulsar PSRB1259-63 in very-high-energy gamma rays aroundperiastron with HESS A&A 442 1 - 10
Aharonian FA , Bogovalov SV , Khangulyan D (2012) Abruptacceleration of a ‘cold’ ultrarelativistic wind from theCrab pulsar Nature 482 507 - 509
Cheng KS , Chernyshov DO , Dogiel VA , Hui CY , Kong AKH (2010) The Origin of Gamma Rays from Globular Clusters ApJ 723 1219 - 1230
Cheng KS , Ho C , Ruderman MA (1986a) Energetic radiation fromrapidly spinning pulsars. I - Outer magnetosphere gaps.II - VELA and Crab ApJ 300 500 - 539
Cheng KS , Ho C , Ruderman MA (1986b) Energetic Radiation fromRapidly Spinning Pulsars. II. VELA and Crab ApJ 300 522 - 539
Cheng KS , Ruderman M , Zhang L (2000) A Three-dimensionalOuter Magnetospheric Gap Model for Gamma-Ray Pulsars: Geometry, Pair Production, EmissionMorphologies, and Phase-resolved Spectra ApJ 537 964 - 976
Cheng KS , Zhang L (1999) Multicomponent X-Ray Emissions fromRegions near or on the Pulsar Surface ApJ 515 337 - 350
Chernyakova M , Neronov A , Aharonian F , Uchiyama Y , Takahashi T (2009) X-ray observations of PSR B1259-63 near the 2007 periastron passage MNRAS 397 2123 - 2132
Chernyakova M , Neronov A , Lutovinov A , Rodriguez J , Johnston S (2006) XMM-Newton observations of PSR B1259-63 near the 2004 periastron passage MNRAS 367 1201 - 1208
Goldreich P , Julian WH (1969) Pulsar Electrodynamics ApJ 157 869 - 880
Hirotani K (2005) High energy emission from pulsars: Outer gapscenario AdSpR 35 1085 - 1091
Hirotani K (2006) Particle Accelerator in Pulsar Magnetospheres:Super-Goldreich-Julian Current with Ion Emission fromthe Neutron Star Surface ApJ 652 1475 - 1493
Hirotani K (2008) Outer-Gap versus Slot-Gap Models forPulsar High-Energy Emissions : The Case of theCrab Pulsar ApJ 688 L25 - L28
Johnston S , Ball L , Wang N , Manchester RN (2005) Radioobservations of PSR B1259-63 through the 2004periastron passage MNRAS 358 1069 - 1075
Johnston S , Manchester RN , Lyne AG , Nicastro L , Spyromilio J (1994) Radio and Optical Observations of the PSR:B1259-63 / SS:2883 Be-Star Binary System MNRAS 268 430 - 436
Johnston S , Manchester RN , Lyne AG , D’Amico N , Bailes M (1996) Radio observations of PSR B1259-63 aroundperiastron MNRAS 279 1026 - 1036
Kennel CF , Coroniti FV (1984a) Confinement of the Crab pulsar'swind by its supernova remnant ApJ 283 694 - 709
Kennel CF , Coroniti FV (1984b) Magnetohydrodynamic modelof Crab nebula radiation ApJ 283 710 - 730
Kong SW , Cheng KS , Huang YF (2012) Modeling the MultiwavelengthLight Curves of PSR B1259-63/LS 2883. II. The Effectsof Anisotropic Pulsar Wind and Doppler Boosting ApJ 753 127 - 135
Kong SW , Yu YW , Huang YF , Cheng KS (2011) Modelling themultiwavelength light curves of PSR B1259-63/SS2883 MNRAS 416 1067 - 1076
Lyne AG , Graham-Smith F 1998 Pulsar Astronomy Cambridge University Press
Manchester RN , Taylor JH 1977 Pulsars Freeman San Francisco
Romani RW (1996) Gamma-Ray Pulsars: Radiation Processes in the Outer Magnetosphere ApJ 470 469 - 478
Takata J , Cheng KS , Taam R (2010) On the Transition fromAccretion-powered to Rotation-powered MillisecondPulsars ApJ 723 L68 - L72
Takata J , Cheng KS , Taam R (2012) X-Ray and Gamma-RayEmissions from Rotation Powered Millisecond Pulsars ApJ 745 100 -
Takata J , Shibata S , Hirotani K (2004) A pulsar outer gap model withtrans-field structure MNRAS 354 1120 - 1132
Takata J , Taam R (2009) Probing the Pulsar Wind in the γ-ray BinarySystem PSR B1259-63/SS 2883 ApJ 702 100 - 113
Takata J , Wang Y , Cheng KS (2010) Pulsar High Energy Emissionsfrom Outer Gap Accelerator Closed by a Magnetic PaircreationProcess ApJ 715 1318 - 1326
Tam PHT , Huang RHH , Takata J , Hui CY , Kong AKH (2011) Discovery of GeV γ-ray Emission from PSR B1259-63/LS 2883 ApJ 736 L10 - L15
Tavani M , Arons J (1997) Theory of High-Energy Emission fromthe Pulsar/Be Star System PSR 1259-63. I. RadiationMechanisms and Interaction Geometry ApJ 477 439 - 464
Uchiyama Y , Tanaka T , Takahashi T , Mori K , Nakazawa K (2009) Suzaku Observations of PSR B1259-63: A NewManifestation of Relativistic Pulsar Wind ApJ 698 911 - 921
Wang Y , Takata J , Cheng KS (2010) Gamma-ray Spectral Propertiesof Mature Pulsars: A Two-layer Model ApJ 720 178 - 190
Wang Y , Takata J , Cheng KS (2013a) Mechanism of the X-Ray andSoft Gamma-Ray Emissions from High Magnetic FieldPulsar: PSR B1509–58 ApJ 764 51 - 66
Wang Y , Takata J , Cheng KS (2013b) in preparation
Wu EMH , Takata J , Cheng KS , Huang RHH , Hui CY (2012) Orbital-phase-dependent γ-Ray Emissions from the Black Widow Pulsar ApJ 761 181 - 187
Zhang L , Cheng KS (1997) High-Energy Radiation from Rapidly Spinning Pulsars with Thick Outer Gaps ApJ 487 370 - 379