A publishing partnership

The following article is Open access

The Third Fermi Large Area Telescope Catalog of Gamma-Ray Pulsars

, , , , , , , , ,

Published 2023 November 27 © 2023. The Author(s). Published by the American Astronomical Society.
, , Citation D. A. Smith et al 2023 ApJ 958 191DOI 10.3847/1538-4357/acee67

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

0004-637X/958/2/191

Abstract

We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems colocated with LAT sources also likely harbor gamma-ray MSPs. This catalog thus reports roughly 340 gamma-ray pulsars and candidates, 10% of all known pulsars, compared to ≤11 known before Fermi. Half of the gamma-ray pulsars are young. Of these, the half that are undetected in radio have a broader Galactic latitude distribution than the young radio-loud pulsars. The others are MSPs, with six undetected in radio. Overall, ≥236 are bright enough above 50 MeV to fit the pulse profile, the energy spectrum, or both. For the common two-peaked profiles, the gamma-ray peak closest to the magnetic pole crossing generally has a softer spectrum. The spectral energy distributions tend to narrow as the spindown power decreases to its observed minimum near 1033 erg s−1, approaching the shape for synchrotron radiation from monoenergetic electrons. We calculate gamma-ray luminosities when distances are available. Our all-sky gamma-ray sensitivity map is useful for population syntheses. The electronic catalog version provides gamma-ray pulsar ephemerides, properties, and fit results to guide and be compared with modeling results.

Export citation and abstractBibTeXRIS

1. Introduction

Fewer than a dozen gamma-ray pulsars were known when the Fermi Gamma-ray Space Telescope was launched on 2008 June 11, and the extent and diversity of the population and its role in Galactic dynamics were subject to debate (Thompson 2008). Fermi's primary instrument, the Large Area Telescope (LAT; Atwood et al. 2009), quickly established that the gamma-ray population is large and varied and is the dominant GeV gamma-ray source class in the Milky Way (Abdo et al. 2010a; The First Fermi LAT Catalog of Gamma-ray Pulsars, hereafter 1PC). The 46 pulsars in 1PC (6 months of data) grew to 132 in 2PC (3 yr of data), the Second Fermi-LAT gamma-ray pulsar catalog (Abdo et al. 2013). This third gamma-ray pulsar catalog (based on 12 yr of data) characterizes 294 confirmed gamma-ray pulsars, and tabulates 33 millisecond pulsars (MSPs) for which gamma-ray pulsations have not yet been seen but likely will be once accurate rotation ephemerides are established. We further tabulate LAT sources likely to reveal new "spider" MSPs, and LAT sources colocated with known pulsars, some of which may ultimately reveal gamma-ray pulsations. Roughly 340 gamma-ray pulsars and candidates are thus now known, or about 10% of the >3400 currently known pulsars (see Table 1).

Table 1. Pulsar Varieties

CategoryCountSubcount
Known rotation-powered pulsars (RPPs) a 3436 
with measured erg s−1  762
 
MSPs ( < 30 ms)681 
with measured erg s−1  250
Field MSPs b  427
MSPs in globular clusters c  254
 
Gamma-ray pulsars in this catalog d 294 
Spectral fits (with free b parameter) f  255 (116)
Profile fits in ≥1, 2, 6 energy bands 236, 167, 28
 
Young gamma-ray pulsars150 
Radio-quiet e  70
Gamma-ray MSPs144 
Isolated, Binary 32, 112
Discovered in LAT blind searches 10
Radio-quiet 6
Black Widows, Redbacks: 32, 13
 
Radio MSPs discovered in LAT sources119 
with gamma-ray pulsations 78
waiting for ephemeris phase-connection d  33

Notes.

a Includes the 3359 pulsars, which are all RPPs, in psrcat, the ATNF Pulsar Catalog (v1.69, Manchester et al. 2005), and as-yet unpublished discoveries. b http://astro.phys.wvu.edu/GalacticMSPs. c http://www.naic.edu/~pfreire/GCpsr. d Table 6 lists 39 MSPs discovered in radio searches of bright 4FGL sources with pulsar-like spectra. At least six were serendipitous. The rest will likely show pulsations once radio timing allows gamma-ray phase-folding. Table 5 lists additional pulsars colocated with 4FGL sources, some of which may reveal pulsations in the future. Table 15 lists 13 "spider" MSP candidates colocated with LAT sources. The number of detected gamma-ray pulsars thus likely exceeds 340, including unpulsed detections. e S1400 < 30 μJy, where S1400 is the radio flux density at 1400 MHz. f Sections 5 and 6 describe the pulse profile fits and energy spectral fits, respectively.

Download table as:  ASCIITypeset image

These results build on much previous work. GeV pulsations from the Crab were glimpsed using a balloon-borne instrument at the start of the 1970s (McBreen et al. 1973), followed by the pulsed detection of Vela by the SAS-2 satellite (Thompson et al. 1975). The COS-B satellite improved the measurements (Swanenburg et al. 1981). In the 1990s, EGRET on the Compton Gamma-Ray Observatory (CGRO) saw six GeV pulsars (Thompson et al. 1999), while a rare pulsar that is brighter below 100 MeV than above, PSR B1509−58 was detected with COMPTEL on CGRO (Kuiper et al. 1999; Abdo et al. 2010b). EGRET data revealed three other strong candidates: PSRs J0659+1414 (Ramanamurthy et al. 1996) and J1048−5832 (Kaspi et al. 2000), and the MSP PSR J0218+4232 (Kuiper et al. 2000). AGILE discovered gamma-ray pulsations from PSR J2021+3651 before Fermi's launch (Halpern et al. 2008). All 11 of these pulsars were quickly confirmed using LAT data, and are noted in Figure 1. Geminga, seen with EGRET, was undetected at radio wavelengths (Bignami & Caraveo 1996; Abdo et al. 2010c) and has turned out to be the prototype of about half of the young gamma-ray pulsars. The 294 pulsars reported here are more numerous than the 271 sources, all object classes combined, in the third EGRET source catalog (Hartman et al. 1999).

Figure 1. Refer to the following caption and surrounding text.

Figure 1. Cumulative number of known gamma-ray pulsars, beginning with the launch of Fermi. The crosses show the numbers included in the first (1PC) and second (2PC) catalogs of LAT pulsars and their publication dates. Some key discoveries are highlighted. See also Table 1.

Standard image High-resolution image

Figure 1 shows that the discovery rate since launch is steady. Table 1 breaks the numbers down by category. Long-term radio observations by the "Pulsar Timing Consortium" (Smith et al. 2008) enabled about half the discoveries and, importantly, also allowed for an unbiased sample of pulsars not seen in gamma-rays (Smith et al. 2019). The "Pulsar Search Consortium" (Ray et al. 2012), later joined by FAST (Li et al. 2018a; Wang et al. 2021) and the TRAPUM 107 project on MeerKAT (Clark et al. 2023a), discovered large numbers of radio MSPs at the positions of unidentified gamma-ray sources, leading to the subsequent detection of gamma-ray pulsations. Gamma-ray blind searches of unidentified sources revealed radio-quiet pulsars that make up a quarter of the current sample (see, e.g., Clark et al. 2017; Wu et al. 2018).

The discovery rate is sustained by innovations in how we detect pulsations. Ever-improving blind search algorithms (Pletsch & Clark 2014; Nieder et al. 2020a) allowed for the first discovery of a radio-quiet MSP (Clark et al. 2018). The increasingly sophisticated use of an optical companion's orbit reduces the parameter space searched using Einstein@Home 108 to discover gamma-ray MSPs in binary systems with perturbed orbits (Nieder et al. 2022). Another example is a method that allows photon weighting (see Section 2) even if the unpulsed gamma-ray source is undetected (Bruel 2019; Smith et al. 2019). Underlying all analysis efforts are the improved sensitivity and low-energy reach afforded by the Pass8 reconstruction method (Atwood et al. 2013; Bruel et al. 2018).

As a result, not just the numbers but the variety of gamma-ray pulsars continues to grow. The minimum spindown power is now 20× lower than the pre-launch expectation of 1034 erg s−1 (Smith et al. 2008). 109 The fastest known field MSP, PSR J0952−0607, was found at the location of a gamma-ray source and subsequently timed with LAT data (Bassa et al. 2017b; Nieder et al. 2019). A third globular cluster, NGC 6652, was found to have gamma-ray emission dominated by a single MSP, PSR J1835−3259B (Gautam et al. 2022b; Zhang et al. 2022), with a fourth recently reported, PSR J1717+4308A in M92 (Zhang et al. 2023). The gamma-ray flux and pulse profile of PSR J2021+4026 in the γ Cygni supernova remnant (SNR) changed during mode transitions in 2011 and 2018 (Allafort et al. 2013; Razzano et al. 2023). The LAT sees more than 40 "spider" MSPs and a dozen candidates (see Section 7), compact binary systems where the pulsar wind ablates its companion. Spiders fall into two categories: "black widows" have companion masses 0.01M < Mc < 0.05 M and orbital periods PB < 10 hr whereas "redbacks" have Mc > 0.2 M and PB < 1 day. Gamma-ray timing of 35 stable MSPs for over 12 yr usefully constrains the intensity of gravitational waves from supermassive black hole binaries in the hearts of distant galaxies. The upper limit should become a measurement in the coming years (Ajello et al. 2022). Following the methods initially applied to PSR J1555−2908 (Nieder et al. 2022), we may be poised to detect planets in multiyear orbits in several compact binary MSP systems.

This heterogeneous population can be classified by comparing the spin period (P) and the period derivative (), shown in Figure 2. Throughout this paper, we call pulsars in the main population "young" to distinguish them from the much older MSPs, thought to be spun up to rapid periods via accretion from a companion (Alpar et al. 1982), although, e.g., the accretion-induced collapse of white dwarfs might also create MSPs (Gautam et al. 2022a).

Figure 2. Refer to the following caption and surrounding text.

Figure 2. Pulsar spindown rate, , vs. the rotation period P. Green dots indicate young, radio-loud (RL) gamma-ray pulsars and blue squares show "radio-quiet" (RQ) pulsars, defined as S1400 < 30 μJy, where S1400 is the radio flux density at 1400 MHz. Red triangles are millisecond gamma-ray pulsars. Black dots indicate pulsars phase-folded in gamma-rays without significant pulsations. Phase-folding was not done for pulsars shown by gray dots. Orange triangles are radio MSPs discovered at the positions of previously unassociated LAT sources, hidden by red triangles when gamma pulsations were subsequently found. The rest are listed in Table 6, and plotted with when is unavailable. The solid black diagonal is the radio deathline of Equation (4) of Zhang et al. (2000). Shklovskii corrections to have been applied only to gamma-ray MSPs with measured proper motion (see Section 4.3).

Standard image High-resolution image

All known gamma-ray pulsars are rotation-powered pulsars (RPPs): LAT has not yet detected accretion-powered pulsars nor the magnetars that populate the upper-right portion of the plane, for which the dominant energy source is magnetic field decay (Parent et al. 2011). An interesting exception is an LAT detection of a few photons for a few minutes from an extragalactic magnetar giant flare (Ajello et al. 2021a). The locations of all 294 gamma-ray pulsars on the sky are shown in Figure 3. The diagram shows diagonal lines of constant , τc , and BS derived from the timing information as follows. For an orthogonal rotator, the magnetic field on the neutron star surface at the magnetic equator (the rotation pole) is . The "characteristic age" assumes that magnetic dipole braking is the only energy-loss mechanism, that the magnetic moment and inclination do not change, and that the initial spin period was much less than the current period. τc thus approximates true age well for some young pulsars, and poorly for MSPs. We set the neutron star radius to RNS = 10 km, and c is the speed of light in a vacuum.

Figure 3. Refer to the following caption and surrounding text.

Figure 3. Pulsar sky map in Galactic coordinates (Hammer projection). Symbols are the same as in Figure 2.

Standard image High-resolution image

The fourth Fermi-LAT source catalog (Abdollahi et al. 2020), and specifically Data Release 3 (DR3; Abdollahi et al. 2022, hereafter 4FGL) 110 characterizes 6658 point and extended sources using 12 yr of LAT data. Half of the sources are various blazar classes of active galactic nuclei, but a third remain unassociated with objects known at other wavelengths. Radio and gamma-ray pulsation searches at the positions of unidentified sources have yielded fully half of the gamma-ray pulsars. The discoveries continue, as detailed in Sections 3 and 7. The 4FGL spectral, flux, and variability measurements are used throughout this work for the pulsar searches and characterization.

We provide the pulsar catalog in FITS 111 and spreadsheet file formats, along with other information here as well as on the Fermi Science Support Center (FSSC) servers. 112 Appendix D provides a description of this material.

2. Observations

Atwood et al. (2009) described the Fermi LAT, and Abdo et al. (2009a), Ackermann et al. (2012a), and Ajello et al. (2021b) reported on-orbit performance. Fermi carries another instrument, the Gamma-ray Burst Monitor (Meegan et al. 2009), which was not used to prepare this catalog.

The LAT is a pair-production telescope composed of a 4 × 4 grid of towers. Each tower consists of a stack of tungsten foil converters interleaved with silicon-strip particle tracking detectors, mated with a hodoscopic cesium-iodide calorimeter. A segmented plastic scintillator anticoincidence detector covers the grid to help discriminate charged particle backgrounds from gamma-ray photons. The LAT field of view is ∼2.4 sr. For most of the Fermi mission, the primary operational mode has been a sky survey where the satellite rocks between a pointing above the orbital plane and one below the plane after each orbit. In this mode, the entire sky is imaged every two orbits (∼3 hr) and any given point on the sky is observed ∼1/6th of the time. For 1 yr, beginning 2013 December, the normal survey mode was changed to favor exposure to the Galactic center region. 113 Survey mode was again modified after a solar panel rotation drive stopped moving on 2018 March 16, detailed in Ajello et al. (2021b).

The LAT is sensitive to gamma-rays with energies E from 20 MeV to over 300 GeV, with an on-axis effective area of ∼8000 cm2 above 1 GeV. Multiple Coulomb scattering of the electron–positron pairs created by converted gamma-rays degrades the per-photon angular resolution, with the average 68% containment radius varying as E−0.8 from 5° at 100 MeV to 0fdg1 at 10 GeV. 114

This energy-dependent point-spread function, the backgrounds from the complex diffuse emission and nearby sources, and the source spectrum are encapsulated in the photon weights (wi , Kerr 2011), which give the probability that a photon originates from a particular source. Weights are used to optimize pulsed signal significance while minimizing event selection trials penalties. This powerful tool has been extended ("simple weights" and "model weights") to sky locations with no point source (Bruel 2019). Weighting is used for all of the discovery techniques (Section 3), and to characterize the pulse profiles (Section 5).

Events recorded by the LAT have time stamps derived from GPS clocks integrated into the satellite's Guidance, Navigation, and Control (GNC) subsystem, accurate to ≈ 300 ns relative to UTC (Abdo et al. 2009a; Ajello et al. 2021b). GNC provides the instantaneous spacecraft position with ≈60 m accuracy. Generally, we compute pulsar rotational phases ϕi using Tempo2 (Hobbs et al. 2006) with the fermi plug-in (Ray et al. 2011), or with PINT (Luo et al. 2021). These use the recorded times and spacecraft positions combined with a pulsar rotational ephemeris (specified in a Tempo2 parameter, or "par," file). The timing chain from the instrument clocks through the barycentering, and epoch folding software is accurate to better than a microsecond (Smith et al. 2008).

We use different data selections and processing for the various analyses presented here, though all cases make use of "Source" class events reconstructed using Pass8 revision 3 (Atwood et al. 2013; Bruel et al. 2018). The pulsation searches and pulsar timing described in Section 3 generally make use of all available data at the time of analysis and employ a variety of data processing methods. The effects of this heterogeneity on, e.g., pulsar discovery efficiency or pulse profile inference are minor, and we do not attempt to provide further details.

For the pulse profiles (see Section 5), the data span 11.8 yr, from MJD 54682 (2008 August 4, when LAT data taking began) to MJD 59000 (2020 May 31). We exclude gamma-rays collected when the LAT was not in nominal science operations mode or when the spacecraft rocking angle exceeded 52°, and we accept photons with reconstructed energies from 0.05 to 100 GeV, within 15° of the pulsar positions. We used Model Weights, requiring weights wi > 0.001.

The spectral measurements (see Section 6) use the same data set as 4FGL-DR3. In brief, the data span 12 yr, to MJD 59063 (2020 August 2) and are selected with energies 50 MeV to 1 TeV. The catalog analysis further uses a heterogeneous zenith angle cut ranging from 80° for energies below 100 MeV to 105° above 1 GeV.

3. Discovery and Timing

Inclusion in the main catalog requires a statistically significant pulsation detection in the Fermi LAT gamma-ray data. The following subsections describe the paths to detection, as well as a brief description of gamma-ray timing that often strengthens the initial signal.

In all cases, detecting and characterizing pulsations requires a rotation ephemeris (or a "timing model") to convert photon arrival times ti to neutron star rotational phases ϕi . This "folding," which "stacks" photons at the same rotational phase, allows a pulse to rise above the background since ≪1 photons are collected per pulse even from the brightest pulsars (Kerr 2022). We detect pulsations with the H-test (de Jager et al. 1989; de Jager & Büsching 2010), a statistical test for discarding the null hypothesis that a set of photon phases is uniformly distributed. For Nγ gamma-rays, the m = 20 harmonic weighted version of the H-test statistic (Kerr 2011) is

with

and αkw and βkw are the empirical trigonometric coefficients and . The w subscripts indicate that this is a photon-weighted version of the test, and the wi are the photon weights, evaluated via spectral analysis (see also 2PC, Abdo et al. 2013). For m = 20, which we adopt universally, the cumulative distribution function for H in the asymptotic limit is (Kerr 2011), and 3σ, 4σ, and 5σ thresholds correspond to = 14, = 24, and = 36. The H-test is unbinned, well suited to the extremely sparse gamma-ray pulsar data: the LAT often detects only one photon in tens of thousands (millions, for MSPs) of pulsar rotations. Bruel (2019) gave corrections for small photon counts. Pulsars with narrow, sharp peaks are easier to detect than pulsars with broad peaks (see Figure 5 of Hou et al. 2014). All pulsars in our sample are detected with = 8. Using = 20 incurs little computational cost and insures sensitivity to putative exotic profile shapes. Kerr (2011) also showed that large m does not cause false-positive detections.

Figure 4 highlights some aspects of gamma-ray phase folding. Figure 9 and Appendix B show other example profiles. The top-most frame shows a weighted phase histogram, duplicated over a second rotation. Section 5 describes the profile fit overlaid in blue. The phase-aligned 1.4 GHz radio pulse overlaid in red comes from the radio timing observations used to create the rotation ephemeris, in this specific case by Parthasarathy et al. (2019). The horizontal dashed line shows the gamma-ray background level, estimated from the photon weights as , the sum of the expected contribution to the weights of background photons not associated with the pulsar. A phase histogram baseline exceeding the background level may indicate the presence of unpulsed magnetospheric emission. Section 5.1 gives details.

Figure 4. Refer to the following caption and surrounding text.

Figure 4. Top: gamma-ray phase histogram for PSR J1648-4611 discovered by Kramer et al. (2003), overlaid with the 1.4 GHz profile (red) obtained during Parkes radio telescope timing (Parthasarathy et al. 2019). The blue curve shows a fit to the histogram, and the horizontal dashed line is an estimate of the background level (see Section 5.1). Bottom-left: phaseogram over the course of the mission. Dots indicate photons, with the grayscale set according to the photon weight. Right: H-test significance accumulated over the course of the mission. The green horizontal dashed line shows when the ephemeris validity finishes (see Section 3.1).

Standard image High-resolution image

The next frame below it shows the phase drifting after the last radio time of arrival (ToA) used to model the neutron star rotation, indicated by the green horizontal dashed line. The start of this timing model's "validity" range, before Fermi's launch, is not shown. Pulsars with irregular spindown, as for this young high- pulsar, require extra parameters to model the rotation, and accuracy of the extrapolation past the validity range rapidly degrades. Stable pulsars can be modeled with few parameters, often accurately predicting the neutron star rotation for years before and/or after validity. The right-hand frame shows the weighted H-test increasing as data accumulated over the years, a nearly straight line for most pulsars. Changes in slope can result from phase drifts, as in this case, or from increased background due to, e.g., a nearby flaring blazar (see Section 6.6), or from changes in the LAT's exposure to the pulsar's sky position. Exposure per unit time increased for this pulsar during 2014 (mid-year was MJD 56810) when LAT pointed more frequently toward the Galactic center, visible in the time versus phase plot, without however affecting the H-test growth. Smith et al. (2019) showed simulations for which Poisson fluctuations in the very low rate of photon arrivals for pulsars near detection threshold significantly perturb the H-test time evolution. Slope variations due to pulsar flux changes nearly never happen: rare exceptions are the young pulsar PSR J2021+4026 (Allafort et al. 2013; Razzano et al. 2023) and transitional MSPs like PSR J1023+0038 (see Stappers et al. 2014, and Appendix A) and J1227−4853 (Johnson et al. 2015; Roy et al. 2015), which lack detectable gamma-ray pulsations during their accretion states.

If different analyses were tried, the number of trials must be accounted for in the pulsation significance calculation, requiring a higher H-test value to claim a detection. Before the advent of photon weighting, we sometimes varied the minimum energy and/or the angular extent of our data set, to explore the pulsar's spectral hardness and local background level. Using a LAT source's measured spectrum to calculate weights allows one single trial for this parameter space.

When the source spectrum is too faint to measure, we use "simple weights" with a parameter μw of unknown optimal value, and thus, possible additional trials. The μw parameter is the logarithm of the energy at which the log-Gaussian weighting function peaks. Since Smith et al. (2019), we have been using six trials for simple weight gamma-ray pulsation searches using radio or X-ray ephemerides, and an H-test > 25 ( = 4.7 × 10−5, or > 4.1σ) detection threshold. We fold the entire data set three times, using μw = (3.2, 3.6, 4.0), and three more times restricting the data to the ephemeris validity period. Although these six trials are not independent, a conservative estimate for the chance probability for a false-positive detection with an H-test > 25 threshold is 6 × 4.7 × 10−5, or 0.3 for a sample of 1000 pulsars. In 2PC we required an H-test > 36 (>5σ). Smith et al. (2019) discovered 16 new gamma-ray pulsars, in part because of this refined, relaxed threshold. Several more pulsars in this catalog were later found in the same way.

This prescription may miss some pulsars if the ephemeris extrapolates well enough that the accumulation of additional data beyond the validity range yields H > 25, but poorly enough that the full data yield H < 25. Such cases also arise from positive fluctuations, so we avoid further trials factors by limiting our search to the two combinations just described.

As further explained in Section 3.4, when the ephemeris is established independently of the LAT, using radio or X-ray observations (Section 3.1), the pulsar can be up to 20 times fainter in gamma-rays than if an unidentified LAT point source (Section 3.2) guides a "blind" gamma-ray search. Once pulsations are found during some epoch of LAT data, gamma-ray timing can extend and improve the ephemeris (Section 3.5).

Using the above detection criteria, we report 294 gamma-ray pulsars in the breakdown tabulated in Table 1 with their distribution on the sky shown in Figure 3. The 150 young pulsars are named in Table 2, and the 144 MSPs are in Table 3. The new pulsars in 3PC were mostly more "difficult" to discover than those reported in the earlier catalogs. For pulsars found with just a few foldings using a radio ephemeris, discussed in Section 3.1, difficulty stems from their extreme faintness (<1 photon per month): an analysis using many trials (for example, from exploring different data intervals) yields significance too low to distinguish from statistical fluctuations. Furthermore, we continue to discover radio MSPs at the positions of unidentified FGL sources, described in Section 3.3. These are difficult both in radio, because of eclipses, scintillation, and intrinsic faintness, and in gamma-rays if their signal-to-noise ratio (S/N) is low and the ephemeris validity is brief, that is, if the radio observations cover only a few years of the LAT mission. The blind searches, discussed in Section 3.4, are now finding MSPs, including those with short, varying orbital periods, feared as perhaps impossible before launch (Ransom 2007). These searches require many orders of magnitude more computing power than for young pulsars. The catalog includes 44 "spider" systems, and in Section 7.1 we tabulate several candidates colocated with gamma-ray sources, for which pulsations will likely be seen in the coming years. Timing models for "noisy" pulsars require many parameters, which presents difficulty in maintaining a coherent timing solution, e.g., with frequent radio monitoring. Our methodical pursuit of ever more difficult gamma-ray pulsars in the LAT's unbiased all-sky data set means that our notion of gamma-ray pulsars is greatly enriched compared to 1PC.

Table 2. Some Parameters of Young LAT-detected Pulsars

PSRCodes l b P S1400
  (°)(°)(ms)(10−15)(erg s−1)(mJy)
J0002+6216GUr117.33−0.07115.46.0153.00.02
J0007+7303Gq119.6610.46315.9355.92445.0<0.005
J0106+4855GUr125.47−13.8783.20.429.40.01
J0117+5914Rr126.28−3.46101.45.8221.00.30
J0139+5814Rr129.22−4.04272.510.720.94.60
J0205+6449Xrx130.723.0865.7192.126688.00.05
J0248+6021Rr136.900.70217.155.2212.013.70
J0357+3205GUq162.76−16.01444.113.105.90<0.004
J0359+5414GUq148.230.8879.416.71317.0 
J0514−4408Rr249.51−35.36320.32.02.450.71
J0534+2200RErx184.56−5.7833.7420.2435319.014.00
J0540−6919Xrx279.72−31.5250.6478.9145529.00.10
J0554+3107GUq179.062.70465.0142.6056.0<0.066
J0622+3749GUq175.8810.96333.225.4227.1<0.012
J0631+0646GUr204.68−1.24111.03.6104.00.025a
J0631+1036Rr201.220.45287.8102.7170.01.11
J0633+0632GUq205.09−0.93297.479.57119.0<0.003
J0633+1746XExq195.134.27237.110.9732.5<0.507
J0659+1414RErx201.118.26384.955.038.02.70
J0729−1448Rr230.391.42251.7113.4280.00.83
J0729−1836Rr233.76−0.34510.218.95.631.90
J0734−1559GUq232.062.02155.112.51132.0<0.005
J0742−2822Rr243.77−2.44166.816.7141.026.00
J0744−2525GUq241.35−0.7392.01.048.3 
J0802−5613GUq269.98−13.19274.12.85.30 
J0834−4159Rr260.89−1.04121.14.395.10.28
J0835−4510RErx263.55−2.7989.4122.36763.01050.00
J0908−4913Rr270.27−1.02106.815.1490.020.00
J0922+0638Rr225.4236.39430.613.76.7710.00
J0940−5428Rr277.51−1.2987.632.81928.00.66
J1016−5857Rr284.08−1.88107.480.42563.00.90
J1019−5749Rr283.84−0.68162.520.1184.03.80
J1023−5746GUq284.17−0.41111.5379.8910820.0<0.030
J1028−5819Rr285.06−0.5091.414.2734.00.24
J1044−5737GUq286.571.16139.054.57801.0<0.020
J1048−5832REr287.430.58123.795.51992.09.10
J1055−6028Rr289.13−0.7499.729.51176.00.95
J1057−5226RErx285.986.65197.15.830.14.40
J1057−5851GUq288.610.80620.4100.616.6 
J1105−6037GUq290.24−0.40194.921.8116.0 
J1105−6107Rr290.49−0.8563.215.82475.01.20
J1111−6039GUq291.02−0.11106.7195.26346.0 
J1112−6103Rr291.22−0.4665.031.54537.02.30
J1119−6127Rrx292.15−0.54409.14042.42330.01.09
J1124−5916Rrx292.041.75135.5751.511914.00.08
J1135−6055GUq293.790.58114.578.232057.0<0.030
J1139−6247GUq294.79−1.06120.44.191.6 
J1151−6108Rr295.810.91101.610.3386.00.06
J1203−6242GUq297.52−0.34100.644.11709.0 
J1208−6238GUq297.99−0.18440.73309.571526.0<0.017
J1224−6407Rr299.98−1.41216.55.019.38.90
J1231−5113GUq299.7611.52206.40.10.525 
J1231−6511GUq300.87−2.40247.428.474.0 
J1253−5820Rr303.204.53255.52.14.984.10
J1341−6220Rr308.73−0.03193.4253.01379.02.70
J1350−6225GUq309.73−0.34138.28.9132.0 
J1357−6429Rrx309.92−2.51166.2354.43047.00.52
J1358−6025GUq311.111.3760.53.0536.0 
J1410−6132Rr312.19−0.0950.131.810000.01.90
J1413−6205GUq312.37−0.74109.727.39818.0<0.024
J1418−6058GUq313.320.13110.6171.004992.0<0.029
J1420−6048Rrx313.540.2368.282.410250.01.19
J1422−6138GUq313.52−0.66341.096.7996.4<0.060
J1429−5911GUq315.261.30115.823.88606.0<0.021
J1447−5757GUq317.851.51158.711.7115.0 
J1459−6053GUqx317.89−1.79103.225.26908.0<0.037
J1509−5850Rr319.97−0.6288.99.2516.00.21
J1513−5908XErx320.32−1.16151.61526.217290.01.43
J1522−5735GUq322.05−0.41204.362.46289.0<0.034
J1528−5838GUq322.17−1.75355.724.821.7 
J1531−5610Rr323.900.0384.213.8911.00.87
J1614−5048Rr332.210.17231.9492.01557.04.10
J1615−5137GUq331.76−0.54179.310.672.8 
J1620−4927GUq333.890.41171.910.4981.5<0.040
J1623−5005GUq333.72−0.3185.14.2266.0 
J1624−4041GUq340.566.15167.94.739.4 
J1641−5317GUq333.29−4.56175.13.727.2 
J1646−4346Rr341.110.98231.7111.8355.01.25
J1648−4611Rr339.44−0.79165.023.7208.00.61
J1650−4601GUq339.78−0.95127.115.1291.0 
J1702−4128Rr344.740.12182.252.3341.01.17
J1705−1906Rr3.1913.03299.04.16.115.66
J1709−4429RErx343.10−2.69102.594.83475.012.10
J1714−3830GUq348.440.1484.170.34657.0 
J1718−3825Rr348.95−0.4374.713.21248.01.70
J1730−3350Rr354.130.09139.584.11222.04.30
J1731−4744Rr342.56−7.67829.9163.511.327.00
J1732−3131GUr356.311.01196.528.0145.00.05
J1736−3422GUq354.33−1.18346.965.562.0 
J1739−3023Rr358.090.34114.411.4300.01.01
J1740+1000Rr34.0120.27154.121.3230.02.70
J1741−2054GUrx6.424.91413.717.09.480.16
J1742−3321GUq355.85−1.69143.31.317.0 
J1746−3239GUq356.96−2.17199.56.5632.6<0.034
J1747−2958Rrx359.31−0.8498.861.32506.00.25
J1748−2815GUq0.91−0.19100.23.5138.0 
J1757−2421Rr5.280.06234.112.739.27.20
J1801−2451Rr5.27−0.87125.089.51810.01.46
J1803−2149GUq8.140.19106.319.50640.0<0.024
J1809−2332Gq7.39−2.00146.834.39429.0<0.025
J1813−1246GUq17.242.4448.117.566238.0<0.017
J1816−0755Rr21.874.09217.66.524.80.17
J1817−1742GUq13.34−0.70149.720.6241.0 
J1826−1256Gxq18.56−0.38110.2120.953564.0<0.013
J1827−1446GUq17.08−1.50499.245.314.4 
J1828−1101Rr20.490.0472.114.81562.02.30
J1831−0952Rr21.90−0.1367.38.31078.00.35
J1833−1034Rr21.50−0.8961.9202.033623.00.07
J1835−1106Rr21.22−1.51165.920.6178.02.50
J1836+5925Gq88.8825.00173.31.5011.4<0.004
J1837−0604Rr25.960.2796.344.91982.00.75
J1838−0537GUq26.510.21145.8450.715746.0<0.017
J1841−0524Rr27.02−0.33445.8233.2103.00.20
J1844−0346GUq28.79−0.19112.9154.74249.0 
J1846+0919GUq40.695.34225.69.9334.2<0.005
(J18460258)b Xxq29.71−0.24326.67107.18055.0 
J1853−0004Rr33.09−0.47101.45.6210.00.70
J1856+0113Rr34.56−0.50267.5205.9424.00.19
J1857+0143Rr35.17−0.57139.831.0448.00.74
J1906+0722GUq41.220.03111.535.861020.0<0.021
J1907+0602Gr40.18−0.89106.686.52814.00.00
J1913+0904Rr43.50−0.68163.317.6159.00.40
J1913+1011Rr44.48−0.1735.93.42882.00.90
J1925+1720Rr52.180.5975.710.5954.00.07
J1928+1746Rr52.930.1168.713.21602.00.28
J1932+1916GUq54.660.08208.293.17407.0<0.075
J1932+2220Rr57.361.55144.557.0746.01.20
J1935+2025Rr56.05−0.0580.060.44638.00.50
J1952+3252RErx68.772.8239.55.83723.01.00
J1954+2836GUq65.240.3892.721.161048.0<0.005
J1954+3852Rr74.045.70352.96.65.931.07
J1957+5033GUq84.5811.01374.87.085.31<0.010
J1958+2846GUq65.88−0.35290.4211.89341.0<0.006
J2006+3102Rr68.67−0.53163.724.9223.00.27
J2017+3625GUq74.510.39166.71.411.6 
J2021+3651Rr75.220.11103.794.83351.00.10
J2021+4026Gq78.232.09265.355.60117.0<0.020
J2022+3842Xrx76.890.9648.686.229707.00.11
J2028+3332GUq73.36−3.01176.74.8634.8<0.005
J2030+3641RUr76.12−1.44200.16.532.00.15
J2030+4415GUq82.342.89227.15.0517.0<0.008
J2032+4127GUbr80.221.03143.211.6156.00.23
J2043+2740Rr70.61−9.1596.11.255.00.34
J2055+2539GUq70.69−12.52319.64.104.96<0.007
J2111+4606GUq88.31−1.45157.8162.631632.0<0.013
J2139+4716GUq92.63−4.02282.81.783.11<0.014
J2208+4056Rr92.57−12.11637.05.30.8070.04
J2229+6114Rrx106.652.9551.775.321557.00.25
J2238+5903GUq106.560.48162.796.85887.0<0.011
J2240+5832Rr106.57−0.11139.915.3219.02.70

Note. Column 2 gives the discovery and detection codes: G = discovered in Fermi-LAT gamma-ray data, R = discovered in the radio and/or gamma-ray pulsations detected using the radio ephemeris, X = discovered in the X-ray and/or gamma-ray pulsations detected using the X-ray ephemeris, E = pulsar was detected in gamma-rays by EGRET/COMPTEL, P = discovered by the Pulsar Search Consortium, U = discovered using a Fermi-LAT seed position, r = pulsations detected in the radio band, x = pulsations detected in the X-ray band, b = binary, and q = no radio detection. In Pletsch et al. (2013) and Kerr et al. (2015), PSR J1522−5735 erroneously has half the spin period shown here and 4× higher . Our larger gamma-ray photon sample made clear that the fundamental spin frequency is half of the discovery value, reflected in the ephemeris we distribute. Columns 3 and 4 give Galactic coordinates for each pulsar. Columns 5 and 6 list the period (P) and its first derivative (), and Column 7 gives the spindown luminosity . The Shklovskii correction to and is negligible for these young pulsars (see Section 4.3). Column 8 gives the radio flux density (or upper limit) at 1400 MHz (S1400, see Section 4.1). PSR J1509−5850 should not be confused with PSR B1509−58 (=J1513−5908) studied with the Compton Gamma-Ray Observatory. (a) J. Wu (2023, private communication). (b) Kuiper et al. (2018) and Section 7.2 describe the detection of pulsed <100 MeV gamma-rays from PSR J1846−0258.

Download table as:  ASCIITypeset images: 1 2 2

Table 3. Some Parameters of the LAT-detected Millisecond Pulsars

PSRCodes l b P S1400
  (°)(°)(ms)(10−20)(erg s−1)(mJy)
J0023+0923RUPbwr111.38−52.853.051.0815.90.73
J0030+0451Rrx113.14−57.614.871.023.491.09
J0034−0534Rbr111.49−68.071.880.5029.70.21
J0101−6422RUPbr301.19−52.722.570.4811.80.16
J0102+4839RUPbr124.87−14.172.961.1717.30.06
J0154+1833Rr143.18−41.812.360.288.470.11
J0218+4232REbrx139.51−17.532.327.74243.00.90
J0248+4230RUPr144.88−15.322.601.6837.80.06
J0251+2606RUPbwr153.88−29.492.540.7618.20.02
J0307+7443RUPbr131.7014.223.161.7121.70.04
J0312−0921RUPbrw191.51−52.383.701.9715.3 
J0318+0253RUr178.46−43.605.191.764.980.01
J0340+4130RUPr153.78−11.023.300.597.740.31
J0418+6635RUPr141.5211.542.911.3721.91.84
J0437−4715Rbrx253.39−41.965.765.7311.8150.20
J0533+6759RUPr144.7818.184.391.265.900.10
J0605+3757RUPbr174.198.022.730.479.170.14
J0610−2100Rbwr227.75−18.183.861.238.450.65
J0613−0200Rbr210.41−9.303.060.9613.22.25
J0614−3329RUPbrx240.50−21.833.151.7822.00.68
J0621+2514RUPbr187.125.072.722.4948.70.08
J0636+5128Rbwrx163.9118.642.860.345.761.00
J0653+4706RUPbr169.2619.784.752.087.610.10
J0737−3039ARbr245.24−4.5022.70176.005.943.06
J0740+6620Rbrx149.7329.602.881.2220.01.10
J0751+1807Rbrx202.7321.093.480.787.301.35
J0931−1902Rr251.0023.054.640.361.430.52
J0952−0607RUPbwr243.6535.381.410.4866.50.02
(J09553947)1 RUbrk269.9311.542.023.73178.0 
J0955−6150RUPbr283.68−5.741.991.4370.40.64
J1012−4235RUr274.2211.223.100.668.680.26
(J1023+0038)2 Rbrk243.4945.781.690.6855.915.01
J1024−0719Rrx251.7040.525.161.865.331.50
J1035−6720GUr290.37−7.842.874.6577.40.06
J1036−8317RUPbr298.94−21.503.413.0630.50.45
J1048+2339RUPbrk213.1762.144.673.0111.70.17
J1124−3653RUPbwr284.0922.762.410.5817.00.04
J1125−5825Rbr291.892.603.106.0980.61.00
J1125−6014Rbr292.500.892.630.378.121.32
J1137+7528RUPbr129.0040.772.510.327.960.10
J1142+0119RUPbr267.5459.405.081.504.520.05
J1207−5050RUPr295.8611.424.840.612.120.39
J1221−0633Rbr289.6855.531.931.0959.20.09
J1227−4853RUPbkr298.9713.801.691.33109.01.56
J1231−1411RUPbrx295.5348.393.682.1217.90.29
(J12598148)1 RUbrw303.24−18.952.090.3314.5 
J1301+0833RUPbwr310.8171.281.841.0566.5 
J1302−3258RUPbr305.5929.843.770.664.830.17
(J13066043)3 Rbr304.752.095.673.046.58 
J1311−3430GUPbwr307.6828.182.562.0949.00.11
J1312+0051RUPbr314.8463.234.231.719.150.19
J1327−0755Rbr318.3853.852.680.183.720.19
J1335−5656GUq308.895.433.241.2114.0 
J1400−1431Rbr326.9945.093.080.729.740.17
(J1402+1306)4 RUPbr356.4268.225.891.352.60 
J1431−4715Rbkr320.0512.252.011.4168.40.67
J1446−4701Rbwr322.5011.432.190.9836.60.46
J1455−3330Rbr330.7222.567.992.431.880.73
J1513−2550RUPbwr338.8226.962.122.1589.10.31
J1514−4946RUPbr325.256.813.591.8715.90.25
(J15262744)5 RUb340.2323.672.490.359.05 
J1536−4948RUPbr328.204.793.082.1228.60.09
J1543−5149Rbr327.922.482.061.6273.30.82
J1544+4937RUPbwr79.1750.172.160.3111.02.13
J1552+5437Rr85.5947.212.430.287.790.04
J1555−2908RUPbrw344.4818.501.794.45307.00.20
J1600−3053Rbr344.0916.453.600.958.052.44
J1614−2230Rbrx352.6420.193.150.9612.11.14
J1622−0315RUPbkr10.7130.683.851.147.93 
(J16236936)5 RUbr319.61−13.962.410.9125.5 
J1625−0021RUPbr13.8931.832.832.1337.00.19
(J1627+3219)6 RUb52.9743.212.180.5520.8 
J1628−3205RUPbkr347.4311.483.211.4814.2 
J1630+3734RUPbr60.2443.213.321.0711.60.02
J1640+2224Rbr41.0538.273.160.283.510.46
J1641+8049Rbwr113.8431.762.010.9846.80.15
J1649−3012GUq351.969.213.421.3213.0 
J1653−0158GUbwq16.6124.941.970.2412.4<0.008
J1658−5324RUPr334.87−6.632.441.1030.40.43
J1713+0747Rbr28.7525.224.570.853.538.30
J1730−2304Rr3.146.028.122.021.494.00
J1732−5049Rbr340.03−9.455.311.423.742.11
J1741+1351Rbr37.8921.643.753.0222.70.29
J1744−1134Rr14.799.184.070.895.212.60
J1744−7619GUq317.11−22.464.690.973.71<0.023
J1745+1017RUPbwr34.8719.252.650.256.450.51
J1747−4036RUPr350.21−6.411.651.33116.01.51
(J17576032)5 RUbr332.98−17.182.910.304.77 
(J18036707)5 RUbk326.85−20.342.131.8575.0 
J1805+0615RUPbwr33.3513.012.132.2893.20.36
J1810+1744RUPbwr44.6416.811.660.4638.40.28
J1811−2405Rbr7.07−2.562.661.3428.01.33
J1816+4510RUbkr72.8324.743.194.3152.20.04
J1823−3021ARr2.79−7.915.44337.62827.00.72
J1824+1014RUPbr39.0810.654.070.553.210.01
(J18240621)7 Rrb24.163.103.230.9110.7 
J1824−2452ARrx7.80−5.583.05161.892243.02.30
J1827−0849GUq22.361.222.241.1038.4<0.800
J1832−0836Rr23.110.262.720.8316.20.92
J1833−3840RUPbrw356.01−13.261.871.77107.0 
(J18353259B)8 Rr1.53−11.381.834.34279.0 
J1843−1113Rr22.05−3.401.850.9660.00.10
(J18521310)2 RUPr21.27−6.154.311.025.01 
J1855−1436RUPbr20.36−7.573.591.099.290.05
(J1857+0943)2 Rbr42.293.065.361.784.575.00
J1858−2216RUPbr13.58−11.392.380.3911.30.06
(J18585422)5 RUbr342.07−22.692.360.4112.4 
J1901−0125GUr32.82−2.902.793.5864.82.50
J1902−5105RUPbr345.65−22.381.740.9068.71.01
J1903−7051RUPbr324.39−26.513.601.048.800.96
J1908+2105RUPbwr53.695.782.561.3832.40.04
J1909−3744Rbr359.73−19.602.951.4021.61.80
J1921+0137RUPbr37.83−5.942.501.8847.60.10
J1921+1929Rbr53.622.452.653.8281.40.20
J1939+2134Rrx57.51−0.291.5610.511097.013.90
J1946+3417Rbr69.294.713.170.323.900.90
J1946−5403RUPbwr343.88−29.582.710.275.330.35
J1959+2048Rbwrx59.20−4.701.611.68159.00.29
J2006+0148RUPbr43.40−15.762.160.3312.80.21
J2017+0603RUPbr48.62−16.032.900.8313.00.18
J2017−1614RUPbwr27.31−26.222.310.247.670.10
(J20294239)1 RUr358.20−35.515.310.942.48 
J2034+3632GUq76.60−2.343.650.171.40 
J2039−3616Rbr6.33−36.523.270.849.470.50
J2039−5617GUbkr341.27−37.152.651.4230.00.58
J2042+0246RUbr48.99−23.024.531.415.980.06
J2043+1711RUPbr61.92−15.312.380.5715.40.12
J2047+1053RUPbwr57.05−19.684.292.1010.4 
J2051−0827Rbwr39.19−30.414.511.275.492.80
J2052+1219RUPbwr59.14−19.991.990.6733.90.44
J2115+5448RUPbwr95.044.112.617.49167.00.46
(J2116+1345)9 RUb64.15−23.822.220.269.57 
J2124−3358Rrx10.92−45.444.932.066.774.50
J2129−0429RUPbkr48.91−36.947.6123.7029.30.00
J2205+6012Rbr103.693.702.411.9855.40.49
J2214+3000RUPbwrx86.86−21.673.121.5019.20.53
J2215+5135RUPbkr99.87−4.162.612.3462.70.16
J2234+0944RUPbwr76.28−40.443.631.9616.61.90
J2241−5236RUPbwrx337.46−54.932.190.8726.01.83
J2256−1024Rbwr59.23−58.292.291.1437.10.73
J2302+4442RUPbr103.40−14.005.191.333.911.40
J2310−0555RUPbr69.70−57.912.610.5011.00.07
J2317+1439Rbr91.36−42.363.450.242.350.60
J2339−0533RUPbkr81.35−62.482.881.4123.2 

Note. In Column 1, pulsars with names in italics, in parentheses, revealed gamma-ray pulsations after the initial sample was defined, and are not analyzed in this catalog. Superscript numbers denote the following references: (1) TRAPUM Collaboration (2023, in preparation), (2) Appendix A, (3) Padmanabh et al. (2023), (4) Cromartie (2020), (5) Clark et al. (2023a), (6) P. Wang et al. (2023, in preparation), (7) Miao et al. (2023), (8) Gautam et al. (2022b), Zhang et al. (2022), and (9) Lewis et al. (2023). Column 2 gives the discovery and detection codes, as in Table 2, with in addition w = black widow and k = redback. Columns 3 and 4 give the Galactic coordinates, with the rotation period P in column 5. The first period time derivative and the spindown luminosity in Columns 6 and 7 are uncorrected for the Shklovskii effect. Corrected values are in Table 8, Section 4.3. Column 9 gives the radio flux density (or upper limit) at 1400 MHz (see Section 4.1). Some MSPs without an S1400 flux measurement suffer intense scintillation, making few-epoch measurements unreliable.

Download table as:  ASCIITypeset images: 1 2 3

3.1. Using Known Rotation Ephemerides

The most straightforward way to search for gamma-ray pulsations is to fold the LAT data using an existing rotation ephemeris that is valid for the full 12+ yr on orbit. To that end, radio and X-ray astronomers shared >1400 rotation ephemerides with the LAT team, ∼40% of the over 3400 known RPPs (mostly from the ATNF Pulsar Catalog v1.69 115 ; Manchester et al. 2005; as well as some unpublished pulsars; see Table 1). Table 4 includes radio telescopes that contributed timing models used in the LAT pulsar searches. These timing models yielded the discovery of ∼110 gamma-ray pulsars. A special effort was made to search for LAT pulsations from the energetic ( erg s−1) subset of the pulsar population, representing about 10% of all pulsars. We have folded >90% of the high- pulsars (see Figure 5), revealing gamma-ray pulsations in about two-thirds of them. Importantly, the ephemerides provided by the radio community sample the entire plane. The >1200 pulsars that were phase-folded without seeing gamma-ray pulsations are shown as black dots in Figure 2. Folding so many pulsars gives a largely unbiased view of known pulsars, which may emit gamma-rays not visible from Earth, discussed by Johnston et al. (2020).

Figure 5. Refer to the following caption and surrounding text.

Figure 5. Left-hand plots: the lines show the fractions of known field pulsars that have been gamma-folded, vs. spindown power (top) and, when a distance estimate exists, the heuristic gamma-ray flux Gh defined in Section 7.1 (bottom). The points show the fractions of folded pulsars for which pulsations have been detected. The vertical dotted line approximates the minimum detectable integral energy flux G100. Right-hand plots: numbers of known and gamma-ray detected field pulsars vs. . (A quarter of ATNF psrcat pulsars have unmeasured .) An earlier version of this figure appears in Laffon et al. (2015).

Standard image High-resolution image

Table 4. Radio Telescopes Searching for New Pulsars in Unidentified LAT Sources, and Timing Discoveries

TelescopeFrequenciesBeam HWHMDecl. Range a Ndisk/N3PC b References
 (MHz)(arcminutes)   
Pulsar Searching
Parkes14007−90° < δ < 33°18/13(1)
Jodrell Bank14005−37° < δ < 90°0/0(2)
Nançay14004 × 22−39° < δ < 90°3/3(3)
Green Bank350, 820, 200018.5, 7.9, 3.1−46° < δ < 90°46/35(4)
Effelsberg14005−31° < δ < 90°1/1(5)
GMRT325, 61085, 44−55° < δ < 90°9/5(6)
Arecibo327, 140014, 3.3 × 3.8−1.3° < δ < 38°16/9(7)
Molonglo8400.8 × 84−90° < δ < 18°0/0(8)
LOFAR c 14010−7° < δ < 90°3/3(9)
MeerKAT d 544–1088, 856–17120.25–0.5, 0.15–0.3−90° < δ < 44°21/2(10)
FAST14003−15° < δ < 65°3/2(11)

Notes.

a The declination ranges for long-duration deep searches are typically 20° narrower than the maximum pointing ranges listed here. b Total numbers of pulsars discovered in LAT-targeted searches and the subset confirmed as gamma-ray pulsars and included in this catalog. Table 6 lists discoveries yet to reveal gamma-ray pulsations. c The value reported here is typical of synthesized beams using the "Superterp" stations (Sanidas et al. 2019). d Beam widths for MeerKAT are those of a coherent tied-array beam (TAB), and are position dependent; hence, we quote a range of values corresponding to elevations above 30°. Several hundred TABs can be recorded simultaneously, so LAT source regions with semimajor axes up to can be covered in single pointings.

References. (1) Weltevrede et al. (2010b), Camilo et al. (2015, 2016), Kerr et al. (2012), Keith et al. (2011); (2) Hobbs et al. (2004c); (3) Guillemot et al. (2012a), Cognard et al. (2011); (4) Ransom et al. (2011), Tabassum et al. (2021), Ray et al. (2020), Bangale et al. (2023), Sanpa-Arsa (2016); (5) Barr et al. (2013); (6) Bhattacharyya et al. (2013, 2021), Roy et al. (2015); (7) Cromartie et al. (2016); (8) UTMOST, Jankowski et al. (2019), Lower et al. (2020); (9) Bassa et al. (2017a);(10) Stappers & Kramer (2016); (11) Wang et al. (2021).

Download table as:  ASCIITypeset image

The availability of so many ephemerides is possible because of the long-term timing campaigns of several radio telescopes, summarized in Smith et al. (2019). Astronomers originally organized support for the LAT as a "Pulsar Timing Consortium" (Smith et al. 2008) and have continued as the Fermi mission proceeds. We aim for P/50 ephemeris accuracy (0.02 in phase), which we have generally achieved: possible gamma-ray pulse width smearing at this level is unlikely to impede pulsation discovery, given that the narrowest known pulse is P/33 (0.03 in phase) for PSR J1959+2048 (Guillemot et al. 2012b). Gamma-ray pulse widths and background rates are such that phase histograms generally require far fewer than 100 bins (Figure 4 uses 50 bins). Of the ∼110 gamma-ray pulsars found by phase-folding using radio ephemerides, 20 are too faint for the 4FGL-DR3 catalog: they rise above the background only when the photons accumulate in a restricted phase range, and thus probe a population beyond that accessible by the source catalog methods. We do not include additional spectral analysis of these sources, so they are absent from the analysis of Section 6.

Fully describing the rotation of young pulsars, which have rotational irregularities known as timing noise, over a timespan of many years often requires many noise-smoothing parameters, which may be covariant with parameters that have a physical interpretation. That is to say, "whitened" solutions used for our analysis and distributed with this paper accurately calculate rotational phase, but should be used with care when studying braking indices, proper motions, or other pulsar properties.

For a few pulsars, X-ray timing was either combined with radio timing to construct an ephemeris (e.g., PSR J2022+3842; Limyansky 2022), or was used exclusively (e.g., PSR J0540−6919 in the Large Magellanic Cloud; Marshall et al. 2016). Section 7.1 further addresses searches for GeV pulsations in X-ray bright pulsars.

Finally, the timing models enable phase alignment between observations at different wavelengths from different observatories. The absolute phase reference is given by the Tempo2 parameters for the arrival time (TZRMJD) and location (TZRSITE) of a reference pulse of a particular frequency (TZRFREQ). The dispersion measure (DM) determines the frequency-dependent delays in the interstellar medium and is necessary to align radio and high-energy data. Carefully aligned pulse profiles provide information about the relative geometry of the different emission regions.

3.2. Pulsar Search Targets

About half of the gamma-ray pulsars were discovered in searches around LAT sources with pulsar-like properties. In addition, many LAT sources are colocated with previously known radio and/or X-ray pulsars (see Table 5). Section 7 discusses the prospects for finding gamma-ray pulsations from the latter, as well as from "spider"-like binary systems found in LAT unidentified sources, likely to yield still more MSP discoveries. This section describes how we select and improve candidate targets for the deep radio searches and the gamma-ray blind searches described in Sections 3.3 and 3.4.

Table 5. Galactic Pulsars Colocated with LAT Sources

PSR4FGL P Dist epsilon = Epeak f = Association
  (ms)(erg s−1)(kpc) (GeV)Δθ/r95  
J0301+35*J0300.4+3450147.063.1 0.121.19 
J0921-5202J0924.1-5202c9.680.7340.40.51 0.66 
J1054-5943 *J1054.0-5938346.913.852.62.311.970.77 
J1146-6610J1147.7-66183.726.091.80.730.561.00 
J1154-6250 *J1155.6-6245c282.010.9841.47.570.600.75SNR G296.8-00.3
J1302-6350 VJ1302.9-634947.76826.02.60.0155.530.32PSR B1259-63
J1306-4035 *J1306.8-40352.204.7 0.100.45sic
J1332-03 *VJ1331.7-03431106.403.5 0.691.06PKS 1328-034
J1430-6623J1431.5-6627785.440.2261.38.440.470.78 
J1439-5501J1440.2-550528.640.2350.73.010.310.35sic
J1455-59J1456.4-5923c176.2058.56.73.87 0.811RXS J145540.4-591320
J1535-5848 *J1534.7-5842307.183.703.02.613.901.08 
J1545-4550J1545.2-45533.5845.32.20.150.561.19 
J1550-5418 *J1550.8-5424c2069.83103.04.00.461.380.89SNR G327.2-00.1
J1604-44 *VJ1604.5-44411389.207.8  0.96PMN J1604-4441
J1616-5017J1616.6-5009491.3815.53.56.890.780.82 
J1632-4818J1631.7-4826c813.6847.65.34.22 1.18 
J1731-1847J1731.7-18502.3477.84.80.570.650.81sic
J1741-34 *J1740.6-3430875.144.6 0.011.08 
J1743-3153J1743.0-3201193.1157.98.84.370.811.05 
J1743-35 *J1743.9-3539569.984.0 0.710.37 
J1801-1417J1801.6-14183.624.291.10.750.510.27sic
J1806+2819J1807.1+282215.080.4311.33.59 0.83sic
J1811-1925J1811.5-192564.716334.05.00.018.940.26sic
J1813-08 *J1812.2-08564.233.3 0.470.68 
J1838-0549J1838.4-0545235.31101.04.11.401.331.02 
J1850-0026J1850.3-0031166.64333.06.70.910.870.57SNR G032.4+00.1
J1852+0158g*J1852.6+0203185.737.6 0.421.13 
J1852-0002g*J1851.8-0007c245.105.6 0.800.49SNR G032.8-00.1
J1855+0455g*J1855.2+0456101.0110.2 0.000.25 
J1858+0310g*J1857.9+0313c372.756.7 0.660.77LQAC 284+003
J1859+0126g*J1900.8+0118957.709.6 0.601.20NVSS J190146+011301
J1904+0603g*J1904.7+06151974.936.1 0.660.72 
J1906+0646g*J1906.2+0631355.525.3 0.650.80SNR G040.5-00.5
J1907+0631J1906.2+0631323.65526.03.40.110.650.66SNR G040.5-00.5
J1908+0811g *J1908.7+0812181.645.6 0.580.36 
J1911+1051J1911.3+1055190.8769.110.15.020.570.43 
J1915+1150J1915.3+1149100.04539.014.03.150.630.25TXS 1913+115
J1917+1121g*J1916.3+1108510.317.1 0.541.19SNR G045.7-00.4
J1928+1725 VJ1929.0+1729289.843.7 0.500.81 
J1929+1731g *VJ1929.0+17293995.409.2 0.500.35 
J1930+1852J1930.5+1853137.0411528.07.00.0110.460.57PWN G54.1+0.3
J1950+2414J1950.6+24164.309.297.317.090.700.25sic
J1957+2516*J1957.3+25173.9617.42.70.611.030.46sic
J2015+0756*J2015.3+07584.332.1  0.91 
J2051+4434g*J2052.3+44371303.1613.5 0.800.42 
J2052+4421g*J2052.3+4437375.3113.5 0.800.77 
J2055+1545J2055.8+15452.1679.23.60.150.570.06sic
J2327+62*J2325.9+6206c266.004.3  0.91NVSS J232543+620829

Note. Colocation means that the pulsar is within f = Δθ/r95 < 1.2 of a 4FGL-DR3 Source, and has efficiency epsilon=Lγ/E˙<10, assuming the pulsar distance (73 matches were overluminous). For association classes "psr" or "msp," we list the pulsar regardless of the efficiency. These pulsars show no gamma-ray pulsations, although many have not been gamma-ray folded, indicated by no value or by * in the PSR column. Δθ is the angular separation between the pulsar and the source, and r95 is the 95% confidence level semimajor axis of the LAT error ellipse. Epeak is the peak of the 4FGL spectral energy distribution (SED) when the PLEC fit was reported in 4FGL-DR3. The last column lists the 4FGL-DR3 association for the LAT source; "sic" means that it matches the pulsar name. "V" indicates Variability_Index > 24.7. Pulsars with a "g" suffix come from the FAST Galactic plane pulsar survey (GPPS; Han et al. 2021).

Download table as:  ASCIITypeset image

3.2.1. Pulsar-like LAT Sources

The LAT source catalog lists "associations" and "identifications" for two-thirds of the sources. "Association" means that the probability that an object known at some other wavelength is responsible for the LAT emission is estimated to be >80%, where the prior for the Bayesian probability is the sky distribution for each population class (Abdollahi et al. 2020). "Identification" requires either detection of gamma-ray pulsations (in the case of pulsars) or an additional match with observations at another wavelength. Examples are simultaneous blazar flares or supernova remnant morphology. Pulsars are the largest class of identified sources, assigned "PSR" or "MSP" classes in the 4FGL catalog (capital letters indicate identification). A spatial match between a known pulsar and a catalog source without gamma-ray pulsations yields association classes "psr" and "msp" (lowercase lettering indicates association without confirmed identification). The 4FGL catalog table has a second association column for low-probability (<20%) associations. Table 5 is complementary to the association process, listing simple colocations without regard for nonpulsar populations.

Additional observables are used to rank candidates in the long list of unidentified sources. Pulsar spectral energy distributions (SEDs), , have sharp cutoffs in the GeV range (Section 6.1) that distinguish them from other categories of gamma-ray sources. The 4FGL catalog includes a variability index for each source, which can help distinguish between pulsars, whose fluxes are stable (see Section 6.6), and the 10× more common blazars, which can exhibit gamma-ray flares. Early successes in selecting pulsar candidates from among the unidentified LAT catalog sources exploited SED curvature versus flux variability correlations (Ackermann et al. 2012b). Subsequent works used machine-learning methods to refine the selections (Lee et al. 2012; Mirabal et al. 2012; Saz Parkinson et al. 2016; Wu et al. 2018; Luo et al. 2020; Finke et al. 2021), or visual inspection and ranking of the LAT source spectra (Camilo et al. 2015). Regardless of the ranking scheme, these lists of pulsar-like unassociated LAT sources have provided a large number of candidate pulsar positions that have been targeted by radio, X-ray, and gamma-ray searches. One recent example, PSR J1653−0158 (Nieder et al. 2020b), was recently found in the brightest remaining pulsar candidate of Saz Parkinson et al. (2016).

In 4FGL-DR3, about 80% of the unidentified sources have error ellipses with semimajor axes between 0fdg05 and 0fdg15. Improved localization can be valuable for the high-frequency radio observations needed to pierce the thick electron column density (i.e., high DM values) in, for example, the direction of the Galactic center. Improved localization also reduces the computing cost of blind period searches, and it can allow a match with optical or X-ray sources necessary to motivate time-consuming observations and analyses (see Section 3.2.2). Radio telescope arrays can coherently combine data from multiple antennas to form tied array beams, which are typically much smaller than the primary beam. These also benefit from improved localization, since computational resources generally limit the number of synthesized beams that can be formed and searched.

3.2.2. Optical, X-Ray, and Radio Studies to Refine Targets

Fewer than 1% of known neutron stars pulse in the optical band, so optical observations of gamma-ray sources might seem unproductive. However, relative to pulsars discovered in radio surveys, the LAT pulsar sample includes many more members of the subclass of interacting binaries—black widows, redbacks, transitional MSPs, and related systems. The companions in these systems can have relatively bright optical magnitudes and (importantly) often exhibit modulation at the orbital period in brightness, color, and radial velocity. Consequently, searches for periodic optical sources within LAT source error regions have identified many candidate gamma-ray pulsars. Some of these have been followed up with radio detections, while others have enabled direct gamma-ray periodicity detections by using the precise position and orbital parameters to reduce the vast parameter space to be searched, as detailed in Section 3.4. Others remain strong candidates that are likely powering the gamma-ray sources, but have yet to be confirmed as pulsars (see Section 7). Some examples of the power of optical studies of LAT source regions include:

  • 1.  
    Optical and X-ray studies strongly indicated that 1FGL J2339.7−0531 hosts a black-widow MSP (Romani & Shaw 2011; Kong et al. 2012). After LAT blind searches found no pulsar, Ray et al. (2020) discovered radio pulsations. PSR J2339−0533 turned out to be a redback with a companion mass outside the parameter space covered by the LAT searches. The resulting ephemeris yielded gamma-ray pulsations as well (Pletsch & Clark 2015).
  • 2.  
    Optical studies by Romani (2012) led to the first discovery of an MSP, PSR J1311−3430, in a blind search of LAT data (Pletsch et al. 2012c). Detection of radio pulsations followed shortly (Ray et al. 2013). Two further binary MSPs have been discovered in this way: the black-widow PSR J1653−0158 (Nieder et al. 2020a) from an optical candidate discovered by Romani et al. (2014) and Kong et al. (2014), and the redback PSR J2039−5617 (Clark et al. 2021) from an optical candidate discovered by Romani (2015) and Salvetti et al. (2015).
  • 3.  
    Four candidate gamma-ray redbacks identified in optical and X-ray searches were recently identified as radio MSPs by TRAPUM (Halpern et al. 2017; Li et al. 2018b; Swihart et al. 2020; Au et al. 2023), listed in Table 6.
  • 4.  
    Optical studies of an unidentified LAT source by Strader et al. (2015) revealed an accretion disk in a binary system with a 5.4 day period. Camilo et al. (2016) subsequently discovered PSR J1417-4402 using the Parkes radio telescope. First thought to be a redback, its orbit is wide, and the wind interactions differ from other spiders. Swihart et al. (2018) argued that it is a possible progenitor of normal field MSP binaries and dubbed it the "huntsman," as the archetype of an emerging class. Unfortunately, it is difficult to time, and a phase-connected ephemeris that would yield gamma-ray pulsations is not anticipated soon.

Table 6. Thirty-nine MSPs Discovered in Radio Searches of LAT Unidentified Sources, Yet to Show Gamma-Ray Pulsations

PSR4FGL P Dist10−33 Lγ DiscoveryReferences
  (ms)(kpc)(erg s−1)  
J0329+50J0330.1+50383.060.30.052GBT 2019Tabassum et al. (2021)
J0506+50J0506.1+50283.391.40.510FAST 2020P. Wang et al. (2023, in preparation)
J0646−54J0646.4−54552.510.40.039PKS 2017as in Camilo et al. (2015)
J0657−4657J0657.4−46583.950.50.085MKT 2022TRAPUM collab. (2023, in preparation)
J0838−2827J0838.7−28273.620.40.147MKT 2021TRAPUM collab. (2023, in preparation)
J0843+67J0843.3+67122.841.61.13GBT 2017Tabassum et al. (2021)
J1008−46X a 2.720.4GMRT 2021J. Roy et al. (2023, in preparation)
J1036−4353J1036.6−43491.680.40.069MKT 2021Clark et al. (2023a)
J1102+02J1102.4+02464.053.72.53GBT 2019Tabassum et al. (2021)
J1103−5403X3.391.7PKS 2009Keith et al. (2011)
J1120−3618X5.551.0GMRT 2011Roy & Bhattacharyya (2013)
J1304+12J1304.4+12034.181.40.442AO 2017Cromartie (2020)
J1346−2610J1345.9−26122.771.20.452MKT 2022TRAPUM collab. (2023, in preparation)
J1356+0230J1356.6+02342.831.80.757MKT 2021TRAPUM collab. (2023, in preparation)
J1417−4402J1417.6−44032.664.418.4PKS 2015"Huntsman," Swihart et al. (2018)
J1551−0658X7.091.3GBT 2010P. Bangale et al. (2023) in preparation
J1624−39J1624.3−39522.962.66.42GBT 2020Tabassum et al. (2021)
J1634+02J1634.7+02352.120.90.094AO 2020S. Tabassum et al. (2023, in preparation)
J1646−2142X5.851.0GMRT 2011Roy & Bhattacharyya (2013)
J1709−0333J1709.4−03283.521.00.902MKT 2021Clark et al. (2023a)
J1716+34J1716.3+34342.112.50.520AO 2021S. Tabassum et al. (2023, in preparation)
J1727−1609J1728.1−16102.453.76.07GBT 2017Tabassum et al. (2021)
J1802−4718J1802.8−47193.671.20.575PKS 2016as in Camilo et al. (2015)
J1803+1358J1803.2+14021.524.68.44AO 2017Cromartie (2020)
J1814+31J1814.4+31322.093.04.87AO 2017Cromartie (2020)
J1823+1208J1823.2+12095.211.70.646MKT 2021TRAPUM collab. (2023, in preparation)
J1823−3543J1823.8−35442.373.73.22MKT 2020Clark et al. (2023a)
J1828+0625X3.631.0GMRT 2011Roy & Bhattacharyya (2013)
J1831−6503J1831.1−65031.851.00.451MKT 2022TRAPUM collab. (2023, in preparation)
J1845+0200J1845.0+01594.311.70.977AO 2016Cromartie (2020)
J1904−11J1904.4−11292.624.74.80GBT 2020Tabassum et al. (2021)
J1906−1754J1906.4−17572.886.818.4MKT 2020Clark et al. (2023a)
J1910−5320J1910.7−53202.331.00.360MKT 2022Au et al. (2023); Dodge et al. (submitted)
J1919+23J1919.1+23544.632.01.36AO 2017Cromartie (2020)
J1947−11J1947.6−11212.243.13.76GBT 2022J. Strader et al. (2023, in preparation)
J2036−02J2036.2−02071.917.710.9GBT 2020Tabassum et al. (2021)
J2045−6837J2045.9−68352.961.30.421PKS 2017as in Camilo et al. (2015)
J2051+50 1.681.4GBT 2020Tabassum et al. (2021)
J2333−5526J2333.1−55272.102.53.13MKT 2021TRAPUM collab. (2023, in preparation)

Notes. Six pulsars with "X" in the 4FGL column are unrelated to the gamma-ray source that guided the radio search. Localization of the two pulsars with no 4FGL name is as yet inadequate for source matching. Distances use DM and the YMW16 model (Yao et al. 2017). Three unpublished pulsars were found using methods as in Camilo et al. (2015). The Discovery column takes "Fermi" entries from http://astro.phys.wvu.edu/GalacticMSPs/GalacticMSPs.txt. See also http://www.trapum.org/discoveries/.

a Search targeted uncatalogued LAT source 605N-0275 (T. Burnett 2023, private communication).

Download table as:  ASCIITypeset image

Neutron stars can also power X-ray emission either directly from their surfaces or magnetospheres, or via shock emission when the pulsar wind interacts with a companion or the interstellar medium (Gentile et al. 2014; Marelli et al. 2015; Coti Zelati et al. 2020). The nonthermal shock spectra, sometimes modulated at the orbital period, provide strong signatures. Consequently, X-ray surveys also provide valuable information as demonstrated by varied studies of the error ellipses of unidentified LAT sources. For example, the Swift XRT systematically studied a large number of LAT unassociated source regions. 116 This data set has been mined using machine-learning techniques to identify promising pulsar candidates for radio follow-up (Kaur et al. 2019; Kerby et al. 2021).

Pulsars are steep-spectrum radio point sources (Jankowski et al. 2018). This was used to identify the source (4C 21.53) that was later found to be the first millisecond pulsar (PSR B1937+21; Backer et al. 1982). The availability of radio surveys covering nearly the whole sky was recently harnessed to look for pulsar candidates associated with LAT source regions (Frail et al. 2016, 2018; Bruzewski et al. 2021; Ray et al. 2022). These candidate lists have high-precision positions from radio interferometric imaging and well-measured fluxes, enabling both gamma-ray blind searches with reduced parameter spaces and radio pulsation searches that discovered several new MSPs.

3.3. Targeted Radio Searches

Shortly after LAT operations began, the Fermi Pulsar Search Consortium (PSC) was organized, bringing together the LAT team and many international collaborators with pulsar searching expertise and access to the world's largest radio telescopes (Ray et al. 2012). The aim was to coordinate radio pulsation searches of LAT sources, and to see which pulsars found in gamma-ray blind searches are truly radio-quiet. The PSC has been exceptionally successful at efficiently discovering MSPs in searches of LAT sources. Beginning in 2010, the TRAPUM collaboration formed to exploit MeerKAT (Stappers & Kramer 2016). It is not formally part of the PSC but includes searches of LAT sources. The FAST radio telescope in China also joined the searches under a separate arrangement. Radio astronomers are provided with early access to LAT source lists as well as the parameters of new pulsars discovered in gamma-ray blind searches (see Section 3.4).

LAT gamma-ray source localizations are much better than was possible with EGRET, making it often possible to search an entire 95% confidence error region with a single radio telescope pointing. Of course, the primary beam of a radio telescope depends on the observation frequency and the dish diameter. Thus, radio astronomers could choose telescope/frequency combinations optimized for the error regions being searched. Initially, the telescopes were all single-dish (Green Bank Telescope, GBT, Parkes, Effelsberg, Nançay Radio Telescope, NRT, Arecibo Telescope, and the Lovell Telescope at Jodrell Bank), but array telescopes were soon added, the first being the Giant Metrewave Radio Telescope (GMRT), followed by LOFAR and MeerKAT (see Table 4).

Array-based telescopes offer increasingly powerful and flexible methods for pulsar searches. Incoherently combining the beams allows for very efficient but less sensitive searches of large error boxes. Some modern arrays have backends powerful enough to search hundreds of coherent (tied-array) beams arranged to cover a large primary beam, enabling extremely efficient full sensitivity searches of large error boxes (Clark et al. 2023a). Gated imaging (Roy & Bhattacharyya 2013) or tied-array beamforming (Stappers & Kramer 2016; Sanidas et al. 2019) can localize newly discovered pulsars extremely well, without the need for long-term timing campaigns.

A key feature of radio pulsar searches, compared to LAT blind searches, is that it is computationally feasible to maintain sensitivity to short-period binary pulsars even when none of the system parameters are known (as described in Section 3.4, binary MSPs can be found in gamma-ray blind searches when approximate orbital parameters are known from optical studies). Because the observation durations are generally short compared to the orbital period, this can be done through Fourier search techniques that assume a constant acceleration (Ransom et al. 2002), or a constant jerk (Andersen & Ransom 2018). Recently, an analysis of PSC search data using the jerk search method resulted in the discovery of 11 new MSPs, at least one of which would likely have been missed in a standard acceleration search (Tabassum et al. 2021).

At least 110 new pulsars have been discovered in the radio searches targeting LAT unassociated sources. All except eight are MSPs—MSPs dominate the high Galactic latitudes where radio surveys have been less thorough. At present, only one of the slow pulsars has been confirmed as a gamma-ray pulsar (PSR J2030+3641; Camilo et al. 2012). Given the many observations, some chance coincidences are expected, especially for telescopes with wide fields of view (particularly arrays of smaller telescopes combined incoherently, like the GMRT; e.g., Bhattacharyya et al. 2022). Table 6 lists 39 MSPs found in the targeted radio searches for which gamma-ray pulsations have not yet been detected. An "X" indicates six, where it has been concluded that the MSP is not associated with the target source. Most of the others were discovered after 2016, so the lack of gamma-ray pulsations from these sources is likely due to the lag in obtaining a rotation ephemeris with which to fold LAT data.

Developing such an ephemeris is an involved process, as the timing model must extrapolate well over the long span of LAT data required for the accumulated signal to become significant. Especially if the position is not well constrained by radio imaging or an X-ray or optical counterpart, then the radio campaign requires a dense set of observations to measure the spin period well and longer intervals (0.5–1.0 yr) to measure the position and . Pulsars with orbital period variability may require years of monitoring. Such follow-up campaigns can require multiple proposal cycles.

In most cases, after getting a good timing solution, LAT pulsations have been detected from the MSPs. If pulsations are detected in a fraction of the LAT data, the solution can often be extended to cover the full mission using LAT timing methods (see Section 3.5). We expect that a large fraction of the recent MSPs in Table 6 will be identified as gamma-ray pulsars in the future. Some of the pulsars in Table 5 may also someday reveal gamma-ray pulsations.

3.4. Blind Periodicity Searches

In addition to the detection of gamma-ray pulsations from known radio or X-ray pulsars via folding, the LAT data also allow for the discovery of entirely new pulsars through direct searches for their gamma-ray pulsations, without prior knowledge of their spin periods or spin-down rates. These searches complement and contrast the efforts described above in Sections 3.1 and 3.3, providing access to a population of otherwise undetectable pulsars. However, these searches are less sensitive than following-up a known radio or X-ray pulsar, and are extremely computationally demanding: searching for isolated young pulsars/MSPs in a single Fermi-LAT source typically takes a few days/weeks, respectively, using tens of thousands of CPU cores, while a search for a single binary MSP can run for months on thousands of GPUs.

Of the pulsars in this catalog, 81 were first discovered in the LAT data, compared to 36 in 2PC. The majority are young pulsars, located near the Galactic plane. Notably, this includes 10 LAT-discovered MSPs of which three are in binary systems, whereas 2PC contained just one such object. That the discovery rate from direct gamma-ray pulsation searches has not fallen off is perhaps surprising when one considers how the computational cost and sensitivity of these searches increases with the data duration.

The computational cost depends on the dimensionality and volume of the parameter space that must be searched, and the required density of search points covering this volume. At minimum, these searches are two-dimensional, covering the range of spin frequencies (or periods) and spin-down rates spanned by the known pulsar population. A subarcsecond precision position—much better than can be provided by the LAT—is required to correct for the Earth's orbital Römer delay. Unless one is available from multiwavelength observations, an additional two-dimensional grid of sky locations covering the localization region must be searched. The required positional precision in each dimension is proportional to the spin frequency, f. This would make searches for MSPs thousands of times more costly than searching for young pulsars, were this not partially mitigated by the far smaller range of spin-down rate over which MSPs are found.

When searching for pulsars in circular binaries, at least three further parameters are required (the orbital period, phase, and radius), with two additional parameters required for eccentric systems. The search grid density in all binary parameters depends on f, again making rapidly spinning MSPs substantially more costly to find. The binary search parameter space is in fact so large that entirely uninformed searches are as yet infeasible, with searches only becoming possible when a candidate binary system is identified and its orbital parameters constrained via multiwavelength observations.

The required search-grid density in most dimensions depends on the duration Tobs of the data over which the signal is integrated (Pletsch & Clark 2014). In f, and , the number of grid points scales linearly and quadratically, respectively, with Tobs. For Tobs < 1 yr, the number of sky positions that must be searched scales quadratically with Tobs, but this saturates at longer timescales. For binary pulsar searches, the required density of orbital period search points also scales linearly with Tobs (Nieder et al. 2020a). The result is that the computing cost scales with at least for isolated pulsars and for binaries. This rapid scaling leads to fully coherent searches of data lasting more than a few months being entirely infeasible.

A compromise is achieved using "semicoherent" search methods, where the signal powers from many shorter windows (of duration Tcoh, typically a few weeks long) are combined incoherently (Atwood et al. 2006; Pletsch et al. 2012a). For a semicoherent search, the numbers of grid points in the frequency and spindown dimensions are each reduced by a factor of Tobs/Tcoh, while the number of sky positions that must be searched is reduced by a factor of (assuming Tcoh < 1 yr < Tobs). This semicoherent search is performed as an initial scan over the entire parameter space, with more sensitive fully coherent stages being used to "follow-up" a smaller number of semicoherent candidates to increase their significance.

The compromise is that the sensitivity of a semicoherent search is far worse than that of a fully coherent search (e.g., for pulsations from a known pulsar). When folding with a known ephemeris, the S/N is proportional to . In a semicoherent search, this scaling becomes . Combining this with the computing cost scaling, we see that the recovered S/N (a proxy for sensitivity) increases very slowly as the computing cost, C(Tobs, Tcoh), increases. For fixed Tobs > 1 yr, , where a = 1/8 if positional information is available, or a = 1/16 otherwise. For fixed Tcoh < 1 yr, where b = 1/4 for isolated pulsars, and b = 1/8 for binary pulsars.

Two additional effects further limit this sensitivity. First, the vast parameter space that must be searched over introduces a large "trials" factor, meaning that pulsed signals need to be detected with far higher test statistic values in order to be statistically significant. Where pulsars can be found in single-trial ephemeris folding searches at H = 25, the associated trials factor of a semicoherent search raises this significance threshold. To be statistically significant, these searches require the coherent power at the fundamental frequency to be P1 ≳ 120. For typical gamma-ray pulse profile shapes, this corresponds to H ≳ 200. Tcoh can be chosen such that signals above this level will be detectable in the initial semicoherent stage.

Second, the accumulation of a pulsed signal can be disrupted by more complex timing behavior such as timing noise, glitching, or orbital period variations in the case of redback binary systems. However, statistical noise still accumulates over the full data duration. For timing behavior that results in a signal being "well behaved" only over a timescale τ, this means that the S/N is further reduced by a factor of . This latter effect does not occur in isolated MSPs, whose spindown rates are extremely stable over time, whereas energetic young pulsars, which tend to exhibit the most timing noise and glitches, will be the worst affected. Clark et al. (2017) predicted that the detection rates of isolated MSPs would likely increase relative to those of young pulsars as a result.

In summary, as the LAT data duration increases, these searches become more expensive, and the minimum detectable flux density decreases more slowly than that of an ephemeris-folding search. The ratio between these sensitivities depends on the background photon flux, and hence on the location of the targeted source in the sky, but varies between a factor of ∼3 and 4 for the ongoing Einstein@Home pulsar surveys (with Tobs = 10 yr, Tcoh = 222 s). The difference in sensitivities between these two methods is emphasized by the difference in the gamma-ray fluxes of the radio-loud and radio-quiet populations shown in Figure 6.

Figure 6. Refer to the following caption and surrounding text.

Figure 6. The radio flux density at 1400 MHz vs. the integral energy flux above 100 MeV. Gamma-ray and radio fluxes are essentially uncorrelated. Dark borders on the triangle symbols indicate MSPs discovered in deep radio searches of previously unidentified LAT sources by the PSC (see Section 3.3). The horizontal line at 30 μJy shows our convention to distinguish radio "loud" from radio "quiet" pulsars.

Standard image High-resolution image

To tackle these growing challenges, search efforts require increasing computing resources and methodological developments to increase efficiency and sensitivity. Fortunately, the available computing power does increase over time as technology progresses, and as search efforts have shifted to larger computing systems. The first successful searches in the early Fermi mission either searched the center of the localization region for each source, or targeted the location of a plausible X-ray counterpart (Abdo et al. 2009b; Saz Parkinson et al. 2010). After the first 2 yr of the Fermi mission, these discoveries became less frequent. Newer searches (Pletsch et al. 2012a, 2012b) mitigated dependence on event selection criteria and source localization by weighting events and shifting to large computing clusters to enable searching over a grid of positions. These searches were later migrated to the Einstein@Home distributed volunteer computing system, allowing for longer Tcoh and hence more sensitivity (Pletsch et al. 2013; Clark et al. 2015, 2016, 2017, 2018).

Additional sensitivity can be gained not just through computing power, but also through careful study. While semicoherent search methods were developed prior to Fermi's launch (Atwood et al. 2006), improvements to their sensitivity and efficiency have been made throughout the LAT mission (Pletsch & Clark 2014; Nieder et al. 2020a). Improvements to the LAT event reconstruction and background rejection increase the recoverable signal strengths (Atwood et al. 2013; Bruel et al. 2018), as do higher-level processing developments such as photon weighting (Kerr 2011; Bruel 2019). An unbinned likelihood method improved LAT's sensitivity at all time scales, useful for a broad range of astrophysical targets in addition to pulsars (Kerr 2019).

Extensive multiwavelength searches for possible counterparts at which to target pulsation searches continue to lead to new discoveries. To give some recent examples: targeted searches of X-ray sources discovered energetic gamma-ray pulsars within two supernova remnants (PSRs J1111−6039 and J1714−3830), and searches targeting steep-spectrum radio sources associated with LAT sources discovered two isolated MSPs (Frail et al. 2018).

Searches for periodic optical and X-ray counterparts to pulsar-like LAT sources have identified several promising black-widow and redback candidates, detailed in Section 7. The precise orbital constraints from long-term optical observations are key inputs for binary searches, although they remain the most computationally demanding projects. Nevertheless, the development of efficient search grid designs (Nieder et al. 2020a) and a GPU-accelerated 117 search algorithm deployed on Einstein@Home has resulted in two binary MSP discoveries (Nieder et al. 2020b; Clark et al. 2021), in addition to PSR J1311−3430 (Pletsch et al. 2012a), which was included in 2PC.

The complementary nature of the blind-search and ephemeris-folding populations is highlighted by the number of radio detections of blind-search pulsars. In 2PC, there were radio detections of just four pulsars discovered in gamma-ray searches. This fraction remains small in this catalog: only 11 radio detections for 81 gamma-ray blind-search pulsars, despite deep radio follow-ups. Of those detected, most have radio flux densities lying below typical radio survey detection thresholds, which are above the 30 μJy limit we use to define "radio-quiet," shown in Figures 6 and 7. This is perhaps unsurprising: most pulsars discovered in gamma-ray searches are young pulsars located in the Galactic plane, which has been extensively surveyed by radio telescopes, and therefore we would expect most radio-loud pulsars in this region to be known. Most gamma-ray pulsars found off the Galactic plane are MSPs, of which most are found in binary systems, and therefore inaccessible to gamma-ray pulsation searches unless orbital constraints from an optical counterpart can be derived. A notable addition since 2PC is the six MSPs that remain undetected in radio observations, including the black-widow PSR J1653-0158 (Nieder et al. 2020b). In such spider binaries, the dense ionized wind can obscure the pulsed radio signal for much, or at times all, of the orbit, hindering radio detection.

Figure 7. Refer to the following caption and surrounding text.

Figure 7. Radio flux density at 1400 MHz vs. pulsar distance. The horizontal line at 30 μJy shows our convention to distinguish radio "loud" from radio "quiet" pulsars. The diagonal line shows a threshold in pseudo-luminosity of 100 μJy-kpc2. The pulsars at the lower-right are assigned distance limits along the Milky Way's rim in Figure 8. DMs for the pulsars on the line at 25 kpc saturate the YMW16 model.

Standard image High-resolution image

3.5. Timing Pulsars Using LAT Data

The timing validity range of a recently discovered radio or X-ray pulsar will likely not cover the entire LAT data set, limiting the statistics available for profile characterization, and limiting the discovery potential for glitches, profile changes, spin-ups (Marshall et al. 2015), or other rare behavior. Often, pulsar timing using LAT data allows for an accurate phase-connected ephemeris covering the entire mission duration.

LAT pulsar timing has significantly improved, with new techniques to handle timing noise and faint pulsars developed by Ray et al. (2011) and Kerr et al. (2015). They use pulse ToA estimation and time-domain models of timing noise. For faint pulsars, ToA estimates for short data segments are not possible, and we adopted the unbinned methods used in Ajello et al. (2022). They optimize the timing model by maximizing the likelihood, using PINT (Luo et al. 2021) to compute the rotational phase for each photon. Parameter uncertainties are derived from the likelihood curvature near the optimum. Particularly useful for MSPs, we have been able to estimate some precise proper motions directly from LAT data, necessary for the Doppler corrections presented in Section 4.3.

When both radio and LAT timing of comparable quality are available, we use the radio ephemeris. LAT timing was used for 142 of the pulsars in this catalog. Two-thirds of those were made with the ToA method. The fraction of timing solutions using the unbinned method is increasing, as analysts adopt the tool, as more difficult pulsars are added to the sample, and as the data sets lengthen. The ephemerides we used to fold the 294 pulsars are included in the supplementary data, at https://fermi.gsfc.nasa.gov/ssc/data/access/lat/3rd_PSR_catalog/, and at https://fermi.gsfc.nasa.gov/ssc/data/access/lat/ephems/. See Appendix D.

4. The Gamma-Ray Pulsars

The search methods yielded gamma-ray pulsations for 294 pulsars. Table 1 breaks down the numbers by different types. Here we compile their radio fluxes, distances, and, for the MSPs, their proper motions, used to Doppler-correct their spindown properties (Shklovskii 1970). Figure 3 shows the sky distribution of the Fermi pulsars. The LAT's uniform sky coverage is important for detecting MSPs, as these relatively nearby sources are approximately uniformly distributed in Galactic latitude.

4.1. Radio Intensities

We tabulate measured or derived 1400 MHz flux density values (S1400) in Tables 2 and 3. ATNF psrcat provides flux densities Sν from among 27 different radio frequencies ν; thus, as a starting point, we used psrqpy (Pitkin 2018) to extract all Sν values from psrcat, as well as the radio spectral index, SPINDX. When available, we use S1400. Otherwise, we calculate with α = SPINDX when present, and α = − 1.7 otherwise. We use the Sν value corresponding to the frequency ν closest to 1400 MHz. We used α = − 1.7 in 2PC, a compromise among published values: −1.6 from Lorimer et al. (1995), −1.8 ± 0.2 from Maron et al. (2000), and −1.60 ± 0.03 from Jankowski et al. (2018). For MSPs, Spiewak et al. (2022) found −1.92 ± 0.06 with a standard deviation of 0.6, and Dai et al. (2015) found −1.81 ± 0.01. The small formal uncertainties are potentially misleading: for example, the value from Jankowski et al. (2018) is for a sample of 400 pulsars with a log-normal distribution of standard deviation 0.54, the distribution being sensitive to, e.g., the frequency range used, yielding the spread in average indices. Using multiple Sν values to fit a power law gives S1400 values generally compatible with those found above, and possibly more accurate in some cases. However, discrepant Sν values lead to unrealistic results in some cases; hence, we use the simple prescription described above.

The psrcat Sν values have accumulated over decades, using a variety of instruments and analysis conventions. Frequency extrapolation introduces additional biases. Fortunately, v1.68 psrcat includes S1400 values from two recent studies. Spiewak et al. (2022) used MeerKAT to obtain S1400 for 189 MSPs, 118 and Parent et al. (2022) analyzed Arecibo data from the PALFA survey for 93 MSPs. We also replace the Hessels et al. (2011) values with the updated results reported by Bangale et al. (2023). For a small number of pulsars, we obtained as-yet-unpublished Sν measurements, noted in the Tables.

In addition to pulsars that are well characterized from radio surveys, those discovered in blind searches are typically searched deeply for radio pulsations (Camilo et al. 2009; Abdo et al. 2010d; Saz Parkinson et al. 2010; Ray et al. 2011, 2012; Pletsch et al. 2012a; Nieder et al. 2020b). The 81 blind discoveries were followed by 11 radio detections, six of which qualify as "radio-loud," with S1400 > 30 μJy (see below). Radio measurements are not yet available for recently discovered LAT pulsars.

Figure 6 shows that radio flux densities are essentially uncorrelated with the integrated gamma-ray energy fluxes, G100. Population studies aim to constrain a power-law model of radio luminosity in which . Johnston & Karastergiou (2017) favored s = 1/4 but recognized that a range of s values match observations, because of the large dispersion. Posselt et al. (2023) found an even smaller index. Since (see Figure 23), flux correlation between the two wavelengths is naturally expected. However, the fraction of neutron star spindown power carried off by radio waves is thousands of times lower, or more, than that of gamma-rays, and radio emission is so weakly connected to the braking mechanism that the correlation is invisible.

Figure 7 shows how the radio flux densities of the gamma-ray pulsars compare with the overall pulsar population, as well as their lack of distance dependence. As in 2PC, we define a pulsar as "radio-loud" if S1400 > 30 μJy, and "radio-quiet" if the measured flux density is lower. The horizontal lines in Figures 6 and 7 show the threshold. Also as in 2PC, the diagonal line in Figure 7 shows an alternate threshold at pseudo-luminosity 100 μJy-kpc2.

Whether a radio or a gamma-ray beam (or both) sweeps through our line of sight from Earth depends on the beam shapes and orientations. Sieber et al. (1975) already showed that radio beam size often increases with decreasing frequency. Grießmeier et al. (2021) searched 27 radio-quiet northern gamma-ray pulsars at 150 MHz but obtained no new detections. This result is compatible with the current understanding of these pulsars' orientations based on their emission geometries, which are derived from LAT light curves and plotted in Johnston et al. (2020), Grießmeier et al. (2021), and which are consistent with the radio beam not crossing the Earth's line of sight.

4.2. Distances

Converting integrated energy fluxes G100 to emitted luminosities Lγ = 4π d2 fΩ G100 requires the source distances d. (We set the beaming fraction fΩ = 1; see Section 6.7.) Knowing pulsar distances also allows us to map neutron star distributions relative to the Galaxy's spiral arms, as in Figure 8, and evaluate their scale height above the plane. The various methods to estimate pulsar distances offer widely differing accuracies. Tables 7 and 8 list the distances that we use, the estimation methods, and the references.

Figure 8. Refer to the following caption and surrounding text.

Figure 8. Gamma-ray pulsar positions projected onto the Milky Way model of Reid et al. (2009). The MSPs that appear coincident with the Galactic center, PSRs J1823−3021A and J1835−3259B in the globular clusters NGC 6624 and 6652, lie well above the Galactic plane. Distance uncertainties are not shown, but can be large for the more distant objects. The squares with arrows indicate the lines of sight toward pulsars for which no distance estimates exist, placed at the distances where 95% of the electron column density has been integrated in the YMW16 model. Other symbols are the same as in Figure 2.

Standard image High-resolution image

Table 7. Distance Estimates for Young LAT-detected Pulsars

Pulsar NameDistance (kpc)MethodReference
J0002+62162.0 ± 0.4OSchinzel et al. (2019)
J0007+73031.4 ± 0.3KPineault et al. (1993)
J0106+48553.1 ± 1.2DMPletsch et al. (2012a)
J0117+59141.8 ± 0.7DMBilous et al. (2016)
J0139+5814 PChatterjee et al. (2009)
J0205+64493.2 ± 0.2KXu et al. (2006)
J0248+60212.0 ± 0.20Theureau et al. (2011)
J0357+3205<7.1DMMYao et al. (2017)
J0359+5414<9.2DMMYao et al. (2017)
J0514−44081.0 ± 0.4DMBhattacharyya et al. (2019)
J0534+2200 PBrown et al. (2021)
J0540−691949.7 ± 1.1OWalker (2012)
J0554+3107<9.0DMMYao et al. (2017)
J0622+3749<8.0DMMYao et al. (2017)
J0631+06464.6 ± 1.8DMWu et al. (2018)
J0631+10362.1 ± 0.8DMWeltevrede et al. (2010a)
J0633+0632<9.7DMMYao et al. (2017)
J0633+17460.2 ± 0.1PCaraveo et al. (1996)
J0659+1414 PBrisken et al. (2003)
J0729−14482.7 ± 1.1DMMcEwen et al. (2020)
J0729−1836 PDeller et al. (2019)
J0734−1559<9.7DMMYao et al. (2017)
J0742−2822 KPLKVerbiest et al. (2012)
J0744−2525<10.6DMMYao et al. (2017)
J0802−5613<11.0DMMYao et al. (2017)
J0834−41595.5 ± 2.2DMKramer et al. (2003)
J0835−4510 PDodson et al. (2003)
J0908−4913 KPLKVerbiest et al. (2012)
J0922+0638 PChatterjee et al. (2001)
J0940−54280.4 ± 0.2DMPetroff et al. (2013)
J1016−58573.2 ± 1.3DMPetroff et al. (2013)
J1019−574910.9 ± 4.4DMPetroff et al. (2013)
J1023−5746<11.9DMMYao et al. (2017)
J1028−58191.4 ± 0.6DMKeith et al. (2008)
J1044−5737<11.2DMMYao et al. (2017)
J1048−5832 KPLKVerbiest et al. (2012)
J1055−60283.8 ± 1.5DMPetroff et al. (2013)
J1057−52260.3 ± 0.2OMignani et al. (2010)
J1057−5851<11.6DMMYao et al. (2017)
J1105−6037<12.7DMMYao et al. (2017)
J1105−61072.4 ± 0.9DMPetroff et al. (2013)
J1111−6039<12.8DMMYao et al. (2017)
J1112−61034.5 ± 1.8DMOswald et al. (2021)
J1119−61278.4 ± 0.4KCaswell et al. (2004)
J1124−5916 XGonzalez & Safi-Harb (2003)
J1135−6055<13.5DMMYao et al. (2017)
J1139−6247<14.4DMMYao et al. (2017)
J1151−61082.2 ± 0.9DMSerylak et al. (2021)
J1203−6242<15.0DMMYao et al. (2017)
J1208−6238<15.1DMMYao et al. (2017)
J1224−64074.0 ± 2.0KPLKVerbiest et al. (2012)
J1231−5113<10.2DMMYao et al. (2017)
J1231−6511<17.0DMMYao et al. (2017)
J1253−58201.6 ± 0.7DMD'Amico et al. (1998)
J1341−622012.6 ± 5.0DMPetroff et al. (2013)
J1350−6225<17.4DMMYao et al. (2017)
J1357−64293.1 ± 1.2DMLorimer et al. (2006)
J1358−6025<16.6DMMYao et al. (2017)
J1410−613213.5 ± 5.4DMPetroff et al. (2013)
J1413−6205<18.0DMMYao et al. (2017)
J1418−60581.6 ± 0.7OYadigaroglu & Romani (1997)
J1420−60485.6 ± 2.3DMPetroff et al. (2013)
J1422−6138<18.2DMMYao et al. (2017)
J1429−5911<17.4DMMYao et al. (2017)
J1447−5757<18.2DMMYao et al. (2017)
J1459−6053<19.7DMMYao et al. (2017)
J1509−58503.4 ± 1.3DMPetroff et al. (2013)
J1513−5908 KPLKVerbiest et al. (2012)
J1522−5735<19.3DMMYao et al. (2017)
J1528−5838<20.5DMMYao et al. (2017)
J1531−56102.8 ± 1.1DMPetroff et al. (2013)
J1614−50485.1 ± 2.1DMPetroff et al. (2013)
J1615−5137<18.6DMMYao et al. (2017)
J1620−4927<18.4DMMYao et al. (2017)
J1623−5005<18.4DMMYao et al. (2017)
J1624−4041<18.1DMMYao et al. (2017)
J1641−5317<22.3DMMYao et al. (2017)
J1646−43466.2 ± 2.5DMHobbs et al. (2004a)
J1648−46114.5 ± 1.8DMPetroff et al. (2013)
J1650−4601<22.6DMMYao et al. (2017)
J1702−41284.0 ± 1.6DMOswald et al. (2021)
J1705−19060.7 ± 0.3DMHobbs et al. (2004b)
J1709−4429 KPLKVerbiest et al. (2012)
J1714−3830<20.5DMMYao et al. (2017)
J1718−38253.5 ± 1.4DMPetroff et al. (2013)
J1730−33503.5 ± 1.4DMPetroff et al. (2013)
J1731−4744 OShternin et al. (2019)
J1732−31310.6 ± 0.3DMMaan et al. (2012, 2017)
J1736−3422<21.8DMMYao et al. (2017)
J1739−30233.1 ± 1.2DMPetroff et al. (2013)
J1740+10001.2 ± 0.5DMBilous et al. (2016)
J1741−2054 OBrownsberger & Romani (2014)
J1742−3321<22.1DMMYao et al. (2017)
J1746−3239<22.2DMMYao et al. (2017)
J1747−29584.8 ± 0.8XGaensler et al. (2004)
J1748−2815<20.0DMMYao et al. (2017)
J1757−24213.1 ± 1.2DMHobbs et al. (2004b)
J1801−24513.8 ± 1.5DMPetroff et al. (2013)
J1803−2149<21.8DMMYao et al. (2017)
J1809−23321.7 ± 1.0KOka et al. (1999)
J1813−1246<22.5DMMYao et al. (2017)
J1816−07553.2 ± 1.3DMMcEwen et al. (2020)
J1817−1742<22.0DMMYao et al. (2017)
J1826−1256<19.0DMMYao et al. (2017)
J1827−1446<22.0DMMYao et al. (2017)
J1828−11014.8 ± 1.9DMPetroff et al. (2013)
J1831−09523.7 ± 1.5DMLorimer et al. (2006)
J1833−10344.1 ± 0.3KPLKVerbiest et al. (2012)
J1835−11063.2 ± 1.3DMD'Amico et al. (1998)
J1836+59250.5 ± 0.3XHalpern et al. (2002)
J1837−06044.8 ± 1.9DMOswald et al. (2021)
J1838−0537<18.2DMMYao et al. (2017)
J1841−05244.1 ± 1.7DMPetroff et al. (2013)
J1844−0346<17.4DMMYao et al. (2017)
J1846+0919<21.7DMMYao et al. (2017)
J1846−0258 KPLKVerbiest et al. (2012)
J1853−00045.3 ± 2.1DMPetroff et al. (2013)
J1856+01133.3 ± 0.6KCaswell et al. (1975)
J1857+01434.6 ± 1.8DMSerylak et al. (2021)
J1906+0722<18.9DMMYao et al. (2017)
J1907+06022.6 ± 1.0DMAbdo et al. (2010d)
J1913+09043.0 ± 1.2DMLorimer et al. (2006)
J1913+10114.6 ± 1.8DMMorris et al. (2002)
J1925+17205.0 ± 2.0DMSerylak et al. (2021)
J1928+17464.3 ± 1.7DMNice et al. (2013)
J1932+1916<17.7DMMYao et al. (2017)
J1932+2220 KPLKVerbiest et al. (2012)
J1935+20254.6 ± 1.8DMLorimer et al. (2013)
J1952+32523.0 ± 2.0KPLKVerbiest et al. (2012)
J1954+2836<16.4DMMYao et al. (2017)
J1954+38524.7 ± 1.9DMMcEwen et al. (2020)
J1957+5033<14.7DMMYao et al. (2017)
J1958+2846<16.1DMMYao et al. (2017)
J2006+31026.0 ± 2.4DMNice et al. (2013)
J2017+3625<15.0DMMYao et al. (2017)
J2021+3651 atnfAKirichenko et al. (2015)
J2021+40261.5 ± 0.4KLandecker et al. (1980)
J2022+3842 OArzoumanian et al. (2011)
J2028+3332<14.5DMMYao et al. (2017)
J2030+36413.0 ± 1.0OCamilo et al. (2012)
J2030+4415 Ode Vries & Romani (2022)
J2032+41271.8 ± 0.1PBrown et al. (2021)
J2043+27401.5 ± 0.6DMBilous et al. (2016)
J2055+2539<10.8DMMYao et al. (2017)
J2111+4606<12.8DMMYao et al. (2017)
J2139+4716<12.0DMMYao et al. (2017)
J2208+40560.8 ± 0.3DMLynch et al. (2018)
J2229+61143.0 ± 1.2XHalpern et al. (2001)
J2238+5903<11.3DMMYao et al. (2017)
J2240+58327.3 ± 2.9DMTheureau et al. (2011)

Note. The best known distances of the 150 young pulsars detected by the Fermi LAT. For the DM method, the reference is for the dispersion measurement, and the distance was obtained with the YMW16 model for the Galactic electron density (Yao et al. 2017). "DMM" is the distance at which the YMW16 DM reaches 95% of its maximum value for that line of sight. Method "P" is parallax, "K" is kinetic (H i absorption), "KPLK" combines the two and corrects for Lutz–Kelker bias (Verbiest et al. 2012), "X" indicates an X-ray measurement of the hydrogen column density, and "O" englobes optical and other methods.

Download table as:  ASCIITypeset images: 1 2 3

Table 8. Millisecond Pulsar Distances and Spindown Doppler Corrections

PSR d MethodReference μ Reference (1 − ξ)
 (pc)  (mas yr−1) (10−21)(10−21)(1033 erg s−1) 
J0023+0923 PAlam et al. (2021)13.9 ± 0.1Same2.57−0.59 0.82
J0030+0451320 ± 10PAlam et al. (2021)6.4 ± 0.1Same0.15−0.523.602 ± 0.0031.0
J0034−05341340 ± 530DMAbdo et al. (2010e)12.1 ± 0.5Desvignes et al. (2016)0.90−0.4026.6 ± 2.20.90
J0101−64221000 ± 400DMSpiewak et al. (2022)17.3 ± 0.5LAT1.87−0.518.0 ± 1.80.72
J0102+48392310 ± 920DMHessels et al. (2011)4.1 ± 0.3NRT0.28−0.1417.6 ± 0.20.99
J0154+18331620 ± 640DMMartinez et al. (2019)13.6 ± 1.5Same1.72−0.173.8 ± 2.40.45
J0218+4232 PLKVerbiest & Lorimer (2014)6.5 ± 0.1Du et al. (2014)0.760.11 0.99
J0248+42301850 ± 740DMBhattacharyya et al. (2021)(6.0 ± 3.3 )Same<1.80.09>33>0.89
J0251+26061170 ± 460DMDeneva et al. (2021)17.8 ± 3.6LAT2.30−0.0312.8 ± 3.20.70
J0307+7443380 ± 150DMHessels et al. (2011)(2.6 ± 1.1 )LAT<0.1−0.04>22>1.00
J0312−0921810 ± 320DMTabassum et al. (2021)32.2 ± 0.6LAT7.62−0.359.7 ± 2.40.63
J0318+02531320 ± 530DMWang et al. (2021)(10.3 ± 13.1 )LAT<22.4−0.21> − 1.3> − 0.26
J0340+41301600 ± 640DMArzoumanian et al. (2018)3.4 ± 0.3Same0.150.226.1 ± 0.10.94
J0418+66352200 ± 880DMTabassum et al. (2021)(2.3 ± 1.0 )LAT<0.30.14>21>0.97
J0437−4715150 ± 10PReardon et al. (2016)140.91 ± 0.00Same43.57−0.282.9 ± 0.60.24
J0533+67592390 ± 950DMSanpa-Arsa (2016)(1.6 ± 0.7 )LAT<0.20.24>5.7>0.96
J0605+3757210 ± 80DMSanpa-Arsa (2016)   0.04≤9.1≤0.99
J0610−21002200 ± 700Pvan der Wateren et al. (2022)18.80 ± 0.03Same7.29−0.073.5 ± 1.60.41
J0613−0200780 ± 70PDesvignes et al. (2016)10.51 ± 0.01Reardon et al. (2021)0.640.0712.2 ± 0.10.93
J0614−3329630 ± 100OBassa et al. (2016a)2.0 ± 0.1Guillemot et al. (2016)0.02−0.1822.69 ± 0.021.0
J0621+25141640 ± 650DMSpiewak et al. (2022)(10.9 ± 52.0 )LAT<142.90.32> − 232> − 4.76
J0636+5128 PAlam et al. (2021)3.63 ± 0.03NRT0.070.06 0.96
J0653+4706910 ± 360DMBassa et al. (2018)(1.6 ± 2.0)JBO<0.30.17>7.4>0.98
J0737−3039A735 ± 60PKramer et al. (2021)3.31 ± 0.03Same0.44−0.885.942 ± 0.0011.0
J0740+66201150 ± 160atnfAFonseca et al. (2021)32.55 ± 0.04Cromartie et al. (2020)8.53−0.046.1 ± 2.00.30
J0751+18071110 ± 200PDesvignes et al. (2016)13.7 ± 0.3Same1.760.095.6 ± 0.30.76
J0931−19023720 ± 1480DMSpiewak et al. (2022)5.0 ± 0.1Alam et al. (2021)1.04−0.731.3 ± 0.20.92
J0952−0607 OKandel & Romani (2022)   −0.14≤68≤1.03
J0955−3947490 ± 190DMYao et al. (2017)   −0.10≤179≤1.00
J0955−61502170 ± 860DMSpiewak et al. (2022)1.8 ± 0.3LAT0.03−0.4073.3 ± 0.81.0
J1012−4235370 ± 140DMSpiewak et al. (2022)6.7 ± 1.0LAT0.13−0.128.7 ± 0.11.00
J1023+00381360 ± 40PBrown et al. (2021)17.98 ± 0.04Deller et al. (2012)1.82−0.2442.8 ± 0.40.77
J1024−07191220 ± 160PBassa et al. (2016b)59.76 ± 0.03Reardon et al. (2021)54.63−0.74-10.2 ± 2.1−1.9
J1035−67201460 ± 580DMSpiewak et al. (2022)11.0 ± 1.4LAT1.22−0.3776.2 ± 1.00.98
J1036−8317 PBrown et al. (2021)11.7 ± 0.2Gaia Collab. et al. (2021)2.44−0.78 0.95
J1048+2339 PBrown et al. (2021)19.3 ± 0.3Gaia Collab. et al. (2021)5.60−0.63 0.83
J1124−3653990 ± 390DMHessels et al. (2011)5.0 ± 0.5Bangale et al. (2023)0.14−0.3016.7 ± 0.31.0
J1125−58251740 ± 690DMSpiewak et al. (2022)10.3 ± 0.3Ng et al. (2014)1.39−0.4379.5 ± 0.80.98
J1125−6014980 ± 390DMSpiewak et al. (2022)17.13 ± 0.01Reardon et al. (2021)1.85−0.194.5 ± 1.60.55
J1137+75283830 ± 1530DMSanpa-Arsa (2016)   −0.16≤8.4≤1.05
J1142+01192160 ± 860DMSanpa-Arsa (2016)(105.0 ± 62.0 )Same<1404.1−1.26> − 418> − 92.52
J1207−50501260 ± 500DMSpiewak et al. (2022)7.0 ± 0.5Bhattacharyya et al. (2021)0.73−0.522.0 ± 0.10.97
J1221−06331250 ± 500DMMcEwen et al. (2020)(44.2 ± 20.0 )LAT<41.5−0.41> − 165> − 2.79
J1227−4853 PBrown et al. (2021)20.1 ± 0.2Jennings et al. (2018)2.92−0.26 0.80
J1231−1411420 ± 100OTesta et al. (2015)62.0 ± 0.1NRT14.41−0.395.7 ± 2.70.34
J1259−81481420 ± 560DMYao et al. (2017)   −0.28≤16≤1.08
J1301+08331770 ± 110OKandel & Romani (2022)22.8 ± 3.5LAT4.11−0.5043.7 ± 8.20.66
J1302−32581430 ± 570DMSpiewak et al. (2022)13.2 ± 0.7NRT2.28−0.633.6 ± 0.70.75
J1306−60431310 ± 520DMYao et al. (2017)   −0.35≤6.7≤1.01
J1311−3430 PBrown et al. (2021)9.1 ± 0.3LAT0.82−0.43 0.98
J1312+00511470 ± 580DMSanpa-Arsa (2016)25.0 ± 0.9Same9.45−1.054.5 ± 2.00.51
J1327−07553700 ± 1480DMBoyles et al. (2013)(6.4 ± 4.2)NRT<5.3−1.15> − 4.9> − 1.31
J1335−5656<15040DMMYao et al. (2017)5.5 ± 1.2C. J. Clark et al. (2023, in preparation)<3.5(-2.9, -0.0)>4.0>0.71
J1400−1431 PSwiggum et al. (2017)42.9 ± 2.2NRT3.82−0.20 0.50
J1402+1306690 ± 270DMYao et al. (2017)   −1.00≤2.8≤1.07
J1431−4715 PBrown et al. (2021)18.8 ± 0.3Jennings et al. (2018)3.06−0.10 0.79
J1446−47011560 ± 620DMSpiewak et al. (2022)5.3 ± 0.1Reardon et al. (2021)0.24−0.0636.2 ± 0.40.98
J1455−3330680 ± 270DMSpiewak et al. (2022)8.2 ± 0.1Alam et al. (2021)0.88−0.231.83 ± 0.030.97
J1513−25503960 ± 1580DMSpiewak et al. (2022)7.2 ± 0.2NRT1.06−0.51 0.97
J1514−4946900 ± 360DMSpiewak et al. (2022)14.0 ± 0.7LAT1.550.0514.6 ± 0.50.91
J1526−27441320 ± 520DMYao et al. (2017)   −0.08≤9.2≤1.02
J1536−4948970 ± 390DMBhattacharyya et al. (2021)7.8 ± 0.3Same0.450.0927.9 ± 0.20.97
J1543−51491140 ± 450DMSpiewak et al. (2022)4.3 ± 0.5PKS0.110.0772.2 ± 0.20.99
J1544+4937 OTang et al. (2014)10.3 ± 0.3NRT1.66−0.68 0.69
J1552+54372630 ± 1050DMPleunis et al. (2017)(3.4 ± 2.3 )LAT<1.0−0.68>6.9>0.88
J1555−2908 OKennedy et al. (2022)(1.1 ± 1.9)LAT<0.5−0.28>306>0.99
J1600−3053 PAlam et al. (2021)7.03 ± 0.03Reardon et al. (2021)0.830.15 0.90
J1614−2230700 ± 300PMatthews et al. (2016)32.4 ± 0.6Alam et al. (2021)5.620.0185.0 ± 3.10.41
J1622−0315 PBrown et al. (2021)13.4 ± 0.3Gaia Collab. et al. (2021)2.06−0.24 0.84
J1623−69361290 ± 510DMYao et al. (2017)   −0.10≤26≤1.01
J1625−0021950 ± 380DMDeneva et al. (2021)(4.0 ± 1.8 )LAT<0.4−0.18>37>0.99
J1627+32194470 ± 1780DMYao et al. (2017)   −1.04≤25≤1.19
J1628−3205 PBrown et al. (2021)22.3 ± 0.4Gaia Collab. et al. (2021)4.620.20 0.67
J1630+37341180 ± 470DMSanpa-Arsa (2016)18.5 ± 0.2JBO3.28−0.628.7 ± 1.40.75
J1640+2224 PVigeland et al. (2018)11.53 ± 0.02Arzoumanian et al. (2018)1.55−0.572.3 ± 0.20.65
J1641+80493030 ± 1210DMLynch et al. (2018)(0.6 ± 1.0 )LAT<0.1−0.28>48>1.02
J1649−3012<15820DMMYao et al. (2017)(10.1 ± 11.4 )LAT (-6.3, 1.3)  
J1653−0158 PBrown et al. (2021)17.6 ± 0.9Gaia Collab. et al. (2021)1.10−0.05 0.56
J1658−5324880 ± 350DMSpiewak et al. (2022)4.9 ± 0.3Camilo et al. (2015)0.130.1029.3 ± 0.20.98
J1713+07471310 ± 30PReardon et al. (2021)6.29 ± 0.00Same0.58−0.293.41 ± 0.010.97
J1730−2304 PReardon et al. (2016)20.3 ± 0.4Reardon et al. (2021)5.040.411.1 ± 0.10.73
J1732−50491870 ± 740DMSpiewak et al. (2022)9.9 ± 0.2Reardon et al. (2021)2.360.453.0 ± 0.30.80
J1741+1351 PAlam et al. (2021)11.64 ± 0.02Arzoumanian et al. (2018)2.69−0.49 0.93
J1744−1134 PReardon et al. (2016)21.03 ± 0.01Reardon et al. (2021)1.730.09 0.80
J1744−7619<1500OClark et al. (2018)18.7 ± 1.2LAT<6.0(−0.4, −0.2)>2.2>0.40
J1745+10171210 ± 480DMSpiewak et al. (2022)8.7 ± 0.2NRT0.59−0.124.4 ± 0.50.81
J1747−40367140 ± 2850DMSpiewak et al. (2022)2.7 ± 0.5Alam et al. (2021)0.210.12 0.98
J1757−60323490 ± 1390DMYao et al. (2017)   −0.30≤5.2≤1.10
J1803−6707 PBrown et al. (2021)10.6 ± 0.2Gaia Collab. et al. (2021)1.08−0.16 0.95
J1805+06153880 ± 1550DMDeneva et al. (2021)15.0 ± 3.0Same4.52−0.3575.9 ± 10.50.82
J1810+1744 PBrown et al. (2021)8.7 ± 0.4LAT0.57−0.18 0.92
J1811−24051830 ± 730DMSpiewak et al. (2022)0.5 ± 0.1Ng et al. (2020)0.000.4727.1 ± 0.50.96
J1816+4510 PBrown et al. (2021)4.4 ± 0.1Gaia Collab. et al. (2021)0.51−1.12 1.0
J1823−3021A7600 ± 400OKuulkers et al. (2003)6.94 ± 0.03Vasiliev & Baumgardt (2021)4.84−0.54 1.00
J1824+10142900 ± 1160DMDeneva et al. (2021)   −0.44≤3.5≤1.08
J1824−06213530 ± 1410DMMiao et al. (2023)   0.47≤10≤0.95
J1824−2452A5370 ± 100OBaumgardt & Vasiliev (2021)8.93 ± 0.03Vasiliev & Baumgardt (2021)3.181.792236.2 ± 0.11.00
J1827−0849<22290DMMAppendix A (1.8 ± 3.9)LAT (−3.4, 0.4)  
J1832−0836 PAlam et al. (2021)22.9 ± 0.1Same7.040.33 0.11
J1833−38404640 ± 1850DMAppendix A    0.43≤105≤0.98
J1835−3259B PBrown et al. (2021)6.95 ± 0.03Vasiliev & Baumgardt (2021)1.80−2.26 1.0
J1843−11131260 ± 500PDesvignes et al. (2016)3.7 ± 0.3Same0.080.1458.2 ± 0.40.98
J1852−13101270 ± 500DMYao et al. (2017)6.8 ± 1.3JBO0.620.314.6 ± 0.20.91
J1855−14365120 ± 2050DMSpiewak et al. (2022)   0.24≤9.1≤0.98
J1857+09431200 ± 350PDesvignes et al. (2016)6.03 ± 0.02Reardon et al. (2021)0.57−0.034.44 ± 0.040.97
J1858−2216920 ± 360DMSanpa-Arsa (2016)6.9 ± 0.4LAT0.260.1110.3 ± 0.30.91
J1858−54221210 ± 480DMYao et al. (2017)   −0.05≤13≤1.01
J1901−01252360 ± 940DMAppendix A (5.8 ± 2.4 )LAT<1.80.11>61>0.95
J1902−51051640 ± 650DMSpiewak et al. (2022)8.0 ± 0.5LAT0.44−0.0264.3 ± 1.40.95
J1903−7051930 ± 370DMSpiewak et al. (2022)18.3 ± 2.0Camilo et al. (2015)2.72−0.266.7 ± 1.10.76
J1908+21052600 ± 1040DMDeneva et al. (2021)8.3 ± 2.4Gaia Collab. et al. (2021)1.13−0.4430.9 ± 1.90.95
J1909−37441140 ± 10PReardon et al. (2016)37.02 ± 0.00Reardon et al. (2021)11.200.064.3 ± 0.20.20
J1921+01375080 ± 2030DMSpiewak et al. (2022)4.3 ± 1.4LAT0.57−0.82 1.0
J1921+19292430 ± 970DMParent et al. (2019)11.5 ± 1.0Same2.07−0.3877.4 ± 2.00.96
J1939+2134 PAlam et al. (2021)0.41 ± 0.00Reardon et al. (2021)0.00−0.351097 ± 11.0
J1946+34176940 ± 2770DMAlam et al. (2021)8.9 ± 0.1Same4.27−1.81 0.22
J1946−54031150 ± 460DMSpiewak et al. (2022)(2.8 ± 2.1 )LAT<0.4−0.17>5.0>0.93
J1959+2048 PBrown et al. (2021)30.4 ± 0.6Arzoumanian et al. (1994)4.59−0.12 0.73
J2006+01482430 ± 970DMDeneva et al. (2021)(3.0 ± 2.5 )LAT<0.8−0.31>11>0.85
J2017+0603 OTesta et al. (2015)2.2 ± 0.1Arzoumanian et al. (2018)0.05−0.26 1.0
J2017−16141440 ± 570DMSanpa-Arsa (2016)(11.0 ± 12.8 )LAT<10.9−0.18> − 26> − 3.44
J2029−4239640 ± 250DMYao et al. (2017)   −0.35≤2.6≤1.04
J2034+3632<14110DMMYao et al. (2017)14.2 ± 0.7Ajello et al. (2022)<25.2(-1.9, -0.1)>20.3>-13.56
J2039−36161700 ± 680DMSpiewak et al. (2022)13.6 ± 0.4Swiggum et al. (2023)2.52−0.437.2 ± 1.20.75
J2039−5617 PBrown et al. (2021)15.5 ± 0.3Clark et al. (2021)2.60−0.40 0.84
J2042+0246640 ± 250DMSanpa-Arsa (2016)6.9 ± 0.9LAT0.34−0.296.0 ± 0.11.00
J2043+1711 PAlam et al. (2021)12.25 ± 0.02Same1.21−0.31 0.84
J2047+10532790 ± 1110DMRay et al. (2012)17.1 ± 1.5JBO8.52−1.196.8 ± 1.90.65
J2051−08271460 ± 580DMShaifullah et al. (2016)6.1 ± 0.1Same0.60−0.645.5 ± 0.11.0
J2052+1219 OKandel & Romani (2022)14.6 ± 0.6Deneva et al. (2021)5.74−1.1110.5 ± 3.70.31
J2115+54483100 ± 1240DMSanpa-Arsa (2016)17.5 ± 0.3NRT6.04−0.56154.2 ± 5.40.93
J2116+13452730 ± 1090DMYao et al. (2017)11 ± 3Lewis et al. (2023)1.8−0.65.4 ± 3.20.56
J2124−3358 PReardon et al. (2016)52.28 ± 0.04Reardon et al. (2021)13.42−0.42 0.37
J2129−04291830 ± 110OBellm et al. (2016)16.0 ± 0.1Jennings et al. (2018)8.66−1.7520.61 ± 0.050.97
J2205+60123530 ± 1410DMDesvignes et al. (2022)3.78 ± 0.01Same0.30−0.4256.2 ± 0.41.0
J2214+3000600 ± 310PGuillemot et al. (2016)21.0 ± 0.1Alam et al. (2021)2.01−0.2717.2 ± 1.40.88
J2215+5135 PBrown et al. (2021)2.2 ± 0.2Gaia Collab. et al. (2021)0.08−0.45 1.0
J2234+0944 PAlam et al. (2021)32.9 ± 0.2Same13.60−0.79 0.35
J2241−52361040 ± 40PReardon et al. (2021)19.61 ± 0.01Same2.13−0.4226.3 ± 0.30.80
J2256−1024 PCrowter et al. (2020)6.0 ± 0.1GBT0.42−0.70 1.0
J2302+4442860 ± 340DMAlam et al. (2021)5.8 ± 0.1Same0.37−0.423.8 ± 0.11.0
J2310−05551550 ± 620DMSanpa-Arsa (2016)   −0.67≤12≤1.14
J2317+1439 PAlam et al. (2021)3.69 ± 0.03Same0.19−0.762.89 ± 0.061.2
J2339−0533 PBrown et al. (2021)11.1 ± 0.4Jennings et al. (2018)1.12−0.67 0.97

Note. Columns 2 and 3 give the distances for the 144 MSPs detected with the Fermi LAT, and the method used to find them: P—parallax; LK—corrected for the Lutz–Kelker effect; DM—from dispersion measure (DM) using the YMW16 model (Yao et al. 2017); O—optical and other methods. For DM, the references in Column 4 give the dispersion measurement. "DMM" is a limit on the pulsar distance, corresponding to 95% of the maximum YMW16 DM for that line of sight. Proper-motion measurements μ are listed in Column 5, taken from the references in Column 6, often the "same" as in Column 4. "LAT" refers to gamma-ray timing, described in Section 3.5. "NRT," "JBO," "PKS," and "GBT" indicate that the measurements were made for this work using data from the Nançay, Lovell, Parkes, and Green Bank radio telescopes. Columns 7 and 8 list the contributions of the Shklovskii effect and of acceleration due to the Galactic rotation, to the intrinsic spindown rate . The last columns list the intrinsic spindown power . See the text for how cases with a missing or insignificant μ measurement, or a distance limit, are handled.

Download table as:  ASCIITypeset images: 1 2 3 4

DM is the column density of free electrons along the path from Earth to the pulsar, in units of pc cm−3. The electrons delay the radio pulse arrival by where ν is the observation frequency in megahertz and p = 2.410 ×10−4 MHz−2 pc cm−3 s−1. Given a model for the electron density ne in the various structures of our Galaxy, integrating DM = along the line of sight dl yields the distance d for which DM matches the radio measurement. We use mainly the YMW16 model (Yao et al. 2017), but often consult the NE2001 model when studying particular lines of sight (Cordes & Lazio 2002). Significant discrepancies between model and true distances occur along some lines of sight, more than three times the probable true distance in the worst cases. Examples are PSR J0614−3329 (Bassa et al. 2016a); PSR J2021+3651 (Abdo et al. 2009c); PSR J0248+6021 (Theureau et al. 2011); and PSR J0908−4913 (Johnston & Lower 2021). Figure 4 of Hou et al. (2014) illustrates one way to identify pulsar lines of sight that have DM model peculiarities.

Radio spectra can show Doppler-shifted neutral hydrogen (H i) absorption or emission lines due to clouds along the line of sight. A Galactic rotation model converts the Doppler-inferred radial velocity to a cloud distance estimate, and the pulsar distance is constrained if it appears to be in or between clouds. This "kinematic" method is tagged with a "K" in the tables. The "X" tag indicates cases where X-ray absorption below 1 keV yields a column density that can be compared with one obtained from 21 cm radio surveys.

To obtain the distances, we first query the ATNF psrcat DIST parameter, which is most often DIST_DM, the DM distance using YMW16. We adopt a 40% uncertainty on DIST_DM as per Figure 12 and Table 5 of Yao et al. (2017). When psrcat has a DIST_A value, obtained by another method, we use that instead and propagate the method and reference to our tables. The exception is that we exclude the distances deduced by Wang (2011) using the LAT G100 values and the correlation between Lγ and shown in Section 6.7, since they would bias our attempts to better determine that correlation. If the psrcat parallax parameter PX is present, DIST_A = 1/PX, and we use the reference for PX. Parallax is measured only for relatively nearby pulsars, using X-ray or optical images, radio interferometric imaging, or accurate timing. 119 When converting parallax measurements to distances, the Lutz–Kelker (LK) effect is an overestimate of parallax values, and hence an underestimate of distances, due to the larger volume of space traced by smaller parallax values (Lutz & Kelker 1973). We use the LK-corrected distance estimates determined by Verbiest et al. (2012). Good estimates of the distance to an optical companion in a binary system are more frequent than in 2PC because we have more spider MSPs than before, and because of the wealth of information provided by Gaia (Gaia Collab et al. 2018). For two of our MSPs, PSRs J1227−4853 and J1431−4715, Jennings et al. (2018) found distances using Gaia companion parallaxes. They used a spatial distribution prior specific to Galactic pulsars, i.e., not isotropic, which means that the LK-like correction brings some pulsars closer to Earth when they are toward the anticenter, as opposed to farther for standard LK. For these two MSPs, we simply use 1/PX, the same within uncertainties as the Jennings et al. (2018) distances.

For several pulsars, measured DM exceeds DMmax, the maximum modeled by Yao et al. (2017) for that line of sight, by between 1% and 50%. The YMW16 model therefore returns the distance to the Galaxy's edge. Modeled DM saturates much closer to Earth, increasing slowly thereafter. For example, toward PSR J1327−0755, 95% of DMmax occurs at 3.7 kpc, compared to 25 kpc returned by the code, and in better agreement with the 1.7 kpc NE2001 distance. Therefore, for pulsars with DM > DMmax, we use the YMW16 distance for 95% of DMmax. Other pulsars cited in this catalog for which this occurs are PSR J1102+02 (Table 6), and PSRs J1302−03, J2051+4434g and J2052+4421G (Table 5).

For 63 pulsars we have no DM measurement or other distance estimator. We again calculate the YMW16 distance for 95% of DMmax, but use it as an upper limit. The method is "DMM" in Tables 7 and 8, and these radio-quiet pulsars are illustrated in Figure 8.

Both Espinoza et al. (2013) and 2PC discussed the problem of PSR J0610−2100 having gamma-ray efficency , using the DM distances of 3.3 and 4 kpc. van der Wateren et al. (2022) used LOFAR timing and optical modeling for this system to find d = 2.2 ± 0.7 kpc. Espinoza et al. (2013) pointed out that the line of sight to the pulsar grazes the electron overdensities bordering the Local Bubble, in both electron models. A slight extension of the modeled bubble yields the increase in modeled DM needed to reconcile the distances. The closer distance for PSR J0610−2100 leads to η ≈ 50%. It also reduces the pulsar's transverse velocity to vT ≈ 220 km s−1, so that this MSP is no longer among the fastest known.

4.3. Doppler Corrections

Many pulsar characteristics depend on the intrinsic spin period Pint and spindown rate . The Doppler shift of the observed period is P = (1 + v R/c)Pint, where v R is the pulsar's radial velocity along the unit vector n from the solar system. The Doppler correction to is obtained by differentiating the equation and separating the effects of the system's proper motion (Shklovskii 1970) from the acceleration due to Galactic rotation (the treatment by Damour & Taylor 1991 is particularly clear). The result is

with

and

where k = 2.43 × 10−21 for pulsar distance d and proper-motion transverse to the line of sight μ. In 2PC, we obtained the accelerations a 1 of the pulsar and a 0 of the Sun using the Galactic potential model of Carlberg & Innanen (1987) and Kuijken & Gilmore (1989). Here, we instead use GalPot 120 with the PJM17_best.Tpot parameter set (McMillan 2017).

The corrections are negligible for the young gamma-ray pulsars, which all have . However, for MSPs can differ noticeably from observed . Expressing the correction magnitude as

leads to . For μ ≲ 10 mas yr−1 and d ≲ 1 kpc, and have comparable magnitudes and opposite signs. For larger μ2 d, dominates.

Table 8 lists d, μ, their uncertainties δ d, δ μ, and the resulting corrections. When we have no μ measurement and only a limit for d, we leave the last columns empty. Parentheses around the μ ± δ μ value indicate that the measurement is insignificant, δ μ/μ > 3, and we express the remaining columns as limits calculated using . In parentheses in the column are shown the minimum and maximum values for that direction for distances from 500 pc to the limit. Given a distance but no proper-motion value, we calculate and limits on ξ, using . There are eight such MSPs, of which PSR J0605+3757 is listed first. In the same spirit, e.g., PSR J1335−5656 has measured proper motion but unknown distance, and we combine the limit for with the extremum from the range.

Table 8 lists a single MSP with well-measured μ but (1 − ξ) < 0, implying implausible spin-up, . It is PSR J1024−0719, with a binary companion in an atypical multicentury orbit that leads to mismeasured (Bassa et al. 2016b; Guillemot et al. 2016; Kaplan et al. 2016). Six other MSPs have erg s−1, some with nominally small uncertainties. The lowest is erg s−1 for PSR J1946+3417. Improved precision and a larger sample will clarify the gamma-ray "deathline" for MSPs, below which emission in the LAT range ceases. A useful tool to understand the uncertainty on are plots of μ versus d showing the constraints suggested by the neutron star's transverse velocity v T = μ d, minimum , and so on. A powerful constraint is to require (Section 6.7 addresses luminosity Lγ ). Examples of μ versus d plots, or similar, are Figure 11 in 2PC (which includes PSR J1024−0719), Figure 1 in Guillemot et al. (2016), and Figure 4 in Clark et al. (2018).

5. Light Curves and Profile Characterization

5.1. Gamma-Ray and Radio Light Curves

Figure 9 shows an example of gamma-ray light curves in different energy bands for the MSP PSR J0030+0451. All 294 sources are available in the Figure Set. Figure 4 shows further examples, overlaid with radio profiles when available.

Figure 9. Refer to the following caption and surrounding text.

Figure 9.

Pulse profiles for the MSP PSR J0030+0451. Weighted LAT pulse profiles in different energy bands are shown in black. The pulse profiles above 10 GeV are in light gray, more apparent for other pulsars shown in Appendix B. The radio pulse profile for this pulsar at 1.4 GHz in shown in red, and a fit of the integrated gamma-ray pulse profile over 0.1 GeV is displayed in magenta when the H-test test significance exceeds 100. Horizontal blue lines correspond to the expected background levels in the gamma-ray pulse profiles, plus or minus one standard deviation. The quantities at the top left correspond to the weighted H-test TS for the LAT pulse profile shown in the top panel, and the total integration time. At top right are listed: the pulsar name, its rotational period and Galactic latitude, and the "PSR" and "CHAR" codes defined in the captions for Tables 2 and 3 (top row), the radio observatory used to acquire the radio profile (if available), and the corresponding radio frequency (second row). For pulsars with radio profiles, a marker at the top right indicates the phase alignment uncertainty stemming from the DM. See Section 5 for more details about light-curve construction and the fitting procedure. (The complete figure set (294 images) is available.)

Standard image High-resolution image

After preparation according to Section 2, the photon arrival times were converted to rotational phase using timing models obtained from pulsar timing observations conducted in radio or X-rays, or by direct timing of the LAT photons, as described in Section 3. In a number of cases, the timing models' validity interval only partially overlapped with the time range of the LAT data set. To avoid contaminating off-pulse phase regions by extrapolating inaccurate models, we restricted the LAT data sets to the validity intervals in the cases of (1) timing solutions containing high-order (≥3) frequency derivatives suggestive of strong timing noise; (2) "IFUNC" or "WAVE" terms (also used to mitigate timing noise); or (3) "ORBIFUNC" or high-order (≥2) orbital frequency derivatives for pulsars in binary systems, or black-widow pulsars that are often subject to timing instabilities. For other objects, generally stable MSPs, all LAT photons were used, and we checked that the pulse phases do not drift outside the ephemeris validity intervals (see Figure 4). Additionally, we conservatively excluded 10 days of data before and after pulsar glitch epochs, to avoid contamination of off-pulse phase intervals caused by imperfect glitch modeling. At the top of all of the plots like Figure 9, tint gives the number of years of LAT data retained.

Pulse profiles as in Figures 4 and 9 were generated by computing histograms of the rotational phases ϕj , weighted by the photon weights wj . The weighted phase histograms repeat over two pulsar rotations to clarify structures near phase 0. The uncertainty for the ith histogram bin containing Ni photons was calculated as . The first term is the sample variance, while the second term adds in quadrature the maximum possible variance for a bin with one photon. This is an ad hoc correction for bins with few counts, where by chance the observed sample variance may be far below the typical value: this prescription provides more realistic error bars than those used in 2PC. The number of bins in the weighted histograms was chosen based on the H-test TS above 0.1 GeV: 25 bins for pulsars with H < 100, 50 bins when 100 ≤ H < 1000, 100 bins when 1000 ≤ H < 10,000, and so on.

Dashed horizontal lines in the figures show the background contribution from neighboring point sources and from diffuse emission. As in 2PC, it was estimated as where N is the number of sources, excluding the pulsar, and Nbin is the number of histogram bins. The background uncertainty is dominated by systematic errors in the diffuse background normalization, estimated by changing the background normalization by ±6%, also shown as the dashed horizontal lines. A pulsar wind nebula (PWN) may surround the pulsar, but with a power-law spectrum extending to high energies. If omitted from the source model, it appears as an excess above background, especially at higher energies. Section 6.2 lists the pulsars for which we added a PWN component to our source model, in addition to the Crab and Vela. 2PC describes a search for off-pulse emission, which measured off-pulse spectra for pulsars with an excess above the background level.

For gamma-ray pulsars with no radio emission, we rotated the pulse profiles such that the phase of the first emission component (as determined from the light-curve fit above 0.1 GeV, described in the next section) is at phase 0.1 (see the profile for PSR J0007+7307 in Figure 9.2 available online for an example). When radio profiles are available, a radio fiducial phase at infinite frequency was determined (see below) and placed at phase 0. The gamma-ray light curves were rotated accordingly, to preserve the relative radio/gamma-ray phase offset, including DM delays (see Section 4.2). We determined the uncertainty on the radio/gamma-ray phase offset stemming from the DM uncertainty neglecting DM time-derivative terms, and show it in the figures with a horizontal error bar at the top-right. Figure 9.104 for PSR J1311-3430 available online  shows a case where a large DM uncertainty leads to a large relative phase uncertainty. For some pulsars, the DM and its uncertainty were not determined during the radio timing analysis; for these objects we used external DM information, generally from the ATNF pulsar catalog or separate analyses. Those cases are highlighted with a marker underneath the phasing error bar, as illustrated in Figure 9.285 available online showing the pulse profile for PSR J2229+6114. In a very few other cases, the radio data analysis did not yield the information necessary to determine the relative radio/gamma-ray phase offset. For these light curves, we indicated the arbitrary relative phase offset with a ⋆ symbol.

Identifying the "fiducial phase" at which the magnetic pole of the neutron star crosses the line of sight is of interest. For some pulsars, this can be inferred from radio polarimetry via the rotating vector model (RVM; Radhakrishnan & Cooke 1969), but we do not consider such models here. As noted in 2PC, for pulsars with a single, symmetric radio peak, the peak maximum can be chosen as an estimator for the fiducial phase. However, many radio profiles in our sample do not consist simply of single, symmetric peaks, so that the choice of the fiducial phase is not always trivial. We adopted the following "Gaussian smoothing" approach, which reduces complex emission structures by replacing them with the best-fitting Gaussian equivalent. Each radio profile was inspected visually, and one peak or peak complex was chosen as the phase reference. This reference peak or set of peaks was then fitted with a single asymmetric Gaussian curve, or with two asymmetric Gaussian curves, depending on the morphology. Single peaks or "patchy" profiles with more than two components were fitted with single asymmetric Gaussian functions. Peaks with two components at their edges were fitted with two asymmetric Gaussians. In the first case, the fiducial phase was simply defined as the best-fit peak phase, and in the second case, the fiducial phase was defined as the average of the two best-fit peak phase values.

While the above procedure for defining the radio fiducial phase works well in most cases, a number of radio profiles in the sample are too complex for the procedure to produce an estimate. We inspected the radio profiles and classified them into four categories:

  • 1.  
    radio profiles for which one can readily define one peak as the main peak, which has a simple, single, or double structure. The procedure for defining the radio fiducial phase should have produced robust results.
  • 2.  
    radio profiles for which one can readily define one peak as the main peak, whose structure is complex (i.e., patchy profiles). The procedure is believed to be robust.
  • 3.  
    one can define a radio peak as the main one, but cannot readily determine whether it is isolated, or part of a broader structure with other peaks. The radio fiducial phase estimate may not be robust for these objects.
  • 4.  
    the profile morphology is too complex for the fiducial phase to be unambiguously chosen, based on the profile alone.

The categories, or "radio classes," assigned to the different pulsars are listed in Tables 9 and 10. The fiducial phases were finally placed at phase 0 in the radio and gamma-ray pulse profiles, and the gamma-ray phases were rotated accordingly.

Table 9. Pulse Shape Parameters of Young LAT-detected Pulsars

PSRPeaksRadio/Gamma Phase Lag, δ Peak Separation, ΔCommentRadio Class
J0002+62162  1
J0007+730330.2383 ± 0.0014 
J0106+485520.057 ± 0.008 u1
J0205+644920.0400 ± 0.00240.5039 ± 0.0019 1
J0248+602110.347 ± 0.005 1
J0357+320520.192 ± 0.005 
J0359+54141 
J0514−440820.4776 ± 0.0034 u3
J0534+220020.067 ± 0.0040.39889 ± 0.00017 2
J0540−691910.96 ± 0.04 1
J0554+310720.339 ± 0.008 
J0622+374920.544 ± 0.027 
J0631+064620.332 ± 0.0140.311 ± 0.015 1
J0631+103610.507 ± 0.007 2
J0633+063220.4776 ± 0.0009 
J0633+174630.50398 ± 0.00015 
J0659+141410.226 ± 0.004 1
J0729−144810.572 ± 0.015 2
J0734−15591 
J0742−28222  2
J0744−252530.468 ± 0.004 
J0802−561320.261 ± 0.023 
J0835−451040.133 ± 0.0040.431865 ± 0.000032 1
J0908−491320.0971 ± 0.00340.500 ± 0.004 2
J0940−542810.454 ± 0.010 1
J1016−58573  1
J1019−574910.566 ± 0.010 1
J1023−574630.4692 ± 0.0005 
J1028−581920.1935 ± 0.00100.4727 ± 0.0004 1
J1044−573730.3604 ± 0.0018 
J1048−583230.1468 ± 0.00120.4221 ± 0.0009 2
J1055−60281  1
J1057−522640.7441 ± 0.00170.3000 ± 0.0010c4
J1057−58511 
J1105−603720.348 ± 0.014 
J1105−61072 0.501 ± 0.006 1
J1111−603920.313 ± 0.018 
J1112−610320.147 ± 0.0180.468 ± 0.020 1
J1119−612710.469 ± 0.008 1
J1124−591620.136 ± 0.0090.4985 ± 0.0011 1
J1135−605520.173 ± 0.013 
J1139−624720.151 ± 0.016 
J1151−61082  2
J1203−624230.468 ± 0.005 
J1208−623820.468 ± 0.007 
J1224−640710.215 ± 0.018 1
J1231−511320.431 ± 0.008 
J1231−65111 
J1350−622530.485 ± 0.005 
J1357−642910.388 ± 0.009 1
J1358−602520.298 ± 0.028 
J1410−613210.537 ± 0.027 1
J1413−620530.3707 ± 0.0014 
J1418−605830.4626 ± 0.0016 
J1420−60482 0.311 ± 0.005 1
J1422−613820.196 ± 0.012 
J1429−591130.4738 ± 0.0030 
J1447−57571 
J1459−60531 
J1509−585020.156 ± 0.0060.258 ± 0.006 1
J1513−590810.321 ± 0.010 1
J1522−573520.465 ± 0.011 
J1528−58381 
J1531−561010.40 ± 0.04 1
J1615−51372  
J1620−492720.246 ± 0.007 
J1623−500520.389 ± 0.010 
J1624−404120.422 ± 0.005 
J1641−53171 
J1646−434610.535 ± 0.014 1
J1648−461120.263 ± 0.0100.310 ± 0.011 1
J1650−460120.358 ± 0.011 
J1702−41281  1
J1705−19061  2
J1709−442930.2364 ± 0.00040.2515 ± 0.0004 1
J1714−38301 
J1718−382510.398 ± 0.005 2
J1730−33502  1
J1732−313130.0876 ± 0.0026  4
J1736−342220.437 ± 0.009 
J1741−205420.035 ± 0.005  2
J1742−33211 
J1746−323920.133 ± 0.009 
J1747−295830.156 ± 0.0050.392 ± 0.004 1
J1748−28151 
J1801−245120.118 ± 0.0070.495 ± 0.007 1
J1803−214930.379 ± 0.005 
J1809−233230.3465 ± 0.0007 
J1813−124630.489 ± 0.004 
J1817−174220.518 ± 0.005 
J1826−125620.4737 ± 0.0004 
J1827−14461 
J1828−110120.110 ± 0.005  1
J1833−103420.142 ± 0.0160.4530 ± 0.0023 1
J1835−11061  1
J1836+592530.5167 ± 0.0016 
J1837−060420.197 ± 0.0090.486 ± 0.021 1
J1838−053720.304 ± 0.009 
J1844−03461 
J1846+091920.227 ± 0.017 
J1856+01132  1
J1857+01431  1
J1906+072220.222 ± 0.023 
J1907+060230.282 ± 0.0080.3850 ± 0.0013 1
J1913+09042  1
J1932+19161 
J1952+325230.1572 ± 0.00070.4752 ± 0.0006 1
J1954+283630.4521 ± 0.0025 
J1957+503320.178 ± 0.009 
J1958+284630.4371 ± 0.0016 
J2006+310210.570 ± 0.007 1
J2017+36252  
J2021+365130.1369 ± 0.00190.47147 ± 0.00011 1
J2021+402630.553 ± 0.005c
J2022+38421  4
J2028+333230.432 ± 0.005 
J2030+364120.300 ± 0.0070.320 ± 0.008 1
J2030+441530.490 ± 0.004 
J2032+412720.1626 ± 0.00090.4831 ± 0.0008u1
J2043+27402  1
J2055+25392  
J2111+460630.3459 ± 0.0034 
J2139+47162  
J2229+61142  1
J2238+590320.5002 ± 0.0005 
J2240+58322  1

Note. Column 2 gives the gamma-ray peak multiplicity. Columns 3 and 4 give the radio/gamma phase lag δ and separation Δ between the two main gamma-ray peaks. In Column 5, "u" indicates uncertain ordering of the first and second gamma-ray peaks, and "c" means the profile is complex. For pulsars with radio emission, Column 6 lists flags quantifying the robustness of the fiducial phase definition (see the text for details).

Download table as:  ASCIITypeset images: 1 2 3

Table 10. Pulse Shape Parameters of LAT-detected MSPs

PSRPeaksRadio/Gamma Phase Lag, δ Peak Separation, ΔCommentRadio Class
J0023+092320.321 ± 0.0220.280 ± 0.027 1
J0030+045150.1266 ± 0.00140.4485 ± 0.0004 2
J0034−053420.839 ± 0.0040.319 ± 0.006 1
J0101−642240.6733 ± 0.00320.4053 ± 0.0031c1
J0102+483920.327 ± 0.0070.444 ± 0.007 1
J0218+423210.691 ± 0.006 1
J0251+260620.157 ± 0.0070.522 ± 0.007 1
J0307+744320.935 ± 0.0080.645 ± 0.013 1
J0312−09212  1
J0340+413020.7501 ± 0.00270.239 ± 0.004 1
J0418+663530.737 ± 0.0170.391 ± 0.005 3
J0437−471510.4299 ± 0.0023 2
J0533+675920.175 ± 0.005  1
J0605+375720.812 ± 0.027 u2
J0610−210030.615 ± 0.004  1
J0613−020020.224 ± 0.0140.159 ± 0.014 2
J0614−332920.1293 ± 0.00050.5558 ± 0.0004 2
J0621+251420.340 ± 0.0300.388 ± 0.032 2
J0737−3039A2 0.506 ± 0.012u3
J0740+66202 0.4622 ± 0.0029 3
J0751+180720.3694 ± 0.00330.272 ± 0.007 1
J0931−19021  4
J0952−060710.053 ± 0.029 1
J0955−615020.883 ± 0.0140.479 ± 0.016 1
J1012−423520.290 ± 0.0050.496 ± 0.008 1
J1024−071910.561 ± 0.024 2
J1035−67202 0.528 ± 0.016 1
J1036−831710.273 ± 0.019 1
J1048+233920.1371 ± 0.0032  1
J1124−365330.259 ± 0.0100.394 ± 0.005c1
J1125−58252  2
J1142+01191  1
J1207−505020.255 ± 0.006  1
J1221−06331  3
J1227−48533 c1
J1231−141140.2341 ± 0.00110.4056 ± 0.0010c2
J1301+083320.104 ± 0.0160.393 ± 0.024 1
J1302−325820.147 ± 0.0040.394 ± 0.008 3
J1311−343030.906 ± 0.0180.3629 ± 0.0017 1
J1312+005120.914 ± 0.0190.690 ± 0.020u3
J1335−565620.209 ± 0.016 
J1400−143110.913 ± 0.013 1
J1446−470120.339 ± 0.0080.236 ± 0.015 1
J1455−33301  2
J1513−25502  3
J1514−494630.273 ± 0.004  1
J1536−494860.429 ± 0.0070.5934 ± 0.0030c3
J1543−514910.959 ± 0.015 1
J1552+54372  1
J1555−290820.029 ± 0.0040.554 ± 0.007 3
J1600−305320.212 ± 0.0070.186 ± 0.015 1
J1614−22303  1
J1625−002120.992 ± 0.005  1
J1628−320520.244 ± 0.016  1
J1630+373420.2544 ± 0.00130.4567 ± 0.0023 3
J1640+222420.332 ± 0.0310.185 ± 0.033 2
J1641+80492  1
J1649−30121 
J1653−015820.420 ± 0.014 
J1658−532420.105 ± 0.0260.270 ± 0.025u4
J1713+074710.318 ± 0.016 2
J1730−23041  2
J1732−50491  2
J1741+13512 u1
J1744−113420.8116 ± 0.00290.727 ± 0.022u1
J1744−761920.343 ± 0.008 
J1745+10171  1
J1747−403610.877 ± 0.013 1
J1805+06152  1
J1810+174420.050 ± 0.0130.243 ± 0.015 2
J1811−24052  3
J1816+451030.5300 ± 0.00290.481 ± 0.005 1
J1823−3021A20.980 ± 0.0040.390 ± 0.005 1
J1824−2452A20.669 ± 0.0310.530 ± 0.034 3
J1827−08491 
J1832−08361  4
J1843−11131  1
J1855−14361  1
J1858−221620.748 ± 0.0060.306 ± 0.014 3
J1901−012530.866 ± 0.0070.378 ± 0.008 1
J1902−510530.015 ± 0.0040.645 ± 0.004 1
J1903−705110.531 ± 0.012 1
J1908+21052 0.419 ± 0.016 1
J1921+013720.134 ± 0.0300.362 ± 0.032 1
J1939+21342  2
J1946−54033 c1
J1959+204820.9953 ± 0.00140.566 ± 0.007 4
J2006+01482  1
J2017+060320.2244 ± 0.00270.2955 ± 0.0035 3
J2017−161420.183 ± 0.014  3
J2034+36323  
J2039−36161  2
J2039−561730.107 ± 0.0070.296 ± 0.009 1
J2042+024610.766 ± 0.008u1
J2043+171130.2503 ± 0.00280.4479 ± 0.0027 3
J2047+10532  1
J2051−08271  2
J2052+121920.675 ± 0.0110.381 ± 0.012 1
J2115+54482  3
J2124−335830.8532 ± 0.00260.744 ± 0.009 3
J2129−042930.6694 ± 0.0032  3
J2214+30002 0.484 ± 0.012u1
J2215+513530.573 ± 0.0040.447 ± 0.004 1
J2234+094420.680 ± 0.0100.187 ± 0.018 3
J2241−52365 0.317 ± 0.024c1
J2256−102420.1510 ± 0.0010  2
J2302+444230.2318 ± 0.00200.3350 ± 0.0021 2
J2310−05552  3
J2339−053330.3036 ± 0.00340.3734 ± 0.0027 2

Note. Same as Table 9, for the gamma-ray MSPs.

Download table as:  ASCIITypeset images: 1 2

5.2. Gamma-Ray Light-curve Fitting above 0.1 GeV

Fits of the LAT weighted histograms above 0.1 GeV are displayed in the Figures as magenta lines (see, e.g., Figure 9). We followed a similar approach as in 2PC for fitting the LAT pulse profiles, by representing the gamma-ray light curves as wrapped probability density functions (PDFs) of , and constructing the PDFs for the individual pulsars as linear combinations of Npeaks unimodal distributions:

In the above expression, ni terms represent normalization factors, and the gi terms represent individually normalized distributions, and the term in parentheses is a uniform distribution representing an unpulsed component. We chose to fit gamma-ray light-curve peaks using asymmetric Gaussian or Lorentzian distributions, wrapped onto a circle:

with defined on the real line, and k being an integer number of rotations. In practice, the sum needed to be truncated, and we summed phases over with Nrot = 2. The PDFs were defined as follows:

with z defined as , with ϕ0 the phase location of the peak and σ = σ1 if ϕ < ϕ0 and σ = σ2 otherwise. Best-fit normalization, phase, and width parameters for each component were determined using maximum likelihood. The logarithm of the likelihood is defined as:

The fit was carried out using the MultiNest Bayesian inference tool (Feroz et al. 2009). We fitted binned weighted histograms above 0.1 GeV with twice as many bins as displayed in the figures. In some cases, the width parameters converged toward the lower edge of the prior; in these cases, we doubled the number of phase bins for the fit. For each gamma-ray component in the pulse profiles, we tested the different PDFs and chose the one maximizing the log likelihood value. The number of gamma-ray components was determined by adding components at possible peak locations (based on visual inspection) and choosing the model that minimized the Bayesian information criterion (BIC). Pulse profile fits were carried out and plotted in the Figures for the 236 pulsars with H-test values above 100.

In most pulsars with more than two gamma-ray components, one can readily distinguish two "main" gamma-ray peaks at phases ϕ1 and ϕ2 whose separation is denoted as Δ = ϕ2ϕ1. In an archetypal gamma-ray pulse profile such as for the Vela pulsar (PSR J0835−4510), a first main gamma-ray peak lags the radio emission peak by δ = ϕ1ϕr ≃ 0.1 − 0.2 (where ϕr denotes the radio fiducial phase, as discussed previously, if radio emission is also detected), and a second main gamma-ray peak lags the first one by ≲0.5 with emission between the peaks. However, in some cases, one cannot straightforwardly determine which of the gamma-ray peaks should be labeled as peak 1 and peak 2. Therefore, for each profile with more than one statistically significant emission component, we chose peak 1 and peak 2 based on similarity with the above-described archetype. Then, for each pulsar, we considered the following criteria:

  • 1.  
    δ < 0.5,
  • 2.  
    Δ < 0.5,
  • 3.  
    σ2 > σ1 for the first peak (i.e., is the first peak oriented toward the second one),
  • 4.  
    σ1 > σ2 for the second peak (i.e., is the second peak oriented toward the first one),
  • 5.  
    presence of statistically significant emission between the two main peaks,

and inspected those meeting fewer than four of the above criteria (i.e., pulse profiles that do not correspond to the "typical" Vela-like profile). In some of these cases, peak 1 and peak 2 were inverted based on visual inspection, and a number of objects were flagged as "uncertain."

In addition to those cases where the ordering of the two main gamma-ray peaks is ambiguous, some other pulsars have "complex" pulse profiles in which the choice of the two main peaks is itself ambiguous, for various reasons. In an attempt to list these objects, we searched the best-fit pulse profiles for examples that contain at least one peak not overlapping with the two main peaks chosen from initial visual inspection, and whose maximum amplitude is at least half that of the smaller of the two main peaks. The handful of objects that matched this criterion were flagged as "complex."

For all of these pulsars for which the choice of peak 1 and peak 2 is not obvious, the values for δ and Δ are also uncertain.

Tables 9 and 10 list the best-fit parameters. Figure 10 shows a plot of the separation of the two main gamma-ray peaks, Δ, versus the radio/gamma-ray phase separation, δ. Figure 11 shows the Δ values plotted against the spin-down power values, . For the pulsars with radio emission and at least two gamma-ray components (and therefore, with both δ and Δ values in Tables 9 and 10 and in Figure 10), and excluding pulsars with profiles considered as complex or with uncertain peak ordering, we find a Spearman rank coefficient of ∼−0.30, indicative of a negative correlation between δ and Δ with a probability of chance correlation of ∼2 × 10−3. The anticorrelation of δ and Δ, already noted in, e.g., 2PC, is a well-known prediction of outer magnetospheric pulsar emission models. No strong relation is seen in Figure 11. These results and other trends visible in these figures are further discussed in Section 8.

Figure 10. Refer to the following caption and surrounding text.

Figure 10. Phase separation Δ between main gamma-ray peaks, vs. phase lag δ of the first gamma-ray peak relative to the radio emission. Values of δ between 0.9 and 1 are plotted between −0.1 and 0 in this graph. Values between 0.5 and 0.9 have been "folded" as . The markers are the same as in Figure 2. Circled symbols are pulsars for which the ordering of the first and second gamma-ray peaks was considered uncertain. Squares correspond to pulsars with "complex" profiles, in which choosing two main peaks is not straightforward. Opaque symbols are pulsars for which the definition of the radio fiducial point is considered robust, other pulsars are shown as transparent symbols (see the text for details on the definition of radio fiducial points). Pulsars with δ and no Δ measurements (i.e., single-peaked gamma-ray pulsars with radio emission) are shown in the gray-shaded region at the bottom, while those with Δ and no δ values (i.e., pulsars with at least two gamma-ray components and no radio emission) are plotted in the gray-shaded area at the right. The artificial staggering of the points in these two regions is to improve clarity. At the top and at the right are plotted histograms of the values of δ and Δ projected on the axes, for all pulsars shown in the central figure (the "All" panels), and restricting the sample of pulsars to those with δ and Δ measurements (the "Both δ and Δ" panels, i.e., excluding pulsars in the gray-shaded areas).

Standard image High-resolution image
Figure 11. Refer to the following caption and surrounding text.

Figure 11. Phase separation Δ between main gamma-ray peaks vs. Shklovskii-corrected spindown power . The markers are the same as in Figure 2. See Figure 10 for the definition of the circle and square symbols.

Standard image High-resolution image

We compared δ and Δ in Tables 9 and 10 with the 2PC values, for pulsars in both catalogs. Table 11 lists pulsars for which δ and/or Δ differ by more than three times the statistical uncertainty, and with an absolute difference >0.03. For Δ, no strong inconsistencies appeared: inverted gamma-ray profile fit components, or differing numbers of components or component phases explain most discrepancies. Similarly, the δ parameters between 2PC and 3PC are qualitatively similar in most cases: differences typically arise from choices of radio pulse reference phase, or choices of the gamma-ray component used to find δ. The only pulsars with significantly discrepant radio/gamma-ray alignment properties in 2PC and 3PC are PSRs J2047+1053 and J2215+5135.

Table 11. Comparisons of Radio/Gamma Phase Lags δ and Peak Separations Δ between 2PC and 3PC, for Pulsars with Apparently Discrepant Results

PSR2PC Value3PC ValueComment
Radio/gamma phase lags, δ
J0030+04510.1600 ± 0.00100.1266 ± 0.0014Different choice of radio reference phase
J0101−64220.145 ± 0.0050.6733 ± 0.0032Different choice of gamma-ray reference peak
J0102+48390.259 ± 0.0040.327 ± 0.007Qualitatively similar results
J0205+64490.075 ± 0.0040.0400 ± 0.0024Qualitatively similar results
J0218+42320.35 ± 0.080.691 ± 0.006Different choice of gamma-ray reference peak
J0534+22000.10900 ± 0.000200.067 ± 0.004Different choice of radio reference phase
J0610−21000.236 ± 0.0060.615 ± 0.004Different choice of radio reference phase
J0742−28220.627 ± 0.005 Different choice of gamma-ray reference peak
J1019−57490.482 ± 0.0100.566 ± 0.010Qualitatively similar results
J1057−52260.3040 ± 0.00290.7441 ± 0.0017Different choice of radio reference phase
J1119−61270.285 ± 0.0150.469 ± 0.008Different choice of gamma-ray reference peak
J1125−58250.6450 ± 0.0017 Different choice of gamma-ray reference peak
J1410−61320.959 ± 0.0230.537 ± 0.027Different choice of gamma-ray reference peak
J1509−58500.271 ± 0.0110.156 ± 0.006Qualitatively similar results
J1514−49460.214 ± 0.0090.273 ± 0.004Different choice of radio reference phase
J1600−30530.147 ± 0.0110.212 ± 0.007Qualitatively similar results
J1658−53240.359 ± 0.0140.105 ± 0.026Different choice of gamma-ray reference peak
J1741−20540.074 ± 0.0060.035 ± 0.005Qualitatively similar results
J1744−11340.189 ± 0.0070.8116 ± 0.0029Different choice of radio reference phase
J1747−40360.031 ± 0.0210.877 ± 0.013Qualitatively similar results
J1801−24510.060 ± 0.0050.118 ± 0.007Qualitatively similar results
J1810+17440.894 ± 0.0180.050 ± 0.013Qualitatively similar results
J1835−11060.139 ± 0.006 Different choice of gamma-ray reference peak
J1907+06020.2090 ± 0.00260.282 ± 0.008Different choice of radio reference phase
J2032+41270.0990 ± 0.00070.1626 ± 0.0009Qualitatively similar results
J2047+10530.567 ± 0.010 Inconsistent results
J2215+51350.257 ± 0.0040.573 ± 0.004Inconsistent results
J2229+61140.187 ± 0.007 Qualitatively similar results
Peak separations, Δ
J0218+42320.39 ± 0.08Only one component in 3PC
J1119−61270.204 ± 0.020Only one component in 3PC
J1124−36530.21 ± 0.040.394 ± 0.005Different component phases
J1410−61320.46 ± 0.04Only one component in 3PC
J1459−60530.063 ± 0.034Only one component in 3PC
J1747−40360.681 ± 0.033Only one component in 3PC
J1823−3021A0.627 ± 0.0100.390 ± 0.005Components inverted
J1835−11060.421 ± 0.011Only one component in 3PC
J2021+40260.687 ± 0.0090.553 ± 0.005Different component phases
J2032+41270.5160 ± 0.00100.4831 ± 0.0008Components inverted
J2229+61140.299 ± 0.008 Qualitatively similar results
J2241−52360.638 ± 0.0310.317 ± 0.024Different component phases

Note. The columns list the pulsar names, 2PC and 3PC values, and a comment on the apparent discrepancy. Pulsars listed here are those for which δ and/or Δ differ between 2PC and 3PC by more than three times the statistical uncertainty, but an absolute difference >0.03.

Download table as:  ASCIITypeset image

We further compared our radio/gamma-ray alignment results to those in the literature but not in 2PC, finding discrepant radio and gamma-ray alignments for PSRs J1341−6220, J1431−4715, J1646−4346, J1731−4744, J1816+4510, J1921+1929, J1935+2025, and J2039−3616. We stress that the δ values listed in Tables 9 and 10, and the relative radio/gamma-ray alignment results in general, should be taken with a grain of salt due to the numerous systematic effects. In the case of PSR J0318+0253, whose δ-parameter is not listed in Table 10, varying rotational phase references during the commissioning of the FAST telescope prevent us from determining the alignment properly. Different choices for the radio fiducial phase, different central frequencies or frequency bandwidths used for the radio observations, or ill-determined (or varying) DM values can also lead to inconsistent δ estimates between different analyses, to name but a few examples. Nevertheless, we have established that for almost all pulsars the relative radio and gamma-ray alignments are qualitatively consistent with published results.

5.3. Energy-resolved Light-curve Fitting

The best-fit profiles above 0.1 GeV were then used to fit gamma-ray pulse profiles in narrower energy ranges. To this end, we divided the LAT data sets into logarithmically spaced energy intervals, and fitted the pulsar pulse profiles in each energy interval independently. The number of energy intervals per decade was determined by the H-test TS above 0.1 GeV: one per decade for pulsars with H-test TS values between 100 and 1000, two per decade when 1000 ≤ H ≤ 10,000, and so on. Profile fits are again done only if H-test > 100 in an energy band: 167 pulsars have at least two such intervals. Profile parameters were allowed to vary within the same ranges as for the fits over the entire energy range. By default, the number of components for a given pulsar was the same as for the fit over the entire energy range. However, profile fits in the narrow energy ranges often require fewer components than for the entire range. We fitted the profiles with Npeaks > 2 components in the individual energy ranges with all possible component combinations, and selected the combination favored by the BIC in most energy intervals. For most objects, the number of peaks remained the same; for the others, the third component (typically, bridge emission) was discarded. PSR J0007+7303, shown in Appendix B, is an example. The peak shapes, Gaussian or Lorentzian, are the same in the sub-bands as in the energy-integrated profile.

Figure 12 summarizes the results for one pulsar. The component amplitudes, phases, and widths vary with energy, as sometimes seen in other pulsars: see, for example, the phase-resolved analyses in Abdo et al. (2010f, 2010g, 2010h). In this Figure, results are shown in energy bins in which the pulsar has an H-test TS greater than 300. Peaks with amplitudes smaller than 1/100 of the maximum peak value in that energy bin, with areas smaller than 5 × 10−3 of the total area of the reconstructed pulsed emission and with widths larger than 20 times the median value, are not plotted and not included in the fit described in the following paragraph. As for the pulse profiles described in Section 5.1, all summary plots from the energy-resolved light-curve fitting are in Figure Set 12.

Figure 12. Refer to the following caption and surrounding text.

Figure 12.

Summary of the energy-resolved light-curve fitting analysis, for PSR J0205+6449. From top to bottom: ratio of the amplitude of the ith pulse profile component to that of the first component as a function of energy; differences between the component phases in the individual energy ranges and those from the global light-curve fit (for energies above 0.1 GeV); component widths as a function of energy. PSR J0205+6449 has two pulse profile components; therefore, each panel has two sets of points, with different colors and symbols. Results are shown in energy bins in which the pulsar is detected with an H-test TS greater than 300. Fit results for components with peak values smaller than 1/100 of the maximum peak value in the energy bin are not shown. Likewise, peaks with areas smaller than 5 × 10−3 of the area of the total pulsed emission and those with widths larger than 20 times the median value are also not plotted. Horizontal bands in the top and bottom panels indicate the reference results from the global light-curve fit above 0.1 GeV. The middle panel lists the reference results. The dashed curve in the top panel represents a fit of vs. energy, with the parameters shown in the bottom-left corner. (The complete figure set (129 images) is available.)

Standard image High-resolution image

We searched for trends in the amplitude ratios between the first and second gamma-ray peaks as a function of energy, for pulsars with at least two identifiable components. Assuming that the first two components have spectra that follow exponentially cutoff power laws, the amplitude ratio can be written as follows:

where C is an energy-independent constant, and Γi and Ec,i , respectively, denote the spectral index and cutoff energy of the ith component. The top panel of Figure 12 shows the fit using the above expression, and the best-fit parameters.

Figures 13 and 14 show the values of the Γ2 − Γ1 and parameters, for pulsars with at least six energy intervals with individual H-test TS > 300. Requiring ≥ 6 data points gives the fit of Equation (12) enough data points and a wide energy range. The minimum detection significance in the energy bins restricts the study to high-quality profile fits. These criteria select 28 pulsars. Figure 12.2 shows that the measured amplitude ratios between the first and second peaks for PSR J0007+7303 are poorly fit; we hence discarded this pulsar, reducing the sample to 27 objects. The choice of BLC in Figures 13 and 14 is arbitrary; we saw little or no correlation with various abscissa parameters, as discussed further in Section 8. Table 12 summarizes the results for the 27 pulsars. The weighted averages of Γ2 − Γ1 and are, respectively, −0.209(2) and −0.044(1) GeV−1, clearly incompatible with 0 but dominated by a few measurements with very small statistical uncertainties. The respective unweighted averages are −0.14(18) and −0.04(7) GeV−1, where the uncertainties here correspond to the standard deviations on the measured values.

Figure 13. Refer to the following caption and surrounding text.

Figure 13. Γ2 − Γ1 from the energy-resolved light-curve fitting, plotted against light cylinder magnetic field, BLC. See Section 5.3 for details regarding the measurement of this parameter. See Figure 10 for the definition of the circle and square symbols.

Standard image High-resolution image
Figure 14. Refer to the following caption and surrounding text.

Figure 14.  from the energy-resolved light-curve fitting, plotted against light cylinder magnetic field, BLC. See Section 5.3 for details regarding the measurement of this parameter. See Figure 10 for the definition of the circle and square symbols.

Standard image High-resolution image

Table 12. Summary of Energy-resolved Fits of the First and Second Gamma-Ray Peak Amplitude Ratios

PSRCodeCommentΓ2 − Γ1 (GeV−1)
J0205+6449rx −0.20 ± 0.06−0.12 ± 0.06
J0534+2200rx −0.105 ± 0.006−0.072 ± 0.006
J0614−3329mbrx −0.461 ± 0.032−0.023 ± 0.011
J0633+0632q 0.23 ± 0.05−0.098 ± 0.022
J0633+1746xq −0.320 ± 0.005−0.1093 ± 0.0034
J0835−4510rx −0.2100 ± 0.0021−0.0352 ± 0.0017
J1023−5746q −0.06 ± 0.040.057 ± 0.027
J1028−5819r −0.004 ± 0.024−0.010 ± 0.008
J1044−5737q −0.16 ± 0.12−0.05 ± 0.06
J1048−5832r −0.14 ± 0.08−0.104 ± 0.034
J1057−5226rxc−0.04 ± 0.11−0.12 ± 0.06
J1124−5916rx −0.13 ± 0.07−0.08 ± 0.06
J1231−1411mbrxc−0.39 ± 0.040.024 ± 0.012
J1413−6205q −0.05 ± 0.080.011 ± 0.020
J1418−6058q −0.12 ± 0.06−0.00 ± 0.04
J1709−4429rx −0.11 ± 0.08−0.06 ± 0.10
J1732−3131r −0.10 ± 0.05−0.104 ± 0.015
J1809−2332q −0.194 ± 0.0280.018 ± 0.008
J1826−1256xq −0.060 ± 0.0250.012 ± 0.011
J1833−1034r −0.02 ± 0.17−0.21 ± 0.14
J1836+5925q 0.19 ± 0.06−0.162 ± 0.028
J1907+0602r −0.23 ± 0.070.111 ± 0.030
J1952+3252rx −0.381 ± 0.0330.029 ± 0.018
J2021+3651r 0.012 ± 0.016−0.047 ± 0.007
J2032+4127bru−0.58 ± 0.070.029 ± 0.017
J2229+6114rx −0.14 ± 0.04−0.014 ± 0.016
J2238+5903q 0.06 ± 0.06−0.019 ± 0.027

Note. In Column 3, "u" indicates that the ordering of the first and second gamma-ray peaks is uncertain, and "c" indicates that the profile is complex. See Section 5.3 for details on the measurement of the parameters listed in Columns 4 and 5.

Download table as:  ASCIITypeset image

6. Phase-averaged Energy Spectra

Two significant improvements enhance our pulsar spectral results relative to 2PC. The first is simply that we have more pulsars, and more photons per pulsar. The second is a change in our expression for the exponentially cutoff power-law spectrum, a characteristic of pulsars in the GeV domain. We replace the classic subexponentially cutoff power-law model ("PLEC") with the new "PLEC4" expression introduced in 4FGL-DR3. Degeneracies between parameters in PLEC4 are reduced, and the likelihood fits converge better, yielding smaller uncertainties. However, some subtleties of the new expression merit the detailed discussion below.

In this Section we first review the motivations behind the exponentially cutoff power-law spectrum, highlighting the differences between PLEC used in 2PC and PLEC4 used here. We then describe our results, using the 4FGL-DR3 phase-integrated spectral analysis available for 255 pulsars, and explain some refinements of the spectral analyses that help characterize the pulsars' emission properties, finally yielding 255 pulsars with spectral analyses in this catalog. Tables 13 and 14 list parameter values. Section 6.5 highlights some features of individual pulsar results. Section 6.6 addresses pulsar flux stability, and in Section 6.7 we combine integrated energy fluxes and distances to obtain gamma-ray luminosities and efficiencies.

Table 13. Spectral Fitting Results for Young LAT-detected Pulsars

PSR a Energy FluxΓ100 Ep dp bfree LuminosityEfficiency b
 (erg cm−2 s−1) (GeV)  (1033 erg s−1)(%)
 (×10−11)      
J0002+62161.9 ± 0.141.6 ± 0.300.82 ± 0.160.65 ± 0.090.85 ± 0.399.01 ± 0.69 ± 3.25.9 ± 0.4 ± 3
J0007+730343.0 ± 0.331.2 ± 0.042.44 ± 0.030.36 ± 0.010.56 ± 0.04101.23 ± 0.76 ± 3922.7 ± 0.2 ± 10
J0106+48551.9 ± 0.071.2 ± 0.202.37 ± 0.110.67 ± 0.061.05 ± 0.26 73 ± 3 ± 70
J0117+5914 †
J0139+5814 †
J0205+64496.4 ± 0.201.6 ± 0.220.26 ± 0.100.45 ± 0.080.20 ± 0.2678.8 ± 2.5 ± 9.50.295 ± 0.009 ± 0.04
J0248+60212.9 ± 0.211.7 ± 0.260.77 ± 0.101.02 ± 0.121.16 ± 0.3814.1 ± 1.0 ± 2.76.6 ± 0.5 ± 1
J0357+32056.0 ± 0.180.73 ± 0.230.87 ± 0.030.95 ± 0.050.64 ± 0.14
J0359+54142.0 ± 0.131.0 ± 0.561.31 ± 0.120.86 ± 0.100.59 ± 0.40
J0514−44080.56 ± 0.040.60 ± 0.061.23 ± 0.19 25 ± 2 ± 20
J0534+2200150.0 ± 1.32.0 ± 0.030.04 ± 0.030.23 ± 0.010.61 ± 0.08
J0540−69192.7 ± 0.150.02 ± 0.080.26 ± 0.057890 ± 440 ± 3505.4 ± 0.3 ± 0.2
J0554+31071.9 ± 0.090.90 ± 0.431.19 ± 0.071.04 ± 0.110.73 ± 0.31
J0622+37491.8 ± 0.090.65 ± 0.500.75 ± 0.041.60 ± 0.150.69 ± 0.32
J0631+06461.6 ± 0.132.65 ± 0.231.05 ± 0.17 39 ± 3 ± 40
J0631+10363.0 ± 0.141.6 ± 0.122.00 ± 0.140.60 ± 0.071.55 ± 0.39 9.4 ± 0.4 ± 9
J0633+06329.6 ± 0.241.1 ± 0.171.43 ± 0.060.60 ± 0.030.48 ± 0.12
J0633+1746420.0 ± 1.11.1 ± 0.011.63 ± 0.010.66 ± 0.000.70 ± 0.0118.23 ± 0.05 ± 9.156.1 ± 0.2 ± 40
J0659+14142.7 ± 0.091.5 ± 0.270.16 ± 0.030.86 ± 0.090.12 ± 0.290.26 ± 0.01 ± 0.050.69 ± 0.02 ± 0.2
J0729−14480.6 ± 0.110.07 ± 0.410.43 ± 0.30 1.8 ± 0.3 ± 2
J0729−1836 †
J0734−15594.6 ± 0.141.4 ± 0.220.62 ± 0.080.36 ± 0.030.29 ± 0.19
J0742−28221.5 ± 0.110.95 ± 0.190.55 ± 0.11
J0744−25250.46 ± 0.121.6 ± 0.301.45 ± 0.160.64 ± 0.101.23 ± 0.58
J0802−56130.66 ± 0.061.26 ± 0.121.10 ± 0.20
J0834−4159 †
J0835−4510930.0 ± 3.01.4 ± 0.011.27 ± 0.010.52 ± 0.000.49 ± 0.01
J0908−49132.2 ± 0.310.60 ± 0.151.07 ± 0.21
J0922+0638 ⋆0.22 ± 0.040.12 ± 0.16
J0940−54281.7 ± 0.151.23 ± 0.180.85 ± 0.15 0.015 ± 0.001 ± 0.01
J1016−58577.0 ± 0.641.5 ± 0.580.89 ± 0.300.63 ± 0.100.58 ± 0.39 3.3 ± 0.3 ± 3
J1019−57492.7 ± 0.590.60 ± 0.740.50 ± 0.16 210 ± 50 ± 200
J1023−574615.0 ± 1.10.62 ± 0.130.67 ± 0.06
J1028−581925.0 ± 0.671.1 ± 0.180.89 ± 0.070.39 ± 0.020.21 ± 0.10 8.1 ± 0.2 ± 8
J1044−573711.0 ± 0.351.8 ± 0.110.65 ± 0.080.49 ± 0.030.83 ± 0.18
J1048−583219.0 ± 0.461.2 ± 0.150.97 ± 0.070.41 ± 0.020.28 ± 0.10
J1055−60281.9 ± 0.331.22 ± 0.490.69 ± 0.20 2.9 ± 0.5 ± 3
J1057−522630.0 ± 0.261.0 ± 0.061.21 ± 0.010.88 ± 0.020.85 ± 0.05
J1057−58511.3 ± 0.161.15 ± 0.131.73 ± 0.36
J1105−60373.1 ± 0.351.14 ± 0.270.76 ± 0.13
J1105−61072.3 ± 0.390.52 ± 0.370.69 ± 0.19 0.6 ± 0.1 ± 0.6
J1111−60395.9 ± 0.340.84 ± 0.532.68 ± 0.250.58 ± 0.060.30 ± 0.24
J1112−61032.2 ± 0.420.55 ± 0.540.55 ± 0.18 1.2 ± 0.2 ± 1
J1119−61274.4 ± 0.32−0.64 ± 2.71.08 ± 0.310.72 ± 0.18−0.26 ± 0.51371 ± 27 ± 3516 ± 1 ± 2
J1124−59166.1 ± 0.201.5 ± 0.220.36 ± 0.090.41 ± 0.040.23 ± 0.20168.3 ± 5.6 ± 741.41 ± 0.05 ± 0.4
J1135−60554.6 ± 0.221.2 ± 0.380.78 ± 0.130.51 ± 0.060.24 ± 0.25
J1139−62471.9 ± 0.211.92 ± 0.540.39 ± 0.09
J1151−61081.0 ± 0.141.39 ± 0.280.90 ± 0.23 1.6 ± 0.2 ± 2
J1203−62423.7 ± 0.361.10 ± 0.260.58 ± 0.08
J1208−62383.7 ± 0.401.00 ± 0.410.49 ± 0.09
J1224−6407 †
J1231−51131.1 ± 0.070.60 ± 0.061.10 ± 0.16
J1231−65111.2 ± 0.140.99 ± 0.190.82 ± 0.18
J1253−58200.43 ± 0.100.51 ± 0.161.38 ± 0.57 28 ± 7 ± 30
J1341−6220 ⋆2.2 ± 0.520.20 ± 0.13 31 ± 7 ± 30
J1350−62253.6 ± 0.411.95 ± 0.330.75 ± 0.13
J1357−64292.9 ± 0.240.35 ± 0.120.70 ± 0.12 1.11 ± 0.09 ± 1
J1358−60253.2 ± 0.360.72 ± 0.300.52 ± 0.11
J1410−61322.3 ± 0.710.75 ± 0.431.08 ± 0.36 5 ± 2 ± 5
J1413−620518.0 ± 0.641.0 ± 0.241.57 ± 0.100.50 ± 0.030.34 ± 0.13
J1418−605830.0 ± 1.81.3 ± 0.201.00 ± 0.130.76 ± 0.050.41 ± 0.13
J1420−604813.0 ± 1.60.73 ± 0.180.78 ± 0.10 4.8 ± 0.6 ± 5
J1422−61385.3 ± 0.250.55 ± 0.523.66 ± 0.180.95 ± 0.090.56 ± 0.24
J1429−591111.0 ± 0.361.6 ± 0.170.54 ± 0.080.51 ± 0.040.52 ± 0.17
J1447−57572.3 ± 0.260.46 ± 0.150.78 ± 0.16
J1459−605312.0 ± 0.371.6 ± 0.160.37 ± 0.070.44 ± 0.030.41 ± 0.15
J1509−585012.0 ± 0.531.5 ± 0.181.22 ± 0.130.53 ± 0.040.65 ± 0.17 32 ± 1 ± 30
J1513−5908‡1.2 ± 0.40
J1522−57357.6 ± 0.401.3 ± 0.290.88 ± 0.110.70 ± 0.070.53 ± 0.21
J1528−58381.8 ± 0.161.69 ± 0.200.97 ± 0.16
J1531−56101.8 ± 0.301.09 ± 0.201.41 ± 0.34 1.9 ± 0.3 ± 2
J1614−50482.0 ± 0.870.97 ± 0.870.94 ± 0.38 4 ± 2 ± 4
J1615−51373.3 ± 0.302.88 ± 0.221.34 ± 0.21
J1620−492713.0 ± 0.521.3 ± 0.202.86 ± 0.140.91 ± 0.060.94 ± 0.18
J1623−50056.1 ± 0.500.55 ± 1.13.03 ± 0.400.63 ± 0.090.29 ± 0.35
J1624−40412.7 ± 0.151.4 ± 0.271.27 ± 0.110.76 ± 0.090.87 ± 0.30
J1641−53171.5 ± 0.121.06 ± 0.130.81 ± 0.13
J1646−4346 †
J1648−46114.8 ± 0.441.1 ± 0.662.48 ± 0.310.81 ± 0.130.67 ± 0.45 55 ± 5 ± 50
J1650−46015.6 ± 0.591.32 ± 0.240.86 ± 0.13
J1702−41282.6 ± 0.640.87 ± 0.321.12 ± 0.30 14 ± 4 ± 10
J1705−19060.36 ± 0.080.88 ± 0.103.06 ± 1.66 3.9 ± 0.8 ± 4
J1709−4429140.0 ± 0.921.3 ± 0.041.43 ± 0.020.37 ± 0.010.43 ± 0.031114.2 ± 7.4 ± 45032.1 ± 0.2 ± 10
J1714−38309.1 ± 0.920.90 ± 0.820.86 ± 0.520.48 ± 0.12−0.13 ± 0.28
J1718−382510.0 ± 0.671.3 ± 0.400.56 ± 0.140.68 ± 0.070.27 ± 0.25 12.1 ± 0.8 ± 10
J1730−33504.4 ± 0.730.07 ± 0.240.47 ± 0.15 5.2 ± 0.9 ± 5
J1731−47440.39 ± 0.120.37 ± 0.073.44 ± 2.49
J1732−313118.0 ± 0.490.37 ± 0.312.09 ± 0.051.04 ± 0.050.47 ± 0.13 6.0 ± 0.2 ± 6
J1736−34221.7 ± 0.271.80 ± 0.390.75 ± 0.18
J1739−30232.1 ± 0.431.18 ± 0.261.34 ± 0.38 8 ± 2 ± 8
J1740+10000.36 ± 0.080.33 ± 0.141.72 ± 0.98 0.28 ± 0.06 ± 0.3
J1741−205412.0 ± 0.260.96 ± 0.150.84 ± 0.021.40 ± 0.050.86 ± 0.13
J1742−33211.5 ± 0.261.76 ± 0.331.18 ± 0.32
J1746−32395.1 ± 0.291.1 ± 0.371.36 ± 0.110.89 ± 0.090.65 ± 0.28
J1747−295816.0 ± 1.10.62 ± 0.540.53 ± 0.120.68 ± 0.06−0.05 ± 0.19430 ± 31 ± 12017 ± 1 ± 6
J1748−2815 †
J1757−2421 ⋆2.5 ± 0.570.09 ± 0.19 70 ± 20 ± 70
J1801−24513.2 ± 0.410.70 ± 0.470.46 ± 0.12 3.1 ± 0.4 ± 3
J1803−21498.9 ± 0.681.2 ± 0.481.34 ± 0.260.58 ± 0.080.31 ± 0.29
J1809−233242.0 ± 0.851.1 ± 0.121.29 ± 0.050.50 ± 0.020.34 ± 0.07
J1813−124625.0 ± 0.621.5 ± 0.120.48 ± 0.050.45 ± 0.020.45 ± 0.11
J1816−0755 †
J1817−17422.8 ± 0.451.12 ± 0.540.50 ± 0.14
J1826−125641.0 ± 1.31.6 ± 0.100.87 ± 0.070.63 ± 0.030.69 ± 0.11
J1827−14461.9 ± 0.182.34 ± 0.142.73 ± 0.44
J1828−1101 ⋆4.8 ± 0.780.38 ± 0.17 8 ± 1 ± 8
J1831−09525.4 ± 0.780.51 ± 0.101.50 ± 0.31 8 ± 1 ± 8
J1833−10348.9 ± 0.611.1 ± 0.760.31 ± 0.340.49 ± 0.15−0.25 ± 0.32179 ± 12 ± 250.53 ± 0.04 ± 0.08
J1835−1106 †
J1836+592562.0 ± 0.281.1 ± 0.031.48 ± 0.010.58 ± 0.010.65 ± 0.03
J1837−06042.8 ± 0.770.33 ± 1.070.41 ± 0.24 4 ± 1 ± 4
J1838−053712.0 ± 1.30.92 ± 0.661.28 ± 0.350.66 ± 0.090.19 ± 0.27
J1841−0524 †
J1844−03464.6 ± 0.531.64 ± 0.231.11 ± 0.18
J1846+09193.6 ± 0.111.2 ± 0.182.04 ± 0.080.86 ± 0.071.05 ± 0.22
J1846−02581.4 ± 0.4655 ± 18 ± 7.40.7 ± 0.2 ± 0.1
J1853−00040.49 ± 0.113.61 ± 0.492.95 ± 1.49 8 ± 2 ± 8
J1856+0113 †
J1857+01431.2 ± 0.262.41 ± 0.372.55 ± 0.86 7 ± 1 ± 6
J1906+07228.6 ± 0.800.58 ± 0.220.63 ± 0.09
J1907+060230.0 ± 0.761.2 ± 0.160.95 ± 0.060.48 ± 0.020.27 ± 0.10 8.4 ± 0.2 ± 8
J1913+09042.0 ± 0.341.69 ± 0.351.19 ± 0.29 13 ± 2 ± 10
J1913+10111.8 ± 0.311.65 ± 5.210.13 ± 0.09 1.6 ± 0.3 ± 2
J1925+17200.95 ± 0.271.18 ± 0.401.24 ± 0.48 3.0 ± 0.9 ± 3
J1928+1746 †
J1932+19166.5 ± 0.451.2 ± 0.400.72 ± 0.130.77 ± 0.080.41 ± 0.25
J1932+22200.52 ± 0.151.06 ± 0.530.75 ± 0.33
J1935+20252.2 ± 0.310.55 ± 0.180.99 ± 0.26 1.2 ± 0.2 ± 1
J1952+325215.0 ± 0.281.3 ± 0.120.99 ± 0.050.42 ± 0.020.34 ± 0.10
J1954+283611.0 ± 0.321.1 ± 0.281.17 ± 0.090.45 ± 0.030.27 ± 0.15
J1954+3852 †
J1957+50332.6 ± 0.080.50 ± 0.400.76 ± 0.031.03 ± 0.070.43 ± 0.20
J1958+284611.0 ± 0.290.72 ± 0.311.15 ± 0.060.60 ± 0.040.21 ± 0.13
J2006+31021.1 ± 0.151.47 ± 0.360.79 ± 0.21 21 ± 3 ± 20
J2017+36257.5 ± 0.461.4 ± 0.300.95 ± 0.100.92 ± 0.090.74 ± 0.25
J2021+365149.0 ± 0.790.99 ± 0.101.06 ± 0.030.50 ± 0.010.25 ± 0.06
J2021+402680.0 ± 1.11.1 ± 0.070.87 ± 0.020.92 ± 0.020.39 ± 0.05215.6 ± 2.9 ± 110183 ± 2 ± 100
J2022+38422.7 ± 0.420.20 ± 0.220.61 ± 0.15
J2028+33325.7 ± 0.19−0.33 ± 0.601.51 ± 0.061.08 ± 0.070.22 ± 0.17
J2030+36414.2 ± 0.170.55 ± 0.491.87 ± 0.081.03 ± 0.090.63 ± 0.2545.5 ± 1.9 ± 25142 ± 6 ± 100
J2030+44154.4 ± 0.220.48 ± 0.571.20 ± 0.071.20 ± 0.110.51 ± 0.27
J2032+412714.0 ± 0.471.0 ± 0.212.61 ± 0.100.56 ± 0.030.55 ± 0.1452.8 ± 1.7 ± 3.034 ± 1 ± 2
J2043+27400.9 ± 0.070.83 ± 0.101.06 ± 0.18 4.3 ± 0.3 ± 4
J2055+25395.3 ± 0.100.69 ± 0.171.21 ± 0.030.89 ± 0.040.65 ± 0.12
J2111+46064.9 ± 0.151.6 ± 0.151.33 ± 0.120.38 ± 0.030.71 ± 0.21
J2139+47162.7 ± 0.101.1 ± 0.260.98 ± 0.051.17 ± 0.090.98 ± 0.25
J2208+40560.46 ± 0.110.32 ± 0.131.72 ± 0.75 39 ± 9 ± 40
J2229+611424.0 ± 0.341.5 ± 0.070.74 ± 0.050.29 ± 0.010.31 ± 0.07 1.20 ± 0.02 ± 1
J2238+59036.6 ± 0.241.3 ± 0.230.67 ± 0.090.47 ± 0.040.30 ± 0.17
J2240+58320.98 ± 0.161.14 ± 0.470.57 ± 0.17 28 ± 4 ± 30

Notes. Unbinned maximum likelihood spectral fit results for the young LAT gamma-ray pulsars, using the PLEC4 model (Equation (15) in Section 6). Column 2 lists the phase-averaged integral energy flux in the 0.1–300 GeV energy band, G100. Column 3 lists Γ100, the photon index at 100 MeV (Equation (23)), reliable only for fits where the b parameter is left free, and not shown when fixed to b = 2/3. Columns 4 and 5 list the energy Ep of the maximum of the SED (Equation (21)), and the SED curvature dp at the maximum (Equation (22)). Column 6 lists the fitted value of b for pulsars with TS > 1000 (fixed to b = 2/3 otherwise). Tabulated values are from the b-free fit when available, and the b = 2/3 fit otherwise. Columns 7 and 8 give the total gamma-ray luminosity Lγ in the 0.1–100 GeV energy band, and the gamma-ray conversion efficiency , assuming fΩ = 1 as described in Section 6.7. The first uncertainties in Lγ and η come from the statistical uncertainties in the spectral fit, whereas the second are due to the distance uncertainty. The strong dependence of these quantities on distance (see Table 7) and beaming factor means that these values should be considered with care.

a A dagger (†) means the point source has TS < 25; hence, we show no spectral fit results. A star (⋆) means that the SED slope decreases over the full LAT energy range; hence, Ep is undefined. PSR J1513−5908 (‡ symbol): the pulsar wind nebula (PWN) MSH 15−52 dominates the colocated 4FGL-DR3 source. We use the on-pulse integral energy flux measurement of Kuiper & Hermsen (2015). b Overestimated distances or the assumed beaming factor, fΩ = 1, can result in an efficiency > 100%.

Download table as:  ASCIITypeset images: 1 2 3

Table 14. Spectral Fitting Results for LAT-detected Millisecond Pulsars

PSR a Energy FluxΓ100 Ep dp bfree LuminosityEfficiency b
 (erg cm−2 s−1) (GeV)  (1033 erg s−1)(%)
 (×10−12)      
J0023+09237.8 ± 0.61.28 ± 0.150.54 ± 0.11
J0030+045160.0 ± 1.10.88 ± 0.151.42 ± 0.030.58 ± 0.030.53 ± 0.110.743 ± 0.013 ± 0.02621.3 ± 0.4 ± 0.8
J0034−053420.0 ± 0.61.5 ± 0.160.99 ± 0.090.32 ± 0.030.65 ± 0.23 14.7 ± 0.5 ± 10
J0101−642213.0 ± 0.40.87 ± 0.311.34 ± 0.070.53 ± 0.050.43 ± 0.20 13.3 ± 0.4 ± 10
J0102+483915.0 ± 0.71.2 ± 0.331.70 ± 0.170.41 ± 0.050.48 ± 0.25 55 ± 3 ± 50
J0154+1833 †
J0218+423248.0 ± 1.21.9 ± 0.070.33 ± 0.080.26 ± 0.020.76 ± 0.17
J0248+42301.9 ± 0.31.80 ± 0.261.39 ± 0.41 2.1 ± 0.3 ± 2
J0251+26064.9 ± 0.41.57 ± 0.180.78 ± 0.17 4.4 ± 0.4 ± 4
J0307+744316.0 ± 0.61.2 ± 0.221.63 ± 0.070.99 ± 0.081.12 ± 0.28 1.35 ± 0.05 ± 1
J0312−09215.5 ± 0.41.53 ± 0.130.95 ± 0.18 2.9 ± 0.2 ± 3
J0318+02536.3 ± 0.51.89 ± 0.180.89 ± 0.16 27 ± 2 ± 30
J0340+413019.0 ± 0.70.76 ± 0.352.91 ± 0.140.56 ± 0.050.55 ± 0.23 77 ± 3 ± 70
J0418+663510.0 ± 0.71.3 ± 0.411.80 ± 0.230.42 ± 0.070.55 ± 0.34 28 ± 2 ± 30
J0437−471518.0 ± 0.51.6 ± 0.110.70 ± 0.030.74 ± 0.061.37 ± 0.240.052 ± 0.001 ± 0.0060.44 ± 0.01 ± 0.06
J0533+67599.3 ± 0.51.2 ± 0.272.29 ± 0.170.53 ± 0.070.80 ± 0.30 108 ± 5 ± 100
J0605+37576.4 ± 0.51.95 ± 0.131.39 ± 0.23 0.39 ± 0.03 ± 0.4
J0610−21007.2 ± 0.51.34 ± 0.120.79 ± 0.124.19 ± 0.27 ± 2.250 ± 3 ± 40
J0613−020038.0 ± 1.11.6 ± 0.161.16 ± 0.070.43 ± 0.030.93 ± 0.272.777 ± 0.078 ± 0.5121.1 ± 0.6 ± 4
J0614−3329110.0 ± 1.40.86 ± 0.092.84 ± 0.060.33 ± 0.010.37 ± 0.075.443 ± 0.069 ± 1.624.7 ± 0.3 ± 8
J0621+25144.1 ± 0.42.69 ± 0.241.83 ± 0.42 2.7 ± 0.3 ± 3
J0636+5128 †
J0653+47062.0 ± 0.31.24 ± 0.151.82 ± 0.58 2.6 ± 0.4 ± 2
J0737−3039A †
J0740+66202.8 ± 0.31.88 ± 0.290.59 ± 0.170.448 ± 0.051 ± 0.122.2 ± 0.3 ± 0.7
J0751+180710.0 ± 0.50.07 ± 1.22.31 ± 0.121.22 ± 0.150.50 ± 0.531.483 ± 0.081 ± 0.4920 ± 1 ± 8
J0931−19021.4 ± 0.21.89 ± 0.192.89 ± 1.16 170 ± 20 ± 200
J0952−06072.4 ± 0.31.83 ± 0.231.29 ± 0.3811.0 ± 1.5 ± 1.417 ± 2 ± 2
J0955−61507.7 ± 0.80.84 ± 0.280.49 ± 0.13 6.2 ± 0.6 ± 6
J1012−42355.6 ± 0.51.30 ± 0.400.40 ± 0.11 1.1 ± 0.1 ± 1
J1024−07194.3 ± 0.41.66 ± 0.240.65 ± 0.160.767 ± 0.075 ± 0.1914 ± 1 ± 4
J1035−672019.0 ± 0.70.55 ± 0.441.65 ± 0.080.87 ± 0.080.46 ± 0.22 6.2 ± 0.2 ± 6
J1036−83174.1 ± 0.42.15 ± 0.250.80 ± 0.19
J1048+23394.9 ± 0.50.89 ± 0.330.31 ± 0.09
J1124−365312.0 ± 0.61.4 ± 0.081.96 ± 0.150.51 ± 0.071.01 ± 0.17 8.6 ± 0.4 ± 8
J1125−58256.5 ± 0.73.02 ± 0.360.93 ± 0.21 2.9 ± 0.3 ± 3
J1125−60143.1 ± 0.91.87 ± 0.511.52 ± 0.69 4 ± 1 ± 4
J1137+75280.6 ± 0.13.20 ± 0.641.72 ± 0.78 13 ± 3 ± 10
J1142+01196.4 ± 0.52.86 ± 0.300.46 ± 0.11 80 ± 6 ± 80
J1207−50504.8 ± 0.51.75 ± 0.201.26 ± 0.28 43 ± 5 ± 40
J1221−06335.8 ± 0.51.48 ± 0.220.48 ± 0.11 1.8 ± 0.2 ± 2
J1227−485319.0 ± 1.61.7 ± 0.450.66 ± 0.350.34 ± 0.070.59 ± 0.53
J1231−1411100.0 ± 1.30.69 ± 0.112.02 ± 0.030.57 ± 0.020.53 ± 0.082.133 ± 0.027 ± 0.8911.9 ± 0.2 ± 6
J1301+08337.7 ± 0.51.30 ± 0.130.59 ± 0.102.88 ± 0.18 ± 0.354.3 ± 0.3 ± 0.5
J1302−325811.0 ± 0.50.78 ± 0.432.54 ± 0.140.84 ± 0.100.86 ± 0.36 55 ± 3 ± 50
J1311−343061.0 ± 1.21.6 ± 0.100.85 ± 0.070.24 ± 0.020.54 ± 0.15
J1312+005115.0 ± 0.71.1 ± 0.271.46 ± 0.080.84 ± 0.090.90 ± 0.30 41 ± 2 ± 40
J1327−0755 †
J1335−56568.2 ± 0.92.09 ± 0.350.52 ± 0.12
J1400−14316.4 ± 0.51.43 ± 0.140.74 ± 0.13
J1431−47154.7 ± 0.70.44 ± 0.240.80 ± 0.30
J1446−47017.7 ± 0.71.38 ± 0.290.47 ± 0.11 6.2 ± 0.6 ± 6
J1455−33301.1 ± 0.31.47 ± 0.381.34 ± 0.66 3.4 ± 0.9 ± 3
J1513−25507.6 ± 0.61.05 ± 0.200.55 ± 0.12 16 ± 1 ± 20
J1514−494642.0 ± 1.21.3 ± 0.162.55 ± 0.110.47 ± 0.040.84 ± 0.19 25.9 ± 0.7 ± 20
J1536−494880.0 ± 1.91.4 ± 0.122.14 ± 0.110.29 ± 0.020.56 ± 0.14 31.9 ± 0.7 ± 30
J1543−514918.0 ± 2.20.64 ± 0.290.58 ± 0.13 3.8 ± 0.5 ± 4
J1544+49372.4 ± 0.32.16 ± 0.390.45 ± 0.14
J1552+54372.7 ± 0.31.71 ± 0.310.50 ± 0.14 29 ± 3 ± 30
J1555−29084.7 ± 0.60.32 ± 0.360.43 ± 0.2114.5 ± 1.9 ± 3.74.7 ± 0.6 ± 1.0
J1600−30537.7 ± 0.52.74 ± 0.230.82 ± 0.14
J1614−223026.0 ± 0.8−0.16 ± 0.562.26 ± 0.080.91 ± 0.070.30 ± 0.18
J1622−03157.2 ± 0.71.51 ± 0.280.48 ± 0.11
J1625−002121.0 ± 0.71.1 ± 0.182.06 ± 0.080.86 ± 0.051.27 ± 0.29 6.1 ± 0.2 ± 6
J1628−320511.0 ± 0.91.07 ± 0.150.81 ± 0.14
J1630+37346.0 ± 0.51.80 ± 0.150.97 ± 0.18 8.7 ± 0.7 ± 8
J1640+22242.3 ± 0.32.09 ± 0.350.88 ± 0.270.621 ± 0.091 ± 0.1118 ± 3 ± 4
J1641+80492.0 ± 0.31.29 ± 0.241.02 ± 0.44 4.7 ± 0.7 ± 4
J1649−30129.0 ± 0.81.71 ± 0.230.71 ± 0.14
J1653−015834.0 ± 1.01.6 ± 0.150.83 ± 0.070.40 ± 0.040.73 ± 0.22
J1658−532421.0 ± 1.10.69 ± 0.610.82 ± 0.080.89 ± 0.110.32 ± 0.30 6.3 ± 0.3 ± 6
J1713+07477.3 ± 0.62.13 ± 0.200.77 ± 0.141.50 ± 0.11 ± 0.07942 ± 3 ± 2
J1730−23048.6 ± 1.70.63 ± 0.181.38 ± 0.53
J1732−50495.6 ± 0.71.34 ± 0.300.69 ± 0.19 63 ± 8 ± 60
J1741+13513.6 ± 0.51.90 ± 0.300.74 ± 0.23
J1744−113437.0 ± 1.40.47 ± 0.471.02 ± 0.051.04 ± 0.080.48 ± 0.24
J1744−761920.0 ± 0.61.0 ± 0.221.78 ± 0.060.87 ± 0.071.04 ± 0.25
J1745+10177.6 ± 0.61.80 ± 0.190.80 ± 0.14 21 ± 2 ± 20
J1747−403613.0 ± 1.10.82 ± 0.180.63 ± 0.12 70 ± 6 ± 70
J1805+06155.3 ± 0.52.05 ± 0.250.90 ± 0.20 10 ± 1 ± 10
J1810+174423.0 ± 0.91.2 ± 0.360.50 ± 0.110.38 ± 0.050.13 ± 0.23
J1811−240514.0 ± 3.50.49 ± 1.250.28 ± 0.14 20 ± 5 ± 20
J1816+451011.0 ± 0.50.98 ± 0.381.76 ± 0.130.49 ± 0.070.46 ± 0.27
J1823−3021A14.0 ± 0.91.84 ± 0.240.49 ± 0.0998.2 ± 6.5 ± 1011.9 ± 0.8 ± 1
J1824+10146.2 ± 0.62.16 ± 0.350.51 ± 0.13 190 ± 20 ± 200
J1824−2452A22.0 ± 1.60.29 ± 1.000.95 ± 0.160.57 ± 0.09−0.06 ± 0.3381.1 ± 5.7 ± 8.63.6 ± 0.3 ± 0.4
J1827−084930.0 ± 3.60.10 ± 0.510.24 ± 0.08
J1832−0836 †
J1833−38402.8 ± 0.41.25 ± 0.291.01 ± 0.34 7 ± 1 ± 6
J1835−3259B5.3 ± 0.745.0 ± 6.0 ± 7.916 ± 2 ± 4
J1843−111323.0 ± 2.20.21 ± 0.200.44 ± 0.10 7.2 ± 0.7 ± 7
J1855−14365.1 ± 0.82.35 ± 0.790.39 ± 0.14 170 ± 30 ± 200
J1858−221611.0 ± 0.70.10 ± 0.941.81 ± 0.121.07 ± 0.160.38 ± 0.38 10.3 ± 0.7 ± 10
J1901−012519.0 ± 1.71.14 ± 0.250.62 ± 0.12 20 ± 2 ± 20
J1902−510524.0 ± 1.00.88 ± 0.470.55 ± 0.120.37 ± 0.05−0.06 ± 0.24 11.1 ± 0.5 ± 10
J1903−70515.3 ± 0.42.57 ± 0.210.92 ± 0.16 6.2 ± 0.5 ± 6
J1908+21054.9 ± 0.81.77 ± 0.660.47 ± 0.18 12 ± 2 ± 10
J1909−3744 †
J1921+013711.0 ± 1.01.32 ± 0.320.50 ± 0.11 70 ± 7 ± 70
J1921+19293.5 ± 0.82.55 ± 0.681.13 ± 0.49 3.0 ± 0.7 ± 3
J1939+213421.0 ± 2.70.68 ± 0.280.67 ± 0.16
J1946+3417 †
J1946−54039.8 ± 0.51.4 ± 0.201.27 ± 0.081.01 ± 0.121.38 ± 0.38 29 ± 1 ± 30
J1959+204816.0 ± 0.90.19 ± 1.11.04 ± 0.100.82 ± 0.110.11 ± 0.41
J2006+01483.0 ± 0.43.37 ± 0.640.65 ± 0.22 17 ± 2 ± 20
J2017+060336.0 ± 1.01.1 ± 0.163.17 ± 0.120.58 ± 0.040.94 ± 0.19
J2017−16146.5 ± 0.61.60 ± 0.310.47 ± 0.12 21 ± 2 ± 20
J2034+363212.0 ± 0.83.16 ± 0.181.38 ± 0.19
J2039−36163.8 ± 0.51.18 ± 0.250.69 ± 0.21 14 ± 2 ± 10
J2039−561715.0 ± 0.61.8 ± 0.081.47 ± 0.170.30 ± 0.041.67 ± 0.43
J2042+02465.8 ± 0.60.84 ± 0.091.36 ± 0.27 4.8 ± 0.5 ± 5
J2043+171128.0 ± 0.81.3 ± 0.171.87 ± 0.100.41 ± 0.030.62 ± 0.196.57 ± 0.19 ± 1.143 ± 1 ± 10
J2047+10534.3 ± 0.62.75 ± 0.690.36 ± 0.14 38 ± 5 ± 40
J2051−08272.5 ± 0.32.33 ± 0.291.27 ± 0.36 12 ± 1 ± 10
J2052+12194.6 ± 0.60.45 ± 0.620.21 ± 0.0917.0 ± 2.1 ± 3.250 ± 6 ± 9
J2115+54487.0 ± 0.73.06 ± 0.301.31 ± 0.27 4.8 ± 0.5 ± 5
J2124−335839.0 ± 0.80.97 ± 0.121.95 ± 0.050.77 ± 0.041.09 ± 0.16
J2129−04296.8 ± 0.51.66 ± 0.170.58 ± 0.112.73 ± 0.19 ± 0.329.3 ± 0.6 ± 1
J2205+60123.7 ± 1.02.40 ± 0.522.90 ± 1.67 10 ± 3 ± 10
J2214+300033.0 ± 0.71.0 ± 0.151.76 ± 0.050.68 ± 0.040.92 ± 0.17
J2215+513518.0 ± 0.81.4 ± 0.301.45 ± 0.170.40 ± 0.050.52 ± 0.28
J2234+094410.0 ± 0.60.74 ± 0.651.87 ± 0.180.59 ± 0.100.36 ± 0.34
J2241−523625.0 ± 1.11.3 ± 0.111.90 ± 0.080.95 ± 0.071.35 ± 0.223.25 ± 0.14 ± 0.2512.5 ± 0.5 ± 1
J2256−10248.2 ± 0.51.7 ± 0.161.41 ± 0.160.43 ± 0.081.44 ± 0.52
J2302+444239.0 ± 0.91.1 ± 0.122.36 ± 0.070.56 ± 0.030.85 ± 0.14 89 ± 2 ± 90
J2310−05555.6 ± 0.41.74 ± 0.150.78 ± 0.14 15 ± 1 ± 10
J2317+14390.6 ± 0.21.58 ± 0.215.78 ± 3.76 9 ± 2 ± 3
J2339−053329.0 ± 0.81.3 ± 0.122.40 ± 0.100.37 ± 0.030.86 ± 0.19

Notes. Unbinned maximum likelihood spectral fit results for the LAT millisecond (MSP) gamma-ray pulsars, using the PLEC4 model (Equation (15) in Section 6). Column 2 lists the phase-averaged integral energy flux in the 0.1–300 GeV energy band, G100. Column 3 lists Γ100, the photon index at 100 MeV (Equation (23)), reliable only for fits where the b parameter is left free, and not shown when fixed to b = 2/3. Columns 4 and 5 list the energy Ep of the maximum of the SED (Equation (21)), and the SED curvature dp at the maximum (Equation (22)). Column 6 lists the fitted value of b for pulsars with TS > 1000 (fixed to b = 2/3 otherwise). Tabulated values are from the b-free fit when available, and the b = 2/3 fit otherwise. Columns 7 and 8 give the total gamma-ray luminosity Lγ in the 0.1–100 GeV energy band, and the gamma-ray conversion efficiency , assuming fΩ = 1 as described in Section 6.7. is not Doppler corrected. The first uncertainties in Lγ and η come from the statistical uncertainties in the spectral fit, whereas the second are due to the distance uncertainty. The strong dependence of these quantities on distance (see Table 7) and beaming factor means that these values should be considered with care.

a A dagger (†) means the point source has TS < 25; hence, we show no spectral fit results. b Overestimated distances or the assumed beaming factor, fΩ = 1, can result in an efficiency > 100%.

Download table as:  ASCIITypeset images: 1 2 3

6.1. Pulsar Spectral Shapes

As we will show, the general GeV pulsar spectrum is peaked around Ep = 1.5 GeV (see also Figure B1 in Appendix B). It is natural to relate this to the synchrotron spectrum radiated by an ultrarelativistic primary particle (Pacholczyk 1970; Dermer & Menon 2009)

Here, the asymptotic spectral index is Γsr = 2/3, and the cutoff energy Ec depends on the specific radiation mechanism, e.g., the field line curvature radius ρ with for curvature radiation within the light cylinder (Romani 1996). Farther out, for example in current sheets, the curvature of the particle trajectory may be quite different. Other emission mechanisms, such as inverse Compton, can also produce peaked spectra. In general, pulsars accelerate particles to a range of energies, and the observed spectrum is a broader superposition of such basic SED shapes. A main goal of gamma-ray spectroscopy is to invert this superposition and reveal properties of the underlying particle energy distribution. In a curvature radiation scenario, e.g., the observed Ep will roughly track the bulk maximum values of Ec , while the width of the peak and the low-energy spectrum (EEp) encode information about the range of Ec . A nearly monoenergetic distribution of Ec will produce a narrow SED peak and a low-energy S(E) ∝ E4/3, while a broad range of Ec will produce a broad gamma-ray peak and Γ0 of ∼1–2 (Romani 1996).

Most previous analyses of LAT pulsar spectra have used the PLEC model:

which is closely related to a monoenergetic curvature emission or synchrotron radiation spectrum, but with the cutoff strength modulated by the parameter b. Blends of particle energies and averaging over pulse phases both tend to broaden the peak, yielding "subexponential" shapes with b < 1 (Abdo et al. 2010g). One innovation since 2PC, adopted beginning with the 8 yr 4FGL catalog and justified therein, is to fix b = 2/3 instead of b = 1 for pulsars with lower TS, as this choice better matches the observed spectra.

A good spectral model should have the flexibility to capture the features evoked above, but in most cases, measuring them all is not possible. The LAT sensitivity peaks at ∼1 GeV and the SED slope at 100 MeV, 2 − Γ100, can only be measured for the brightest pulsars. Furthermore, it is only in the cases of sharply peaked spectra that quickly transition to a low-energy power law that Γ100 may be a good proxy for the asymptotic low-energy slope Γ0. Above 1 GeV, pulsar spectra cut off quickly, making it difficult to measure the cutoff shape and strength. We show below that for most pulsars, the only reliable features that can be measured are the position and width of the SED peak. Fortunately, it is these features that strongly correlate with underlying pulsar properties like , so they are likely to be the most important in capturing information about the pulsar emission mechanism.

The PLEC formula suffers from large covariance between the model parameters (Abdollahi et al. 2022), magnifying the effect of statistical uncertainty. To reduce these correlations, the spectral analysis of 4FGL-DR3 adopted the PLEC4 model (PLSuperExpCutoff4 in the Fermitools, 121 specifically designed to minimize the covariance between the model parameters at a reference energy, E0):

For PLEC4, N0 is the flux density at E0. The PLEC4 spectral shape is exactly the same as for the PLEC model, as can be seen by identifying

In the fits presented below, the reference energies E0 are confined to a relatively narrow energy range with a typical value of 1.6 GeV.

The PLEC4 model parameters Γ and d are related to the logarithmic derivatives of the model. The gamma-ray convention defines the photon spectral index as

from which it clear that Γ is the photon index at the reference energy (x = 1), and Γ0 is the asymptotic spectral index (x → 0). The spectral curvature can be defined as , so d is the SED curvature at E0.

Figure 15 shows some example realizations of the PLEC4 model and in particular highlights the role of d in producing narrow (large d) or broad (small d) peaks. b controls the overall shape of the spectrum, with larger values of b producing greater asymmetry. For b ≪ 1, the spectrum becomes more symmetric, with both a slower cutoff and a slower convergence to a power law at low energies. b > 1 produces asymmetric spectra with sharp, "superexponential" cutoffs and power-law behavior at low energy. When b < 0, the cutoff and power-law regimes are reflected. For many of these examples, the transition to power-law behavior has not yet happened at 100 MeV. In this example, the local spectral index at 100 MeV agrees with the asymptotic index to within 0.1 for b = 0.6, but for b = 0.3 it differs by 0.3 (d = 0.2), 0.9 (d = 0.6), and 1.9 (d = 1.3).

Figure 15. Refer to the following caption and surrounding text.

Figure 15. Realizations of the PLEC4 model for three values of the peak curvature, dp = 0.2, 0.6, and 1.3, and the asymmetry parameter b = 0.3, 0.6, and 1.3. We set Ep = E0 = 1 GeV, which uniquely determines Γ and Γ0. The SEDs are normalized to 1 at Ep but offset by 50% for clarity. Each curve is labeled with the tuple (d, b, Γ0). The shaded regions show the difference between the SED and the asymptotic power law, normalized to 0.1 GeV.

Standard image High-resolution image

The energy at which the SED peaks (if the peak exists, is nonzero, and is finite) is

Note that b < 0 can yield a valid Ep if Γ0 ≥ 2. Ep is a key physical property of interest, and we can obtain a second key property, the width of the SED at its peak, by evaluating the curvature at xp :

For reference, the maximum expected spectral curvature, obtained for monoenergetic synchrotron radiation, is dp = 4/3.

We expect these quantities associated with the SED peak to be well determined and relatively independent because Ep is generally close to E0. However, Equation (22) indicates that we cannot measure the peak curvature and the asymptotic index independently when b is fixed, and since we expect dp to be the better-determined quantity, we cannot use the PLEC4 model to extrapolate the asymptotic spectral index, Γ0. Γ100 could be used as a proxy for Γ0 (when Ep ≫ 100 MeV and the peak is not too broad), but it too suffers from bias: for a typical E0 = 1.6 GeV and with b = 2/3,

This formulation all but maximizes the correlation between the model parameters Γ and d, undoing the benefits of PLEC4, and like Γ0, Γ100 will potentially be biased if b = 2/3 does not describe the true SED.

Thus, to use the PLEC4 model to estimate reliable physical properties, we need to answer two primary questions:

  • 1.  
    For bright pulsars where b can be measured, does the PLEC4 model allow for robust estimation of other physical parameters of interest, like Γ100? That is, can PLEC4 reproduce the true shapes of pulsar SEDs?
  • 2.  
    For fainter pulsars, for which b must be fixed to a canonical value, does the reduced parameter space introduce a bias on physical parameters?

6.2. Spectral Model Fits

We now address these questions using phase-averaged spectral fits. All spectral results in this section rely on the 4FGL-DR3 spectral fitting methodology. Some results are incorporated directly from that catalog, and we have extended these results using the DR3 model as a baseline. In particular, we have

  • 1.  
    performed spectral fits in which b was free to vary for 119 pulsars with TS > 1000 (see 28 pulsars with TS > 10,000 in 4FGL-DR3).
  • 2.  
    fit the 28 b-free pulsars of 4FGL with a b = 2/3 model for comparison.
  • 3.  
    re-fit PSRs J0922+0638, J1757−2421, J1913+1011, and J2205+6012 with a PLEC4 model, as opposed to the power law used in 4FGL-DR3.
  • 4.  
    added a power-law component to the spectral models for pulsars with TeV counterparts (PSRs J0007+7303, J0205+6449, J1016−5857, J1119−6127, J1714−3830, J1833−1034, J1907+0602, and J2032+4127) to encapsulate potential emission from a PWN and re-optimized the PLEC4 model. Crab, Vela, and a half-dozen other pulsars were already colocated with TeV power-law sources in 4FGL.
  • 5.  
    directly measured the low-energy spectral indices with a power-law fit of the 50–300 MeV band, for the 42 pulsars with TS > 25 in that band and TS > 1000 overall.

Of the 255 pulsars, 251 produce a spectral fit with a well-defined Ep.

Here we quickly note a possible point of confusion. The reference energy E0 in Equation (15) is used in fitting models and evaluating covariance between parameters. From these results, a "pivot energy" can be determined at which the uncertainty on the flux density normalization N0 is minimized. The results presented here require the covariance matrices from the likelihood fits and so use the reference energy, and thus parameter values will generally differ from those in the FITS files accompanying 4FGL-DR3, which have been scaled to the pivot energy. A script provided with the supplementary material shows how to correctly unpack the parameters and covariance from the electronic catalog.

6.2.1. Distribution of b

First, we consider the distribution of best-fit b values, plotted against in Figure 16. It is clear that b does not take on a constant value, so we used maximum likelihood to fit a linear model to each class and to the full population. We assumed that each measured value follows a Gaussian distribution with the measurement uncertainty, and we also included an intrinsic scatter, σi , which is added in quadrature to the measurement error.

Figure 16. Refer to the following caption and surrounding text.

Figure 16. The best-fit values of b, with 1σ uncertainties, shown as a function of . Maximum likelihood linear models are shown for each class and for the combined population. The centroid indicates the model fit, while the distance between the dashed lines indicates the intrinsic scatter (±σi , see the main text). The dark shaded region indicates the same quantities for the full population with measured values of b.

Standard image High-resolution image

We excluded the Crab pulsar from the fits. In Figures 1721 its well-measured parameters generally deviate from those of most other pulsars. With its high , including it changes the fit results substantially, such that they no longer describe most other gamma-ray pulsars well. If the Crab is atypical only in having both high and good statistics, then perhaps our simple linear fit breaks down for high . But the Crab is also atypical in being a rare pulsar with much stronger flux density in the MeV range than in the GeV range, as we further discuss in Section 8.1.

Figure 17. Refer to the following caption and surrounding text.

Figure 17. The comparison between spectral indices measured directly with a low-energy fit (y-axis) and the predicted local index at 100 MeV from b-free models.

Standard image High-resolution image

The results, depicted in Figure 16, are (1) there is strong evidence (log likelihood improves by >40) of intrinsic scatter, i.e., σi > 0, and (2) there is no evidence (log likelihood improves only by a few) for a different linear relation among the classes. Therefore, we report a single model for all pulsars, finding b0 = 0.42, mB = − 0.15, and σi = 0.15.

These values are intended as a guide only: the simple model does not fully capture the observed variations, and the fitting approximates the measurement uncertainty on b as a Gaussian rather than using the full likelihood surface. Nonetheless, it is clear that many pulsars fit with a b = 2/3 model will suffer appreciable biases. In particular, estimates of Γ0 and Γ100 for MSPs (with ∼ 1034 erg s−1) will tend to be too soft, and those of RL pulsars (with erg s−1) too hard. We analyze these biases more quantitatively below. With reference to Figure 15, we also see that the distributions of b indicate that MSP spectra hew more closely to the power-law with exponential cutoff model, while young pulsars, particularly high- RL pulsars, have substantial curvature over the full LAT band and more gradual high-energy cutoffs. This is suggestive of MSPs possessing an intrinsically narrower radiating particle distribution than do their RL and RQ counterparts; how this might connect to the smaller magnetospheres of MSPs remains to be understood.

In Appendix C, we study the change in model parameters between the 119 pulsars that have b-free and b = 2/3 fits available. In particular, we find that d and Γ agree well between the models, while Γ100 is irrecoverably biased. The physical properties dp and Ep are robustly measured in both b-free and b = 2/3 fits, and we conclude that we can make inferences on their distribution using the entire population of pulsars. On the other hand, the low-energy spectrum is unreliable, and values inferred from the b = 2/3 model could be very far from the true values.

Figure 17 establishes the converse, that the b-free PLEC4 model yields a good estimate of Γ100, at least for the small sample for which model comparison is possible. The figure shows the indices resulting from the 119 low-energy fits, of which 42 have TS > 25 in the 50–300 MeV band. The geometric mean energy of this band is 122 MeV, but the improvement of LAT sensitivity with energy, along with the use of likelihood weights, means the effective energy is greater. Nevertheless, the agreement is good for all three pulsar classes, 122 and consequently we conclude that the b-free Γ100 measurements are reliable. If additionally Ep ≫ 100 MeV and dp is not too small, Γ100 may approximate Γ0 well. We discuss this further below.

6.3. Physical Properties

We can thus make inferences about three physical properties: dp and Ep for the full population of 251 pulsars, and Γ100 for the 119 b-free pulsars. Unless otherwise noted, results from b-free fits are preferred and presented whenever available.

First, we consider the distribution of spectral curvature, dp, or equivalently, (inverse) SED peak widths, with results shown in Figure 18. dp depends strongly on , varying in width by a factor of about 3 over the main population. The narrowest peaks—primarily of MSPs, but with members of each class—saturate at about 4/3, corresponding to a monoenergetic synchrotron spectrum. Only PSR J1827−1446, a radio-quiet pulsar, lies out of family with dp = 2.6 ± 0.5 and one of the largest Ep values (see below). To quantify the relation, we fit a simple power-law model to the full population (left panel, excluding Crab), finding with an additional scatter of 17% added in quadrature.

Figure 18. Refer to the following caption and surrounding text.

Figure 18. Dependence of the spectral peak curvature dp on , showing a well-defined trend. The horizontal dashed line at 4/3 indicates the expected peak curvature for a monoenergetic curvature radiation spectrum. The left-hand panel shows results for all 251 pulsars (i.e., mixes b-free and b = 2/3 fits), while the right-hand panel shows only the 116 pulsars with a b-free fit (and a well-defined peak energy). The solid line indicates a best-fit relation to the merged sample (left-hand panel).

Standard image High-resolution image

Next, we consider the distribution of peak energies, Ep, shown in Figure 19. As with dp, there is a trend with , which we characterized as with an additional scatter of 30%. This evolution leads to a distinct difference between the pulsar classes, with median values of 1.7 GeV (MSP), 1.2 GeV (RQ), and 0.9 GeV (RL). Other population variables, such as the magnetic field at the surface or light cylinder, do not reduce the larger intrinsic scatter.

Figure 19. Refer to the following caption and surrounding text.

Figure 19. Dependence of the spectral peak energy on , with low- pulsars having a peak energy roughly twice that of high- pulsars. The samples and best-fit model are the same as in the left panel of Figure 18.

Standard image High-resolution image

dp and Ep are also correlated (Figure 20), but this relation appears to be entirely captured by evolution. Using the single-variable models determined above removes the correlation to the level of the intrinsic scatter in Ep. In general, we arrive at a picture of larger peak widths at lower peak energies, a paradigm clearly illustrated by Crab-like pulsars with very broad peaks with, presumably, substantial contributions from secondary particles. The large widths are such that some pulsars, like the Crab and Vela, emit into the TeV range (HESS Collab. et al. 2023), implying very high particle energies, which emission models need to accommodate.

Figure 20. Refer to the following caption and surrounding text.

Figure 20. Dependence of the spectral curvature dp on the peak energy, Ep. Panel selections are the same as in Figure 18.

Standard image High-resolution image

Finally, we consider Γ100, a proxy for the asymptotic spectral index. As an illustration of the possibility of making an incorrect inference about the physical properties of the pulsar population, in the first panel of Figure 21 we show the results of 251 b = 2/3 fits, revealing an apparently significant trend with hardening spectral indices for increasing . However, this result is largely due to the bias relating Γ100 to the spectral curvature and thus reproduces the trend shown in Figure 18. In the right panel, we show the b-free fits, which we demonstrated provide relatively unbiased estimators for Γ100. The trend is reduced or absent: the bulk of the pulsars with erg s−1 have a spectral index near 1.2, and somewhat lower (≈1) for lower . Below this (possible) threshold, harder spectra become apparent, with values approaching the monoenergetic limit of 2/3. These pulsars have the highest values of Ep, and so for them, Γ100 may be closer to Γ0.

Figure 21. Refer to the following caption and surrounding text.

Figure 21. Dependence of the low-energy spectral index Γ100 on . The left panel shows fits in which every pulsar has b = 2/3, revealing a strong but erroneous trend. The right panel shows only the 119 unbiased b-free fits. The horizontal dashed line indicates the monoenergetic limit, 2/3.

Standard image High-resolution image

Both observed SED shapes and pulse profile shapes should depend on the magnetic field strength and configuration, as well as on "geometry," meaning, for example, the inclination α of the magnetic axis relative to the rotation axis; the inclination ζ of the rotation axis relative to the line of sight; and the radius of the light cylinder RLC = cP/2π. We thus searched for correlations between measured profile and spectral parameters. Figure 22 suggests that pulsars with classic Vela- or Crab-like gamma-ray pulse profiles, with Δ ≲ 0.5 rotations, have typical SED peak energies. The highest Ep values seem to occur mainly for pulsars where the second gamma-ray peak follows the first more closely. This is perhaps related to the evidence that the gamma pulse farther from the magnetic pole generally has the harder spectrum (Figure 13).

Figure 22. Refer to the following caption and surrounding text.

Figure 22. Left: the SED peak energy Ep vs. the phase separation Δ between the two principal gamma-ray peaks. The Crab and PSR J1730-3350 are off-scale with Ep ≈ 50 MeV. Right: Ep distributions for pulsars with closely (top) and widely (bottom) spaced peaks. Pulsars with closely spaced peaks perhaps have the highest Ep values.

Standard image High-resolution image

6.4. Summary of Results

We summarize our main results here:

  • 1.  
    The broadband shape, governed by the parameter b, evolves with , yielding slowly decaying spectra for high- pulsars and sharply cutoff spectra for low- pulsars.
  • 2.  
    Pulsar spectral energy densities () are most effectively characterized by the curvature at the spectral energy density peak, dp, which is inversely proportional to the spectral energy density peak width. dp anticorrelates with with low (17%) scatter; see Figure 18.
  • 3.  
    The energy at which the spectral energy density peaks, Ep, also varies inversely with with somewhat larger (30%) scatter; see Figure 19.
  • 4.  
    Consequently, the spectra of MSPs are more sharply peaked than those of young pulsars, and reach the highest peak energies.
  • 5.  
    The asymptotic (E → 0) spectral index, Γ0 lies out of the LAT band. The spectral index near 100 MeV, Γ100, is close to Γ0 when Ep ≫ 100 MeV, and it can be estimated for 116 of the brightest pulsars (Figure 21). These indices appear to take on an almost universal value of about 1.2, though low- pulsars, particularly MSPs, may have harder spectra that approach the "monoenergetic limit" of 2/3.
  • 6.  
    Thus, for our LAT pulsars, nearly all of the rotational energy that is converted to electromagnetic radiation is concentrated in a single decade of the spectrum.
  • 7.  
    The spectral energy density peak energy Ep seems to be highest for pulsars with pulse phase separations Δ < 0.4 (Figure 22). This, coupled with the observation that the second pulse tends to have a harder spectrum than the first (Figure 13), is new empirical evidence that the same underlying magnetosphere configurations affect both profile and spectral shapes.

6.5. Notes on Individual Pulsars

Twenty gamma-ray pulsars in Tables 13 and 14 are flagged with † symbols, meaning that no 4FGL-DR3 source with Test Statistic > 25 is colocated with the pulsar. The pulsed signature is detectable in these cases, in spite of a low gamma-ray flux or large local background.

Guillemot et al. (2013) provided spectral parameters for the double pulsar PSR J0737−3039A via a dedicated analysis. PSR J1513−5908 (B1509−58) is brightest in the MeV domain, with a spectral cutoff below 100 MeV, measured by Kuiper et al. (2018). The colocated extended source 4FGL J1514.2−5909e corresponds to the PWN MSH 15−52, and its spectral parameters do not apply to the pulsar. Similarly, PSR J1023+0038 is a transitional pulsar, and the 4FGL-DR3 spectral results are dominated by the data obtained after it transitioned to its accreting state, where the (unpulsed) gamma-ray flux is greatly enhanced.

PSR J1748-2815 was colocated with a TS ∼ 300 source in the LAT 8 yr catalog. Prior to 4FGL-DR1, the LAT team adopted a new model for the diffuse Galactic gamma-ray emission that is brighter around the Galactic center, and the source disappeared, although its pulsations remain significant with the model weights used here.

For PSR J1531−5610, 2PC mistakenly reported a luminosity Lγ over 10 times smaller than the value in Table 13. An outlier in 2PC, with corrected Lγ , its efficiency becomes quite typical.

6.6. Searches for Pulsar Flux Variability

In 4FGL-DR3, fluxes for sources with Variability_Index < 24.7 are considered stable. Three pulsars exceed this threshold: PSR J2021+4026, cited in Section 3 with Variability_Index = 285, and the two transitional MSPs J1023+0038 and J1227−4853 with Variability_Index = 58 and 1470, respectively, which exhibit large flux changes when switching between accretion- and rotation-powered states, but appear to be stable emitters like all other gamma-ray MSPs when in the rotation-powered state. In contrast, of PSR J0540−6910 in the Large Magellanic Cloud increased by one-third, with no corresponding change in the gamma-ray flux (Marshall et al. 2015), contrary to what the correlation between Lγ and seen in Figure 23 suggests (Variability_Index = 11.1). A search for variability in our pulsar sample yielded no new cases.

Figure 23. Refer to the following caption and surrounding text.

Figure 23. Gamma-ray luminosity Lγ = 4π fΩ d2 G100 in the 0.1–100 GeV energy band vs. spindown power . The vertical error bars from the statistical uncertainty on the energy flux G100 are colored. The vertical error bars due to the distance uncertainties are black-dashed, and generally larger. Doppler corrections (Section 4.3) have been applied to MSPs with known proper motions, leading to visible horizontal error bars in many cases, and when the corrections yield negative , we place the limit at erg s−1 for clarity. Luminosity upper limits use distance upper limits, while limits result from insignificant proper-motion measurements or undetermined distances. We use a light color shade to plot these limits in order to de-emphasize them, and we do not include them in the tallies in the legend. The upper diagonal line indicates 100% conversion of spindown power into gamma-ray flux: pulsars above this line may have smaller distance d, and/or the assumed beam correction fΩ ≡ 1 is wrong. The lower diagonal line indicates the heuristic luminosity , to guide the eye. For the Crab and PSR J0540−6919 in the Large Magellanic Cloud, at far right, the upper points include the X-ray energy flux.

Standard image High-resolution image

The linear increase in pulsed significance (weighted H-test) with time of PSR J2256−1024 changes abruptly near MJD 56600. It has Variability_Index = 18.3. A dedicated analysis confirms that the pulsar flux is constant, but that nearby 4FGL J2250.0+3825, associated with a blazar, has a very bright flare peaking at the same time. The softer significance slope is thus due to increased background rather than decreased signal. Exposure changes due to pointing changes also cause H-test versus MJD slope changes.

PSR J2116+1345 is an MSP discovered using the Arecibo radio telescope in a search of a pulsar-like LAT source. Gamma-ray pulsations were found using Arecibo and GMRT timing after we froze the list of pulsars reported in this catalog. Its corresponding source, 4FGL J2117.0+1344, has Variability_Index = 16.2. Here too, pulsed significance is not linear with time. Analysis showed that the decreased rate of significance growth between MJD 56900 and 57900 is due to the flaring activity of 4FGL J2035.4+1056 (PKS 2032+107), which is 10fdg5 away from the pulsar.

Clark et al. (2023b) discovered dips in the gamma-ray flux of "spider" MSPs at orbital phases when the companion star lies between Earth and the pulsar. The 4FGL-DR3 fluxes are integrated in both rotational and orbital phase; these eclipses hence do not affect Variability_Index. Orbital modulations (see, for example, An et al. 2020, or the references cited in Table 15) similarly do not influence Variability_Index.

Table 15. Candidate "Spider" and Transitional MSPs (tMSPs) Discovered in Deep Searches of Previously Unidentified LAT Catalog Sources

LAT NameType Epeak Class, Assn NameReferences
J0212.1+5321RB1.20bin 1SXPS J021210.6+532136Li et al. (2016), Linares et al. (2017)
J0336.0+7502BW1.71lmb 4FGL J0336.0+7502Li et al. (2021)
J0407.7−5702DStMSPPL(––)Miller et al. (2020)
J0427.8−6704DStMSPPLLMB 1SXPS J042749.2-670434Strader et al. (2016)
J0523.3−2527RB1.77bin CRTS J052316.9-252737Strader et al. (2014)
J0540.0−7552DStMSP0.59(––)Strader et al. (2021)
J0846.0+2820RB  Swihart et al. (2017)
J0935.3+0901RB0.50(––)Wang et al. (2020)
J0940.3−7610RB1.13lmb 2SXPS J094023.5-761001Swihart et al. (2021)
J1120.0−2204RB+1.0(––)Swihart et al. (2022a)
J1408.6−2917BW0.18bin 4FGL J1408.6−2917Swihart et al. (2022b)
J1544.5−1126DStMSP0.36lmb 1RXS J154439.4-112820Britt et al. (2017)
J1702.7−5655RB0.54(––)Corbet et al. (2022)

Note. "BW" means black widow, "RB" means redback, and "DStMSP" means that a "disk state" is presumed to have obscured pulsed detections to date. "PL" in place of the SED peak energy (in GeV) indicates that the best spectral fit uses a power law instead of log parabola. Two tMSPs seen to pulse in gamma-rays are PSRs J1023+0038 and J1227−4853. All LAT names are 4FGL except for 2FGL J0846.0+2820, detected only in 2009.

Download table as:  ASCIITypeset image

6.7. Luminosity

Most of the low-frequency Poynting flux generated by the dipole radiation of the rotating neutron star converts to a wind of electrons and ions toward and/or beyond the light cylinder, carrying away angular momentum and energy, braking the rotation. About 10% of the known pulsar population—those in this catalog—convert much of their spindown power into intense gamma-ray beams. Pulsar energetics contribute a piece to the broader puzzle of Galactic ecology, as recently illustrated by the PeV energy detection by LHAASO of the Crab nebula (Cao et al. 2021) and by the Fermi observation of GeV-band halos around three bright HAWC sources (Di Mauro et al. 2021), and gamma-ray pulsars allow for the most direct look at the acceleration mechanisms.

Integrating the LAT spectral energy flux above 100 MeV yields G100, which is tabulated in 4FGL. The range of G100 values for the gamma-ray pulsars is apparent in Figure 6. We can then calculate the luminosity

using the distance d to the pulsar, and the beaming fraction fΩ.

Figure 23 shows that the luminosity (Equation (24)) is highly correlated with the available braking power. However, at a given , there are two decades of dispersion in the observed efficiency , shown in Figure 24. Distance errors cause some of the dispersion. The total radiated power likely depends on the inclination α between the neutron star's magnetic and rotation axes, adding more pulsar-to-pulsar variation. The factor fΩ in Equation (24) is the ratio of total beam power averaged over the whole sky to the power in the beam slice illuminating the Earth averaged in phase. It is defined in 2PC Equation (16), and by Romani & Watters (2010), and depends on both α and on the inclination ζ between the rotation axis and the line of sight. Different models predict different pulse profile shapes and thus fΩ, and in any case, observational difficulties cause large uncertainties in the α, ζ model inputs. Lacking fΩ(α, ζ), our choice of fΩ = 1 likely contributes to the observed η dispersion. A key test of a successful theory of pulsar radiation is its predictions of Lγ for given pulsars. Kalapotharakos et al. (2023) generally found fΩ < 1 for the LAT pulsar sample, which would make our Lγ values overestimates.

Figure 24. Refer to the following caption and surrounding text.

Figure 24. Gamma-ray efficiency vs. spindown power . The error bars are the same as in Figure 23. The dashed diagonal shows , corresponding to the heuristic luminosity , to guide the eye. The correlation is more striking than in 2PC.

Standard image High-resolution image

The upper diagonal line in Figure 23 shows 100% efficiency, , for converting spindown power into gamma-rays. Progress since 2PC in determining pulsar distances and MSP proper motions has eliminated pulsars more than 1σ above the line. Furthermore, RNS = 10 km is an outdated convention, and is ≈2 times larger than what we use, reducing all η values. The neutron star mass values MNS in our sample probably have variance <20% (see, for example, Figure 5 in Shamohammadi et al. 2023), contributing imperceptibly to the luminosity dispersion.

The lower diagonal shows the "heuristic" gamma-ray pulsar luminosity,

motivated by the idea that above some minimum open field-line voltage volts, electron–positron cascades occur in an acceleration region with a self-governing size that produces a gamma-ray efficiency epsilonγ ∝ 1/V (Arons 1996). Recent particle-in-cell simulations of pulsar magnetospheres also support the idea of a critical pair production rate below which pulsar magnetospheres become "dead" (e.g., Kalapotharakos et al. 2018; Philippov & Spitkovsky 2018; Chen et al. 2020). The scaling of Lγ with depends on whether gamma-rays are radiated primarily via synchrotron or curvature processes (Kalapotharakos et al. 2019). We discuss our results in the context of these simulations further in Section 8.

7. The Pulsars Not Seen

Here we address known pulsars not seen with the LAT, in spite of characteristics shared with gamma-ray emitters. PSR J1846−0258 (Section 7.2) is a unique case of a pulsar seen in LAT data only below 100 MeV (Kuiper et al. 2018), not using the methods of this catalog. Finally, we present the LAT's sensitivity for future pulsar discoveries.

7.1. Radio and X-Ray Pulsars

Figure 5 shows that about 5% and 35% of young pulsars and MSPs, respectively, are seen with the LAT in the erg s−1 decade, climbing to over 60% and 80% at high . Removing pulsars discovered using LAT sources as seeds, the low fractions decrease to about 3% and 25%. 2PC Table 13 listed 28 pulsars with erg s−1 that were not seen at the time, of the 64 then known. Since 2PC, half of these 28 have been detected and are in this catalog, while the total number increased to 71. The majority of those not yet seen are young X-ray pulsars, discussed below.

With a distance estimate, we calculate the pulsar's "heuristic" flux , which converts to the units of the spectral measurements using G100/Gh = 3.33 × 10−27 = 10−26.48. Figure 5 shows the fraction of detected pulsars rising steadily with Gh above the LAT detection threshold near G100 ≈ 10−12 erg cm−2 s−1 (detailed in Section 7.3). Here, removing pulsars discovered in searches of unidentified LAT sources moves the MSP points closer to those of the young pulsars. Smith et al. (2019) discussed whether the nondetection of lower pulsars is due to a deathline, meaning that the emission mechanisms cease below a minimum braking power, or comes from the small number of pulsars close to Earth and the low detection fraction for low Gh .

For the nondetections at high , for the young pulsars, Johnston et al. (2020) found that beam geometry and pulsar inclination dominate (the beam misses the Earth). Rookyard et al. (2017) observed that broad radio pulses are a good proxy for the inclinations unfavorable to gamma detection. An example is the young pulsar PSR J1302-6350 (B1259-63), with large and colocated with an LAT source (see Table 5), yet showing no gamma-ray pulsations (Abdo et al. 2011; Smith et al. 2019). GeV flares occur during periastron every 3.4 yr, when the neutron star moves out of the stellar disk of the companion star (Chang et al. 2018). Here, the absence of gamma-ray pulsations is likely due to the beam orientation. Some gamma-ray pulsars are detected at much lower efficiency than expected from the general trend (e.g., PSRs J0659+1414, J0940−5428 and J1740+1000 near and 190 × 1034 erg s−1, respectively, in Figure 24). If these objects are typical gamma-ray pulsars, then Earth lies well outside of the standard strong beam, and a weak fringe of this beam or an alternate emission component may be responsible for such "subluminous" gamma-ray pulsars (Romani et al. 2011).

Different factors contribute to the larger fraction of MSP gamma-ray detections. For gamma-rays emitted in the inner magnetosphere, the small light cylinders broaden the beams. However, current sheet emission models predict similar beams for MSPs and young pulsars. MSPs have lower , and thus lower Lγ than young pulsars (see Figure 23), but they are also generally closer to Earth. Figure 19 shows that Ep is slightly higher for MSPs, facilitating detection, as was highlighted by Kalapotharakos et al. (2017) using 2PC data. The details of which MSPs are not seen by the LAT remain to be clarified.

Another category of as-yet unseen pulsars is those for which the radio rotation ephemerides do not suffice for multiyear gamma folding. An example is the 33 radio MSPs found in deep radio searches targeting unidentified LAT sources that have not yet been confirmed as gamma-ray pulsars (Tables 1 and 6). It also occurs for pulsars found in a survey, but never included in a long-term timing campaign: psrcat lists 278 pulsars with 30 < P < 800 ms but without . Nearly half were recently discovered by the GPPS 123 and CRAFT 124 surveys, both on the FAST telescope (Han et al. 2021), with improved measurements ongoing. Some of the others have been added to the timing program for Fermi at the Parkes radio telescope, when their decl. and radio brightness allow. Of those found to have large , a fraction may be detected in LAT data in the coming years.

X-ray timing with RXTE enabled the gamma-ray discovery of PSR J0540−6919 in the Large Magellanic Cloud (Ackermann et al. 2015). Subsequent Swift timing revealed a rare 30% increase in and a near-zero braking index (Marshall et al. 2016). XMM-Newton timing assisted the LAT pulsation searches for PSR J2022+3842 (Limyansky 2022).

Several of the high- pulsars unseen in LAT data mentioned above are radio-quiet X-ray pulsars. Timing noise for such young pulsars is generally high, and while X-ray observations made over a few years with XMM or Chandra yield astrometry and the spindown rate, they cannot give the phase-connected ephemeris necessary for multiyear gamma-folding. Thus, for many years, only pulsars that have low timing noise and/or are bright enough for routine observations with RXTE or Swift could be gamma-ray phase-folded with the same level of confidence as we have for our 1000 radio ephemerides.

NICER has provided new capabilities, allowing precision X-ray timing for more pulsars than previously possible. An example is the study of PSR J1813−1749, the third pulsar in 2PC Table 13, by Ho et al. (2020). The pulsar has no 4FGL counterpart and gamma folding yields a null result. Ho et al. (2022) presented NICER timing for five other high- X-ray pulsars, unseen with the LAT.

X-ray and optical studies of unidentified 4FGL sources have led to the pulsar candidates listed in Table 15. Many of these have led to MSP discoveries in radio pulsation searches using the methods described in Section 3.3 or gamma-ray MSP discoveries in restricted blind searches (Section 3.4), once the orbital parameters are better constrained by further multiwavelength studies. Transitional MSPs, of which PSR J1023+0038 is the archetype (see Stappers et al. 2014, and Appendix A), may "turn on" and show pulsations at some future date. In a similar vein, Hare et al. (2019) identified a point source in the heart of an X-ray PWN that merits pulsation searches.

7.2. PSR J1846-0258

PSR J1846-0258 is a young (τ ≈ 723 yr), high magnetic field (5 × 1013 G) pulsar located in the Kes 75 SNR (Gotthelf et al. 2000). Although radio-quiet (Archibald et al. 2008), it has been timed regularly in X-rays by RXTE, Swift, INTEGRAL, and NICER. PSR J1846-0258 behaves largely as an RPP, but has twice exhibited magnetar-like outbursts (Gavriil et al. 2008; Blumer et al. 2021). These outbursts have garnered interest in this source as a "magnetar-pulsar transitional object," possibly able to shed light on the evolutionary relationship between RPPs and magnetars. Additionally, it is one of 18 "soft" gamma-ray pulsars (Kuiper & Hermsen 2015), with strong emission in the MeV range. The X-ray spectrum suggested that the LAT would be able to detect the pulsed high-energy tail, aiding research into the soft pulsar emission mechanism (Torres 2018).

PSR J1846-0258 is the only pulsar with a published gamma-ray detection (pulsed 4.2σ, via unweighted , see Equation (2) in Section 3) that is not detected using the methods of this catalog (Kuiper et al. 2018). It has no 4FGL-DR3 counterpart source, and thus no Lγ measurement. Following Kuiper et al.'s (2018) prescription, we confirmed this detection. It increases to 4.3σ using Pass 8 data. As discussed in Section 3, timing solutions may produce significance versus time curves that peak outside the ephemeris validity range. This is indeed the case with PSR J1846-0258, with a maximum H-Test () of 4.4σ (5.0σ) occurring ∼484 days after the end of the published 3295 day timing solution (including a 153 day loss of coherence). We also used NICER to construct an independent ∼770 days timing solution outside the time range covered by Kuiper et al. (2018), which we applied to similarly filtered LAT data. This did not yield an LAT detection, which, given the 4× shorter span of the NICER data, is consistent with expectations from the previous detection (Limyansky 2022).

7.3. Flux Upper Limits and Sensitivity

Discovery of most of the pulsars in this list started with an initial phase-integrated source detection of the type applied to all Fermi-LAT catalog entries. We summarize the 4FGL procedure: for an all-sky model defined by pointlike (Kerr 2010), which includes all current sources, a new trial source is introduced at the position of each of HEALPix nside = 256 (∼1/4 deg) pixels. Its flux is optimized using each of five fixed trial spectral shapes. Three are power laws, but two are peaked, at 1 and 2 GeV. The one with the highest TS, if greater than 10, is used as a seed to add to the model. The model is re-optimized with now free spectral parameters, including the new source or sources. All sources are also relocalized. All that survive with TS > 10 are used for additional iterations. Finally, the list is used as seeds for the full catalog gtlike analysis. After this stage, only sources with TS > 25 are included in the final catalog. Similarly, 2PC started from the 3FGL catalog, which however used only a single power-law trial spectrum.

Construction of an estimated sensitivity threshold map, shown in Figure 25, followed a similar procedure. For each HEALPix pixel location, a trial source was introduced, now with a pulsar-like spectrum, and the flux likelihood function determined. Instead of optimizing it as before, the value corresponding to its 95% cumulative probability was recorded for that point. A more precise procedure would have instead inserted simulated trial data with varying intensity, applying the full subsequent analysis procedure, to determine the threshold curve.

Figure 25. Refer to the following caption and surrounding text.

Figure 25. Hammer–Aitoff projection of the estimated LAT 12 yr sky-survey energy flux pulsar detection sensitivity.

Standard image High-resolution image

We assume that our procedure yields the correct position dependence. Converting the sensitivity estimates to a 50% completeness threshold requires an ad hoc scale adjustment, which we obtain by comparing the distribution of sensitivity values to the detected pulsar fluxes, shown in Figure 26 as a function of latitude. This detection threshold analysis considers only phase-integrated source detections.

Figure 26. Refer to the following caption and surrounding text.

Figure 26. Integral energy flux from 0.1−100 GeV, G100, vs. Galactic latitude b (scaled as b0.65 for clarity). The shaded area shows the 5th–95th percentile range of the 12 yr pulsar detection sensitivity in the sky directions at the corresponding latitude. The points correspond to the gamma-ray pulsars in this catalog.

Standard image High-resolution image

As discussed in 2PC and by Hou et al. (2014), and further demonstrated by the 15 pulsars in this catalog that were detected via their pulsations but have TS below the 4FGL threshold, narrow gamma-ray pulses can facilitate source discovery. Such pulsars appear in Figure 26 near the shaded area.

Finally, a requirement for inclusion in the catalog is that no new candidate source be within 0fdg5 of an existing one. Therefore, in the map in Figure 25, we set such pixels around each 4FGL source to nan. The all-sky sensitivity threshold map is provided in FITS format.

8. Discussion

8.1. Phase-averaged Spectral Trends

In Section 6, we established that the energy where the SED culminates, Ep, and the SED's curvature at that energy, dp, are robust quantities that can be reliably estimated when using the general b-free SED model or the restricted b = 2/3 case. Ep, dp, and b all evolve with , but the scatter of b and Ep is much larger compared to dp. For 90% of the pulsar sample, 0.5 GeV < Ep < 2.5 GeV, with young, high- pulsars at the bottom of the range and older, unrecycled pulsars and MSPs at the top.

A consequence of the opposed evolution of Ep and dp is a nearly universal cutoff energy, epsilon37,hi, at which the SED has fallen 1/e ≈ 0.37 from the peak (Figure 27). Another way of looking at this universality is to consider individual spectral shapes. Figure 28 shows the best-fit SED shapes for the b-free pulsars, encoded according to . The region between epsilon37,hi and epsilon37,lo (the mirror 1/e point below Ep) is highlighted, allowing for a visual estimate of the effective maximum energy and width. The upper peak bounds cluster in the 3–10 GeV range, independently of . Thus, we have the simple descriptive result: almost all γ-ray pulsars emit their power within a two-decade (0.1–10 GeV) envelope, with the Crab-like, flat-spectrum, high- pulsars nearly filling the envelope, while sharply peaked MSPs concentrate their power into only the higher decade. In contrast with this absolute energy scale, the spectral curvature as estimated by dp varies strongly with . Selection bias, if any, contributes weakly to this result, as we discuss below.

Figure 27. Refer to the following caption and surrounding text.

Figure 27. One realization of the cutoff energy, epsilon37,hi, takes on a nearly constant value, shown here in the estimates for the b-free sample.

Standard image High-resolution image
Figure 28. Refer to the following caption and surrounding text.

Figure 28. The spectral energy density, S(E), for 116 pulsars with an estimate of b. All pulsars have a valid Ep, marked with a circle (b > 0) or a cross (six pulsars with b < 0). The spectra are normalized such that . Curve colors in decades guide the eye. Curves are dashed below epsilon37,lo and end at epsilon37,lo.

Standard image High-resolution image

In Section 6, we also showed that the low-energy spectral index Γ100 could not be measured reliably for b = 2/3 models, and of the population of 116 pulsars for which it can be measured, the values tend to cluster around 1.2 with a scatter of a few tenths and only a very modest evolution with . This result is clearly visible in Figure 28, in which it can further be seen that many of these SEDs have not transitioned to a power law at 100 MeV, making these measurements effectively upper limits on the asymptotic spectral index. The low-energy spectral shape is thus another frequent property of gamma-ray pulsars, taking on values that are well in excess of the monoenergetic synchrotron radiation spectral index (Γ0 = Γsr = 2/3) for all but few of the lowest- MSPs.

As stated, the curvature evolves strongly with , and depends on the spectral steepness beyond Ep and epsilon37,hi. Indeed, Lyutikov et al. (2012) and Richards & Lyutikov (2018) argued that the shape of the high-energy spectral tail is particularly rich in information about the highest energies reached by the particle acceleration mechanisms, and thus, whether and when the radiation-reaction limit is reached. The tail's shape can reveal, for example, the presence of an inverse Compton "rebound," and also informs about the best candidates for >10 GeV emission to be observed by Cherenkov telescopes. Barnard et al. (2022) calculated the very high-energy Vela spectrum radiated by primary particles in the pulsar magnetosphere and current sheet. Measuring the high-energy spectral shape with Fermi alone is, however, challenging: the photon rates decrease rapidly, and possible contributions from a surrounding PWN and/or SNR must be handled carefully. We envisage future studies dedicated specifically to these measurements.

The low- and high-energy spectral shapes, and thus the curvature, depend on the underlying acceleration of leptons; on variations in the local conditions that govern the balance between acceleration and radiation, observational convolutions of which depend on pulse phase; and on the extent and re-acceleration of pair cascades. Disentangling these effects may not be possible with only GeV data, but the robust evolution of spectral shape demonstrated here is an important input in constraining models of pulsar emission. Future successful theories should also reproduce the approximately constant spectral indices at 100 MeV, the weak decline of Ep with pulsar , curvatures which track , and the overall shape that evolves from flat, Crab-like spectra to MSP-like spectra with well-defined power-law indices and sharp peaks. These measured trends also provide new ingredients for population syntheses: by adopting a particular spectral shape for a given synthesized , more reliable estimates of the LAT sensitivity and selection effects can be obtained.

The physical implication of these results may be somewhat diluted since the spectra are phase-averaged, with multiple distinct acceleration regions (whose observational sampling depends on the pulse phase) contributing to and broadening the measured spectral features. On the other hand, we can point to the relatively tight clustering of dp about the scaling relation, which suggests that these accelerators must have broadly similar properties. Indeed, while we established a general trend of harder trailing peaks in Section 5.3, these differences do not seem to be large enough to add appreciable scatter to the observed dp values. Furthermore, at least some MSPs approach the shape of monoenergetic curvature radiation, further limiting the possible variations in spectral shape with phase. On the other hand, the preference of b < 1 for high- pulsars may indicate a richer dependence on phase.

Kalapotharakos et al. (2019) examined the mechanisms and configurations at play in the current sheet just outside the light cylinder, and conclude that curvature radiation dominates over synchrotron radiation yielding a "fundamental plane" (FP) such that . The spectral cutoff energy epsiloncut tracks the pulsar's highest photon energies, which reflect the highest attained particle energies. Kalapotharakos et al. (2022) developed a proxy for epsiloncut called epsilonc1, obtained by imposing a PLEC shape with b = 1 on the 4FGL spectra and then finding the energy epsilonc1 > Ep where the SED falls to 10% of its maximum. This is similar to our epsilon37,hi shown in Figure 27, and we include PLEC_E10_b23 and PLEC_E10_bfr in the 3PC catalog FITS file, with the important difference that we use the b = 2/3 and b-free PLEC4 fits that give the best LAT photon data likelihoods. Kalapotharakos et al. (2022, 2023) thus improved the FP data fits as compared to their earlier work, finding good agreement with their FP exponent predictions. Further improvement could include selecting pulsars with the most reliable distances, and reviewing the choice of the epsiloncut estimator in light of what we established here regarding measurement bias, and improving the uncertainties used in the FP fits. Applying fΩ corrections to Lγ for pulsars with reliable estimates of the inclination angles α, ζ, including Shklovskii corrections for the MSPs, and perhaps correcting for pulsars with neutron star mass measurements may also improve FP tests.

8.2. Selection Bias

Are our results representative of all gamma-ray pulsars? Here we address whether observation bias might skew our Ep distribution, and prospects for discovering pulsars with lower Ep. In 4FGL-DR3, over 2000 sources have Ep < 500 MeV. They were fit with the LogParabola function rather than PLEC4, but that affects the Ep values only marginally. The LAT detection and analysis chain is clearly sensitive to SEDs with maxima below what we see in our sample.

Half of our gamma-ray pulsars were discovered using radio rotation ephemerides of pulsars discovered independently of the LAT, code "R," 125 generally in radio surveys covering large sky areas, with selection biases a priori unrelated to gamma-ray properties. Conceivably, neutron star orientations such that both radio and gamma-ray beams sweep the Earth could correlate with gamma-ray spectral shape. The "X" pulsars in our sample mitigate such a putative radio-selection bias somewhat.

On the other hand, the average Ep of "U" pulsars is higher than for "R" pulsars, with "P" MSPs having higher Ep than the "G" pulsars. MSPs generally have lower than young radio-quiet pulsars, so this reflects the Ep versus trend we report. Our target selection procedures may have tended to favor larger Ep, to be addressed in future work. Because the Ep anticorrelation is fairly modest with substantial scatter, its slope does not change substantially when excluding LAT-induced discoveries.

We thus confirm that our highly diverse pulsar sample is striking in its spectral homogeneity: almost all gamma-ray pulsars are "GeV pulsars," with 0.5 < Ep < 2.5 GeV and following the same trends.

Exceptions exist. Upon its EGRET discovery, Ramanamurthy et al. (1996) highlighted the low gamma-ray efficiency of PSR J0659+1414. Romani et al. (2011) cast it as the archetype of a hypothesized category of pulsars subluminous in gamma-rays. It stands out here as having very low Ep and dp (Figure 20) and the lowest epsilon37,hi (the isolated dot at lower left in Figure 27).

Another clear exception is the Crab: although its measured Ep = 90 ± 80 MeV and spectral shape might be construed as a high- continuation of the trend in Figure 28, its true SED peak occurs in the soft gamma-ray band (Kuiper & Hermsen 2015). A further two pulsars, PSRs B1509−58 and J1846−0258, have spectral shapes similar to GeV pulsars but with values of Ep well outside the continuum of GeV pulsars. Both pulsars have erg s−1 and little or no LAT signal above 100 MeV. Kuiper & Hermsen (2015) highlighted them as archetypes of a pulsar population with SED peaks in the MeV range, and show 18 high- pulsars with SEDs rising in the keV range. Half of the 18 are in this catalog, as are two of five similar pulsars they discuss. PSR B1509−58 has the best-defined MeV peak. They argue that some of the others may not be detected by the LAT because their SEDs also peak below the LAT's range, as for PSR J1846−0258, at the very limit of our instrument's sensitivity (see Section 7.2). However, the SEDs shown by Kuiper & Hermsen (2015) leave open an alternate explanation: PSRs J0205+6449 and J2229+6114 indeed have keV slopes similar to PSR J1846−0258, yet are bright up to a GeV. The as-yet undetected pulsars may simply be too faint, located as they are in very high background regions. Phase-connected rotation ephemerides for a few of the radio-quiet, noisy pulsars are unavailable. Others may be GeV pulsars, but with gamma-ray beams that simply miss the Earth. Harding & Kalapotharakos (2017) further discussed the MeV pulsars.

To see whether the paucity of Ep < 500 MeV pulsars reflects magnetospheric emission processes or stems from observational limits, we start with a tally of LAT candidate sources. There are 238 unidentified nonvariable LAT sources with Ep < 500 MeV, undetected beyond 10 GeV, and significantly curved (well fit with a LogParabola function). Of these, 122 have localization ellipse semimajor axes with 95% C.L. <10', suitable for radio searches (see Table 4). Restricting the search to the 87 sources with ∣b∣ < 10° would favor young, energetic pulsars. Blind searches for radio-quiet pulsars require a large gamma-ray count rate, even more so if timing noise, glitch rates, or background are high. Unfortunately, only four (25) of the 238 targets have TS > 200 (100). The four sources are 4FGL J0340.4+5302, J0426.5+5434, J1749.8-0303, and J1805.7+3401. All were searched several times for radio pulsations but none have been reported to date.

Many high- pulsars have already been found in radio and X-ray searches, effective at finding such sources because high- pulsars tend to have broad radio beams, and because searches can target the supernova remnants near to which such young sources likely remain. But high- pulsars are young, and therefore rare, limited by the Galactic supernova rate. What is rare is far, and thus faint. The subset of high- pulsars for which only the spectral tail reaches the LAT's sensitivity range is thus even rarer. Kuiper et al.'s (2018) detection of PSR J1846−0258 importantly demonstrates the LAT's ability to see such objects. Learning how many more Ep < 500 MeV pulsars exist will likely have to wait for a sensitive, large field-of-view MeV instrument such as AMEGO-X (Caputo et al. 2022) and ASTROMeV (De Angelis et al. 2021). It would allow the missing SED segments for the known and candidate low-Ep pulsars to be filled in, and to search for pulsations from new ones. But again, the ultimate low Ep population is likely to be small compared to that for GeV pulsars; thus, we propose that the trends shown here do capture the spectral character of most gamma-ray pulsars.

8.3. Light-curve Trends

The 237 energy-integrated profile fits in this catalog generally confirm the broad categories of light-curve shapes known since EGRET, consolidated in 1PC and 2PC: a majority of gamma-ray pulsars have two principal peaks separated by Δ ≈ 0.4 ± 0.15 rotations (see Figure 10). For radio-detected pulsars, the first peak generally trails the radio pulse by δ ≈ 0.2 ± 0.2 in phase. Gamma-ray emission from near or beyond the light cylinder can produce the observed anticorrelation between Δ and δ, as discussed in 2PC and explicitly reproduced by the models of Kalapotharakos et al. (2014, 2023). About 40 pulsars have a single broad gamma-ray peak. Very few pulsars deviate from these trends: only two young pulsars in Table 9 are flagged as "complex," appearing to have a third peak half-way between the "standard" two peaks (one is EGRET pulsar B1055-52). This number increases to seven for the MSPs in Table 10. Overall, the MSP and young pulsar profiles are more similar than they differ.

On the other hand, 167 pulsars have enough accumulated photons to allow profile fits in multiple energy sub-bands, fits not done for 2PC. For the first time, we quantify how profiles evolve with energy for a large sample. Figures 13 and 14 highlight the 28 brightest profiles, fit in at least six energy sub-bands. Reminiscent of the spectral results, profile morphology and evolution are remarkably homogeneous given the large diversity of pulsar types and environments. We describe two examples below.

In the second panel from the top, 126 the energy-resolved fit summaries show the phase evolution for P1 and P2, and P3 for the few cases where it is a strong component. Changes in the peak separation Δ are nearly always small, ≲0.1 at most. For <5 pulsars, such as PSR J2229+6114 shown in Figure Set 12, Δ appears to change with energy. However, inspection of the profiles in Figure Set 12 shows that individual profile components vary to accommodate slight pulse shape changes, but Δ in fact varies little.

The summary plots also illustrate the other example. Whereas P3 for the Vela pulsar increases by 0.1 in phase over the LAT energy range, visibly detaching from P1 and moving across the bridge region, this turns out to be the only known pulsar where a gamma-ray peak phase varies substantially with Eγ . Three pulsars, PSRs J0633+1736 (Geminga), J1747-2958, and J1836+5925 ("the next Geminga"), show P3 evolution in the summary plots. Again, inspection of the profiles shows it comes from weak broadening of P2 rather than a peak drift as seen for Vela.

Such weak dependence of pulse phase on energy does not necessarily imply that the region illuminating an observer is uniform. Relativistic aberrations and time delays create caustics such that photons emitted at different distances from the neutron star remain in narrow phase ranges. Thus, whether different electron energies occur at different magnetospheric depths or not, we could see little light-curve evolution with photon energy. Nevertheless, since the SEDs suggest populations of gamma-ray emitting electrons with a narrow range of maximum energies for the lowest pulsars, broadening gradually for systems with more available power, the light curves also encourage vision of emission zones along caustic viewing lines that change little with electron energy. Figure 22 attempts to relate spectral and profile properties: it suggests that when the line of sight cuts across the emission region such that the viewer sees two closely spaced peaks (small Δ), then the line of sight sometimes samples a region with the highest electron energies (large Ep).

Indeed, how the line of sight intersects the emission region is key. Beginning with Watters et al. (2009), and continuing as mastery of magnetospheric processes evolves (e.g., Venter et al. 2009; Bai & Spitkovsky 2010; Contopoulos & Kalapotharakos 2010; Arka & Dubus 2013; Kalapotharakos et al. 2014, 2018, 2023; Cerutti et al. 2016, 2020; Harding 2016; Philippov & Spitkovsky 2018; Pétri 2019), modelers have provided "atlases" allowing for comparison between observed and calculated pulse profiles. Sky maps show beam intensity versus neutron star latitude and longitude, that is, as a function of the inclination ζ between the rotation axis and the line of sight, and the rotational phase ϕ, with other parameters fixed. A cut across a sky map yields the profile sampled by a given line of sight. The dominant "other" parameter is the magnetic inclination α, but, e.g., period or magnetic field strength have more or less importance for different models. Some theories favor emission from plasmoids in current sheets beyond the light cylinder, where magnetic reconnection accelerates electrons to high energies to then radiate synchrotron gamma-rays (Pétri & Lyubarsky 2007; Cerutti et al. 2020). Other models favor electric field acceleration near the current sheet and curvature radiation. In either case, the picture is quite different from the various "gaps" that previously dominated the field (Philippov & Kramer 2022), and the resulting sky maps differ in consequence. The particle-in-cell calculations of the relativistic magnetohydrodynamic processes consume substantial computer resources. Realistic atlases are currently out of reach, with progress in kinetic plasma simulations expected in the coming years.

Figure 29 shows histograms of the gamma-ray duty cycle D, the total phase for which the fitted profile is >25% of the range from its baseline to the maximum. (The supplementary catalog also includes D with 10% and 50% thresholds, and psrcat includes radio pulse width.) Hou et al. (2014) highlighted some pulsars with large duty cycle, and discussed how detection sensitivity depends on D. With two exceptions, the young radio-loud pulsars have D < 0.4. In contrast, three times more young radio-quiet pulsars have 0.3 < D < 0.5 as loud ones.

Figure 29. Refer to the following caption and surrounding text.

Figure 29. Duty cycle D of the gamma-ray pulse profiles (fraction of a rotation that the fitted profile exceeds 25% of the maximum amplitude, measured from the background level). More young radio-quiet pulsars have large duty cycles than do the young radio-loud pulsars.

Standard image High-resolution image

Radio detection occurs when ∣ζα∣ ≲ ρ, that is, along a broad diagonal swath of the plane of neutron star inclinations (α, ζ) (see Figure 3 of Johnston et al. 2020). The opening angle of the radio cone is , for emission height h. Gamma-ray detection favors the α > 30° and ζ > 30° part of the plane. Finally, long duty cycles arise when the line of sight slices a skymap in a ζ range with emission extended in ϕ, but these "flat" beam shapes favor small α (see, e.g., Figure 15 of Kalapotharakos et al. 2018). The tallies in Figure 29 are in good qualitative agreement with this simple geometric beam model. Thus, study of the high duty cycle pulsars may help constrain (α, ζ) and/or discriminate between models. The small light cylinders characteristic of MSPs make for broader beams, the likely reason for the MSP D tail in Figure 29.

If D could constrain the "beaming fraction" fΩ used to calculate luminosity Lγ (Equation (24)), the dispersion in the Lγ versus could be reduced. Large D pulsars, with their distinctive skymaps, could be particularly useful. But they are rare: the α range producing the right beam shape is limited, and the line of sight must fall into a narrow ζ range. PSR J1857+0143, the young radio-loud outlier in Figure 29, may be an example. Serylak et al. (2021) showed its polarization position angle versus phase at 1.4 GHz, but unfortunately scattering excludes an RVM estimate of (α, ζ) at this frequency while its narrow radio pulse width poorly constrains the RVM at higher frequencies. Table 6 points out unidentified LAT sources colocated with known pulsars for which we have not seen gamma-ray pulsations: some may be D ≈ 1 pulsars and could benefit from a (α, ζ) study.

The differences in spectral hardness (Figure 13) between the peaks nearer (P1) and farther (P2) from the magnetic pole, from the observer's point of view, are another powerful discriminant, mentioned above. Harding et al. (2021) and Barnard et al. (2022) addressed this specifically. Early LAT studies measured spectra in as many phase bins as the statistics would allow (Abdo et al. 2010c, 2010g, 2010h; DeCesar 2013). The observed phase evolution of the spectral parameters is smooth, and the range of parameter values is moderate (the cutoff energy can vary by a factor of 2). A simple analysis, measuring on-peak spectral shape in the peak phase intervals we provide, is easier to perform than a full phase-resolved spectral analysis that samples light curves at many points. Furthermore, a larger sample of pulsars is bright enough for on-pulse measurements than for full phase-resolved analysis; yet, a simple analysis is well suited to the current model accuracy. Such future work is encouraged.

8.4. Luminosity and Pulsar Populations

Our large gamma-ray pulsar sample, along with the known radio pulsars, is a starting point for a census of the total Milky Way population. Characterizing the population sheds light on many open questions. These include the supernova rate and neutron star evolution. The large number of "spider" binary MSPs caused a renaissance in the study of binary evolution and pulsar recycling (Strader et al. 2019; van der Merwe et al. 2020; Swihart et al. 2022b). The number and luminosity of high-energy pulsars helps determine the relative contributions of pulsars and putative dark matter annihilation to the observed cosmic positron spectrum (Venter et al. 2015; Cholis & Krommydas 2022) and, via PWNs, to the leptonic cosmic-ray spectrum generally. Applying emission models to a synthetic population and then simply tallying the observed numbers of, for example, radio-loud versus radio-quiet pulsars tests predictions of beam geometry (Johnston et al. 2020).

The strong dependence of Lγ on (Figure 23) and the LAT's sensitivity (Figure 25) are key inputs. Key tests of a synthesis's fidelity are the pulsar latitude distribution (Figure 3) and the cumulatative log N–log S distribution of our pulsar sample, shown in Figure 30. The slopes in the log N–log S plot are related to the spatial distributions of the sources: a flux-limited sample of a population of sources having a narrow luminosity distribution has slope −3/2 if the population is isotropic, and −1 if the distribution is disk-like. Indeed, in Figure 30 the MSPs have a steeper slope, being more isotropic, while the young pulsars have a slope close to −1. The turnover at low fluxes is the detection threshold, lower for objects away from the intense background in the plane. Distance uncertainties do not affect the log N–log S plot. However, an extrapolation from the observed to the unresolved population needs to take selection biases into account. One example is that most of our young radio-loud pulsars were discovered in radio surveys near the plane—results sensitive to details of the high-latitude population may suffer from bias. Another is that target sources for gamma-ray blind searches for radio-quiet pulsars have a higher threshold than shown in Figure 25, explained in Section 3.4.

Figure 30. Refer to the following caption and surrounding text.

Figure 30. The cumulative number of pulsars detected with >100 MeV flux over G100, "log N–log S." The lines show slopes of −1 and −3/2.

Standard image High-resolution image

The scatter around the Lγ versus correlation (or, equivalently, of efficiency versus shown in Figure 24) is a combination of physics and uncertainties. The "physics" bias includes beaming, since fΩ = 1 in our Lγ computation. We use a canonical moment of inertia I = 1045 g cm2 for all pulsars, but the true average is probably twice that, since neutron star radii are closer to 14 km than to 10 km, shifting the X-axis of all plots to the right. Furthermore, MSPs especially have a range of masses, of the order of only 60%, neglected in our single I value. Luminosity may also depend on parameters other than , such as the inclination angles, which would contribute to the scatter. The correlation is stronger at higher spindown powers, and is less pronounced for the MSPs. It remains to be understood whether this is due to high-energy processes "shrinking" as tends toward the "deathline" near 1033 erg s−1, or is because, e.g., of the compactness of the MSP magnetosphere.

Highly polarized radio beams generally occur only for pulsars with erg s−1 (Weltevrede & Johnston 2008; Johnston & Kerr 2018). Smith et al. (2019) suggested it is an unlikely coincidence that both radio polarization and gamma-ray emission require the same minimum spin parameters for independent reasons—that these two phenomena share the same deathline may mean that the same electron populations produce both signals.

8.5. Toward 4PC and 400 Gamma-Ray Pulsars

A small number (∼20) of our gamma-ray pulsars are below the 4FGL catalog detection threshold. Thus, we expect that the bulk of future pulsed discoveries will correspond to the 1000 or more currently unidentified LAT sources. Requiring nonvariability and a pulsar-like spectrum reduces this to a few hundred. It thus appears unlikely that the LAT will ever more than double the size of the current sample.

However, another 100 gamma-ray pulsars is possible, and dozens are likely. Machine-learning methods applied to the 4FGL sources typically find ∼70 pulsar-like unassociated objects with ∣b∣ > 10° (Mirabal & Bonaca 2021). Future high-sensitivity radio surveys with SKA are predicted to find 27,000 young pulsars and 3000 MSPs (Keane et al. 2015). Radio timing solutions for those colocated with LAT source would allow for gamma-ray folding, at least for the years following the radio discovery. A radio position, period, and orbital period for the binaries would improve the sensitivity and decrease the computing cost for gamma-ray blind searches. Without waiting for SKA, TRAPUM at MeerKAT and FAST are already finding new pulsars. Sustained timing of the new discoveries will require thoughtful management of available resources. The insights we have achieved in this work into the shapes of pulsar spectra will help quantify how many of the unidentified LAT sources are likely to be pulsars without pulsed detections.

We intend to maintain the catalog (see the footnote at the end of Section 1). We will add new gamma-ray pulsars as they accumulate. If new analyses are performed, e.g., on-peak spectral analysis, we will add those results.

Pulsar observations enable investigations over a spectacular range of topics, from fundamental physics to variations of the interstellar medium, including several aspects of stellar evolution, and more. Our gamma-ray pulsar sample is so large that even subsamples selected for specific studies have more members than did the entire collection acquired by Compton Gamma-ray Observatory, over 20 yr ago. Discoveries enabled by exploiting the sample will expand for many years to come, at the same time that Fermi LAT keeps finding more, and better characterizes those already known.

Acknowledgments

The Fermi-LAT Collaboration acknowledges generous ongoing support from a number of agencies and institutes that have supported both the development and the operation of the LAT as well as scientific data analysis. These include the National Aeronautics and Space Administration and the Department of Energy in the United States, the Commissariat à l'Energie Atomique and the Centre National de la Recherche Scientifique/Institut National de Physique Nucléaire et de Physique des Particules in France, the Agenzia Spaziale Italiana and the Istituto Nazionale di Fisica Nucleare in Italy, the Ministry of Education, Culture, Sports, Science and Technology (MEXT), High Energy Accelerator Research Organization (KEK) and Japan Aerospace Exploration Agency (JAXA) in Japan, and the K. A. Wallenberg Foundation, the Swedish Research Council and the Swedish National Space Board in Sweden.

Additional support for science analysis during the operations phase is gratefully acknowledged from the Istituto Nazionale di Astrofisica in Italy and the Centre National d'Études Spatiales in France. This work performed in part under DOE contract DE-AC02-76SF00515.

The Parkes Radio Telescope is part of the Australia Telescope, which is funded by the Commonwealth Government for operation as a National Facility managed by CSIRO. The Green Bank Telescope is operated by the National Radio Astronomy Observatory, a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. The Arecibo Observatory is part of the National Astronomy and Ionosphere Center (NAIC), a national research center operated by Cornell University under a cooperative agreement with the National Science Foundation. The Nançay Radio Observatory is operated by the Paris Observatory, associated with the French Centre National de la Recherche Scientifique (CNRS). The Lovell Telescope is owned and operated by the University of Manchester as part of the Jodrell Bank Centre for Astrophysics with support from the Science and Technology Facilities Council of the United Kingdom. The International LOFAR Telescope and the Westerbork Synthesis Radio Telescope are operated by ASTRON, the Netherlands Institute for Radio Astronomy. This work made extensive use of the ATNF pulsar catalog (Manchester et al. 2005). I.G., A.K.H., and D.A.S. thank the ISSI (Bern, Switzerland) for financial support of Team 459 meetings that helped to improve the present work.

Facility: Fermi(LAT) - Fermi Gamma-Ray Space Telescope (formerly GLAST).

Software: Tempo2 (Hobbs et al. 2006), PINT (Luo et al. 2021), astropy (Astropy Collaboration et al. 2018), numpy (Harris et al. 2020).

Note added in proof

Perez et al. (2023) discovered a 2.11 ms pulsar at the position of 4FGL J0212.1+5321.

Appendix A: Radio Timing for This Work

PSR J0653+4706: PSR J0653+4706 was discovered with LOFAR (van Haarlem et al. 2013) at 135 MHz in a targeted pulsar survey toward unidentified Fermi LAT gamma-ray sources (Bassa et al. 2017b; Pleunis et al. 2017). The survey used semicoherent dedispersion to minimize the effects of dispersive smearing at these low frequencies while retaining sensitivity to MSPs (Bassa et al. 2017a). The pulsar has a 4.76 ms spin period and a dispersion measure (DM) = 25.74 ± 0.14 pc cm−3, placing it at a distance of ∼900 pc with both the NE2001 and the YMW16 models. It is in a 5.84 day binary system with an Mc ≳ 0.21 M binary companion (Bassa et al. 2018). The timing solution used in the catalog is obtained from LOFAR timing observations covering a 2.8 yr timespan. Subsequently, a 5.2 yr solution using Jodrell Bank Observatory L-band observations confirmed the model.

The LOFAR data indicate a flux density of 5 ± 2 mJy at 150 MHz. An 1800 s Green Bank Telescope observation at 820 MHz, with 200 MHz bandwidth, yielded a 37.2σ detection. For two polarizations, gain G = 2 K Jy−1, system temperature Tsys = 29 K, and a 12% duty cycle seen in both the 150 and 820 MHz profiles, we obtain a flux density of 260 μJy, calculated using the PRESTO sum_profiles.py routine.

Analysis of 101 positive detections at 1400 MHz with the Lovell telescope at Jodrell Bank Observatory yields flux densities of 89 ± 11 and 10 ± 3 μJy for the main peak and interpulse, respectively. The uncertainties are statistical, and include variations due to interstellar scintillation. The total is S1400 = 99 ± 23 μJy, having added the 20% systematic calibration uncertainty in quadrature. The spectral index obtained from these three measurements agrees well with the average value of 1.7 of the pulsar population.

The JBO ROACH timing data show no significant pulsar proper motion. We found PMDEC = −1.5(2) and PMRA = 0.5(1.8). The Doppler correction is hence 0.98 < (1 − ξ) < 1, with , as in Section 4.3.

PSR J1023+0038: PSR J1023+0038 was the first transitional MSP to be discovered (Archibald et al. 2009), and it underwent a transition from a rotationally powered MSP state to an accretion-powered state in 2013 (Stappers et al. 2014). A faint source coincident with the pulsar was detected by the LAT in the pre-transition data (Nolan et al. 2012), and the flux from this source increased five-fold when accretion began (Stappers et al. 2014).

Using the radio timing ephemeris of Archibald et al. (2013), we folded the pre-transition data (from the beginning of the LAT data to 2013 June), discovering significant pulsations, with H = 62. This ephemeris contains a single orbital period derivative to account for long-term trends in the orbital period, but does not account for shorter-timescale orbital phase variations seen in radio timing observations. We refined this solution using an unbinned timing analysis on the LAT data, also using only one orbital period derivative, resulting in an increased H = 92.

We also folded the post-transition data using a timing model built from X-ray timing by Jaodand et al. (2016) and X-ray and optical timing from Illiano et al. (2023). We fit a polynomial function to the orbital phase deviations measured by Illiano et al. (2023), and folded the LAT data using this orbital model and the spindown model of Jaodand et al. (2016). No pulsations were detected. We also searched for pulsations with slightly different spin periods and spin-down rates, and with slight shifts in orbital phase, using a sliding window approach, but again did not detect any significant pulsations.

This nondetection implies that the flux increase due to accretion is not an enhancement of the pulsed flux (seen in other nonaccreting redbacks; e.g., An et al. 2020; Clark et al. 2021), but is instead an additional unpulsed component. If the five-fold increase in gamma-ray flux is indeed due to unpulsed emission from the accretion process, and the gamma-ray pulsations remain unaffected by the accretion, then the pulsed fraction would be reduced to 20% of the pre-transition level, and this would leave the pulsations undetectable to our analysis. Like PSR J1227−4853 (Johnson et al. 2015), it therefore remains unclear whether or not this system emits gamma-ray pulsations in the accretion-powered state.

PSR J18270849: Applying the methods of Clark et al. (2017) and Section 3.4 to LAT data around the position of a steep-spectrum radio source coincident with unidentified LAT source 3FGL J1827.6-0846, found by Frail et al. (2016), we discovered gamma-ray pulsations with spin period 2.24 ms in 2016. The radio positional uncertainties of ∼2'' have greatly reduced the computing cost and allowed us to increase our sliding coherence window length to 223 s (or ∼97 days) and still perform the search on the ATLAS cluster 127 in about 1 day per source.

The associated radio source has flux density 229.9 ± 14.8 mJy in the GMRT 150 MHz All-Sky Survey, and 777 ± 126 mJy in VLSSr at 74 MHz. NVSS has no detection, giving S1400 < 0.8 mJy per beam. Table 8 lists an insignificant proper-motion estimate from LAT timing as per Section 3.5. C.J. Clark and coworkers are studying other LAT unidentified sources colocated with steep-spectrum radio sources and will detail their methods and results in future work.

PSR J18333840: This 1.87 ms pulsar was discovered in 2015 by F. Camilo, in observations of unidentified LAT sources using the Parkes radio telescope, using the methods reported by Camilo et al. (2015). Later, Frail et al. (2018) independently highlighted it in a search for steep radio spectrum pulsar candidates in LAT unidentified sources using the GMRT 150 MHz all-sky survey (TGSS ADR1). Fewer than 10% of known pulsars have such high spindown power, erg s−1, in the Milky Way or in globular clusters. The orbital period is 0.9 day. In 2020, then-undergraduate Camryn Phillips phase-connected GBT timing observations during PINT code development she was conducting. We phase-folded the gamma-ray data using her ephemeris and obtained the pulsed detection reported here.

PSR J1852-1310: The P = 4.31 ms isolated pulsar PSR J1852−1310 was discovered in a 1900s pointing on MJD 57793 toward an LAT unidentified source with a pulsar-like spectrum, at 800 MHz using the GUPPI backend on the Green Bank Radio Telescope. It was initially named PSR J1852−13. Routine timing observations at 1400 MHz with the Lovell telescope at Jodrell Bank Observatory began on MJD 57818. As of MJD 59927, 124 arrival times have been acquired using ROACH, allowing a phase-connected rotation ephemeris. The spindown power is erg s−1. For DM = 44.95 ± 0.01 pc cm−3 in the direction of this pulsar, the YMW16 model gives a distance of 1.27 ± 0.5 kpc. The integral energy flux of colocated 4FGL J1852.2−1309 thus gives Lγ = (5.8 ± 1.0) × 1032 erg s−1, for a gamma-ray efficiency of 11%.

GBT timing observations at 1400 MHz between MJD 58698 and 59050 yielded 49 times of arrival. Combining these with the JBO data slightly improves the MSP's proper-motion measurement, and we obtain PMDEC = −6.6 ± 1.4 and PMRA = −1.84 ± 0.25. The Doppler correction is hence (1 − ξ) = 0.91, as defined in Section 4.3, such that erg s−1.

Gamma-ray phase folding 14.5 yr of data using the JBO timing model revealed pulsations with test statistic H = 122 using the simple weights parameter μw = 4.3. Model weights increase this to H = 133. H increases steadily for the 8.5 yr before the start of the timing model, indicating that the MSPs spindown rate is stable, well-described using only spin frequency and its first derivative. The radio source is bright (250σ in the discovery observation), and the profile has sharp features; hence, this MSP may rate inclusion in gravitational wave searches.

PSR J1857+0943: This 5.36 ms pulsar has long been known for its stable rotation rate (Segelstein et al. 1986). Gamma-ray phase folding using the rotation ephemeris provided by Reardon et al. (2021) yielded H-test >30, using simple weights. Located in a high background region, (l, b) = (42.2, 3.12) degrees, the gamma-ray profile is ill-defined but appears to be dominated by a single narrow peak. The timing parallax distance from Table 8 and the integral energy flux of the colocated LAT source 4FGL J1857.2+0941 give erg s−1, where the first uncertainty comes from the flux measurement and the second from the distance, about 30% of the spindown power.

PSR J19010125: The same procedure as described for PSR J1827−0849, above, led to the 2016 discovery of 2.79 ms gamma-ray pulsations. F. Camilo then re-analyzed GBT radio data from previous observations of the LAT source position using the LAT rotation ephemeris and discovered radio pulsations and measured the DM. Table 3 lists S1400 = 2.50 mJy, extrapolated from the Frail et al. (2016) GMRT flux density of 362 ± 15 mJy at 150 MHz and radio spectral index −2.23. Table 8 lists a 2.4σ proper-motion measurement using the methods of Section 3.5. Clark and co-workers are studying other LAT unidentified sources colocated with steep-spectrum radio sources and will detail their methods and results in future work.

PSR J22561024: Crowter et al. (2020) found μ = 9 ± 3mas yr−1 using GBT timing data. NRT has also monitored this pulsar, but with a substantial gap between the end of the GBT data and the beginning of the NRT observations. This gap and limited frequency coverage cause a strong degeneracy between the spindown parameters, time-varying DM, and orbital frequency evolution. We therefore perform a hybrid analysis, using the methods of Section 3.5, in which we use LAT data to determine the "long-term" parameters (F0, F1, FB0, FB1, and FB2). We then fix these parameters and fit for the astrometry using the radio data, finding a proper motion of PMRA = 2.80 ± 0.05, PMDEC = −5.30 ± 0.11 from which μ = 6.0 ± 0.1 mas yr−1. A systematic error comparable to or larger than the statistical uncertainty is likely, due to the inhomogeneous data used in the analysis.

Appendix B: Typical Gamma-Ray Spectral Energy Distrbutions

Figure B1 shows sample SEDs provided as part of 4FGL-DR3. 128 The figure set for this pulsar catalog includes a compilation of the SEDs of the gamma-ray pulsars. Additional material is available at https://fermi.gsfc.nasa.gov/ssc/data/access/lat/3rd_PSR_catalog/. This also includes the spectral parameters with their covariance matrices and the data points necessary to reproduce the SEDs, and a python script that performs the task.

Figure B1. Refer to the following caption and surrounding text.

Figure B1.

Typical gamma-ray spectral energy distributions (SEDs). The plots were made using the python script sed_plotter.py and data files. The SEDs for all of the gamma-ray pulsars detected as point sources (TS > 25) are also available. Similar to those provided with 4FGL-DR3, they use the refits done for some pulsars, described in Section 6. When bfr exists, both b23 and bfr fits are shown, with the bfr uncertainty envelope. The complete figure set (294 images) is available. (The complete figure set (294 images) is available.)

Standard image High-resolution image

Appendix C: Additional Notes on Spectral Analysis

We study here the distribution of the model parameters and derived physical quantities as they vary among the 116 pulsars, which have both b-free and b = 2/3 fits available. These results are presented in Figure C1. In brief, they are as follows:

  • 1.  
    Γ, the spectral index at the reference energy E0, is stable, varying by <10% over the effective range in b.
  • 2.  
    The curvature at the reference energy, d, varies more than Γ, by up to ∼50%. If considering b < 1, the total range is only about ±20%.
  • 3.  
    Γ100 varies almost linearly with b before saturating at b = 1.2. The shift is up to 1 unit harder for pulsars with b near 0, and 0.5 units softer for (mostly) MSPs with b ≥ 1.2. Thus, clearly, freeing b profoundly alters the inferred low-energy spectral index.
  • 4.  
    The curvature at the peak, dp, is stable if the intrinsic value of b > 0.5. For pulsars with spectra more consistent with small b, the b-free pulsars become sharper, with a difference in log-width ranging up to 2×. It is not (yet) clear if this difference primarily "fixes" the true width of the peak or delivers a better low-energy index.
  • 5.  
    The peak energy Ep varies by <10%–15% in the "core" b range and <25% overall. The largest deviations are for the MSPs with fit values b > 1, which move to higher Ep.

In summary, the physical properties measured at the peak, dp and Ep, are robustly measured in both b-free and b = 2/3 fits, and so we can make inferences on their distribution using the entire pulsar sample. On the other hand, the low-energy spectrum is unreliable, and values inferred from the b = 2/3 model could be very far from the true values.

Figure C1. Refer to the following caption and surrounding text.

Figure C1. Comparisons of base and derived parameters between the and b-free populations. See the main text for further discussion.

Standard image High-resolution image

Appendix D: Description of the Catalog Files

Supplementary material is available at https://fermi.gsfc.nasa.gov/ssc/data/access/lat/3rd_PSR_catalog/. Here we describe the files posted as of publication. Additional material may be added in the future. Static versions of some data used in this work are also available. The "xxx" in the file names is a multidigit time stamp of the file creation date. Contact the authors to request missing information.

3PC_Catalog+SEDs_xxx.fits, 3PC_Catalog_xxx.xls: The pulsar catalog "flat file," in FITS and spreadsheet formats. Each has 294 rows, one for each gamma-ray pulsar. The 112 columns contain basic pulsar properties (name, position, spin parameters, distance, radio flux density, etc.) as well as the spectral and profile parameters described in this article. The FITS version includes a subtable ("HDU") with the detailed spectral fit parameters and covariance matrix needed to reproduce the SEDs, provided in 3PC_SEDs.tgz (individual plots) and 3PC_AllSEDs.pdf (all plots concatenated into a single file). 3PC_SEDPlotter.tgz contains a python script that reads 3PC_Catalog+SEDs_xxx.fits to generate the SED plots. Most Figures and Tables in this article were created from the FITS and spreadsheet files: the python script Example_3PC_CatalogReader.py is a simple example. README_3PC_Catalog+SEDs_xxx.txt provides more detail. These materials are available with the article in the datacat.tar.gz package.

3PC_ProfileData_xxx.tgz: Contains two files for each pulsar, one in FITS format and the other in simple text (ascii), containing the weighted gamma-ray phase histograms and, when they exist, the fit to the histogram and/or the radio profile. The python script Example_3PC_PlotProfileData.py illustrates how to read the files to reproduce the plots provided in 3PC_ProfilePlotsPDF_xxx.tgz and 3PC_ProfilePlotsPNG_xxx.tgz. The plots of gamma-ray peak evolution with energy, as in Figure Set 12, are provided in 3PC_LC_Eres_pdf_xxx.tgz, 3PC_LC_Eres_png_xxx.tgz. README_ProfileData.txt provides details. These materials are available with the article in the datapro.tar.gz package.

3PC_TimingModels_xxx.tgz: A compressed archive of the 294 rotation ephemerides (".par files") used to calculate the gamma-ray rotational phases using software such as Tempo2 or PINT. README_3PC_TimingModels.txt provides details. These materials are available with in the datapar.tar.gz package.

Phased, weighted FT1 FITS files: available at https://heasarc.gsfc.nasa.gov/FTP/fermi/data/lat/catalogs/3PC/photon/. Two subdirectories called 3deg_50MeV and 15deg_20MeV each contain one file for each gamma-ray pulsar, in the standard Fermi LAT data format, with additional columns PHASE and MODEL_WEIGHT. The 277 files in 3deg_50MeV include photons with energies >50 MeV, within 3° of the pulsar position, and were used to create the gamma-ray phase histograms. They occupy roughly 30 GB of disk space, with a large variation in single file size depending on the background intensity (generally decreasing away from the Galactic plane) and pulsar intensity. The files in 15deg_20MeV are on average 15× larger, including photons with energies >20 MeV within 15° of the pulsar position. They are suitable for spectral analyses (e.g., using the gtlike Fermi tool) or to explore the lowest LAT photon energies. README_3PC_PhotonData_.txt describes the FT1 files in more detail.

3PC_sensitivity_xxx.fits: The all-sky Fermi LAT pulsar sensitivity, in FITS format, is useful for, e.g., population syntheses. Example_3PC_SensitivityMap_xxx.py is a python script that reads the FITS file and generates a plot as in Figure 25. README_3PC_SensitivityMap_xxx.txt provides details.

Finally, we provide HTML files for each pulsar, to allow for a visual summary of the information for each gamma-ray pulsar, named 3PC/J0007+7303_LAT.html and so on.

Footnotes

10.3847/1538-4357/acee67
undefined