Millisecond pulsars (MSPs) are a growing class of gamma-ray emitters. Pulsed gamma-ray signals have been detected from more than 40 MSPs with the Fermi Large Area Telescope (LAT). The wider radio beams and more compact magnetospheres of MSPs enable studies of emission geometries over a broader range of phase space than non-recycled radio-loud gamma-ray pulsars. We have modeled the gamma-ray light curves of 40 LAT-detected MSPs using geometric emission models assuming a vacuum retarded-dipole magnetic field. We modeled the radio profiles using a single-altitude hollow-cone beam, with a core component when indicated by polarimetry; however, for MSPs with gamma-ray and radio light curve peaks occurring at nearly the same rotational phase, we assume that the radio emission is co-located with the gamma rays and caustic in nature. The best-fit parameters and confidence intervals are determined using a maximum likelihood technique. We divide the light curves into three model classes, with gamma-ray peaks trailing (Class I), aligned (Class II), or leading (Class III) the radio peaks. Outer gap and slot gap (two-pole caustic) models best fit roughly equal numbers of Class I and II, while Class III are exclusively fit with pair-starved polar cap models. Distinguishing between the model classes based on typical derived parameters is difficult. We explore the evolution of the magnetic inclination angle with period and spin-down power, finding possible correlations. While the presence of significant off-peak emission can often be used as a discriminator between outer gap and slot gap models, a hybrid model may be needed.

Cosmic rays are particles (mostly protons) accelerated to relativistic speeds. Despite wide agreement that supernova remnants (SNRs) are the sources of galactic cosmic rays, unequivocal evidence for the acceleration of protons in these objects is still lacking. When accelerated protons encounter interstellar material, they produce neutral pions, which in turn decay into gamma rays. This offers a compelling way to detect the acceleration sites of protons. The identification of pion-decay gamma rays has been difficult because high-energy electrons also produce gamma rays via bremsstrahlung and inverse Compton scattering. We detected the characteristic pion-decay feature in the gamma-ray spectra of two SNRs, IC 443 and W44, with the Fermi Large Area Telescope. This detection provides direct evidence that cosmic-ray protons are accelerated in SNRs.

We report on the discovery of gamma-ray pulsations from five millisecond pulsars (MSPs) using the Fermi Large Area Telescope (LAT) and timing ephemerides provided by various radio observatories. We also present confirmation of the gamma-ray pulsations from a sixth source, PSR J2051-0827. Five of these six MSPs are in binary systems: PSRs J1713+0747, J1741+1351, J1600-3053 and the two black widow binary pulsars PSRs J0610-2100 and 2051-0827. The only isolated MSP is the nearby PSR J1024-0719, which is also known to emit X-rays. We present X-ray observations in the direction of PSRs J1600-3053 and J2051-0827. While the latter is firmly detected, we an only give upper limits for the X-ray flux of the former. There are no dedicated X-ray observations available for the other 3 objects. The MSPs mentioned above, together with most of the MSPs detected by Fermi, are used to put together a sample of 30 gamma-ray MSPs which is used to study the morphology and phase connection of radio and gamma-ray pulse profiles. We show that MSPs with pulsed gamma-ray emission which is phase aligned with the radio emission present the steepest radio spectra and the largest magnetic fields at the light cylinder among all MSPs. As well, we also observe a trend towards very low, or undetectable, radio linear polarisation levels. These properties could be attributed to caustic radio emission produced at a range of different altitudes in the magnetosphere. We note that most of these characteristics are also observed in the Crab pulsar, the only other radio pulsar known to exhibit phase-aligned radio and gamma-ray emission.

Using the 100-m Effelsberg radio telescope operating at 1.36 GHz, we have performed a targeted radio pulsar survey of 289 unassociated gamma-ray sources discovered by the Large Area Telescope (LAT) aboard the Fermi satellite and published in the 1FGL catalogue (Abdo et al., 2010). This survey resulted in the discovery of millisecond pulsar J1745+1017, which resides in a short-period binary system with a low-mass companion, Mmin ~ 0.0137 Msun, indicative of `Black Widow' type systems. A two-year timing campaign has produced a refined radio ephemeris, accurate enough to allow for phase-folding of the LAT photons, resulting in the detection of a dual-peaked gamma-ray light-curve, proving that PSR J1745+1017 is the source responsible for the gamma-ray emission seen in 1FGL J1745.5 + 1018 (2FGL J1745.6+1015; Nolan et al., 2012). We find the gamma-ray spectrum of PSR J1745+1017 to be well modelled by an exponentially-cut-off power law with cut-off energy 3.2 GeV and photon index 1.6. The observed sources are known to contain a further 10 newly discovered pulsars which were undetected in this survey. Our radio observations of these sources are discussed and in all cases limiting flux densities are calculated. The reasons behind the seemingly low yield of discoveries are also discussed.

The Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope (formerly GLAST), with its improved sensitivity relative to previous generation gamma-ray telescopes, is significantly increasing the number of known gamma-ray sources in the sky, including pulsars. In addition to searching for gamma-ray emission from known radio pulsars, it is now possible to successfully perform blind searches for pulsars on the gamma-ray data alone, with the goal of uncovering a new population of potentially radio-quiet (Geminga-like) pulsars. We describe our methods and some recent results from our searches, performed using the time-differencing technique, which have resulted in the discovery of a large number of new gamma-ray pulsars.

The Fermi Large Area Telescope (Fermi-LAT, hereafter LAT), the primary instrument on the Fermi Gamma-ray Space Telescope (Fermi) mission, is an imaging, wide field-of-view, high-energy γ-ray telescope, covering the energy range from 20 MeV to more than 300 GeV. During the first years of the mission the LAT team has gained considerable insight into the in-flight performance of the instrument. Accordingly, we have updated the analysis used to reduce LAT data for public release as well as the Instrument Response Functions (IRFs), the description of the instrument performance provided for data analysis. In this paper we describe the effects that motivated these updates. Furthermore, we discuss how we originally derived IRFs from Monte Carlo simulations and later corrected those IRFs for discrepancies observed between flight and simulated data. We also give details of the validations performed using flight data and quantify the residual uncertainties in the IRFs. Finally, we describe techniques the LAT team has developed to propagate those uncertainties into estimates of the systematic errors on common measurements such as fluxes and spectra of astrophysical sources.

The Fermi Large Area Telescope (LAT) First Source Catalog (1FGL) provided spatial, spectral, and temporal properties for a large number of gamma-ray sources using a uniform analysis method. After correlating with the most-complete catalogs of source types known to emit gamma rays, 630 of these sources are "unassociated" (i.e., have no obvious counterparts at other wavelengths). Here, we employ two statistical analyses of the primary gamma-ray characteristics for these unassociated sources in an effort to correlate their gamma-ray properties with the active galactic nucleus (AGN) and pulsar populations in 1FGL. Based on the correlation results, we classify 221 AGN-like and 134 pulsar-like sources in the 1FGL unassociated sources. The results of these source "classifications" appear to match the expected source distributions, especially at high Galactic latitudes. While useful for planning future multiwavelength follow-up observations, these analyses use limited inputs, and their predictions should not be considered equivalent to "probable source classes" for these sources. We discuss multiwavelength results and catalog cross-correlations to date, and provide new source associations for 229 Fermi-LAT sources that had no association listed in the 1FGL catalog. By validating the source classifications against these new associations, we find that the new association matches the predicted source class in similar to 80% of the sources.

Gamma-ray binaries are stellar systems containing a neutron star or black hole, with gamma-ray emission produced by an interaction between the components. These systems are rare, even though binary evolution models predict dozens in our Galaxy. A search for gamma-ray binaries with the Fermi Large Area Telescope (LAT) shows that 1FGL J1018.6-5856 exhibits intensity and spectral modulation with a 16.6-day period. We identified a variable x-ray counterpart, which shows a sharp maximum coinciding with maximum gamma-ray emission, as well as an O6V((f)) star optical counterpart and a radio counterpart that is also apparently modulated on the orbital period. 1FGL J1018.6-5856 is thus a gamma-ray binary, and its detection suggests the presence of other fainter binaries in the Galaxy.

Millisecond pulsars, old neutron stars spun-up by accreting matter from a companion star, can reach high rotation rates of hundreds of revolutions per second. Until now, all such "recycled" rotation-powered pulsars have been detected by their spin-modulated radio emission. In a computing-intensive blind search of gamma-ray data from the Fermi Large Area Telescope (with partial constraints from optical data), we detected a 2.5-millisecond pulsar, PSR J1311-3430. This unambiguously explains a formerly unidentified gamma-ray source that had been a decade-long enigma, confirming previous conjectures. The pulsar is in a circular orbit with an orbital period of only 93 minutes, the shortest of any spin-powered pulsar binary ever found.

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

A young and energetic pulsar powers the well-known Crab Nebula. Here we describe two separate gamma-ray (photon energy >100 MeV) flares from this source detected by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. The first flare occurred in February 2009 and lasted approximately 16 days. The second flare was detected in September 2010 and lasted approximately 4 days. During these outbursts the gamma-ray flux from the nebula increased by factors of four and six, respectively. The brevity of the flares implies that the gamma rays were emitted via synchrotron radiation from PeV (10^15 eV) electrons in a region smaller than 1.4 10^-2 pc. These are the highest energy particles that can be associated with a discrete astronomical source, and they pose challenges to particle acceleration theory.

We report the discovery of nine previously unknown gamma-ray pulsars in a blind search of data from the Fermi Large Area Telescope (LAT). The pulsars were found with a novel hierarchical search method originally developed for detecting continuous gravitational waves from rapidly rotating neutron stars. Designed to find isolated pulsars spinning at up to kHz frequencies, the new method is computationally efficient, and incorporates several advances, including a metric-based gridding of the search parameter space (frequency, frequency derivative and sky location) and the use of photon probability weights. The nine pulsars have spin frequencies between 3 and 12 Hz, and characteristic ages ranging from 17 kyr to 3 Myr. Two of them, PSRs J1803–2149 and J2111+4606, are young and energetic Galactic-plane pulsars (spin-down power above 6×1035 erg s−1 and ages below 100 kyr). The seven remaining pulsars, PSRs J0106+4855, J0622+3749, J1620–4927, J1746–3239, J2028+3332, J2030+4415, J2139+4716, are older and less energetic; two of them are located at higher Galactic latitudes (|b|> 10◦). PSR J0106+4855 has the largest characteristic age (3 Myr) and the smallest surface magnetic field (2×1011G) of all LAT blind-search pulsars. PSR J2139+4716 has the lowest spin-down power (3×1033 erg s−1) among all non-recycled gamma-ray pulsars ever found. Despite extensive multi-frequency observations, only PSR J0106+4855 has detectable pulsations in the radio band. The other eight pulsars belong to the increasing population of radio-quiet gamma-ray pulsars. Subject headings: gamma rays: stars, pulsars: general, pulsars: individual (PSR J0106+4855, PSR J0622+3749, PSR J1620–4927, PSR J1746–3239, PSR J1803–2149, PSR J2028+3332, PSR J2030+4415, PSR J2111+4606, PSR J2139+4716) 1 Albert-Einstein-Institut, Max-Planck-Institut für Gravitationsphysik, D-30167 Hannover, Germany 2 Institut für Gravitationsphysik, Leibniz Universität Hannover, D30167 Hannover, Germany 3 email: holger.pletsch@aei.mpg.de 4 Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany 5 email: guillemo@mpifr-bonn.mpg.de 6 Physics Department, University of Wisconsin – Milwaukee, Milwaukee, WI 53211, USA 7 Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, The University of Manchester, M13 9PL, UK 8 Space Science Division, Naval Research Laboratory, Washington, DC 20375-5352, USA 9 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 10 Università degli Studi di Pavia, 27100 Pavia, Italy 11 INAF-Istituto di Astrofisica Spaziale e Fisica Cosmica, I-20133 Milano, Italy 12 Columbia Astrophysics Laboratory, Columbia University, New York, NY 10027, USA 13 NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 14 Center for Research and Exploration in Space Science and Technology (CRESST) and NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 15 Department of Physics and Center for Space Sciences and Technology, University of Maryland Baltimore County, Baltimore, MD 21250, USA 16 Astronomical Institute “Anton Pannekoek,” University of Amsterdam, Postbus 94249, 1090 GE Amsterdam, Netherlands 17 Netherlands Institute for Radio Astronomy (ASTRON), Postbus 2, 7990 AA Dwingeloo, Netherlands 18 CSIRO Astronomy and Space Science, Australia Telescope National Facility, Epping NSW 1710, Australia 19 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 20 Istituto Universitario di Studi Superiori (IUSS), I-27100 Pavia, Italy 21 Department of Physics, West Virginia University, Morgantown, WV 26506, USA 22 Center for Earth Observing and Space Research, College of Science, George Mason University, Fairfax, VA 22030, resident at Naval Research Laboratory, Washington, DC 20375, USA 23 National Radio Astronomy Observatory (NRAO), Charlottesville, VA 22903, USA 24 Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, I-56127 Pisa, Italy ar X iv :1 11 1. 05 23 v1 [ as tr oph .H E ] 2 N ov 2 01 1

The high sensitivity of the Fermi-LAT (Large Area Telescope) offers the first opportunity to study faint and extended GeV sources such as pulsar wind nebulae (PWNe). After one year of observation the LAT detected and identified three PWNe: the Crab Nebula, Vela-X, and the PWN inside MSH 15–52. In the meantime, the list of LAT detected pulsars increased steadily. These pulsars are characterized by high energy loss rates () from ∼3 × 1033 erg s−1 to 5 × 1038 erg s−1 and are therefore likely to power a PWN. This paper summarizes the search for PWNe in the off-pulse windows of 54 LAT-detected pulsars using 16 months of survey observations. Ten sources show significant emission, seven of these likely being of magnetospheric origin. The detection of significant emission in the off-pulse interval offers new constraints on the γ-ray emitting regions in pulsar magnetospheres. The three other sources with significant emission are the Crab Nebula, Vela-X, and a new PWN candidate associated with the LAT pulsar PSR J1023−5746, coincident with the TeV source HESS J1023−575. We further explore the association between the HESS and the Fermi source by modeling its spectral energy distribution. Flux upper limits derived for the 44 remaining sources are used to provide new constraints on famous PWNe that have been detected at keV and/or TeV energies.

We report on the first year of Fermi γ-ray observations of pulsed high-energy emission from the old PSR J2043 + 2740. The study of the γ-ray efficiency of such old pulsars gives us an insight into the evolution of pulsars’ ability to emit in γ rays as they age. The γ-ray light curve of this pulsar above 0.1 GeV is clearly defined by two sharp peaks, 0.353 ± 0.035 periods apart. We have combined the γ-ray profile characteristics of PSR J2043 + 2740 with the geometrical properties of the pulsar's radio emission, derived from radio-polarization data, and constrained the pulsar-beam geometry in the framework of a two-pole caustic (TPC) and an outer gap (OG) model. The ranges of magnetic inclination and viewing angle were determined to be {α, ζ} ∼ {52°–57°, 61°–68°} for the TPC model, and {α, ζ} ∼ {62°–73°, 74°–81°} and {α, ζ} ∼ {72°–83°, 60°–75°} for the OG model. Based on this geometry, we assess possible birth locations for this pulsar and derive a likely proper motion, sufficiently high to be measurable with VLBI. At a characteristic age of 1.2 Myr, PSR J2043 + 2740 is the third oldest of all discovered, non-recycled, γ-ray pulsars: it is twice as old as the next oldest, PSR J0357 + 32, and younger only than the recently discovered PSR J1836 + 5925 and PSR J2055 + 25, both of which are at least five and ten times less energetic, respectively.

We report on the gamma-ray activity of the blazar Mrk 501 during the first 480 days of Fermi operation. We find that the average LAT gamma-ray spectrum of Mrk 501 can be well described by a single power-law function with a photon index of 1.78 +/- 0.03. While we observe relatively mild flux variations with the Fermi-LAT (within less than a factor of 2), we detect remarkable spectral variability where the hardest observed spectral index within the LAT energy range is 1.52 +/- 0.14, and the softest one is 2.51 +/- 0.20. These unexpected spectral changes do not correlate with the measured flux variations above 0.3GeV. In this paper, we also present the first results from the 4.5-month-long multifrequency campaign (2009 March 15 - August 1) on Mrk 501, which included the VLBA, Swift, RXTE, MAGIC and VERITAS, the F-GAMMA, GASP-WEBT, and other collaborations and instruments which provided excellent temporal and energy coverage of the source throughout the entire campaign. The average spectral energy distribution of Mrk 501 is well described by the standard one-zone synchrotron self-Compton model. In the framework of this model, we find that the dominant emission region is characterized by a size <~ 0.1 pc (comparable within a factor of few to the size of the partially-resolved VLBA core at 15-43 GHz), and that the total jet power (~10^{44} erg s^{-1}) constitutes only a small fraction (~10^{-3}) of the Eddington luminosity. The energy distribution of the freshly-accelerated radiating electrons required to fit the time-averaged data has a broken power-law form in the energy range 0.3GeV-10TeV, with spectral indices 2.2 and 2.7 below and above the break energy of 20GeV. We argue that such a form is consistent with a scenario in which the bulk of the energy dissipation within the dominant emission zone of Mrk 501 is due to relativistic, proton-mediated shocks.

The BL Lacertae object 3C 66A was detected in a flaring state by the Fermi Large Area Telescope (LAT) and VERITAS in 2008 October. In addition to these gamma-ray observations, F-GAMMA, GASP-WEBT, PAIRITEL, MDM, ATOM, Swift, and Chandra provided radio to X-ray coverage. The available light curves show variability and, in particular, correlated flares are observed in the optical and Fermi-LAT gamma-ray band. The resulting spectral energy distribution can be well fitted using standard leptonic models with and without an external radiation field for inverse Compton scattering. It is found, however, that only the model with an external radiation field can accommodate the intra-night variability observed at optical wavelengths.

Contrary to expectations, the γ -rays from a distant cluster of stars are dominated by emission from a single neutron star.