# Einstein@Home discovers a radio-quiet gamma-ray millisecond pulsar

Colin J. Clark, Holger J. Pletsch, Jason Wu, Lucas Guillemot, Matthew Kerr, Tyrel J. Johnson, Fernando Camilo, David Salvetti, Bruce Allen, David Anderson, et al.

► To cite this version:

Colin J. Clark, Holger J. Pletsch, Jason Wu, Lucas Guillemot, Matthew Kerr, et al. Einstein@Home discovers a radio-quiet gamma-ray millisecond pulsar. Sc.Adv., 2018, 4 (2), pp.eaao7228. 10.1126/sci-adv.aao7228. hal-01797367

HAL Id: hal-01797367

https://hal.science/hal-01797367

Submitted on 20 Dec 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.
# Einstein@Home discovers a radio-quiet gamma-ray millisecond pulsar

Colin J. Clark,¹,²,³* Holger J. Pletsch,¹,² Jason Wu,⁴ Lucas Guillemot,⁴,⁵,⁶ Matthew Kerr,⁷

Tyrel J. Johnson,⁸† Fernando Camilo,⁹ David Salvetti,¹⁰ Bruce Allen,¹,²,¹¹ David Anderson,¹²

Carsten Aulbert,¹,² Christian Beer,¹,² Oliver Bock,¹,² Andres Cuéllar,¹,² Heinz-Bernd Eggenstein,¹,²,¹³

Henning Fehrmann,¹,² Michael Kramer,³,⁴ Shawn A. Kwang,¹¹ Bernd Machenschalk,¹,² Lars Nieder,¹,²

Markus Ackermann,¹⁴ Marco Ajello,¹⁵ Luca Baldini,¹⁶ Jean Ballet,¹⁷ Guido Barbiellini,¹⁸,¹⁹

Denis Bastieri,²⁰,²¹ Ronaldo Bellazzini,²² Elisabetta Bissaldi,²³,²⁴ Roger D. Blandford,²⁵

Elliott D. Bloom,²⁵ Raffaella Bonino,²⁶,²⁷ Eugenio Bottacini,²⁵ Terri J. Brandt,²⁸ Johan Bregeon,²⁹

Philippe Bruel,³⁰ Rolf Buehler,¹⁴ Toby H. Burnett,³¹ Sara Buson,²⁸† Rob A. Cameron,²⁵

Regina Caputo,³² Patrizia A. Caraveo,¹⁰ Elisabetta Cavazzuti,³³ Claudia Cecchi,³⁴,³⁵ Eric Charles,²⁵

Alexandre Chekhtman,⁸‡ Stefano Ciprini,³⁴,³⁶ Lynn R. Cominsky,³⁷ Denise Costantin,²¹ Sara Cutini,³⁴,³⁶

Filippo D'Ammando,³⁸,³⁹ Andrea De Luca,¹⁰,⁴⁰ Rachele Desiante,²⁶,⁴¹ Leonardo Di Venere,²³,²⁴

Mattia Di Mauro,²⁵ Niccolò Di Lalla,¹⁶ Seth W. Digel,²⁵ Cecilia Favuzzi,²³,²⁴ Elizabeth C. Ferrara,²⁸,⁴²

Anna Franckowiak,¹⁴ Yasushi Fukazawa,⁴³ Stefan Funk,⁴⁴ Piergiorgio Fusco,²³,²⁴ Fabio Gargano,²⁴

Dario Gasparrini,³⁴,³⁶ Nico Giglietto,²³,²⁴ Francesco Giordano,²³,²⁴ Marcello Giroletti,³⁸

Germán A. Gomez-Vargas,⁴⁵,⁴⁶ David Green,²⁸,⁴⁷ Isabelle A. Grenier,¹⁷ Sylvain Guiriec,²⁸,⁴⁸

Alice K. Harding,²⁸ John W. Hewitt,⁴⁹ Deirdre Horan,³⁰ Guðlaugur Jóhannesson,⁵⁰,⁵¹ Shiki Kensei,⁴³

Michael Kuss,²² Giovanni La Mura,²¹ Stefan Larsson,⁵²,⁵³ Luca Latronico,²⁶ Jian Li,⁵⁴

Francesco Longo,¹⁸,¹⁹ Francesco Loparco,²³,²⁴ Michael N. Lovellette,⁷ Pasquale Lubrano,³⁴

Jeffrey D. Magill,⁴⁷ Simone Maldera,²⁶ Alberto Manfreda,¹⁶ Mario N. Mazziotta,²⁴ Julie E. McEnery,²⁸,⁴⁷

Peter F. Michelson,²⁵ Nestor Mirabal,²⁸† Warit Mitthumsiri,⁵⁵ Tsunefumi Mizuno,⁵⁶

Maria Elena Monzani,²⁵ Aldo Morselli,⁴⁶ Igor V. Moskalenko,²⁵ Eric Nuss,²⁹ Takashi Ohsugi,⁵⁶

Nicola Omodei,²⁵ Monica Orienti,³⁸ Elena Orlando,²⁵ Michele Palatiello,¹⁸,¹⁹ Vaidehi S. Paliya,¹⁵

Francesco de Palma,²⁴,⁵⁷ David Paneque,⁵⁸ Jeremy S. Perkins,²⁸ Massimo Persic,¹⁸,⁵⁹

Melissa Pesce-Rollins,²² Troy A. Porter,²⁵ Giacomo Principe,⁴⁴ Silvia Rainò,²³,²⁴ Riccardo Rando,²⁰,²¹

Paul S. Ray,⁷ Massimiliano Razzano,²² Anita Reimer,²⁵,⁶⁰ Olaf Reimer,²⁵,⁶⁰ Roger W. Romani,²⁵

Pablo M. Saz Parkinson,⁶¹,⁶²,⁶³ Carmelo Sgrò,²² Eric J. Siskind,⁶⁴ David A. Smith,⁶⁵ Francesca Spada,²²

Gloria Spandre,²² Paolo Spinelli,²³,²⁴ Jana B. Thayer,²⁵ David J. Thompson,²⁸ Diego F. Torres,⁵⁴,⁶⁶

Eleonora Troja,²⁸,⁴⁷ Giacomo Vianello,²⁵ Kent Wood,⁶⁷† Matthew Wood²⁵

Millisecond pulsars (MSPs) are old neutron stars that spin hundreds of times per second and appear to pulsate as their emission beams cross our line of sight. To date, radio pulsations have been detected from all rotation-powered MSPs. In an attempt to discover radio-quiet gamma-ray MSPs, we used the aggregated power from the computers of tens of thousands of volunteers participating in the Einstein@Home distributed computing project to search for pulsations from unidentified gamma-ray sources in Fermi Large Area Telescope data. This survey discovered two isolated MSPs, one of which is the only known rotation-powered MSP to remain undetected in radio observations. These gamma-ray MSPs were discovered in completely blind searches without prior constraints from other observations, raising hopes for detecting MSPs from a predicted Galactic bulge population.

## INTRODUCTION

Since its launch in 2008, the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope (1) has detected pulsations from more than 200 gamma ray-emitting rotation-powered pulsars (see http://tinyurl.com/fermipulsars). Nearly half of these are millisecond pulsars (MSPs), old neutron stars thought to have been spun up by accreting matter from a companion star (2, 3), reaching fast rotation rates of hundreds of revolutions per second (4). The first gamma-ray MSPs were found shortly after Fermi's launch by folding gamma-ray photon arrival times using rotational ephemerides from concurrent radio telescope observations (5). Dedicated radio searches targeting unidentified LAT sources (6) have also led to the discovery of more than 80 new MSPs.

However, radio pulsations from many gamma-ray pulsars are not detectable from Earth because of their low intrinsic luminosity or their radio beams not intersecting our line of sight. These radio-quiet pulsars can only be found by directly searching for pulsations at unknown frequencies in gamma-ray data. These “blind” searches in LAT data (7, 8) have discovered 57 young pulsars, 53 of which are radio quiet (9–11). However, the blind search detection of gamma-ray MSPs is challenging. Only one gamma-ray MSP (12), also detected later in radio (13), has been found in previous blind searches, and this was only possible with additional positional and orbital constraints from optical and x-ray observations (14). With this exception, all rotationally powered MSPs have been discovered through their radio pulsations, limiting knowledge of the

Copyright © 2018
The Authors, some
rights reserved;
exclusive licensee
American Association
for the Advancement
of Science. No claim to
original U.S. Government
Works. Distributed
under a Creative
Commons Attribution
NonCommercial
License 4.0 (CC BY-NC).
population of MSPs to nearby bright pulsars with radio emission beamed toward Earth (15).

Thus, in contrast to the large number of radio-quiet young pulsars seen by the LAT, no radio-quiet gamma-ray MSP has been found. Some disparity between the numbers observed in each of these two groups is expected; MSPs have wider radio beams that are visible from a larger range of viewing angles, making radio-quiet MSPs less common (16). However, the apparent absence of this source class so far is more likely to be due to the inherent difficulty of detecting the signal from an MSP at an unknown location. At higher pulsation frequencies, higher precision is required for the position-dependent “barycentering” corrections applied to photon arrival times to account for the Doppler shift due to Fermi’s motion through the solar system. The localization region of a typical unidentified gamma-ray source, limited by the LAT’s angular resolution to a few arc minutes (17), is far larger than the arc second precision required to detect gamma-ray pulsations from MSPs. Hundreds of thousands of sky locations covering the source localization region must therefore be searched, incurring a large computational cost.

In a recent study (18), we improved gamma-ray pulsar blind search methods to reduce their computational cost and increase their sensitivity. These methods are based on an efficient “semicoherent” first stage, in which only photons arriving within a short time of one another (within a “lag window”) are combined coherently. More sensitive, but more computationally expensive, follow-up stages are then used to pick out and refine weak candidate signals. The efficient search-grid spacings derived from our study (18) allowed us to use a longer

lag window and, hence, achieve higher sensitivity than was possible in previous surveys.

Using these techniques, we performed a survey (10) on the distributed volunteer computing system Einstein@Home (19) to search for pulsations from 152 pulsar-like unidentified sources (which have curved spectra and low flux variability) from the Fermi-LAT Third Source Catalog (3FGL) (17). The survey searched 5.5 years of improved “Pass 8” Fermi-LAT data (20). We searched the full gamma-ray source localization region for each source and over a range of spin frequencies and spin-down rates broad enough to cover all known young pulsars and MSPs. This search volume for each source was split into hundreds of thousands of smaller regions, each of which could be searched by a typical personal computer within a few hours. These are then distributed among the computers of tens of thousands of volunteers located across the globe, which run the search when otherwise idle. Using more than 10,000 years of volunteered central processing unit (CPU) time, this gamma-ray pulsar survey has discovered 17 young pulsars (10, 11).

RESULTS

Among the sources searched by the Einstein@Home survey were two of the three unassociated 3FGL sources with the highest detection significance, 3FGL J1035.7–6720 and 3FGL J1744.1–7619 (18). Both of these have high Galactic latitudes and pulsar-like spectra, marking them as promising MSP candidates (21). These sources were searched

<footnote n="1">¹Albert-Einstein-Institut, Max-Planck-Institut für Gravitationsphysik, D-30167 Hannover, Germany. ²Leibniz Universität Hannover, D-30167 Hannover, Germany. ³Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, University of Manchester, Manchester M13 9PL, UK. ⁴Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, D-53121 Bonn, Germany. ⁵Laboratoire de Physique et Chimie de l'Environnement et de l'Espace-Université d'Orléans/CNRS, F-45071 Orléans Cedex 02, France. ⁶Station de radioastronomie de Nançay, Observatoire de Paris, CNRS/INSU, F-18330 Nançay, France. ⁷Space Science Division, Naval Research Laboratory, Washington, DC 20375–5352, USA. ⁸College of Science, George Mason University, Fairfax, VA 22030, USA. ⁹Square Kilometre Array South Africa, Pinelands 7405, South Africa. ¹⁰Istituto Nazionale di Astrofisica (INAF)–Istituto di Astrofisica Spaziale e Fisica Cosmica Milano, via E. Bassini 15, I-20133 Milano, Italy. ¹¹Department of Physics, University of Wisconsin-Milwaukee, P.O. Box 413, Milwaukee, WI 53201, USA. ¹²Space Sciences Laboratory, University of California at Berkeley, Berkeley, CA 94720, USA. ¹³Albert-Einstein-Institut, Max-Planck-Institut für Gravitationsphysik, D-14476 Potsam-Golm, Germany. ¹⁴Deutsches Elektronen Synchrotron (DESY), D-15738 Zeuthen, Germany. ¹⁵Department of Physics and Astronomy, Clemson University, Kinard Laboratory of Physics, Clemson, SC 29634–0978, USA. ¹⁶Università di Pisa and Istituto Nazionale di Fisica Nucleare, Sezione di Pisa I-56127 Pisa, Italy. ¹⁷Laboratoire AIM, CEA-IRFU/CNRS/Université Paris Diderot, Service d'Astrophysique, CEA Saclay, F-91191 Gif sur Yvette, France. ¹⁸Istituto Nazionale di Fisica Nucleare, Sezione di Trieste, I-34127 Trieste, Italy. ¹⁹Dipartimento di Fisica, Università di Trieste, I-34127 Trieste, Italy. ²⁰Istituto Nazionale di Fisica Nucleare, Sezione di Padova, I-35131 Padova, Italy. ²¹Dipartimento di Fisica e Astronomia “G. Galilei,” Università di Padova, I-35131 Padova, Italy. ²²Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, I-56127 Pisa, Italy. ²³Dipartimento di Fisica “M. Merlin” dell’Università e del Politecnico di Bari, I-70126 Bari, Italy. ²⁴Istituto Nazionale di Fisica Nucleare, Sezione di Bari, I-70126 Bari, Italy. ²⁵W. W. Hansen Experimental Physics Laboratory, Kavli Institute for Particle Astrophysics and Cosmology, Department of Physics and SLAC National Accelerator Laboratory, Stanford University, Stanford, CA 94305, USA. ²⁶Istituto Nazionale di Fisica Nucleare, Sezione di Torino, I-10125 Torino, Italy. ²⁷Dipartimento di Fisica, Università degli Studi di Torino, I-10125 Torino, Italy. ²⁸NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA. ²⁹Laboratoire Univers and Particules de Montpellier, Université Montpellier, CNRS/IN2P3, F-34095 Montpellier, France. ³⁰Laboratoire Lepince-Ringuet, École polytechnique, CNRS/IN2P3, F-91128 Palaiseau, France. ³¹Department of Physics, University of Washington, Seattle, WA 98195–1560, USA. ³²Center for Research and Exploration in Space Science and Technology (CRESST) and NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA. ³³Italian Space Agency, Via del Politecnico snc, 00133 Roma, Italy. ³⁴Istituto Nazionale di Fisica Nucleare, Sezione di Perugia, I-06123 Perugia, Italy. ³⁵Dipartimento di Fisica, Università degli Studi di Perugia, I-06123 Perugia, Italy. ³⁶Space Science Data Center–Agenzia Spaziale Italiana, Via del Politecnico snc, I-00133, Roma, Italy. ³⁷Department of Physics and Astronomy, Sonoma State University, Rohnert Park, CA 94928–3609, USA. ³⁸INAF Istituto di Radioastronomia, I-40129 Bologna, Italy. ³⁹Dipartimento di Astronomia, Università di Bologna, I-40127 Bologna, Italy. ⁴⁰Istituto Universitario di Studi Superiori (IUS), I-27100 Pavia, Italy. ⁴¹Università di Udine, I-33100 Udine, Italy. ⁴²Center for Research and Exploration in Space Sciences and Technology II, University of Maryland, College Park, MD 20742, USA. ⁴³Department of Physical Sciences, Hiroshima University, Higashi-Hiroshima, Hiroshima 739-8526, Japan. ⁴⁴Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen Centre for Astroparticle Physics, Erwin-Rommel-Strasse 1, 91058 Erlangen, Germany. ⁴⁵Istituto de Astrofisica, Facultad de Fisica, Pontificia Universidad Católica de Chile, Casilla 306, Santiago 22, Chile. ⁴⁶Istituto Nazionale di Fisica Nucleare, Sezione di Roma “Tor Vergata,” I-00133 Roma, Italy. ⁴⁷Department of Physics and Department of Astronomy, University of Maryland, College Park, MD 20742, USA. ⁴⁸Department of Physics, The George Washington University, 725 21st Street, NW, Washington, DC 20052, USA. ⁴⁹Department of Physics, University of North Florida, 1 UNF Drive, Jacksonville, FL 32224, USA. ⁵⁰Science Institute, University of Iceland, IS-107 Reykjavik, Iceland. ⁵¹NORDITA, Roslagstullsbacken 23, 106 91 Stockholm, Sweden. ⁵²Department of Physics, KTH Royal Institute of Technology, AlbaNova, SE-106 91 Stockholm, Sweden. ⁵³The Oskar Klein Centre for Cosmoparticle Physics, AlbaNova, SE-106 91 Stockholm, Sweden. ⁵⁴Institute of Space Sciences (IEEC-CSIC), Campus UAB, Carrer de Magrans s/n, E-08193 Barcelona, Spain. ⁵⁵Department of Physics, Faculty of Science, Mahidol University, Bangkok 10400, Thailand. ⁵⁶Hiroshima Astrophysical Science Center, Hiroshima University, Higashi-Hiroshima, Hiroshima 739-8526, Japan. ⁵⁷Università Telematica Pegaso, Piazza Trieste e Trento, 48, I-80132 Napoli, Italy. ⁵⁸Max-Planck-Institut für Physik, D-80805 München, Germany. ⁵⁹Osservatorio Astronomico di Trieste, Istituto Nazionale di Astrofisica, I-34143 Trieste, Italy. ⁶⁰Institut für Astro- und Teilchenphysik and Institut für Theoretische Physik, Leopold-Franzens-Universität Innsbruck, A-6020 Innsbruck, Austria. ⁶¹Santa Cruz Institute for Particle Physics, Department of Physics and Department of Astronomy and Astrophysics, University of California at Santa Cruz, Santa Cruz, CA 95064, USA. ⁶²Department of Physics, The University of Hong Kong, Pokfulam Road, Hong Kong, China. ⁶³Laboratory for Space Research, The University of Hong Kong, Hong Kong, China. ⁶⁴NYCB Real-Time Computing Inc., Lattingtown, NY 11560–1025, USA. ⁶⁵Centre d'Etudes Nucléaires de Bordeaux Gradignan, IN2P3/CNRS, Université Bordeaux 1, BP120, F-33175 Gradignan Cedex, France. ⁶⁶Institució Catalana de Recerca i Estudis Avançats (ICREA), E-08010 Barcelona, Spain. ⁶⁷Praxis Inc., Alexandria, VA 22303, USA.</footnote>

<footnote>*Corresponding author. Email: colin.clark-2@manchester.ac.uk</footnote>

<footnote>†NASA Postdoctoral Program Fellow, USA.</footnote>

<footnote>‡Resident at Naval Research Laboratory, Washington, DC 20375, USA.</footnote>
9 and 10 times, respectively, in a survey of unidentified Fermi-LAT sources performed with the Parkes radio telescope (6), but no pulsations were detected. Previous blind gamma-ray searches of these sources, including an earlier Einstein@Home survey (22) using less data and a shorter lag window, were also unsuccessful.

The latest Einstein@Home survey revealed gamma-ray pulsations from these sources, at frequencies of $v = 348$ Hz and $v = 213$ Hz, identifying them as isolated MSPs now known as PSR J1035-6720 and PSR J1744-7619. The gamma-ray pulsations had *H* test [a significance test for pulsations in unbinned event data (23)] values of *H* = 389 and *H* = 494, corresponding to single-trial false-alarm probabilities of $< 5 \times 10^{-68}$ and $< 4 \times 10^{-86}$, respectively. For comparison, Einstein@Home can perform fewer than $10^{20}$ independent trials in 1 year.

Following their discoveries, we refined the rotational and astrometric parameters of each pulsar by timing their gamma-ray pulsations. The resulting estimates and uncertainties on their timing and spectral parameters are given in Table 1, and the resulting integrated pulse profiles are shown in Fig. 1.

These MSPs remained undetected in more sensitive analyses of the earlier Parkes survey data (6), using the now-known pulsation periods. However, a weak signal [with average flux density $S_{1400} \approx 40$ µJy (microjanskys) at a radio frequency of 1400 MHz] was detected from PSR J1035-6720 in four follow-up Parkes searches. Its estimated distance derived from its dispersion measure (DM = $84.16 \pm 0.22$ pc cm$^{-3}$) is $d \approx 1.46$ kpc or $d \approx 2.24$ kpc, according to the YMW16 (24) and NE2001 (25) Galactic electron-density models, respectively. PSR J1035-6720 therefore has a pseudo-luminosity of $L_{1400} = S_{1400}d^2 \sim 0.085$ to $0.20$ mJy (millijanskys) kpc$^2$. This is less luminous than the majority of known radio MSPs with measured pseudo-luminosities in the Australia Telescope National Facility (ATNF) Catalogue (26), as shown in Fig. 2.

In contrast, PSR J1744-7619 has remained undetected in two dedicated 3-hour follow-up searches on 19 March 2017 and 10 April 2017. These give an upper limit of $\sim 23$ µJy on its radio flux density, below the 30 µJy threshold used by Abdo et al. (9) to classify gamma-ray pulsars as radio quiet. This threshold is chosen to represent a limiting sensitivity for pulsar surveys; a new generation of radio telescopes will be needed to discover new Galactic field pulsars with flux densities below this level (27).

In the absence of a dispersion measurement, the distance to PSR J1744-7619 is hard to constrain. A weak upper limit of $d \sim 1$ kpc can be obtained by assuming that its loss in rotational energy as it spins down is converted entirely into the observed gamma-ray emission. This distance limit implies a maximum $L_{1400} \approx 0.023$ mJy kpc$^2$. Only two known MSPs have lower reliably measured pseudo-luminosities (28, 29). Even if the lack of detection in these follow-up searches was due to temporary interstellar scintillation, the upper limit from the multiple shorter observations from the previous survey suggests a maximum $L_{1400} \approx 0.15$ mJy kpc$^2$, weaker than 90% of radio MSPs. PSR J1744-7619 therefore has (at least) an unusually low radio luminosity.

The fields of PSRs J1035-6720 and J1744-7619 were observed in x-rays by XMM-Newton as part of a complementary effort to discover radio-quiet MSPs by searching for possible x-ray counterparts to unidentified LAT sources (21). Positionally coincident x-ray counterparts were detected for both pulsars. Their gamma-ray-to-x-ray flux ratios ($F_γ/F_X \approx 700$ and 1100, respectively) are consistent with those of other known isolated gamma-ray MSPs (30). The gamma-ray spectral properties of these pulsars are also consistent with those of the

<table><caption>**Table 1. Measured and derived properties of the newly discovered MSPs.**
For timing parameters, we report the mean values from a Monte Carlo timing analysis and $1\sigma$ uncertainties on the final digits in parentheses. The first set of uncertainties on spectral parameters is statistical and the second set estimates the effects of systematic uncertainties in the LAT's effective collecting area and the Galactic diffuse emission model. Timing parameters are in barycentric dynamical time (TDB) units. The spin-down power of PSR J1035-6720 has been corrected for the Doppler shift induced by its proper motion; only an upper limit is given for PSR J1744-7619, because this correction cannot be applied as a result of its uncertain distance. mas, milli-arc second.</caption><thead><tr><th>Parameter</th><th colspan="2">PSR J1035-6720 PSR J1744-7619</th></tr></thead><tbody><tr><td colspan="3"><strong>Timing parameters</strong></td></tr><tr><td>Reference time (MJD)</td><td colspan="2">55716</td></tr><tr><td>Data span (MJD)</td><td colspan="2">54682-57828</td></tr><tr><td>Right ascension (R.A.) (J2000.0), α (hh:mm:ss)</td><td>10:35:27.478(1)</td><td>17:44:00.488(2)</td></tr><tr><td>Declination (Decl.) (J2000.0), δ (dd:mm:ss)</td><td>-67:20:12.692(6)</td><td>-76:19:14.710(9)</td></tr><tr><td>Proper motion in R.A., µ<sub>α</sub> cos δ (mas year<sup>-1</sup>)</td><td>-12(3)</td><td>-21(3)</td></tr><tr><td>Proper motion in Decl. µ<sub>δ</sub> (mas year<sup>-1</sup>)</td><td>1(3)</td><td>-7(3)</td></tr><tr><td>Spin frequency, v (Hz)</td><td>348.18864014054 (8)</td><td>213.33223675351 (5)</td></tr><tr><td>Spin-down rate, -v (10<sup>-15</sup> Hz s<sup>-1</sup>)</td><td>5.633(1)</td><td>0.4405(8)</td></tr><tr><td>Second frequency derivative, |v| (10<sup>-25</sup> Hz s<sup>-2</sup>)</td><td>&lt; 1.1</td><td>&lt; 0.7</td></tr><tr><td>Spin period, P (ms)</td><td>2.8720063916972(7)</td><td>4.687524094895(1)</td></tr><tr><td>Period derivative, Ṗ (10<sup>-20</sup> s s<sup>-1</sup>)</td><td>4.647(1)</td><td>0.968(2)</td></tr><tr><td colspan="3"><strong>Derived parameters</strong></td></tr><tr><td>Galactic longitude, l (°)</td><td>290.37</td><td>317.11</td></tr><tr><td>Galactic latitude, b (°)</td><td>-7.84</td><td>-22.46</td></tr><tr><td>Spin-down power, Ė (10<sup>33</sup> erg s<sup>-1</sup>)</td><td>75.0</td><td>&lt; 3.7</td></tr><tr><td>Characteristic age, τ<sub>c</sub> (10<sup>9</sup> years)</td><td>1.0</td><td>7.7</td></tr><tr><td>Surface magnetic field, B<sub>s</sub> (10<sup>8</sup> G)</td><td>3.7</td><td>2.2</td></tr><tr><td>Light-cylinder magnetic field, B<sub>LC</sub> (10<sup>5</sup> G)</td><td>1.4</td><td>0.2</td></tr><tr><td colspan="3"><strong>Phase-averaged gamma-ray spectral parameters above 100 MeV</strong></td></tr><tr><td>Test statistic, TS</td><td>1839.2</td><td>2492.2</td></tr><tr><td>Photon index, Γ</td><td>1.46 ± 0.07 ± 0.09</td><td>1.07 ± 0.10 ± 0.02</td></tr><tr><td>Cutoff energy, E<sub>c</sub> (GeV)</td><td>2.76 ± 0.26 ± 0.36</td><td>1.82 ± 0.19 ± 0.01</td></tr><tr><td>Photon flux (10<sup>-9</sup> cm<sup>-2</sup> s<sup>-1</sup>)</td><td>24.4 ± 1.7 ± 1.5</td><td>19.2 ± 1.5 ± 1.0</td></tr><tr><td>Energy flux, F<sub>γ</sub> (10<sup>-12</sup> erg cm<sup>-2</sup> s<sup>-1</sup>)</td><td>21.5 ± 0.8 ± 1.1</td><td>20.8 ± 0.8 ± 1.1</td></tr><tr><td colspan="3"><strong>Off-pulse spectral parameters above 100 MeV</strong></td></tr><tr><td>Test statistic, TS</td><td>7.4</td><td>33.8</td></tr><tr><td>TS of exponential cutoff, TS<sub>cut</sub></td><td>—</td><td>8.3</td></tr><tr><td>Photon index, Γ</td><td>—</td><td>1.35 ± 0.70 ± 0.09</td></tr><tr><td>Cutoff energy, E<sub>c</sub> (GeV)</td><td>—</td><td>1.06 ± 0.76 ± 0.13</td></tr><tr><td>Photon flux (10<sup>-9</sup> cm<sup>-2</sup> s<sup>-1</sup>)</td><td>—</td><td>1.9 ± 0.9 ± 0.1</td></tr><tr><td>Energy flux, F<sub>γ</sub> (10<sup>-12</sup> erg cm<sup>-2</sup> s<sup>-1</sup>)</td><td>—</td><td>1.2 ± 0.3 ± 0.1</td></tr></tbody></table>
<figure><img src="image_1.png" /><figcaption>**Fig. 1.** Pulsed signals from the two newly discovered MSPs. Rotational phases of individual gamma-ray photons (bottom) and integrated pulse profiles (top) of the newly detected MSPs. Each photon has been assigned a weight, determined by its energy and arrival direction, representing the probability of it having come from the gamma-ray source in question. These weights are indicated by the color bar. In the top panels, the black dashed line shows the estimated background level, derived from the photon weights as in study of Abdo et al. (9). For PSR J1035-6720, the blue line shows the radio profile as measured by the Parkes radio telescope at 1400 MHz, averaged over four observations after subtracting an estimated baseline (dotted line). Two identical rotations are shown for clarity. MJD, Modified Julian Date.</figcaption></figure>

MSPs studied in the Second Fermi-LAT Catalog of Gamma-ray Pulsars (9). Their x-ray and gamma-ray spectra therefore give us no reason to suspect that these radio-faint MSPs belong to another class of objects different from that of other gamma-ray emitting MSPs.

One explanation for their weak radio pulsations could therefore be that their radio beams either clip our line of sight or, in the case of PSR J1744-7619, do not cross it at all. To investigate this scenario, we modeled the gamma-ray emission geometries by fitting simulated pulse profiles to the observed gamma-ray photon phases, using the fitting technique and models described in the study of Johnson et al. (31) and references therein. We considered three emission models: an outer gap (OG) model, a two-pole caustic (TPC) model, and a pair-starved polar cap (PSPC) model.

For PSR J1035-6720, only the PSPC model is consistent with the detection of radio pulsations. This model predicts that the gamma-ray pulse should lead the radio pulse, an unusual characteristic shared by only six other MSPs studied by Johnson et al. (31). The observed lag between the radio and gamma-ray pulses, shown in Fig. 1, is consistent with this prediction. However, the uncertainty of DM leads to an additional uncertainty of the arrival phase of the radio pulse of 0.15 rotations, larger than the observed lag between the two pulses. Because of this phase uncertainty and the faintness of the radio pulse, we have not attempted a joint fit of the radio and gamma-ray data.

The best-fitting models (TPC and PSPC) for PSR J1744-7619 have our line of sight cutting across the predicted radio emission cone and not merely clipping it. The lack of radio emission from PSR J1744-7619 is therefore not easily explained by small offsets of a few degrees from the best-fitting geometry. The best-fitting TPC model for this pulsar is further supported by a spectral analysis of off-pulse flux from PSR J1744-7619, which revealed likely magnetospheric unpulsed emission. Therefore, either the nondetection of radio pulsations is due to an extremely low intrinsic radio luminosity or it is at odds with the best-fitting emission models. More elaborate pulsar emission models, including emission from the current sheet (32, 33), may be required to address this tension.

## DISCUSSION

At least 12 high-confidence pulsar-like 3FGL sources at high Galactic latitudes have been identified but have remained undetected in repeated radio pulsation searches. Several convincing MSP candidates have also been highlighted by Saz Parkinson et al. (21). The most natural explanation for these sources is that they are either radio-quiet MSPs or MSPs in eclipsing binary systems, which are also difficult to detect in radio searches. Our results demonstrate that blind gamma-ray pulsation searches of LAT data are now able to find isolated MSPs
<figure><img src="image_2.png" /><figcaption>Fig. 2. Pseudo-luminosities ($L_{1400} = S_{1400} d^2$) of known MSPs (including both Galactic field MSPs and those found in globular clusters) with flux density and distance measurements in the ATNF Catalogue (26), version 1.57. The measured pseudo-luminosities of PSR J1035–6720 according to the NE2001 and YMW16 Galactic dispersion models are shown by solid and dashed vertical lines, respectively. The dotted and dashed-dotted vertical lines show the upper limits on the pseudo-luminosity of PSR J1744–7619, assuming a maximum distance of 1 kpc, from two dedicated follow-up radio observations and (a conservative estimate) from the 10 shorter observations of Camilo et al. (6), respectively.</figcaption></figure>

among these sources without requiring any additional constraints from multiwavelength observations and therefore provide an entirely new method through which to investigate the Galactic MSP population.

In particular, blind gamma-ray searches may partially address the severe biases in the observed radio MSP luminosity distribution because of the difficulty of detecting faint MSPs (15). A sizeable sample of gamma ray–selected MSPs would also provide an unbiased estimate of the relative numbers of radio-loud and radio-quiet gamma-ray MSPs, constraining the ratio of the radio and gamma-ray beam solid angles, an important discriminator for pulsar emission models (16, 34). These models suggest that a large fraction of gamma-ray MSPs could be radio quiet; the contrast with the present situation of just three gamma ray–selected MSPs compared to almost 100 radio-loud gamma-ray MSPs illustrates both the technical barrier faced by blind searches for gamma-ray MSPs until now and their large discovery potential.

Furthermore, the rotational stability and low flux variability of MSPs mean that their recoverable signals will accumulate steadily over time. Blind gamma-ray searches will therefore become sensitive to fainter and more distant MSPs as the Fermi mission continues. The rotation rates of MSPs are also unaffected by glitches or timing noise, which often interrupt the signals from young pulsars. We may therefore expect future blind gamma-ray searches to find a larger proportion of MSPs compared to previous surveys, whose discoveries have been dominated by young pulsars.

The ability to detect radio-quiet MSPs in blind gamma-ray searches is also an important tool for the search for gamma-ray emission from dark matter. A number of unidentified LAT sources at high Galactic latitudes have curved spectra that are consistent with the annihilation of a candidate for dark matter in ultra-faint dwarf spheroidal galaxies or dark matter subhalos (35, 36). However, these sources are perhaps more likely to be MSPs, whose spectra are similar. If these sources continue to remain undetected at other wavelengths, then blind gamma-ray searches are necessary to confirm or rule out at least an isolated,

radio-quiet MSP explanation. The observed GeV (billion electron volts) excess toward the Galactic center (37), another possible dark matter annihilation signal, is also consistent with an unresolved population of thousands of MSPs in the Galactic bulge (38–40). No more than a handful of these are expected to be detected with current radio telescopes because of the large distance to the bulge and the strong scattering in that region of the Galaxy (40). This latter issue does not affect gamma-ray pulsations, meaning that gamma-ray searches perhaps have the highest potential to uncover the brightest members of this population.

## MATERIALS AND METHODS

### LAT data preparation

The data in which the two MSPs were first detected covered approximately 5.5 years of LAT observations between 4 August 2008 and 6 April 2014. They were produced using preliminary internal versions of the Pass 8 instrument response functions (IRFs) and background models. Using gtselect, which is part of the Fermi Science Tools (http://fermi.gsfc.nasa.gov/ssc/data/analysis/software/), we selected SOURCE-class photons, within an 8° region of interest (ROI) around the positions of the corresponding 3FGL sources, with energies > 100 MeV and zenith angles < 100°. We excluded photons recorded when the LAT’s rocking angle was > 52° and when the LAT was not operating in normal science mode.

A photon-weighting scheme developed by Kerr (41) greatly increases the sensitivity when searching for gamma-ray pulsations by weighting the contribution of each photon to a detection statistic according to the probability of it having been emitted by the target source, derived from its reconstructed energy and arrival direction. To assign these weights, we first carried out a spectral analysis of the target sources. The initial source model consisted of all sources in the 3FGL catalog within 13° of the target source and Galactic diffuse emission and isotropic diffuse background models. With this source model, we performed a binned maximum likelihood analysis using the pointlike package (42). The spectrum of the target source was modeled with an exponentially cutoff power law (PL) with pre-factor $N_0$, spectral index $\Gamma$, and cutoff energy $E_c$

$$ \frac{dN}{dE} = N_0 \left( \frac{E}{1 \text{ GeV}} \right)^{-\Gamma} \exp \left( -\left( \frac{E}{E_c} \right) \right) \quad (1) $$

The spectral parameters of the target and all sources within 5° were free to vary in the fit, as were the normalizations of the background models. From the results of this spectral fitting, probability weights were calculated for each photon using gtsrcprob.

After discovering the pulsars, we refined their spin parameters (see the “Pulsar parameter estimation” section) using extended data sets, which included LAT photons recorded from 4 August 2008 to 16 March 2017. These data consisted of publicly available Pass 8 R2 SOURCE-class photons, processed with the P8R2_SOURCE_V6 IRFs. When producing the extended data sets, we extended the size of the ROI to 15° and reduced the zenith angle cutoff to < 90°. To produce a source model for use with gtsrcprob, another binned spectral analysis was performed using gtlike, including sources within 20° of the pulsars, and with the latest gll_iem_v06.fits map cube (43) and iso_P8R2_SOURCE_V6_v06.txt template (http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html) used to model the diffuse
<figure><img src="image_3.png" /><figcaption>**Fig. 3.** Pulse profiles of the two new gamma-ray MSPs, showing the phase definitions of the off-pulse regions (blue dashed-dotted lines) used to search for unpulsed gamma-ray emission. The estimated background levels are shown by the black dashed lines. No significant off-pulse emission was detected from PSR J1035-6720. Off-pulse emission was detected with a significance of ~5σ from PSR J1744-7619.</figcaption></figure>

The source positions were fixed at the preliminary timing positions of the pulsars.

The estimates of the spectral parameters contain systematic uncertainties due to our choice of IRFs and diffuse background models. To estimate these systematic uncertainties, we performed the same spectral analysis with rescaled effective areas (http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/Aeff_Systematics.html) and with the normalization of the Galactic diffuse emission rescaled to ± 6 % of the best-fit value.

The results of this spectral analysis for both pulsars are summarized in Table 1. The spectral properties (Γ and $E_c$) of PSRs J1035-6720 and J1744-7619 are similar to those of the MSP population seen in the Second Fermi-LAT Catalog of Gamma-ray Pulsars (9).

## LAT off-pulse analysis

The integrated pulse profile for PSR J1744-7619 includes unpulsed emission at all phases above the estimated background level. We performed a spectral analysis of this “off-pulse” emission to determine whether this is likely to be magnetospheric emission or contamination from a nearby source. We performed a similar analysis to search for unpulsed emission from PSR J1035-6720.

We defined off-pulse regions for the two pulsars using the Bayesian block decomposition method described by Scargle *et al.* (44). The definitions of the off-pulse regions are shown in Fig. 3. We first computed residual test statistic maps ($TS = 2 \Delta \log L$, where $\Delta \log L$ is the difference in log-likelihood between models with and without a putative source) with gamma-ray photons coming solely from the off-pulse region to look for putative sources around the pulsar position. No significant off-pulse emission was detected from PSR J1035-6720. This could be due to the low photon statistics in the defined off-pulse region; because the duty cycle for PSR J1035-6720 is close to 100%, a phase interval covering only ~14% of a rotation was selected for the analysis. Off-pulse emission at the position of PSR J1744-7619 was detected with $TS = 33.8$. Under the null hypothesis, the $TS$ is approximately $\chi^2$ distributed with three degrees of freedom (from the normalization, spectral index, and cutoff that are free in this test). This $TS$ therefore corresponds to a significance of ~5.0 $\sigma$. We computed the $TS$ value with both PL and exponentially cutoff PL (PLEC) models to test for curvature of the gamma-ray spectrum. The off-

pulse emission from PSR J1744-7619 shows marginal evidence of a cutoff, with $TS_{cut} = TS_{PLEC} - TS_{PL} = 8.3$ and one degree of freedom, corresponding to a significance of ~2.6$\sigma$. The positional consistency between the off-pulse source and the pulsar and the hint of a spectral cutoff suggest a magnetospheric origin for the unpulsed emission (9).

## Pulsar parameter estimation

After discovering the pulsars, we performed dedicated timing analyses using the extended data sets described above. To speed up our parameter estimation procedures, we applied a photon probability weight cutoff to remove the lowest-weighted photons. The cutoff values for each pulsar were chosen such that only 1% of the pulsation signal-to-noise ratio would be lost by this but around 90% of the lowest-weight photons would be removed. The analyses followed the timing procedure described by Clark *et al.* (10), with photon arrival times being corrected to the solar system barycenter using the JPL DE405 solar system ephemeris.

Starting with the signal parameters detected by the search (sky position, rotational frequency, and spin-down rate), we “phase-folded” the photon arrival times. We then fit template pulse profiles by minimizing the Bayesian information criterion (BIC) (45), the sum of the negative log-likelihood and a penalty factor proportional to the number of model parameters that prevents overfitting. We then varied the signal parameters using a parallelized affine-invariant Monte Carlo sampling algorithm (46, 47) to maximize the likelihood given our template pulse profile. This process continued iteratively until the likelihood stopped increasing. After each Monte Carlo step, we refolded the data with the most likely set of signal parameters, constructed a new template pulse profile, and repeated the Monte Carlo sampling.

We then added additional parameters to our phase model one by one, performed a further Monte Carlo step, and kept the new parameter in our model if its inclusion led to a decrease in the BIC. By this criterion, we were unable to measure a second frequency derivative or parallax from either pulsar. The inclusion of proper motion in our timing models led to a moderate decrease in the BIC ($\Delta$BIC = -7.8) for PSR J1035-6720 and a significant decrease ($\Delta$BIC = -40.6) for PSR J1744-7619.

Because we measured proper motion from each pulsar, a portion of the observed spin-down rate for each must be caused by the Doppler shift introduced by their velocity gaining an increasing radial component (9, 48), known as the Shklovskii effect. Correcting the observed
<figure><img src="image_4.png" /><figcaption>**Fig. 4.** Distance and spin-down constraints for the two new MSPs obtained from their gamma-ray flux, radio dispersion measure, and timing measurements of their proper motion. (Top) constraints on MSP distances, as in figure 11 of the work by Abdo et al. (9). Assuming certain gamma-ray efficiencies, constraints on the pulsar distance are inferred from the gamma-ray flux and spin-down power, after correcting the observed $\dot{v}$ for the Shklovskii effect due to the measured proper motion (1σ range given by the gray shaded region). The physically realistic region is to the left of the 100% efficiency line ($L_γ = \dot{E}$), although higher apparent efficiencies are possible, depending on the beam correction factor. Contours of constant transverse velocity are shown in gray. For PSR J1035–6720, the green dotted and dashed-dotted vertical lines show the DM distance according to the YMW16 and NE2001 models, respectively. For PSR J1744–7619, the $\dot{v}_{INT} = 0$ line shows the required distance for the observed spin down to be purely due to the Shklovskii effect, giving a hard upper limit. (Bottom) Spin-down powers after correcting $\dot{v}$ for the Shklovskii effect. The shaded region shows the allowed range for $\dot{E}$ at each distance, given the observed proper motion. The solid curves show $\dot{E}$ as a function of distance at fixed gamma-ray efficiencies. The dashed line shows the maximum $\dot{E}$, that is, if there was no proper motion.</figcaption></figure>

spin-down rates ($\dot{v}_{OBS}$) of the pulsars for the Doppler-induced apparent spin-down ($\dot{v}_{SHK}$) due to their proper motion to retrieve their intrinsic spin-down rate ($\dot{v}_{INT}$) requires knowledge of their distances, which are uncertain without additional dispersion or parallax measurements. However, we can use the observed gamma-ray fluxes ($F_γ$) and proper motions ($\mu$), and assume certain gamma-ray efficiencies ($\eta$), to retrieve constraints on the distances ($d$) by solving

$$ \dot{v}_{OBS} = \dot{v}_{INT} + \dot{v}_{SHK} = -\frac{F_{\gamma}f_{\Omega}d^2}{\eta\pi I v} - \frac{\mu^2 v d}{c} \quad (2) $$

Figure 4 shows the results of this, assuming a canonical moment of inertia $I = 10^{45}$ g cm$^2$, a gamma-ray beam correction factor of $f_Ω = 1$ (49), and a realistic range of gamma-ray efficiencies $0.01 < \eta < 1$ (converting loss of rotational energy into gamma-ray emission) (9). From the inferred distances, we can retrieve Shklovskii-corrected spin-down powers ($\dot{E}$) via $\dot{E} = F_γ4\pi f_Ω d^2/\eta$. These ranges are shown in the lower panels of Fig. 4.

From its DM distance, the gamma-ray efficiency of PSR J1035–6720 is around 7% (YMW16 model) to 17% (NE2001 model). For PSR J1744–7619, the Shklovskii effect can account for up to 65% of the apparent spin-down rate of PSR J1744–7619. The observed gamma-ray flux constrains the distance to PSR J1744–7619 to be less than ~1 kpc, assuming $\eta = 1$ and $f_Ω = 1$.

If the gamma-ray efficiency of PSR J1744–7619 is high, then its spin-down power could be as low as $1.5 \times 10^{33}$ erg s$^{-1}$, which would make this one of the least energetic gamma-ray MSPs and close to the empirical gamma-ray “death line,” that is, the minimum spin-down power an MSP must have to emit gamma rays (50).

An additional comparison can be made between the transverse velocities of known isolated MSPs and the transverse velocities given by the proper motion measurements and distance constraints described above. Lines of constant transverse velocity are shown in Fig. 4 to illustrate this. The mean (median) transverse velocity of isolated MSPs (outside of globular clusters) in the ATNF Catalogue is 84 (59) km s$^{-1}$. For PSR J1035–6720, the DM distances correspond to transverse velocities of
84 ± 20 km s⁻¹ (YMW16) and 130 ± 30 km s⁻¹ (NE2001). For PSR J1744–7619, assuming 100% gamma-ray efficiency leads to a transverse velocity of around 90 km s⁻¹, suggesting that the 1-kpc distance limit derived above is likely to be realistic for this pulsar. However, larger distances (for example, due to a lower beam correction factor) cannot be excluded on the basis of the resulting transverse velocity alone.

## Radio observations

The radio pulsar survey performed by Camilo et al. (6) targeted 56 unidentified Fermi-LAT sources with the Parkes radio telescope and detected 11 MSPs, 10 of which were new discoveries. As part of this survey, we observed the LAT sources now known to be associated with the two MSPs multiple times between 2009 and 2012, at a center frequency of 1390 MHz using an analog filter-bank system. The total power from the central beam of the Parkes multibeam receiver (full width at half maximum = 14'), filtered through 512 frequency channels spanning a bandwidth of 256 MHz, was sampled every 125 µs and recorded for off-line analysis.

Typical integration times were around 1 hour, although they ranged between 41 and 136 min for these two sources [see table 1 of Abdo et al. (7)]. The LAT source associated with PSR J1035–6720 was observed a total of nine times, at a typical offset from the actual pulsar position of 4.3'. This reduced the sensitivity of those observations to 80% of the maximum (on axis) sensitivity. The LAT source associated with PSR J1744–7619 was observed 10 times at a typical offset of 1', with no material impact on the sensitivity.

The data were analyzed using standard pulsar search techniques implemented in PRESTO (51). Before the new gamma-ray pulsars were known, the data were searched as described by Camilo et al. (6), considering possible DMs up to approximately twice the maximum DM predicted by the NE2001 model (25) for the corresponding line of sight [about 500 pc cm⁻³ for PSR J1035–6720 and 200 pc cm⁻³ for PSR J1744–7619; table 1 of Abdo et al. (7)]. Subsequently, we analyzed the data sets using the timing ephemerides obtained from gamma-ray data, folding the radio data while searching only in DM. None of the observations yielded significant radio pulsations.

The nominal threshold for these searches depends on the putative radio duty cycle of these pulsars and their DM. For an assumed duty cycle of 25% and DM < 50 pc cm⁻³, an indicative 1.4-GHz flux density threshold for PSR J1035–6720 is about 0.25 mJy normalized to a 1-hour integration. For PSR J1744–7619, the equivalent figure is approximately 0.15 mJy.

Some of the MSPs discovered by Camilo et al. (6) were not detected in all their search observations. This can be caused by a combination of orbital acceleration, eclipses, or interstellar scintillation. The first two effects are not relevant for the two new isolated MSPs. Depending on the distance of the new MSPs, especially if they are nearby objects, interstellar scintillation could very much modulate any putative radio flux density recorded by our Parkes observing system. In extreme cases, such as PSR J1514–4946 (6), this can lead to radio-loud MSPs being discoverable for less than one-third of the time. However, we know of no established isolated radio-loud MSP in the Galactic disk that is undetectable as much as 90% of the time.

Following the discoveries of the pulsars, we performed five long observations toward PSR J1035–6270 with the Parkes radio telescope. All observations were performed at a center frequency of 1369 MHz and sampled voltages were processed online with a digital polyphase filter bank into 2048 channels across a bandwidth of 256 MHz, for an effective sampling time of 4 µs. The filter-bank outputs were detected

and folded online into 256-bin profiles using the known spin ephemeris. Before each observation, we observed a pulsed noise diode to facilitate polarization and flux density (via archival observation of the radio galaxy Hydra A) calibration. The first four observations used the wideband “H-OH” receiver, whereas the latter used the center pixel of the “multibeam” receiver; for consistency, we report results only from the first four observations.

We searched the first observation over a range of DMs and found a promising signal with a signal-to-noise ratio of ~10 and a duty cycle of about 10% at a DM of 84.2 pc cm⁻³. The signal appears at the same DM and pulse phase in the observations on 24 August 2016 and, particularly, 26 September 2016, confirming the detection. Using a Gaussian template, we maximized the likelihood for the data from 09 June 2016, 24 August 2016, and 26 September 2016 to estimate an optimal DM of 84.16 ± 0.22 pc cm⁻³. With this template and DM, we measured the average flux density in each observation. The results are summarized in Table 2. (The pulsed emission in the observation on 07 August 2016 is too weak to be detected independently, but with knowledge of the phase and profile shape, the flux density can be estimated.)

For PSR J1744–7619, we performed two long observations using an identical observational setup and the multibeam system. We performed a similar search over DM (up to 800 pc cm⁻³) and found no pulsed emission. Restricting consideration to DM < 200 pc cm⁻³, we smoothed the fully averaged profile at each trial DM with a top hat of width 12.5% and recorded the maximum flux density. The results depend modestly on the baseline subtraction method used, and taking the average of two methods gives 95% confidence flux density upper limits of 0.031 and 0.032 mJy, respectively, for the two observations and 0.023 mJy for their coadded sum.

Combined with our estimated distance upper limit of *d* < 1 kpc, we found an upper limit on the pseudo-luminosity of PSR J1744–7619 of *L*<sub>1400</sub> < 0.023 mJy kpc². Only two known MSPs with both reliable flux density measurements and DM-independent distance estimates have lower pseudo-luminosities, PSR J1400–1431 with *L*<sub>1400</sub> = 0.01 mJy kpc² (28) and PSR J2322–2650 with *L*<sub>1400</sub> = 0.008 mJy kpc² (29). Of the 164 remaining MSPs in the ATNF Catalogue with reported flux density measurements and distance estimates (including those inferred from DM measurements), only two have apparent pseudo-luminosities below that of PSR J1744–7619. For one of these, PSR J0922–52, the reported flux density is only a lower bound because of an unknown pointing offset in its discovery observation (52). For the other, PSR J1745–0952, the YMW16 distance estimate (0.23 kpc) strongly disagrees with that from NE2001

<table><caption>Table 2. Summary of dedicated follow-up radio observations with the Parkes radio telescope.</caption><thead><tr><th>Target</th><th>Date</th><th>Duration (min)</th><th>Receiver</th><th>S<sub>1400</sub> (mJy)</th></tr></thead><tbody><tr><td rowspan="5">PSR J1035–6720</td><td>09 June 2016</td><td>188</td><td>H-OH</td><td>0.052 ± 0.0064</td></tr><tr><td>07 August 2016</td><td>164</td><td>H-OH</td><td>0.023 ± 0.0064</td></tr><tr><td>24 August 2016</td><td>169</td><td>H-OH</td><td>0.039 ± 0.0061</td></tr><tr><td>26 September 2016</td><td>232</td><td>H-OH</td><td>0.045 ± 0.005</td></tr><tr><td>19 March 2017</td><td>180</td><td>Multibeam</td><td>—</td></tr><tr><td rowspan="2">PSR J1744–7619</td><td>19 March 2017</td><td>180</td><td>Multibeam</td><td>&lt; 0.031</td></tr><tr><td>10 April 2017</td><td>163</td><td>Multibeam</td><td>&lt; 0.032</td></tr></tbody></table>
(1.8 kpc). This pulsar lies toward the Galactic Loop I bubble, within which the YMW16 distance for the closest pulsar with an independently estimated distance (PSR J1744–1134, angular separation $Δθ = 1.7°$ from PSR J1745–0952) was seen to be underestimated by a factor of 2.6 [parallax distance $d = 0.395$ kpc and DM distance $d_{DM} = 0.148$ kpc; (24)]. In combination with the extremely low apparent pseudo-luminosity, this suggests that the DM distance estimate is unreliable in this case.

## X-ray observations

The fields of PSRs J1035–6720 and J1744–7619 were observed by XMM-Newton for 25 ks (obsids 692830101 and 692830201) with the aim of detecting x-ray counterparts to possible radio-quiet MSPs in unidentified LAT sources (21). These gamma-ray sources were selected before the pulsation discoveries, because their gamma-ray properties indicated a probable pulsar nature.

Plausible x-ray counterparts for both pulsars were detected with a significance greater than 10σ at locations consistent with the newly discovered pulsars' timing positions. Their x-ray unabsorbed flux in the 0.3- to 10-keV energy range is $3.06_{-0.50}^{+0.96} \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$ for PSR J1035–6720 and $1.92_{-0.39}^{+0.59} \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$ for PSR J1744–7619. See table 11 of Saz Parkinson et al. (21) for additional information. We computed the probability that the association between the x-ray source and the MSP is due to a chance coincidence using $P = 1 - \exp(\pi r^2)$, where $r$ is the matching radius (in our case, the x-ray source error radius) and $\rho$ is the density of x-ray objects in the XMM-Newton EPIC (European Photon Imaging Camera) field, regardless of their flux. We estimated $P \sim 1.3 \times 10^{-4}$ for PSR J1035–6720 ($r = 2''$ and $\rho \sim 0.038 \text{ arc min}^{-2}$) and $P \sim 1.4 \times 10^{-4}$ for PSR J1744–7619 ($r = 1.9''$ and $\rho \sim 0.046 \text{ arc min}^{-2}$), which make a chance positional coincidence unlikely. In addition, the gamma-ray-to-x-ray flux ratios of the two likely x-ray counterparts (~700 and ~1100) were consistent with an MSP nature (30), confirming their associations with PSRs J1035–6720 and J1744–7619.

## Gamma-ray pulse profile modeling

To model the gamma-ray emission geometry of the MSPs, we fit simulated pulse profiles to the observed photon phases, using the fitting technique described by Johnson et al. (31). We considered three emission models: an OG model, a TPC model, and a PSPC model. These are also described in the study of Johnson et al. (31) and references therein, and briefly summarized here. In the first two models, particle acceleration takes place in narrow “gaps” in the magnetosphere, where the plasma charge density deviates from the force-free configuration. In both models, these gaps border the last-closed magnetic field lines. In the OG model, the lower bound of the gap is defined by the “null-charge surface,” where the plasma charge density changes sign. In the TPC model, the gap begins at the pulsar surface and extends to the light cylinder. The PSPC model is valid for low-$\dot{E}$ pulsars, where pair creation may be insufficient to reach the force-free configuration, allowing for particle acceleration throughout the regions of open field lines. We assumed a hollow-cone model (16) for the radio beams.

We used simulated rotation periods of 2.5 and 4.5 ms, respectively, and a simulated period derivative of $1 \times 10^{-20} \text{ s s}^{-1}$ for both. Instead of a Poisson likelihood, we used a $\chi^2$ statistic to fit the weighted-counts light curves, using 60 bins for both MSPs, with background levels and uncertainties calculated as in the study of Abdo et al. (9). The best-fitting parameters (the angle between spin axis and magnetic dipole, $\alpha$; the angle between the spin axis and the line-of-sight, $\zeta$; and the estimated beam correction factor, $f_\Omega$) for each model are given in Table 3, with 95% confidence-level uncertainties estimated as described by Johnson et al.

(31). The best-fitting light curves from all models for each pulsar are shown in Fig. 5. It should be noted that, as pointed out by Pierbattista et al. (53), fitting only the gamma-ray light curves with these toy models can lead to systematic biases in the best-fit parameters of ~10°. We also note that no model considered here successfully reproduces the broad double-peaked pulse profile of PSR J1035–6720.

<table><caption>**Table 3. Best-fit parameters from gamma-ray pulse profile modeling.**
For each pulsar, we report the best-fitting magnetic inclination angles ($\alpha$), viewing angles ($\zeta$), and beam correction factors ($f_\Omega$), according to the OG, TPC, and PSPC models.</caption><thead><tr><th rowspan="2">Parameter</th><th colspan="3">PSR J1035–6720</th><th colspan="3">PSR J1744–7619</th></tr><tr><th>OG</th><th>TPC</th><th>PSPC</th><th>OG</th><th>TPC</th><th>PSPC</th></tr></thead><tbody><tr><td>Model</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>-ln(L)</td><td>84.5</td><td>78.9</td><td>72.0</td><td>84.6</td><td>58.1</td><td>72.7</td></tr><tr><td>α (deg)</td><td>9<sup>+1</sup><sub>-2</sub></td><td>7<sup>+2</sup><sub>-1</sub></td><td>51<sup>+4</sup><sub>-3</sub></td><td>67<sup>+2</sup><sub>-2</sub></td><td>62<sup>+1</sup><sub>-2</sub></td><td>62<sup>+5</sup><sub>-1</sub></td></tr><tr><td>ζ (deg)</td><td>75<sup>+1</sup><sub>-1</sub></td><td>70<sup>+1</sup><sub>-1</sub></td><td>78<sup>+3</sup><sub>-3</sub></td><td>32<sup>+2</sup><sub>-2</sub></td><td>44<sup>+3</sup><sub>-3</sub></td><td>74<sup>+3</sup><sub>-3</sub></td></tr><tr><td>f<sub>Ω</sub></td><td>0.22<sup>+0.25</sup><sub>-0.04</sub></td><td>0.62<sup>+0.08</sup><sub>-0.11</sub></td><td>0.98<sup>+0.12</sup><sub>-0.04</sub></td><td>0.94<sup>+0.05</sup><sub>-0.10</sub></td><td>0.89<sup>+0.03</sup><sub>-0.09</sub></td><td>1.14<sup>+0.15</sup><sub>-0.09</sub></td></tr></tbody></table>

<figure><img src="image_5.png" /><figcaption>**Fig. 5. Gamma-ray pulse profiles of the newly detected MSPs.** The overlaying curves are the best-fitting pulse profiles predicted by fits to OG, TPC, and PSPC gamma-ray emission models. The black dashed line is the estimated background level, derived from the photon weights as in the study of Abdo et al. (9). Predicted radio pulse profiles (with arbitrary baseline and normalization) are shown by dashed-dotted lines. The three models are shown in separate rotations for clarity.</figcaption></figure>
The PSPC model for PSR J1035-6720 predicted a phase offset between the gamma-ray and radio pulses that was reasonably consistent with that observed. However, the uncertainty in DM led to an additional uncertainty in the radio pulse phase of 0.15 rotations, precluding a joint fit of the radio and gamma-ray pulse profiles as is normally possible for other gamma-ray pulsars.

For the models that did not predict radio emission (TPC and OG for PSR J1035-6720 and OG for PSR J1744-7619), the lines of sight barely missed the edges of the emission cones at 300 MHz, suggesting that these MSPs might appear brighter at even lower frequencies.

## REFERENCES AND NOTES

1. W. B. Atwood, A. A. Abdo, M. Ackermann, W. Althouse, B. Anderson, M. Axelsson, L. Baldini, J. Ballet, D. L. Band, G. Barbiellini, J. Bartelt, D. Bastieri, B. M. Baughman, K. Bechtol, D. Bedéréde, F. Bellardi, R. Bellazzini, B. Bereni, G. F. Bignami, D. Bisello, E. Bissaldi, R. D. Blandford, E. D. Bloom, J. R. Bogart, E. Bonamente, J. Bonnell, A. W. Borgland, A. Bouvier, J. Bregeon, A. Brez, M. Brigida, P. Bruel, T. H. Burnett, G. Busetto, G. A. Caliandro, R. A. Cameron, P. A. Caraveo, S. Carius, P. Carlson, J. M. Casandjian, E. Cavazzuti, M. Ceccanti, C. Cecchi, E. Charles, A. Chekhtman, C. C. Cheung, J. Chiang, R. Chipaux, A. N. Gillis, S. Ciprini, R. Claus, J. Cohen-Danugi, S. Condamoor, J. Conrad, R. Corbet, L. Corucci, L. Costamonic, S. Cutini, D. S. Davis, D. Decotigny, M. DeKlotz, C. D. Dermer, A. de Angelis, S. W. Digel, E. do Couto e Silva, P. S. Drell, R. Dubois, D. Dumora, Y. Edmonds, D. Fabiani, C. Farnier, C. Favuzzi, D. L. Flath, P. Fleury, W. B. Focke, S. Funk, P. Fusco, F. Gargano, D. Gasparrini, N. Gehrels, F.-X. Gentit, S. Germani, B. Giebels, N. Giglietto, P. Giommi, F. Giordano, T. Glanzman, G. Godfrey, I. A. Grenier, M.-H. Grondin, J. E. Grove, L. Guillemot, S. Guiriec, G. Haller, A. K. Harding, P. A. Hart, E. Hays, S. E. Healey, M. Hirayama, L. Hjalmsdorfer, R. Horn, R. E. Hughes, G. Johansson, G. Johansson, A. S. Johnson, R. P. Johnson, T. J. Johnson, W. N. Johnson, T. Kamae, H. Katagiri, J. Kataoka, A. Kavelaars, N. Kawai, H. Kelly, M. Kerr, W. Klamra, J. Knödlseder, M. L. Kocian, N. Komin, F. Kuehn, M. Kuss, D. Landriu, L. Latronico, B. Lee, S.-H. Lee, M. Lemoine-Goumard, A. M. Lionetto, F. Longo, F. Loparco, B. Lott, M. M. Lovellette, P. Lubrano, G. M. Madejski, A. Makeev, B. Marangelli, M. M. Massai, M. N. Mazziotta, J. E. McEnery, N. Menon, C. Meurer, P. F. Michelson, M. Minuti, N. Mirizzi, W. Mitthumsiri, T. Mizuno, A. A. Moiseev, Monte, P. E. Monzani, E. Moretti, A. Morselli, I. V. Moskalenko, S. Murgia, T. Nakamori, S. Nishino, P. L. Nolan, J. P. Norris, E. Nuss, M. Ohno, T. Ohsugi, N. Omodei, E. Orlando, J. F. Ormes, M. Paccagnella, D. Picque, J. H. Panetta, D. Parent, M. Pearce, M. Pepe, A. Pezzo, M. PaccRollins, P. Picozza, L. Pieri, M. Pinchera, F. Piron, T. A. Porter, L. Poupard, S. Rainò, R. Rando, E. Rapposelli, M. Razzano, A. Reimer, O. Reimer, T. Reposeur, L. C. Reyes, S. Ritz, L. S. Rochester, A. Y. Rodriguez, R. W. Romano, M. Roth, J. J. Russell, F. Ryde, S. Sabatini, F.-W. Sadrozinski, D. Sanchez, A. Sander, L. Sapoznikov, P. M. Saz Parkinson, J. D. Scargle, T. L. Schalk, G. Scolieri, C. Sgrò, G. H. Share, M. Shaw, T. Shimokawabe, C. Shrader, A. Sierpowska-Tartosik, E. J. Siskind, D. A. Smith, P. D. Smith, G. Spandre, P. Spinelli, J.-L. Starck, T. E. Stephens, M. S. Strickman, A. W. Strong, D. J. Suson, H. Tajima, H. Takahashi, T. Takahashi, T. Tanaka, A. Tenze, S. Tether, J. B. Thayer, J. G. Thayer, D. J. Thompson, L. Tibaldo, O. Tibolla, D. F. Torres, G. Tosti, A. Tramacere, Y. Uchiyama, T. L. Usher, A. Van Etten, V. Vasileiou, N. Vilchez, V. Vitale, A. P. Waite, P. Wang, K. Watters, B. L. Winer, K. S. Wood, T. Ylinen, M. Ziegler, The Large Area Telescope on the Fermi Gamma-ray Telescope mission. Astrophys. J. 697, 1071 (2009).

2. M. A. Alpar, A. F. Cheng, M. A. Ruderman, J. Shaham, A new class of radio pulsars. Nature 300, 728-730 (1982).

3. D. Bhattacharyya, E. P. J. van den Heuvel, Formation and evolution of binary and millisecond radio pulsars. Phys. Rep. 203, 1-124 (1991).

4. J. W. T. Hessels, S. M. Ransom, I. H. Stairs, P. C. C. Freire, V. M. Kaspi, F. Camilo, A radio pulsar spinning at 716 Hz. Science 311, 1901-1904 (2006).

5. A. A. Abdo, M. Ackermann, M. Ajello, W. B. Atwood, M. Axelsson, L. Baldini, J. Ballet, G. Barbiellini, M. G. Baring, D. Bastieri, B. M. Baughman, K. Bechtol, R. Bellazzini, B. Bereni, G. F. Bignami, R. D. Blandford, E. D. Bloom, E. Bonamente, A. W. Borgland, J. Bregeon, A. Brez, M. Brigida, P. Bruel, T. H. Burnett, G. A. Caliandro, R. A. Cameron, F. Camilo, P. A. Caraveo, P. Carlson, J. M. Casandjian, C. Cecchi, Ö. Çelik, E. Charles, A. Chekhtman, C. C. Cheung, J. Chiang, S. Ciprini, R. Claus, I. Cognard, J. Cohen-Tanugi, L. R. Cominsky, J. Conrad, R. Corbet, S. Cutini, C. D. Dermer, G. Desvignes, A. de Angelis, A. de Luca, F. de Palma, S. W. Digel, M. Dormody, E. do Couto e Silva, P. S. Drell, R. Dubois, D. Dumora, Y. Edmonds, C. Farnier, C. Favuzzi, S. J. Fegan, W. B. Focke, M. Frailis, P. C. C. Freire, Y. Fukazawa, S. Funk, P. Fusco, F. Gargano, D. Gasparrini, N. Gehrels, S. Germani, B. Giebels, N. Giglietto, F. Giordano, T. Glanzman, G. Godfrey, I. A. Grenier, M. H. Grondin, J. E. Grove, L. Guillemot, S. Guiriec, T. Y. Hanabata, A. K. Harding, M. Hayashida, E. Hays, G. Hobbs, R. E. Hughes, G. Jóhannesson, A. S. Johnson, R. P. Johnson, T. J. Johnson, W. N. Johnson, S. Johnston, T. Kamae, H. Katagiri, J. Kataoka, N. Kawai, M. Kerr, J. Knödlseder, M. L. Kocian, M. Kramer, M. Kuss, J. Lande, L. Latronico, M. Lemoine-Goumard, F. Longo, F. Loparco,

B. Lott, M. N. Lovellette, P. Lubrano, G. M. Madejski, A. Makeev, R. N. Manchester, M. Marelli, M. N. Mazziotta, W. McConville, J. E. McEnery, M. A. McLaughlin, C. Meurer, P. F. Michelson, M. W. Mitthumsiri, T. Mizuno, A. A. Moiseev, C. Monte, M. E. Monzani, A. Morselli, I. V. Moskalenko, S. Murgia, P. L. Nolan, J. P. Norris, E. Nuss, T. Ohsugi, N. Omodei, E. Orlando, J. F. Ormes, D. Paneque, J. H. Panetta, D. Parent, V. Pelassa, M. Pepe, M. Pesce-Rollins, F. Piron, T. A. Porter, S. Rainò, R. Rando, S. M. Ransom, P. S. Ray, M. Razzano, N. Rea, O. Reimer, T. Reposeur, S. Ritz, L. S. Rochester, A. Y. Rodriguez, R. W. Romano, M. Roth, F. Ryde, H. F. W. Sadrozinski, D. Sanchez, A. Sander, P. M. Saz Parkinson, J. D. Scargle, T. L. Schalk, C. Sgrò, E. J. Siskind, D. A. Smith, P. D. Smith, G. Spandre, P. Spinelli, B. W. Stappers, J. L. Starck, E. Striani, M. S. Strickman, D. J. Suson, H. Tajima, H. Takahashi, T. Tanaka, J. B. Thayer, J. G. Thayer, G. Theureau, D. J. Thompson, S. E. Thorsett, L. Tibaldo, D. F. Torres, G. Tosti, A. Tramacere, Y. Uchiyama, T. L. Usher, A. Van Etten, V. Vasileiou, C. Venter, N. Vilchez, V. Vitale, A. P. Waite, E. Wallace, P. Wang, K. Watters, N. Webb, P. Weltevrede, B. L. Winer, K. S. Wood, T. Ylinen, M. Ziegler, A population of gamma-ray millisecond pulsars seen with the Fermi Large Area Telescope. Science 325, 848-852 (2009).

6. F. Camilo, M. Kerr, P. S. Ray, S. M. Ransom, J. Sarkissian, H. T. Cromartie, S. Johnston, J. E. Reynolds, M. T. Wolff, P. C. C. Freire, B. Bhattacharyya, E. C. Ferrara, M. Keith, P. F. Michelson, P. M. Saz Parkinson, K. S. Wood, P. Carver, SPARCS radio searches of Fermi gamma-ray sources and millisecond pulsar discoveries. Astrophys. J. 810, 85 (2015).

7. A. A. Abdo, M. Ackermann, M. Ajello, B. Anderson, W. B. Atwood, M. Axelsson, L. Baldini, J. Ballet, G. Barbiellini, M. G. Baring, D. Bastieri, B. M. Baughman, K. Bechtol, R. Bellazzini, B. Bereni, G. F. Bignami, R. D. Blandford, E. D. Bloom, E. Bonamente, A. W. Borgland, J. Bregeon, A. Brez, M. Brigida, P. Bruel, T. H. Burnett, G. A. Caliandro, R. A. Cameron, P. A. Caraveo, J. M. Casandjian, C. Cecchi, Ö. Çelik, A. Chekhtman, C. C. Cheung, J. Chiang, S. Ciprini, R. Claus, J. Cohen-Tanugi, J. Conrad, S. Cutini, C. D. Dermer, A. de Angelis, A. de Luca, F. de Palma, S. W. Digel, M. Dormody, E. do Couto e Silva, P. S. Drell, R. Dubois, D. Dumora, C. Farnier, C. Favuzzi, S. J. Fegan, Y. Fukazawa, S. Funk, P. Fusco, F. Gargano, D. Gasparrini, N. Gehrels, S. Germani, B. Giebels, N. Giglietto, P. Giommi, F. Giordano, T. Glanzman, G. Godfrey, I. A. Grenier, M.-H. Grondin, J. E. Grove, L. Guillemot, S. Guiriec, C. Gwon, Y. Hanabata, A. K. Harding, M. Hayashida, E. Hays, R. E. Hughes, G. Jóhannesson, R. P. Johnson, T. J. Johnson, W. N. Johnson, T. Kamae, H. Katagiri, J. Kataoka, N. Kawai, M. Kerr, J. Knödlseder, M. L. Kocian, M. Kuss, J. Lande, L. Latronico, M. Lemoine-Goumard, F. Longo, F. Loparco, B. Lott, M. M. Lovellette, P. Lubrano, G. M. Madejski, A. Makeev, M. Marelli, M. N. Mazziotta, W. McConville, J. E. McEnery, C. Meurer, P. F. Michelson, M. W. Mitthumsiri, T. Mizuno, Monte, M. E. Monzani, A. Morselli, I. V. Moskalenko, S. Murgia, P. L. Nolan, J. P. Norris, E. Nuss, T. Ohsugi, N. Omodei, E. Orlando, J. F. Ormes, D. Paneque, D. Parent, V. Pelassa, M. Pepe, M. Pesce-Rollins, M. Pierbattista, F. Piron, T. A. Porter, J. R. Primack, S. Rainò, R. Rando, P. S. Ray, M. Razzano, N. Rea, A. Reimer, O. Reimer, T. Reposeur, S. Ritz, L. S. Rochester, A. Y. Rodriguez, R. W. Romano, F. Ryde, H. F.-W. Sadrozinski, D. Sanchez, A. Sander, P. M. Saz Parkinson, J. D. Scargle, C. Sgrò, E. J. Siskind, D. A. Smith, P. D. Smith, G. Spandre, P. Spinelli, J.-L. Starck, M. S. Strickman, D. J. Suson, H. Tajima, H. Takahashi, T. Takahashi, T. Tanaka, J. G. Thayer, D. J. Thompson, L. Tibaldo, O. Tibolla, D. F. Torres, G. Tosti, A. Tramacere, Y. Uchiyama, T. L. Usher, A. Van Etten, V. Vasileiou, N. Vilchez, V. Vitale, A. P. Waite, P. Wang, K. Watters, B. L. Winer, M. T. Wolff, K. S. Wood, T. Ylinen, M. Ziegler, Detection of 10 gamma-ray pulsars through blind frequency searches using the Fermi LAT. Science 325, 840-844 (2009).

8. H. J. Pletsch, L. Guillemot, B. Allen, M. Kramer, C. Aulbert, F. H. Fehrmann, P. S. Ray, E. D. Barr, A. Belfiore, F. Camilo, P. A. Caraveo, Ö. Çelik, D. J. Champion, M. Dormody, R. P. Eatough, E. C. Ferrara, P. C. C. Freire, M. Keith, M. Kerr, A. de Luca, A. G. Lyne, M. Marelli, M. A. McLaughlin, D. Parent, S. M. Ransom, M. Razzano, W. Reich, P. M. Saz Parkinson, B. W. Stappers, M. T. Wolff, Discovery of nine gamma-ray pulsars in Fermi Large Area Telescope data using a new blind search method. Astrophys. J. 744, 105 (2012).
9. A. A. Abdo, M. Ajello, A. Allafort, L. Baldini, J. Ballet, G. Barbiellini, M. G. Baring, D. Bastieri, A. Belfiore, R. Bellazzini, B. Bhattacharyya, E. Bissaldi, E. D. Bloom, E. Bonamente, E. Bottacini, T. J. Brandt, J. Bregeon, M. Brigida, P. Bruel, R. Buehler, M. Burgay, T. H. Burnett, G. Busetto, S. Buson, G. A. Caliandro, R. A. Cameron, F. Camilo, P. A. Caraveo, J. M. Casandjian, C. Cecchi, Ö. Çelik, E. Charles, S. Chaty, R. C. G. Chaves, A. Chekhtman, A. W. Chen, J. Chiang, G. Chiaro, S. Ciprini, R. Claus, I. Cognard, J. Cohen-Tanugi, L. R. Cominsky, J. Conrad, S. Cutini, F. D'Ammando, A. de Angelis, M. E. DeCesar, A. De Luca, P. R. den Hartog, F. de Palma, C. D. Dermer, G. Desvignes, S. W. Digel, L. Di Venere, P. S. Drell, A. Drlica-Wagner, R. Dubois, D. Dumora, C. M. Espinoza, L. Falletti, C. Favuzzi, E. C. Ferrara, W. B. Focke, A. Franckowiak, P. C. C. Freire, S. Funk, P. Fusco, F. Gargano, D. Gasparrini, S. Germani, N. Giglietto, P. Giommi, F. Giordano, M. Giroletti, T. Glanzman, G. Godfrey, E. V. Gotthelf, I. A. Grenier, M.-H. Grondin, J. E. Grove, L. Guillemot, S. Guiriec, D. Hadasch, Y. Hanabata, A. K. Harding, M. Hayashida, E. Hays, J. Hessels, J. Hewitt, A. B. Hill, D. Horan, X. Hou, R. E. Hughes, M. S. Jackson, G. H. Janssen, T. Jogler, G. Jóhannesson, R. P. Johnson, A. S. Johnson, T. J. Johnson, W. N. Johnson, S. Johnston, T. Kamae, J. Kataoka, M. Keith, M. Kerr, J. Knödlseder, M. Kramer, M. Kuss, J. Lande, S. Larsson, L. Latronico, M. Lemoine-Goumard, F. Longo, F. Loparco, M. N. Lovellette, P. Lubrano, A. G. Lyne, R. N. Manchester, M. Marelli, F. Massaro, M. Mayer, M. N. Mazziotta, J. E. McEnery, M. A. McLaughlin, J. Mehault, P. F. Michelson, R. P. Mignani, W. Mitthumsiri, T. Mizuno, A. A. Moiseev, M. E. Monzani, A. Morselli, I. V. Moskalenko, S. Murgia, T. Nakamori,
R. Nemmen, E. Nuss, M. Ohno, T. Ohsugi, M. Orienti, E. Orlando, J. F. Ormes, D. Paneque, J. H. Panetta, D. Parent, J. S. Perkins, M. Pesce-Rollins, M. Pierbattista, F. Piron, G. Pivato, H. J. Pletsch, T. A. Porter, A. Possenti, S. Rainò, R. Rando, S. M. Ransom, S. P. Ray, M. Razzano, N. Rea, A. Reimer, O. Reimer, N. Renault, T. Reposeur, S. Ritz, R. W. Romani, M. Roth, R. Rousseau, J. Roy, J. Ruan, A. Sartori, P. M. Saz Parkinson, J. D. Scargle, A. Schulz, C. Sgrò, R. Shannon, E. J. Siskind, D. A. Smith, G. Spandre, P. Spinelli, B. W. Stappers, A. W. Strong, D. J. Suson, H. Takahashi, J. G. Thayer, J. B. Thayer, G. Theureau, D. J. Thompson, S. E. Thorsett, L. Tibaldo, O. Tibolla, M. Tinivella, D. F. Torres, G. Tosti, E. Troja, Y. Uchiyama, T. L. Usher, J. Vandenbroucke, V. Vasileiou, C. Venter, G. Vianello, V. Vitale, N. Wang, P. Weltevrede, B. L. Winer, M. T. Wolff, D. L. Wood, K. S. Wood, M. Wood, Z. Yang, The second Yang Large Area Telescope catalog of gamma-ray pulsars. Astrophys. J. Suppl. Ser. 208, 17 (2013).

10. C. J. Clark, J. Wu, H. J. Pletsch, L. Guillemot, B. Allen, C. Aulbert, C. Beer, O. Bock, A. Cuéllar, H. B. Eggenstein, H. Fehrmann, M. Kramer, B. Machenschalk, L. Nieder, The Einstein@Home gamma-ray pulsar survey. I. Search methods, sensitivity, and discovery of new young gamma-ray pulsars. Astrophys. J. 834, 106 (2017).

11. C. J. Clark, "Blind-search Survey for Gamma-ray Pulsars," thesis, Leibniz Universität Hannover (2017).

12. H. J. Pletsch, L. Guillemot, H. Fehrmann, B. Allen, M. Kramer, C. Aulbert, M. Ackermann, M. Ajello, A. de Angelis, W. B. Atwood, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, K. Bechtol, R. Bellazzini, A. W. Borgland, E. Bottacini, T. J. Brandt, J. Bregeon, M. Brigida, P. Bruel, R. Buehler, S. Buson, G. A. Caliandro, R. A. Cameron, P. A. Caraveo, J. M. Casandjian, C. Cecchi, Ö. Çelik, E. Charles, R. C. G. Chaves, C. C. Cheung, J. Chiang, S. Ciprini, R. Claus, J. Cohen-Tanugi, J. Conrad, S. Cutini, F. D'Ammando, C. D. Dermer, S. W. Digel, P. S. Drell, A. Drlica-Wagner, R. Dubago, N. Dumora, C. Favuzzi, E. C. Ferrara, A. Franckowiak, Y. Fukazawa, P. Fusco, F. Gargano, N. Gehrels, S. Germani, N. Giglietto, F. Giordano, M. Giroletti, G. Godfrey, I. A. Grenier, M.-H. Grondin, J. E. Grove, S. Guiriec, D. Hadasch, Y. Hanabata, A. K. Harding, P. R. den Grind, M. Hayashida, E. Hays, A. B. Hill, X. Hou, R. E. Hughes, G. Johannesson, M. S. Jackson, T. Jogler, A. S. Johnson, W. N. Johnson, J. Kataoka, M. Kerr, J. Knödlseder, M. Kuss, J. Lande, S. Larsson, L. Latronico, M. Lemoine-Goumard, F. Longo, F. Loparco, M. N. Lovellette, P. Lubrano, F. Massaro, M. M. Mazziotta, J. E. McEnery, J. Mehault, P. F. Michelson, W. Mitthumsiri, T. Mizuno, M. E. Monzani, A. Morselli, I. V. Moskalenko, S. Murgia, T. Nakamori, R. Nemmen, E. Nuss, M. Ohno, T. Ohsugi, N. Omodei, M. Orienti, E. Orlando, F. de Palma, D. Paneque, J. S. Perkins, F. Piron, G. Pivato, T. A. Porter, S. Rainò, R. Rando, P. S. Ray, M. Razzano, A. Reimer, O. Reimer, T. Reposeur, S. Ritz, R. W. Romani, C. Romoli, D. A. Sanchez, P. M. Saz Parkinson, A. Schulz, C. Sgrò, E. do Couto e Silva, E. J. Siskind, D. A. Smith, G. Spandre, P. Spinelli, D. J. Suson, H. Takahashi, T. Tanaka, J. B. Thayer, J. G. Thayer, D. J. Thompson, L. Tibaldo, M. Tinivella, E. Troja, T. L. Usher, J. Vandenbroucke, V. Vasileiou, G. Vianello, V. Vitale, A. P. Waite, B. L. Winer, K. S. Wood, M. Wood, Z. Yang, S. Zimmer, Y. millisecond pulsar discovery via gamma-ray pulsations. Science 338, 1314–1317 (2012).

13. P. S. Ray, S. M. Ransom, C. C. Cheung, M. Giroletti, I. Cognard, F. Camilo, B. Bhattacharyya, J. Roy, R. W. Romani, E. C. Ferrara, L. Guillemot, S. Johnston, M. Keith, M. Kerr, M. Kramer, H. J. Pletsch, P. M. Saz Parkinson, K. S. Wood, Radio detection of the Fermi-LAT blind search millisecond pulsar J1311–3430. Astrophys. J. 763, L13 (2013).

14. R. W. Romani, 2FGL J1311.7–3429 joins the black widow club. Astrophys. J. 754, L25 (2012).

15. D. R. Lorimer, Binary and millisecond pulsars. Living Rev. Relativ. 11, 8 (2008).

16. S. A. Story, P. L. Gonthier, A. K. Harding, Population synthesis of radio and γ-ray millisecond pulsars from the galactic disk. Astrophys. J. 671, 713 (2007).

17. F. Acero, M. Ackermann, M. Ajello, A. Albert, W. B. Atwood, M. Axelsson, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, A. Belfiore, R. Bellazzini, E. Bissaldi, R. D. Blandford, E. D. Bloom, J. R. Bogart, R. Bonino, E. Bottacini, J. Bregeon, R. J. Britto, P. Bruel, R. Buehler, T. H. Burnett, S. Buson, G. A. Caliandro, R. A. Cameron, R. Caputo, P. A. Caraveo, J. M. Casandjian, E. Cavazzuti, E. Charles, R. C. G. Chaves, A. Chekhtman, C. C. Cheung, J. Chiang, G. Chiaro, S. Ciprini, R. Claus, J. Cohen-Tanugi, L. R. Cominsky, J. Conrad, S. Cutini, F. D'Ammando, A. de Angelis, M. Deklotz, F. de Palma, R. Desiante, S. W. Digel, L. Di Venere, P. S. Drell, R. Dubois, D. Dumora, C. Favuzzi, S. J. Fegan, E. C. Ferrara, J. Finke, A. Franckowiak, Y. Fukazawa, S. Funk, P. Fusco, F. Gargano, D. Gasparrini, B. Giebels, N. Giglietto, P. Giommi, F. Giordano, M. Giroletti, T. Glanzman, G. Godfrey, I. A. Grenier, M.-H. Grondin, J. E. Grove, L. Guillemot, S. Guiriec, D. Hadasch, A. K. Harding, E. Hays, J. W. Hewitt, A. B. Hill, D. Horan, G. Irfate, T. Jogler, G. Johannesson, R. P. Johnson, A. S. Johnson, T. J. Johnson, W. N. Johnson, T. Kamae, J. Kataoka, J. Katsuta, M. Kuss, G. La Mura, D. Landriu, S. Larsson, L. Latronico, M. Lemoine-Goumard, J. Li, L. Li, F. Longo, F. Loparco, B. Lott, M. N. Lovellette, P. Lubrano, G. M. Madejski, F. Massaro, M. Mayer, M. N. Mazziotta, J. E. McEnery, P. F. Michelson, N. Mirabal, T. Mizuno, A. A. Moiseev, M. Mongelli, M. E. Monzani, A. Morselli, I. V. Moskalenko, S. Murgia, E. Nuss, M. Ohno, T. Ohsugi, N. Omodei, M. Orienti, E. Orlando, J. F. Ormes, D. Paneque, J. H. Panetta, J. S. Perkins, M. Pesce-Rollins, F. Piron, G. Pivato, T. A. Porter, J. L. Racusin, R. Rando, M. Razzano, S. Razzaque, A. Reimer, O. Reimer, T. Reposeur, L. S. Rochester, R. W. Romani, D. Salvetti, M. Sánchez-Conde, P. M. Saz Parkinson, A. Schulz, E. J. Siskind, A. A. Smith, F. Spada, G. Spandre, P. Spinelli, T. E. Stephens, A. W. Strong, D. J. Suson, H. Takahashi, T. Takahashi, Y. Tanaka, J. G. Thayer, J. B. Thayer, D. J. Thompson, L. Tibaldo, O. Tibolla, D. F. Torres, E. Torres, G. Tosti, E. Troja, B. Van Klaveren, G. Vianello, B. L. Winer, K. S. Wood,

M. Wood, S. Zimmer, Fermi Large Area Telescope third source catalog. Astrophys. J. Suppl. Ser. 218, 23 (2015).

18. H. J. Pletsch, C. J. Clark, Optimized blind gamma-ray pulsar searches at fixed computing budget. Astrophys. J. 795, 75 (2014).

19. B. Knispel, B. Allen, J. M. Cordes, J. S. Deneva, D. Anderson, C. Aulbert, N. D. R. Bhat, O. Bock, S. Bogdanov, A. Brazier, F. Camilo, D. J. Champion, S. Chatterjee, F. Crawford, P. B. Demorest, H. Fehrmann, P. C. Freire, M. E. Gonzalez, D. Hammer, J. W. T. Hessels, F. A. Jenet, L. Kasian, V. M. Kaspi, M. Kramer, P. Lazarus, J. van Leeuwen, D. R. Lorimer, A. G. Lyne, B. Machenschalk, M. A. McLaughlin, M. Keith, M. Kramer, B. Machenschalk, C. Ng, M. A. Papa, P. S. Ray, X. Siemens, Einstein@Home discovery of four young gamma-ray pulsars in Fermi LAT data. Astrophys. J. 779, L11 (2013).

20. W. Atwood, A. Albert, L. Baldini, M. Tinivella, J. Bregeon, M. Pesce-Rollins, C. Sgrò, P. Bruel, E. Charles, A. Drlica-Wagner, A. Franckowiak, T. Jogler, L. Rochester, T. Usher, M. Wood, J. Cohen-Tanugi, S. Zimmer, Pass 8: Toward the Full Realization of the Fermi-LAT Scientific Potential, in Proceedings of the 4th Fermi Symposium, T. J. Brandt, N. Omodei, C. Wilson-Hodge, Eds., eConf C121028 (2012), arXiv:1303.3514.

21. P. M. Saz Parkinson, H. Xu, P. L. H. Yu, D. Salvetti, M. Marelli, A. D. Falcone, Classification and ranking of Fermi LAT gamma-ray sources from the 3FGL catalog using machine learning techniques. Astrophys. J. 820, 8 (2016).

22. H. J. Pletsch, L. Guillemot, B. Allen, D. Anderson, C. Aulbert, O. Bock, D. J. Champion, H. B. Eggenstein, H. Fehrmann, D. Hammer, R. Karupuspyas, M. Keith, M. Kramer, B. Machenschalk, C. Ng, M. A. Papa, P. S. Ray, X. Siemens, Einstein@Home discovery of four young gamma-ray pulsars in Fermi LAT data. Astrophys. J. 779, L11 (2013).

23. O. C. De Jager, J. W. H. Swanepoel, B. C. Raubenheimer, A powerful test for weak periodic signals with unknown light curve shape in sparse data. Astron. Astrophys. 221, 180–190 (1989).

24. J. M. Yao, R. N. Manchester, N. Wang, A new electron-density model for estimation of pulsar and FRB distances. Astrophys. J. 835, 29 (2017).

25. J. M. Cordes, T. J. W. Lazio, NE2001.1. A new model for the galactic distribution of free electrons and its fluctuations. arXiv:astro-ph/0207156 (2002).

26. R. N. Manchester, G. B. Hobbs, A. Teoh, M. Hobbs, The Australia Telescope National Facility Catalogue. Astron. J. 129, 1993 (2005).

27. M. Kramer, B. Stappers. Pulsar Science with the SKA, in *Advancing Astrophysics with the Square Kilometre Array (AASKA14)* (Proceedings of Science, Trieste, 2014), p. 36.

28. J. K. Swiggum, D. L. Kaplan, M. A. McLaughlin, D. R. Lorimer, S. Bogdanov, P. S. Ray, R. Lynch, P. Gentile, R. Rosen, S. A. Heatherly, B. N. Barlow, R. J. Hegedus, A. Vasquez Soto, P. Clancy, V. I. Kondratiev, K. Stovall, A. Istrate, B. Penprase, E. C. Bellm, A multiwavelength study of nearby millisecond pulsar PSR J1400–1431: Improved astrometry and an optical detection of its cool white dwarf companion. Astrophys. J. 847, 25 (2017).

29. R. Spiewak, M. Bailes, E. D. Barr, N. D. R. Bhat, M. Burgay, A. D. Cameron, D. J. Champion, C. M. L. Flynn, A. Jameson, S. Johnston, M. J. Keith, M. Kramer, S. R. Kulkarni, L. Levin, A. G. Lyne, V. Morello, C. Ng, A. Possenti, V. Ravi, B. W. Stappers, W. van Straten, C. Tiburzi, PSR J2322–2650—A low-luminosity millisecond pulsar with a planetary-mass companion. Mon. Not. R. Astron. Soc. 475, 469–477 (2018).

30. M. Marelli, A. De Luca, P. A. Caraveo, A multiwavelength study on the high-energy behavior of the Fermi/LAT pulsars. Astrophys. J. 733, 82 (2011).

31. T. J. Johnson, C. Venter, A. K. Harding, L. Guillemot, D. A. Smith, M. Kramer, Ö. Çelik, P. R. den Hartog, E. C. Ferrara, X. Hou, J. Lande, P. S. Ray, Constraints on the emission geometries and spin evolution of gamma-ray millisecond pulsars. Astrophys. J. Suppl. Ser. 213, 6 (2014).

32. C. Kalapotharakos, A. K. Harding, D. Kazanas. Gamma-ray emission in dissipative pulsar magnetospheres: From theory to Fermi observations. Astrophys. J. 793, 97 (2014).

33. B. Cerutti, A. A. Philippov, A. Spitkovsky, Modelling high-energy pulsar light curves from first principles. Mon. Not. R. Astron. Soc. 457, 2401–2414 (2016).

34. J. Takata, Y. Wang, K. S. Cheng, Population study for γ-ray pulsars with the outer gap model—II. Millisecond pulsars. Mon. Not. R. Astron. Soc. 414, 2173–2185 (2011).

35. N. Mirabal, E. Charles, E. C. Ferrara, P. L. Gonthier, A. K. Harding, M. A. Sánchez-Conde, D. J. Thompson, 3FGL demographics outside the galactic plane using supervised machine learning: Pulsar and dark matter subhalo interpretations. Astrophys. J. 825, 69 (2016).

36. M. Ackermann, A. Albert, B. Anderson, W. B. Atwood, L. Baldini, G. Barbiellini, D. Bastieri, K. Bechtol, R. Bellazzini, E. Bissaldi, R. D. Blandford, E. D. Bloom, R. Bonino, E. Bottacini, T. J. Brandt, J. Bregeon, P. Bruel, R. Buehler, G. A. Caliandro, R. A. Cameron, R. Caputo, M. Caragiulo, P. A. Caraveo, C. Cecchi, E. Charles, A. Chekhtman, J. Chiang, G. Chiaro, S. Ciprini, R. Claus, J. Cohen-Tanugi, J. Conrad, A. Cuoco, S. Cutini, F. D'Ammando, A. de Angelis, F. de Palma, R. Desiante, S. W. Digel, L. Di Venere, P. S. Drell, A. Drlica-Wagner, R. Essig, C. Favuzzi, S. J. Fegan, E. C. Ferrara, W. B. Focke, A. Franckowiak, Y. Fukazawa, S. Funk, P. Fusco, F. Gargano, D. Gasparrini, N. Giglietto, F. Giordano, M. Giroletti, T. Glanzman, G. Godfrey, G. A. Gomez-Vargas, I. A. Grenier, S. Guiriec, M. Gustafsson, E. Hays, J. W. Hewitt, D. Horan, T. Jogler, G. Johannesson, M. Kuss, S. Larsson, L. Latronico, J. Li, L. Li, M. Llena Garde, D. Longo, F. Loparco, P. Lubrano, D. Malyshev, M. Mayer, M. N. Mazziotta, J. E. McEnery, M. Meyer, P. F. Michelson, T. Mizuno, A. A. Moiseev, M. E. Monzani, A. Morselli, S. Murgia, E. Nuss, T. Ohsugi, M. Orienti, E. Orlando, J. F. Ormes, D. Paneque, J. H. Panetta, D. Salvetti, M. Sánchez-Conde, P. M. Saz Parkinson, A. Schulz, E. J. Siskind, A. A. Smith, F. Spada, G. Spandre, P. Spinelli, T. E. Stephens, A. W. Strong, D. J. Suson, H. Takahashi, T. Takahashi, Y. Tanaka, J. G. Thayer, J. B. Thayer, D. J. Thompson, L. Tibaldo, O. Tibolla, D. F. Torres, E. Torres, G. Tosti, E. Troja, B. Van Klaveren, G. Vianello, B. L. Winer, K. S. Wood,
M. Sánchez-Conde, A. Schulz, N. Sehgal, C. Sgrò, E. J. Siskind, F. Spada, G. Spandre, P. Spinelli, L. Strigari, H. Tajima, H. Takahashi, J. B. Thayer, L. Tibaldo, D. F. Torres, E. Troja, G. Vianello, M. Werner, B. L. Winer, K. S. Wood, M. Wood, G. Zaharijas, S. Zimmer, Fermi-LAT Collaboration, Searching for dark matter annihilation from Milky Way dwarf spheroidal galaxies with six years of Fermi Large Area Telescope data. *Phys. Rev. Lett.* **115**, 231301 (2015).

37. M. Ackermann, M. Ajello, A. Albert, W. B. Atwood, R. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, R. D. Blandford, E. D. Bloom, R. Bonino, E. Bottaci, T. J. Brandt, J. Bregeon, P. Bruel, R. Buehler, T. H. Burnett, R. A. Cameron, R. Caputo, M. Caragiulo, P. A. Caraveo, E. Cavazzuti, C. Cecchi, E. Charles, A. Chekhtman, J. Chiang, A. Chiappo, G. Chiaro, S. Ciprinì, J. Conrad, F. Costanza, A. Cuoco, S. Cutini, F. D'Ammando, F. de Palma, R. Desiante, S. W. Digel, N. Di Lalla, M. Di Mauro, L. Di Venere, P. S. Drell, C. Favuzzi, S. J. Fegan, E. C. Ferrara, W. B. Focke, A. Franckowiak, Y. Fukazawa, S. Funk, P. Fusco, F. Gargano, D. Gasparrini, N. Giglietto, F. Giordano, M. Giroletti, T. Glanzman, G. A. Gomez-Vargas, D. Green, I. A. Grenier, J. E. Grove, L. Guillemot, S. Guiriec, M. Gustafsson, A. K. Harding, E. Hays, J. W. Hewitt, D. Horan, T. Jogler, A. S. Johnson, T. Kamae, D. Kocevski, M. Kuss, G. La Mura, S. Larsson, L. Latronico, J. Li, F. Longo, F. Loparco, M. N. Lovellette, P. Lubrano, J. D. Magill, S. Maldera, D. Malyshev, A. Manfreda, P. Martin, M. N. Mazziotta, P. F. Michelson, N. Mirabal, W. Mitthumsiri, T. Mizuno, A. A. Moiseev, M. E. Monzani, A. Morselli, M. Negro, E. Nuss, T. Ohsugi, M. Orienti, E. Orlando, J. F. Ormes, D. Paneque, M. Perkins, M. Persic, M. Pesce-Rollins, F. Piron, G. Principe, S. Rainò, R. Rando, M. Razzano, S. Razzaque, A. Reimer, O. Reimer, M. Sánchez-Conde, C. Sgrò, D. Simone, E. J. Siskind, F. Spada, G. Spandre, P. Spinelli, D. J. Suson, H. Tajima, K. Tanaka, J. B. Thayer, L. Tibaldo, D. F. Torres, E. Troja, Y. Uchiyama, G. Vianello, K. S. Wood, M. Wood, G. Zaharijas, S. Zimmer, The Fermi Galactic center GeV excess and implications for dark matter. *Astrophys. J.* **840**, 43 (2017).

38. R. Bartels, S. Krishnamurthy, C. Weniger, Strong support for the millisecond pulsar origin of the galactic center GeV excess. *Phys. Rev. Lett.* **116**, 051102 (2016).

39. S. K. Lee, M. Lisanti, B. R. Safdi, T. R. Slatyer, W. Xue, Evidence for unresolved γ-ray point sources in the inner galaxy. *Phys. Rev. Lett.* **116**, 051103 (2016).

40. F. Calore, M. Di Mauro, F. Donato, J. W. T. Hessels, C. Weniger, Radio detection prospects for a bulge population of millisecond pulsars as suggested by Fermi-LAT observations of the inner galaxy. *Astrophys. J.* **827**, 143 (2016).

41. M. Kerr, Improving sensitivity to weak pulsations with photon probability weighting. *Astrophys. J.* **732**, 38 (2011).

42. M. Kerr, "Likelihood methods for the detection and characterization of gamma-ray pulsars with the Fermi Large Area Telescope," thesis, University of Washington (2010).

43. F. Acero, M. Ackermann, M. Ajello, A. Albert, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, E. D. Bloom, R. Bonino, E. Bottaci, T. J. Brandt, J. Bregeon, P. Bruel, R. Buehler, S. Buson, G. A. Caliandro, R. A. Cameron, M. Caragiulo, P. A. Caraveo, J. M. Casandjia, E. Cavazzuti, C. Cecchi, E. Charles, A. Chekhtman, J. Chiang, G. Chiaro, S. Ciprinì, R. Claus, J. Cohen-Tanugi, J. Conrad, A. Cuoco, S. Cutini, F. D'Ammando, A. de Angelis, F. de Palma, R. Desiante, S. W. Digel, L. Di Venere, P. S. Drell, C. Favuzzi, S. J. Fegan, E. C. Ferrara, W. B. Focke, A. Franckowiak, S. Funk, P. Fusco, F. Gargano, D. Gasparrini, N. Giglietto, F. Giordano, M. Giroletti, T. Glanzman, G. Godfrey, I. A. Grenier, S. Guiriec, D. Hadasch, A. K. Harding, K. Hayashi, E. Hays, J. W. Hewitt, A. B. Hill, D. Horan, X. Hou, T. Jogler, G. Jóhannesson, T. Kamae, M. Kuss, D. Landriu, S. Larsson, L. Latronico, J. Li, L. Li, F. Longo, F. Loparco, M. N. Lovellette, P. Lubrano, S. Maldera, D. Malyshev, A. Manfreda, P. Martin, M. Mayer, M. N. Mazziotta, J. E. McEnery, P. F. Michelson, N. Mirabal, T. Mizuno, M. E. Monzani, A. Morselli, E. Nuss, T. Ohsugi, N. Omodei, M. Orienti, E. Orlando, J. F. Ormes, D. Paneque, M. Pesce-Rollins, F. Piron, G. Pivato, S. Rainò, R. Rando, M. Razzano, S. Razzaque, A. Reimer, O. Reimer, Q. Remy, N. Renault, M. Sánchez-Conde, M. Schaal, A. Schulz, C. Sgrò, E. J. Siskind, F. Spada, G. Spandre, P. Spinelli, A. W. Strong, D. J. Suson, H. Tajima, H. Takahashi, J. B. Thayer, D. J. Thompson, L. Tibaldo, M. Tinivella, D. F. Torres, G. Tosti, E. Troja, G. Vianello, M. Werner, K. S. Wood, M. Wood, G. Zaharijas, S. Zimmer, Development of the model of galactic interstellar emission for standard point-source analysis of Fermi large area telescope data. *Astron. Astrophys. Suppl. Ser.* **223**, 26 (2016).

44. J. D. Scargle, J. P. Norris, B. Jackson, J. Chiang, Studies in astronomical time series analysis. VI. Bayesian block representations. *Astrophys. J.* **764**, 167 (2013).

45. G. Schwarz, Estimating the dimension of a model. *Ann. Stat.* **6**, 461–464 (1978).

46. D. Foreman-Mackey, D. W. Hogg, D. Lang, J. Goodman, emcee: The MCMC Hammer. *Publ. Astron. Soc. Pac.* **125**, 306 (2013).

47. J. Goodman, Weare, Ensemble samplers with affine invariance. *Comm. App. Math. Comp. Sci.* **5**, 65–80 (2010).

48. I. S. Shklovskii, Possible causes of the secular increase in pulsar periods. *Soviet Astron.* **13**, 562 (1970).

49. K. P. Watters, R. W. Romani, P. Weltevrede, S. Johnston, An atlas for interpreting γ-ray pulsar light curves. *Astrophys. J.* **695**, 1289 (2009).

50. L. Guillemot, D. A. Smith, H. Laffon, G. H. Janssen, I. Cognard, G. Theureau, G. Desvignes, E. C. Ferrara, P. S. Ray, The gamma-ray millisecond pulsar deathline, revisited—New velocity and distance measurements. *Astron. Astrophys.* **587**, A109 (2016).

51. S. M. Ransom, "New search techniques for binary pulsars," thesis, Harvard University (2001).

52. M. B. Mickaliger, D. R. Lorimer, J. Boyles, M. A. McLaughlin, A. Collins, L. Hough, N. Tehrani, C. Tenney, A. Liska, J. Swiggum, Discovery of five new pulsars in archival data. *Astrophys. J.* **759**, 127 (2012).

53. M. Pierbattista, A. K. Harding, I. A. Grenier, T. J. Johnson, P. A. Caraveo, M. Kerr, P. L. Gonthier, Light-curve modelling constraints on the obliquities and aspect angles of the young Fermi pulsars. *Astron. Astrophys.* **575**, A3 (2015).

**Acknowledgments:** We are grateful to the thousands of volunteers who have supported Einstein@Home. We also thank those volunteers whose computers discovered the new MSPs: PSR J1035−6720: “Ws5”; K. Kovacs, of Seattle, WA, USA, and the ATLAS Cluster, AEI, Hannover, Germany; PSR J1744−7619: D. Hoberer, of Gainesville, TX, USA, and the ATLAS Cluster, AEI, Hannover, Germany; I. Yakushin of Chicago, IL, USA, and the LIGO Laboratory, USA; and K. Pickstone of Oldham, UK. **Funding:** This work was supported by the Max-Planck-Gesellschaft, the Deutsche Forschungsgemeinschaft through an Emmy Noether Research Grant [no. PL 710/1-1 (H.J.P.)], and NSF award 1104902. C.J.C. acknowledges support from the European Research Council under the European Union’s Horizon 2020 research and innovation programme (grant agreement no. 715051; Spiders). W.M. was partially supported by the Thailand Research Fund (grants TRG5880173 and RTA5980003). M.R. was funded by the Italian Ministry of Education, University and Research (contract no. FIRB-2012-RBRF12PM1F). Work at the Naval Research Laboratory is supported by NASA. The Parkes Observatory is part of the Australia Telescope, funded by the Commonwealth of Australia for operation as a National Facility managed by the Commonwealth Scientific and Industrial Research Organisation (CSIRO). 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 e 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-76FS00515. **Author contributions:** C.J.C., H.J.P., J.W., and L.G. designed and implemented the Einstein@Home gamma-ray pulsar searches. M. Kerr performed follow-up searches with the Parkes radio telescope. T.J.J. carried out the gamma-ray pulse profile modeling. F.C. reanalyzed archival Parkes observations. D.S. performed analyses of x-ray observations. The manuscript was written by the aforementioned authors. The remaining authors contributed to the development and operation of the Einstein@Home project, to the design and operation of the Fermi LAT and analysis of its observations, or to the Fermi Pulsar Search Consortium. **Competing interests:** The authors declare that they have no competing interests. **Data and materials availability:** Fermi-LAT data are available from the Fermi Science Support Center (http://fermi.gsc.nasa.gov/ssc). XMM-Newton data are available from the XMM-Newton Science Archive (http://nxsa.esac.esa.int). Parkes radio telescope data recorded for this paper will be made available from the CSIRO Data Access Portal (https://data.csiro.au/) 18 months after the observation dates and are available upon request before that. All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

Submitted 18 August 2017
Accepted 24 January 2018
Published 28 February 2018
10.1126/sciadv.aa07228

**Citation:** C. J. Clark, H. J. Pletsch, J. Wu, L. Guillemot, M. Kerr, T. J. Johnson, F. Camilo, D. Salvetti, B. Allen, D. Anderson, C. Aulbert, C. Beer, O. Bock, A. Cuéllar, H.-B. Eggenstein, H. Fehrmann, M. Kramer, S. A. Kwang, B. Machenschalk, L. Nieder, M. Ackermann, M. Ajello, L. Baldini, J. Ballet, G. Barbiellini, D. Bastieri, R. Bellazzini, E. Bissaldi, R. D. Blandford, E. D. Bloom, R. Bonino, E. Bottaci, T. J. Brandt, J. Bregeon, P. Bruel, R. Buehler, T. H. Burnett, S. Buson, R. A. Cameron, R. Caputo, P. A. Caraveo, E. Cavazzuti, C. Cecchi, E. Charles, A. Chekhtman, S. Ciprinì, L. R. Cominsky, D. Costantin, S. Cutini, F. D'Ammando, A. De Luca, R. Desiante, L. Di Venere, M. Di Mauro, N. Di Lalla, S. W. Digel, C. Favuzzi, E. C. Ferrara, A. Franckowiak, Y. Fukazawa, S. Funk, P. Fusco, F. Gargano, D. Gasparrini, N. Giglietto, F. Giordano, M. Giroletti, G. A. Gomez-Vargas, D. Green, I. A. Grenier, S. Guiriec, A. K. Harding, J. W. Hewitt, D. Horan, G. Jóhannesson, S. Kensei, M. Kuss, G. La Mura, S. Larsson, L. Latronico, J. Li, F. Longo, F. Loparco, M. N. Lovellette, P. Lubrano, J. D. Magill, S. Maldera, A. Manfreda, P. S. Paliya, F. de Palma, D. Paneque, J. S. Perkins, M. Persic, M. Pesce-Rollins, T. A. Porter, G. Principe, S. Rainò, R. Rando, P. S. Ray, M. Razzano, A. Reimer, O. Reimer, R. W. Romani, P. M. Saz Parkinson, C. Sgrò, E. J. Siskind, D. A. Smith, F. Spada, G. Spandre, P. Spinelli, J. B. Thayer, D. J. Thompson, D. F. Torres, E. Troja, G. Vianello, K. Wood, M. Wood, Einstein@Home discovers a radio-quiet gamma-ray millisecond pulsar. *Sci. Adv.* **4**, ea007228 (2018).
ScienceAdvances

Einstein@Home discovers a radio-quiet gamma-ray millisecond pulsar

Colin J. Clark, Holger J. Pletsch, Jason Wu, Lucas Guillemot, Matthew Kerr, Tyrel J. Johnson, Fernando Camilo, David Salvetti, Bruce Allen, David Anderson, Carsten Aulbert, Christian Beer, Oliver Bock, Andres Cuéllar, Heinz-Bernd Eggenstein, Henning Fehrmann, Michael Kramer, Shawn A. Kwang, Bernd Machenschalk, Lars Nieder, Markus Ackermann, Marco Ajello, Luca Baldini, Jean Ballet, Guido Barbiellini, Denis Bastieri, Ronaldo Bellazzini, Elisabetta Bissaldi, Roger D. Blandford, Elliott D. Bloom, Raffaella Bonino, Eugenio Bottacini, Terri J. Brandt, Johan Bregeon, Philippe Bruel, Rolf Buehler, Toby H. Burnett, Sara Buson, Rob A. Cameron, Regina Caputo, Patrizia A. Caraveo, Elisabetta Cavazzuti, Claudia Cecchi, Eric Charles, Alexandre Chekhtman, Stefano Ciprini, Lynn R. Cominsky, Denise Costantin, Sara Cutini, Filippo D'Ammando, Andrea De Luca, Rachele Desiante, Leonardo Di Venere, Mattia Di Mauro, Niccolò Di Lalla, Seth W. Digel, Cecilia Favuzzi, Elizabeth C. Ferrara, Anna Franckowiak, Yasushi Fukazawa, Stefan Funk, Piergiorgio Fusco, Fabio Gargano, Dario Gasparrini, Nico Giglietto, Francesco Giordano, Marcello Giroletti, Germán A. Gomez-Vargas, David Green, Isabelle A. Grenier, Sylvain Guiriec, Alice K. Harding, John W. Hewitt, Deirdre Horan, Guðlaugur Jóhannesson, Shiki Kensei, Michael Kuss, Giovanni La Mura, Stefan Larsson, Luca Latronico, Jian Li, Francesco Longo, Francesco Loparco, Michael N. Lovellette, Pasquale Lubrano, Jeffrey D. Magill, Simone Maldera, Alberto Manfreda, Mario N. Mazziotta, Julie E. McEnery, Peter F. Michelson, Nestor Mirabal, Warit Mitthumsiri, Tsunefumi Mizuno, Maria Elena Monzani, Aldo Morselli, Igor V. Moskalenko, Eric Nuss, Takashi Ohsugi, Nicola Omodei, Monica Orienti, Elena Orlando, Michele Palatiello, Vaidehi S. Paliya, Francesco de Palma, David Paneque, Jeremy S. Perkins, Massimo Persic, Melissa Pesce-Rollins, Troy A. Porter, Giacomo Principe, Silvia Rainò, Riccardo Rando, Paul S. Ray, Massimiliano Razzano, Anita Reimer, Olaf Reimer, Roger W. Romani, Pablo M. Saz Parkinson, Carmelo Sgrò, Eric J. Siskind, David A. Smith, Francesca Spada, Gloria Spandre, Paolo Spinelli, Jana B. Thayer, David J. Thompson, Diego F. Torres, Eleonora Troja, Giacomo Vianello, Kent Wood and Matthew Wood

Sci Adv 4 (2), eaao7228.

DOI: 10.1126/sciadv.aao7228

ARTICLE TOOLS

http://advances.sciencemag.org/content/4/2/eaao7228

REFERENCES

This article cites 47 articles, 5 of which you can access for free
http://advances.sciencemag.org/content/4/2/eaao7228#BIBL

PERMISSIONS

http://www.sciencemag.org/help/reprints-and-permissions