Arxiv:1302.0869V1 [Astro-Ph.HE] 4 Feb 2013

The Neutron Star Zoo

Alice K. Harding∗ Astrophysics Science Division NASA Goddard Space Flight Center Greenbelt, MD 20771 (Dated: February 6, 2013) Neutron stars are a very diverse population, both in their observational and their physical properties. They prefer to radiate most of their energy at X-ray and gamma-ray wavelengths. But whether their emission is powered by rotation, accretion, heat, magnetic fields or nuclear reactions, they are all different species of the same animal whose magnetic field evolution and interior composition remain a mystery. This article will broadly review the properties of inhabitants of the neutron star zoo, with emphasis on their high-energy emission.

PACS numbers: 97.60.Jd,97.60.Gb Keywords: Stars: Neutron stars, Pulsars, Binary stars

CONTENTS have been fascinating celestial objects, both for study of their exotic interiors and environments and for their im- I. Introduction 1 portant place in stellar evolution. Neutron stars (NSs) all typically have detectable pulsa- II. Rotation-Powered Pulsars 3 tions, since they are rapidly spinning and their emission patterns are highly anisotropic. This may be one of the III. Magnetars 6 few observational properties that they have in common, because they otherwise show an amazing variety of pul- IV. Compact Central Objects 9 sating and bursting behaviors. Neutron star types are classified according to the primary power source for their V. Isolated Neutron Stars 9 emission and spin evolution. Rotation-powered pulsars (RPP) derive their energy primarily from the rotation VI. Accreting Neutron Stars 10 of the NS, magnetars from magnetic field energy, iso- A. Low-Mass X-Ray Binaries 10 lated NSs (INS) from the latent heat of the NS matter, B. Intermediate- and High-Mass X-Ray and accretion-powered NSs from the energy released by Binaries 11 matter accreting onto the NS from a binary companion. C. Microquasars and Gamma-Ray Binaries 12 A subclass of accreting NSs are X-ray bursters whose bursts are powered by thermo-nuclear explosions. An ad- VII. Summary 13 ditional class, Central Compact Objects (CCO), are seen References 13 as soft X-ray point sources inside supernova remnants and seem to be quiet at all other wavelengths. Figure 1 shows where these different NS types roughly fall in pe- I. INTRODUCTION riod and surface magnetic field space. Although period is a measured quantity, surface magnetic field is derived in different ways for the various NS types. Magnetic Neutron stars are the remnants of massive stars whose fields in RPPs, magnetars, INS and CCO are derived cores collapse during the supernova explosions at the end from measured period derivatives, P˙ , assuming magnetic of their nuclear fusion lifetimes. Conservation of both the dipole radiation spin down, although in magnetars this angular momentum and the magnetic flux of the progen-

arXiv:1302.0869v1 [astro-ph.HE] 4 Feb 2013 can only give an approximate value since their spin be- itor star during the collapse gives the neutron star an ex- havior is complicated by magnetically-driven NS heating tremely high spin rate and magnetic field. The collapse and bursting. The surface dipole field at the pole deter- ends when the degeneracy pressure of neutrons balances mined assuming the observed P˙ is from magnetic dipole the gravitational forces of the matter. At that point the radiation is core radius is about 10 km and with a mass ∼ 1.4M the 1/2 center of the star has already reached nuclear densities. 3 3Ic P P˙ 12 1/2 Neutron stars thus possess the highest spin frequencies, Bd = 2 × 10 G(P P˙15) , (1) magnetic fields and densities of any known objects in the 2π2R6 Universe. Since their theoretical conception by [20] they −15 −1 where P˙15 ≡ P/˙ (10 ss ), P is in units of seconds, and I ( 1045 gcm2)andR ( 106 cm) are the neutron star moment of inertia and radius. Magnetic fields of ∗ [email protected] accretion-powered NSs cannot be measured from their 2

1016

1015 Magnetars 1014 Rotation Powered INS 13 10 Pulsars Accreting X-Ray 1012 Pulsars / HMXB

11 10 CCO

1010 MSP Surface magnetic field (Gauss) 109 LMXB

108 10-3 10-2 10-1 100 101 102 Period (s)

FIG. 1. Schematic diagram of neutron star populations with respect to their periods and derived surface magnetic ﬁeld strengths.

P˙ since their spin evolution is governed by accretion types generally occupy different areas of P − B space, torques. For accreting X-ray pulsars, the cyclotron lines although there is a good deal of overlap. The LMXB seen in their spectra give good estimates of surface field population, being progenitors of rotation-powered MSPs, strength: overlap the MSP population. CCOs lie at the lower end of the RPPs, coinciding with some of the MSP population. E c 12 The lower end of the magnetar population overlaps the Bcyc 10 G(2) 11.6keV upper end of the RPPs, and the INS overlaps both. The accreting X-ray pulsars have magnetic fields and periods E where c is the cyclotron line energy. For low-mass X- similar to some of the older RPPs. The magnetars have ray binaries (LMXB) and bursting sources, the Alfven the highest surface magnetic fields, while the rotation- radius where the NS magnetic pressure balances that of powered millisecond pulsars and LMXBs have very low the accretion flow gives an estimate of the surface mag- surface fields but the shortest periods. netic field strength: Figure 2 shows the distribution of non-accreting NS in 1/2 measured period and period derivative. The radio pul- M˙ B ∼ 12 P sars occupy the largest region of this phase space with A 10 G eq −9 −1 (3) 10 M yr their population extending from the very short period, low P˙ MSPs up to the high P˙ , high magnetic field pul- [46] where Peq (in seconds) is the equilibrium spin period sars that border the magnetar range. The magnetars of the NS and M˙ is the mass accretion rate. The NS have the highest P˙ and some of the longest periods in 3

38 36

10 erg/s 10 10

10 yr

Radio pulsar

Period Derivative Magnetar

X, -ray pulsar 11

10 yr

INS

RRAT

30 32

10 erg/s 10 CCO

Period (s)

FIG. 2. Plot of period vs. period derivative for the presently known rotation-powered pulsars, Isolated Neu- tron Stars (INS), Compact Central Objects (CCO), Rotating Radio Transients (RRATs) and magnetars (from http://www.atnf.csiro.au/people/pulsar/psrcat/). Lines of constant characteristic age, P/2P˙ , and dipole spin-down luminosity, E˙ d, are superposed. the NS zoo. with some having P ∼ 11 s (the longest pe- emission are known as rotation-powered pulsars (RPP). riod of a radio pulsar is about 8 s). The INS have periods The energy from their spin down appears as broad- very similar to those of magnetars, also reaching up to 11 band pulsations from radio to gamma-ray wavelengths s, but with P˙ and magnetic fields about a factor of ten and as a wind of energetic particles flowing into their lower than those of the lower field magnetars. The CCOs surrounding pulsar wind nebulae. Since the discov- have very low P˙ , almost as low as the MSPs, but their ery of RPP through their radio pulsations in 1967 [62], spin periods are more like those of young RPPs. High en- more than 2000 radio pulsars are now known with pe- ergy pulsars (RPPs with X-ray or gamma-ray pulsations) riods ranging from a few ms to several seconds [82], typically have the highest spin-down power (see Eq [4]) http://www.atnf.csiro.au/people/pulsar/psrcat/). X- but are not necessarily the youngest, as many of the very ray, gamma-ray and optical pulsations were soon discov- old MSPs are efficient X-ray and gamma-ray pulsars. ered in a few of these pulsars by folding the time series obtained at these wavelengths at the radio periods. At present, there are over 100 RPP detected at X-ray ener- II. ROTATION-POWERED PULSARS gies and over 130 gamma-ray pulsars [11]; most were discovered using known radio ephemerides, but many were also discovered through their X-ray or gamma-ray pulsa- Neutron stars that are spinning down as a result tions and are radio quiet. The spin down of RPPs is typi- of torques from magnetic dipole radiation and particle 4

FIG. 3. Skymap in Galactic coordinates showing the gamma-ray pulsars of different types detected by Fermi: Blue squares: radio quiet pulsars. Red triangles: millisecond gamma-ray pulsars. Green circles: radio loud gamma-ray pulsars. Black dots: Gamma-rays were phase-folded using a rotation ephemeris. Gray dots: Pulsars for which no rotation ephemeris was available. From [11]. cally smooth and predictable, but they have occasionally RPPs gives estimates of their main electrodynamic prop- been observed to undergo sudden changes in spin called erties. The Poynting flux of a rotating dipole in vacuum ‘glitches’, where the period decreases and then recovers gives an estimate of the spin-down power: back to its normal spin-down rate on a timescale of days 2 2 4 4π IP˙ 2B0 Ω 31 −1 2 −4 to weeks. E˙ d = = =10 erg s B P (4) P 3 3R6c5 12 There are two main populations of RPPs: normal pulsars having characteristic ages τ = P/2P<˙ 100 Myr, and where Ω = 2π/P is the spin angular velocity, R the τ > NS radius, B0 is the surface magnetic field and B12 ≡ millisecond pulsars (MSP) with ∼ 100 Myr. The peri- 12 ods of the normal pulsars are thought to have increased B0/10 G. The rotating dipole induces strong electric at a steady rate from their birth periods at the time of su- fields, E = Ω × r × B that in vacuum are orders of mag- pernovacorecollapse.ThebirthperiodsofRPPsarenot nitude larger than the gravitational force, pulling charges well known, and inferences of the period range from tens out of the NS surface to fill the magnetosphere [47]. The of ms to hundreds of seconds [74][49]. MSPs are thought maximum charge density, which screens the electric field to have originally been members of the normal RPP pop- along magnetic field lines, is ulation, spun down for tens of Myr and then spun up by Ω · B ρGJ − (5) accretion from a binary companion [14]. MSPs make 2πc up about 10% of the RPP population and about 80% of them are in binary systems. A large number of new which gives a scale for the current density flowing along J ρ c MSPs have recently been found through radio searches at the open field lines, GJ = GJ . The particles and mag- positions of Fermi unidentified gamma-ray sources (e.g. netic field corotate with the NS out to the light cylinder R c/ [98]), doubling the known numbers of radio MSPs in the radius, LC = Ω, where the corotation velocity reaches Galactic disk. Their radio ephemerides can then be used the speed of light. The magnetosphere is divided into re- R to find the gamma-ray pulsations. They are extremely gions of open field lines that cross LC and regions of R good clocks since they rarely glitch, they have very little field lines that close within LC. The footpoints on the of the timing noise seen in young RPPs, and their spin NS surface of the open/closed boundary define the polar down is very stable. These characteristics make them po- caps. The induced potential drop in vacuum across the tentially valuable sources for use in celestial navigation open field lines is [31] and gravitational wave detection [68]. 2 1 2π 3 12 −2 Vpc B0R =6× 10 Volts B12 P (6) The rotating dipole model that seems to well describe 2 cP 5

FIG. 4. Gamma-ray (black) and radio (red) pulse proﬁles for a selection of pulsars detected by Fermi, from [11].

The open field lines and the particles that flow along high-energy exponential cutoffs in the range 1 - 10 GeV. them become the pulsar wind, which carries the bulk of The pulsed gamma-ray emission of RPPs results from the Poynting flux beyond RLC. radiation of particles accelerated to energies of order 10 The pulsed radiation of RPPs carries less than about TeV by electric fields parallel to the magnetic fields in 10% of the total spin-down power. Most of the power the open magnetosphere. Acceleration could take place in pulsed emission is in gamma rays around a GeV for near the polar caps [17], in the outer magnetosphere near all but a few RPPs, notably the Crab pulsar whose the light cylinder [33] in slot gaps along the last open pulsed emission power peaks in hard X rays. The ra- field lines [18, 89] or from reconnection in the striped dio pulsations make up typically only 10−4 or less of wind outside the light cylinder [92]. The Fermi measure- the spin-down power, but are easier to detect since the ment of an exponential shape for the spectral cutoff of photon flux is much higher. Since the gamma-ray pul- the Vela pulsar [1] ruled out emission in the intense mag- sations typically carry the highest percentage of E˙ d, netic field near the polar caps, whose spectra would have they can reveal the most information about the parti- much sharper, super-exponential cutoffs [37], as being cle acceleration in pulsar magnetospheres. The discov- the primary source of the observed gamma-rays. Instead, ery of gamma-ray pulsations from over 130 RPPs by the thesoftercutoffsaswellastheobservedprofileshapes Fermi Gamma-Ray Space Telescope [7, 11], where there and peak separation vs. radio phase-lag correlation [7] were only seven known previously [110], has thus revo- predicted by outer magnetosphere models [102], strongly lutionized the study of this type of NS (see Figure 3). suggest that the gamma-ray emission comes from high Two new populations of gamma-ray pulsars have been altitudes up to and possibly beyond the light cylinder. established with Fermi: radio quiet pulsars discovered In these models, the narrow gamma-ray peaks are caus- through blind searches for gamma-ray pulsations [2, 93], tics formed by phase bunching of the emission along the and gamma-ray MSPs whose pulsations are found us- open field boundary through aberration and light travel ing radio ephemerides [3]. These are presently equal time effects [42, 52]. The main radiation mechanism at in number to the radio-loud gamma-ray pulsars. The GeV energies is thought to be curvature radiation by the gamma-ray pulse profiles show recurring patterns of two primary particles that have been accelerated to 10 TeV narrow peaks separated by phase intervals of 0.1 - 0.5, energies, but inverse-Compton emission may also con- with occasional single peaks. The most striking pattern tribute radiation at the higher energies (up to 400 GeV) is that the gamma-ray peaks are not in phase with the seen from the Crab pulsar [13]. radio pulses, but typically arrive later in phase, with a Many RPPs (over 100 at present) also show X-ray lag that is correlated with the gamma-ray peak separa- emission, with pulsations detected in many of these. This tion (smaller lags for larger peak separation, see Figure emission is usually made up of two components: non- 4). The pulsed gamma-ray spectra are power laws with thermal emission that is probably magnetospheric, and 6

0.05 - 0.1 keV thermal emission from surface cooling or nova shell that is moving more slowly than the relativistic heated polar caps (see [74] for a detailed review). Some, pulsar wind, or by the bow shock driven by the supersonic like the Crab pulsar, show only non-thermal X-ray emis- motion of the pulsar through the interstellar medium. sion, but in these cases the thermal components are prob- The particles are mostly electron-positron pairs from the ably buried under the strong non-thermal emission. Oth- pair cascades [36] that occur in the magnetosphere, but ers, like MSP J2124-3358, show only thermal components there could be some protons or positively-charged ions with statistically insignificant non-thermal emission. The as well. Most of the Poynting flux of the RPP must at thermal and non-thermal peaks in the pulse profiles are some point, either in the wind or at the wind termination usually not in phase, but the non-thermal peaks are shock, be transferred to the particles. In the Crab neb- sometimes in phase with one or two gamma-ray peaks. ula, the highest energy particles receive nearly the entire In the youngest RPPs, emission from cooling of the NS open-field voltage, as deduced from the ∼ 100 MeV cut- surface is thought to dominate the thermal radiation. off in the synchrotron spectrum [19, 39]. Many PWNe In contrast, middle-aged RPPs like Vela, Geminga and are detected at energies up to several or even tens of PSR J0659+1414 show a two-component thermal spec- TeV [38]. This VHE emission is thought to be inverse- trum from both heating and cooling, and a non-thermal Compton emission from the nebular particles, scattering power-law component. MSPs show all of this behavior either their own synchrotron radiation as in the Crab but, since these sources are too old to have detectable PWN or, in most cases, the microwave and infrared back- emission from cooling, the thermal components, some of ground radiation. Powerful gamma-ray flares from the which are multiple blackbodies, must be due to polar cap Crab nebula that appear only in the synchrotron emis- heating. The non-thermal emission is in most models due sion component have recently been observed by Fermi to synchrotron radiation from electron-positron pairs, ei- [8] and AGILE [109]. Their short timescales of hours ther from the polar cap or the outer gap, emitting at to days indicate that they are coming from a very small high altitude [56, 63]. The non-cooling thermal emission region (< 1 arcsec) and their energies of up to several could come from polar caps heated by high-energy par- GeV challenge traditional models of acceleration in pul- ticles flowing toward the NS from the polar cap [55] or sar wind nebulae [28]. outer gap [50] accelerators. There is an unusual and growing subpopulation of RPPs known as Rotating Radio Transients (RRATS) III. MAGNETARS (see [77], for review). They were discovered only very recently [84], through detections of their single, isolated NSs whose primary power source is the tapping of en- radio pulses. They show a variety of transient radio be- ergy stored in their magnetic fields are known as magne- havior, ranging from nulling (or turning off for up to 4 tars (see [119], for review). There are two sub-classes of ∼ 10 s) for long time periods to steady pulsations that magnetars, Anomalous X-Ray Pulsars (AXPs) and Soft are highly modulated. While nulling behavior in normal Gamma-Ray Repeaters (SGRs), that were thought for radio pulsars has been known for many years, the nulling many years to be separate and unrelated objects. Today, of RRATS is extreme. Over 70 RRATS have been iden- we know that SGRs and AXPs are both NSs possessing ˙ 14 16 tified and P s have been measured for about 20 of these magnetic fields of unprecedented strength of 10 − 10 ˙ sources. Their P and P are scattered throughout the nor- G, and that show both steady X-ray pulsations as well as mal RPP population (see Figure 2), but they tend to be soft γ-ray bursts. Their inferred steady X-ray luminosi- > 5 older (τ ∼ 10 yr) and a number of them have higher than ties are about one hundred times higher than their spin- average surface magnetic fields. X-ray pulsations have down luminosities, requiring a source of power way be- been detected from one RRAT, PSR J1819-1458 [85], re- yond the magnetic dipole spin-down that powers RPPs. vealing a thermal spectrum and an absorption line at 0.5 New high-energy components discovered in the spectra of keV, possibly due to proton cyclotron resonant scattering 14 a number of AXPs and SGRs require non-thermal par- inafieldof2×10 G. The causes of such radio transient ticle acceleration and look very similar to high-energy behavior is unknown, but global changes in the currents spectral components of young rotation-powered pulsars or charge density in the pulsar magnetosphere has been [41] (see Figure 6). suggested [70, 81]. It is not yet clear how RRATS fit SGRs were first detected around 1979 as γ-ray tran- into the normal RPP population. They share some sim- sients and were thought to be a type of classical γ-ray ilar properties with INS (see Section V), but if they are burst. They undergo repeated bursts with several tenths a separate evolutionary group they would significantly of second duration and average energy 1040 − 1041 erg, raise the birthrate of NSs in the Galaxy [76]. and their bursting often occurs in episodes spaced years In addition to pulsed emission, broadband, un-pulsed apart. They more rarely undergo giant superflares of to- emission from radio to high-energy (TeV) gamma-rays is tal energy 1045 − 1047 erg, consisting of an initial spike of detected from the pulsar wind nebulae (PWNe) associ- duration several hundred ms followed by a longer decay ated with RPPs. The PWN is a repository of high-energy of duration several hundred seconds showing pulsations particles and fields from the RPP, accumulated over the (Figure 5). Such superflares have been observed in three pulsar’s lifetime and trapped by the surrounding super- SGR sources, SGR0526-66 (the famous 5th March 1979 7

FIG. 5. Light curve of the 2005 superburst from SGR 1806-20, seen by the Burst Alert Telescope on Swift. Strong pulsations are seen in the decaying tail of the burst. (from [90]) event) [34], SGR1900+14 [69] and SGR1806-20 [90]. In of the same type of object, very strongly magnetized, 1998, SGR1806-20 was discovered to have 7.47 s pulsa- isolated NSs possibly powered by magnetic field decay. tions in its quiescent X-ray emission [78] and a large P˙ In both sources, the high-state quiescent luminosities of 15 35 −1 that implies a huge surface magnetic field of 10 G LX ∼ 10 erg s are much higher than their spin-down 33 −1 if due to dipole spin-down. Quiescent periodicities of 8 luminosities of E˙ d ∼ 2 − 6 × 10 erg s , demanding an s and 5.16 s and large P˙ were subsequently detected in alternative power source. SGR0526-66 and SGR1900+14, implying similarly high The quiescent spectra of AXPs and SGRs (see Figure surface magnetic fields. In all three sources, the qui- 6) consist of a thermal component fit by ∼ 0.5-1 keV escent periods are the same periods seen in the decay blackbodies and one or more non-thermal components. phases of their superflares. The quiescent pulse profiles Until recently, non-thermal spectra were seen only be- are very broad and undergo dramatic changes just below 10 keV and were fit with steep power laws having fore and after superflares. The profiles are often more photon indices ∼ 3 − 4. When INTEGRAL and RXTE complex, with multiple peaks before flares, changing to recently measured the spectra above 10 keV for the first more simple single peaked profiles following the flares. time hard, non-thermal components were discovered in Since the modulation is thought to result from beam- three AXPs, and also SGR 1806-20. In two of the AXPs, ing along magnetic field lines, the profile changes signal the differential spectra between 10 keV and 50 keV are a re-arrangement of the magnetic field structure during extremely flat: 1E 1841-045 [79] has a power-law pho- the flares. All of the SGRs lie near the Galactic plane ton index of s =0.94 and 4U 0142+61 [40] displays a and are thought to have distances around 10-15 kpc (ex- photon index of s =0.45, both much flatter than the cept for SGR0526-66, which is in the LMC). Recently, a steep non-thermal components in the classic X-ray band. transient SGR was discovered with an apparent surface RXS J1708-40 possesses a slightly steeper continuum 12 magnetic field strength < 7.5×10 G [99]. This could be with s =1.18. The non-thermal tail of quiescent emis- an aging magnetar that has experienced significant field sion in SGR 1806-20 is similarly pronounced [86, 88], but decay over its lifetime. somewhat steeper, with a photon index of s =1.6 − 1.9 The first AXPs were discovered as pulsating X-ray extending to 100 keV. Such hard non-thermal compo- sources in the early 1980s by EXOSAT and were thought nents require continuous particle acceleration during the to be a strange (anomalous) type of accreting X-ray pul- quiescent state. Transient, highly variable radio emis- sar. They are bright X-ray sources possessing luminosi- sion has been detected from several magnetars [29, 30] 35 −1 ties (in their highest states) of LX ∼ 10 erg s , but and may be correlated with their X-ray luminosity and show no sign of any companion or accretion disks that outbursts[100]. would be required to support the accretion hypothesis. It was proposed early on [75, 97] that SGRs were NSs TheAXPshavepulsationperiodsinarelativelynarrow with strong magnetic fields in the range 1011 − 1013 Gto range of 5 - 11 s and are observed to be spinning down confine the super-Eddington burst radiation. That SGRs with large period derivatives [115]. Their pulse profiles and AXPs had much higher magnetic fields above 1014 are broad and very similar to those of SGR sources. The G was proposed by Duncan & Thompson [DT] [43] be- very high surface magnetic fields of 1014 − 1015 Gim- fore their existence was observationally verified. In this plied by dipole spin-down were originally controversial, model, some NSs generate huge magnetic fields by dy- but have come to be accepted after the quiescent periods namo action soon after their birth in supernova explo- were found in SGRs and especially following the recent sions. Such high fields have several properties that can discovery of SGR-like bursts from several AXPs [73]. It cause the NSs to behave differently from NS with lower is now believed that SGRs and AXPs are two varieties fields. These fields can decay on much shorter timescales, 8

FIG. 6. Magnetar spectral components. Left: quiescent phase-averaged spectrum of AXP 4U 0142+61 showing lower energy thermal and power law components below 10 keV and high-energy component above 10 keV (from [41]). Right: typical burst spectra from AXP/SGR, compared with burst spectra from X-ray and gamma-ray bursts (from [119]). due to ambipolar diffusion [60]. According to the DT model [113], the magnetar su- −2 perflares result from reconnection of sheared or twisted 5 Bcore external field lines, leading to particle acceleration and t =∼ 10 yr (7) amb 1015G radiation of hard emission. The estimated luminosity of such events, Diffusion of magnetic flux out of the NS core on these 2 2 timescales provides the power to magnetars in the B Bcore core R3 ≈ 4 × 1046erg , (8) DT model. Magnetar-strength fields also apply higher 8π 1015G stresses to the stellar crust, so that the yield strain can exceed the crustal strength. This property is responsi- is similar to observed luminosities of superflares. The ble for the small SGR and AXP bursts in the DT model smaller bursts result from cracking of the crust, which 15 [112]. If a toroidal component of the field Bcore > 10 is continually overstressed by diffusion of magnetic flux G develops in the interior of the star, it can twist the ex- from the NS interior. The shaking of magnetic footpoints ternal field [45]. Such action can cause the superflares if then excites Alfven waves that accelerate particles. The the twisted field lines reconnect. Finally, due to the much energy radiated in such events would be faster heat transport in very strong magnetic fields, there −2 2 2 is a greater heat flux through the crust [60]. Such a prop- ∼ 41 B0 l θmax ESGR = 10 erg , (9) erty may explain the much hotter surface temperatures 1015G 1km 10−3 of magnetars and the high quiescent X-ray emission. l B Magnetar fields also produce a variety of different be- where is the length scale of the displacements, 0 is θ havior of radiative processes (see [44] and [53] for re- the crustal field and max is the yield strain of the crust. views). In general, in magnetic fields approaching and The quiescent emission in the DT model is powered by 13 exceeding the quantum critical field Bcr =4.4 × 10 magnetic field decay through conduction of heat from the G, for radiative processes such as Compton scatter- core. The NS crust is heated to a temperature of ing, cyclotron and synchrotron emission and absorp- T 5/9 tion, and pair production and annihilation, classical de- T ∼ . × 6K core crust = 1 3 10 8 (10) scriptions are largely inaccurate and QED descriptions 10 K must be used. In addition, new processes become possi- where Tcore is the core temperature, and luminosity ble in strong fields, such as one-photon pair production 4.4 and annihilation, vacuum polarization and photon split- ∼ 35 −1 Bcore Lx = 6 × 10 ergs . (11) ting, that cannot take place in field-free environments. 1016G These processes, in particular vacuum polarization [80], strongly influence the propagation of photons in the NS [111] have proposed that the hard, non-thermal qui- atmosphere and the spectrum of the emergent radiation. escent component is due to the creation of a strong 9

E induced by twisting of field in the closed magneto- predict sudden cooling at around the age of CasA (300 sphere (see also [26]), producing synchrotron radiation yr). These measurements could thus also be the strongest from electron acceleration at high altitude. An alterna- evidence for superfluidity in NS cores. Some CCOs show tive model for magnetar activity and emission has been absorption lines in their spectra. 1E 1207.4 -5209 has 2-4 discussed by [61]. The burst and quiescent radiation are harmonically spaced lines at 0.7, 0.14, 0.21 and 0.28 keV a result of shocks from fast-mode plasma waves and the [27] and the spectrum of J0822.0-4300 may have an ab- hard quiescent component is due to a pair-synchrotron sorption feature at 0.8 keV [48]. These could be electron cascade. On the other hand, Baring [21] and Baring & cyclotron lines in a field around 8 × 1010 G. Harding [22] propose that resonant Compton upscatter- Several mysteries surround these objects. Were they ing of thermal X-rays by accelerated particles produces really born with such low magnetic fields or were higher the quiescent hard emission. interior fields submerged by mass fallback accretion to emerge later, enabling them to turn on as normal RPPs [64]? This hypothesis would fit with the need for C at- IV. COMPACT CENTRAL OBJECTS mospheres that make their spectra consistent with full- surface NS cooling. Compact Central Objects (CCOs) are X-ray sources detected close to the centers of young supernova remnants (SNRs) that have no apparent emission in other V. ISOLATED NEUTRON STARS wavebands and no binary companions. Although these sources have been known and studied for several decades Another class of NSs that appear to be thermally cool- without much understanding of their nature, exciting re- ing with no emission outside the soft X-ray band, except sults over the past few years have brought them into for faint optical/UV counterparts, are known as Isolated the forefront of NS studies. They have soft, exclusively Neutron Stars (INS) (see review by [71]). Although these thermal spectra in the few hundred eV range and X- properties are similar to those of CCOs, they are a dis- 33 34 −1 ray luminosities around 10 − 10 erg s .Aboutten tinct class because they lack any observable associated CCOs are presently known, including the central sources supernova remnant or nebula. There are presently seven of CasA (Figure 7), Puppis A and Kes 79 supernova confirmed INS (sometimes referred to as The Magnificent remnants. Several, J1852+0040 in Kes79 [49], J0822.0- Seven), six of which have measured weakly modulated X- 4300 in Puppis A [48] and 1E 1207.4 -5209 in PKS 1209- ray pulsations with periods between 3 s and 11 s, much 51/52 [121] have detected pulsations in the hundreds of longer than those of CCOs. Two have high measured P˙ ms range. J1852+0040 has a detected P˙ [49], indicat- near 10−13, giving surface magnetic fields around 1013 G, ing that it is spinning down like a RPP. The measured characteristic ages τ ∼ 2 Myr and spin-down luminosity P and either measurements or constraints on P˙ indicate around 4 × 1030 erg s−1 assuming dipole spindown. They that these sources have very low magnetic fields in the have significant proper motions, indicating distances of range 1010 − 1011 G assuming magnetic dipole braking, less than 500 pc. The INS thus appear to be very nearby, so that their birth periods were close to their present val- cooling middle-aged NSs. ues. Since their SNRs are all young, ∼ 103 − 104 yr, they The X-ray spectra of INS are soft blackbodies with were probably born with unusually low magnetic fields, kT in the range 40 - 100 eV, and X-ray luminosities in which makes them ‘anti-magnetars’ [49]. the range 1030 − 1032 erg s−1, with no measurable power The observed spectra are purely thermal, but some law or non-thermal components. Four have observed have multiple blackbody components. The measured optical and/or UV counterparts that appear at first to radii of the emitting areas, assuming blackbody spectra be Rayleigh-Jeans extensions of the X-ray blackbodies. or even H atmospheres, are less that 1 km, much smaller However the optical spectra are about a factor of 10 than a NS radius. This is surprising if the emission is above the extrapolation of the X-ray blackbody. Alterna- from cooling of the NS, especially since the very low mag- tivestosimpleblackbodyfitshavebeenexplored,suchas netic fields do not preferentially conduct the cooling to a thin magnetic, ionized H atmosphere for the spectrum small hotspots at the poles [120]. Recently though, fits of J1856.5–3754, the original INS [65], or high B-field at- of the spectrum of the CCO in CasA with C atmosphere mospheres [96]. Most INS show multiple or complex ab- models, that are harder than those of H or He, give emit- sorption lines in their spectra at energy between 0.3 and 1 ting radii that are consistent with a NS radius of 10 - keV (e.g. [72]). It is not clear what causes these features, 12 km [66]. This could be an indication of recent ac- but some possibilities are proton cyclotron, and neutral cretion activity that would deposit some heavy elements or molecular H absorption. The X-ray luminosities of on the NS surface. Even more recently, the temperature several INS are larger than their spin down luminosities, of CasA is observed to be decreasing, by 4% in just 10 by a factor of 60 in the case of RX J0720, so that rotation years [59]. If this is in fact the first measured cooling cannot power the observed X-ray emission. Furthermore, of a NS, the rate of cooling is much more rapid than thetemperaturesandthermalluminositiesofINSaretoo standard NS cooling models predict. The measurements high (by about a factor of 10) for conventional cooling require proton or neutron superfluid models [105] that at their characteristic ages. This could imply that their 10

FIG. 7. X-ray image of the supernova remnant Cassiopeia A taken with the Chandra X-Ray Observatory, color coded for energy (Red: 0.5-1.5 keV; Green: 1.5-2.5 keV; Blue 4.0-6.0 keV). The CCO is the central white point source. Credit: NASA/CXC/MIT/UMass Amherst/M.D.Stage et al. magnetic fields have not been constant since birth but VI. ACCRETING NEUTRON STARS have decayed [60], providing another power source in addition to cooling and making them younger than their NSs in binary systems can accrete matter from the spin-down ages. companion stars, either from the stellar winds or from an It is not clear how INS are connected to the rest of accretion disk that forms if the companion overflows its the RPP population, other than being NSs with thermal Roche Lobe. The gravitational energy from the infalling emission. Do they lack detectable radio emission because matter provides at least part of the energy for the ob- our viewing angle does not cross the radio beam? A fit served radiation and the accretion torques dominate the of the X-ray profile of J1856.5–3754 [65] constrains mag- spin evolution. Despite these common properties, accret- netic inclination and viewing angles to be far apart. Or ing NSs display a wide variety of behaviors, depending on are they dead radio pulsars or extreme RRATS? They the NS magnetic field strength, mass of the companion do lie very near or beyond the empirical radio death and properties of the accretion. line. But why do they not have detectable non-thermal high energy emission? Other RPPs of similar ages or older have detectable X-ray or gamma-ray pulsations. If their magnetic fields are decaying or have decayed from A. Low-Mass X-Ray Binaries magnetar-strength fields, then they could be old magnetars [60] that wound up in the normal pulsar part of P -P˙ LMXBs are binary systems in which one member is space [95]. At the moment the population of INS is too a NS or black hole and the other star is a low-mass small to draw any firm conclusions about their nature main sequence star, white dwarf or red giant that fills its and origin, but searches for new INS are underway. Roche Lobe, transferring matter onto the compact object 11 through an accretion disk. Almost all of the radiation is surface. The true rotation frequency of the NS is then the emitted as X rays with a very small amount (about 1%) observed asymptotic frequency. The NS spin frequency in optical light. The NS in LMXB are thought to be can also appear as the frequency difference between pairs old, have weak magnetic fields (∼ 108 G) and are being of QPO that move up and down in frequency together. spun up by torques from the accretion disk. One expla- Since the difference in the frequencies of the QPO pair nation for their very weak fields is that the accretion has stays the same and is in the range 200 - 400 Hz, it is reduced or submerged the magnetic flux [35, 101, 103]. thought to be the NS rotation frequency. In some cases Depending on the accretion rate M˙ , they will reach an though, the frequency difference appears to be closer to equilibrium where the pressure of the NS magnetic field 1/2 of the spin frequency, and the relation between QPO balances the pressure of the accretion flow after about and spin frequency has been questioned [83]. 1 Myr, at which point the NS period will have reached The distribution of NS spin frequencies from these an equilibrium period Peq.ThevalueofPeq depends on various measurements cuts off suddenly at a maximum theNSfieldandM˙ given by Eq (3). For magnetic fields frequency around 620 Hz [32], which is well below the 8 ∼ around 10 GandM˙ near the Eddington limit, Peq will 1000 Hz break-up frequency of a NS. This indicates reach millisecond periods. X-ray pulsations at millisec- that there are some mechanisms that halt or counteract ond periods were in fact found in a number of LMXBs the spin up process. Several possibilities have been sug- with the Rossi X-Ray Timing Explorer (RXTE), the first gested, including gravitational wave damping [51, 116], being SAX J1808.4-3658 [118]. Around 24 of the more making these good sources of kHz gravitational radiation than 100 LMXB sources have shown ms X-ray pulsa- for possible detection with Advanced LIGO. tions, with frequencies in the range 100 - 700 Hz (see review by [32]). Discovery of these long-sought ms X- ray pulsations in LMXB finally established them as the B. Intermediate- and High-Mass X-Ray Binaries progenitors of the rotation-powered MSPs. The transi- tion between an LMXB and a MSP was actually observed Binary systems consisting of a neutron star or black > when ms radio pulsations were discovered from the source hole and high-mass O or B star (M ∼ 5M )areknown J102347.6+003841 which less than ten years earlier had as High-Mass X-Ray Binaries (HMXB), an example be- shown optical emission indicating an accretion disk [16]. ingVelaX-1.InalargesubclassofHMXBs,theNS Some LMXBs have steady X-ray emission, with occa- orbits a Be star in a very eccentric orbit, accreting mate- sional bursting behavior while others are detectable only rial only occasionally in outbursts when the NS crosses a during their bursts. A good fraction are transient, with disk of material surrounding the Be star. An example of outbursts on month to year timescales. The outbursts a Be X-ray binary is A0535+26. If the donor star is an < < typically last weeks to months. The X-ray bursts are intermediate-mass star (1.0M ∼ M ∼ 5M), the sys- thought to be thermonuclear explosions on the NS sur- tem is an intermediate-mass X-ray binary (IMXB) [94], face, (see [108] for review) when the density and temper- a well-known example being Her X-1. In HMXBs, the ature of the accreted material reach the critical point compact object accretes material from the wind of the for igniting nuclear reactions. The bursts have typi- companion star or from a Be star disk, while in IMXBs cal timescales of around 10 - 100 s and total energy of the NS accretes from a disk as in LMXBs. Both of these 39 40 −1 10 − 10 erg s . LMXBs also show kHz quasiperiodic classes, numbering about several hundred, are bright X- oscillations (QPO) in their X-ray emission (see [114] for ray sources, and a large number show either persistent review), thought to originate in the inner accretion disk or transient X-ray pulsations (in the case of Be star bi- flow. In most models of QPO, the oscillation frequen- naries). In contrast to the ms pulsations of LMXBs, pul- cies are the orbital frequencies of accreting matter. Since sation periods of HMXBs are much longer, in the range these are stable orbits around a NS, they must be outside 1 - 1000 s. The relation in Eq (3) then requires that the the innermost stable circular orbit (ISCO) determined by NS magnetic field in these sources are near 1012 G. In General Relativity, so their frequencies can be used to fact, many X-ray pulsars in IMXBs and HMXBs (about place limits on the NS mass [87]. Since the maximum 20) show absorption features in their spectra [58] that measured QPO frequencies are about 1200 Hz, using the are almost certainly due to electron cyclotron resonance 2 ISCO of 6GM/c in the Schwarzschild metric gives the scattering. The energies of these cyclotron resonant scat- NS mass M<2M [122]. tering features (CRSF) are in the range 10 - 50 keV, also A subset (∼ 13) of LMXB show ms pulsations in their suggesting (Eq. [2]) that the surface magnetic of the NSs persistent emission. The rest (11) show the pulsations are around 1012 − 1013 G. only during X-ray bursts [106, 117] or in QPO behavior. Because the NS magnetic fields are high in these sys- During the X-ray bursts, oscillations have been observed, tems, the Alfven radius usually during the decay phase (see Fig 8), whose fre- −2/7 quency rises to an asymptotic limit as the intensity drops M˙ [107]. The interpretation of this behavior is that the ro- r × 8 B4/7 M/M −1/7 , A =3 10 cm 12 17 −1 ( ) tation of the NS is modulating the burst hot spot which 10 gs expands and rises in altitude, then drops back onto the (12) 12

FIG. 8. X-ray burst oscillations in 4U 1702-43 (from [107]). The histogram shows the X-ray intensity of the X-ray burst. The contours show the Fourier power level as a function of frequency and time, indicating a drifting oscillation starting at 328 Hz (t = 7 s) and ending at 330.5 Hz. lies at a much larger distance from the NS than in sorption lines and as many as four harmonics also appear LMXBs, and the accreting material is strongly chan- in some X-ray pulsars such as 4U 0115+63 [57]. They are neled via the magnetic field lines onto a small polar formed through resonant scattering of continuum pho- cap. The infalling material then heats the NS atmo- tons by electrons that have a thermal distribution of mo- sphere above the polar cap to temperatures around 108 menta along the strong magnetic field lines (see [53], for K . The proposed mechanisms by which the kinetic en- a detailed description of the physics of CRSFs). Since ergy of the infalling material is decelerated and trans- the resonant scattering cross section has a strong depen- ferred to heat depends on the mass accretion rate, which dence on angle, CRSFs can provide a good diagnostic in turn determines its pressure and the source luminos- of the geometry and physical conditions of the emitting ity LX = GMM˙ . If the accretion pressure exceeds the atmosphere [15, 104]. radiation pressure, which is much higher than the spher- ical Eddington limit due to the channeling of the flow, a radiative shock will form at some distance above the C. Microquasars and Gamma-Ray Binaries NS surface [23]. If the pressure of the accretion flow does not exceed the radiation pressure, the atmosphere is heated either by Coulomb collisions [91] or a collision- A subclass of accreting X-ray binaries are called micro- less shock. In the case where a stand-off shock forms, a quasars since they display properties similar to those of slowly sinking accretion column radiates mostly perpen- quasars, including rapid variability of their X-ray emis- dicular to the magnetic axis in a fan beam [24], while in sion and radio jets. Strong, broadband emission as well the Coulomb-heated case, the atmosphere is a thin slab as broad emission lines are also observed from their ac- at the NS surface and radiates primarily along the mag- cretion disks. They consist of a compact object, either netic axis [54]. The continuum spectrum of the X-ray neutron star or black hole, with a normal companion. pulsar radiation is thermal bremsstrahling and Compton These sources undergo repeated and sometimes periodic scattering emission (e.g. [25]). The CRSFs appear as ab- radio, optical and X-ray flaring with associated formation of relativistic jets. Cygnus X-3 is a well-known micro- 13 quasar that could be a black hole or neutron star in a 4.8 existence requires some kind of early dynamo activity. hr orbit around a massive hot star. One of the strongest The intermediate-strength magnetic fields of RPPs de- X-ray sources in the sky, it is also detected from radio to termines their spin-down behavior and radiation prop- GeV gamma-ray wavelengths [4]. erties. However, the overlap of the magnetar and RPP The term gamma-ray binary has come to include bi- populations is not understood at present. Why don’t nary systems that contain a compact object and emit high-field RPPs burst or flare? Why don’t lower-field high energy gamma rays. Since several microquasars magnetars sometimes behave more like RPPs? INS may have been detected in gamma rays, they are one of the be old magnetars whose high fields have decayed, but gamma-ray binary source types. The other gamma-ray they do not account for the existence of younger RPPs binary source type is a system containing a RPP orbit- with magnetar-strength fields. Not only the strength of ing a massive star. The prototype of this source class the magnetic field but also its configuration may be im- is PSR B1259-63, which is a 48 ms radio pulsar in a portant in making a NS a magnetar or a RPP. Magnetic 3.4 yr eccentric orbit around a Be star. It appears as a field decay is a critical link between other NS popula- normal radio pulsar at apastron, but at periastron the tions as well. “Decay” of the magnetic field is necessary pulsar is thought to make several crossings through the for normal RPPs to evolve into MSPs through accretion Be star disk, producing X-rays and gamma-rays at GeV andspinupinLMXBs.Somekindofaccretion-driven [9] and TeV [12] energies. The X-ray binaries LSI+61 303 field reduction is the most likely mechanism, but it is and LS5039 have been detected by Fermi in high-energy controversial since it is not clear how effective it is or on gamma rays [5, 6], have characteristics similar to PSR what timescale a buried field might re-emerge. One piece B1259-63 and are thought to contain undetected pulsars. of evidence in favor of accretion-driven field reduction is The Fermi source 1FGL J1018.65856 is the first Galactic the fact that NSs in LMXBs, which are older systems 8 binary system to have been discovered at gamma-ray en- (> 10 yr), have mostly low fields and NSs in HMXBs, 7 8 ergies [10] but appears similar to LS5039, as both X-ray which are younger systems (10 − 10 yr), have higher and radio modulation at the orbital period were subse- fields. This may be an indication that accretion-driven quently detected. field reduction or decay has not had enough time to operate in HMXBs but has in LMXBs. However, there does not seem to be any evidence of decaying fields in either the LMXB or HMXB populations; e.g. smaller magnetic VII. SUMMARY fields in older systems. On the other hand, CCOs are very young so if they acquired their low fields through mass Neutron stars are found in a wide variety of sources, fallback accretion, the field submergence would have had displaying an amazing array of behavior. They can be to operate on much faster timescales than it apparently isolated or in binary systems, accreting, heating, cooling, does in LMXBs. But as we continue to find new species spinning down, spinning up, pulsing, flaring and burst- in the NS zoo, one of these may someday be the “Rosetta ing. The one property that seems to determine their Stone” that will give us the clues for solving these puz- behavior most strongly is their magnetic field strength, zles. structure and evolution. The hot polar caps, bursts and I would like to thank Megan Decesar, Isabelle Gre- flares of magnetars are likely due to the rapid decay and nier, Peter Meszaros, Tod Strohmayer and Bing Zhang twisting of their superstrong magnetic fields, whose very for providing helpful comments on the manuscript.

[1] Abdo, A. A. et al. 2009a, ApJ, 696, 1084. [17] Arons, J., & Scharlemann, E. T. 1979, ApJ, 231, 854. [2] Abdo, A. A. et al. 2009b, Science, 325, 840. [18] Arons, J. 1983, ApJ, 266, 215. [3] Abdo, A. A. et al. 2009c, Science, 325, 848. [19] Arons, J. 2011, Space Sci Rev, High-Energy Emission [4] Abdo, A. A. et al. 2009d, Science 326, 1512. from Pulsars and their Systems, Astrophysics and Space [5] Abdo, A. A. et al. 2009e, ApJ, 701, L123. Science Proceedings, Springer-Verlag Berlin Heidelberg, [6] Abdo, A. A. et al. 2009f, ApJ, 706, L56. p. 165. [7] Abdo, A. A. et al. 2010, ApJS, 187, 460. [20] Baade, W. & Zwicky, F. 1934, PhRv 46 76. [8] Abdo, A. A. et al. 2011a, Science, 331, 739. [21] Baring, M. G. 2004, in Young Neutron Stars and Their [9] Abdo, A. A. et al. 2011b, ApJ, 736, L11. Environments eds F. Camilo and B. M. Gaensler (IAU [10] Abdo, A. A. et al. 2012, Science, 335, 189. Symposium 218, ASP Conference Proceedings), p. 267. [11] Abdo, A. A. et al. 2013, in preparation. [22] Baring, M. G. & Harding A. K. 2007, Astr. Space Sci., [12] Aharonian, F. A. et al. 2005, A & A, 442, 1. 308, 109. [13] Aliu, E. et al., 2011, Science, 334, 69. [23] Basko, M. M. & Sunyaev, R. A. 1976, MNRAS 175 395. [14] Alpar, M. A.; Cheng, A. F.; Ruderman, M. A.; Shaham, [24] Becker, P. A. & Wolﬀ, M. T. 2005, ApJ, 630, 465 J. 1982, Nature 300 28. [25] Becker, P. A. & Wolﬀ, M. T. 2007, ApJ, 654, 435 [15] Araya, R. & Harding, A. K. 1999, ApJ 517 334. [26] Beloborodov, A. 2009, ApJ, 703, 1044. [16] Archibald, A. et al. 2009, Science, 324, 1411. [27] Bignami, G. F.; Caraveo, P. A.; De Luca, A.; 14

Mereghetti, S. 2003, Nature, 423, 725 [65] Ho, W. C. 2007, MNRAS, 380, 71. [28] Buehler, R. et al. 2012, ApJ, 749, 26. [66] Ho, W. C. & Heinke, C. O. 2009, Nature, 462, 71. [29] Camilo, F., Ransom, S. M., Halpern, J. P., & Reynolds, [67] Ho W. C. G., Kaplan D. L., Chang P., van Adelsberg J. 2007, ApJ, 666, L93 M., Potekhin A. Y., 2007, MNRAS, 375, 821 [30] Camilo, F., Ransom, S. M., Halpern, J. P., et al. 2006, [68] Hobbs, G. et al. 2010, Classical and Quantum Gravity, Nature, 442, 892 27:8. [31] Coll, B. & Tarantola, A. 2009, [arXiv:0905.4121]. [69] Hurley, K. et al. 1999, ApJ, 510, L110. [32] Chakrabarty, D. 2005, in Binary Radio Pulsars, ASP [70] Kalapotharakos, C., Kazanas, D., Harding, A., & Con- Conference Series, Vol. 328, ed. F. A. Rasio and I. H. topoulos, I. 2012, ApJ, 749, 2 Stairs, p. 279. [71] Kaplan, D. L. 2008, 40 YEARS OF PULSARS: Mil- [33] Cheng, K. S., Ho, C., & Ruderman, M. A. 1986, ApJ, lisecond Pulsars, Magnetars and More. AIP Conference 300, 500. Proceedings, Volume 983, pp. 331-339 [34] Cline, T. L. et al. 1980, ApJ Letters, 237, L1. [72] Kaplan, D. L. & van Kerkwijk, M. H. 2009, ApJ, 705, [35] Cumming, A., Zweibel, E., & Bildsten, L. 2001, ApJ, 798. 557, 958. [73] Kaspi, V. et al. 2003, ApJ, 588L, 93. [36] Daugherty, J. K., & Harding, A. K. 1982, ApJ 252 337. [74] Kaspi, V.M., Roberts, M.S.E., & Harding, A.K. 2005, [37] Daugherty, J. K., & Harding A. K. 1996, ApJ, 458, 278. in Compact Stellar X-ray Sources, ed. W.H.G. Lewin & [38] DeJager, O. C. 2005, AIP Conference Proc. 801, 298. M. van der Klis, [39] DeJager, O. C. & Harding, A. K. 1992, ApJ, 396, 161. [75] Katz, J. I. 1982, ApJ, 260, 371. [40] denHartog, P. R., Kuiper,L., Hermsen,W. & Vink,J [76] Keane E. F., Kramer M., 2008, MNRAS, 391, 2009. 2004, Astron. Tel. 293. [77] Keane, E. F. & McLaughlin, M. A. 2011, Bull. Astr. [41] den Hartog, P. R.; Kuiper, L.; Hermsen, W.; Kaspi, V. Soc. India, 39, 333. M.; Dib, R.; Kndlseder, J.; Gavriil, F. P. 2008, A&A, [78] Kouveliotou, C. et al. 1998, Nature, 393, 235. 489, 245. [79] Kuiper,L., Hermsen,W. & Mendez, 2004, ApJLett, 613, [42] Dyks, J., Harding, A. K. & Rudak, B. 2004, ApJ, 606, 1173. 1125. [80] Lai, D. 2005, ASTROPHYSICAL SOURCES OF HIGH [43] Duncan, R. C. & Thompson, C. 1992, ApJ, 392,9. ENERGY PARTICLES AND RADIATION. AIP Con- [44] Duncan, R. C. 2000, Fifth Huntsville Gamma-Ray Burst ference Proceedings, Volume 801, pp. 259-264. Symposium, Eds: R. Marc Kippen, Robert S. Mallozzi, [81] Li, J., Spitkovsky, A., & Tchekhovskoy, A. 2012, ApJ, Gerald J. Fishman. AIP Vol. 526 (American Institute of 746, 60. Physics, Melville, New York) p.830 [astro-ph/0002442]. [82] Manchester, R. N., Hobbs, G. B., Teoh, A. & Hobbs, [45] Duncan, R. C & Thompson, C. 2001, ApJ, 561, 980. M., AJ, 129, 1993-2006 (2005) [46] Ghosh, P., & Lamb, F. K. 1979, ApJ, 234, 296. [83] Mendez, M. & Belloni, T. 2007, MNRAS 381, 790. [47] Goldreich, P., & Julian, W. H. 1969, ApJ, 157, 869. [84] McLaughlin M. A. et al., 2006, Nature, 439, 817. [48] Gotthelf, E. V. & Halpern, J. P. 2009, ApJ, 695, L35. [85] McLaughlin M. A. et al., 2007, ApJ, 670, 1307. [49] Gotthelf, e. V. & Halpern. J. P. 2010, ApJ, 709, 436. [86] Mereghetti, S., et al. 2005, A& A Lett., 433, L9. [50] Halpern, J. P. & Ruderman, M. A. 1993, ApJ, 415, 286. [87] Miller, M. C., Lamb, F. K. & Psaltis, D. 1998, ApJ, [51] Harding, A. K. 1983, Nature, 303, 23. 508, 791. [52] Harding 2005, Proc. 22nd Texas Symposium on Rel- [88] Molkov, S., et al. 2005, A& A Lett., 433, L13. ativistic Astrophysics, Edited by Pisin Chen, Elliott [89] Muslimov, A. G., & Harding, A. K. 2004, ApJ, 606, Bloom, Greg Madejski, and Vahe Patrosian, p.149-158 1143 [arXiv:astro-ph/0503300] [90] Palmer, D. M., et al. 2005, Nature, 434, 1107. [53] Harding, A. K. & Lai, D. 2006, Rep. Prog. Phys., 69, [91] Pavlov, G. G. & Yakovlev, D. G. 1976, Sov. Phys. JETP 2631. 43 389. [54] Harding, A. K., Meszaros, P., Kirk, J. G. & Galloway, [92] Petri, J. & Kirk, J. G. 2005, ApJ, 627, 37. D. J. 1984, ApJ, 278 369. [93] Pletsch, H. J. et al. 2012, ApJ, 744, 105 [55] Harding, A. K. & Muslimov, A. G. 2001, ApJ, 556, 987. [94] Podsiadlowski, Ph.; Rappaport, S.; Pfahl, E. D. 2002, [56]Harding,A.K.,J.V.Stern,J.Dyks&M.Frackowiak, ApJ, 565, 1107 2008, ApJ, 680, 1376. [95] Pons et al. 2009, A & A, 496, 207. [57] Heindl, W. A., Coburn, W., Gruber, D. E., et al., 1999, [96] Potekhin, A. Y.; Suleimanov, V. F.; van Adelsberg, M.; ApJ, 521, L49 Werner, K. 2012, [arXiv1208.6582]. [58] Heindl, W. A. et al. 2004, X-ray Timing 2003: Rossi [97] Ramaty, R., Bonazzola, S., Cline, T. L., Kazanas, D., and Beyond. AIP Conference Proceedings, Vol. 714, M´esz´aros, P. & Lingenfelter, R. E. 1980, Nature 287 held 3-5 November, 2003 in Cambridge, MA. Edited 122. by Philip Kaaret, Frederick K. Lamb, and Jean H. [98] Ransom, S. M. et al. 2011, ApJ, 727, L16. Swank. Melville, NY: American Institute of Physics, [99] Rea, N. et al. 2010, Science, 330, 944. 2004., p.323-330 [100] Rea, N.; Pons, Jos A.; Torres, D.o F.; Turolla, R. 2012, [59] Heinke, C. O. & Ho, W. C. 2010, ApJ, 719, L167. ApJ, 748, L12. [60] Heyl, J. S. & Kulkarni, S. R. 1998, ApJ 506, L61. [101] Romani, R. 1993, in Isolated Pulsars, van Riper, K., [61] Heyl, J. & Hernquist, L. 2005, ApJ, 618, 463. Epstein, R., & Ho, C., eds., 66 [62] Hewish, A., Bell, S. J., Pilkington, J. D., Scott, P. F. & [102] Romani, R. W. & Yadigaroglu, I. A. 1995, ApJ, 438, Collins, R. A. 1968, Nature 217 709. 314. [63] Hirotani, K. 2006, ApJ, 652, 1475. [103] Ruderman, M. 2004, in The Electromagnetic Spec- [64] Ho, W. C. 2011, MNRAS, 414, 2567 trum of Neutron Stars (NATO-ASI Proceedings) eds. 15

A. Baykal et al. (astro-ph/0410607). [113] Thompson, C. & Duncan, R. C. 1995, MNRAS 275 255. [104] Schonherr, G., et al. 2007, A&A, 472, 353 [114] van der Klis, M. 2000, ARA& A, 38, 717 [105] Shternin, Peter S.; Yakovlev, Dmitry G.; Heinke, Craig [115] Vasisht, G. & Gotthelf, E. V. 1997, ApJ, 486, L129. O.; Ho, Wynn C. G.; Patnaude, Daniel J. 2011, MN- [116] Wagoner, R. V. 1984, ApJ, 278, 345 RAS, 412, L108 [117] Watts, A. 2012, Annual Review of Astronomy and As- [106] Strohmayer, T. E., et al. 1996, ApJ, 469, L9 trophysics, vol. 50, p.609-640. [107] Strohmayer, T. E., & Markwardt, C. B. 1999, ApJ, 516, [118] Wijnands, R., & van der Klis, M., 1998, Nature, 394, L81 344. [108] Strohmayer, T., & Bildsten, L. 2004, in Compact Stellar [119] Woods, P. M. & Thompson, C. 2006, in Compact Stellar X-Ray Sources, ed. W. H. G. Lewin & M. van der Klis X-ray Sources, eds. W.H.G. Lewin and M. van der Klis, (Cambridge: Cambridge U. Press), (astro- ph/0301544) Cambridge Astrophysics Series, No. 39. Cambridge, UK: [109] Tavani, M. et al. 2011, Science, 331, 736. Cambridge University Press, p. 547. [astro-ph/0406133]. [110] Thompson, D. J. 2001, HIGH ENERGY GAMMA-RAY [120] Zavlin, V. E.; Pavlov, G. G.; Shibanov, Y. A.; Ventura, ASTRONOMY, AIP Conference Proceedings, Volume J. 1995, A & A, 297, 441 558, pp. 103-114 [121] Zavlin, V. E.; Pavlov, G. G.; Sanwal, D.; Trmper, J. [111] Thompson & Beloborodov 2005, ApJ, 634, 565. 2000, ApJ, 540, L25. [112] Thompson, C. & Duncan, R. C. 1996, ApJ 473 332. [122] Zhang, W.; Strohmayer, T. E.; Swank, J. H. 1997, ApJ, 482, L167.