Skip to main content

Full text of "Three Millisecond Pulsars in FERMI LAT Unassociated Bright Sources"

See other formats

Accepted for publication in ApJ Letters 

Preprint typeset using I^'T^]X style emulatoapj v. 11/10/09 


S. M. RansomI-2, p. S. RAY''•^ F. Camilo", M. S. E. Roberts", 6. Qelik'-^-'', M. T. Wolff'\ C. C. Cheung"', 
M. Kerr", T. Pennucci", M. E. DeCesar' "^, I. Cognard^-', A. G. LYNE^^ B. W. Stapfers^''^ P. C. C. Freire^ 
J. E. Grove^, A. A. Abdo^", G. Desvignes^''^**, D. Donato*"', E. C. Ferrara^, N. Gehrels', L. Guillemot"', 
C. GwoN^ A. K. Harding^ S. JoHNSTON"^ M. Keith"\ M. Kramer"'' "', P. F. Michelson", D. Parent^", ^ 
P. M. Saz Parkinson^\ R. W. Romani", D. A. Smith^^, G. Theureau^^, D. J. Thompson', P. Weltevrede^^, 

K. S. Wood-', M. Ziegler2\ 

Accepted for publication in ApJ Letters 


We searched for radio pulsars in 25 of the non-variable, unassociated sources in the Fermi LAT 
Bright Source List with the Green Bank Telescope at 820 MHz. We report the discovery of three radio 
and 7-ray millisecond pulsars (MSPs) from a high Galactic latitude subset of these sources. All of the 
pulsars are in binary systems, which would have made them virtually impossible to detect in blind 
7-ray pulsation searches. They seem to be relatively normal, nearby (<2kpc) millisecond pulsars. 
These observations, in combination with the Fermi detection of 7-rays from other known radio MSPs, 
imply that most, if not all, radio MSPs are efficient 7-ray producers. The 7-ray spectra of the pulsars 
are power-law in nature with exponential cutoffs at a few GeV, as has been found with most other 
pulsars. The MSPs have all been detected as X-ray point sources. Their soft X-ray luminosities of 
^-[^q30-31 gj-gg-i are typical of the rare radio MSPs seen in X-rays. 

Subject headings: pulsars: general — pulsars: individual (J0614— 3329, J1231 — 1411, J2214+3000) 

^ National Radio Astronomy Observatory (NRAO), Char- 
lottesville, VA 22903, USA 
^ email: 

^ Space Science Division, Naval Research Laboratory, Wash- 
ington, DC 20375, USA 

* email: 

^ Columbia Astrophysics Laboratory, Columbia University, 
New York, NY 10027, USA 

6 Eureka Scientific, Oakland, CA 94602, USA 

7 NASA Goddard Space Flight Center, Greenbelt, MD 20771, 

® Center for Research and Exploration in Space Science and 
Technology (CRESST) and NASA Goddard Space Flight Cen- 
ter Greenbelt, MD 20771, USA 

^ Department of Physics and Center for Space Sciences and 
Technology, University of Maryland Baltimore County, Balti- 
more, MD 21250, USA 

National Research Council Research Associate, National 
Academy of Sciences, Washington, DC 20001, resident at Naval 
Research Laboratory, Washington, DC 20375, USA 

W. W. Hansen Experimental Physics Laboratory, Kavli In- 
stitute for Particle Astrophysics and Cosmology, Department of 
Physics and SLAG National Accelerator Laboratory, Stanford 
University Stanford, CA 94305, USA 

12 University of Virginia, Charlottesville, VA 22904, USA 

Department of Physics and Department of Astronomy, Uni- 
versity of Maryland, College Park, MD 20742, USA 

^'^ Laboratoire de Physique et Chimie de I'Environnement, 
LPCE UMR 6115 CNRS, F-45071 Orleans Cedex 02, and 
Station de radioastronomie de Nangay, Observatoire de Paris, 
CNRS/INSU, F-18330 Nangay France 

Jodrell Bank Centre for Astrophysics, School of Physics and 
Astronomy, The University of Manchester, M13 9PL, UK 

Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 
69, 53121 Bonn, Germany 

Department of Astronomy, University of California, Berke- 
ley CA 94720-3411, USA 

^ Radio Astronomy Laboratory, University of California, 
Berkeley, CA 94720, USA 

CSIRO Astronomy and Space Science, Australia Telescope 
National Facility, Epping NSW 1710, Australia 

College of Science, George Mason University, Fairfax, VA 
22030, resident at Naval Research Laboratory, Washington, DC 
20375, USA 

Santa Cruz Institute for Particle Physics, Department of 

Physics and Department of Astronomy and Astrophysics, Uni- 
versity of California at Santa Cruz, Santa Cruz, CA 95064, USA 
22 Universite Bordeaux 1, CNRS/IN2p3, Centre d'Etudes 
Nucleaires de Bordeaux Gradignan, 33175 Gradignan, France 


Ransom et al. 


Before the launch of the Fermi Gamma-ray Space Tele- 
scope, the only pulsars with definitive detections in 7- 
rays (using EGRET on CGRO) were young and very en- 
ergetic {E > 10^^ ergs""'^) or nearby older systems {E > 
lO^'^ergs"^) (jThoinpsonI 120041 ) . A possible detection 
of pulsed 7-rays fr om the energetic MSP J0218-H4232 
(jKuiper et al.|[2000(), sparked interest in modeling MSP 
7-ray emission fe.g. lZhang fc Chengjl2003l:lHarding et al.l 
mm\ . and encouraged one group (|Storv et al.ll2007D to 
predict that many new MSPs might be detected in 7- 
rays or discovered in radio follow-up of unidentified Fermi 

The launch of Fermi and the ext raordinary sensitivit y 
of the Large Area Telescope (LAT, lAtwood et aIll2009D . 
confirmed those predictions of 7-ray-bright MSPs with 
detections of eight relatively normal radio MSPs using 
only the first few months of Fermi events (jAbdo et al.l 
[2009a). Those MSPs were detected via the folding of 
7-rays modulo the known spi n and orbital ephe merides 
from radio timing campaigns (ISmith et al.l[2008l) . 

In order to best utilize radio telescope time to search 
either for radio counterparts to new 7-ray-selected pul- 
sars or to search blindly for radio pulsations from 7-ray 
sources that might contain pulsars, we formed the Pul- 
sar Search Consortium (PSC), a group of approximately 
20 LAT-team members and/or pulsar experts associated 
with large radio telescopes around the world. This pa- 
per describes one of the PSC's first programs, which used 
the Green Bank Telescope (GBT) to search 25 unasso- 
ciated sources from the Fermi LAT Bright Source List 
(lAbdo et al. l l2009bD . 


We selected 25 sources from the Fermi LAT Bright 
Source List that were a) unassociated with known pul- 
sars or active galactic nuclei (AGN) b) unassociated 
with X-ray counterparts that had been previously deeply 
searched for radio pulsations (e.g. IC443, Camilo et al., 
in prep.) c) statistically non-variable and d) at decli- 
nations >— 35°. We observed each of the sourceS for 
approximately 45— 50 minutes using the prime focus re- 
ceiver at the GBT centered at 820 MHz with 200 MHz 
of bandwidth. The GBT pointings, all taken between 
July and October 2009, were actually centered on the 
positions from an internal LAT source list using nine 
months of sky-survey data prepared in a similar fashion 
to the Fermi LAT First Source Catalog (i.e. "IFGL" , 
lAbdo^ et all [2(31(1 . The individual GBT pointings at 
820 MHz had FWHM=0.25° and covered either all or a 
substantial fraction of the 95% error regions for the vast 
majority of the sources. 

We sampled the summed power from two polarizations 
in 2048 frequency channels with 8-bits every 61.44 /is us- 
ing the GUPPI pulsar backentQ- Each pointing gener- 
is The LAT Bright Sources observed with the GBT were OFGLs 
J0614.3-3330, J1231. 5-1410, J1311.9-3419, J1653. 4-0200, 
J1741. 4-3046, J1746.0-2900, J1801.6-2327, J1805. 3-2138, 
J1814.3-1739, J1821.4-1444, J1834.4-0841, J1836. 1-0727. 
J1839. 0-0549, J1844.1-0335, J1848.6-0138, J1855. 9-1-0126. 
J1900. 0-1-0356, J1911.0-I-0905, J1923.0+1411, J2001. 0+4352, 
J2027.5-I-3334, J2110. 8-1-4608, J2214.8+3002, J2302. 9-1-4443, and 

J2339. 8-0530. 

'^'^ https : //safe . nrao . edu/wiki/bin/view/CICADA/NGNPPl 

ated approximately 100 GB of data, which were recorded 
to hard drives for processing off-site. The 820 MHz cen- 
ter frequency was chosen as a compromise between the 
competing effects of sky temperature (from the Galactic 
synchrotron background) and beamsize, as well as steep 
pulsar spectra and the effects of interstellar dispersion 
and scattering. 

For a pulsar with a pulse width of '^10% of the pe- 
riod, the search sensitivity was approximately 0.06 x 
(29 K -I- Tsky)/ (32 K) mJy, where Tgky is the contribution 
at 820 MHz of the Galactic synchrotron background. The 
majority of the sources (17 of them) were within three 
degrees of the Galactic plane where Tgky ^ 12 — 40 K, 
although two were very near the Galactic center with 
Jsky ^ 100 — 150 K. Eight of the sources were well off the 
Galactic plane (|6| > 5°) and had T^ky - 3 - 10 K. For 
those sources, our search sensitivity was 0.06— 0.08 mJy 
for normal pulsars at all reasonable dispersion measures 
(DMs) and MSPs up to DM-lOOpccm-^. In general, 
the observations were factors of 2—12 deeper than the 
best pulsar surveys that have previously covered thes e 
regions (e.g. iManchester eran[200lL iCordes et al1 [200g) . 

We processed the data, after de-dispersing into ^9000 
DMs over the range 0— 1055pccm~^, using both acceler- 
ation searches (to improve sensitivity to pulsars in binary 
systems) and single pulse searches (to provide sensitiv- 
ity to pulsars with sporadic or giant-pulse-like emission) 
using standard tools found in PRESTCE!I(R a nsom et al.l 
[2OOI . No new pulsar-like signals were found in any of the 
low Galactic latitude {\b\ < 5°) sources. However, four 
new MSPs were detected amid the eight high Galactic 
latitude sources, in OFGLs J0614.3-3330, J1231. 5-1410, 
J2214.8+3002, and J2302.9-I-4443. The MSP in OFGL 
J2302. 9+4443 was detected first in an independent PSC 
survey by the Nangay telescope and will be reported else- 
where (Cognard et al. in prep.). The rest of this paper 
details the properties of the other three MSPs. 

2.1. The New MSPs 

The first two pulsars detected, J2214+3000 and 
J1231 — 1411, were undergoing substantial accelerations 
due to orbital motion during the discovery observations. 
PSR J0614— 3329 was initially uncovered in an unacceler- 
ated search, although orbital motion was detected in the 
discovery observations via a more precise timing analysis. 
The fact that all three MSPs were in unknown binaries 
demanded a radio timing program to determine precise 
orbital parameters and constrain their astrometric po- 
sitions before detailed 7-ray timing and analysis could 

Each MSP was observed with several different observ- 
ing setups at the GBT, the Lovell Telescope at Jodrell 
Bank, and the Nangay radio telescope. The Arecibo tele- 
scope also observed PSR J2214+3000 several times. At 
the GBT, GUPPI was used with bandwidths of 100, 200 
and 800 MHz centered at 350, 820, and 1500 MHz re- 
spectively. At Jodrell Bank and Nangay, observations 
were made with bandwidths of 200— 300 MHz centered 
near 1400 MHz. Standard radio timing procedures were 
used (Lorimer & Kramer 2005) and the or bital parame- 
ters w ere fit to high precision with TEMPD2 (|Hobbs et al.l 
[2006h . 

Ihttp : //www . cv . nrao . edu/- sransom/presto/ 1 

MSPs in Fermi Sources 


After several months of radio timing, and using trial 
X-ray positions based on point sources from Swift, 
XMM-Newton, and/or Chandra, 7-ray pulsations from 
each MSP were detected with the LAT using orbitally- 
demodulated events (see also H2.2p . We determined aver- 
age 7-ray pulse Times of Arrival (TOAs) using the max- 
imum like l ihood 7-ray timing techniques described by 
I Ray et all (|201Clf ) after integrating source photons mod- 
ulo the predicted pulse period for between 22 and 36 
days per TOA. The resulting joint timing solutions, us- 
ing both radio and 7-ray TOAs, as well as the derived 
physical parameters of the MSPs, are presented in Ta- 
ble 1. The radio and 7-ray pulse profiles are shown in 
Figure [2l 

The three new pulsars, besides being three of 
the brightest 7-ray MSPs in the sky, appear to 
be relatively normal, nearby (<2kpc), radio MSPs, 
with ~3 ms spin periods, surface magnetic field 
strengths of (2— 3)xlO^G, and spin-down luminosi- 
ties of ~2x 10^4 erg s-i. J0614-3329 and J1231-1411 
have orbital periods of 53.6 and 1.9 days respectively, 
with companions of mass ^0.2— 0.3 Mq, consistent with 
the orbital peri o d - white dwarf mass relation of 
iRappaport etall (fl995l ). J2214-H3000 is a so-called 
"black-widow" system with a very low-mass compan- 
ion ('^0.02 M0) and likely timing irregula rities, similar 
to pulsars B 1 957+20 fFruchter et a l.|[T98l) J205 1-0827 
JStappers et al]|199^, and J0610-2100 (|Burgav et all 
[2006I) . and only the fourth such system known in the 
Galactic disk. While we currently have no evidence for 
radio eclipses from the pulsar (at least at frequencies 
>1.4GHz), its formatio n was li kely similar to that of 
those other systems (e.g. iKing^t al. 2005). 

2.2. J -ray Analysis 

The Fermi L AT is sensitive to 7 -rays with energies 
0.02-300GeV (jAtwood et all [20091 ). The Fermi LAT 
sky-survey data set used here for spectral analysis spans 
from 2008 August 4 to 2010 February 4. We selected 
"Pass 6 Diffuse" class events - i.e. events passing the 
most stringent background rejection cuts - with ener- 
gies above 0.1 GeV and rejected events with zenith an- 
gles > 105° to limit contamination from 7-rays from the 
Earth's limb. We used "Pass6 v3" instrument response 
functions (IRFs). 

The 7-ray light curves shown in Figure [T] are con- 
structed from events as described above, although using 
seven additional months of data through 2010 Septem- 
ber 14, and with energy and radius cuts to optimizes 
the signal to noise for each pulsar. The energy and ra- 
dius cuts used for PSRs J0614-3329, J1231-1411, and 
J2214-h3000 are (1.0°, 0.35 GeV), (1.2°, 0.35 GeV), and 
(1.0°, 0.7 GeV), respectively. 

We derived the 7-ray spectrum of each pulsar us- 
ing a maximum-likelihood method implemented in the 
LAT Science Tool gtlikj^. We analyzed a region of 
10° radius centered on the radio position of each pul- 
sar, and modeled each region by i ncluding all sources 
from the IFGL (|Abdo et al.l l2010f) within 17° of the 
pulsar along with Galactic and isotropic diffuse emis- 
sion (models gll_iem_v02 and isotropic_iem_v02, re- 

FlG. 1. — Radio and 7-ray pulse profiles for the three new MSPs. 
The red lines correspond to the 820 MHz discovery pulse profiles 
from the GBT with the 820 MHz flux density scale on the right. 
The blue lines are the Fermi pulse profiles with the photon counts 
and low energy cut used for the 7-ray selections listed on the left. 
There are 2043, 2341, and 621 photons in the 7-ray profiles for 
PSRs J0614-3329, J1231-1411, and J2214+3000, respectively 

spectiveljO)- The power-law spectral parameters for all 
sources within 10° of the pulsar and a normalizing scale 
factor for the diffuse emission spectrum were allowed to 
be free in the fit. 

We modeled the spectrum of each MSP using a power- 
law with an exponential cutoff where the three parame- 
ters, the differential flux K, the photon index F, and the 
cutoff energy i?cutoff, were allowed to vary in the fit. The 
phase-averaged 7-ray spectra obtained for each pulsar are 
shown in Figure [2] and the spectral parameters are given 
in Table 1. The uncertainty in the LAT effective area 
is estimated to be <5% near 1 GeV, 10% below 0.1 GeV 
and 20% over 10 GeV. The resulting systematic errors 
on the three spectral parameters, propagated from the 
uncertainties on the LAT effective area, were calculated 
using a set of "modified IRFs" bracketing the nominal 
(Pass6 v3) one. 

We verified the significance of the exponential c utoff in 
each spectrum with a likelihood ratio test (Mattox et al.l 
119960 . A simple power-law model is rejected significantly 
for all three pulsars relative to an exponentially cutoff 
power-law, as indicated by the test statistic ATScutoff 
listed in Table 1 for the addition of one free parameter. 

The flux points in Figure [5] were obtained by repeating 
the likelihood analysis in each energy band, assuming a 
power-law spectrum with a photon index fixed at 2 and 
a free flux normalization parameter for all sources. 

2.3. X-ray Analysis 
To search for X-ray counterparts. 


field of each pulsar with the Swift-XKT 
120051 ) with exposures ranging from 2.6— 15.9 ks. For 

observed the 
Burrows et al.l 

^ |http: //ferml .gsfc. analysis/document at ion/ I ^^ |http: //ferml .gsfc .h1 


Ransom et al. 

10"' 1 10 

Energy (GeV) 

Fig. 2. — The 7-ray spectra for the three MSPs. The flux points 
on the curve were obtained from independent fits in each energy 
bin, as explained in the text. The curves represent the spectrum 
fit with a simple exponentially cut off power-law in the full energy 
range of 0.1 - 100 GeV. 

J0614-3329 and J1231-1411 we detected X-ray point 
sources at the locations of the pulsars. In the Swift obser- 
vation of PSR J2214-I-3000, there is no significant source 
at the location of the pulsar, however we also obtained a 
Chandra ACIS-I observation of this region and detected 
the pulsar with it. A detailed analysis of this Chandra 
observation will be presented elsewhere. 

The X-ray counterparts to the MSPs are soft sources 
and were fit to black body spectra using XSPEC, fix- 
ing the absorption values to zero. The resultant tem- 
peratures for the three X-ray sources were in the range 
~0.21-0.25 keV, which in the case of J1231-1411, is 
consistent with that derived from deeper Suzaku data 
(Maeda et al. 2010, in prep.). Allowing for additional 
Galactic absorpt ion fixed to the valu es of (3.4 — 5.7) x 
10 cm' -2 from (IKalberla et al.l[2005l) . we found negligi- 
ble differences in the fitted parameters. The Swift and 
Chandra positional localizations, as well as approximate 
fluxes and luminosities from the black body fits, are in 
Table 1. 

For PSR J1231 — 1411 we were able to do more detailed 
spectral analysis using XMM-Newton. On 2009 July 15 
XMM-Newton observed the field of this as- yet unidenti- 
fied BSL source with all three EPIC instruments. Data 
from each instrument were analyzed utilizing the Science 
Analysis System software version 10.0.0 and the calibra- 
tion update of May 8, 2010. We filtered the data for bad 
events and excluded times of high particle background, 
yielding 29.4 ks and 29.5 ks of good time for the EPIC- 
MOSl and -M0S2 instruments and 24.1 ks for EPIC-PN. 
All three instruments utilized Full Frame mode so none 
had sufficient time resolution to allow searches for X-ray 

The XMM-Newton X-ray images of the J1231-1411 
field reveal a relatively isolated, moderately bright, point 
source that we name XMMUJ123112-141146 at the best 
radio position of the pulsar to the accuracy of the EPIC 
instruments. We generated spectra of this source by ex- 
tracting events from around the pulsar in 35" regions for 
the MOS images and a 25" region (due to a chip gap) 
for the PN image. This yielded spectra consisting of 
853, 812, and 2164 events from the MOSl, M0S2, and 
PN cameras, respectively. We group the counts with at 

least 30 events per spectral bin for the MOS spectra and 
45 events per bin for the PN spectra using the FTOOL 
grppha. Background spectra were extracted from nearby 
regions 100" and 55" in radius from the same CCD chips 
for the MOS and PN instruments, respectively. 

Using XSPEC (vl2.6.0) we simultaneously fit the 
MOSl, M0S2, and PN spectra in the energy range 
0.4—3.0 keV. An absorbed power-law model is formally 
acceptable with a reduced x^=0.977 (61 dof) but with 
an extremely steep photon index of F = 4.23lQ3g (90% 
confidence) and column density of uh ~ {l.8t.Q'^) x 
lO^-'^cm"^ (90% confidence). Such a column density is 
significantly higher than that expected for this direc- 
tion based on galact ic surveys {nn = 3.45 x lO^'^cm^^, 
IKalberla et al.ll2005D . A fit to these data utilizing a 
model of an absorbed neutron star non-magnetic hy- 
drogen atmosphere (phabsxnsatmos AHeuike et al.l[2006[ l 
with the neutron star mass and radius held fixed at 
1.4Mq and 10km, and the source distance fixed at 
the dispersion measure value of 0.4 kpc, yields a re- 
duced x^=1.41 (61 dof). In this model there are sig- 
nificant residuals above 1.5 keV so we add a power- law 
component with photon index fixed at 1.8 and obtain 
an improved fit with (60 dof). This latter 

fit yields an atmospheric temperature (seen at infinity) 
Toff = 6lt^2eV (90% confidence), a best fit Uh con- 
sistent with zero and a 90% confidence upper limit of 
5 X 10^°cm~^, and flux in the 0.5— 3keV energy band 
(1.15 ± 0.05) X 10"^'^ ergscm"^ s^^. In this model the 
power-law spectral component accounts for roughly 25% 
of the total flux in the 0.5— 3keV energy band. 


We have identified three new nearby radio MSPs as 
the counterparts of bright and previously unassociated 
Fermi LAT sources at high Galactic latitude. Our non- 
detection of young pulsars or MSPs in the more numer- 
ous sources searched at low Galactic latitude is likely due 
to our only moderate sensitivity improvements (typically 
2— 3x) over the best sur veys of those regions to date (e.g. 
[Manchester et all 1200 It ) due to higher sky temperatures 
resulting from our lower observing frequency. Addition- 
ally, the complicated and confused nature of the Galactic 
plane in 7-rays makes the positive identification of point 
sources difficult. Several of the bright sources may be 
blends of other sources or the result of insufficient mod- 
elling of the Galactic background. Nonetheless, deeper 
surveys at frequencies of 1.5— 2 GHz of these sources may 
prove more fruitful in the future. 

The new pulsars are very typical radio MSPs in 
terms of spin period, binary parameters, magnetic field 
strength, spin-down luminosity, and characteristic age, 
and their unusual brightness in 7-rays is likely due more 
to their proximity than to especially energetic emission 
processes in their magnetospheres. The very high im- 
plied 7-ray efficiency for PSR J0614— 3329 suggests it is 
likely closer, by up to a fac tor of 2 or m ore, than pre- 
dicted by the NE2001 model (jCordes k Laz io 2002). The 
line-of-sight to PSR J0614— 3329 is nearly tangent to the 
Gum Nebula where NE2001 shows an exceptionally steep 
DM gradient. Additionally, the pulsar's 7-ray emission 
is likely not isotropic, but only covers tens of percent 
of the sky. These large efficiencies in general, though. 

MSPs in Fermi Sources 


are consistent with the tens of percent values found by 
lAbdo et all (|2009al) for radio MSPs detected in 7-rays 
and imply that MSPs are very efficient producers of 7- 

We do not have proper motion measurements for pul- 
sars J0614— 3329 or J1231 — 1411 and so their measured 
spin-down rates are contamina ted at some leve l (likely 
<10%) by the Shklovskii effect (IShklovskiilll970l) . There 
is a statistically significant proper motion measurement 
of ^lOOmasyr"^ for J1231 — 1411, though, which implies 
a Shklovskii effect at 400 pc larger than the measured 
spin-down rate for the pulsar. If the proper motion is 
confirmed at this level, the requirement to have the pul- 
sar intrinsically spinning down gives an upper limit for 
the pulsar's distance of ~240pc. Timing observations 
over the next several years will determine the proper mo- 
tions and possibly the timing parallaxes for each of the 

In all three cases we identified X-ray counterparts to 
the pulsars which substantially aided in the rapid estab- 
lishment of timing solutions. The three MSPs appear 
to h ave fairly typic al X-ray properties for radio MSPs 
(e.g. iBogdanov et a l. 2006) with primarily soft thermal- 
like spectra and X-ray luminosities in the 10^°"^^ ergs~^ 
range, approximately 10~^ to 10"'^ of their 7-ray lumi- 

The radio flux densities of ^ImJy near 1 GHz are 
large enough to make the MSPs potentially useful for 
a wide variety of timing projects, such as the detection 
of gravitational waves via long-term pulsar timing (e.g. 
NAN0Gra\|3), yet they are small enough to explain why 
earlier large-area surve ys for pu lsars missed them (e.g. 
[Manchester et al. 1996'^ iLorimer ct al. 2005). In addi- 
tion, the fact that many of the nearby radio MSPs are 

being detected in 7-rays and vice-versa argues that the 
sizes of the radio and 7-ray beams are comparable for 
MSPs (likely within a factor of ~2), and that deep radio 
and 7-ray surveys may allow us to eventually detect a 
large percentage of the local population of these sources. 
In the short-term, the fact that Fermi can point us to 
nearby radio MSPs is already causing a large increase in 
the number of known systems, with much less effort than 
is required by sensitive large-area radio surveys. If most 
radio MSPs produce 7-rays as these early results seem 
to indicate, MSPs may contribute to the diffuse iso tropic 
7-ray background (|Faucher-Giguere fc Loebir2010l ). 

We acknowledge helpful discussions with Natalie Webb 
and Lynne Valencic. The National Radio Astronomy 
Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Asso- 
ciated Universities, Inc. This work was partially sup- 
ported by NASA Grant No. NNG09EE57I. The Fermi 
LAT Collaboration acknowledges support from a num- 
ber of agencies and institutes for both development and 
the operation of the LAT as well as scientific data analy- 
sis. These include NASA and DOE in the US, CEA/Irfu 
and IN2P3/CNRS in France, ASI and INFN in Italy, 
MEXT, KEK, and JAXA in Japan, and the K. A. Wal- 
lenberg Foundation, the Swedish Research Council and 
the National Space Board in Sweden. Additional sup- 
port from INAF in Italy and CNES in France for science 
analysis during the operations phase is also gratefully 

Facilities: GET (GUPPI), Fermi (LAT), XMM 
(EPIC), CXO (ACIS), Swift (XRT) 


Abdo, A. A., et al. 2009a, Science, 325, 848 

— . 2009b, ApJS, 183, 46 

— . 2010, ApJS, 188, 405 

Atwood, W. B., et al. 2009, ApJ, 697, 1071 

Bogdanov, S., Grindlay, J. E., Heinke, C. O., Camilo, P., Freire, 

P. C. C, & Becker, W. 2006, ApJ, 646, 1104 
Burgay, M., et al. 2006, MNRAS, 368, 283 
Burrows, D. N., et al. 2005, Space Sci. Rev., 120, 165 
Cordes, J. M., & Lazio, T. J. W. 2002, astro-ph/0207156 
Cordes, J. M., et al. 2006, ApJ, 637, 446 

Paucher-Giguere, C, & Loeb, A. 2010, Journal of Cosmology and 

Astro-Particle Physics, 1, 5 
Pruchter, A. S., Stinebring, D. R., & Taylor, J. H. 1988, Nature, 

333, 237 

Harding, A. K., Usov, V. V., & Muslimov, A. G. 2005, ApJ, 622, 

Heinke, C. O., Rybicki, G. B., Narayan, R., & Grindlay, J. E. 

2006, ApJ, 644, 1090 
Hobbs, G. B., Edwards, R. T., & Manchester, R. N. 2006, 

MNRAS, 369, 655 
Kalberla, P. M. W., Burton, W. B., Hartmann, D., Arnal, E. M., 

Bajaja, E., Morras, R., & Poppel, W. G. L. 2005, A&A, 440, 


King, A. R., Beer, M. E., Rolfe, D. J., Schenker, K., & Skipp, 

J. M. 2005, MNRAS, 358, 1501 
Kuiper, L., Hermsen, W., Verbunt, P., Thompson, D. J., Stairs, 

I. H., Lyne, A. G., Strickman, M. S., k, Cusumano, G. 2000, 

A&A, 359, 615 
Lorimer, D. R., et al. 2005, MNRAS, 359, 1524 
Lorimer, D. R. and Kramer, M. 2005, Handbook of Pulsar 

Astronomy (Cambridge University Press) 
Manchester, R. N., et al. 1996, MNRAS, 279, 1235 
— . 2001, MNRAS, 328, 17 
Mattox, J. R., et al. 1996, ApJ, 461, 396 

Ransom, S. M., Eikenberry, S. S., k. Middleditch, J. 2002, AJ, 
124, 1788 

Rappaport, S., Podsiadlowski, P., Joss, P. C, DiStefano, R., & 

Han, Z. 1995, MNRAS, 273, 731 
Ray, P. S., et al. 2010, ApJ submitted. larXiv:1011.2468l 
Shklovskii, I. S. 1970, Soviet Ast., 13, 562 
Smith, D. A., et al. 2008, A&A, 492, 923 
Stappers, B. W., et al. 1996, ApJ, 465, L119 
Story, S. A., Gonthier, P. L., & Harding, A. K. 2007, ApJ, 671, 


Thompson, D. J. 2004, in Astrophysics and Space Science 
Library, Vol. 304, Cosmic Gamma-Ray Sources, ed. 
K. S. Cheng &; G. E. Romero, 149 

Zhang, L., & Cheng, K. S. 2003, A&A, 398, 639 

* |http : //nanograv ■ org | 


Ransom et al. 

Parameters for the New MSPs 

Parameter PSR J0614-3329 PSR J1231-1411 PSR J2214+3000 

Fermi BSL Association (OFGL) J0614.3-3330 J1231. 5-1410 J2214.8+3002 

Fermi 1-Year Source (IFGL) J0614.1-3328 J1231. 1-1410 J2214.8+3002 

Timing Parameters 

Right Ascension (RA, J2000) 06'^ 14™ 10!3478(3) 12^ 31™ 11!3132(7) 221^ 14™ 38!8460(1) 

Declination (DEC, J2000) -33° 29' 54'.'118(4) -14° 11' 43^63(2) +30° 00' 38'.'234(4) 

Proper Motion in RA (mas yr~l) ... ... -1.0(2)xl02 

Proper Motion in DEC (masyr-i) . . ... -3(4)xl0^ 

Pulsar Period (ms) 3.148669579439(9) 3.683878711077(3) 3.119226579079(4) 

Pulsar Frequency (Hz) 317.5944552995(9) 271.4530196103(2) 320.5922925597(4) 

Frequency Derivative (Hz s-i) -1.77(7)x lO'^^ -1.68(1) X lO^^^ -1.44(3) X 10~15 

Frequency 2nd Deriv. (Hzs"^) ... . .. 1.7(4)xl0~^^ 

Reference Epoch (MJD) 55100 55100 55100 

Dispersion Measure (pc cm'^) 37.049(1) 8.090(1) 22.557(1) 

Orbital Period (days) 53.5846127(8) 1.860143882(9) 0.416632943(5) 

Projected Semi-Major Axis (It-s) . . . 27.638787(2) 2.042633(3) 0.0590800(9) 

Orbital Eccentricity 0.0001801(1) 4(3)xl0-'' < 2 X lO"* 

Longitude of Periastron (deg) 15.92(4) 3.2(4) xlO^ 

Epoch of Periastron (MJD) 55146.821(7) 55016.8(2) 

Epoch of Ascending Node (MJD) .. . ... ... 55094.137854(2) 

Span of Timing Data (MJD) 54683-55422 54683-55430 54683-55415 

Number of 7-ray TOAs 24 32 20 

RMS 7-ray TOA Residual (fis) 99.1 24.4 110.2 

Number of radio TOAs 328 136 437 

RMS radio TOA Residual {(is) 7.1 9.3 5.0 

Derived Parameters 

Mass Function (M0) 0.007895133(3) 0.00264460(2) 1.2755(l)xl0"'' 

Min Companion Mass (Mq) >0.28 >0.19 > 0.014 

Galactic Longitude (deg) 240.50 295.53 86.86 

Galactic Latitude (deg) -21.83 48.39 -21.67 

DM-derived Distance (kpc) 1.9 0.4 1.5 

Flux Density at 820 MHz (mjy) .... 1.5 0.4 2.1 

Surface Magnetic Field (10* G) 2.4 2.9 2.1 

Characteristic Age (Gyr) 2.8 2.6 3.5 

Spin-down Lumin, E (10^"' ergss-i) 2.2 1.8 1.8 

7-ray Spectral Fit Parameters 

K (10-" phcm-2 s-i MeV-i) 2.12 ± 0.10 ± 0.13 2.62 ± 0.16 ± 0.18 0.94 ± 0.11 ± 0.05 

Spectral Index F 1.44 ± 0.05 ± 0.07 1.40 ± 0.07 ± 0.05 1.44 ± 0.13 ± 0.11 

Ecutoff (GeV) 4.49 ± 0.54 2.98 ± 0.33 tofg 2.53 ± 0.50 g 

Fioo (10-s phcm-2 s-i) 9.52 ± 0.46 ± 0.45 10.57 ± 0.62 ± 0.39 3.83 ± 0.44 ± 0.06 

Gioo (10-11 ergscm-2s-i) 10.86 ± 0.35 ± 1.06 10.33 ± 0.35 ± 0.87 3.32 ± 0.21 ± 0.24 

TS 5270.3 4798.4 958.0 

ATSeutoff 184.7 203.3 63.3 

7? (%) 210 11 49 

X-ray Parameters 

X-ray Source RA (J2000) Oe*" 14™ 10!3(3) 121^ 31™ 11!3(4) 22^ 14™ 38';84(3) 

X-ray Source Dec (J2000) -33° 29' 54(5) -14° 11' 43(6) 30° 00' 38'.'2(6) 

BB Temperature (keV) 0.23(5) 0.21(5) 0.25(4) 

i^BB,o.5-8kcV (10-1* ergs cm-2s-i) 8.7^^-4 15+5.3 2.9 

iBB,o.5-8kcV (10^0 ergss-i) 38tlr '^■^ -H 7.8 

Note. — Numbers in parentheses represent 2-a uncertainties in the last digit as determined by TEMP02 using the DE405 Solar System 
Ephemeris for the timing parameters and l-a uncertainties for the other parameters. The time system used is Barycentric Dynamical 
Time (TDB). Minimum companion masses were calculated assuming a pulsar mass of 1.4 M p) . The DM-distances were estimated using the 
NE2001 Galactic electron density model and likely have ~20% uncertainties l lCordes fc La zio 2002). The gamma-ray spectral parameters 
are from fits of exponentially cutoff power-laws as described in §2.2. Fioo and Gioo give the integrated photon or energy iiux above 0.1 GeV, 
respectively, while the last two parameters are gamma-ray detection significance of the source and significance of an exponential cutoff (as 
compared to a simple power-law), where the approximate Gaussian significance is given by ~ VTS, and TS is the Test Statistic TS = 
2A log(likelihood) between models with and without the source. The first errors are statistical and the second errors are systematic errors 
calculated from the bracketing IRFs. The 7-ray efficiency, rj = L-j/q/E = Air D^G 100/ E, assumes a beaming correction factor, /q = 1. The 
X-ray results are from Swift for PSRs J0614-3329 and J1231-1411, and Chandra for PSR J2214-I-3000. The quoted positional uncertainty 
for J2214-I-3000 is dominated by a systematic error of 0.6" in the absolute Chandra astrometry (statistical error =0.2"). X-ray fluxes and 
luminosities from 0.5— 8 keV, ^bb, 0.5-8 kcV and ^BB, 0.5-8 kcVi are based on black body fits to the point source counts using kT = 0.16 keV.